Genetic Population History Research Articles

Sum of parts is greater than the whole: inference of common genetic history of populations

BackgroundReconstructability of population history, from genetic information of extant individuals, is studied under a simulation setting. We do not address the issue of accuracy of the reconstruction algorithms: we assume the availability of the theoretical best algorithm. On the other hand, we focus on the fraction (1 - f) of the common genetic history that is irreconstructible or impenetrable. Thus the fraction, f, gives an upper bound on the extent of estimability. In other words, there exists no method that can reconstruct a fraction larger than f of the entire common genetic history. For the realization of such a study, we first define a natural measure of the amount of genetic history. Next, we use a population simulator (from literature) that has at least two features. Firstly, it has the capability of providing samples from different demographies, to effectively reflect reality. Secondly, it also provides the underlying relevant genetic history, captured in its entirety, where such a measure is applicable. Finally, to compute f, we use an information content measure of the relevant genetic history. The simulator of choice provided the following demographies: Africans, Europeans, Asians and Afro-Americans.ResultsWe observe that higher the rate of recombination, lower the value of f, while f is invariant over varying mutation rates, in each of the demographies. The value of f increases with the number of samples, reaching a plateau and suggesting that in all the demographies at least about one-third of the relevant genetic history is impenetrable. The most surprising observation is that the the sum of the reconstructible history of the subsegments is indeed larger than the reconstructible history of the whole segment. In particular, longer the chromosomal segment, smaller the value of f, in all the demographies.ConclusionsWe present the very first framework for measuring the fraction of the relevant genetic history of a population that is mathematically elusive. Our observed results on the tested demographies suggest that it may be better to aggregate the analysis of smaller chunks of chromosomal segments than fewer large chunks. Also, no matter the richness of samples in a population, at least one-third of the population genetic history is impenetrable. The framework also opens up possible new lines of investigation along the following. Given the characteristics of a population, possibly derived from observed extant individuals, to estimate the (1) optimal sample size and (2) optimal sequence length for the most informative analysis.

Isolation and characterization of microsatellite DNA loci from Naihai cutlassfish (Trichiurus nanhaiensis)

The Naihai cutlassfish (Trichiurus nanhaiensis) is an important commercial fish species. In this study, we report 25 microsatellite DNA markers for T. nanhaiensis. The number of alleles per locus ranged from 2 to 10 in a sample of 37 individuals. Observed and expected heterozygosities per locus varied from 0.300 to 1.000 and from 0.403 to 0.845, respectively. Thirteen of 25 loci showed significant departure from Hardy–Weinberg equilibrium (HWE) after sequential Bonferroni correction (P < 0.002). No significant linkage disequilibrium between pairs of loci was found. These microsatellite markers will provide powerful tools for investigating population genetic structure, population history and conservation management of T. nanhaiensis. The family Trichiuridae, including 32 species in nine genera, is generally benthopelagic on continental shelves and slopes from tropical to warm temperate waters (Nakamura and Parin 1993). Three species of genus Trichiurus were identified in South China Sea based on morphological characteristics and results of isozyme electrophoresis (Wang et al. 1992, 1993). T. nanhaiensis ranks second only behind T. haumela in tonnage of fish caught and both are important commercial species in South China Sea. So far, the fisheries biology, growth or mortality parameters and genetic structure of T. haumela have been studied, but few studies on T. nanhaiensis have been reported (Zhang et al. 1998; Chakraborty et al. 2005; Zhu and Qiu 2005; Hsu et al. 2007; Yang et al. 2007). In the present study, we developed 25 microsatellite loci isolated from a dinucleotide-enriched genomic library of T. nanhaiensis using the fast isolation by AFLP of sequences containing repeats (FIASCO) according to previous report (Zane et al. 2002). These markers should provide a tool for

Microsatellite markers developed for Corallodiscus lanuginosus (Gesneriaceae) and their cross‐species transferability

Corallodiscus is widespread across China and in adjacent countries to the west and south, with C. lanuginosus covering the entire range of the genus. Microsatellite markers will be useful to address within-complex taxonomic and biogeographic structures of this species. • Fourteen markers were developed using the Fast Isolation by AFLP of Sequences COntaining repeats (FIASCO) protocol. Polymorphism was assessed in 12 individuals each from two populations from West Yunnan, China. Nine primers displayed polymorphisms. The number of alleles per locus ranged from one to six (mean: 2.7). The observed and expected heterozygosities ranged from 0 to 1 and from 0 to 0.772 (average: 0.487 and 0.448). Only one of the polymorphic loci deviated from Hardy-Weinberg equilibrium. Among the four congeneric species, cross-amplification success ranged from two to eight markers. • The markers will be useful for population genetic and evolutionary history studies for C. lanuginosus and the allied species C. bhutanicus, C. cooperi, and C. kingianus.

PCAdmix: Principal Components-Based Assignment of Ancestry Along Each Chromosome in Individuals with Admixed Ancestry from Two or More Populations

Identifying ancestry along each chromosome in admixed individuals provides a wealth of information for understanding the population genetic history of admixture events and is valuable for admixture mapping and identifying recent targets of selection. We present PCAdmix (available at https://sites.google.com/site/pcadmix/home ), a Principal Components-based algorithm for determining ancestry along each chromosome from a high-density, genome-wide set of phased single-nucleotide polymorphism (SNP) genotypes of admixed individuals. We compare our method to HAPMIX on simulated data from two ancestral populations, and we find high concordance between the methods. Our method also has better accuracy than LAMP when applied to three-population admixture, a situation as yet unaddressed by HAPMIX. Finally, we apply our method to a data set of four Latino populations with European, African, and Native American ancestry. We find evidence of assortative mating in each of the four populations, and we identify regions of shared ancestry that may be recent targets of selection and could serve as candidate regions for admixture-based association mapping.

PCAdmix: Principal Components-Based Assignment of Ancestry along Each Chromosome in Individuals with Admixed Ancestry from Two or More Populations

Identifying ancestry along each chromosome in admixed individuals provides a wealth of information for understanding the population genetic history of admixture events and is valuable for admixture mapping and identifying recent targets of selection. We present PCAdmix (available at https://sites.google.com/site/pcadmix/home), a Principal Components-based algorithm for determining ancestry along each chromosome from a high-density, genome-wide set of phased single-nucleotide polymorphism (SNP) genotypes of admixed individuals. We compare our method to HAPMIX on simulated data from two ancestral populations, and we find high concordance between the methods. Our method also has better accuracy than LAMP when applied to three-population admixture, a situation as yet unaddressed by HAPMIX. Finally, we apply our method to a data set of four Latino populations with European, African, and Native American ancestry. We find evidence of assortative mating in each of the four populations, and we identify regions of shared ancestry that may be recent targets of selection and could serve as candidate regions for admixture-based association mapping.

Population genetic structure and demographic history of small yellow croaker, Larimichthys polyactis (Bleeker, 1877), from coastal waters of China

Population genetic structure and demographic history of small yellow croaker, Larimichthys polyactis (Bleeker, 1877), from coastal waters of China

The relative efficacy of functional and developmental cranial modules for reconstructing global human population history

This study tests the relative efficacy of human cranial modules, defined on the basis of developmental and functional criteria, for reconstructing neutral genetic population history. Specifically, two hypotheses were tested: 1) The "basicranial hypothesis" predicts that the endochondrally ossifying basicranium will be more reliable for reconstructing population history than intramembranously ossifying regions of the human cranium. This is based on the assumption that early ossification of the basicranium and its distinct functional constraints produce a cranial structure that is relatively immune to non-neutral evolutionary forces. 2) The "single function hypothesis" predicts that cranial regions associated with a single (sensory) function are less reliable indicators of neutral genetic history. Here the prediction is based on the logic that complex, multi-functional, integrated cranial regions are less likely toexhibit homoplasy and, therefore, provide a more accurate morphological proxy for genetic relationships. The congruence between craniometric affinity matrices and neutral genetic population matrices based on autosomal microsatellite and classical markers was assessed using a series of Mantel and Dow-Cheverud tests. The results did not support the predictions of the "basicranial hypothesis," as the endochondrally ossifying basicranium was not significantly more congruent with the genetic data than intramembraneously ossifying modules. Moreover, although the results provided some support for the "single function hypothesis," defining cranial modules on the basis of anatomical or functional complexity did not provide a consistent means of predicting their phylogenetic efficacy. These results have important implications for building an accurate inference model of cranial evolution in the human fossil record.

Pedigree analysis and inbreeding depression on growth traits in Brazilian Marchigiana and Bonsmara breeds1

The study of population structure by pedigree analysis is useful to identify important circumstances that affect the genetic history of populations. The intensive use of a small number of superior individuals may reduce the genetic diversity of populations. This situation is very common for the beef cattle breeds. Therefore, the objectives of the present study were to analyze the pedigree and possible inbreeding depression on traits of economic interest in the Marchigiana and Bonsmara breeds and to test the inclusion of the individual inbreeding coefficient (F(i)) or individual increases in inbreeding coefficient (ΔF(i)) in the genetic evaluation model for the quantification of inbreeding depression. The complete pedigree file of the Marchigiana breed included 29,411 animals born between 1950 and 2003. For the Bonsmara breed, the pedigree file included 18,695 animals born between 1988 and 2006. Only animals with at least 2 equivalent generations of known pedigree were kept in the analyses of inbreeding effect on birth weight, weaning weight measured at about 205 d, and BW at 14 mo in the Marchigiana breed, and on birth weight, weaning weight, and scrotal circumference measured at 12 mo in the Bonsmara breed. The degree of pedigree knowledge was greater for Marchigiana than for Bonsmara animals. The average generation interval was 7.02 and 3.19 for the Marchigiana and Bonsmara breed, respectively. The average inbreeding coefficient was 1.33% for Marchigiana and 0.26% for Bonsmara. The number of ancestors explaining 50% of the gene pool and effective population size computed via individual increase in coancestry were 13 and 97.79 for Marchigiana and 41 and 54.57 for Bonsmara, respectively. These estimates indicate reduction in genetic variability in both breeds. Inbreeding depression was observed for most of the growth traits. The model including ΔF(i) can be considered more adequate to quantify inbreeding depression. The inclusion of F(i) or ΔF(i) in the genetic evaluation model may not result in better fit to the data. A genetic evaluation with simultaneous estimation of inbreeding depression can be performed in Marchigiana and Bonsmara breeds, providing additional information to producers and breeders.

Evidence from Cameroon reveals differences in the genetic structure and histories of chimpanzee populations

The history of the genus Pan is a topic of enduring interest. Chimpanzees (Pan troglodytes) are often divided into subspecies, but the population structure and genetic history of chimpanzees across Africa remain unclear. Some population genetics studies have led to speculation that, until recently, this species constituted a single population with ongoing gene flow across its range, which resulted in a continuous gradient of allele frequencies. Chimpanzees, designated here as P. t. ellioti, occupy the Gulf of Guinea region that spans southern Nigeria and western Cameroon at the center of the distribution of this species. Remarkably, few studies have included individuals from this region, hindering the examination of chimpanzee population structure across Africa. Here, we analyzed microsatellite genotypes of 94 chimpanzees, including 32 designated as P. t. ellioti. We find that chimpanzees fall into three major populations: (i) Upper Guinea in western Africa (P. t. verus); (ii) the Gulf of Guinea region (P. t. ellioti); and (iii) equatorial Africa (P. t. troglodytes and P. t. schweinfurthii). Importantly, the Gulf of Guinea population is significantly different genetically from the others, sharing a last common ancestor with the populations in Upper Guinea ~0.46 million years ago (mya) and equatorial Africa ~0.32 mya. Equatorial chimpanzees are subdivided into up to three populations occupying southern Cameroon, central Africa, and eastern Africa, which may have constituted a single population until ~0.10-0.11 mya. Finally, occasional hybridization may be occurring between the Gulf of Guinea and southern Cameroon populations.

Population genetic structure and history of fragmented remnant populations of the New England cottontail (Sylvilagus transitionalis)

The New England cottontail (Sylvilagus transitionalis) has suffered from extensive loss and fragmentation of its habitat and is now a species of conservation priority in the northeastern United States. Remnant New England cottontail populations currently occur in five geographically disjunct locations: southern Maine and southeastern New Hampshire (MENH); the Merrimack Valley in central New Hampshire (NH-MV); Cape Cod, Massachusetts (CC); parts of eastern Connecticut and Rhode Island (CTRI); and western Connecticut, southeastern New York and southwestern Massachusetts (CTNY). We used microsatellite genotyping to discern patterns of population structure, genetic variability, and demographic history across the species’ range and to assess whether the observed patterns are a consequence of recent habitat loss and fragmentation. Our findings show that the geographic populations are highly differentiated (overall FST = 0.145; P < 0.001). Using Bayesian clustering analyses, we identified five genetic clusters, which corresponded closely to the geographic populations, but grouped MENH & NH-MV together (ME/NH) and identified an isolated population in eastern Connecticut (Bluff Point). The genetic clusters showed little evidence of recent gene flow and are highly influenced by genetic drift. The CC and Bluff Point populations show signs they experienced a genetic bottleneck, whereas the ME/NH population shows evidence of ongoing decline. Populations in Bluff Point, CC, and ME/NH also show significantly reduced genetic variation relative to the other clusters (CTNY and CTRI without Bluff Point). Without immediate human intervention, the short-term persistence of New England cottontail populations in Maine, New Hampshire and Cape Cod is at great risk. Conservation efforts at this time should focus on within-population sustainability and eventually restoring connectivity among these isolated populations.

Population genetic structure and evolutionary history of grass carp Ctenopharyngodon idella in the Yangtze River, China

Ample studies have been conducted to investigate the population genetic structure of grass carp Ctenopharyngodon idella in the Yangtze River, China. However, samples from the upper reaches were not included. In this study, we collected samples from the entire river, including three locations in the upper reaches: Yibin, Banan and Yunyang, two locations in the middle reaches: Shishou and Ruichang and one location in the lower reaches: Hanjiang, and sequenced three mitochondrial coding genes (ND5, ND6 and Cytb) and one control region (i.e., the D-loop). Nineteen haplotypes were observed in grass carp of the Yangtze River through the analysis of combined sequence data sets (around 4428 bp). Haplotype diversity indices (0.6000 ∼ 0.9333) and nucleotide diversity indices (0.0002 ∼ 0.0020) demonstrated low genetic diversity in the Yangtze grass carp. The analysis of molecular variance and the fixation index (F ST = 0.0202) revealed insignificant genetic difference between samples from different reaches. Two monophyletic lineages of haplotypes were identified, with the lineage A experiencing potential expansion events. Along with previous findings, this study provides a better understanding of genetic diversity and variation of grass carp in the Yangtze River and will be served as an important baseline to evaluate the long-term impact of the Three Gorges Dam and other hydroelectric facilities on fish biodiversity.

Population genetic structure and history of a generalist parasite infecting multiple sympatric host species

Host specificity is predicted to shape patterns of parasite gene flow between host species; specialist parasites should have low gene flow between host species, while generalists are predicted to have high gene flow between species. However, even for generalist parasites external forces, including ecological differences between host species may sometimes intervene to limit gene flow and create genetic structure. To investigate the potential for cryptic parasite genetic structure to arise under such circumstances, we examined the population genetic structure and history of the generalist nematode, Trichostrongylus axei, infecting six sympatric wild ungulate species in North America. Using genotypes for 186 T. axei larvae at two mitochondrial genes, cox1 and nad4, we found that T. axei was completely panmictic across host species, with 0% of genetic variation structured between host species and 97% within individual hosts. In addition, T. axei showed no evidence of recent genetic bottlenecks, had high nucleotide diversities (above 2%), and an effective population size estimated to be in the tens of millions. Our result that T. axei maintains high rates of gene flow between multiple sympatric host species adds to a growing body of information on trichostrongylid population genetic structure in different ecological contexts. Furthermore, the high rates of gene flow, coupled with high levels of genetic diversity and large effective population size which we observed in T. axei, point to a potentially broad capacity for rapid evolutionary change in this parasite.

Identification of population substructure among Jews using STR markers and dependence on reference populations included.

BackgroundDetecting population substructure is a critical issue for association studies of health behaviors and other traits. Whether inherent in the population or an artifact of marker choice, determining aspects of a population's genetic history as potential sources of substructure can aid in design of future genetic studies. Jewish populations, among which association studies are often conducted, have a known history of migrations. As a necessary step in understanding population structure to conduct valid association studies of health behaviors among Israeli Jews, we investigated genetic signatures of this history and quantified substructure to facilitate future investigations of these phenotypes in this population.ResultsUsing 32 autosomal STR markers and the program STRUCTURE, we differentiated between Ashkenazi (AJ, N = 135) and non-Ashkenazi (NAJ, N = 226) Jewish populations in the form of Northern and Southern geographic genetic components (AJ north 73%, south 23%, NAJ north 33%, south 60%). The ability to detect substructure within these closely related populations using a small STR panel was contingent on including additional samples representing major continental populations in the analyses.ConclusionsAlthough clustering programs such as STRUCTURE are designed to assign proportions of ancestry to individuals without reference population information, when Jewish samples were analyzed in the absence of proxy parental populations, substructure within Jews was not detected. Generally, for samples with a given grandparental country of birth, STRUCTURE assignment values to Northern, Southern, African and Asian clusters agreed with mitochondrial DNA and Y-chromosomal data from previous studies as well as historical records of migration and intermarriage.

Population genetic history of the dreissenid mussel invasions: expansion patterns across North America

This study tests population genetic patterns across the Eurasian dreissenid mussel invasions of North America—encompassing the zebra mussel Dreissena polymorpha (1986 detection) and the quagga mussel D. rostriformis bugensis (detected in 1990, which now has largely displaced the former in the Great Lakes). We evaluate their source-spread relationships and invasion genetics using 9–11 nuclear microsatellite loci for 583 zebra mussels (21 sites) and 269 quagga mussels (12 sites) from Eurasian and North American range locations, with the latter including the Great Lakes, Mississippi River basin, Atlantic coastal waterways, Colorado River system, and California reservoirs. Additionally, mtDNA cytochrome b gene sequences are used to verify species identity. Our results indicate that North American zebra mussels originate from multiple non-native northern European populations, whereas North American quagga mussels trace to native estuaries in the Southern Bug and Dnieper Rivers. Invasive populations of both species show considerable genetic diversity and structure (zebra F ST = 0.006–0.263, quagga F ST = 0.008–0.267), without founder effects. Most newer zebra mussel populations have appreciable genetic diversity, whereas quagga mussel populations from the Colorado River and California show some founder effects. The population genetic composition of both species changed over time at given sites; with some adding alleles from adjacent populations, some losing them, and all retaining closest similarity to their original composition. Zebra mussels from Kansas and California appear genetically similar and assign to a possible origin from the St. Lawrence River, whereas quagga mussels from Nevada and California assign to a possible origin from Lake Ontario. These assignments suggest that overland colonization pathways via recreational boats do not necessarily reflect the most proximate connections. In conclusion, our microsatellite results comprise a valuable baseline for resolving present and future dreissenid mussel invasion pathways.

Do taxonomic divisions reflect genetic differentiation? A comparison of morphological and genetic data in Coenonympha tullia (Müller), Satyrinae

Historically, morphological variation has been used to classify many species (and subspecies) of Lepidoptera. However, some of this variation may be unsuitable for inferring the recent evolutionary history of populations. Genetic data provide an alternative. We examine the morphological and genetic variation within and between British subspecies of Coenonympha tullia (Muller 1764) to test the hypothesis that neutral genetic variation corresponds to morphological variation. We find that most morphological and genetic variation occurs within populations and that those populations designated as subspecies based on morphological characters are not necessarily most closely related for mitochondrial DNA or nuclear DNA (amplified fragment length polymorphisms and allozymes). Thus, the notion that wing spot variation reflects population isolation and therefore genetic differentiation does not hold. The present study highlights the need for genetic data where taxonomy may be based on environmentally plastic or locally adapted characters because such characters will not reflect the true population genetic history.

The value of some Corsican sub-populations for genetic association studies

BackgroundGenetic isolates with a history of a small founder population, long-lasting isolation and population bottlenecks represent exceptional resources in the identification of disease genes. In these populations the disease allele reveals Linkage Disequilibrium (LD) with markers over significant genetic intervals, therefore facilitating disease locus identification. In a previous study we examined the LD extension on the Xq13 region in three Corsican sub-populations from the inner mountainous region of the island. On the basis of those previous results we have proposed a multistep procedure to carry out studies aimed at the identification of genes involved in complex diseases in Corsica. A prerequisite to carry out the proposed multi-step procedure was the presence of different degrees of LD on the island and a common genetic derivation of the different Corsican sub-populations. In order to evaluate the existence of these conditions in the present paper we extended the analysis to the Corsican coastal populations.MethodsSamples were analyzed using seven dinucleotide microsatellite markers on chromosome Xq13-21: DXS983, DXS986, DXS8092, DXS8082, DXS1225, DXS8037 and DXS995 spanning approximately 4.0 cM (13.3 Mb). We have also investigated the distribution of the DXS1225-DXS8082 haplotype which has been recently proposed as a good marker of population genetic history due to its low recombination rate.Resultsthe results obtained indicate a decrease of LD on the island from the central mountainous toward the coastal sub-populations. In addition the analysis of the DXS1225-DXS8082 haplotype revealed: 1) the presence of a particular haplotype with high frequency; 2) the derivation from a common genetic pool of the sub-populations examined in the present study.ConclusionThese results indicate the Corsican sub-populations useful for the fine mapping of genes contributing to complex diseases.

Complex population genetic and demographic history of the Salangid, Neosalanx taihuensis, based on cytochrome b sequences.

BackgroundThe Salangid icefish Neosalanx taihuensis (Salangidae) is an economically important fish, which is endemic to China, restricted to large freshwater systems (e.g. lakes, large rivers and estuaries) and typically exhibit low vagility. The continuous distribution ranges from the temperate region of the Huai and Yellow River basins to the subtropical region of the Pearl River basin. This wide ranging distribution makes the species an ideal model for the study of palaeoclimatic effects on population genetic structure and phylogeography. Here, we aim to analyze population genetic differentiation within and between river basins and demographic history in order to understand how this species responded to severe climatic oscillations, decline of the sea levels during the Pleistocene ice ages and tectonic activity.ResultsWe obtained the complete mtDNA cytochrome b sequences (1141 bp) of 354 individuals from 13 populations in the Pearl River, the Yangze River and the Huai River basin. Thirty-six haplotypes were detected. Haplotype frequency distributions were strongly skewed, with most haplotypes (n = 24) represented only in single samples each and thus restricted to a single population. The most common haplotype (H36) was found in 49.15% of all individuals. Analysis of molecular variance (AMOVA) revealed a random pattern in the distribution of genetic diversity, which is inconsistent with contemporary hydrological structure. Significant levels of genetic subdivision were detected among populations within basins rather than between the three basins. Demographic analysis revealed that the population size in the Pearl River basin has remained relatively constant whereas the populations in the Yangze River and the Huai River basins expanded about 221 and 190 kyr ago, respectively, with the majority of mutations occurring after the last glacial maximum (LGM).ConclusionThe observed complex genetic pattern of N. taihuensis is coherent with a scenario of multiple unrelated founding events by long-distance colonization and dispersal combined with contiguous population expansion and locally restricted gene flow. We also found that this species was likely severely impacted by past glaciations. More favourable climate and the formation of large suitable habitations together facilitated population expansion after the late Quaternary (especially the LGM). We proposed that all populations should be managed and conserved separately, especially for habitat protection.

Fixation Probabilities When Generation Times Are Variable: The Burst–Death Model

Estimating the fixation probability of a beneficial mutation has a rich history in theoretical population genetics. Typically, to attain mathematical tractability, we assume that generation times are fixed, while the number of offspring per individual is stochastic. However, fixation probabilities are extremely sensitive to these assumptions regarding life history. In this article, we compute the fixation probability for a "burst-death" life-history model. The model assumes that generation times are exponentially distributed, but the number of offspring per individual is constant. We estimate the fixation probability for populations of constant size and for populations that grow exponentially between periodic population bottlenecks. We find that the fixation probability is, in general, substantially lower in the burst-death model than in classical models. We also note striking qualitative differences between the fates of beneficial mutations that increase burst size and mutations that increase the burst rate. In particular, once the burst size is sufficiently large relative to the wild type, the burst-death model predicts that fixation probability depends only on burst rate.

Sampling genetic diversity in the sympatrically and allopatrically speciating Midas cichlid species complex over a 16 year time series

BackgroundSpeciation often occurs in complex or uncertain temporal and spatial contexts. Processes such as reinforcement, allopatric divergence, and assortative mating can proceed at different rates and with different strengths as populations diverge. The Central American Midas cichlid fish species complex is an important case study for understanding the processes of speciation. Previous analyses have demonstrated that allopatric processes led to species formation among the lakes of Nicaragua as well as sympatric speciation that is occurring within at least one crater lake. However, since speciation is an ongoing process and sampling genetic diversity of such lineages can be biased by collection scheme or random factors, it is important to evaluate the robustness of conclusions drawn on individual time samples.ResultsIn order to assess the validity and reliability of inferences based on different genetic samples, we have analyzed fish from several lakes in Nicaragua sampled at three different times over 16 years. In addition, this time series allows us to analyze the population genetic changes that have occurred between lakes, where allopatric speciation has operated, as well as between different species within lakes, some of which have originated by sympatric speciation. Focusing on commonly used genetic markers, we have analyzed both DNA sequences from the complete mitochondrial control region as well as nuclear DNA variation at ten microsatellite loci from these populations, sampled thrice in a 16 year time period, to develop a robust estimate of the population genetic history of these diversifying lineages.ConclusionThe conclusions from previous work are well supported by our comprehensive analysis. In particular, we find that the genetic diversity of derived crater lake populations is lower than that of the source population regardless of when and how each population was sampled. Furthermore, changes in various estimates of genetic diversity within lakes are minimal and provide no evidence for drastic changes during the last 20 years, supporting the hypothesis that the processes which have resulted in rapid speciation are primarily historical. In contrast, there is some evidence for ongoing evolution, particularly selection, in all lakes except crater Lake Masaya, perhaps reflecting the persistence of speciational processes. Importantly, we find that the crater Lake Apoyo population, for which strong evidence of sympatric speciation has been demonstrated, has lower genetic diversity than other crater lakes and the strongest evidence for ongoing selection.

Divergent Haplotypes and Human History as Revealed in a Worldwide Survey of X-Linked DNA Sequence Variation

The population genetic history of a 10.1-kbp noncoding region of the human X chromosome was studied using the males of the HGDP-CEPH Human Genome Diversity Panel (672 individuals from 52 populations). The geographic distribution of patterns of variation was roughly consistent with previous studies, with the major exception that 1 highly divergent haplotype (haplotype X, hX) was observed at low frequency in widely scattered non-African populations and not at all observed in sub-Saharan African populations. Microsatellite (short tandem repeat) variation within the sequenced region was low among copies of hX, even though the estimated time of ancestry of hX and other sequences was 1.44 Myr. The estimated age of the common ancestor of all hX copies was 5,230 years (95% consistency index: 2,000-75,480 years). To further address the presence of hX in Africa, additional samples from Chad and Tanzania were screened. Five additional copies of hX were observed, consistent with a history in which hX was present in Africa prior to the migration of modern humans out of Africa and with eastern Africa being the source of non-African modern human populations. Taken together, these features of hX-that it is much older than other haplotypes and uncommon and patchily distributed throughout Africa, Europe, and Asia-present a cautionary tale for interpretations of human history.