The pedigree rate of sequence divergence in the human mitochondrial genome: there is a difference between phylogenetic and pedigree rates

Whole genome identity-by-descent determination

High-throughput single nucleotide polymorphism genotyping assays conveniently produce genotype data for genome-wide genetic linkage and association studies. For pedigree datasets, the unphased genotype data is used to infer the haplotypes for individuals, according to Mendelian inheritance rules. Linkage studies can then locate putative chromosomal regions based on the haplotype allele sharing among the pedigree members and their disease status. Most existing haplotyping programs require rather strict pedigree structures and return a single inferred solution for downstream analysis. In this research, we relax the pedigree structure to contain ungenotyped founders and present a cubic time whole genome haplotyping algorithm to minimize the number of zero-recombination haplotype blocks. With or without explicitly enumerating all the haplotyping solutions, the algorithm determines all distinct haplotype allele identity-by-descent (IBD) sharings among the pedigree members, in linear time in the total number of haplotyping solutions. Our algorithm is implemented as a computer program iBDD. Extensive simulation experiments using 2 sets of 16 pedigree structures from previous studies showed that, in general, there are trillions of haplotyping solutions, but only up to a few thousand distinct haplotype allele IBD sharings. iBDD is able to return all these sharings for downstream genome-wide linkage and association studies.

The history of Slavs inferred from complete mitochondrial genome sequences

To shed more light on the processes leading to crystallization of a Slavic identity, we investigated variability of complete mitochondrial genomes belonging to haplogroups H5 and H6 (63 mtDNA genomes) from the populations of Eastern and Western Slavs, including new samples of Poles, Ukrainians and Czechs presented here. Molecular dating implies formation of H5 approximately 11.5–16 thousand years ago (kya) in the areas of southern Europe. Within ancient haplogroup H6, dated at around 15–28 kya, there is a subhaplogroup H6c, which probably survived the last glaciation in Europe and has undergone expansion only 3–4 kya, together with the ancestors of some European groups, including the Slavs, because H6c has been detected in Czechs, Poles and Slovaks. Detailed analysis of complete mtDNAs allowed us to identify a number of lineages that seem specific for Central and Eastern Europe (H5a1f, H5a2, H5a1r, H5a1s, H5b4, H5e1a, H5u1, some subbranches of H5a1a and H6a1a9). Some of them could possibly be traced back to at least ∼4 kya, which indicates that some of the ancestors of today's Slavs (Poles, Czechs, Slovaks, Ukrainians and Russians) inhabited areas of Central and Eastern Europe much earlier than it was estimated on the basis of archaeological and historical data. We also sequenced entire mitochondrial genomes of several non-European lineages (A, C, D, G, L) found in contemporary populations of Poland and Ukraine. The analysis of these haplogroups confirms the presence of Siberian (C5c1, A8a1) and Ashkenazi-specific (L2a1l2a) mtDNA lineages in Slavic populations. Moreover, we were able to pinpoint some lineages which could possibly reflect the relatively recent contacts of Slavs with nomadic Altaic peoples (C4a1a, G2a, D5a2a1a1).

Inferring recent historic abundance from current genetic diversity

Recent historic abundance is an elusive parameter of great importance for conserving endangered species and understanding the pre-anthropogenic state of the biosphere. The number of studies that have used population genetic theory to estimate recent historic abundance from contemporary levels of genetic diversity has grown rapidly over the last two decades. Such assessments often yield unexpectedly large estimates of historic abundance. We review the underlying theory and common practices of estimating recent historic abundance from contemporary genetic diversity, and critically evaluate the potential issues at various estimation steps. A general issue of mismatched spatio-temporal scales between the estimation itself and the objective of the estimation emerged from our assessment; genetic diversity-based estimates of recent historic abundance represent long-term averages, whereas the objective typically is an estimate of recent abundance for a specific population. Currently, the most promising approach to estimate the difference between recent historic and contemporary abundance requires that genetic data be collected from samples of similar spatial and temporal duration. Novel genome-enabled inference methods may be able to utilize additional information of dense genome-wide distributions of markers, such as of identity-by-descent tracts, to infer recent historic abundance from contemporary samples only.

Purifying selection in mammalian mitochondrial protein-coding genes is highly effective and congruent with evolution of nuclear genes

The mammalian mitochondrial genomes differ from the nuclear genomes by maternal inheritance, absence of recombination, and higher mutation rate. All these differences decrease the effective population size of mitochondrial genome and make it more susceptible to accumulation of slightly deleterious mutations. It was hypothesized that mitochondrial genes, especially in species with low effective population size, irreversibly degrade leading to decrease of organismal fitness and even to extinction of species through the mutational meltdown. To interrogate this hypothesis, we compared the purifying selections acting on the representative set of mitochondrial (potentially degrading) and nuclear (potentially not degrading) protein-coding genes in species with different effective population size. For 21 mammalian species, we calculated the ratios of accumulation of slightly deleterious mutations approximated by Kn/Ks separately for mitochondrial and nuclear genomes. The 75% of variation in Kn/Ks is explained by two independent variables: type of a genome (mitochondrial or nuclear) and effective population size of species approximated by generation time. First, we observed that purifying selection is more effective in mitochondria than in the nucleus that implies strong evolutionary constraints of mitochondrial genome. Mitochondrial de novo nonsynonymous mutations have at least 5-fold more harmful effect when compared with nuclear. Second, Kn/Ks of mitochondrial and nuclear genomes is positively correlated with generation time of species, indicating relaxation of purifying selection with decrease of species-specific effective population size. Most importantly, the linear regression lines of mitochondrial and nuclear Kn/Ks's from generation times of species are parallel, indicating congruent relaxation of purifying selection in both genomes. Thus, our results reveal that the distribution of selection coefficients of de novo nonsynonymous mitochondrial mutations has a similar shape with the distribution of de novo nonsynonymous nuclear mutations, but its mean is five times smaller. The harmful effect of mitochondrial de novo nonsynonymous mutations triggers highly effective purifying selection, which maintains the fitness of the mammalian mitochondrial genome.

An alternative model for the early peopling of southern South America revealed by analyses of three mitochondrial DNA haplogroups

After several years of research, there is now a consensus that America was populated from Asia through Beringia, probably at the end of the Pleistocene. But many details such as the timing, route(s), and origin of the first settlers remain uncertain. In the last decade genetic evidence has taken on a major role in elucidating the peopling of the Americas. To study the early peopling of South America, we sequenced the control region of mitochondrial DNA from 300 individuals belonging to indigenous populations of Chile and Argentina, and also obtained seven complete mitochondrial DNA sequences. We identified two novel mtDNA monophyletic clades, preliminarily designated B2l and C1b13, which together with the recently described D1g sub-haplogroup have locally high frequencies and are basically restricted to populations from the extreme south of South America. The estimated ages of D1g and B2l, about ~15,000 years BP, together with their similar population dynamics and the high haplotype diversity shown by the networks, suggests that they probably appeared soon after the arrival of the first settlers and agrees with the dating of the earliest archaeological sites in South America (Monte Verde, Chile, 14,500 BP). One further sub-haplogroup, D4h3a5, appears to be restricted to Fuegian-Patagonian populations and reinforces our hypothesis of the continuity of the current Patagonian populations with the initial founders. Our results indicate that the extant native populations inhabiting South Chile and Argentina are a group which had a common origin, and suggest a population break between the extreme south of South America and the more northern part of the continent. Thus the early colonization process was not just an expansion from north to south, but also included movements across the Andes.

Phylogenetic analysis of four nuclear protein-encoding genes largely corroborates the traditional classification of Bivalvia (Mollusca)

Revived interest in molluscan phylogeny has resulted in a torrent of molecular sequence data from phylogenetic, mitogenomic, and phylogenomic studies. Despite recent progress, basal relationships of the class Bivalvia remain contentious, owing to conflicting morphological and molecular hypotheses. Marked incongruity of phylogenetic signal in datasets heavily represented by nuclear ribosomal genes versus mitochondrial genes has also impeded consensus on the type of molecular data best suited for investigating bivalve relationships. To arbitrate conflicting phylogenetic hypotheses, we evaluated the utility of four nuclear protein-encoding genes-ATP synthase β, elongation factor-1α, myosin heavy chain type II, and RNA polymerase II-for resolving the basal relationships of Bivalvia. We sampled all five major lineages of bivalves (Archiheterodonta, Euheterodonta [including Anomalodesmata], Palaeoheterodonta, Protobranchia, and Pteriomorphia) and inferred relationships using maximum likelihood and Bayesian approaches. To investigate the robustness of the phylogenetic signal embedded in the data, we implemented additional datasets wherein length variability and/or third codon positions were eliminated. Results obtained include (a) the clade (Nuculanida+Opponobranchia), i.e., the traditionally defined Protobranchia; (b) the monophyly of Pteriomorphia; (c) the clade (Archiheterodonta+Palaeoheterodonta); (d) the monophyly of the traditionally defined Euheterodonta (including Anomalodesmata); and (e) the monophyly of Heteroconchia, i.e., (Palaeoheterodonta+Archiheterodonta+Euheterodonta). The stability of the basal tree topology to dataset manipulation is indicative of signal robustness in these four genes. The inferred tree topology corresponds closely to those obtained by datasets dominated by nuclear ribosomal genes (18S rRNA and 28S rRNA), controverting recent taxonomic actions based solely upon mitochondrial gene phylogenies.

The influence of rate heterogeneity among sites on the time dependence of molecular rates

Molecular evolutionary rate estimates have been shown to depend on the time period over which they are estimated. Factors such as demographic processes, calibration errors, purifying selection, and the heterogeneity of substitution rates among sites (RHAS) are known to affect the accuracy with which rates of evolution are estimated. We use mathematical modeling and Bayesian analyses of simulated sequence alignments to explore how mutational hotspots can lead to time-dependent rate estimates. Mathematical modeling shows that underestimation of molecular rates over increasing time scales is inevitable when RHAS is ignored. Although a gamma distribution is commonly used to model RHAS, we show that when the actual RHAS deviates from a gamma-like distribution, rates can either be under- or overestimated in a time-dependent manner. Simulations performed under different scenarios of RHAS confirm the mathematical modeling and demonstrate the impacts of time-dependent rates on estimates of divergence times. Most notably, erroneous rate estimates can have narrow credibility intervals, leading to false confidence in biased estimates of rates, and node ages. Surprisingly, large errors in estimates of overall molecular rate do not necessarily generate large errors in divergence time estimates. Finally, we illustrate the correlation between time-dependent rate patterns and differential saturation between quickly and slowly evolving sites. Our results suggest that data partitioning or simple nonparametric mixture models of RHAS significantly improve the accuracy with which node ages and substitution rates can be estimated.

High mitochondrial mutation rates estimated from deep-rooting Costa Rican pedigrees

Estimates of mutation rates for the noncoding hypervariable Region I (HVR-I) of mitochondrial DNA vary widely, depending on whether they are inferred from phylogenies (assuming that molecular evolution is clock-like) or directly from pedigrees. All pedigree-based studies so far were conducted on populations of European origin. In this article, we analyzed 19 deep-rooting pedigrees in a population of mixed origin in Costa Rica. We calculated two estimates of the HVR-I mutation rate, one considering all apparent mutations, and one disregarding changes at sites known to be mutational hot spots and eliminating genealogy branches which might be suspected to include errors, or unrecognized adoptions along the female lines. At the end of this procedure, we still observed a mutation rate equal to 1.24 × 10(-6) , per site per year, i.e., at least threefold as high as estimates derived from phylogenies. Our results confirm that mutation rates observed in pedigrees are much higher than estimated assuming a neutral model of long-term HVRI evolution. We argue that until the cause of these discrepancies will be fully understood, both lower estimates (i.e., those derived from phylogenetic comparisons) and higher, direct estimates such as those obtained in this study, should be considered when modeling evolutionary and demographic processes.

Statistical guidelines for detecting past population shifts using ancient DNA

Populations carry a genetic signal of their demographic past, providing an opportunity for investigating the processes that shaped their evolution. Our ability to infer population histories can be enhanced by including ancient DNA data. Using serial-coalescent simulations and a range of both quantitative and temporal sampling schemes, we test the power of ancient mitochondrial sequences and nuclear single-nucleotide polymorphisms (SNPs) to detect past population bottlenecks. Within our simulated framework, mitochondrial sequences have only limited power to detect subtle bottlenecks and/or fast post-bottleneck recoveries. In contrast, nuclear SNPs can detect bottlenecks followed by rapid recovery, although bottlenecks involving reduction of less than half the population are generally detected with low power unless extensive genetic information from ancient individuals is available. Our results provide useful guidelines for scaling sampling schemes and for optimizing our ability to infer past population dynamics. In addition, our results suggest that many ancient DNA studies may face power issues in detecting moderate demographic collapses and/or highly dynamic demographic shifts when based solely on mitochondrial information.

Postglacial recolonization of eastern Blacknose Dace,Rhinichthys atratulus(Teleostei: Cyprinidae), through the gateway of New England

During the last ice age, much of North America far south as 40°N was covered by glaciers (Hewitt 2000). About 20,000 years ago, as the glaciers retreated, the hydrologic landscape changed dramatically creating waterways for fish dispersal. The number of populations responsible for recolonization and the regions from which they recolonized are unknown for many freshwater fishes living in New England and southeastern Canada. The Blacknose Dace,Rhinichthys atratulus, is one of the freshwater fish species that recolonized this region. We hypothesize that the earliest deglaciated region, modern-day Connecticut, was recolonized byR. atratulusvia a single founding event by a single population. In this paper, we test this hypothesis phylogenetically with regard to the major drainage basins within Connecticut. The mitochondrial DNA exhibits low nucleotide diversity, high haplotype diversity, and a dominant haplotype found across the state. A small percentage of individuals in the Housatonic drainage basin, however, share a haplotype with populations in New York drainage basins, a haplotype not found elsewhere in Connecticut's drainage basins. We calculated a range for the rate of divergence for NADH dehydrogenase subunit 2 (nd2) and control region (ctr) of 4.43-6.76% and 3.84-8.48% per million years (my), respectively. While this range is higher than the commonly accepted rate of 2% for mitochondrial DNA, these results join a growing list of publications finding high rates of divergence for various taxa (Peterson and Masel 2009). The data support the conclusion that Connecticut as a whole was recolonized initially by a single founding event that came from a single refugium. Subsequently, the Housatonic basin alone experienced a secondary recolonization event.

Genetic differences among Vietnamese Haplorchis taichui populations using the COI genetic marker

Adults of the fish-borne intestinal trematode species Haplorchis taichui were collected from humans in three provinces of Vietnam: Ha Giang, Thanh Hoa and Quang Tri. Genetic analysis revealed three groups of the parasite from clustering dendrograms, correlating with the localities in which they were collected. Measurements of evolutionary divergence over sequence pairs were greater between the different populations than within them, which indicated that the three populations were genetically different. The significance (Fst= 0.73; P value < 0.05) of the genetic variation of the three studied populations implied that genetic separation of the populations had already occurred, which may have been caused by a low gene flow among the different H. taichui populations. Factors contributing to the low gene flow may include isolation resulting from the intermediate-host fish rarely being sold outside of the rural commune where they are raised and the enclosed aquacultural areas themselves.

Reconciling deep calibration and demographic history: bayesian inference of post glacial colonization patterns in Carcinus aestuarii (Nardo, 1847) and C. maenas (Linnaeus, 1758)

A precise inference of past demographic histories including dating of demographic events using Bayesian methods can only be achieved with the use of appropriate molecular rates and evolutionary models. Using a set of 596 mitochondrial cytochrome c oxidase I (COI) sequences of two sister species of European green crabs of the genus Carcinus (C. maenas and C. aestuarii), our study shows how chronologies of past evolutionary events change significantly with the application of revised molecular rates that incorporate biogeographic events for calibration and appropriate demographic priors. A clear signal of demographic expansion was found for both species, dated between 10,000 and 20,000 years ago, which places the expansions events in a time frame following the Last Glacial Maximum (LGM). In the case of C. aestuarii, a population expansion was only inferred for the Adriatic-Ionian, suggestive of a colonization event following the flooding of the Adriatic Sea (18,000 years ago). For C. maenas, the demographic expansion inferred for the continental populations of West and North Europe might result from a northward recolonization from a southern refugium when the ice sheet retreated after the LGM. Collectively, our results highlight the importance of using adequate calibrations and demographic priors in order to avoid considerable overestimates of evolutionary time scales.

Phylogeographic patterns of decapod crustaceans at the Atlantic-Mediterranean transition

Comparative multispecies studies allow contrasting the effect of past and present oceanographic processes on phylogeographic patterns. In the present study, a fragment of the COI gene was analyzed in seven decapod crustacean species from five families and with different bathymetric distributions. A total of 769 individuals were sampled along the Atlantic-Mediterranean transition area in order to test the effect of three putative barriers to gene flow: Strait of Gibraltar, Almeria-Oran Front and Ibiza Channel. A significant effect of the Strait of Gibraltar was found in the crabs Liocarcinus depurator and Macropipus tuberculatus. The Ibiza Channel had a significant effect for L. depurator. However, the Almeria-Oran front was not found to have a significant effect on any of the studied species. Higher levels of population structure were found in shallow-water species, although the number of species sampled should be increased to obtain a conclusive pattern. The haplotypes within the different species coalesced at times that could be related with past climatic events occurring before, during and after the last glacial maximum. Given the large diversity of phylogeographic patterns obtained within decapods, it is concluded that both historical and contemporary processes (marine current patterns, bathymetry and life-history traits) shape the phylogeographic patterns of these crustaceans.

Faunal histories from Holocene ancient DNA

Recent studies using ancient DNA have been instrumental in advancing understanding of the impact of Holocene climate change on biodiversity. Ancient DNA has been used to track demography, migration and diversity, and is providing new insights into the long-term dynamics of species and population distributions. The Holocene is key to understanding how the past has impacted on the present, as it bridges the gap between contemporary phylogeographic studies and those with inference on Pleistocene patterns, based on ancient DNA studies. Here, we examine the major patterns of Holocene faunal population dynamics and connectivity; highlighting the dynamic nature of species and population responses to Holocene climatic change, thereby providing an 'analogue' for understanding potential impacts of future change.

Time-dependent rates of molecular evolution

For over half a century, it has been known that the rate of morphological evolution appears to vary with the time frame of measurement. Rates of microevolutionary change, measured between successive generations, were found to be far higher than rates of macroevolutionary change inferred from the fossil record. More recently, it has been suggested that rates of molecular evolution are also time dependent, with the estimated rate depending on the timescale of measurement. This followed surprising observations that estimates of mutation rates, obtained in studies of pedigrees and laboratory mutation-accumulation lines, exceeded long-term substitution rates by an order of magnitude or more. Although a range of studies have provided evidence for such a pattern, the hypothesis remains relatively contentious. Furthermore, there is ongoing discussion about the factors that can cause molecular rate estimates to be dependent on time. Here we present an overview of our current understanding of time-dependent rates. We provide a summary of the evidence for time-dependent rates in animals, bacteria and viruses. We review the various biological and methodological factors that can cause rates to be time dependent, including the effects of natural selection, calibration errors, model misspecification and other artefacts. We also describe the challenges in calibrating estimates of molecular rates, particularly on the intermediate timescales that are critical for an accurate characterization of time-dependent rates. This has important consequences for the use of molecular-clock methods to estimate timescales of recent evolutionary events.

Dynamics of mitochondrial heteroplasmy in three families investigated via a repeatable re-sequencing study

Background

Originally believed to be a rare phenomenon, heteroplasmy - the presence of more than one mitochondrial DNA (mtDNA) variant within a cell, tissue, or individual - is emerging as an important component of eukaryotic genetic diversity. Heteroplasmies can be used as genetic markers in applications ranging from forensics to cancer diagnostics. Yet the frequency of heteroplasmic alleles may vary from generation to generation due to the bottleneck occurring during oogenesis. Therefore, to understand the alterations in allele frequencies at heteroplasmic sites, it is of critical importance to investigate the dynamics of maternal mtDNA transmission.

Results

Here we sequenced, at high coverage, mtDNA from blood and buccal tissues of nine individuals from three families with a total of six maternal transmission events. Using simulations and re-sequencing of clonal DNA, we devised a set of criteria for detecting polymorphic sites in heterogeneous genetic samples that is resistant to the noise originating from massively parallel sequencing technologies. Application of these criteria to nine human mtDNA samples revealed four heteroplasmic sites.

Conclusions

Our results suggest that the incidence of heteroplasmy may be lower than estimated in some other recent re-sequencing studies, and that mtDNA allelic frequencies differ significantly both between tissues of the same individual and between a mother and her offspring. We designed our study in such a way that the complete analysis described here can be repeated by anyone either at our site or directly on the Amazon Cloud. Our computational pipeline can be easily modified to accommodate other applications, such as viral re-sequencing.

Purifying selection of mtDNA and its implications for understanding evolution and mitochondrial disease

Mutations of mitochondrial DNA (mtDNA) are frequent in humans and are implicated in many different types of pathology. The high substitution rate and the maternal, asexual mode of transmission of mtDNA make it more likely to accumulate deleterious mutations. Here, we discuss recent evidence that mtDNA transmission is subject to strong purifying selection in the mammalian female germ line, limiting the accumulation of such mutations. This process shapes mitochondrial sequence diversity and is therefore probably of fundamental importance for animal evolution and in human mitochondrial disease.

Post-glacial partitioning of mitochondrial genetic variation in the field vole

Genetic markers are often used to examine population history. There is considerable debate about the behaviour of molecular clock rates around the population-species transition. Nevertheless, appropriate calibration is critical to any inference regarding the absolute timing and scale of demographic changes. Here, we use a mitochondrial cytochrome b gene genealogy, based entirely on modern sequences and calibrated from recent geophysical events, to date the post-glacial expansion of the Eurasian field vole (Microtus agrestis), a widespread temperate mammal species. The phylogeographic structure reflects the subsequent expansion of populations that went through bottlenecks at the time of the Younger Dryas (ca 12,000 years BP) rather than the Last Glacial Maximum (LGM, ca 24,000 years BP), which is usually seen as the time when present-day patterns were determined. The nucleotide substitution rate that was estimated here, ca 4 × 10(-7) substitutions/site/year, remains extremely high throughout the relevant time frame. Calibration with similarly high population-based substitution rates, rather than long-term rates derived from species divergence times, will show that post-LGM climatic events generated current phylogeographic structure in many other organisms from temperate latitudes.

Mitochondrial DNA reveals distinct evolutionary histories for Jewish populations in Yemen and Ethiopia

Southern Arabia and the Horn of Africa are important geographic centers for the study of human population history because a great deal of migration has characterized these regions since the first emergence of humans out of Africa. Analysis of Jewish groups provides a unique opportunity to investigate more recent population histories in this area. Mitochondrial DNA is used to investigate the maternal evolutionary history and can be combined with historical and linguistic data to test various population histories. In this study, we assay mitochondrial control region DNA sequence and diagnostic coding variants in Yemenite (n = 45) and Ethiopian (n = 41) Jewish populations, as well as in neighboring non-Jewish Yemeni (n = 50) and Ethiopian (previously published Semitic speakers) populations. We investigate their population histories through a comparison of haplogroup distributions and phylogenetic networks. A high frequency of sub-Saharan African L haplogroups was found in both Jewish populations, indicating a significant African maternal contribution unlike other Jewish Diaspora populations. However, no identical haplotypes were shared between the Yemenite and Ethiopian Jewish populations, suggesting very little gene flow between the populations and potentially distinct maternal population histories. These new data are also used to investigate alternate population histories in the context of historical and linguistic data. Specifically, Yemenite Jewish mitochondrial diversity reflects potential descent from ancient Israeli exiles and shared African and Middle Eastern ancestry with little evidence for large-scale conversion of local Yemeni. In contrast, the Ethiopian Jewish population appears to be a subset of the larger Ethiopian population suggesting descent primarily through conversion of local women.

Mitochondrial pathogenic mutations are population-specific

Background

Surveying deleterious variation in human populations is crucial for our understanding, diagnosis and potential treatment of human genetic pathologies. A number of recent genome-wide analyses focused on the prevalence of segregating deleterious alleles in the nuclear genome. However, such studies have not been conducted for the mitochondrial genome.

Results

We present a systematic survey of polymorphisms in the human mitochondrial genome, including those predicted to be deleterious and those that correspond to known pathogenic mutations. Analyzing 4458 completely sequenced mitochondrial genomes we characterize the genetic diversity of different types of single nucleotide polymorphisms (SNPs) in African (L haplotypes) and non-African (M and N haplotypes) populations. We find that the overall level of polymorphism is higher in the mitochondrial compared to the nuclear genome, although the mitochondrial genome appears to be under stronger selection as indicated by proportionally fewer nonsynonymous than synonymous substitutions. The African mitochondrial genomes show higher heterozygosity, a greater number of polymorphic sites and higher frequencies of polymorphisms for synonymous, benign and damaging polymorphism than non-African genomes. However, African genomes carry significantly fewer SNPs that have been previously characterized as pathogenic compared to non-African genomes.

Conclusions

Finding SNPs classified as pathogenic to be the only category of polymorphisms that are more abundant in non-African genomes is best explained by a systematic ascertainment bias that favours the discovery of pathogenic polymorphisms segregating in non-African populations. This further suggests that, contrary to the common disease-common variant hypothesis, pathogenic mutations are largely population-specific and different SNPs may be associated with the same disease in different populations. Therefore, to obtain a comprehensive picture of the deleterious variability in the human population, as well as to improve the diagnostics of individuals carrying African mitochondrial haplotypes, it is necessary to survey different populations independently.

Reviewers

This article was reviewed by Dr Mikhail Gelfand, Dr Vasily Ramensky (nominated by Dr Eugene Koonin) and Dr David Rand (nominated by Dr Laurence Hurst).

Recovering population parameters from a single gene genealogy: an unbiased estimator of the growth rate

We show that the number of lineages ancestral to a sample, as a function of time back into the past, which we call the number of lineages as a function of time (NLFT), is a nearly deterministic property of large-sample gene genealogies. We obtain analytic expressions for the NLFT for both constant-sized and exponentially growing populations. The low level of stochastic variation associated with the NLFT of a large sample suggests using the NLFT to make estimates of population parameters. Based on this, we develop a new computational method of inferring the size and growth rate of a population from a large sample of DNA sequences at a single locus. We apply our method first to a sample of 1,212 mitochondrial DNA (mtDNA) sequences from China, confirming a pattern of recent population growth previously identified using other techniques, but with much smaller confidence intervals for past population sizes due to the low variation of the NLFT. We further analyze a set of 63 mtDNA sequences from blue whales (BWs), concluding that the population grew in the past. This calls for reevaluation of previous studies that were based on the assumption that the BW population was fixed.

Pedigree likelihood ratio for lineage markers

Lineage-based haplotype markers (e.g., Y chromosome STRs and mitochondrial DNA sequences) are important adjunct tools to the autosomal markers for kinship analysis and for specialized kinship applications such as database searching. Traditionally, the prosecution or kinship hypothesis considers the haplotypes in the same lineage and the probability of genotype data given the lineage hypothesis is simply set at 1 if the number of mismatched loci or nucleotides between the questioned person and the references is less than a predefined threshold. In this study, a kinship hypothesis based on a fixed relationship of the questioned person in the reference family is introduced. A graphical model is proposed to calculate the probability of the genotype data given the kinship hypothesis, which is the product of haplotype frequency of the founder in the pedigree and the transmission probability from the founder to all descendants. Proper mutation models are suggested for Y chromosome STRs and mitochondrial DNA sequence variants (i.e., SNPs) to calculate the transmission probability. The methods to infer the genotypes of the untyped individuals in the pedigree and the computational complexity of handling these untyped individuals are also addressed. Lastly, numerical examples of the applications are given to demonstrate the kinship hypothesis and the algorithms.

Spatial-temporal stratifications in natural populations and how they affect understanding and estimation of effective population size

The concept of effective population size (N(e) ) is based on an elegantly simple idea which, however, rapidly becomes very complex when applied to most real-world situations. In natural populations, spatial and temporal stratifications create different classes of individuals with different vital rates, and this in turn affects (generally reduces) N(e) in complex ways. I consider how these natural stratifications influence our understanding of effective size and how to estimate it, and what the consequences are for conservation and management of natural populations. Important points that emerge include the following: 1  The relative influences of local vs metapopulation N(e) depend on a variety of factors, including the time frame of interest. 2  Levels of diversity in local populations are strongly influenced by even low levels of migration, so these measures are not reliable indicators of local N(e) . 3  For long-term effective size, obtaining a reliable estimate of mutation rate is the most important consideration; unless this is accomplished, estimates can be biased by orders of magnitude. 4  At least some estimators of contemporary N(e) appear to be robust to relatively high (approximately 10%) equilibrium levels of migration, so under many realistic scenarios they might yield reliable estimates of local N(e) . 5  Age structure probably has little effect on long-term estimators of N(e) but can strongly influence contemporary estimates. 6  More research is needed in several key areas: (i) to disentangle effects of selection and drift in metapopulations connected by intermediate levels of migration; (ii) to elucidate the relationship between N(b) (effective number of breeders per year) and N(e) per generation in age-structured populations; (iii) to perform rigorous sensitivity analyses of new likelihood and coalescent-based methods for estimating demographic and evolutionary histories.

Genetic and historic evidence for climate-driven population fragmentation in a top cetacean predator: the harbour porpoises in European water

Recent climate change has triggered profound reorganization in northeast Atlantic ecosystems, with substantial impact on the distribution of marine assemblages from plankton to fishes. However, assessing the repercussions on apex marine predators remains a challenging issue, especially for pelagic species. In this study, we use Bayesian coalescent modelling of microsatellite variation to track the population demographic history of one of the smallest temperate cetaceans, the harbour porpoise (Phocoena phocoena) in European waters. Combining genetic inferences with palaeo-oceanographic and historical records provides strong evidence that populations of harbour porpoises have responded markedly to the recent climate-driven reorganization in the eastern North Atlantic food web. This response includes the isolation of porpoises in Iberian waters from those further north only approximately 300 years ago with a predominant northward migration, contemporaneous with the warming trend underway since the ‘Little Ice Age’ period and with the ongoing retreat of cold-water fishes from the Bay of Biscay. The extinction or exodus of harbour porpoises from the Mediterranean Sea (leaving an isolated relict population in the Black Sea) has lacked a coherent explanation. The present results suggest that the fragmentation of harbour distribution range in the Mediterranean Sea was triggered during the warm ‘Mid-Holocene Optimum’ period (approx. 5000 years ago), by the end of the post-glacial nutrient-rich ‘Sapropel’ conditions that prevailed before that time.

Measurements of spontaneous rates of mutations in the recent past and the near future

The rate of spontaneous mutation in natural populations is a fundamental parameter for many evolutionary phenomena. Because the rate of mutation is generally low, most of what is currently known about mutation has been obtained through indirect, complex and imprecise methodological approaches. However, in the past few years genome-wide sequencing of closely related individuals has made it possible to estimate the rates of mutation directly at the level of the DNA, avoiding most of the problems associated with using indirect methods. Here, we review the methods used in the past with an emphasis on next generation sequencing, which may soon make the accurate measurement of spontaneous mutation rates a matter of routine.

Evaluating the Farming/Language Dispersal Hypothesis with genetic variation exhibited by populations in the Southwest and Mesoamerica

The Farming/Language Dispersal Hypothesis posits that prehistoric population expansions, precipitated by the innovation or early adoption of agriculture, played an important role in the uneven distribution of language families recorded across the world. In this case, the most widely spread language families today came to be distributed at the expense of those that have more restricted distributions. In the Americas, Uto-Aztecan is one such language family that may have been spread across Mesoamerica and the American Southwest by ancient farmers. We evaluated this hypothesis with a large-scale study of mitochondrial DNA (mtDNA) and Y-chromosomal DNA variation in indigenous populations from these regions. Partial correlation coefficients, determined with Mantel tests, show that Y-chromosome variation in indigenous populations from the American Southwest and Mesoamerica correlates significantly with linguistic distances (r = 0.33-0.384; P < 0.02), whereas mtDNA diversity correlates significantly with only geographic distance (r = 0.619; P = 0.002). The lack of correlation between mtDNA and Y-chromosome diversity is consistent with differing population histories of males and females in these regions. Although unlikely, if groups of Uto-Aztecan speakers were responsible for the northward spread of agriculture and their languages from Mesoamerica to the Southwest, this migration was possibly biased to males. However, a recent in situ population expansion within the American Southwest (2,105 years before present; 99.5% confidence interval = 1,273-3,773 YBP), one that probably followed the introduction and intensification of maize agriculture in the region, may have blurred ancient mtDNA patterns, which might otherwise have revealed a closer genetic relationship between females in the Southwest and Mesoamerica.

Compensatory evolution in mitochondrial tRNAs navigates valleys of low fitness

A long-standing controversy in evolutionary biology is whether or not evolving lineages can cross valleys on the fitness landscape that correspond to low-fitness genotypes, which can eventually enable them to reach isolated fitness peaks. Here we study the fitness landscapes traversed by switches between different AU and GC Watson-Crick nucleotide pairs at complementary sites of mitochondrial transfer RNA stem regions in 83 mammalian species. We find that such Watson-Crick switches occur 30-40 times more slowly than pairs of neutral substitutions, and that alleles corresponding to GU and AC non-Watson-Crick intermediate states segregate within human populations at low frequencies, similar to those of non-synonymous alleles. Substitutions leading to a Watson-Crick switch are strongly correlated, especially in mitochondrial tRNAs encoded on the GT-nucleotide-rich strand of the mitochondrial genome. Using these data we estimate that a typical Watson-Crick switch involves crossing a fitness valley of a depth of about 10(-3) or even about 10(-2), with AC intermediates being slightly more deleterious than GU intermediates. This compensatory evolution must proceed through rare intermediate variants that never reach fixation. The ubiquitous nature of compensatory evolution in mammalian mitochondrial tRNAs and other molecules implies that simultaneous fixation of two alleles that are individually deleterious may be a common phenomenon at the molecular level.

Explaining the imperfection of the molecular clock of hominid mitochondria

The molecular clock of mitochondrial DNA has been extensively used to date various genetic events. However, its substitution rate among humans appears to be higher than rates inferred from human-chimpanzee comparisons, limiting the potential of interspecies clock calibrations for intraspecific dating. It is not well understood how and why the substitution rate accelerates. We have analyzed a phylogenetic tree of 3057 publicly available human mitochondrial DNA coding region sequences for changes in the ratios of mutations belonging to different functional classes. The proportion of non-synonymous and RNA genes substitutions has reduced over hundreds of thousands of years. The highest mutation ratios corresponding to fast acceleration in the apparent substitution rate of the coding sequence have occurred after the end of the Last Ice Age. We recalibrate the molecular clock of human mtDNA as 7990 years per synonymous mutation over the mitochondrial genome. However, the distribution of substitutions at synonymous sites in human data significantly departs from a model assuming a single rate parameter and implies at least 3 different subclasses of sites. Neutral model with 3 synonymous substitution rates can explain most, if not all, of the apparent molecular clock difference between the intra- and interspecies levels. Our findings imply the sluggishness of purifying selection in removing the slightly deleterious mutations from the human as well as the Neandertal and chimpanzee populations. However, for humans, the weakness of purifying selection has been further exacerbated by the population expansions associated with the out-of Africa migration and the end of the Last Ice Age.

Coalescent simulations of Yakut mtDNA variation suggest small founding population

The Yakuts are a Turkic-speaking population from northeastern Siberia who are believed to have originated from ancient Turkic populations in South Siberia, based on archaeological and ethnohistorical evidence. In order to better understand Yakut origins, we modeled 25 demographic scenarios and tested by coalescent simulation whether any are consistent with the patterns of mtDNA diversity observed in present-day Yakuts. The models consist of either two simulated demes that represent Yakuts and a South Siberian ancestral population, or three demes that also include a regional Northeast Siberian population that served as a source of local gene flow into the Yakut deme. The model that produced the best fit to the observed data defined a founder group with an effective female population size of only 150 individuals that migrated northwards approximately 1,000 years BP and who experienced significant admixture with neighboring populations in Northeastern Siberia. These simulation results indicate a pronounced founder effect that was primarily kin-structured and reconcile reported discrepancies between Yakut mtDNA and Y chromosome diversity levels.

Quantitative prediction of molecular clock and ka/ks at short timescales

Recent empirical studies of taxa including humans, fish, and birds have shown elevated rates of molecular evolution between species that diverged recently. Using the Moran model, we calculate expected divergence as a function of time. Our findings suggest that the observed phenomenon of elevated rates at short timescales is consistent with standard population genetics theory. The apparent acceleration of the molecular clock at short timescales can be explained by segregating polymorphisms present at the time of the ancestral population, both neutral and slightly deleterious, and not newly arising slightly deleterious mutations as has been previously hypothesized. Our work also suggests that the duration of the rate elevation depends on the effective population size, providing a method to correct time estimates of recent divergence events. Our model concords with estimates of divergence obtained from African cichlid fish and humans. As an additional application of our model, we calculate that K(a)/K(s) is elevated within a population before decaying slowly to its long-term value. Similar to the molecular clock, the duration and magnitude of K(a)/K(s) elevation depend on the effective population size. Unlike the molecular clock, however, K(a)/K(s) elevation is caused by newly arising slightly deleterious mutations. This elevation, although not as severe in magnitude as had been previously predicted in models neglecting ancestral polymorphism, persists slightly longer.

Rapid response of a marine mammal species to holocene climate and habitat change

Environmental change drives demographic and evolutionary processes that determine diversity within and among species. Tracking these processes during periods of change reveals mechanisms for the establishment of populations and provides predictive data on response to potential future impacts, including those caused by anthropogenic climate change. Here we show how a highly mobile marine species responded to the gain and loss of new breeding habitat. Southern elephant seal, Mirounga leonina, remains were found along the Victoria Land Coast (VLC) in the Ross Sea, Antarctica, 2,500 km from the nearest extant breeding site on Macquarie Island (MQ). This habitat was released after retreat of the grounded ice sheet in the Ross Sea Embayment 7,500–8,000 cal YBP, and is within the range of modern foraging excursions from the MQ colony. Using ancient mtDNA and coalescent models, we tracked the population dynamics of the now extinct VLC colony and the connectivity between this and extant breeding sites. We found a clear expansion signal in the VLC population ∼8,000 YBP, followed by directional migration away from VLC and the loss of diversity at ∼1,000 YBP, when sea ice is thought to have expanded. Our data suggest that VLC seals came initially from MQ and that some returned there once the VLC habitat was lost, ∼7,000 years later. We track the founder-extinction dynamics of a population from inception to extinction in the context of Holocene climate change and present evidence that an unexpectedly diverse, differentiated breeding population was founded from a distant source population soon after habitat became available.

In order to understand how biodiversity is generated and maintained over time, we need to understand the process by which populations form and diverge. Natural variation within species is typically partitioned among populations, which sometimes forms the basis for speciation events. One mechanism for the establishment of novel variation at the population level is through a response to emerging habitat. Here we use data from ancient DNA to show how elephant seal populations responded when new breeding habitat was gained and then lost over the course of approximately 7,000 years. We show that the seals quickly took advantage of newly available breeding habitat, far from the nearest extant breeding site, and that a highly diverse and genetically differentiated population was established over a matter of generations. The key factors were likely the abundant local food resource and extensive physical habitat that allowed rapid expansion after the initial founder event and a tendency for females to return to annual breeding sites in this species. Tracking the founder-extinction dynamics of historical populations provides insight into the likely implications of future environmental change. This is an important tool in our efforts to mitigate the impact of human-induced change.

On the origins and admixture of Malagasy: new evidence from high-resolution analyses of paternal and maternal lineages

The Malagasy have been shown to be a genetically admixed population combining parental lineages with African and South East Asian ancestry. In the present paper, we fit the Malagasy admixture history in a highly resolved phylogeographic framework by typing a large set of mitochondrial DNA and Y DNA markers in unrelated individuals from inland (Merina) and coastal (Antandroy, Antanosy, and Antaisaka) ethnic groups. This allowed performance of a multilevel analysis in which the diversity among main ethnic divisions, lineage ancestries, and modes of inheritance could be concurrently evaluated. Admixture was confirmed to result from the encounter of African and Southeast Asian people with minor recent male contributions from Europe. However, new scenarios are depicted about Malagasy admixture history. The distribution of ancestral components was ethnic and sex biased, with the Asian ancestry appearing more conserved in the female than in the male gene pool and in inland than in coastal groups. A statistic based on haplotype sharing (D(HS)), showing low sampling error and time linearity over the last 200 generations, was introduced here for the first time and helped to integrate our results with linguistic and archeological data. The focus about the origin of Malagasy lineages was enlarged in space and pushed back in time. Homelands could not be pinpointed but appeared to comprise two vast areas containing different populations from sub-Saharan Africa and South East Asia. The pattern of diffusion of uniparental lineages was compatible with at least two events: a primary admixture of proto-Malay people with Bantu speakers bearing a western-like pool of haplotypes, followed by a secondary flow of Southeastern Bantu speakers unpaired for gender (mainly male driven) and geography (mainly coastal).

Correcting for purifying selection: an improved human mitochondrial molecular clock

There is currently no calibration available for the whole human mtDNA genome, incorporating both coding and control regions. Furthermore, as several authors have pointed out recently, linear molecular clocks that incorporate selectable characters are in any case problematic. We here confirm a modest effect of purifying selection on the mtDNA coding region and propose an improved molecular clock for dating human mtDNA, based on a worldwide phylogeny of > 2000 complete mtDNA genomes and calibrating against recent evidence for the divergence time of humans and chimpanzees. We focus on a time-dependent mutation rate based on the entire mtDNA genome and supported by a neutral clock based on synonymous mutations alone. We show that the corrected rate is further corroborated by archaeological dating for the settlement of the Canary Islands and Remote Oceania and also, given certain phylogeographic assumptions, by the timing of the first modern human settlement of Europe and resettlement after the Last Glacial Maximum. The corrected rate yields an age of modern human expansion in the Americas at approximately 15 kya that-unlike the uncorrected clock-matches the archaeological evidence, but continues to indicate an out-of-Africa dispersal at around 55-70 kya, 5-20 ky before any clear archaeological record, suggesting the need for archaeological research efforts focusing on this time window. We also present improved rates for the mtDNA control region, and the first comprehensive estimates of positional mutation rates for human mtDNA, which are essential for defining mutation models in phylogenetic analyses.