The role of selection in the evolution of human mitochondrial genomes

Yeast models of mutations in the mitochondrial ATP6 gene found in human cancer cells

Since the discovery of somatic mtDNA mutations in tumor cells, multiple studies have focused on establishing a causal relationship between those changes and alterations in energy metabolism, a hallmark of cancer cells. Yet the consequences of these mutations on mitochondrial function remain largely unknown. In this study, Saccharomyces cerevisiae has been used as a model to investigate the functional consequences of four cancer-associated missense mutations (8914C>A, 8932C>T, 8953A>G, 9131T>C) found in the mitochondrial MT-ATP6 gene. This gene encodes the a-subunit of F1FO-ATP synthase, which catalyzes the last steps of ATP production in mitochondria. Although the four studied mutations affected well-conserved residues of the a-subunit, only one of them (8932C>T) had a significant impact on mitochondrial function, due to a less efficient incorporation of the a-subunit into ATP synthase. Our findings indicate that these ATP6 genetic variants found in human tumors are neutral mitochondrial genome substitutions with a limited, if any, impact on the energetic function of mitochondria.

Ulcerative Colitis Is Under Dual (Mitochondrial and Nuclear) Genetic Control

BACKGROUND

Cellular oxidative stress and genetic susceptibility have been implicated in the multifactorial etiology of ulcerative colitis (UC). The nuclear genome association with UC has been intensely investigated, but the role of the mitochondrial DNA (mtDNA) has received far less attention and may account for part of the missing heritability. This study is a comprehensive analysis of the mtDNA contribution to UC susceptibility.

METHODS

The association of mitochondrial single-nucleotide polymorphisms (mtSNPs) and haplogroups with UC was tested in 488 cases and 833 controls of European ancestry from the NIDDK IBD Genetics Consortium Ulcerative Colitis Genome-Wide Association Study available through dbGaP and from the Illumina Genotype Control Database (studies 64 and 65).

RESULTS

No evidence of population stratification could be detected using 218 ancestry informative markers for European Americans. Seven of the 58 tested mtSNPs were nominally associated with UC, and A10550G in MT-ND4L withstands the Bonferroni correction (P = 1.29E-06, odds ratio [ORG] [95% confidence interval (CI)] = 4.80 [2.54-9.05], 10550G allele: 8.1% of patients and 1.9% of controls). A10550G remains equally associated after conditional analyses on the 11 UC genome-wide association studies (GWAS) top SNPs (6.35E-07 < Pcond < 4.58E-06), which suggests that it constitutes an independent risk factor from nuclear-encoded susceptibility loci. We detected additive (but not multiplicative) epistatic interactions between A10550G and all 11 top GWAS hits. Subhaplogroup K1 (P = 0.021, OR [95% CI] = 1.71 [1.08-2.69]) increased the risk for UC, whereas the U5b lineage conferred protection (P = 0.016, OR [95% CI] = 0.34 [0.14-0.82]).

CONCLUSIONS

These results suggest that UC has a dual mitochondrial and nuclear genetic control that warrants further replication in independent data sets and reinforces its etiopathogenic complexity.

Massively parallel sequencing of the entire control region and targeted coding region SNPs of degraded mtDNA using a simplified library preparation method

The application of next-generation sequencing (NGS) to forensic genetics is being explored by an increasing number of laboratories because of the potential of high-throughput sequencing for recovering genetic information from multiple markers and multiple individuals in a single run. A cumbersome and technically challenging library construction process is required for NGS. In this study, we propose a simplified library preparation method for mitochondrial DNA (mtDNA) analysis that involves two rounds of PCR amplification. In the first-round of multiplex PCR, six fragments covering the entire mtDNA control region and 22 fragments covering interspersed single nucleotide polymorphisms (SNPs) in the coding region that can be used to determine global haplogroups and East Asian haplogroups were amplified using template-specific primers with read sequences. In the following step, indices and platform-specific sequences for the MiSeq(®) system (Illumina) were added by PCR. The barcoded library produced using this simplified workflow was successfully sequenced on the MiSeq system using the MiSeq Reagent Nano Kit v2. A total of 0.4 GB of sequences, 80.6% with base quality of >Q30, were obtained from 12 degraded DNA samples and mapped to the revised Cambridge Reference Sequence (rCRS). A relatively even read count was obtained for all amplicons, with an average coverage of 5200 × and a less than three-fold read count difference between amplicons per sample. Control region sequences were successfully determined, and all samples were assigned to the relevant haplogroups. In addition, enhanced discrimination was observed by adding coding region SNPs to the control region in in silico analysis. Because the developed multiplex PCR system amplifies small-sized amplicons (<250 bp), NGS analysis using the library preparation method described here allows mtDNA analysis using highly degraded DNA samples.

Endothelial Cell Bioenergetics and Mitochondrial DNA Damage Differ in Humans Having African or West Eurasian Maternal Ancestry

BACKGROUND

We hypothesized that endothelial cells having distinct mitochondrial genetic backgrounds would show variation in mitochondrial function and oxidative stress markers concordant with known differential cardiovascular disease susceptibilities. To test this hypothesis, mitochondrial bioenergetics were determined in endothelial cells from healthy individuals with African versus European maternal ancestries.

METHODS AND RESULTS

Bioenergetics and mitochondrial DNA (mtDNA) damage were assessed in single-donor human umbilical vein endothelial cells belonging to mtDNA haplogroups H and L, representing West Eurasian and African maternal ancestries, respectively. Human umbilical vein endothelial cells from haplogroup L used less oxygen for ATP production and had increased levels of mtDNA damage compared with those in haplogroup H. Differences in bioenergetic capacity were also observed in that human umbilical vein endothelial cells belonging to haplogroup L had decreased maximal bioenergetic capacities compared with haplogroup H. Analysis of peripheral blood mononuclear cells from age-matched healthy controls with West Eurasian or African maternal ancestries showed that haplogroups sharing an A to G mtDNA mutation at nucleotide pair 10398 had increased mtDNA damage compared with those lacking this mutation. Further study of angiographically proven patients with coronary artery disease and age-matched healthy controls revealed that mtDNA damage was associated with vascular function and remodeling and that age of disease onset was later in individuals from haplogroups lacking the A to G mutation at nucleotide pair 10398.

CONCLUSIONS

Differences in mitochondrial bioenergetics and mtDNA damage associated with maternal ancestry may contribute to endothelial dysfunction and vascular disease.

Whole mitochondrial genomes unveil the impact of domestication on goat matrilineal variability

Background

The current extensive use of the domestic goat (Capra hircus) is the result of its medium size and high adaptability as multiple breeds. The extent to which its genetic variability was influenced by early domestication practices is largely unknown. A common standard by which to analyze maternally-inherited variability of livestock species is through complete sequencing of the entire mitogenome (mitochondrial DNA, mtDNA).

Results

We present the first extensive survey of goat mitogenomic variability based on 84 complete sequences selected from an initial collection of 758 samples that represent 60 different breeds of C. hircus, as well as its wild sister species, bezoar (Capra aegagrus) from Iran. Our phylogenetic analyses dated the most recent common ancestor of C. hircus to ~460,000 years (ka) ago and identified five distinctive domestic haplogroups (A, B1, C1a, D1 and G). More than 90 % of goats examined were in haplogroup A. These domestic lineages are predominantly nested within C. aegagrus branches, diverged concomitantly at the interface between the Epipaleolithic and early Neolithic periods, and underwent a dramatic expansion starting from ~12–10 ka ago.

Conclusions

Domestic goat mitogenomes descended from a small number of founding haplotypes that underwent domestication after surviving the last glacial maximum in the Near Eastern refuges. All modern haplotypes A probably descended from a single (or at most a few closely related) female C. aegagrus. Zooarchaelogical data indicate that domestication first occurred in Southeastern Anatolia. Goats accompanying the first Neolithic migration waves into the Mediterranean were already characterized by two ancestral A and C variants. The ancient separation of the C branch (~130 ka ago) suggests a genetically distinct population that could have been involved in a second event of domestication. The novel diagnostic mutational motifs defined here, which distinguish wild and domestic haplogroups, could be used to understand phylogenetic relationships among modern breeds and ancient remains and to evaluate whether selection differentially affected mitochondrial genome variants during the development of economically important breeds.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2342-2) contains supplementary material, which is available to authorized users.

Mitochondrial DNA variation in human radiation and disease

Environmental adaptation, predisposition to common diseases, and, potentially, speciation may all be linked through the adaptive potential of mitochondrial DNA (mtDNA) alterations of bioenergetics. This Perspective synthesizes evidence that human mtDNA variants may be adaptive or deleterious depending on environmental context and proposes that the accrual of mtDNA variation could contribute to animal speciation via adaptation to marginal environments.

Fine Dissection of Human Mitochondrial DNA Haplogroup HV Lineages Reveals Paleolithic Signatures from European Glacial Refugia

Genetic signatures from the Paleolithic inhabitants of Eurasia can be traced from the early divergent mitochondrial DNA lineages still present in contemporary human populations. Previous studies already suggested a pre-Neolithic diffusion of mitochondrial haplogroup HV*(xH,V) lineages, a relatively rare class of mtDNA types that includes parallel branches mainly distributed across Europe and West Asia with a certain degree of structure. Up till now, variation within haplogroup HV was addressed mainly by analyzing sequence data from the mtDNA control region, except for specific sub-branches, such as HV4 or the widely distributed haplogroups H and V. In this study, we present a revised HV topology based on full mtDNA genome data, and we include a comprehensive dataset consisting of 316 complete mtDNA sequences including 60 new samples from the Italian peninsula, a previously underrepresented geographic area. We highlight points of instability in the particular topology of this haplogroup, reconstructed with BEAST-generated trees and networks. We also confirm a major lineage expansion that probably followed the Late Glacial Maximum and preceded Neolithic population movements. We finally observe that Italy harbors a reservoir of mtDNA diversity, with deep-rooting HV lineages often related to sequences present in the Caucasus and the Middle East. The resulting hypothesis of a glacial refugium in Southern Italy has implications for the understanding of late Paleolithic population movements and is discussed within the archaeological cultural shifts occurred over the entire continent.

Estimate of within population incremental selection through branch imbalance in lineage trees

Incremental selection within a population, defined as limited fitness changes following mutation, is an important aspect of many evolutionary processes. Strongly advantageous or deleterious mutations are detected using the synonymous to non-synonymous mutations ratio. However, there are currently no precise methods to estimate incremental selection. We here provide for the first time such a detailed method and show its precision in multiple cases of micro-evolution. The proposed method is a novel mixed lineage tree/sequence based method to detect within population selection as defined by the effect of mutations on the average number of offspring. Specifically, we propose to measure the log of the ratio between the number of leaves in lineage trees branches following synonymous and non-synonymous mutations. The method requires a high enough number of sequences, and a large enough number of independent mutations. It assumes that all mutations are independent events. It does not require of a baseline model and is practically not affected by sampling biases. We show the method's wide applicability by testing it on multiple cases of micro-evolution. We show that it can detect genes and inter-genic regions using the selection rate and detect selection pressures in viral proteins and in the immune response to pathogens.

The complete mitogenome of a 500-year-old Inca child mummy

In 1985, a frozen mummy was found in Cerro Aconcagua (Argentina). Archaeological studies identified the mummy as a seven-year-old Inca sacrifice victim who lived >500 years ago, at the time of the expansion of the Inca Empire towards the southern cone. The sequence of its entire mitogenome was obtained. After querying a large worldwide database of mitogenomes (>28,000) we found that the Inca haplotype belonged to a branch of haplogroup C1b (C1bi) that has not yet been identified in modern Native Americans. The expansion of C1b into the Americas, as estimated using 203 C1b mitogenomes, dates to the initial Paleoindian settlements (~18.3 thousand years ago [kya]); however, its internal variation differs between Mesoamerica and South America. By querying large databases of control region haplotypes (>150,000), we found only a few C1bi members in Peru and Bolivia (e.g. Aymaras), including one haplotype retrieved from ancient DNA of an individual belonging to the Wari Empire (Peruvian Andes). Overall, the results suggest that the profile of the mummy represents a very rare sub-clade that arose 14.3 (5–23.6) kya and could have been more frequent in the past. A Peruvian Inca origin for present-day C1bi haplotypes would satisfy both the genetic and paleo-anthropological findings.

Early Holocenic and Historic mtDNA African Signatures in the Iberian Peninsula: The Andalusian Region as a Paradigm

Determining the timing, identity and direction of migrations in the Mediterranean Basin, the role of "migratory routes" in and among regions of Africa, Europe and Asia, and the effects of sex-specific behaviors of population movements have important implications for our understanding of the present human genetic diversity. A crucial component of the Mediterranean world is its westernmost region. Clear features of transcontinental ancient contacts between North African and Iberian populations surrounding the maritime region of Gibraltar Strait have been identified from archeological data. The attempt to discern origin and dates of migration between close geographically related regions has been a challenge in the field of uniparental-based population genetics. Mitochondrial DNA (mtDNA) studies have been focused on surveying the H1, H3 and V lineages when trying to ascertain north-south migrations, and U6 and L in the opposite direction, assuming that those lineages are good proxies for the ancestry of each side of the Mediterranean. To this end, in the present work we have screened entire mtDNA sequences belonging to U6, M1 and L haplogroups in Andalusians--from Huelva and Granada provinces--and Moroccan Berbers. We present here pioneer data and interpretations on the role of NW Africa and the Iberian Peninsula regarding the time of origin, number of founders and expansion directions of these specific markers. The estimated entrance of the North African U6 lineages into Iberia at 10 ky correlates well with other L African clades, indicating that U6 and some L lineages moved together from Africa to Iberia in the Early Holocene. Still, founder analysis highlights that the high sharing of lineages between North Africa and Iberia results from a complex process continued through time, impairing simplistic interpretations. In particular, our work supports the existence of an ancient, frequently denied, bridge connecting the Maghreb and Andalusia.

The Genomic Legacy of the Transatlantic Slave Trade in the Yungas Valley of Bolivia

During the period of the Transatlantic Slave Trade (TAST) some enslaved Africans were forced to move to Upper Peru (nowadays Bolivia). At first they were sent to Potosí, but later to the tropical Yungas valley where the Spanish colonizers established a so-called “hacienda system” that was based on slave labor, including African-descendants. Due to their isolation, very little attention has been paid so far to ‘Afro-Bolivian’ communities either within the research field of TAST or in genetic population studies. In this study, a total of 105 individuals from the Yungas were sequenced for their mitochondrial DNA (mtDNA) control region, and mitogenomes were obtained for a selected subset of these samples. We also genotyped 46 Ancestry Informative Markers (AIM) in order to investigate continental ancestry at the autosomal level. In addition, Y-chromosome STR and SNP data for a subset of the same individuals was also available from the literature. The data indicate that the partitioning of mtDNA ancestry in the Yungas differs significantly from that in the rest of the country: 81% Native American, 18% African, and 1% European. Interestingly, the great majority of ‘Afro-descendant’ mtDNA haplotypes in the Yungas (84%) concentrates in the locality of Tocaña. This high proportion of African ancestry in the Tocaña is also manifested in the Y-chromosome (44%) and in the autosomes (56%). In sharp contrast with previous studies on the TAST, the ancestry of about 1/3 of the ‘Afro-Bolivian’ mtDNA haplotypes can be traced back to East and South East Africa, which may be at least partially explained by the Arab slave trade connected to the TAST.

A study of genetic polymorphisms in mitochondrial DNA hypervariable regions I and II of the five major ethnic groups and Vedda population in Sri Lanka

Diversity of the hypervariable regions (HV) I and II of the mitochondrial genome was studied in maternally unrelated Sri Lankans (N=202) from six ethnic groups (i.e.: Sinhalese, Sri Lankan Tamil, Muslim, Malay, Indian Tamil and Vedda). DNA was extracted from blood and buccal swabs and HVI and HVII regions were PCR amplified and sequenced. Resulting sequences were aligned and edited between 16024-16365 and 73-340 regions and compared with revised Cambridge reference sequences (rCRS). One hundred and thirty-five unique haplotypes and 22 shared haplotypes were observed. A total of 145 polymorphic sites and 158 polymorphisms were observed. Hypervariable region I showed a higher polymorphic variation than hypervariable region II. Nucleotide diversities were quite low and similar for all ethnicities apart from a slightly higher value for Indian Tamils and a much lower value for the Vedda population compared to the other groups. When the total population was considered South Asian (Indian) haplogroups were predominant, but there were differences in the distribution of phylo-geographical haplogroups between ethnic groups. Sinhalese, Sri Lankan Tamil and Vedda populations had a considerable presence of West Eurasian haplogroups. About 2/3rd of the Vedda population comprised of macro-haplogroup N or its subclades R and U, whereas macro-haplogroup M was predominant in all other populations. The Vedda population clustered separately from other groups and Sri Lankan Tamils showed a closer genetic affiliation to Sinhalese than to Indian Tamils. Thus this study provides useful information for forensic analysis and anthropological studies of Sri Lankans.

Carriers of Mitochondrial DNA Macrohaplogroup N Lineages Reached Australia around 50,000 Years Ago following a Northern Asian Route

Background

The modern human colonization of Eurasia and Australia is mostly explained by a single-out-of-Africa exit following a southern coastal route throughout Arabia and India. However, dispersal across the Levant would better explain the introgression with Neanderthals, and more than one exit would fit better with the different ancient genomic components discovered in indigenous Australians and in ancient Europeans. The existence of an additional Northern route used by modern humans to reach Australia was previously deduced from the phylogeography of mtDNA macrohaplogroup N. Here, we present new mtDNA data and new multidisciplinary information that add more support to this northern route.

Methods

MtDNA hypervariable segments and haplogroup diagnostic coding positions were analyzed in 2,278 Saudi Arabs, from which 1,725 are new samples. Besides, we used 623 published mtDNA genomes belonging to macrohaplogroup N, but not R, to build updated phylogenetic trees to calculate their coalescence ages, and more than 70,000 partial mtDNA sequences were screened to establish their respective geographic ranges.

Results

The Saudi mtDNA profile confirms the absence of autochthonous mtDNA lineages in Arabia with coalescence ages deep enough to support population continuity in the region since the out-of-Africa episode. In contrast to Australia, where N(xR) haplogroups are found in high frequency and with deep coalescence ages, there are not autochthonous N(xR) lineages in India nor N(xR) branches with coalescence ages as deep as those found in Australia. These patterns are at odds with the supposition that Australian colonizers harboring N(xR) lineages used a route involving India as a stage. The most ancient N(xR) lineages in Eurasia are found in China, and inconsistently with the coastal route, N(xR) haplogroups with the southernmost geographical range have all more recent radiations than the Australians.

Conclusions

Apart from a single migration event via a southern route, phylogeny and phylogeography of N(xR) lineages support that people carrying mtDNA N lineages could have reach Australia following a northern route through Asia. Data from other disciplines also support this scenario.

Mitochondrial DNA haplogroup phylogeny of the dog: Proposal for a cladistic nomenclature

Canis lupus familiaris mitochondrial DNA analysis has increased in recent years, not only for the purpose of deciphering dog domestication but also for forensic genetic studies or breed characterization. The resultant accumulation of data has increased the need for a normalized and phylogenetic-based nomenclature like those provided for human maternal lineages. Although a standardized classification has been proposed, haplotype names within clades have been assigned gradually without considering the evolutionary history of dog mtDNA. Moreover, this classification is based only on the D-loop region, proven to be insufficient for phylogenetic purposes due to its high number of recurrent mutations and the lack of relevant information present in the coding region. In this study, we design 1) a refined mtDNA cladistic nomenclature from a phylogenetic tree based on complete sequences, classifying dog maternal lineages into haplogroups defined by specific diagnostic mutations, and 2) a coding region SNP analysis that allows a more accurate classification into haplogroups when combined with D-loop sequencing, thus improving the phylogenetic information obtained in dog mitochondrial DNA studies.

West Eurasian mtDNA lineages in India: an insight into the spread of the Dravidian language and the origins of the caste system

There is no indication from the previous mtDNA studies that west Eurasian-specific subclades have evolved within India and played a role in the spread of languages and the origins of the caste system. To address these issues, we have screened 14,198 individuals (4208 from this study) and analyzed 112 mitogenomes (41 new sequences) to trace west Eurasian maternal ancestry. This has led to the identification of two autochthonous subhaplogroups--HV14a1 and U1a1a4, which are likely to have originated in the Dravidian-speaking populations approximately 10.5-17.9 thousand years ago (kya). The carriers of these maternal lineages might have settled in South India during the time of the spread of the Dravidian language. In addition to this, we have identified several subsets of autochthonous U7 lineages, including U7a1, U7a2b, U7a3, U7a6, U7a7, and U7c, which seem to have originated particularly in the higher-ranked caste populations in relatively recent times (2.6-8.0 kya with an average of 5.7 kya). These lineages have provided crucial clues to the differentiation of the caste system that has occurred during the recent past and possibly, this might have been influenced by the Indo-Aryan migration. The remaining west Eurasian lineages observed in the higher-ranked caste groups, like the Brahmins, were found to cluster with populations who possibly arrived from west Asia during more recent times.

Maternal ancestry and population history from whole mitochondrial genomes

MtDNA has been a widely used tool in human evolutionary and population genetic studies over the past three decades. Its maternal inheritance and lack of recombination have offered the opportunity to explore genealogical relationships among individuals and to study the frequency differences of matrilineal clades among human populations at continental and regional scales. The whole mtDNA genome sequencing delivers molecular resolution that is sufficient to distinguish patterns that have arisen over thousands of years. However, mutation rate is highly variable among the functional and non-coding domains of mtDNA which makes it challenging to obtain accurate split dates of the mitochondrial clades. Due to the shallow coalescent time of mitochondrial TMRCA at approximately 100 to 200 thousand years (ky), mtDNA data have only limited power to inform us about the more distant past and the early stages of human evolutionary history. The variation shared by mitochondrial genomes of individuals drawn from different continents outside Africa has been used to illuminate the details of the colonization process of the Old World, whereas regional patterns of variation have been at the focus of studies addressing questions of a more recent time scale. In the era of whole nuclear genome sequencing, mitochondrial genomes are continuing to be informative as a unique tool for the assessment of female-specific aspects of the demographic history of human populations.

A mitogenomic phylogeny and genetic history of sable (Martes zibellina)

We assessed phylogeny of sable (Martes zibellina, Linnaeus, 1758) by sequence analysis of nearly complete, new mitochondrial genomes in 36 specimens from different localities in northern Eurasia (Primorye, Khabarovsk and Krasnoyarsk regions, the Kamchatka Peninsula, the Kuril Islands and the Urals). Phylogenetic analysis of mtDNA sequences demonstrates that two clades, A and BC, radiated about 200-300 thousandyears ago (kya) according to results of Bayesian molecular clock and RelTime analyses of different mitogenome alignments (nearly complete mtDNA sequences, protein-coding region, and synonymous sites), while the age estimates of clades A, B and C fall within the Late Pleistocene (~50-140 kya). Bayesian skyline plots (BSPs) of sable population size change based on analysis of nearly complete mtDNAs show an expansion around 40 kya in the warm Karganian time, without a decline of population size around the Last Glacial Maximum (21 kya). The BSPs based on synonymous clock rate indicate that M. zibellina experienced demographic expansions later, approximately 22 kya. The A2a clade that colonized Kamchatka ~23-50 kya (depending on the mutation rate used) survived the last glaciation there as demonstrated by the BSP analysis. In addition, we have found evidence of positive selection acting at ND4 and cytochrome b genes, thereby suggesting adaptive evolution of the A2a clade in Kamchatka.

Improved calibration of the human mitochondrial clock using ancient genomes

Reliable estimates of the rate at which DNA accumulates mutations (the substitution rate) are crucial for our understanding of the evolution and past demography of virtually any species. In humans, there are considerable uncertainties around these rates, with substantial variation among recent published estimates. Substitution rates have traditionally been estimated by associating dated events to the root (e.g., the divergence between humans and chimpanzees) or to internal nodes in a phylogenetic tree (e.g., first entry into the Americas). The recent availability of ancient mitochondrial DNA sequences allows for a more direct calibration by assigning the age of the sequenced samples to the tips within the human phylogenetic tree. But studies also vary greatly in the methodology employed and in the sequence panels analyzed, making it difficult to tease apart the causes for the differences between previous estimates. To clarify this issue, we compiled a comprehensive data set of 350 ancient and modern human complete mitochondrial DNA genomes, among which 146 were generated for the purpose of this study and estimated substitution rates using calibrations based both on dated nodes and tips. Our results demonstrate that, for the same data set, estimates based on individual dated tips are far more consistent with each other than those based on nodes and should thus be considered as more reliable.

The episode of genetic drift defining the migration of humans out of Africa is derived from a large east African population size

Human genetic variation particularly in Africa is still poorly understood. This is despite a consensus on the large African effective population size compared to populations from other continents. Based on sequencing of the mitochondrial Cytochrome C Oxidase subunit II (MT-CO2), and genome wide microsatellite data we observe evidence suggesting the effective size (Ne) of humans to be larger than the current estimates, with a foci of increased genetic diversity in east Africa, and a population size of east Africans being at least 2-6 fold larger than other populations. Both phylogenetic and network analysis indicate that east Africans possess more ancestral lineages in comparison to various continental populations placing them at the root of the human evolutionary tree. Our results also affirm east Africa as the likely spot from which migration towards Asia has taken place. The study reflects the spectacular level of sequence variation within east Africans in comparison to the global sample, and appeals for further studies that may contribute towards filling the existing gaps in the database. The implication of these data to current genomic research, as well as the need to carry out defined studies of human genetic variation that includes more African populations; particularly east Africans is paramount.

Dog mitochondrial genome sequencing to enhance dog mtDNA discrimination power in forensic casework

A Belgian dog population sample and several population studies worldwide have confirmed that only a limited number of mtDNA control region haplotypes is observed in the majority of dogs. The high population frequency of these haplotypes negatively impacts both the exclusion probability of dog mtDNA analysis and the evidential value of a match with one of these haplotypes in casework. Variation within the mtDNA coding region was explored to improve the discrimination power of dog mtDNA analysis. In the current study, the entire mitochondrial genome of 161 dogs was sequenced applying a quality assured strategy and resulted in a total of 119 different mitochondrial genome sequences. Our research was focused on those dogs with the six most common control region haplotypes from a previous Belgian population study. We identified 33 informative SNPs that successfully divide the six most common control region haplotypes into 32 clusters of mitochondrial genome sequences. Determining the identity of these 33 polymorphic sites in addition to control region sequencing in case of a match with one of these 6 control region haplotypes could augment the exclusion probability of forensic dog mtDNA analysis from 92.5% to 97.5%.

Complete mitogenome analysis of indigenous populations in Mexico: its relevance for the origin of Mesoamericans

Mesoamerica has an important role in the expansion of Paleoamericans as the route to South America. In this study, we determined complete mitogenome sequences of 113 unrelated individuals from two indigenous populations of Mesoamerica, Mazahua and Zapotec. All newly sequenced mitogenomes could be classified into haplogroups A2, B2, C1 and D1, but one sequence in Mazahua was D4h3a, a subclade of haplogroup D4. This haplogroup has been mostly found in South America along the Pacific coast. Haplogroup X2a was not found in either population. Genetic similarity obtained using phylogenetic tree construction and principal component analysis showed that these two populations are distantly related to each other. Actually, the Mazahua and the Zapotec shared no sequences (haplotypes) in common, while each also showed a number of unique subclades. Surprisingly, Zapotec formed a cluster with indigenous populations living in an area from central Mesoamerica to Central America. By contrast, the Mazahua formed a group with indigenous populations living in external areas, including southwestern North America and South America. This intriguing genetic relationship suggests the presence of two paleo-Mesoamerican groups, invoking a scenario in which one group had expanded into South America and the other resided in Mesoamerica.

A phylogeny for the pomatiopsidae (Gastropoda: Rissooidea): a resource for taxonomic, parasitological and biodiversity studies

BACKGROUND

The Pomatiopsidae are reported from northern India into southern China and Southeast Asia, with two sub-families, the Pomatiopsinae (which include freshwater, amphibious, terrestrial and marine species) and the freshwater Triculinae. Both include species acting as intermediate host for species of the blood-fluke Schistosoma which cause a public health problem in East Asia. Also, with around 120 species, triculine biodiversity exceeds that of any other endemic freshwater molluscan fauna. Nevertheless, the origins of the Pomatiopsidae, the factors driving such a diverse radiation and aspects of their co-evolution with Schistosoma are not fully understood. Many taxonomic questions remain; there are problems identifying medically relevant species. The predicted range is mostly unsurveyed and the true biodiversity of the family is underestimated. Consequently, the aim of the study was to collect DNA-sequence data for as many pomatiopsid taxa as possible, as a first step in providing a resource for identification of epidemiologically significant species (by non-malacologists), for use in resolving taxonomic confusion and for testing phylogeographical hypotheses.

RESULTS

The evolutionary radiation of the Triculinae was shown to have been rapid and mostly post late Miocene. Molecular dating indicated that the radiation of these snails was driven first by the uplift of the Himalaya and onset of a monsoon system, and then by late-Pliocene global warming. The status of Erhaia as Anmicolidae is supported. The genera Tricula and Neotricula are shown to be non-monophyletic and the tribe Jullieniini may be polyphyletic (based on convergent characters). Triculinae from northern Vietnam could be derived from Gammatricula of Fujian/Yunnan, China.

CONCLUSIONS

The molecular dates and phylogenetic estimates in this study are consistent with an Australasian origin for the Pomatiopsidae and an East to West radiation via Oligocene Borneo-Philippines island hopping to Japan and then China (Triculinae arising mid-Miocene in Southeast China), and less so with a triculine origin in Tibet. The lack of monophyly in the medically important genera and indications of taxonomic inaccuracies, call for further work to identify epidemiologically significant taxa (e.g., Halewisia may be potential hosts for Schistosoma mekongi) and highlight the need for surveys to determine the true biodiversity of the Triculinae.

The impact of mitochondrial DNA and nuclear genes related to mitochondrial functioning on the risk of Parkinson's disease

Mitochondrial dysfunction and oxidative stress are the major factors implicated in Parkinson’s disease (PD) pathogenesis. The maintenance of healthy mitochondria is a very complex process coordinated bi-genomically. Here, we review association studies on mitochondrial haplogroups and subhaplogroups, discussing the underlying molecular mechanisms. We also focus on variation in the nuclear genes (NDUFV2, PGC-1alpha, HSPA9, LRPPRC, MTIF3, POLG1, and TFAM encoding NADH dehydrogenase (ubiquinone) flavoprotein 2, peroxisome proliferator-activated receptor gamma coactivator 1-alpha, mortalin, leucine-rich pentatricopeptide repeat containing protein, translation initiation factor 3, mitochondrial DNA polymerase gamma, and mitochondrial transcription factor A, respectively) primarily linked to regulation of mitochondrial functioning that recently have been associated with PD risk. Possible interactions between mitochondrial and nuclear genetic variants and related proteins are discussed.

Large-scale mitochondrial DNA analysis in Southeast Asia reveals evolutionary effects of cultural isolation in the multi-ethnic population of Myanmar

Background

Myanmar is the largest country in mainland Southeast Asia with a population of 55 million people subdivided into more than 100 ethnic groups. Ruled by changing kingdoms and dynasties and lying on the trade route between India and China, Myanmar was influenced by numerous cultures. Since its independence from British occupation, tensions between the ruling Bamar and ethnic minorities increased.

Results

Our aim was to search for genetic footprints of Myanmar’s geographic, historic and sociocultural characteristics and to contribute to the picture of human colonization by describing and dating of new mitochondrial DNA (mtDNA) haplogroups. Therefore, we sequenced the mtDNA control region of 327 unrelated donors and the complete mitochondrial genome of 44 selected individuals according to highest quality standards.

Conclusion

Phylogenetic analyses of the entire mtDNA genomes uncovered eight new haplogroups and three unclassified basal M-lineages. The multi-ethnic population and the complex history of Myanmar were reflected in its mtDNA heterogeneity. Population genetic analyses of Burmese control region sequences combined with population data from neighboring countries revealed that the Myanmar haplogroup distribution showed a typical Southeast Asian pattern, but also Northeast Asian and Indian influences. The population structure of the extraordinarily diverse Bamar differed from that of the Karen people who displayed signs of genetic isolation. Migration analyses indicated a considerable genetic exchange with an overall positive migration balance from Myanmar to neighboring countries. Age estimates of the newly described haplogroups point to the existence of evolutionary windows where climatic and cultural changes gave rise to mitochondrial haplogroup diversification in Asia.

The case for the continuing use of the revised Cambridge Reference Sequence (rCRS) and the standardization of notation in human mitochondrial DNA studies

Since the determination in 1981 of the sequence of the human mitochondrial DNA (mtDNA) genome, the Cambridge Reference Sequence (CRS), has been used as the reference sequence to annotate mtDNA in molecular anthropology, forensic science and medical genetics. The CRS was eventually upgraded to the revised version (rCRS) in 1999. This reference sequence is a convenient device for recording mtDNA variation, although it has often been misunderstood as a wild-type (WT) or consensus sequence by medical geneticists. Recently, there has been a proposal to replace the rCRS with the so-called Reconstructed Sapiens Reference Sequence (RSRS). Even if it had been estimated accurately, the RSRS would be a cumbersome substitute for the rCRS, as the new proposal fuses--and thus confuses--the two distinct concepts of ancestral lineage and reference point for human mtDNA. Instead, we prefer to maintain the rCRS and to report mtDNA profiles by employing the hitherto predominant circumfix style. Tree diagrams could display mutations by using either the profile notation (in conventional short forms where appropriate) or in a root-upwards way with two suffixes indicating ancestral and derived nucleotides. This would guard against misunderstandings about reporting mtDNA variation. It is therefore neither necessary nor sensible to change the present reference sequence, the rCRS, in any way. The proposed switch to RSRS would inevitably lead to notational chaos, mistakes and misinterpretations.

Evolutionary perspectives on the links between mitochondrial genotype and disease phenotype

BACKGROUND

Disorders of the mitochondrial respiratory chain are heterogeneous in their symptoms and underlying genetics. Simple links between candidate mutations and expression of disease phenotype typically do not exist. It thus remains unclear how the genetic variation in the mitochondrial genome contributes to the phenotypic expression of complex traits and disease phenotypes.

SCOPE OF REVIEW

I summarize the basic genetic processes known to underpin mitochondrial disease. I highlight other plausible processes, drawn from the evolutionary biological literature, whose contribution to mitochondrial disease expression remains largely empirically unexplored. I highlight recent advances to the field, and discuss common-ground and -goals shared by researchers across medical and evolutionary domains.

MAJOR CONCLUSIONS

Mitochondrial genetic variance is linked to phenotypic variance across a variety of traits (e.g. reproductive function, life expectancy) fundamental to the upkeep of good health. Evolutionary theory predicts that mitochondrial genomes are destined to accumulate male-harming (but female-friendly) mutations, and this prediction has received proof-of-principle support. Furthermore, mitochondrial effects on the phenotype are typically manifested via interactions between mitochondrial and nuclear genes. Thus, whether a mitochondrial mutation is pathogenic in effect can depend on the nuclear genotype in which is it expressed.

GENERAL SIGNIFICANCE

Many disease phenotypes associated with OXPHOS malfunction might be determined by the outcomes of mitochondrial-nuclear interactions, and by the evolutionary forces that historically shaped mitochondrial DNA (mtDNA) sequences. Concepts and results drawn from the evolutionary sciences can have broad, but currently under-utilized, applicability to the medical sciences and provide new insights into understanding the complex genetics of mitochondrial disease. This article is part of a Special Issue entitled Frontiers of Mitochondrial Research.

Mitochondrial DNA content contributes to climate adaptation using Chinese populations as a model

Maintaining a balance between ATP synthesis and heat generation is crucial for adapting to changes in climate. Variation in the mitochondrial DNA (mtDNA), which encodes 13 subunits of the respiratory chain complexes, may contribute to climate adaptation by regulating thermogenesis and the use of bioenergy. However, studies looking for a relationship between mtDNA haplogroups and climate have obtained mixed results, leaving unresolved the role of mtDNA in climate adaptation. Since mtDNA content can regulate human bioenergy processes and is known to influence many physiological traits and diseases, it is possible that mtDNA content contributes to climate adaptation in human populations. Here, we analyze the distribution of mtDNA content among 27 Chinese ethnic populations residing across China and find a significant association between mtDNA content and climate, with northern populations having significantly higher mtDNA content than southern populations. Functional studies have shown that high mtDNA content correlates with an increase in the expression of energy metabolism enzymes, which may accelerate thermogenesis. This suggests that the significantly higher mtDNA content observed in northern populations may confer a selective advantage in adapting to colder northern climates

Frequency and pattern of heteroplasmy in the complete human mitochondrial genome

Determining the levels of human mitochondrial heteroplasmy is of utmost importance in several fields. In spite of this, there are currently few published works that have focused on this issue. In order to increase the knowledge of mitochondrial DNA (mtDNA) heteroplasmy, the main goal of this work is to investigate the frequency and the mutational spectrum of heteroplasmy in the human mtDNA genome. To address this, a set of nine primer pairs designed to avoid co-amplification of nuclear DNA (nDNA) sequences of mitochondrial origin (NUMTs) was used to amplify the mitochondrial genome in 101 individuals. The analysed individuals represent a collection with a balanced representation of genders and mtDNA haplogroup distribution, similar to that of a Western European population. The results show that the frequency of heteroplasmic individuals exceeds 61%. The frequency of point heteroplasmy is 28.7%, with a widespread distribution across the entire mtDNA. In addition, an excess of transitions in heteroplasmy were detected, suggesting that genetic drift and/or selection may be acting to reduce its frequency at population level. In fact, heteroplasmy at highly stable positions might have a greater impact on the viability of mitochondria, suggesting that purifying selection must be operating to prevent their fixation within individuals. This study analyses the frequency of heteroplasmy in a healthy population, carrying out an evolutionary analysis of the detected changes and providing a new perspective with important consequences in medical, evolutionary and forensic fields.

The other genome: a systematic review of studies of mitochondrial DNA haplogroups and outcomes of HIV infection and antiretroviral therapy

Mitochondrial toxicity is implicated in some treatment-limiting antiretroviral therapy complications, and reports of mitochondrial dysfunction in untreated HIV infection suggest antiretroviral therapy independent effects of HIV. Several studies have explored associations between mtDNA haplogroups (patterns of mtDNA polymorphisms) and outcomes of HIV infection and/or antiretroviral therapy, but findings have been inconsistent. We systematically reviewed published studies examining mtDNA haplogroups in HIV-infected persons to summarize reported outcome associations, and to highlight potential future research directions. We identified 21 articles published from 2005-2013. Multiple different phenotypes were studied; most were antiretroviral therapy associated metabolic outcomes (e.g. lipodystrophy, insulin resistance, and dyslipidemia). Haplogroup H was associated with the most outcomes, including AIDS progression, CD4 T-cell recovery, cirrhosis (in hepatitis C coinfection), and metabolic outcomes. This review is the first to focus on the emerging area of mtDNA haplogroups in HIV, and summarizes the published literature on associations between mtDNA haplogroups and clinical outcomes in populations of European and African descent. Several reported associations require replication and ideally biological verification before definitive conclusions can be drawn, but research in this area has the potential to explain outcome disparities and impact clinical management of patients.

Origin of HBV and its arrival in the Americas--the importance of natural selection on time estimates

BACKGROUND

The strong geographic structure shown by the global pattern of HBV lineages suggests an ancient origin for this virus; however, estimates based on the molecular clock suggest a very recent origin for the Native American genotypes F and H. In this study, we contribute to this debate by estimating the divergence times of genotypes F and H and by discussing how evolutionary rates estimated from recent samples may underestimate the divergence time of more ancient nodes in HBV phylogenies.

METHODS

A total of 108 complete HBV genotype F and H genomes were compared to 44 reference genomes from other genotypes. Time estimates were based on a Bayesian method with evolutionary rates taken from the literature. To assess the pattern of substitutions in recent versus old branches we mapped the phylogenetic distribution of all mutations occurring in genotypes F and H using a maximum likelihood approach and compared the number of synonymous and non-synonymous mutations in young and old branches of HBV genotype F and H phylogeny using a χ² test.

RESULTS

Estimated divergence times between genotypes F and H depend heavily on the evolutionary rate. While fast rates suggest a recent separation of these genotypes (approximately 800 years ago), slow rates suggest an earlier divergence (up to approximately 13,000 years ago). There is a clear excess of non-synonymous substitutions in the most recent branches of HBV phylogeny (P=4.87×10⁻¹⁵), most likely suggesting the action of purifying selection.

CONCLUSIONS

These results suggest that rates estimated based on recent samples will overestimate the evolutionary rate and underestimate the coalescence times for ancient nodes in HBV phylogeny.

Selective constraints determine the time dependency of molecular rates for human nuclear genomes

In contrast to molecular rates for neutral mitochondrial sequences, rates for constrained sites (including nonsynonymous sites, D-loop, and RNA) in the mitochondrial genome are known to vary with the time frame used for their estimation. Here, we examined this issue for the nuclear genomes using single-nucleotide polymorphisms (SNPs) from six complete human genomes of individuals belonging to different populations. We observed a strong time-dependent distribution of nonsynonymous SNPs (nSNPs) in highly constrained genes. Typically, the proportion of young nSNPs specific to a single population was found to be up to three times higher than that of the ancient nSNPs shared between diverse human populations. In contrast, this trend disappeared, and a uniform distribution of young and old nSNPs was observed in genes under relaxed selective constraints. This suggests that because mutations in constrained genes are highly deleterious, they are removed over time, resulting in a relative overabundance of young nSNPs. In contrast, mutations in genes under relaxed constraints are nearly neutral, which leads to similar proportions of young and old SNPs. These results could be useful to researchers aiming to select appropriate genes or genomic regions for estimating evolutionary rates and species or population divergence times.

Evaluating purifying selection in the mitochondrial DNA of various mammalian species

Mitochondrial DNA (mtDNA), the circular DNA molecule inside the mitochondria of all eukaryotic cells, has been shown to be under the effect of purifying selection in several species. Traditional testing of purifying selection has been based simply on ratios of nonsynonymous to synonymous mutations, without considering the relative age of each mutation, which can be determined by phylogenetic analysis of this non-recombining molecule. The incorporation of a mutation time-ordering from phylogeny and of predicted pathogenicity scores for nonsynonymous mutations allow a quantitative evaluation of the effects of purifying selection in human mtDNA. Here, by using this additional information, we show that purifying selection undoubtedly acts upon the mtDNA of other mammalian species/genera, namely Bos sp., Canis lupus, Mus musculus, Orcinus orca, Pan sp. and Sus scrofa. The effects of purifying selection were comparable in all species, leading to a significant major proportion of nonsynonymous variants with higher pathogenicity scores in the younger branches of the tree. We also derive recalibrated mutation rates for age estimates of ancestors of these various species and proposed a correction curve in order to take into account the effects of selection. Understanding this selection is fundamental to evolutionary studies and to the identification of deleterious mutations.

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).

The Q2 mitochondrial haplogroup in Oceania

Many details surrounding the origins of the peoples of Oceania remain to be resolved, and as a step towards this we report seven new complete mitochondrial genomes from the Q2a haplogroup, from Papua New Guinea, Fiji and Kiribati. This brings the total to eleven Q2 genomes now available. The Q haplogroup (that includes Q2) is an old and diverse lineage in Near Oceania, and is reasonably common; within our sample set of 430, 97 are of the Q haplogroup. However, only 8 are Q2, and we report 7 here. The tree with all complete Q genomes is proven to be minimal. The dating estimate for the origin of Q2 (around 35 Kya) reinforces the understanding that humans have been in Near Oceania for tens of thousands of years; nevertheless the Polynesian maternal haplogroups remain distinctive. A major focus now, with regard to Polynesian ancestry, is to address the differences and timing of the ‘Melanesian’ contribution to the maternal and paternal lineages as people moved further and further into Remote Oceania. Input from other fields such as anthropology, history and linguistics is required for a better understanding and interpretation of the genetic data.

Chicken domestication: an updated perspective based on mitochondrial genomes

Domestic chickens (Gallus gallus domesticus) fulfill various roles ranging from food and entertainment to religion and ornamentation. To survey its genetic diversity and trace the history of domestication, we investigated a total of 4938 mitochondrial DNA (mtDNA) fragments including 2843 previously published and 2095 de novo units from 2044 domestic chickens and 51 red junglefowl (Gallus gallus). To obtain the highest possible level of molecular resolution, 50 representative samples were further selected for total mtDNA genome sequencing. A fine-gained mtDNA phylogeny was investigated by defining haplogroups A-I and W-Z. Common haplogroups A-G were shared by domestic chickens and red junglefowl. Rare haplogroups H-I and W-Z were specific to domestic chickens and red junglefowl, respectively. We re-evaluated the global mtDNA profiles of chickens. The geographic distribution for each of major haplogroups was examined. Our results revealed new complexities of history in chicken domestication because in the phylogeny lineages from the red junglefowl were mingled with those of the domestic chickens. Several local domestication events in South Asia, Southwest China and Southeast Asia were identified. The assessment of chicken mtDNA data also facilitated our understanding about the Austronesian settlement in the Pacific.

Divorcing the Late Upper Palaeolithic demographic histories of mtDNA haplogroups M1 and U6 in Africa

Background

A Southwest Asian origin and dispersal to North Africa in the Early Upper Palaeolithic era has been inferred in previous studies for mtDNA haplogroups M1 and U6. Both haplogroups have been proposed to show similar geographic patterns and shared demographic histories.

Results

We report here 24 M1 and 33 U6 new complete mtDNA sequences that allow us to refine the existing phylogeny of these haplogroups. The resulting phylogenetic information was used to genotype a further 131 M1 and 91 U6 samples to determine the geographic spread of their sub-clades. No southwest Asian specific clades for M1 or U6 were discovered. U6 and M1 frequencies in North Africa, the Middle East and Europe do not follow similar patterns, and their sub-clade divisions do not appear to be compatible with their shared history reaching back to the Early Upper Palaeolithic. The Bayesian Skyline Plots testify to non-overlapping phases of expansion, and the haplogroups’ phylogenies suggest that there are U6 sub-clades that expanded earlier than those in M1. Some M1 and U6 sub-clades could be linked with certain events. For example, U6a1 and M1b, with their coalescent ages of ~20,000–22,000 years ago and earliest inferred expansion in northwest Africa, could coincide with the flourishing of the Iberomaurusian industry, whilst U6b and M1b1 appeared at the time of the Capsian culture.

Conclusions

Our high-resolution phylogenetic dissection of both haplogroups and coalescent time assessments suggest that the extant main branching pattern of both haplogroups arose and diversified in the mid-later Upper Palaeolithic, with some sub-clades concomitantly with the expansion of the Iberomaurusian industry. Carriers of these maternal lineages have been later absorbed into and diversified further during the spread of Afro-Asiatic languages in North and East Africa.