Inter- and intragenerational transmission of a human mitochondrial DNA heteroplasmy among 13 maternally-related individuals and differences between and within tissues in two family members

The transmission of a C16,291C/T heteroplasmy in the HV1 region of human mitochondrial DNA (mtDNA) was examined in buccal cells from 13 maternally-related individuals across three generations and in additional tissues (hair, blood, or finger nails) from three members of this family. The ratio of C:T at nucleotide position (np) 16,291 showed wide intra- and intergenerational variation as well as tissue variation within individuals. Our results demonstrate that one or two sequence differences between samples in the mtDNA does not warrant an exclusion. To avoid false exclusions especially when comparing mtDNA from hair samples, we recommend the analysis of as many samples as possible in order to minimize the possibility that the detection of a rare polymorphism in a single sample would be considered an exclusion when it is really a match. The observation that the transmission of a mtDNA heteroplasmy from one individual to her offspring is likely to differ among the first-generation offspring and between that generation and subsequent generations lends further credence to the bottleneck theory of inheritance of human mtDNA.     

Human ribosomal RNA variants from a single individual and their expression in different tissues.

We have investigated the extent of sequence variation in human ribosomal RNA (rRNA) genes and the expression of specific rRNA gene variants in different tissues of an individual. Focusing on the fifth variable region (V5; nt 2065-2244) of the 28S rRNA gene, we find that sequence differences between rRNA genes of a single individual are characterized by differences in number of repeats of simple sequences at four specific sites. These data support and extend previous findings which show similar V5 sequence variation in rRNA genes from a group of individuals. We performed experiments to determine if there is differential gene expression within the rRNA multigene family. From the analysis of data of six variant V5 probes protected from RNase digestion by rRNAs isolated from different tissues of the individual, we conclude that each variant rRNA is present in a similar proportion in these tissues, whereas the actual contributions of variants differ, their relative proportion is maintained from tissue to tissue in an individual. We favor the explanation of a gene dosage effect over that of a regulated gene effect to account for this pattern of rRNA gene expression. In addition, computer generated secondary structure models of each V5 clone structure predict the same three helix structure with the regions of sequence variation contained in one stem-loop structure.     

Direct selection of mutations in the human mitochondrial tRNAThr gene: reversion of an 'uncloneable' phenotype

Several regions of the human mitochondrial genome are refractory to cloning in plasmid and bacteriophage DNA vectors. For example, recovery of recombinant M13 clones containing a 462 basepair MboI-Kpn I restriction fragment that spans nucleotide positions 15591 to 16053 of HeLa cell mitochondrial DNA was as much as 100-fold lower than the recovery of M13 clones containing other regions of the human mitochondrial genome. All of 50 recombinant M13 clones containing this 'uncloneable' fragment had one or more changes in nucleotide sequence. Each clone contained at least one alteration in two nucleotide positions within the tRNAThr gene that encode portions of the anticodon loop and D-stem of the HeLa mitochondrial tRNAThr. These results imply that the HeLa mitochondrial tRNAThr gene is responsible for the 'uncloneable' phenotype of this region of human mitochondrial (mt) DNA. A total of 61 nucleotide sequence alterations were identified in 50 independent clones containing the HeLa mt tRNAThr gene. 56 mutations were single-base substitutions; 5 were deletions. Approximately 80% of the base substitution mutations were A:T----G:C transitions. A preference for A:T----G:C transition mutations also characterizes polymorphic base substitution variants in the mitochondrial DNA of unrelated individuals. This similarity suggests that human mitochondrial DNA sequence variation within and between individuals may have a common origin.     

Mitochondrial DNA Variation and the Evolution of Robertsonian Chromosomal Races of House Mice, Mus Domesticus

The house mouse, Mus domesticus, includes many distinct Robertsonian (Rb) chromosomal races with diploid numbers from 2n = 22 to 2n = 38. Although these races are highly differentiated karyotypically, they are otherwise indistinguishable from standard karyotype (i.e., 2n = 40) mice, and consequently their evolutionary histories are not well understood. We have examined mitochondrial DNA (mtDNA) sequence variation from the control region and the ND3 gene region among 56 M. domesticus from Western Europe, including 15 Rb populations and 13 standard karyotype populations, and two individuals of the sister species, Mus musculus. mtDNA exhibited an average sequence divergence of 0.84% within M. domesticus and 3.4% between M. domesticus and M. musculus. The transition/transversion bias for the regions sequenced is 5.7:1, and the overall rate of sequence evolution is approximately 10% divergence per million years. The amount of mtDNA variation was as great among different Rb races as among different populations of standard karyotype mice, suggesting that different Rb races do not derive from a single recent maternal lineage. Phylogenetic analysis of the mtDNA sequences resulted in a parsimony tree which contained six major clades. Each of these clades contained both Rb and standard karyotype mice, consistent with the hypothesis that Rb races have arisen independently multiple times. Discordance between phylogeny and geography was attributable to ancestral polymorphism as a consequence of the recent colonization of Western Europe by mice. Two major mtDNA lineages were geographically localized and contained both Rb and standard karyotype mice. The age of these lineages suggests that mice have moved into Europe only within the last 10,000 years and that Rb populations in different geographic regions arose during this time.     

Direct selection of mutations in the human mitochondrial tRNAThr gene: reversion of an ‘uncloneable’ phenotype

Several regions of the human mitochondrial genome are refractory to cloning in plasmid and bacteriophage DNA vectors. For example, recovery of recombinant M13 clones containing a 462 basepair MboI-Kpn I restriction fragment that spans nucleotide positions 15591 to 16053 of HeLa cell mitochondrial DNA was as much as 100-fold lower than the recovery of M13 clones containing other regions of the human mitochondrial genome. All of 50 recombinant M13 clones containing this ‘uncloneable’ fragment had one or more changes in nucleotide sequence. Each clone contained at least one alteration in two nucleotide positions within the tRNA^Thr gene that encode portions of the anticodon loop and D-stem of the HeLa mitochondrial tRNA^Thr. These results imply that the HeLa mitochondrial tRNA^Thr gene is responsible for the ‘uncloneable’ phenotype of this region of human mitochondrial (mt) DNA.A total of 61 nucleotide sequence alterations were identified in 50 independent clones containing the HeLa mt tRNA^Thr gene. 56 mutations were single-base substitutions; 5 were deletions. Approximately 80% of the base substitution mutations were A:T → G:C transitions. A preference for A:T → G:C transition mutations also characterizes polymorphic base substitution variants in the mitochondrial DNA of unrelated individuals. This similarity suggests that human mitochondrial DNA sequence variation within and between individuals may have a common origin.     

小麂、赤麂、黑麂mtDNA序列变异性及反映的进化关系

利用已测定的鹿科麂亚科动物小麂、赤麂、黑麂的线粒体全基因组序列,统计它们各自连接在一起的13个蛋白编码基因、22个tRNA基因、2个rRNA基因和1个控制区序列的碱基长度和组成,计算rRNA基因遗传距离,估算分歧时间,比较蛋白编码基因的碱基水平和氨基酸水平上的差异,基于连接在一起的13个氨基酸序列,以羊为外群,通过邻位相连法和最大简约性法构建进化树,探讨小麂、赤麂、黑麂的进化关系。结果表明,小麂是较原始的物种,赤麂和黑麂较为近缘,是从类似小麂的祖先演化而来。     

Length heteroplasmy in the first hypervariable segment of the human mtDNA control region.

The first hypervariable segment of the human mtDNA control region contains a homopolymeric tract of cytosines between nt 16184 and 16193, interrupted at position 16189 by a thymine, according to the Cambridge reference sequence. A variant commonly found in population screening is a T-to-C transition at nt 16189, resulting in an uninterrupted homopolymeric tract. Direct sequencing of individuals with this variant produces a characteristic blurred sequence in nucleotides beyond the tract. Sequencing clones from these individuals revealed that this is caused by high levels of length heteroplasmy in the homopolymeric tract and low levels of length heteroplasmy in the four adenines following the tract. We have developed a rapid method involving densitometry of sequencing gels to quantify the relative proportions of different length variants present in an individual. We have used this to study the proportions of length variants in individuals from three twin pairs and two maternal lineages. While unrelated individuals usually have different proportions of length variants, all maternally related individuals studied have the same proportions, even if they are only distantly related. It is not obvious how identical heteroplasmic profiles are maintained in maternally related individuals, but some possible mechanisms are suggested.     

Recurrent Tissue-Specific mtDNA Mutations Are Common in Humans

Mitochondrial DNA (mtDNA) variation can affect phenotypic variation; therefore, knowing its distribution within and among individuals is of importance to understanding many human diseases. Intra-individual mtDNA variation (heteroplasmy) has been generally assumed to be random. We used massively parallel sequencing to assess heteroplasmy across ten tissues and demonstrate that in unrelated individuals there are tissue-specific, recurrent mutations. Certain tissues, notably kidney, liver and skeletal muscle, displayed the identical recurrent mutations that were undetectable in other tissues in the same individuals. Using RFLP analyses we validated one of the tissue-specific mutations in the two sequenced individuals and replicated the patterns in two additional individuals. These recurrent mutations all occur within or in very close proximity to sites that regulate mtDNA replication, strongly implying that these variations alter the replication dynamics of the mutated mtDNA genome. These recurrent variants are all independent of each other and do not occur in the mtDNA coding regions. The most parsimonious explanation of the data is that these frequently repeated mutations experience tissue-specific positive selection, probably through replication advantage.     

Denaturing gradient gel method for mapping single base changes in human mitochondrial DNA.

A denaturing gradient gel electrophoresis (DGGE) method is described that detects even single base pair changes in mitochondrial DNA (mtDNA). In this method, restriction fragments of mtDNA are electrophoresed in a urea/formamide gradient gel at 60 degrees C. Migration distance of each mtDNA fragment in the gel depends on melting behavior which reflects base composition. Fragments are located by Southern blotting with specific mtDNA probes. With just four carefully chosen restriction enzymes and as little as 50-100 ng of mtDNA, the method covers almost the entire human mitochondrial genome. To demonstrate the method, human mtDNA was analyzed. In six normal individuals, DGGE revealed melting behavior polymorphisms (MBPs) in mtDNA fragments that were not detected by restriction fragment length polymorphism (RFLP) analysis in agarose gels. Another individual, shown to have a melting behavior polymorphism in the cytochrome b coding region, was studied in detail. By mapping, the mutation was deduced to lie between nt 14905 and 15370. The affected fragment was amplified by PCR and sequenced. Specific base changes were identified in the region predicted by the gel result. This method will be especially useful as a diagnostic tool in mitochondrial disease for rapid localization of mtDNA mutations to specific regions of the genome, but DGGE also could complement RFLP analysis as a more sensitive method to follow maternal lineage in human and animal populations in a variety of research fields.     

Length Mutations in Human Mitochondrial DNA

By high-resolution, restriction mapping of mitochondrial DNAs purified from 112 human individuals, we have identified 14 length variants caused by small additions and deletions (from about 6 to 14 base pairs in length). Three of the 14 length differences are due to mutations at two locations within the D loop, whereas the remaining 11 occur at seven sites that are probably within other noncoding sequences and at junctions between coding sequences. In five of the nine regions of length polymorphism, there is a sequence of five cytosines in a row, this sequence being comparatively rare in coding DNA. Phylogenetic analysis indicates that, in most of the polymorphic regions, a given length mutation has arisen several times independently in different human lineages. The average rate at which length mutations have been arising and surviving in the human species is estimated to be many times higher for noncoding mtDNA than for noncoding nuclear DNA. The mystery of why vertebrate mtDNA is more prone than nuclear DNA to evolve by point mutation is now compounded by the discovery of a similar bias toward rapid evolution by length mutation.     

On the estimation of the rate of nucleotide substitution for the control region of human mitochondrial DNA

To apply molecular clock for studying human evolution, the pattern of nucleotide substitution for the control region of human mtDNA was analyzed in detail. It is well known that the rate of nucleotide substitution for the control region is much higher than that for any other part of mtDNA. In this study, the higher substitution rate was attributed to the higher rate of transition-type substitution between pyrimidines within the D-loop part, whereas the rates of other types of substitution were essentially the same over the entire mtDNA molecule. Even within the control region, the rate and pattern of nucleotide substitution were different between the D-loop part and the rest. The rate and pattern for the non-D-loop part were very similar to those for fourfold-degenerate sites in the protein-coding region. In contrast, the D-loop and non-D-loop parts showed similarities in the base composition, whereas the base composition of fourfold-degenerate sites slightly different from that of the both parts of the control region. It is concluded, therefore, that the nucleotide frequencies of the control region should be used to estimate the number of substitutions (d) between the control region sequences. However, a method to verify the accuracy of the estimation of d by means of the transition/transversion (s/v) ratio was theoretically studied. It was suggested that the s/v ratio becomes constant over a wide range of d values only when the estimation of d is unbiased. On the basis of this result, the estimates of d previously obtained between human sequences were evaluated.     

Homopolymeric tract heteroplasmy in mtDNA from tissues and single oocytes: support for a genetic bottleneck.

While mtDNA polymorphisms at single base positions are common, the overwhelming majority of the mitochondrial genomes within a single individual are usually identical. When there is a point-mutation difference between a mother and her offspring, there may be a complete switching of mtDNA type within a single generation. It is generally assumed that there is a genetic bottleneck whereby a single or small number of founder mtDNA(s) populate the organism, but it is not known at which stages the restriction/amplification of mtDNA subtype(s) occur, and this uncertainty impedes antenatal diagnosis for mtDNA disorders. Length polymorphisms in homopolymeric tracts have been demonstrated in the large noncoding region of mtDNA. We have developed a new method, T-PCR (trimmed PCR), to quantitate heteroplasmy for two of these tracts (D310 and D16189). D310 variation is sufficient to indicate clonal origins of tissues and single oocytes. Tissues from normal individuals often possessed more than one length variant (heteroplasmy). However, there was no difference in the pattern of the length variants between somatic tissues in any control individual when bulk samples were taken. Oocytes from normal women undergoing in vitro fertilization were frequently heteroplasmic for length variants, and in two cases the modal length of the D310 tract differed in individual oocytes from the same woman. These data suggest that a restriction/amplification event, which we attribute to clonal expansion of founder mtDNA(s), has occurred by the time oocytes are mature, although further segregation may occur at a later stage. In contrast to controls, the length distribution of the D310 tract varied between tissues in a patient with heteroplasmic mtDNA rearrangements, suggesting that these mutants influence segregation. These findings have important implications for the genetic counselling of patients with pathogenic mtDNA mutations.     

The Use of Restriction Endonucleases to Measure Mitochondrial DNA Sequence Relatedness in Natural Populations. I. Population Structure and Evolution in the Genus Peromyscus

In this study we introduce to natural population analysis a molecular technique that involves the use of restriction endonucleases to compare mitochondrial DNA (mtDNA) sequences. We have examined the fragment patterns produced by six restriction endonucleases acting upon mtDNA isolated from 23 samples of three species of the rodent Peromyscus. Our observations confirm the following conclusions derived from previous experiments with laboratory animals: (1) mtDNA within an individual homogeneous; (2) at least the majority of mtDNA present in an individual is inherited from the female parent. Our experiments demonstrate for the first time that there is detectable heterogeneity in mtDNA sequences within and among natural geographic populations of a species and that this heterogeneity can readily be used to estimate relatedness between individuals and populations. Individuals collected within a single locale show less than 0.5% sequence divergence, while those collected from conspecific populations separated by 50 ti 500 miles differ by approximately 1.5%. The mtDNAs of the closely related sibling species P. polionotus and P. maniculatus differ from each other by 13 to 17%; nonsibling species differ by more than 20%. Qualitative and quantitative approaches to analysis of digestion patterns are suggested. The results indicate that restriction analysis of mtNDA may become the most sensitive and powerful technique yet available for reconstructing evolutionary relationships among conspecific organisms.     

Sequence variations of mitochondrial DNA D-loop region are highly frequent events in familial breast cancer

Mitochondrial DNA (mtDNA) is known for its high frequencies of polymorphisms and mutations. The non-coding displacement (D)-loop, especially a mononucleotide repeat (poly-C) between 303 and 315 nucleotides (D310), has been recently identified as a frequent hotspot of mutations in human neoplasia, including breast cancer. To further explore the sequence variations of mitochondrial D-loop region in familial breast cancer and their possible associations with breast cancer risk, PCR-SSCP and direct DNA sequencing methods were used to detect the variants of the mtDNA D-Loop in 23 familial breast cancer patients as well as three high-risk cancer families. Compared to that in sporadic breast tumors (53.3%, 16/30) and healthy blood donors (6.7%, 2/30), we identified a total of 126 sequence alterations in 23/23 (100%) of familial breast cancer patients, including eight novel nucleotide variants. Among these changes, A to G at nt.263, T to C at nt.489, T to C at nt.310, TC insertion at nt.311, CA deletion at nt.522, and C to G at nt.527 were highly frequent ones. In addition, among three high-risk cancer families, we found that individuals affected with breast cancer harbored more mtDNA sequence variants in mtDNA D310 area than other affected family members. Together, our data indicate that sequence variants within the mtDNA D-Loop region are frequent events in Chinese familial breast cancer patients. Some of these nucleotide abnormalities, particularly those in D310 segment, might be involved in the breast carcinogenesis and could be included in a panel of molecular biomarkers for cancer susceptibility early-detection strategy.     

A new disease-related mutation for mitochondrial encephalopathy lactic acidosis and strokelike episodes (MELAS) syndrome affects the ND4 subunit of the respiratory complex I.

The molecular lesions in two patients exhibiting classical clinical manifestations of MELAS (mitochondrial encephalopathy, lactic acidosis, and strokelike episodes) syndrome have been investigated. A recently reported disease-related A----G base substitution at nt 3243 of the mtDNA, in the DHU loop of tRNA(Leu), was detected by restriction-enzyme analysis of the relevant PCR-amplified segment of the mtDNA of one patient but was not observed, by either restriction-enzyme analysis or nucleotide sequencing, in the other. To define the molecular lesion in the patient who does not have the A----G base substitution at nt 3243, the total mitochondrial genome of the patient has been sequenced. An A----G base substitution at nt 11084, leading to a Thr-to-Ala amino acid replacement in the ND4 subunit of the respiratory complex I, is suggested to be a disease-related mutation.     

Mutational mechanisms by which an inactive replication origin of bacteriophage M13 is turned on are similar to mechanisms of activation of ras proto-oncogenes.

M13 viral strand synthesis is initiated by nicking of the viral strand of the duplex replicative form by the M13 gene II initiator protein at a specific site within a sequence of about 40 base pairs having dyad symmetry. Efficient replication of the M13 viral strand also requires the presence of an adjacent sequence of ca. 100 base pairs. Together these sequences constitute the minimal origin for M13 viral strand synthesis. A pBR322 derivative having a 182-base-pair insert of M13 DNA contains a functional M13 viral strand origin and, when provided with M13 gene functions in trans, replicates under conditions nonpermissive for the parent plasmid. Chimeric plasmids containing deletions within the sequence flanking the viral strand origin are unable to replicate under these conditions. We isolated spontaneous mutants of M13 based on their ability to activate replication of such plasmids. The mutations found in these strains, as well as several produced by oligonucleotide-directed mutagenesis, all result in the substitution of any of at least four different amino acids for a specific glycine residue near the amino-terminal end of the initiator protein. Other studies have shown that overproduction of the wild-type initiator protein also restores replication. These alternate mechanisms are discussed in terms of their striking similarity to the mechanisms of activation of the ras proto-oncogenes which can be activated either by increased expression of the wild-type protein or by substitution of any of several amino acids for a glycine residue near the amino terminus.     

Evidence for mitochondrial DNA recombination in a human population of island Melanesia

Mitochondrial DNA (mtDNA) analysis has proved useful in studies of recent human evolution and the genetic affinities of human groups of different geographical regions. As part of an extensive survey of mtDNA diversity in present-day Pacific populations, we obtained sequence information of the hypervariable mtDNA control region of 452 individuals from various localities in the western Pacific. The mtDNA types fell into three major groups which reflect the settlement history of the area. Interestingly, we detected an extremely rare point mutation at high frequency in the small island of Nguna in the Melanesian archipelago of Vanuatu. Phylogenetic analysis of the mtDNA data indicated that the mutation was present in individuals of separate mtDNA lineages. We propose that the multiple occurrence of a rare mutation event in one isolated locality is highly improbable, and that recombination between different mtDNA types is a more likely explanation for our observation. If correct, this conclusion has important implications for the use of mtDNA in phylogenetic and evolutionary studies.     

Detection and quantification of mitochondrial DNA deletions in individual cells by real-time PCR

Defects of mitochondrial DNA (mtDNA) are an important cause of disease and play a role in the ageing process. There are multiple copies of the mitochondrial genome in a single cell. In many patients with acquired or inherited mtDNA mutations, there exists a mixture of mutated and wild type genomes (termed heteroplasmy) within individual cells. As a biochemical and clinical defect is only observed when there are high levels of mutated mtDNA, a crucial investigation is to determine the level of heteroplasmic mutations within tissues and individual cells. We have developed an assay to determine the relative amount of deleted mtDNA using real-time fluorescence PCR. This assay detects the vast majority of deleted molecules, thus eliminating the need to develop specific probes. We have demonstrated an excellent correlation with other techniques (Southern blotting and three- primer competitive PCR), and have shown this technique to be sensitive to quantify the level of deleted mtDNA molecules in individual cells. Finally, we have used this assay to investigate patients with mitochondrial disease and shown in individual skeletal muscle fibres that there exist different patterns of abnormalities between patients with single or multiple mtDNA deletions. We believe that this technique has significant advantages over other methods to quantify deleted mtDNA and, employed alongside our method to sequence the mitochondrial genome from single cells, will further our understanding of the role of mtDNA mutations in human disease and ageing.     

Population variation of human mtDNA control region sequences detected by enzymatic amplification and sequence-specific oligonucleotide probes.

A method for detecting sequence variation of hypervariable segments of the mtDNA control region was developed. The technique uses hybridization of sequence-specific oligonucleotide (SSO) probes to DNA sequences that have been amplified by PCR. The nucleotide sequences of the two hypervariable segments of the mtDNA control region from 52 individuals were determined; these sequences were then used to define nine regions suitable for SSO typing. A total of 23 SSO probes were used to detect sequence variants at these nine regions in 525 individuals from five ethnic groups (African, Asian, Caucasian, Japanese, and Mexican). The SSO typing revealed an enormous amount of variability, with 274 mtDNA types observed among these 525 individuals and with diversity values, for each population, exceeding .95. For each of the nine mtDNA regions significant differences in the frequencies of sequence variants were observed between these five populations. The mtDNA SSO-typing system was successfully applied to a case involving individual identification of skeletal remains; the probability of a random match was approximately 0.7%. The potential useful applications of this mtDNA SSO-typing system thus include the analysis of individual identity as well as population genetic studies.