Mitochondrial DNA-like sequences in the human nuclear genome. Characterization and implications in the evolution of mitochondrial DNA

Human mtDNA-Encoded Long ncRNAs: Knotty Molecules and Complex Functions

Until a few decades ago, most of our knowledge of RNA transcription products was focused on protein-coding sequences, which were later determined to make up the smallest portion of the mammalian genome. Since 2002, we have learnt a great deal about the intriguing world of non-coding RNAs (ncRNAs), mainly due to the rapid development of bioinformatic tools and next-generation sequencing (NGS) platforms. Moreover, interest in non-human ncRNAs and their functions has increased as a result of these technologies and the accessibility of complete genome sequences of species ranging from Archaea to primates. Despite not producing proteins, ncRNAs constitute a vast family of RNA molecules that serve a number of regulatory roles and are essential for cellular physiology and pathology. This review focuses on a subgroup of human ncRNAs, namely mtDNA-encoded long non-coding RNAs (mt-lncRNAs), which are transcribed from the mitochondrial genome and whose disparate localisations and functions are linked as much to mitochondrial metabolism as to cellular physiology and pathology.

The Mighty NUMT: Mitochondrial DNA Flexing Its Code in the Nuclear Genome

Nuclear-mitochondrial DNA segments (NUMTs) are mitochondrial DNA (mtDNA) fragments that have been inserted into the nuclear genome. Some NUMTs are common within the human population but most NUMTs are rare and specific to individuals. NUMTs range in size from 24 base pairs to encompassing nearly the entire mtDNA and are found throughout the nuclear genome. Emerging evidence suggests that the formation of NUMTs is an ongoing process in humans. NUMTs contaminate sequencing results of the mtDNA by introducing false positive variants, particularly heteroplasmic variants present at a low variant allele frequency (VAF). In our review, we discuss the prevalence of NUMTs in the human population, the potential mechanisms of de novo NUMT insertion via DNA repair mechanisms, and provide an overview of the existing approaches for minimizing NUMT contamination. Apart from filtering known NUMTs, both wet lab-based and computational methods can be used to minimize the contamination of NUMTs in analyses of human mtDNA. Current approaches include: (1) isolating mitochondria to enrich for mtDNA; (2) applying basic local alignment to identify NUMTs for subsequent filtering; (3) bioinformatic pipelines for NUMT detection; (4) k-mer-based NUMT detection; and (5) filtering candidate false positive variants by mtDNA copy number, VAF, or sequence quality score. Multiple approaches must be applied in order to effectively identify NUMTs in samples. Although next-generation sequencing is revolutionizing our understanding of heteroplasmic mtDNA, it also raises new challenges with the high prevalence and individual-specific NUMTs that need to be handled with care in studies of mitochondrial genetics.

Insertions of mitochondrial DNA into the nucleus—effects and role in cell evolution

We review the insertion of mitochondrial DNA (mtDNA) fragments into nuclear DNA (NUMTS) as a general and ongoing process that has occurred many times during genome evolution. Fragments of mtDNA are generated during the lifetime of organisms in both somatic and germinal cells, by the production of reactive oxygen species in the mitochondria. The fragments are inserted into the nucleus during the double-strand breaks repair via the non-homologous end-joining machinery, followed by genomic instability, giving rise to the high variability observed in NUMT patterns among species, populations, or genotypes. Some de novo produced mtDNA insertions show harmful effects, being involved in human diseases, carcinogenesis, and ageing. NUMT generation is a non-stop process overpassing the Mendelian transmission. This parasitic property ensures their survival even against their harmful effects. The accumulation of mtDNA fragments mainly at pericentromeric and subtelomeric regions is important to understand the transmission and integration of NUMTs into the genomes. The possible effect of female meiotic drive for mtDNA insertions at centromeres remains to be studied. In spite of the harmful feature of NUMTs, they are important in cell evolution, representing a major source of genomic variation.

Numts help to reconstruct the demographic history of the ocellated lizard (Lacerta lepida) in a secondary contact zone

In northwestern Iberia, two largely allopatric Lacerta lepida mitochondrial lineages occur, L5 occurring to the south of Douro River and L3 to the north, with a zone of putative secondary contact in the region of the Douro River valley. Cytochrome b sequence chromatograms with polymorphisms at nucleotide sites diagnostic for the two lineages were detected in individuals in the region of the Douro River and further north within the range of L3. We show that these polymorphisms are caused by the presence of four different numts (I-IV) co-occurring with the L3 genome, together with low levels of heteroplasmy. Two of the numts (I and II) are similar to the mitochondrial genome of L5 but are quite divergent from the mitochondrial genome of L3 where they occur. We show that these numts are derived from the mitochondrial genome of L5 and were incorporated in L3 through hybridization at the time of secondary contact between the lineages. The additional incidence of these numts to the north of the putative contact zone is consistent with an earlier postglacial northward range expansion of L5, preceding that of L3. We show that genetic exchange between the lineages responsible for the origin of these numts in L3 after secondary contact occurred prior to, or coincident with, the northward expansion of L3. This study shows that, in the context of phylogeographic analysis, numts can provide evidence for past demographic events and can be useful tools for the reconstruction of complex evolutionary histories.

Time and time again: unisexual salamanders (genus Ambystoma) are the oldest unisexual vertebrates

Background

The age of unisexual salamanders of the genus Ambystoma is contentious. Recent and ancient evolutionary histories of unisexual Ambystoma were proposed by a few separate studies that constructed phylogenies using mitochondrial DNA markers (cytochrome b gene vs. non-coding region). In contrast to other studies showing that unisexual Ambystoma represent the most ancient unisexual vertebrates, a recent study by Robertson et al. suggests that this lineage has a very recent origin of less than 25,000 years ago.

Results

We re-examined the phylogenetic relationship of the unisexuals to A. barbouri from various populations using both mitochondrial markers as well as the complete mitochondrial genomes of A. barbouri and a unisexual individual from Kentucky. Lineage dating was conducted using BEAST and MultiDivTime on a complete mitochondrial genome phylogeny. Our results support a monophyletic lineage for unisexual Ambystoma that shares its most recent common ancestor with an A. barbouri lineage from western Kentucky. In contrast to the Robertson et al.'s study, no A. barbouri individual shared an identical or almost identical cytochrome b haplotype with any unisexual. Molecular dating supports an early Pliocene origin for the unisexual linage (~5 million years ago). We propose that a unisexual-like cytochrome b numt (or pseudogene) exists in the controversial A. barbouri individuals from Kentucky, which was likely the cause of an erroneous phylogeny and time estimate in Robertson et al.'s study.

Conclusion

We reject a recent origin of unisexual Ambystoma and provide strong evidence that unisexual Ambystoma are the most ancient unisexual vertebrates known to exist. The likely presence of an ancient cytochrome b numt in some Kentucky A. barbouri represents a molecular "fossil" reinforcing the hypothesis that these individuals are some of the closest extant relatives to unisexual Ambystoma.

Beyond linear sequence comparisons: the use of genome-level characters for phylogenetic reconstruction

The first whole genomes to be compared for phylogenetic inference were those of mitochondria, which provided the first sets of genome-level characters for phylogenetic reconstruction. Most powerful among these characters has been the comparisons of the relative arrangements of genes, which has convincingly resolved numerous branch points, including those that had remained recalcitrant even to very large molecular sequence comparisons. Now the world faces a tsunami of complete nuclear genome sequences. In addition to the tremendous amount of DNA sequence that is becoming available for comparison, there is also a potential for many more genome-level characters to be developed, including the relative positions of introns, the domain structures of proteins, gene family membership, the presence of particular biochemical pathways, aspects of DNA replication or transcription, and many others. These characters can be especially convincing owing to their low likelihood of reverting to a primitive condition or occurring independently in separate lineages, thereby reducing the occurrence of homoplasy. The comparisons of organelle genomes pioneered the way for using such features for phylogenetic reconstructions, and it is almost certainly true, as ever more genomic sequence becomes available, that further use of genome-level characters will play a big role in outlining the relationships among major animal groups.

Organization and evolution of mitochondrial gene clusters in human

Currently, the spatial patterns of mitochondrial genes and how the genomic localization of (pseudo)genes originated from mitochondrial DNA remain largely unexplained. The aim of this study was to elucidate the organization of mitochondrial (pseudo)genes given their evolutionary origin. We used a keyword finding method and a bootstrapping method to estimate parameter values that represent the distribution pattern of mitochondrial genes in the nuclear genome. Almost half of mitochondrial genes showing physical clusters were located in the pericentromeric and subtelomeric regions of the chromosome. Most interestingly, the size of these clusters ranged from 0.085 to 3.2 Mb (average+/-SD 1.3+/-0.73 Mb), which coincides with the size of the evolutionary pocket, or the average size of evolutionary breakpoint regions. Our findings imply that the localization of mitochondrial genes in the human genome is determined independent of adaptation.

Extensive mitochondrial DNA transfer in a rapidly evolving rodent has been mediated by independent insertion events and by duplications

Mitochondrial DNA translocations to the nucleus (numt pseudogenes) are pervasive among eukaryotes, but copy number within the nuclear genome varies widely among taxa. As an increasing number of genomes are sequenced in their entirety, the origins, transfer mechanisms and insertion sites of numts are slowly being characterized. We investigated mitochondrial transfers within a genetically diverse rodent lineage and here report 15 numts totaling 21.8 kb that are harbored within the nuclear genome of the vole Microtus rossiaemeridionalis. The 15 numts total 21.8 kb and range from 0.39 to over 3.0 kb in length. Phylogenetic analyses revealed that these numts resulted from three independent insertions to the nucleus, two of which were followed by subsequent nuclear duplication events. The dates of the two translocations that led to subsequent duplications were estimated at 1.97 and 1.19 MYA, which coincide with the origin and radiation of the genus Microtus. Numt sequence data from five Microtus species were used to estimate an average rate of nucleotide substitution as 2.6x10(-8) subs/site/yr. This substitution rate is higher than in many other mammals, but is concordant with the elevated rate of mtDNA substitution in this lineage. Our data suggest that numt translocation in Microtus is more extensive than in either Mus or in Rattus, consistent with the elevated rate of speciation, karyotypic rearrangement, and mitochondrial DNA evolution in Microtus.

The complete mitochondrial genome of Flustrellidra hispida and the phylogenetic position of Bryozoa among the Metazoa

The complete mitochondrial genome of Flustrellidra hispida (Bryozoa, Ctenostomata, Flustrellidridae) was sequenced using a transposon-mediated approach. All but one of the 36 genes were identified (trnS2). The genome is 13,026 bp long, being one of the smallest metazoan mitochondrial genomes sequenced to date with a unique gene order when compared to other Metazoa. The genome has an overall AT richness of 59.4%. We found seven regions of overlaps between tRNAs and protein-coding genes ranging from 2 to 11 nt, and seven regions of overlap between tRNAs, ranging from 1 to 8 nt, resulting in a total number of 46 overlapping nucleotides. Genes nad4, cox2, atp8, and nad3 are terminated by the abbreviated stop codon T and cytb is suggested to terminate on (ACT)AA; we postulate that mRNA editing is required to remove AC for TAA to be functional in terminating translation. Phylogenetic analysis of nucleotide and amino acid data place Flustrellidra in the Lophotrochozoa. DNA for this study originated from two populations resulting in a contig consisting of multiple haplotypes. Twenty-seven SNP sites were detected, the majority occurring in cox1 and nad5. With cox1 already established as a marker in bryozoan studies, we advocate the further testing of nad5.

Adaptive responses to mitochondrial dysfunction in the rho degrees Namalwa cell

Mitochondrial function depends on the synchronised expression of nuclear and mitochondrial-encoded genes. The expression of many nuclear genes is sensitive to the functional state of mitochondria. Mitochondria containing mutated mtDNA populations have been shown to withstand stress to a high level of tolerance. The prolonged survival of cells that are devoid of mitochondrial DNA is further evidence for the existence of this adaptive response. Yeast cells lacking the mitochondrial genome configure their pattern of nuclear gene expression to accommodate changes in nutrient availability and in response to different stress conditions. Little is known about the signalling pathways that mediate the adaptive response in mammalian cells with mitochondrial dysfunction. To gain a global perspective on adaptive responses in cells with mitochondrial dysfunction, we used Affymetrix microarray technology to compare the whole genome transcription profile of Namalwa cells (rho(+)) to that of a Namalwa cell depleted of mitochondrial DNA (rho degrees ). Our results demonstrate that rho degrees Namalwa induce genes encoding hypoxia target proteins, mitochondrial ribosomal proteins, transport channels, tRNA synthetases and enzymes of the glycolytic pathway and the fatty oxidation pathway suggesting that hypoxia and the ensuing cellular adaptations could play an important adaptive role in cells with mitochondrial dysfunction.

Forty million years of independent evolution: a mitochondrial gene and its corresponding nuclear pseudogene in primates

Sequences from nuclear mitochondrial pseudogenes (numts) that originated by transfer of genetic information from mitochondria to the nucleus offer a unique opportunity to compare different regimes of molecular evolution. Analyzing a 1621-nt-long numt of the rRNA specifying mitochondrial DNA residing on human chromosome 3 and its corresponding mitochondrial gene in 18 anthropoid primates, we were able to retrace about 40 MY of primate rDNA evolutionary history. The results illustrate strengths and weaknesses of mtDNA data sets in reconstructing and dating the phylogenetic history of primates. We were able to show the following. In contrast to numt-DNA, (1) the nucleotide composition of mtDNA changed dramatically in the different primate lineages. This is assumed to lead to significant misinterpretations of the mitochondrial evolutionary history. (2) Due to the nucleotide compositional plasticity of primate mtDNA, the phylogenetic reconstruction combining mitochondrial and nuclear sequences is unlikely to yield reliable information for either tree topologies or branch lengths. This is because a major part of the underlying sequence evolution model--the nucleotide composition--is undergoing dramatic change in different mitochondrial lineages. We propose that this problem is also expressed in the occasional unexpected long branches leading to the "common ancestor" of orthologous numt sequences of different primate taxa. (3) The heterogeneous and lineage-specific evolution of mitochondrial sequences in primates renders molecular dating based on primate mtDNA problematic, whereas the numt sequences provide a much more reliable base for dating.

Nuclear mitochondrial pseudogenes as molecular outgroups for phylogenetically isolated taxa: a case study in Sphenodon

'Living fossil' taxa, by definition, have no close relatives, and therefore no outgroup to provide a root to phylogenetic trees. We identify and use a molecular outgroup in the sole extant lineage of sphenodontid reptiles, which separated from other reptiles 230 million years ago. We isolated and sequenced a partial nuclear copy of the mitochondrial cytochrome b gene. We confirm the copy is indeed not mitochondrial, is older than all extant mitochondrial copies in Sphenodon (tuatara), and is therefore useful as a molecular outgroup. Under phylogenetic analysis, the nuclear copy places the root of the tuatara mitochondrial gene tree between the northern and the southern (Cook Strait) groups of islands of New Zealand that are the last refugia for Sphenodon. This analysis supports a previous mid-point rooted mitochondrial gene tree. The mitochondrial DNA tree conflicts with allozyme analyses which place a Cook Strait population equidistant to all northern and other Cook Strait populations. This population on North Brother Island is the only natural population of extant S. guntheri; thus, we suggest that the current species designations of tuatara require further investigation.

Molecular systematics and origin of sociality in mongooses (Herpestidae, Carnivora)

The Herpestidae are small terrestrial carnivores comprising 18 African and Asian genera, currently split into two subfamilies, the Herpestinae and the Galidiinae. The aim of this work was to resolve intra-familial relationships and to test the origin of sociality in the group. For this purpose we analysed sequences of the complete cytochrome b gene for 18 species of Herpestidae. The results showed that the mongooses were split into three clades: (1) the Malagasy taxa (Galidiinae and Cryptoprocta), (2) the true social mongooses and (3) the solitary mongooses, each group being also supported by morphological and chromosomal data. Our results suggested unexpected phylogenetic relationships: (1) the genus Cynictis is included in the solitary mongoose clade, (2) the genera Liberiictis and Mungos are sister-group, and (3) the genus Herpestes is polyphyletic. We examined the evolution of the sociality in mongooses by combining behavioural traits with the cytochrome b data. Some of the behavioural traits provided good synapomorphies for characterizing the social species clade, showing the potential benefit of using such characters in phylogeny. The mapping of ecological and behavioural features resulted in hypothesizing solitary behavior and life in forest as the conditions at the base of the mongoose clade.

Low number of mitochondrial pseudogenes in the chicken (Gallus gallus) nuclear genome: implications for molecular inference of population history and phylogenetics

Background

Mitochondrial DNA has been detected in the nuclear genome of eukaryotes as pseudogenes, or Numts. Human and plant genomes harbor a large number of Numts, some of which have high similarity to mitochondrial fragments and thus may have been inadvertently included in population genetic and phylogenetic studies using mitochondrial DNA. Birds have smaller genomes relative to mammals, and the genome-wide frequency and distribution of Numts is still unknown. The release of a preliminary version of the chicken (Gallus gallus) genome by the Genome Sequencing Center at Washington University, St. Louis provided an opportunity to search this first avian genome for the frequency and characteristics of Numts relative to those in human and plants.

Results

We detected at least 13 Numts in the chicken nuclear genome. Identities between Numts and mitochondrial sequences varied from 58.6 to 88.8%. Fragments ranged from 131 to 1,733 nucleotides, collectively representing only 0.00078% of the nuclear genome. Because fewer Numts were detected in the chicken nuclear genome, they do not represent all regions of the mitochondrial genome and are not widespread in all chromosomes. Nuclear integrations in chicken seem to occur by a DNA intermediate and in regions of low gene density, especially in macrochromosomes.

Conclusion

The number of Numts in chicken is low compared to those in human and plant genomes, and is within the range found for most sequenced eukaryotic genomes. For chicken, PCR amplifications of fragments of about 1.5 kilobases are highly likely to represent true mitochondrial amplification. Sequencing of these fragments should expose the presence of unusual features typical of pseudogenes, unless the nuclear integration is very recent and has not yet been mutated. Metabolic selection for compact genomes with reduced repetitive DNA and gene-poor regions where Numts occur may explain their low incidence in birds.

Pattern of organization of human mitochondrial pseudogenes in the nuclear genome

Mitochondrial pseudogenes in the human nuclear genome have been previously described, mostly as a source of artifacts during the analysis of the mitochondrial genome. With the availability of the complete human genome sequence, we performed a comprehensive analysis of mtDNA insertions into the nucleus. We found 612 independent integrations that are evenly distributed among all chromosomes as well as within each individual chromosome. The identified pseudogenes account for a content of at least 0.016% of the human nuclear DNA. Up to 30% of a chromosome's mtDNA pseudogene content is composed of fragments that encompass two or more adjacent mitochondrial genes, and we found no correlation between the abundance of mitochondrial transcripts and the multiplicity of integrations. These observations indicate that the migrations of mitochondrial DNA sequences to the nucleus were predominantly DNA mediated. Phylogenetic analysis of the mtDNA pseudogenes and mtDNA sequences of primates indicate a continuous transfer into the nucleus. Because of the limited window of opportunity for mtDNA transfer to the germline, sperm mtDNA, which is released from degenerating mitochondria after fertilization, could be an important source of nuclear mtDNA pseudogenes.

Using secondary structure to identify ribosomal numts: cautionary examples from the human genome

The identification of inadvertently sequenced mitochondrial pseudogenes (numts) is critical to any study employing mitochondrial DNA sequence data. Failure to discriminate numts correctly can confound phylogenetic reconstruction and studies of molecular evolution. This is especially problematic for ribosomal mtDNA genes. Unlike protein-coding loci, whose pseudogenes tend to accumulate diagnostic frameshift or premature stop mutations, functional ribosomal genes are not constrained to maintain a reading frame and can accumulate insertion-deletion events of varying length, particularly in nonpairing regions. Several authors have advocated using structural features of the transcribed rRNA molecule to differentiate functional mitochondrial rRNA genes from their nuclear paralogs. We explored this approach using the mitochondrial 12S rRNA gene and three known 12S numts from the human genome in the context of anthropoid phylogeny and the inferred secondary structure of primate 12S rRNA. Contrary to expectation, each of the three human numts exhibits striking concordance with secondary structure models, with little, if any, indication of their pseudogene status, and would likely escape detection based on structural criteria alone. Furthermore, we show that the unwitting inclusion of a particularly ancient (18-25 Myr old) and surprisingly cryptic human numt in a phylogenetic analysis would yield a well-supported but dramatically incorrect conclusion regarding anthropoid relationships. Though we endorse the use of secondary structure models for inferring positional homology wholeheartedly, we caution against reliance on structural criteria for the discrimination of rRNA numts, given the potential fallibility of this approach.

Mitochondrial pseudogenes are pervasive and often insidious in the snapping shrimp genus Alpheus

Here we show that multiple DNA sequences, similar to the mitochondrial cytochrome oxidase I (COI) gene, occur within single individuals in at least 10 species of the snapping shrimp genus Alpheus. Cloning of amplified products revealed the presence of copies that differed in length and (more frequently) in base substitutions. Although multiple copies were amplified in individual shrimp from total genomic DNA (gDNA), only one sequence was amplified from cDNA. These results are best explained by the presence of nonfunctional duplications of a portion of the mtDNA, probably located in the nuclear genome, since transfer into the nuclear gene would render the COI gene nonfunctional due to differences in the nuclear and mitochondrial genetic codes. Analysis of codon variation suggests that there have been 21 independent transfer events in the 10 species examined. Within a single animal, differences between the sequences of these pseudogenes ranged from 0.2% to 20.6%, and those between the real mtDNA and pseudogene sequences ranged from 0.2% to 18.8% (uncorrected). The large number of integration events and the large range of divergences between pseudogenes and mtDNA sequences suggest that genetic material has been repeatedly transferred from the mtDNA to the nuclear genome of snapping shrimp. Unrecognized pseudogenes in phylogenetic or population studies may result in spurious results, although previous estimates of rates of molecular evolution based on Alpheus sister taxa separated by the Isthmus of Panama appear to remain valid. Especially worrisome for researchers are those pseudogenes that are not obviously recognizable as such. An effective solution may be to amplify transcribed copies of protein-coding mitochondrial genes from cDNA rather than using genomic DNA.

Mitochondrial pseudogenes: evolution's misplaced witnesses

Nuclear copies of mitochondrial DNA (mtDNA) have contaminated PCR-based mitochondrial studies of over 64 different animal species. Since the last review of these nuclear mitochondrial pseudogenes (Numts) in animals, Numts have been found in 53 of the species studied. The recent evidence suggests that Numts are not equally abundant in all species, for example they are more common in plants than in animals, and also more numerous in humans than in Drosophila. Methods for avoiding Numts have now been tested, and several recent studies demonstrate the potential utility of Numt DNA sequences in evolutionary studies. As relics of ancient mtDNA, these pseudogenes can be used to infer ancestral states or root mitochondrial phylogenies. Where they are numerous and selectively unconstrained, Numts are ideal for the study of spontaneous mutation in nuclear genomes.

Recombination among multiple mitochondrial pseudogenes from a Passerine genus

PCR products of a fragment of the mitchondrial protein coding subunit 5 of NADH-dehydrogenase (ND5) from eight individuals representing five species of the South American bird genus Conirostrum were cloned. The 130 clones, which were subsequently sequenced, constituted 55 different sequences. Due to the observed differences in substitution patterns 58% of the cloned sequences were identified as pseudogenes. Recombination could be traced in 19% of the inferred nuclear pseudogenes, but this figure probably represents a significant underestimation of the factual recombination events. The nonrecombined pseudogenes consisted of multiple haplotypes found to diverge from 1 to 16% from the mitochondrial gene. The number of mitochondrial nuclear copies and their apparent frequent recombination suggest that pseudogenes constitute a serious potential risk in confounding phylogenetic studies and population genetic analysis.

Genomic gigantism: DNA loss is slow in mountain grasshoppers

Several studies have shown DNA loss to be inversely correlated with genome size in animals. These studies include a comparison between Drosophila and the cricket, Laupala, but there has been no assessment of DNA loss in insects with very large genomes. Podisma pedestris, the brown mountain grasshopper, has a genome over 100 times as large as that of Drosophila and 10 times as large as that of Laupala. We used 58 paralogous nuclear pseudogenes of mitochondrial origin to study the characteristics of insertion, deletion, and point substitution in P. pedestris and Italopodisma. In animals, these pseudogenes are "dead on arrival"; they are abundant in many different eukaryotes, and their mitochondrial origin simplifies the identification of point substitutions accumulated in nuclear pseudogene lineages. There appears to be a mononucleotide repeat within the 643-bp pseudogene sequence studied that acts as a strong hot spot for insertions or deletions (indels). Because the data for other insect species did not contain such an unusual region, hot spots were excluded from species comparisons. The rate of DNA loss relative to point substitution appears to be considerably and significantly lower in the grasshoppers studied than in Drosophila or Laupala. This suggests that the inverse correlation between genome size and the rate of DNA loss can be extended to comparisons between insects with large or gigantic genomes (i.e., Laupala and Podisma). The low rate of DNA loss implies that in grasshoppers, the accumulation of point mutations is a more potent force for obscuring ancient pseudogenes than their loss through indel accumulation, whereas the reverse is true for Drosophila. The main factor contributing to the difference in the rates of DNA loss estimated for grasshoppers, crickets, and Drosophila appears to be deletion size. Large deletions are relatively rare in Podisma and Italopodisma.

A novel chimeric mitochondrial RNA localized in the nucleus of mouse sperm

Six identical cDNA clones corresponding to an RNA of 1685 nucleotides that is enriched in mouse sperm compared with testis were isolated from a mouse testis cDNA library. The sequence of these clones corresponds to the 16S mitochondrial RNA plus an inverted repeat of 120 bp covalently joined to the 5' end of the RNA. By RT-PCR, it was demonstrated that this transcript, referred to as chimeric RNA, was present in mouse sperm, testis, liver, kidney, brain, and spleen. The absence of an equivalent sequence in mitochondrial DNA or as a mitochondrial pseudogene in total DNA extracted from sperm, testis, and somatic tissues suggests that the chimeric RNA is a post-transcriptional product, maybe resulting from a trans splicing reaction. The chimeric RNA was found by RT-PCR in total RNA extracted from purified sperm heads. This result was confirmed by in situ hybridization, which showed clear staining of the sperm nucleus with probes corresponding to sequences of the mitochondrial 16S RNA and the inverted repeat.

Multiple nuclear insertions of mitochondrial cytochrome b sequences in callitrichine primates

We report the presence of four nuclear paralogs of a 380-bp segment of cytochrome b in callitrichine primates (marmosets and tamarins). The mitochondrial cytochrome b sequence and each nuclear paralog were obtained from several species, allowing multiple comparisons of rates and patterns of substitution both between mitochondrial and nuclear sequences and among nuclear sequences. The mitochondrial DNA had high overall rates of molecular evolution and a strong bias toward substitutions at third codon positions. Rates of molecular evolution among the nuclear sequences were low and constant, and there were small differences in substitution patterns among the nuclear clades which were probably attributable to the small number of sites involved. A novel method of phylogenetic reconstruction based on the large difference in rates of evolution at different codon positions among mitochondrial and nuclear clades was used to determine whether different nuclear paralogs represent independent transposition events or duplications following a single insertion. This method is generally applicable in cases where differences in pattern of molecular evolution are known, and it showed that at least three of the four nuclear clades represent independent transposition events. The insertion events giving rise to two of the nuclear clades predate the divergence of the callitrichines, whereas those leading to the other two nuclear clades may have occurred in the common ancestor of marmosets.

Mitochondrial DNA-like sequence in the nuclear genome of Saguinus (Callitrichinae, Primates): transfer estimation

Mitochondrial DNA-like sequences have been found in the nuclei of a variety of organisms. These nuclear pseudogenes can be used to estimate relative evolutionary rates of mitochondrial genes, and can be used as outgroups in phylogenetic analyses. In this study, mitochondrial sequences with pseudogene-like characteristics, including deletions and/or insertions and stop codons, were found in tamarins (Saguinus spp., Callitrichinae, Primates). Phylogenetic analysis allowed estimation of the timing of the migration of these sequences to the nuclear genome, and also permitted inferences on the phylogeny of the genus. The choice of an inadequate outgroup (Aotus infulatus) prevented a good phylogenetic resolution of the subfamily Callitrichinae. The relatively ancient divergence of the Cebidae (Callitrichinae, Aotinae and Cebinae) may have favored confounding homoplasies.

Frequent assimilation of mitochondrial DNA by grasshopper nuclear genomes

Multiple copies of mitochondrial-like DNA were found in the brown mountain grasshopper, Podisma pedestris (Orthoptera: Acrididae), paralogous to COI and ND5 regions. The same was discovered using the ND5 regions of nine other grasshopper species from four separate subfamilies (Podisminae, Calliptaminae, Cyrtacanthacridinae, and Gomphocerinae). The extra ND5-like sequences were shown to be nuclear in the desert locust, Schistocerca gregaria (Cyrtacanthacridinae), and probably so in P. pedestris and an Italopodisma sp. (Podisminae). Eighty-seven different ND5-like nuclear mitochondrial pseudogenes (Numts) were sequenced from 12 grasshopper individuals. Different nuclear mitochondrial pseudogenes, if descended from the same mitochondrial immigrant, will have diverged from each other under no selective constraints because of their loss of functionality. Evidence of selective constraints in the differences between any two Numt sequences (e.g., if most differences are at third positions of codons) implies that they have separate mitochondrial origins. Through pairwise comparisons of pseudogene sequences, it was established that there have been at least 12 separate mtDNA integrations into P. pedestris nuclear genomes. This is the highest reported rate of horizontal transfer between organellar and nuclear genomes within a single animal species. The occurrence of numerous mitochondrial pseudogenes in nuclear genomes derived from separate integration events appears to be a common phenomenon among grasshoppers. More than one type of mechanism appears to have been involved in generating the observed grasshopper Numts.

Nuclear integrations of mitochondrial DNA in primates: inference of associated mutational events

To infer the possible mutational events taking place along the interorganellar transfer of genetic material from mitochondria to the nucleus, four integrations of mitochondrial DNA (mtDNA) in the human genome were characterized together with their flanking nuclear sequences. By determining their presence/absence status in different primate species, these integrations were inferred to have occurred on the lineages leading to catarrhines (Old World monkeys and hominoids), to hominoids and to humans, respectively. In case of a polymorphic state, with respect to its presence in a certain species, each preintegration sequence was either cloned in the same species or in a primate taxon that branched off before the transfer of the mtDNA to the nucleus took place. For the four mtDNA integrations presented here, random mobilization of the mtDNA and differing mechanisms for generating free ends in the nuclear target sequences can be inferred. Additionally, no common sequence features at the preintegration sites could be observed for these integrations. Moreover, the comparisons of the sites before and after integration suggest different ways of integration. Thus, mtDNA integrations represent unique molecular recombinations in the evolutionary history and can, according to their presence/absence status in different species, help to determine the branching order in phylogenetic trees.

Nuclear pseudogenes of mitochondrial DNA as a variable part of the human genome

Novel pseudogenes homologous to the mitochondrial (mt) 16S rRNA gene were detected via different approaches. Eight pseudogenes were sequenced. Copy number polymorphism of the mtDNA pseudogenes was observed among randomly chosen individuals, and even among siblings. A mtDNA pseudogene in the Y-chromosome was observed in a YAC clone carrying only repetitive sequence tag site (STS). PCR screening of human yeast artificial chromosome (YAC) libraries showed that there were at least 5.7 x 10(5) bp of the mtDNA pseudogenes in each haploid nuclear genome. Possible involvement of the mtDNA pseudogenes in the variable part of the human nuclear genome is discussed.

Introduction of heteroplasmic mitochondrial DNA (mtDNA) from a patient with NARP into two human rho degrees cell lines is associated either with selection and maintenance of NARP mutant mtDNA or failure to maintain mtDNA

Mitochondria from a patient heteroplasmic at nucleo-tide position 8993 of mitochondrial DNA (mtDNA) were introduced into two human tumour cell lines lacking mtDNA. The donor mitochondria contained between 85 and 95% 8993G:C mtDNA. All detectable mtDNA in the mitochondrially transformed cells contained the pathological 8993G:C mutation 3 months after transformation. These results suggest that 8993G:C mtDNA had a selective advantage over 8993T:A mtDNA in both lung carcinoma and osteo-sarcoma cell backgrounds. In contrast, two other presumed pathological mtDNA variants were lost in favour of 'wild-type' mtDNA molecules in the same lung carcinoma cell background. Taken together, these findings suggest that the transmission bias of mtDNA variants is dependent upon a combination of nuclear background and mtDNA genotype. A second phenomenon observed was a marked decrease in the growth rate of many putative transformed cell lines after 6 weeks of culturing in selective medium, and in these cell lines mtDNA was not readily detectable by Southern blotting. Restriction endonuclease analysis and sequencing of amplified mtDNA demonstrated that the slow growing cells contained little or no mtDNA. It is concluded that these cells represented transient mitochondrial transformants.

A novel mitochondrial DNA-like sequence in the human nuclear genome

We describe here a nuclear mitochondrial DNA-like sequence (numtDNA) that is nearly identical in sequence to a continuous 5842 bp segment of human mitochondrial DNA (mtDNA) that spans nucleotide positions 3914 to 9755. On the basis of evolutionary divergence among modern primates, this numtDNA molecule appears to represent mtDNA from a hominid ancestor that has been translocated to the nuclear genome during the recent evolution of humans. This numtDNA sequence harbors synonymous and nonsynonymous nucleotide substitutions relative to the authentic human mtDNA sequence, including an array of substitutions that was previously found in the cytochrome c oxidase subunit 1 and 2 genes. These substitutions were previously reported to occur in human mtDNA, but subsequently contended to be present in a nuclear pseudogene sequence. We now demonstrate their exclusive association with this 5842-bp numtDNA, which we have characterized in its entirety. This numtDNA does not appear to be expressed as a mtDNA-encoded mRNA. It is present in nuclear DNA from human blood donors, in human SH-SY5Y and A431 cell lines, and in rho(0) SH-SY5Y and rho(0) A431 cell lines that were depleted of mtDNA. The existence of human numtDNA sequences with great similarities to human mtDNA renders the amplification of pure mtDNA from cellular DNA very difficult, thereby creating the potential for confounding studies of mitochondrial diseases and population genetics.

Applications of constant denaturant capillary electrophoresis/high-fidelity polymerase chain reaction to human genetic analysis

Constant denaturant capillary electrophoresis (CDCE) permits high-resolution separation of single-base variations occurring in an approximately 100 bp isomelting DNA sequence based on their differential melting temperatures. By coupling CDCE for highly efficient enrichment of mutants with high-fidelity polymerase chain reaction (hifi PCR), we have developed an analytical approach to detecting point mutations at frequencies equal to or greater than 10(-6) in human genomic DNA. In this article, we present several applications of this approach in human genetic studies. We have measured the point mutational spectra of a 100 bp mitochondrial DNA sequence in human tissues and cultured cells. The observations have led to the conclusion that the primary causes of mutation in human mitochondrial DNA are spontaneous in origin. In the course of studying the mitochondrial somatic mutations, we have also identified several nuclear pseudogenes homologous to the analyzed mitochondrial DNA fragment. Recently, through developments of the means to isolate the desired target sequences from bulk genomic DNA and to increase the loading capacity of CDCE, we have extended the CDCE/hifi PCR approach to study a chemically induced mutational spectrum in a single-copy nuclear sequence. Future applications of the CDCE/hifi PCR approach to human genetic analysis include studies of somatic mitochondrial mutations with respect to aging, measurement of mutational spectra of nuclear genes in healthy human tissues and population screening for disease-associated single nucleotide polymorphisms (SNPs) in large pooled samples.

Expression of atrC - encoding a novel member of the ATP binding cassette transporter family in Aspergillus nidulans - is sensitive to cycloheximide

A new member of the ABC superfamily of transmembrane proteins in Aspergillus nidulans has been cloned and characterized. The topology of conserved motifs subgroups AtrC in the P-glycoprotein cluster of ABC permeases, the members of this subfamily, are known to participate in multidrug resistance (MDR) in diverse organisms. Alignment results display significant amino acid similarity to AfuMDR1 and AflMDR1 from Aspergillus fumigatus and flavus, respectively. Northern analysis reveals that atrC mRNA levels are 10-fold increased in response to cycloheximide. Evidence for the existence of eight additional hitherto unpublished ABC transporter proteins in A. nidulans is provided.

Long-extension PCR to detect deleted mitochondrial DNA molecules is compromized by technical artefacts

Long-extension PCR (LX-PCR), followed by Southern hybridization to probes for two different regions of the mitochondrial genome, was used to evaluate the presence of deleted mtDNA molecules in heart muscle samples from alcoholic cardiomyopathy patients compared with age-matched controls. Two different primer pairs capable of amplifying the entire genome, as well as a variety of other primer pairs predicted to amplify the genome in large, overlapping fragments, were tested. Products indicating the presence of a variety of subgenomic, deleted molecules were detected in variable amounts from patient and control myocardial samples alike. Most of these hybridized with a probe for the 16S/ND1 region, but not with a probe for the ND4/ND5 region that is commonly deleted. Dilution of a given template DNA in which deleted products were prominent resulted in the disappearance of the subgenomic bands in favour of the full-length, undeleted product. Therefore, the appearance and amount of such products is subject to template concentration or quality. The results indicate that the application of LX-PCR to the detection and quantitation of deleted mtDNAs is inherently unreliable, and findings using this technique should be treated with caution unless supported by an independent method.

Selected nuclear LINE elements with mitochondrial-DNA-like inserts are more plentiful and mobile in tumor than in normal tissue of mouse and rat

The nuclear DNA of normal and tumor mouse and rat tissue was examined for mitochondrial-DNA-like inserts by means of the Southern blot technique. The two probes were 32P-labeled cloned mitochondrial DNA. KpnI, which doesn't cut either mitochondrial DNA, was one of the restriction enzymes, while the enzymes that fragment mitochondrial DNA were for mouse and rat PstI and BamHI, respectively. When KpnI alone was used in the procedure a nuclear LINE family whose elements had mitochondrial-DNA-like insertions was selected. Such elements were much more abundant in tumor than in normal tissue. The results with PstI alone and BamHI alone and each combined with KpnI indicated that there were mobile LINE elements with mitochondrial-DNA-like inserts in the nuclear genome of tumor. The mouse tissues were normal liver and a transplantable lymphoid leukemic ascites cell line L1210 that had been carried for 40 years. The rat tissues were normal liver and a hepatoma freshly induced by diethylnitrosoamine in order to minimize the role of 40 years of transplantation. Our unitary hypothesis for carcinogenesis of 1971, which suggested these experiments, has been augmented to include mobile nuclear elements with inserts of mitochondrial-DNA-like sequences. Such elements have been related to diseases of genetic predisposition such as breast cancer and Huntington's disease.

Ancient mtDNA sequences in the human nuclear genome: a potential source of errors in identifying pathogenic mutations

Nuclear-localized mtDNA pseudogenes might explain a recent report describing a heteroplasmic mtDNA molecule containing five linked missense mutations dispersed over the contiguous mtDNA CO1 and CO2 genes in Alzheimer's disease (AD) patients. To test this hypothesis, we have used the PCR primers utilized in the original report to amplify CO1 and CO2 sequences from two independent rho degrees (mtDNA-less) cell lines. CO1 and CO2 sequences amplified from both of the rho degrees cells, demonstrating that these sequences are also present in the human nuclear DNA. The nuclear pseudogene CO1 and CO2 sequences were then tested for each of the five "AD" missense mutations by restriction endonuclease site variant assays. All five mutations were found in the nuclear CO1 and CO2 PCR products from rho degrees cells, but none were found in the PCR products obtained from cells with normal mtDNA. Moreover, when the overlapping nuclear CO1 and CO2 PCR products were cloned and sequenced, all five missense mutations were found, as well as a linked synonymous mutation. Unlike the findings in the original report, an additional 32 base substitutions were found, including two in adjacent tRNAs and a two base pair deletion in the CO2 gene. Phylogenetic analysis of the nuclear CO1 and CO2 sequences revealed that they diverged from modern human mtDNAs early in hominid evolution about 770,000 years before present. These data would be consistent with the interpretation that the missense mutations proposed to cause AD may be the product of ancient mtDNA variants preserved as nuclear pseudogenes.

Apparent mtDNA heteroplasmy in Alzheimer's disease patients and in normals due to PCR amplification of nucleus-embedded mtDNA pseudogenes

In an unprecedented finding, Davis et al. [Davis, R. E., Miller, S., Herrnstadt, C., Ghosh, S. S., Fahy, E., Shinobu, L. A., Galasko, D., Thal, L. J., Beal, M. F., Howell, N. & Parker, W. D., Jr. (1997) Proc. Natl. Acad. Sci. USA 94, 4526-4531] used an unusual DNA isolation method to show that healthy adults harbor a specific population of mutated mitochondrial cytochrome c oxidase (COX) genes that coexist with normal mtDNAs. They reported that this heteroplasmic population was present at a level of 10-15% in the blood of normal individuals and at a significantly higher level (20-30%) in patients with sporadic Alzheimer's disease. We provide compelling evidence that the DNA isolation method employed resulted in the coamplification of authentic mtDNA-encoded COX genes together with highly similar COX-like sequences embedded in nuclear DNA ("mtDNA pseudogenes"). We conclude that the observed heteroplasmy is an artifact.

Molecular evolution of a portion of the mitochondrial 16S ribosomal gene region in scleractinian corals

Relationships among families and suborders of scleractinian corals are poorly understood because of difficulties 1) in making inferences about the evolution of the morphological characters used in coral taxonomy and 2) in interpreting their 240-million-year fossil record. Here we describe patterns of molecular evolution in a segment of the mitochondrial (mt) 16S ribosomal gene from taxa of 14 families of corals and the use of this gene segment in a phylogenetic analysis of relationships within the order. We show that sequences obtained from scleractinians are homologous to other metazoan 16S ribosomal sequences and fall into-two distinct clades defined by size of the amplified gene product. Comparisons of sequences from the two clades demonstrate that both sets of sequences are evolving under similar evolutionary constraints: they do not differ in nucleotide composition, numbers of transition and transversion substitutions, spatial patterns of substitutions, or in rates of divergence. The characteristics and patterns observed in these sequences as well as the secondary structures, are similar to those observed in mt 16S ribosomal DNA sequences from other taxa. Phylogenetic analysis of these sequences shows that they are useful for evaluating relationships within the order. The hypothesis generated from this analysis differs from traditional hypotheses for evolutionary relationships among the Scleractinia and suggests that a reevaluation of evolutionary affinities in the order is needed.

Multiple independent transpositions of mitochondrial DNA control region sequences to the nucleus

Transpositions of mtDNA sequences to the nuclear genome have been documented in a wide variety of individual taxa, but little is known about their taxonomic frequency or patterns of variation. We provide evidence of nuclear sequences homologous to the mtDNA control region in seven species of diving ducks (tribe Aythyini). Phylogenetic analysis places each nuclear sequence as a close relative of the mtDNA haplotypes of the specie(s) in which it occurs, indicating that they derive from six independent transposition events, all occurring within the last approximately 1.5 million years. Relative-rate tests and comparison of intraspecific variation in nuclear and mtDNA sequences confirm the expectation of a greatly reduced rate of evolution in the nuclear copies. By representing mtDNA haplotypes from ancestral populations, nuclear insertions may be valuable in some phylogenetic analyses, but they also confound the accurate determination of mtDNA sequences. In particular, our data suggest that the presumably nonfunctional but more slowly evolving nuclear sequences often will not be identifiable by changes incompatible with function and may be preferentially amplified by PCR primers based on mtDNA sequences from related taxa.

Human evolution and the mitochondrial genome

The use of mitochondrial DNA (mtDNA) continues to dominate studies of human genetic variation and evolution. Recent work has re-affirmed the strict maternal inheritance of mtDNA, yielded new insights into the extent and nature of intra-individual variation, supported a recent African origin of human mtDNA, and amply demonstrated the utility of mtDNA in tracing population history and in analyses of ancient remains.

Evidence for association of mitochondrial DNA sequence amplification and nuclear localization in human low-grade gliomas

Gliomas are tumors which have been found to exhibit consistent genetic changes. Recent studies have shown mitochondrial DNA is also altered in these tumors, and include large deletions and gene amplification. Other studies of the mitochondrial genome in cancer have revealed a variety of different alterations, including the localization and insertion of mitochondrial DNA into the nucleus and nuclear genome in HeLa cells and diethylnitrosurea-induced hepatoma cells. Whether these changes are ontogenically early in the multistep pathway to the development of malignancy, or if this phenomenon occurs in human glial tumors is unknown. I sought to study these questions in a panel of unselected primary glial tumors of pathologically low grade. Fifteen tumors were assessed with a mitochondrial cDNA probe with homology to positions 1679-1948, and 2017-2057. All low-grade tumors revealed increases in copy number when compared to a normal brain control. Nuclear suspensions of these tumors were evaluated by fluorescent in situ hybridization (FISH), using the entire mitochondrial genome as a probe after labeling with rhodamine. All tumors showed evidence of mitochondrial sequence localization within the nuclei. A corresponding glioblastoma and two normal brain specimens were also evaluated which did not have amplification of the mitochondrial genome; FISH with the mitochondrial probe revealed minimal hybridization signal within the nuclei of these samples. Mitochondrial DNA nuclear localization can be found in primary low-grade brain neoplasms, and is correlated to increases in mitochondrial DNA.

Highly conserved nuclear copies of the mitochondrial control region in the desert locust Schistocerca gregaria: some implications for population studies

Animal mitochondrial DNA has proved a valuable marker in intraspecific systematic studies. However, if nucleotide sequence heterogeneity exists at the individual level, its usefulness will be much reduced. The study demonstrates that the presence of highly conserved non-coding mitochondrial sequences in the nuclear genome of Schistocerca gregaria greatly impairs the use of mtDNA in population genetic studies. Caution is called for in other organisms; and it seems necessary to check for conserved nuclear copies of mitochondrial sequences before launching into a large scale analysis of populations using mtDNA as a genetic marker. Experimental procedures are suggested for this purpose.

Mitochondrial DNA: molecular fossils in the nucleus

The nuclear genomes of many animals contain non-functional copies of mitochondrial genes that provide new opportunities for evolutionary analysis.

Rotenone, a mitochondrial NADH dehydrogenase inhibitor, induces cell surface expression of CD13 and CD38 and apoptosis in HL-60 cells

We previously demonstrated that the mitochondrial NADH dehydrogenase subunit 2 (ND2) gene was overexpressed in human acute myelogenous leukemia (AML) cells. Since this finding suggested that ND2 gene expression was related to myeloid differentiation, we here investigated the effects of rotenone, a specific NADH dehydrogenase inhibitor, on HL-60 cell growth, differentiation and death. Fifty nM rotenone inhibited the growth of HL-60 cells and caused an increase in the cell population in the G(2) +M phase. In the quantitative comparison of myeloid antigen, the expression of CD13 and CD38 were relatively increased in the rotenone-treated cells. These findings suggest that the inhibition of NADH dehydrogenase changes the cell cycle and induces some specific surface antigens of HL-60 cells. On the other hand, the expression of ND2 gene remained unchanged after the rotenone treatment, suggesting the rotenone-mediated mitochondrial inhibition did not affect the mitochondrial gene expression. Five mu M rotenone strongly inhibited the cellular proliferation. Electron microscopy and an electrophoretic analysis of DNA showed that the majority of the HL-60 cells were induced into typical apoptosis within 24-48 hours. On the basis of this and other studies, we believe that mitochondrial function is directly involved in both cellular differentiation and apoptotic cell death.

Escape and migration of nucleic acids between chloroplasts, mitochondria, and the nucleus

The escape and migration of genetic information between mitochondria, chloroplasts, and nuclei have been an integral part of evolution and has a continuing impact on the biology of cells. The evolutionary transfer of functional genes and fragments of genes from chloroplasts to mitochondria, from chloroplasts to nuclei, and from mitochondria to nuclei has been documented for numerous organisms. Most documented instances of genetic material transfer have involved the transfer of information from mitochondria and chloroplasts to the nucleus. The pathways for the escape of DNA from organelles may include transient breaches in organellar membranes during fusion and/or budding processes, terminal degradation of organelles by autophagy coupled with the subsequent release of nucleic acids to the cytoplasm, illicit use of nucleic acid or protein import machinery, or fusion between heterotypic membranes. Some or all of these pathways may lead to the escape of DNA or RNA from organellar compartments with subsequent uptake of nucleic acids from the cytoplasm into the nucleus. Investigations into the escape of DNA from mitochondria in yeast have shown the rate of escape for gene-sized fragments of DNA from mitochondria and its subsequent migration to the nucleus to be roughly equivalent to the rate of spontaneous mutation of nuclear genes. Smaller fragments of mitochondrial DNA may appear in the nucleus even more frequently. Mutations of nuclear genes that define gene products important in controlling the rate of DNA escape from mitochondria in yeast also have been described. The escape of genetic material from mitochondria and chloroplasts has clearly had an impact on nuclear genetic organization throughout evolution and may also affect cellular metabolic processes.