Lack of Structural Variation but Extensive Length Polymorphisms and Heteroplasmic Length Variations in the Mitochondrial DNA Control Region of Highly Inbred Crested Ibis,Nipponia nippon

The animal mitochondrial DNA (mtDNA) length polymorphism and heteroplasmy are accepted to be universal. Here we report the lack of structural variation but the presence of length polymorphism as well as heteroplasmy in mtDNA control region of an endangered avian species – the Crested Ibis (Nipponia nippon). The complete control region was directly sequenced while the distribution pattern and inheritance of the length variations were examined using both direct sequencing and genotyping of the PCR fragments from captive birds with pedigrees, wild birds and a historical specimen. Our results demonstrated that there was no structural variation in the control region, however, different numbers of short tandem repeats with an identical motif of CA₃CA₂CA₃ at the 3′-end of the control region determined the length polymorphisms among and heteroplasmy within individual birds. There were one to three predominant fragments in every bird; nevertheless multiple minor fragments coexist in all birds. These extremely high polymorphisms were suggested to have derived from the ‘replication slippage’ of a perfect microsatellite evolution following the step-wise mutational model. The patterns of heteroplasmy were found to be shifted between generations and among siblings but rather stable between blood and feather samples. This study provides the first evidence of a very extensive mtDNA length polymorphism and heteroplasmy in the highly inbred Crested Ibis which carries an mtDNA genome lack of structural genetic diversity. The analysis of pedigreed samples also sheds light on the transmission of mtDNA length heteroplasmy in birds following the genetic bottleneck theory. Further research focusing on the generation and transmission of particular mtDNA heteroplasmy patterns in single germ line of Crested Ibis is encouraged by this study.     

Analysis of heteroplasmy in the major noncoding region of mitochondrial DNA in the blood and atherosclerotic plaques of carotid arteries

For identification of somatic mitochondrial DNA (mtDNA) mutations, the mtDNA major noncoding region (D-loop) sequence in blood samples and carotid atherosclerosis plaques from patients with atherosclerosis was analyzed. Five point heteroplasmic positions were observed in 4 of 23 individuals (17%). Only in two cases could heteroplasmy have resulted from somatic mutation, whereas three heteroplasmic positions were found in both vascular tissue and blood. In addition, length heteroplasmy in a polycytosine stretches was registered at nucleotide positions 303–315 in 16 individuals, and also in the 16184–16193 region in four patients. The results suggest that somatic mtDNA mutations can occur during atherosclerosis, but some heteroplasmic mutations may appear in all tissues, possibly being inherited.     

Frequency and Pattern of Heteroplasmy in the Control Region of Human Mitochondrial DNA

In this work, we present the results of the screening of human mitochondrial DNA (mtDNA) heteroplasmy in the control region of mtDNA from 210 unrelated Spanish individuals. Both hypervariable regions of mtDNA were amplified and sequenced in order to identify and quantify point and length heteroplasmy. Of the 210 individuals analyzed, 30% were fully homoplasmic and the remaining presented point and/or length heteroplasmy. The prevalent form of heteroplasmy was length heteroplasmy in the poly(C) tract of the hypervariable region II (HVRII), followed by length heteroplasmy in the poly(C) tract of hypervariable region I (HVRI) and, finally, point heteroplasmy, which was found in 3.81% of the individuals analyzed. Moreover, no significant differences were found in the proportions of the different kinds of heteroplasmy in the population when blood and buccal cell samples were compared. The pattern of heteroplasmy in HVRI and HVRII presents important differences. Moreover, the mutational profile in heteroplasmy seems to be different from the mutational pattern detected in population. The results suggest that a considerable number of mutations and, particularly, transitions that appear in heteroplasmy are probably eliminated by drift and/or by selection acting at different mtDNA levels of organization. Taking as a whole the results reported in this work, it is mandatory to perform a broad-scale screening of heteroplasmy to better establish the heteroplasmy profile which would be important for medical, evolutionary, and forensic proposes.     

MtDNA substitution rate and segregation of heteroplasmy in coding and noncoding regions

The mitochondrial DNA (mtDNA) substitution rate and segregation of heteroplasmy were studied for the non-coding control region (D-loop) and 500 bp of the coding region between nucleotide positions 5550 and 6050, by sequence analysis of blood samples from 194 individuals, representing 33 maternal lineages. No homoplasmic nucleotide substitutions were detected in a total of 292 transmissions. The estimated substitution rate per nucleotide per million years for the control region (micro>0.21, 95% CI 0-0.6) was not significantly different from that for the coding region (micro>0.54, 95% CI 0-1.0). Variation in the length of homopolymeric C streches was observed at three sites in the control region (positions 65, 309 and 16,189), all of which were in the heteroplasmic state. Segregation of heteroplasmic genotypes between generations was observed in several maternal pedigrees. At position 309, a longer poly C tract length was strongly associated with a higher probability for heteroplasmy and rapid segregation between generations. The length heteroplasmy at positions 65 and 16,189 was found at low frequency and was confined to a few families.     

Polymorphism of mitochondrial genome noncoding regions in the three Kazakh populations inhabited different areas of Kazakhstan and in the samples of DNA from ancient people of Kazakhstan Altai

Polymorphism of major noncoding region of mitochondrial DNA (mtDNA D-loop, 528 bp in length) from the three modem kazakh populations and from DNA samples of ancient people inhabited modern Kazakhstani Altai were studied. PCR and RFLP analysis of 13 sites of restriction--BamHI, EcoRV, Sau3AI (1 restriction site), KpnI (2 sites), HaeIII (3 sites), RsaI (5 restriction sites), were carried out. The distribution of each site frequencies was determined. Nucleotide diversity (h) and genetic distance between different kazakh population and other populations of world were estimated. The same RFLP analysis of the mitochondrial DNA control region was carried out for the paleogenomic samples. It was shown that two samples of ancient mitochondrial DNA were monomorphous throughout all analyzed restriction sites.     

Mitochondrial genome size variation in New World and Old World populations of Drosophila melanogaster

Drosophila melanogaster originated in Africa, spread to Europe and Asia, and is believed to have colonized the New World in the past few hundred years. Levels of genetic variation are typically reduced in New World populations, consistent with a founder event following range expansion out of Africa and the Old World. We describe the patterns of mtDNA length variation within and among several populations of Drosophila melanogaster from the Old and New World. MtDNA length variation is due to insertion and deletion of tandem repeats in the control region (D-loop) of D. melanogaster mitochondrial genome. The distinct mutational dynamics of this system provide an opportunity to compare the patterns of variation in this marker to those of other markers with different mutational pressures and linkage relationships. The data show significantly more length variation in African and Asian samples than in New World samples. New World samples also show more pronounced skew of the length distribution. Our results are distinct from an earlier study that showed significantly higher levels of length variation and heteroplasmy. The level of heteroplasmy is highly correlated with the number of years that samples have been maintained in laboratory culture, suggesting that relaxed selection in small populations permits the accumulation of mtDNA length variation and heteroplasmy. Together, the data indicate that mtDNA length variants retain a signature of founder events and selection, and suggest that further investigation into the mutation-selection dynamics of the D-loop region of mtDNA would provide a distinct and informative marker for analysis of the recent history of populations.     

Heteroplasmy of Short Tandem Repeats in Mitochondrial DNA of Atlantic Cod, Gadus Morhua

The mitochondrial DNA of the Atlantic cod (Gadus morhua) contains a tandem array of 40-bp repeats in the D-loop region of the molecule. Variation among molecules in the copy number of these repeats results in mtDNA length variation and heteroplasmy (the presence of more than one form of mtDNA in an individual). In a sample of fish collected from different localities around Iceland and off George's Bank, each individual was heteroplasmic for two or more mtDNAs ranging in repeat copy number from two (common) to six (rare). An earlier report on mtDNA heteroplasmy in sturgeon (Acipenser transmontanus) presented a competitive displacement model for length mutations in mtDNAs containing tandem arrays and the cod data deviate from this model. Depending on the nature of putative secondary structures and the location of D-loop strand termination, additional mechanisms of length mutation may be needed to explain the range of mtDNA length variants maintained in these populations. The balance between genetic drift and mutation in maintaining this length polymorphism is estimated through a hierarchical analysis of diversity of mtDNA length variation in the Iceland samples. Eighty percent of the diversity lies within individuals, 8% among individuals and 12% among localities. An estimate of theta = 2N(eo) mu greater than 1 indicates that this system is characterized by a high mutation rate and is governed primarily by deterministic dynamics. The sequences of repeat arrays from fish collected in Norway, Iceland and George's Bank show no nucleotide variation suggesting that there is very little substructuring to the North Atlantic cod population.     

Selection against pathogenic mtDNA mutations in a stem cell population leads to the loss of the 3243A-->G mutation in blood

The mutation 3243A-->G is the most common heteroplasmic pathogenic mitochondrial DNA (mtDNA) mutation in humans, but it is not understood why the proportion of this mutation decreases in blood during life. Changing levels of mtDNA heteroplasmy are fundamentally related to the pathophysiology of the mitochondrial disease and correlate with clinical progression. To understand this process, we simulated the segregation of mtDNA in hematopoietic stem cells and leukocyte precursors. Our observations show that the percentage of mutant mtDNA in blood decreases exponentially over time. This is consistent with the existence of a selective process acting at the stem cell level and explains why the level of mutant mtDNA in blood is almost invariably lower than in nondividing (postmitotic) tissues such as skeletal muscle. By using this approach, we derived a formula from human data to correct for the change in heteroplasmy over time. A comparison of age-corrected blood heteroplasmy levels with skeletal muscle, an embryologically distinct postmitotic tissue, provides independent confirmation of the model. These findings indicate that selection against pathogenic mtDNA mutations occurs in a stem cell population.     

Simultaneous A8344G heteroplasmy and mitochondrial DNA copy number quantification in myoclonus epilepsy and ragged-red fibers (MERRF) syndrome by a multiplex molecular beacon based real-time fluorescence PCR

The association of a particular mitochondrial DNA (mtDNA) mutation with different clinical phenotypes is a well-known feature of mitochondrial diseases. A simple genotype–phenotype correlation has not been found between mutation load and disease expression. Tissue and intercellular mosaicism as well as mtDNA copy number are thought to be responsible for the different clinical phenotypes. As disease expression of mitochondrial tRNA mutations is mostly in postmitotic tissues, studies to elucidate disease mechanisms need to be performed on patient material. Heteroplasmy quantitation and copy number estimation using small patient biopsy samples has not been reported before, mainly due to technical restrictions. In order to resolve this problem, we have developed a robust assay that utilizes Molecular Beacons to accurately quantify heteroplasmy levels and determine mtDNA copy number in small samples carrying the A8344G tRNA^Lys mutation. It provides the methodological basis to investigate the role of heteroplasmy and mtDNA copy number in determining the clinical phenotypes.     

The use of next generation sequencing technology to study the effect of radiation therapy on mitochondrial DNA mutation

The human mitochondrial genome has an exclusively maternal mode of inheritance. Mitochondrial DNA (mtDNA) is particularly vulnerable to environmental insults due in part to an underdeveloped DNA repair system, limited to base excision and homologous recombination repair. Radiation exposure to the ovaries may cause mtDNA mutations in oocytes, which may in turn be transmitted to offspring. We hypothesized that the children of female cancer survivors who received radiation therapy may have an increased rate of mtDNA heteroplasmy mutations, which conceivably could increase their risk of developing cancer and other diseases. We evaluated 44 DNA blood samples from 17 Danish and 1 Finnish families (18 mothers and 26 children). All mothers had been treated for cancer as children and radiation doses to their ovaries were determined based on medical records and computational models. DNA samples were sequenced for the entire mitochondrial genome using the Illumina GAII system. Mother's age at sample collection was positively correlated with mtDNA heteroplasmy mutations. There was evidence of heteroplasmy inheritance in that 9 of the 18 families had at least one child who inherited at least one heteroplasmy site from his or her mother. No significant difference in single nucleotide polymorphisms between mother and offspring, however, was observed. Radiation therapy dose to ovaries also was not significantly associated with the heteroplasmy mutation rate among mothers and children. No evidence was found that radiotherapy for pediatric cancer is associated with the mitochondrial genome mutation rate in female cancer survivors and their children.     

Genetic status of the wood stork (Mycteria americana) from the southeastern United States and the Brazilian Pantanal as revealed by mitochondrial DNA analysis

The wood stork (Mycteria americana) is a colonial wading bird that inhabits the Neotropical region from the southeastern United States (US) to northern Argentina. The species is considered to be endangered in the US due to degradation of its foraging and breeding habitat. In other parts of its range, such as in the Brazilian Pantanal region, breeding populations of this species appear to be stable. We compared the levels of genetic variability and population structuring of the US and the Pantanal breeding populations using mitochondrial DNA (mtDNA) control region sequences. Twenty-seven haplotypes were identified among 88 wood stork samples collected from eight breeding colonies in the US and eight in the Pantanal. Patterns indicative of heteroplasmy were observed in 35.3% of the mtDNA sequences that were examined. Significantly higher levels of haplotype diversity were observed in the Pantanal samples compared to those from the US, suggesting that during the last century, demographic declines or a recent evolutionary bottleneck reduced the levels of mtDNA variability of the US population. Analyses of genetic structuring revealed non-significant genetic differentiation between these regions, indicating that either the populations were only recently separated or that gene flow continues to occur at low levels. Haplotype network analysis indicated low current levels of gene flow between populations that were closely related in the past.     

Variation of the mitochondrial genome in the evolution of Drosophila

The evidence on mitochondrial genome variation and its role in evolution of the genus Drosophila are reviewed. The mitochondrial genome is represented by a circular double-stranded DNA molecule 16 to 19 kb in length. The genome contains no introns involved in recombination. The entire mitochondrial genome can be arbitrarily divided into three parts: (1) protein-coding genes; (2) genes encoding rRNA and tRNA; and (3) the noncoding regulatory region (A + T region). The selective importance of mutations within different mtDNA regions is therefore unequal. In Drosophila, the content of the A + T pairs in mtDNA is extremely low and a pattern of nucleotide substitution is characterized by a low transition/transversion ratio (and a low threshold of mutation saturation). The deletions and duplications are of common occurrence in the mitochondrial genome. However, this genome lacks such characteristic for the nuclear genome aberrations as the inversions and transpositions. The phenomena of introgression and heteroplasmy provide an opportunity to study the adaptive role of the mitochondrial genome and its role in speciation. Analysis of evidence concerning mtDNA variation in different species of the genus Drosophila made it possible to ascertain data on phylogenetic relationships among species obtained by studying nuclear genome variation. In some species, mtDNA variation may serve as a reliable marker for population differentiation within a species, although evidence on the population dynamics of the mtDNA variation is very scarce.     

Polymorphism of the Noncoding Region of the Mitochondrial Genome in Three Kazakh Populations Inhabiting Different Areas of Kazakhstan and in DNA Samples from Ancient Populations of the Kazakhstani Altai

The polymorphism of the major noncoding region of mitochondrial DNA (mtDNA D loop, 528 bp) has been studied in samples from three modern Kazakh populations (from Almaty, the Semipalatinsk Region, and the Altai Mountains) and in DNA samples of ancient human populations of the Kazakhstani Altai. Polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) analysis for 13 restriction sites, including BamHI, EcoRV, Sau3AI (one site each), KpnI (two sites), HaeIII (three sites), and RsaI (five sites) were used. The frequency distributions of all sites have been determined. The gene diversity (h) and the genetic distances between different Kazakh populations and other populations of the world have been calculated. The RFLP analysis of the mtDNA control region of fossil samples has been performed similarly to the analysis of modern mtDNA samples. Two fossil mtDNA samples from burial mound 11 are monomorphic with respect to all restriction sites analyzed.     

线粒体DNA的异质性与检测方法

线粒体DNA的异质性在生物学、医学、法医、生物技术等领域有重要的应用。本文从自然发生(包括体细胞突变、父系渗入和杂交)和人工产生(包括转基因、细胞融合、核移植和转线粒体)两大方面系统地介绍了线粒体DNA异质性的产生机制及其遗传。并介绍了线粒体DNA异质性的检测方法,已知突变位点mtDNA异质性的检测方法有原位PCR技术、PCR-RFLP和实时荧光定量PCR技术,未知突变位点mtDNA异质性的检测方法有长PCR技术、时相温度梯度凝胶电泳(TTGE)和变性高效液相色谱(DHPLC)等。最后对克隆生物中线粒体异质性检测的应用实例作了介绍。     

Variation of the Mitochondral Genome in the Evolution of Drosophila

The evidence on mitochondrial genome variation and its role in evolution of the genus Drosophila are reviewed. The mitochondrial genome is represented by a circular double-stranded DNA molecule 16 to 19 kb in length. Mitochondrial genes lack introns and recombination. The entire mitochondrial genome can be arbitrarily divided into three parts: (1) protein-coding genes; (2) genes encoding rRNA and tRNA; and (3) the noncoding regulatory region (A + T region). The selective importance of mutations within different mtDNA regions is therefore unequal. In Drosophila, the content of the A + T pairs in mtDNA is extremely high and a pattern of nucleotide substitution is characterized by a low transition/transversion ratio (and a low threshold of mutation saturation). The deletions and duplications are of common occurrence in the mitochondrial genome. However, this genome lacks such characteristic for the nuclear genome aberrations as inversions and transpositions. The phenomena of introgression and heteroplasmy provide an opportunity to study the adaptive role of the mitochondrial genome and its role in speciation. Analysis of evidence concerning mtDNA variation in different species of the genus Drosophilamade it possible to ascertain data on phylogenetic relationships among species obtained by studying nuclear genome variation. In some species, mtDNA variation may serve as a reliable marker for population differentiation within a species, although evidence on the population dynamics of the mtDNA variation is very scarce.     

Postglacial population expansion of Japanese macaques (<Emphasis Type="Italic">Macaca fuscata</Emphasis>) inferred from mitochondrial DNA phylogeography

We investigated the diversity and phylogeography of mitochondrial DNA (mtDNA) in Japanese macaques (Macaca fuscata), an endemic species in Japan that has the northernmost distribution of any non-human primate species. DNA samples from 135 localities representing the entire range of this species were compared. A total of 53 unique haplotypes were observed for the 412-bp partial mtDNA control region sequence, with length variation distinguishing the two subspecies. Clustering analyses suggested two putative major haplogroups, of which one was geographically distributed in eastern Japan and the other in western Japan. The populations in the east showed lower mtDNA diversity than those in the west. Phylogeographical relationships of haplotypes depicted with minimum spanning network suggested differences in population structure. Population expansion was significant for the eastern but not the western population, suggesting establishment of the ancestral population was relatively long ago in the west and recent in the east. Based on fossil evidence and past climate and vegetation changes, we inferred that the postulated population expansion may have taken place after the last glacial period (after 15,000 years ago). Mitochondrial DNA showed contrasting results in both variability and phylogenetic status of local populations to those of previous studies using protein variations, particularly for populations in the periphery of the range, with special inference on habitat change during the glacial period in response to cold adaptation. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Stock Structure and Homing Fidelity in Gulf of Mexico Sturgeon (Acipenser Oxyrinchus Desotoi) Based on Restriction Fragment Length Polymorphism and Sequence Analyses of Mitochondrial DNA

Efforts have been proposed worldwide to restore sturgeon populations through the use of hatcheries to supplement natural reproduction and to reintroduce sturgeon where they have become extinct. We examined the population structure and inferred the extent of homing in the anadromous Gulf of Mexico (Gulf) sturgeon (Acipenser oxyrinchus desotoi). Restriction fragment length polymorphism and control region sequence analyses of mitochondrial DNA (mtDNA) were used to identify haplotypes of Gulf sturgeon specimens obtained from eight drainages spanning the subspecies' entire distribution from Louisiana to Florida. Significant differences in haplotype frequencies indicated substantial geographic structuring of populations. A minimum of four regional or river-specific populations were identified (from west to east): (1) Pearl River, LA and Pascagoula River, MS, (2) Escambia and Yellow rivers, FL, (3) Choctawhatchee River, FL, and (4) Apalachicola, Ochlockonee, and Suwannee rivers, FL. Estimates of maternally mediated gene flow between any pair of the four regional or river-specific stocks ranged between 0.15 to 1.2. Tandem repeats in the mtDNA control region of Gulf sturgeon were not perfectly conserved. This result, together with an absence of heteroplasmy and length variation in Gulf sturgeon mtDNA, indicates that the molecular mechanisms of mtDNA control region sequence evolution differ among acipenserids.     

Length Variation, Heteroplasmy and Sequence Divergence in the Mitochondrial DNA of Four Species of Sturgeon (Acipenser)

The extent of mtDNA length variation and heteroplasmy as well as DNA sequences of the control region and two tRNA genes were determined for four North American sturgeon species: Acipenser transmontanus, A. medirostris, A. fulvescens and A. oxyrhnychus. Across the Continental Divide, a division in the occurrence of length variation and heteroplasmy was observed that was concordant with species biogeography as well as with phylogenies inferred from restriction fragment length polymorphisms (RFLP) of whole mtDNA and pairwise comparisons of unique sequences of the control region. In all species, mtDNA length variation was due to repeated arrays of 78-82-bp sequences each containing a D-loop strand synthesis termination associated sequence (TAS). Individual repeats showed greater sequence conservation within individuals and species rather than between species, which is suggestive of concerted evolution. Differences in the frequencies of multiple copy genomes and heteroplasmy among the four species may be ascribed to differences in the rates of recurrent mutation. A mechanism that may offset the high rate of mutation for increased copy number is suggested on the basis that an increase in the number of functional TAS motifs might reduce the frequency of successfully initiated H-strand replications.     

A twin study of mitochondrial DNA polymorphisms shows that heteroplasmy at multiple sites is associated with mtDNA variant 16093 but not with zygosity

The mitochondrial theory of ageing proposes that damage to mitochondria and diminished mitochondrial DNA (mtDNA) repair are major contributors to cellular dysfunction and age-related diseases. We investigate the prevalence of heteroplasmy in the mtDNA control region in buccal swab and blood derived samples for 178 women from the TwinsUK cohort (41 DZ pair 39 MZ pairs, 18 singletons, mean age 57.5 range 28–82) and its relationship to age, BMI and fasting insulin and glucose serum levels. The overall estimated prevalence of heteroplasmy for both tissues in the control region measured for 37 sites was 17%. The prevalence of heteroplasmy was higher among the older half of the study subjects than in the younger half (23% vs 10% p<0.03), primarily reflecting the increase in the prevalence of a heteroplasmic dinucleotide CA repeat in variable region II (VRII) with age. The VRII 523–524 heteroplasmic site (heteroplasmic in 25 subjects) was also associated with a decrease in BMI. In addition, concordance rates for common heteroplasmy were observed to be near complete for both dizygotic (DZ = 94%) and monozygotic twin pairs (MZ = 100%), consistent with previous reports that suggest variation in heteroplasmy rates between generations are determined by bottlenecks in maternal transmission of mitochondria. Differences in the prevalence of heteroplasmy were observed overall between samples derived from buccal swabs (19%) and blood (15%, p<0.04). These were particularly marked at position 16093 of hypervariable region I (HVI, 7% vs 0%, respectively, p<4×10⁻¹¹). The presence of the C allele at position 16093 in blood was associated with the presence of heteroplasmy in buccal swabs at this position (p = 3.5×10⁻¹⁴) and also at VRII (p = 2×10⁻⁴) suggesting a possible predisposing role for this site in the accumulation of heteroplasmy. Our data indicate that BMI is potentially associated with control region heteroplasmy.