Mitochondrial DNA Population Research Articles

Bayesian Emulation and Calibration of a Stochastic Computer Model of Mitochondrial DNA Deletions in Substantia Nigra Neurons

This article considers the problem of parameter estimation for a stochastic biological model of mitochondrial DNA population dynamics using experimental data on deletion mutation accumulation. The stochastic model is an attempt to describe the hypothesized link between deletion accumulation and neuronal loss in the substantia nigra region of the human brain. Inference for the parameters of the model is complicated by the fact that the model is both analytically intractable and slow to sample from. We show how the stochastic model can be approximated using a simple parametric statistical model with smoothly varying parameters. These parameters are treated as unknown functions and modeled using Gaussian process priors. Several simplifications of our Bayesian model are implemented to ease the computational burden. Throughout the article, we validate our models using predictive simulations. We demonstrate the validity of our fitted model on an independent dataset of substantia nigra neuron survival.

Evidence for Variation in the Effective Population Size of Animal Mitochondrial DNA

BackgroundIt has recently been shown that levels of diversity in mitochondrial DNA are remarkably constant across animals of diverse census population sizes and ecologies, which has led to the suggestion that the effective population of mitochondrial DNA may be relatively constant.ResultsHere we present several lines of evidence that suggest, to the contrary, that the effective population size of mtDNA does vary, and that the variation can be substantial. First, we show that levels of mitochondrial and nuclear diversity are correlated within all groups of animals we surveyed. Second, we show that the effectiveness of selection on non-synonymous mutations, as measured by the ratio of the numbers of non-synonymous and synonymous polymorphisms, is negatively correlated to levels of mitochondrial diversity. Finally, we estimate the effective population size of mitochondrial DNA in selected mammalian groups and show that it varies by at least an order of magnitude.ConclusionsWe conclude that there is variation in the effective population size of mitochondria. Furthermore we suggest that the relative constancy of DNA diversity may be due to a negative correlation between the effective population size and the mutation rate per generation.

Mitochondrial DNA diversity and population structure of a forest-dependent rodent, Praomys taitae (Rodentia: Muridae) Heller 1911, in the fragmented forest patches of Taita Hills, Kenya

The population genetic structure of the forest-dependent rodent, Praomys taitae, sampled from nine indigenous forest fragments distributed over three ranges of the Taita Hills in Kenya, was determined using mitochondrial DNA (mtDNA) control region sequence variation. Thirty-three unique haplotypes were observed in a total sample of 132 individuals, whereas the number of different haplotypes per population ranged from two to 10. An overall low nucleotide diversity of 0.9% was observed in the total sample but varied widely between populations (0.2–1.3%). Significant genetic differentiation was observed in 30 of the 36 possible pair-wise comparisons between populations, while a hierarchical AMOVA revealed significant genetic subdivision between groups of populations on the three hill ranges of Dabida, Mbololo and Kyulu (FCT = 0.404, P less than 0.01), among populations on each of the hill ranges (FSC = 0.112, P less than 0.01) and among populations in the total sample (FST = 0.471, P less than 0.001). Demographic history analyses based on pair-wise nucleotide sequence mismatch distributions revealed that all the populations were in mutation-drift disequilibrium except the populations of the Kyulu and Ronge forest fragments.

Mitochondrial DNA population data of HV1 and HV2 sequences from Japanese individuals

Mitochondrial DNA sequences of the hypervariable regions HV1 and HV2 were determined for 1204 unrelated Japanese individuals. A total of 741 different mtDNA haplotypes were found, 157 of which were seen in multiple individuals. Twenty-seven of these individuals showed point heteroplasmy. The most frequent haplotype (16223T-16362C-73G-263G-315.1C) was found in 31 individuals and the second most frequent haplotypes (16129A-16223T- 16362C-73G-152C-263G-309.1C-315.1C) was found in 24 individuals. The haplotypes diversity and random match probability were calculated to be 0.9969 and 0.0040, respectively.

Mitochondrial DNA Transmission and Transcription After Somatic Cell Fusion to One or More Cytoplasts

Following fertilization, mitochondrial DNA is inherited from the oocyte and transmitted homoplasmically. However, following nuclear transfer, mitochondrial DNA can be transmitted from both the donor cell and recipient oocyte, resulting in a state of heteroplasmy. To determine whether the genetic diversity between donor cell and recipient cytoplast mitochondrial DNA influences development, we generated bovine embryos by fusing a donor cell to one or more enucleated cytoplasts. Analysis of mitochondrial DNA from embryos, fetal tissues, and blood samples from offspring revealed that early preimplantation embryos from two or three cytoplasts had significantly more mitochondrial DNA variants than fetal tissues. Phylogenic analysis of embryos generated using single cytoplasts divided the mitochondrial DNA sequence variants into three separate groups with various amounts of genetic divergence from the donor cell line. In heteroplasmic tissue and blood samples, the predominant mitochondrial DNA population was significantly more divergent from the donor cell than the less frequent allele. Furthermore, analysis of the mitochondrially encoded cytochrome B gene showed that two heteroplasmic alleles encoded for different amino acids, and the ratios of mitochondrial DNA/mRNA for each allele differed significantly between tissues. The degree of evolutionary distance between the donor cell and the cytoplast and the variability in heteroplasmy between tissues may have an impact on more divergent intergeneric nuclear transfer and the use of this approach for the generation of embryonic stem cells.

Mitochondrial Resequencing Arrays Detect Tumor-Specific Mutations in Salivary Rinses of Patients with Head and Neck Cancer

Alterations of the mitochondrial genome have been identified in multiple solid tumors and in many head and neck squamous cell carcinomas (HNSCC). Identification of mitochondrial mutations in the salivary rinses of patients with HNSCC has potential application in disease detection. In this study, we used the MitoChip v2.0 mitochondrial genome resequencing array to detect minor populations of mitochondrial DNA in salivary rinses of patients with HNSCC. Salivary rinses from 13 patients with HNSCC, whose tumors carried mitochondrial mutations, were collected before surgical resection. DNA isolated from salivary rinses and serial dilutions of DNA derived from HNSCC-derived cell lines with known mitochondrial mutations were sequenced using the MitoChip, and analyzed using a quantitative algorithm which we developed to detect minor populations of mitochondrial DNA from MitoChip probe intensity data. We detected heteroplasmic populations of mitochondrial DNA up to a 1:200 dilution using MitoChip v2.0 and our analysis algorithm. A logarithmic relationship between the magnitude of assay intensity and concentration of minor mitochondrial populations was shown. This technique was able to identify tumor-specific mitochondrial mutations in salivary rinses from 10 of 13 (76.9%) patients with head and neck cancer. Minor populations of mitochondrial DNA and disease-specific mitochondrial mutations in salivary rinses of patients with HNSCC can be successfully identified using the MitoChip resequencing array and the algorithm which we have developed. This technique has potential application in the surveillance of patients after resection and may have applicability in the surveillance of body fluids in other tumor types.

Use of Y Chromosome and Mitochondrial DNA Population Structure in Tracing Human Migrations

Well-resolved molecular gene trees illustrate the concept of descent with modification and exhibit the opposing processes of drift and migration, both of which influence population structure. Phylogenies of the maternally inherited mtDNA genome and the paternally inherited portion of the nonrecombining Y chromosome retain sequential records of the accumulation of genetic diversity. Although knowledge regarding the diversity of the entire human genome will be needed to completely characterize human genetic evolution, these uniparentally inherited loci are unique indicators of gender in modulating the extant population structure. We compare and contrast these loci for patterns of continuity and discreteness and discuss how their phylogenetic diversity and progression provide means to disentangle ancient colonization events by pioneering migrants from subsequent overlying migrations. We introduce new results concerning Y chromosome founder haplogroups C, DE, and F that resolve their previous trifurcation and improve the harmony with the mtDNA recapitulation of the out-of-Africa migration.

Mitochondrial and chloroplast phylogeography of Picea crassifolia Kom. (Pinaceae) in the Qinghai‐Tibetan Plateau and adjacent highlands

The disjunct distribution of forests in the Qinghai-Tibetan Plateau (QTP) and adjacent Helan Shan and Daqing Shan highlands provides an excellent model to examine vegetation shifts, glacial refugia and gene flow of key species in this complex landscape region in response to past climatic oscillations and human disturbance. In this study, we examined maternally inherited mitochondrial DNA (nad1 intron b/c and nad5 intron 1) and paternally inherited chloroplast DNA (trnC-trnD) sequence variation within a dominant forest species, Picea crassifolia Kom. We recovered nine mitotypes and two chlorotypes in a survey of 442 individuals from 32 populations sampled throughout the species' range. Significant mitochondrial DNA population subdivision was detected (G(ST) = 0.512; N(ST) = 0.679), suggesting low levels of recurrent gene flow through seeds among populations and significant phylogeographical structure (N(ST) > GST, P < 0.05). Plateau haplotypes differed in sequence from those in the adjacent highlands, suggesting a long period of allopatric fragmentation between the species in the two regions and the presence of independent refugia in each region during Quaternary glaciations. On the QTP platform, all but one of the disjunct populations surveyed were fixed for the same mitotype, while most populations at the plateau edge contained more than one haplotype with the mitotype that was fixed in plateau platform populations always present at high frequency. This distribution pattern suggests that present-day disjunct populations on the QTP platform experienced a common recolonization history. The same phylogeographical pattern, however, was not detected for paternally inherited chloroplast DNA haplotypes. Two chlorotypes were distributed throughout the range of the species with little geographical population differentiation (G(ST) = N(ST) = 0.093). This provides evidence for highly efficient pollen-mediated gene flow among isolated forest patches, both within and between the QTP and adjacent highland populations. A lack of isolation to pollen-mediated gene flow between forests on the QTP and adjacent highlands is surprising given that the Tengger Desert has been a geographical barrier between these two regions for approximately the last 1.8 million years.

Mitochondrial DNA diversity and population differentiation in southern East Asia

Mitochondrial DNA (mtDNA) polymorphism has been studied systematically in the Han, Tibeto-Buman, and Hmong-Mien ethnic families of southern East Asia. Only two families in this region, Daic and Austro-Asiatic, were still uninvestigated. Daic is a major ethnic family in South China and Southeast Asia and has a long history. To study mtDNA polymorphism within this family, all the Daic populations of China and some of Vietnam (774 individuals from 30 populations) were typed by HVS-1 region sequencing and by PCR-RFLP assays. The observed high Southern type frequencies (B, F, M7, R) confirmed Daic as a typical Southern group. mtDNAs of other populations (126 individuals from 14 populations) from Austro-Asiatic ethnic families neighboring the Daic were also typed. Networks of mtDNA haplogroups in South China were traced from these new data and those from the literature. Ethnic families share many haplogroups, indicating their common origin. However, the two largest families in South China, Daic, and Hmong-Mien, polarized into several ethnic family specific haplogroups. Haplogroup ages were estimated in the networks of high-frequency haplogroups (B, F, M7, R), and they were found to originate about 50,000 years ago. In contrast, ethnic family specific haplogroups all originated around 20,000 years ago. We therefore conclude that modern humans have lived in South China for a long time, inside-ethnogenesis was a rather late event, and frequent inmixing was taking place throughout. MtDNA data of Daic, Austro-Asiatic and other populations in South China has therefore proven pivotal for studying the human history of East Asia.

Mitochondrial DNA control region population data from Macedonia

Mitochondrial DNA sequences of the entire control region were analyzed in 200 unrelated individuals from Macedonia. A total of 163 different haplotypes were found as determined by 177 polymorphic sites. The probability of a random match was calculated as 1:121 (0.83%). The basic phylogenetic structure of the Macedonian population as derived from its haplogroup distribution is in agreement with other West-Eurasian populations. Upon publication, the population data are going to be available in the EMPOP database ( www.empop.org) [W. Parson, A. Dür, EMPOP—a forensic mtDNA database, FSI:Genetics 1 (2) (2007) 88–92; W. Parson, A. Brandstätter, A. Alonso, N. Brandt, B. Brinkmann, A. Carracedo, et al., The EDNAP mitochondrial DNA population database (EMPOP) collaborative exercises: organisation, results and perspectives, Forensic Sci. Int. 139 (2–3) (2004) 215–226.].

Mitochondrial DNA population data of HVS-I and HVS-II sequences from a northeast German sample

Mitochondrial DNA sequences of the control region's two hypervariable regions HVS-I and HVS-II were determined for 213 unrelated west Eurasian individuals from northeast Germany (Mecklenburg). A total of 174 different mtDNA haplotypes were found, 25 of which were shared by more than 1 individual. The most frequent haplotypes were 263G–309.1C–315.1C, found in seven individuals, 263G–309.1C–309.2C–315.1C, found in six individuals and 263G–315.1C, found in five individuals. These sequences are also the most common haplotypes in other published European data sets. The sequence polymorphisms consisting of 150 polymorphic nucleotide positions were compared with other European databases. The genetic diversity and random match probability were calculated. Our results corroborate certain features which are characteristic for west Eurasian mtDNA population samples.

Comment on "Population Size Does Not Influence Mitochondrial Genetic Diversity in Animals"

Bazin et al. (Reports, 28 April, 2006, p. 570) found no relationship between mitochondrial DNA (mtDNA) diversity and population size when comparing across large groups of animals. We show empirically that species with smaller populations, as represented by eutherian mammals, exhibit a positive correlation between mtDNA and allozyme variation, suggesting that mtDNA diversity may correlate with population size in these animals.

Mitochondrial DNA variation and population structure of the yarkand hareLepus yarkandensis

The yarkand hareLepus yarkandensis Gunther, 1875 is one of the endemic mammals in China. We used mitochondrial control region (D-loop) sequences to investigate the genetic variation, population structure and phylogeography of the Yarkand hare from Xinjiang Uighur Autonomous Region. The Yarkand hare holds a high level of haplotype diversity in most populations. In contrast, we found a relatively low level of gene diversity in Niya and Shule populations. The pairwise mismatch distribution suggested a pattern of population expansion. Fu’s test of neutrality was significant for the total population (F s=−10.458), which also suggested a recent population expansion. The pattern of diversification among populations was not consistent with its geographic distributions. Significant divergence was found among some populations, however, no significant divergence among other populations. The results indicate some degree of dispersal of individuals may exist among populations in recent history. Due to past isolation and habitat fragmentation, genetic variation was reduced within populations and division among populations increased. In order to conserve genetic diversity of the Yarkand hare in isolated populations, in the future, it may be advisable to translocate individuals between isolated populations to maintain current levels of genetic variation.

Decoupled mitochondrial and chloroplast DNA population structure reveals Holocene collapse and population isolation in a threatened Mexican‐endemic conifer

Chihuahua spruce (Picea chihuahuana Martínez) is a montane subtropical conifer endemic to the Sierra Madre Occidental in northwestern México. Range-wide variation was investigated using maternally inherited mitochondrial (mtDNA) and paternally inherited chloroplast (cpDNA) DNA markers. Among the 16 mtDNA regions analysed, only two mitotypes were detected, while the study of six cpDNA microsatellite markers revealed eight different chlorotypes. The average cpDNA diversity (H = 0.415) was low but much higher than that for mtDNA (H = 0). The distribution of mitotypes revealed two clear nonoverlapping areas (G(ST) = N(ST) = 1), one including northern populations and the second one including the southern and central stands, suggesting that these two regions may represent different ancestral populations. The cpDNA markers showed lower population differentiation (G(ST) = 0.362; R(ST) = 0.230), implying that the two ancestral populations continued to exchange pollen after their initial geographic separation. A lack of a phylogeographic structure was revealed by different spatial analyses of cpDNA (G(ST) > R(ST); and samova), and reduced cpDNA gene flow was noted among populations (Nm = 0.873). Some stands deviated significantly from the mutation-drift equilibrium, suggesting recent bottlenecks. Altogether, these various trends are consistent with the hypothesis of a population collapse during the Holocene warming and suggest that most of the modern P. chihuahuana populations are now effectively isolated with their genetic diversity essentially modelled by genetic drift. The conservation efforts should focus on most southern populations and on the northern and central stands exhibiting high levels of genetic diversity. Additional mtDNA sequence analysis confirmed that P. martinezii (Patterson) is not conspecific with P. chihuahuana, and thus deserves separate conservation efforts.

Mitochondrial DNA diversity and population structure of four Chinese donkey breeds

Mitochondrial DNA diversity and population structure of four Chinese donkey breeds

Heteroplasmy as a common state of mitochondrial genetic information in plants and animals

Plant and animal mitochondrial genomes, although quite distinct in size, structure, expression and evolutionary dynamics both may exhibit the state of heteroplasmy--the presence of more than one type of mitochondrial genome in an organism. This review is focused on heteroplasmy in plants, but we also highlight the most striking similarities and differences between plant and animal heteroplasmy. First we summarize the information on heteroplasmy generation and methods of its detection. Then we describe examples of quantitative changes in heteroplasmic populations of mitochondrial DNA (mtDNA) and consequences of such events. We also summarize the current knowledge about transmission and somatic segregation of heteroplasmy in plants and animals. Finally, factors which influence the stoichiometry of heteroplasmic mtDNA variants are discussed. Despite the apparent differences between the plant and animal heteroplasmy, the observed similarities allow one to conclude that this condition must play an important role in the mitochondrial biology of living organisms.

Genetic structure of Brazilian populations of white mouth croaker Micropogonias furnieri (Perciformes : Sciaenidae)

Mitochondrial (mt) DNA population structure was examined by using restriction fragment length polymorphism analysis (RFLP) among 149 white mouth croakers (Micropogonias furnieri) sampled from six locations in the Brazilian coast. Heterogeneity tests revealed no differentiation in mtDNA haplotype frequencies within the region between 23°S and 34°S (P=0.263), but significant heterogeneity was detected between north and south of 23°S (P&lt;0.003). Analysis of molecular variance (AMOVA) indicated a low F ST value (0.008, P=0.180) among south-central localities (23-34°S), but a significant level of population subdivision (F ST = 0.322, P&lt; 0.003) between north and south-central localities. UPGMA analysis of mtDNA sequence divergence revealed differentiation between white mouth croakers collected from north and south of 23°S. Mantel test found significant association between matrices of genetic distance and geographic distance among samples. Collectively, these data were consistent with a single genetic stock of Micropogonias furnieri in the Brazilian coast, with semi-isolated populations occurring in the north and south of 23° S.

Mitochondrial DNA variation and population structure of the Japanese mitten crab Eriocheir japonica

The Japanese mitten crab Eriocheir japonica is a common grapsid species found throughout freshwater and estuarine regions in Japan. In order to obtain information on the genetic variation and population structure of this species, a polymerase chain reaction (PCR) restriction fragment length polymorphism (RFLP) analysis was conducted on the cytochrome oxidase subunit I (CCI) of mitochondrial DNA, on 666 individuals from 19 sample sites covering the three main geographic regions of Japan (Main Islands, Okinawa, and Ogasawara). Genetic analysis using seven restriction enzymes produced an array of 61 composite haplotypes. Three regional groups corresponding to the three geographic regions were clearly identified by cluster and molecular variance model (AMOVA) analyses. Each of the three groups showed dominant haplotypes that were almost completely absent in populations from the other geographic areas. Comparison with published information for other species indicates that the degree of genetic divergence between these three main groups is equivalent to the genetic distance between congeneric species. Thus, the population structure of the Japanese mitten crab, as inferred from mtDNA analysis, is formed by genetically distinct groups that closely reflect their geographic distribution in the Japanese archipelago as well as restricted gene flow.

A microarray system for genotyping 150 single nucleotide polymorphisms in the coding region of human mitochondrial DNA

We established a genotyping system for a panel of 150 SNPs in the coding regions of mitochondrial DNA based on multiplex tag-array minisequencing. We show the feasibility of this system for simultaneous identification of individuals and prediction of the geographical origin of the mitochondrial DNA population lineage of the sample donors by genotyping the panel of SNPs in 265 samples representing nine different populations from Africa, Europe, and Asia. Nearly 40,000 genotypes were produced in the study, with an overall genotyping success rate of 95% and accuracy close to 100%. The gene diversity value of the panel of 150 SNPs was 0.991, compared to 0.995 for sequencing 500 nucleotides of the hypervariable regions I and II of mtDNA. For 17 individuals with identical sequences in the hypervariable regions of mtDNA, our panel of SNPs increased the power of discrimination. We observed 144 haplotypes that correspond to previously determined mitochondrial “haplogroups,” and they allowed prediction of the origin of the maternal population lineage of 97% of the analyzed samples.

Mitochondrial DNA variation and population structure of the island endemic Azorean bat (Nyctalus azoreum).

The Azorean bat Nyctalus azoreum is the only endemic mammal native to the remote archipelago of the Azores. It evolved from a continental ancestor related to the Leisler's bat Nyctalus leisleri and is considered threatened because of its restricted and highly fragmented distribution. We studied the genetic variability in 159 individuals from 14 colonies sampled throughout the archipelago. Sequences of the D-loop region revealed moderate but highly structured genetic variability. Half of the 15 distinct haplotypes were restricted to a single island, but the most common was found throughout the archipelago, suggesting a single colonization event followed by limited interisland female gene flow. All N. azoreum haplotypes were closely related and formed a star-like structure typical of expanded populations. The inferred age of demographic expansions was consistent with the arrival of founder animals during the Holocene, well before the first humans inhabited the Azores. Comparisons with a population of N. leisleri from continental Portugal confirmed not only that all N. azoreum lineages were unique to the archipelago, but also that the current levels of genetic diversity were surprisingly high for an insular species. Our data imply that the Azorean bat has a high conservation value. We argue that geographical patterns of genetic structuring indicate the existence of two management units.