Journal of Heredity Advance Access originally published online on December 7, 2006
Journal of Heredity 2007 98(2):158-164; doi:10.1093/jhered/esl058
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Population Growth Confounds Phylogeographic Inference in Namaqua Sandgrouse
From the DST-NRF Centre of Excellence at the Percy FitzPatrick Institute, Molecular Ecology and Evolution Programme, Department of Genetics, University of Pretoria, Pretoria 0002, South Africa (Delport and Bloomer); and the DST-NRF Centre of Excellence at the Percy FitzPatrick Institute of African Ornithology, University of Cape Town, Rondebosch 7701, South Africa (Crowe and Lloyd). Wayne Delport is now at room 2.04, John Day Zoology Building, University Avenue, University of Cape Town, Rondebosch, Cape Town, 7701, South Africa
Address correspondence to W. Delport at the address above, or e-mail: wdelport{at}botzoo.uct.ac.za.
The Namaqua sandgrouse, Pterocles namaqua, is a highly nomadic granivore of semiarid to arid habitats. As a result of nomadic movements in response to rainfall, the size of the breeding population in any one area fluctuates dramatically between breeding seasons. This high mobility in response to spatial and temporal abundance of food resources is expected to result in little population genetic structuring. Namaqua sandgrouse also shows a seasonally predictable partial migration between the southeast and northwest regions of South Africa, and a further possible north–south migration between southwestern South Africa and central Namibia. It is unclear whether birds migrating between these regions breed in only one or both regions. If populations breed in only one region of their migratory range, then population genetic structuring is predicted to occur. This study addresses Namaqua sandgrouse movements with the analysis of mitochondrial DNA control region sequences. In general, little population genetic structure was evident, yet strong signals of population growth were detected. Several populations have private alleles, which is in direct contradiction to the spatial genetic pattern expected under high levels of gene flow. We suggest that the inference of high levels of female gene flow could be an artifact of population growth and that additional loci will allow a greater understanding of Namaqua sandgrouse movements.
The inference of the demographic processes of migration, admixture, isolation by distance, bottlenecks, and population growth has become a central challenge in phylogeographic studies (Emerson et al. 2001; Knowles and Maddison 2002). Such demographic processes, if acting in isolation, would each be expected to generate different spatial genetic patterns. However, each of these processes may occur during limited time periods, or multiple processes may occur simultaneously, and thus, the untangling of multiple historical factors contributing to spatial genetic pattern remains a major challenge (Wakeley and Hey 1998; Knowles 2004). Recent developments in coalescent approaches (Kuhner et al. 1998; Bahlo and Griffiths 2000; Beerli and Felsenstein 2001; Nielsen and Wakeley 2001) have, however, allowed the coestimation of demographic parameters of interest, such as population growth, effective population size, migration, time of isolation, among other parameters. In this article, we use such available computational methods to identify demographic processes that have contributed to observed spatial genetic patterns in the partially migratory and nomadic Namaqua sandgrouse (Pterocles namaqua).
Namaqua sandgrouse inhabit arid to semiarid regions of southern Africa (Figure 1). As obligate granivores, they must drink regularly and are known to congregate in flocks during a synchronized flight to the nearest surface water each morning (Ward 1972). These flocks may cover a distance of more than 50 km a day in their movement between water, feeding, and roosting sites (Lloyd et al. 2000). Namaqua sandgrouse are highly nomadic over most of their distribution, exhibiting large interannual fluctuations in abundance in any one area (Lloyd, Little, Crowe, and Simmons 2001). In addition, this species has a partial migration, restricted to the southern African sub-region, between the southwestern and northeastern regions of its distribution (Malan et al. 1994; Lloyd, Little, Crowe, and Simmons 2001). Populations in the southwest (Namaqualand/Karoo) are generally larger in the summer months (September–April) and often disappear during the winter months (May–August), whereas populations in the northeast (southern Kalahari) are larger in the winter months (Malan et al. 1994; Lloyd, Little, Crowe, and Simmons 2001). This suggests a partial migration between these regions, with a large proportion of birds from both regions residing in the northeast and southwest, during winter and summer, respectively. Breeding occurs in both regions, but it is unclear whether birds moving between these regions breed in only one or both regions (Lloyd, Little, Crowe, and Simmons 2001).
|
The high mobility of Namaqua sandgrouse is expected to facilitate a high level of gene flow between populations, resulting in little genetic structuring, unless these movements are structured in some manner. Clancey (1979) suggested that, during their nonbreeding season, individuals from the southern Karoo in southwestern South Africa may range up to 1000 km north to overwinter alongside breeding individuals in northern and central Namibia. This behavior could potentially result in the genetic mixing of populations if individuals from the southern Karoo breed again in central Namibia. This study addresses intraspecific mtDNA population structure and gene flow between populations of Namaqua sandgrouse in southern Africa. Our results, in general, contribute to the understanding of population processes in birds inhabiting arid zones of southern Africa and more broadly to the inference of demographic processes from molecular data.
 |
Methods |
|---|
 Top Methods ResultsDiscussionReferences |
|---|
Study Area and Sampling
The Namaqua sandgrouse, a near endemic to southern Africa, occurs throughout the arid western region of the subcontinent (see Figure 1). We sampled 10 populations during wing shooting hunts from 29 October 1994 to 23 June 2004. The 10 populations sampled (Table 1) were Benfontein (BF), Droëgrond (DG), Graaff Reinet (GR), Kgalagadi Transfrontier Park (K), Papkuil (PK), Springbok (SB), Groblershoop (GH), Soetdoring (SD), Tontelbos (TB), and Tswalu Kalahari Desert Reserve (TS). For the remainder of the article, we refer to these collecting localities as populations. Collecting covered the range of 3 morphologically described subspecies (Figure 1), yet preliminary molecular analysis indicated no corresponding genetically defined lineages (Delport W, unpublished data). Heart muscle was removed from a total of 81 Namaqua sandgrouse from these 10 localities. Tissue samples were stored in a dimethyl sulfoxide (DMSO)–saturated salt solution (5 M NaCl and 20% DMSO), at room temperature, with a ratio of 1 part tissue to 3 parts solution.
|
DNA Extraction, Polymerase Chain Reaction Amplification, and Cycle Sequencing
Genomic DNA was extracted from all samples using a standard phenol-chloroform DNA extraction procedure. Approximately 50 ng of genomic DNA was used as template in a polymerase chain reaction (PCR) (Saiki et al. 1988). Reactions were performed in a total volume of 50 µl containing, in addition to the genomic DNA, 2.5 mM MgCl2, 1x reaction buffer, 0.2 mM of each of 4 nucleotides, 1.5 U of Super-therm DNA polymerase (Southern Cross Biotechnology, Cape Town, South Africa) and 12.5 picamol of each primer. The primers L16075 (5' ACATTGGTCTTGTAAACC 3'—Delport et al. 2002) and ND6H (5' GCAGGTTGAGGTTTGTTGT 3'—Grant 1999) were used to amplify the complete mitochondrial control region of the 81 individuals. PCR conditions were as follows: denaturing at 94 °C for 2 min; followed by 35 cycles of denaturing at 94 °C for 30 s, primer annealing at 58 °C for 30 s, and elongation at 72 °C for 45 s; and finally an extended elongation period of 10 min at 72 °C. The PCR products were purified using a High Pure PCR Product Purification kit (Roche Diagnostics, Basel, Switzerland), according to the manufacturer's instructions. Purified PCR products were cycle sequenced using BigDye (Applied Biosystems, Foster City, CA), under conditions specified by the manufactured. Primers, H522 (5'GCC TGA CCG AGG AAC CAG AG 3'—modified from Quinn and Wilson 1993) and THR-EXT (5' CCC CTT CCC CCC CCG GCT 3'—Grant 1999), were used to yield a 318-bp sequence, from the 5' end of the control region, of both the heavy and light strands for each of the 81 individuals. Grant (1999) extracted purified mtDNA and sequenced the entire control region for several sandgrouse species. The sequences generated in this study from genomic DNA were comparable to those generated in the aforementioned study, and thus, we could be certain that we had not sequenced nuclear copies. All sequences used in this study were deposited in GenBank under the accession numbers AY181949–AY181975 and DQ025759–DQ025766.
Sequence Analysis
Sequences were proofread in Sequence Navigator (Applied Biosystems) and aligned using CLUSTAL X version 1.74 (Thompson et al. 1997). Analysis was conducted to gain an understanding of the demographic processes that could have produced the observed data. First, we tested for the presence of selection in the data (Tajima 1989; Fu and Li 1993). Outgroup-based tests of selection (D*/F* of Fu and Li) utilized chestnut-bellied sandgrouse (Pterocles exustus). Significance was evaluated from approximations to the underlying null distribution (Tajima's D) or from simulations of empirical distributions for theta in the interval [2, 20] (Fu and Li 1993). Second, we tested for the existence of population structure with 1) permutations of the chi-square test (Roff and Bentzen 1989) using custom software, 2) the calculation of pairwise FST's, using Arlequin (Schneider et al. 2000) and their deviation from the random association of alleles among localities (10 000 permutations), and 3) the construction of allele networks using TCS (Clement et al. 2000). This preliminary assessment of population structure was used to determine the next appropriate analysis.
Because little population structure was detected, we assumed a panmictic model and coestimated population growth and theta using FLUCTUATE (Kuhner et al. 1998). The Markov Chain Monte Carlo search strategy comprised 10 short chains and 2 long chains both with sampling increments of 20, yet of chain lengths 200 and 20 000, respectively. Replicate runs with alternate random seeds produced similar results, thus ensuring convergence on the correct parameter estimates. Because we assumed a panmictic population model, we used GENETREE (Bahlo and Griffiths 2000) to estimate symmetric migration rates between groups of populations. GENETREE was chosen over MIGRATE-n (Beerli and Felsenstein 2001) because the infinitely many alleles model of GENETREE allows for faster estimation times and less homoplasic inconsistencies, especially, in light of the high levels of migration between the specified groups. Typically, MIGRATE-n outperforms GENETREE in the estimation of migration when populations are structured, yet the inverse is true for populations that approach panmixia (Beerli P, personal communication). Populations were grouped by proximity for the GENETREE analysis (Figure 1). Because GENETREE calculates the likelihood surface in the vicinity of the true parameter estimates, it is necessary to provide upper and lower bounds for estimates on both theta and population growth. Here we used estimates of population growth and theta from FLUCTUATE to define the region in which to estimate the theta-migration-growth likelihood surface in GENETREE (theta: 0–20, g: 0–1000, 1 000 000 repetitions). This process of first estimating theta and growth parameters in FLUCTUATE, and subsequently, defining upper and lower bounds on the likelihood surface in GENETREE, allowed for faster estimation times and smoother likelihood surfaces from the GENETREE analysis.
|
Results |
|---|
|
TopMethodsResultsDiscussionReferences |
|---|
Mitochondrial control region sequences of 318 bp were obtained from the 81 individuals. Mean nucleotide frequencies (A = 31.8, C = 28.1, G = 14.8, and T = 25.3) were similar to that of other 5' avian control regions sequenced (Baker and Marshall 1997). Of the 318 positions sequenced, 22 were variable and defined 31 alleles, of which 20 were private alleles, only occurring in their respective populations (Table 1). Allelic diversity was high (0.926) relative to nucleotide diversity (0.008), an indication of either a selective sweep or rapid population growth (Slatkin and Hudson 1991). Neutrality tests indicated that the mitochondrial region was selectively neutral (Tajima's D = –1.35, D* of Fu and Li = –1.65, F* of Fu and Li = –1.84, D of Fu and Li = –1.02, F of Fu and Li = –1.36, and all P > 0.05). Preliminary tests of population structure (Roff and Bentzen 1989) indicated that the observed distribution of alleles among populations was not significantly different from random for populations (
2 = 22.81, P = 0.27) and for the Namaqualand/Karoo and Kalahari regions (2 = 37.3, P = 0.075). Little geographic structuring is evident in the allele network (Figure 2), yet many populations had unique alleles (Table 1). Pairwise FST's calculated in Arlequin again showed little evidence of geographic structure; however, some population pairs had FST values significantly different from zero and show significant evidence of population structure (Table 2).
|
|
We coestimated theta and population growth using FLUCTUATE. The maximum likelihood estimates of theta and population growth were 0.1484 and 587.6, respectively (Figure 3a). The larger estimate of theta under the maximum likelihood model (0.1484), compared with the unbiased estimate (average number of pairwise differences or nucleotide diversity 0.008), is consistent with longer times to coalescence given the observed population growth. The growth parameter in this case simply indicates that growth is positive. However, if inferences of the timing of growth are to be made, an estimate of the mutation rate is required. We calculated changes in population size, given the estimated growth parameter and upper and lower bounds on mtDNA mutation rates. Because much variability exists in mitochondrial control region mutation rate estimates (Zink and Blackwell 1998; Ruokonen and Kvist 2002; Zink and Weckstein 2003), we applied upper and lower bounds of 10–5 and 10–8, respectively, that encompass both protein-coding and control region rates of evolution. Given the bounds on mutation rate, the population has either undergone a recent rapid radiation (within the last 2500 generations) or a steady increase in population size over approximately 200 000 generations. Which of these hypotheses is most likely is difficult to infer without independent information on mutation rate.
|
Following the FLUCTUATE analysis, we used GENETREE to infer migration rates between groups of populations. GENETREE requires the data to be fully consistent with an infinitely many sites model (Bahlo and Griffiths 2000), which has the advantages of reducing both homoplasy and computation time. Eight alleles and 2 polymorphic sites not consistent with the infinitely many sites model were removed from the data set following suggestions in the GENETREE manual (Griffiths 2000). The coalescent-based tests of population structure using GENETREE showed high levels of gene flow between all populations (Figure 3b), where an overall level of symmetrical gene flow between all pairs of populations is represented. The GENETREE estimated theta (log theta = 4.5) is substantially larger than for that estimated using FLUCTUATE (log theta = –1.9), when assuming a panmictic population (Figure 3). Although the estimated migration rates between subpopulations in GENETREE are high and suggest panmixia, the specification of population structure in the GENETREE analysis clearly has an effect on the maximum likelihood estimate of theta (Figure 3). Theta is expected to be larger in subdivided populations because coalescent genealogies tend to be longer when populations are structured (Hudson 1990), and the alleles available for coalescence are fewer than what would be available given a panmictic population structure.
|
Discussion |
|---|
|
TopMethodsResultsDiscussionReferences |
|---|
Demographics and Population Structure
Analysis of population structure showed little evidence of the association of alleles with geographic localities. Neutrality tests indicate the absence of selection, which is furthermore supported by several studies that show mtDNA control region polymorphisms are usually neutral (Fry and Zink 1998; Milot et al. 2000; Griswold and Baker 2002). Although little geographic structure is evident, statistical permutations of FST indicated that at least some of the populations are significantly differentiated. This result may be in part due to high frequencies of alleles that are unique to particular populations. The significant differences detected with summary statistics have north to south (SD–PK) and southwest to northeast (PK–TS and SB–TS) orientations. Of these orientations, the north to south divergences could be the result of isolation by distance and the geographic isolation of the SD population. Differences between southwest (PK, SB) and northeast (TS) populations would be consistent with the partial migration accompanied by little mixing of flocks at the overwintering sites. Piertney et al. (2000) in their analysis of red grouse (Lagopus lagopus scoticus) mtDNA population structure similarly identified common alleles that were shared among the majority of populations, yet most other alleles occurred in at least 2 of the 12 populations. Similarly, Ball et al. (1988) found a wide geographic dispersion of alleles in red-winged blackbirds (Agelaius phoeniceus). Both these studies infer panmixia because it is unlikely that the same alleles have arisen in different populations by parallel mutations. A southern African migratory species, the red-billed quelea (Quelea quelea), also shows a lack of geographic population structure, in spite of an apparent migratory divide (Dallimer et al. 2003) in central southern Africa. The authors attribute the lack of genetic structure to the potential for interbreeding between populations on either side of the migratory divide during a second breeding attempt after returning to the interior of southern Africa.
It is likely that a high level of gene flow in Namaqua sandgrouse is causing the populations to approach panmixia. However, Namaqua sandgrouse differs in one respect from these former studies that inferred panmixia. The majority of Namaqua sandgrouse alleles (20 of 31) occur in only a single population. These 20 private alleles occur in 9 of the 10 populations, with BF as the only population without private alleles. This high degree of private alleles within populations should define the existence of population structure; however, the population growth signal overpowers its detection. The more recent an expansion event the more difficult it will to be detect population subdivision. Assuming bounds on mutation rate of 10–5 and 10–8, we have calculated that the expansion may have been rapid and fairly recent (2500 generations ago) or slow and old starting approximately 200 000 generations ago. Ray et al. (2003) have used simulations to make inferences of the process of range expansion and its effect on population genetic data. Specifically, they have shown that patterns of expansion observed in populations are not only the result of population size, and the age of the expansion event, but also depend on the level of gene flow between neighboring populations at the time of expansion. Range expansions with large levels of gene flow between populations tend to have most of their coalescent events occurring close to the time of expansion and, therefore, produce star-like genealogies. In contrast, range expansions with low levels of gene flow between populations have most coalescent events occurring before migration events. The crucial result of the former, where spatial expansions occurs rapidly due to high levels of gene flow between populations, is that recent demographic expansions with high levels of initial gene flow lead to patterns of molecular diversity that are equivalent to populations that have always been exchanging a large number of individuals with their neighbors (Ray et al. 2003). This would account for the extremely high levels of gene flow inferred using the maximum likelihood methods in this study, even though some populations have private alleles and significant differences detected with summary statistics.
The causes of such an expansion are uncertain, yet Grant and Leslie (1993) provide support for similar trends in mtDNA patterns from southern African arid zone taxa. Low levels of mtDNA diversity and high levels of inferred gene flow between demes are attributed to population cycles that are the result of extreme climatic variation in the region (Grant and Leslie 1993), where organellar diversity is lost more rapidly than nuclear DNA under severe population bottlenecks (Wilson et al. 1985). Grant and Leslie (1993), however, further demonstrate using simulation modeling that metapopulation structure combined with only low to moderate population bottlenecks can achieve similar patterns in loss of organellar diversity. The population structure of arid zone birds, such as Namaqua sandgrouse, combined with intermediate population size fluctuations, and a life history susceptible to climatic effects could account for the observed patterns. Namaqua sandgrouse has a nest that is a simple scrape in the ground (Lloyd et al. 2000) and is highly vulnerable to the elements. Furthermore, breeding success is highly variable with productivity estimates of 0.07–0.15 young per pair per year to 0.26–0.58 over 4 years in the Namaqualand, but 1.13–4.64 young per pair per year over a 14-month breeding season in the Kgalagadi National Park (Lloyd, Little, and Crowe 2001). Southern Africa is characterized by episodic droughts, which potentially could influence the reproductive success of Namaqua sandgrouse and cause population crashes. These factors most likely contribute to the observed patterns of mtDNA diversity in Namaqua sandgrouse. Higher resolution molecular markers have the potential to provide more power for the detection of population structure in Namaqua sandgrouse. In particular, multilocus approaches for the detection of cryptic population structure have been developed (Pritchard et al. 2000) and would certainly assist in the understanding of population movements within Namaqua sandgrouse.
|
Acknowledgments |
|---|
We thank Peter Bassingthwaighte, Rory Kroon, and Andrew Stainthorpe for collecting the samples from SD, GR, and TS, respectively. Furthermore, we would like to thank all the wing shooters who assisted in the collection of material for this research. The research reported in this article was supported by National Research Foundation grants to T.M.C (grant number: 2053582) and P.B (grant number: 2053653). The opinions and views presented in this article are, however, not necessarily those of the National Research Foundation. We thank 3 anonymous reviewers and subject editor, Robert C. Fleischer, for suggestions that allowed us to improve the manuscript.
|
Footnotes |
|---|
Corresponding Editor: Robert Fleischer
Received August 25, 2005
Accepted March 11, 2006
|
References |
|---|
|
TopMethodsResultsDiscussionReferences |
|---|
-
Bahlo M and Griffiths RC. (2000) Inference from gene trees in a subdivided population. Theor Popul Biol 57:79–95.[CrossRef][Web of Science][Medline]
Baker AJ and Marshall HD. (1997) Mitochondrial control region sequences as tools for understanding evolution. In Mindell DP (Ed.). Avian molecular evolution and systematics(Academic Press, San Diego (CA)) pp. 51–83.
Ball RM, Freeman S, James FC, Bermingham E, Avise JC. (1988) Phylogeographic population structure of redwinged blackbirds assessed by mitochondrial DNA. Proc Natl Acad Sci USA 85:1558–1562.
Beerli P and Felsenstein J. (2001) Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach. Proc Natl Acad Sci USA 98:4563–4568.
Clancey PA. (1979) On Pterocles namaqua (Gmelin) in South West Africa. Madoqua 2:261265.
Clement M, Posada D, Crandall KA. (2000) TCS: a computer program to estimate gene genealogies. Mol Ecol 9:1657–1659.[CrossRef][Medline]
Dallimer M, Jones PJ, Pemberton JM, Cheke RA. (2003) Lack of genetic and plumage differentiation in the red-billed quelea Quelea quelea across a migratory divide in southern Africa. Mol Ecol 12:345–353.[CrossRef][Medline]
Delport W, Ferguson JWH, Bloomer P. (2002) Characterization and evolution of the mitochondrial DNA control region in Hornbills (Bucerotiformes). J Mol Evol 54:794–806.[CrossRef][Web of Science][Medline]
Emerson BC, Paradis E, Thébaud C. (2001) Revealing the demographic histories of species using DNA sequences. Trends Ecol Evol 16:707–716.[CrossRef]
Fry AJ and Zink RM. (1998) Geographic analysis of nucleotide diversity and song sparrow (Aves: Emberizidae) population history. Mol Ecol 7:1303–1313.[CrossRef][Medline]
Fu Y-X and Li W-H. (1993) Statistical tests of neutrality of mutations. Genetics 133:693709.
Grant TJ. (1999) Characterisation of the Namaqua Sandgrouse mitochondrial DNA control region [MSc dissertation]. (University of Pretoria, Pretoria (South Africa)).
Grant WS and Leslie RW. (1993) Effect of metapopulation structure on nuclear and organellar DNA variability in semi-arid environments of southern Africa. S Afr J Sci 89:287–293.[Web of Science]
Griffiths RC. (2000) Genetree version 9.0 documentation [Internet]. [cited 2005 Aug 1]. Available from: http://www.stats.ox.ac.uk/
griff/software.html.
Griswold CK and Baker AJ. (2002) Time to most recent common ancestor and divergence times of populations of common chaffinches (Fringilla coelebs) in Europe and North Africa: insights into Pleistocene refugia and current levels of migration. Evolution 56:143153.
Hudson RR. (1990) Gene genealogies and the coalescent process. In Futuyama D and Antonovics J (Eds.). Oxford surveys in evolutionary biology(Oxford University Press, Oxford) Vol. 7:.
Knowles LL and Maddison WP. (2002) Statistical phylogeography. Mol Ecol 11:2623–2635.[CrossRef][Medline]
Knowles LL. (2004) The burgeoning field of statistical phylogeography. J Evol Biol 17:1–10.[CrossRef][Web of Science][Medline]
Kuhner MK, Yamato J, Felsenstein J. (1998) Maximum likelihood estimation of population growth rates based on the coalescent. Genetics 149:429–434.
Lloyd P, Little RM, Crowe TM. (2000) Daily activity pattern of a Namaqua Sandgrouse, Pterocles namaqua, population. Ostrich 71:425437.
Lloyd P, Little RM, Crowe TM. (2001) The breeding biology of the Namaqua Sandgrouse, Pterocles namaqua. Ostrich 72:169178.
Lloyd P, Little RM, Crowe TM, Simmons RE. (2001) Rainfall and food availability as factors influencing the migration and breeding activity of Namaqua Sandgrouse, Pterocles namaqua.. Ostrich 72:5062.
Malan G, Little RM, Crowe TM. (1994) Temporal and spatial patterns of abundance and breeding activity of Namaqua Sandgrouse in South Africa. S Afr J Zool 29:162–167.
Milot E, Gibbs LH, Hobson KA. (2000) Phylogeography and genetic structure of northern populations of the Yellow Warbler (Oendroica petechia). Mol Ecol 9:667–681.[CrossRef][Medline]
Nielsen R and Wakeley J. (2001) Distinguishing migration from isolation: a Markov Chain Monte Carlo approach. Genetics 158:885–896.
Piertney SB, MacColl ADC, Bacon PJ, Racey PA, Lambin X, Dallas JF. (2000) Matrilineal genetic structure and female-mediated gene flow in Red Grouse (Lagopus lagopus scoticus): an analysis using mitochondrial DNA. Evolution 54:279–289.[CrossRef][Web of Science][Medline]
Pritchard JK, Stephens M, Donnelly P. (2000) Inference of population structure using multilocus genotype data. Genetics 55:945959.
Quinn TW and Wilson AC. (1993) Sequence evolution in and around the mitochondrial control region in birds. J Mol Evol 37:417–425.[CrossRef][Web of Science][Medline]
Ray N, Currat M, Excoffier L. (2003) Intra-deme molecular diversity in spatially expanding expanding populations. Mol Biol Evol 20:76–86.
Roff DA and Bentzen P. (1989) The statistical analysis of mitochondrial DNA polymorphisms: chi-square and the problem of small samples. Mol Biol Evol 6:539–545.[Abstract]
Ruokonen M and Kvist L. (2002) Structure and evolution of the avian mitochondrial control region. Mol Phylogenet Evol 23:422–432.[CrossRef][Web of Science][Medline]
Saiki RK, Gelfand DH, Stoffel S, Scharf SJ, Higuchi R, Horn GT, Mullis KB, Erlich HA. (1988) Primer-directed enzymatic amplification of DNA with thermostable DNA polymerase. Science 239:487–491.
Schneider S, Roessli D, Excoffier L. (2000) Arlequin version. 2.000: a software for population genetics data analysis. (Genetics Biometry Laboratory, University of Geneva, Geneva (Switzerland)).
Slatkin M and Hudson RR. (1991) Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics 129:555–562.[Abstract]
Tajima F. (1989) Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123:585–595.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. (1997) The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 24:4876–4882.
Wakeley J and Hey J. (1998) Testing speciation models with DNA sequence data. In DeSalle R and Schierwater B (Eds.). Molecular approaches to ecology(Birkhäuser Verlag, Basel (Switzerland)) pp. 157–175.
Ward P. (1972) The functional significance of mass drinking flights by sandgrouse: Pteroclidae. Ibis 114:533536.
Wilson AC, Cann RL, Carr SM, George M, Gyllensten UB, Helm-Bychowski KM, Higuchi RG, Palumbi SR, Prager EM, Sage RD, et al. (1985) Mitochondrial DNA and two perspectives on evolutionary genetics. Biol J Linn Soc 26:375–400.[Web of Science]
Zink RM and Blackwell RC. (1998) Molecular systematics and biogeography of aridland gnatcatchers (genus Polioptila) and evidence supporting species status of the California Gnatcatcher (Polioptila californien). Mol Phylogenet Evol 9:26–32.[CrossRef][Web of Science][Medline]
Zink RM and Weckstein JD. (2003) Recent evolutionary history of the Fox Sparrows (genus: Passerella). Auk 12:522–527.[CrossRef]
CiteULike 
Connotea 
Del.icio.us What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||