Mitochondrial DNA: molecular fossils in the nucleus

Comparative Genome Microsynteny Illuminates the Fast Evolution of Nuclear Mitochondrial Segments (NUMTs) in Mammals

The escape of DNA from mitochondria into the nuclear genome (nuclear mitochondrial DNA, NUMT) is an ongoing process. Although pervasively observed in eukaryotic genomes, their evolutionary trajectories in a mammal-wide context are poorly understood. The main challenge lies in the orthology assignment of NUMTs across species due to their fast evolution and chromosomal rearrangements over the past 200 million years. To address this issue, we systematically investigated the characteristics of NUMT insertions in 45 mammalian genomes and established a novel, synteny-based method to accurately predict orthologous NUMTs and ascertain their evolution across mammals. With a series of comparative analyses across taxa, we revealed that NUMTs may originate from nonrandom regions in mtDNA, are likely found in transposon-rich and intergenic regions, and unlikely code for functional proteins. Using our synteny-based approach, we leveraged 630 pairwise comparisons of genome-wide microsynteny and predicted the NUMT orthology relationships across 36 mammals. With the phylogenetic patterns of NUMT presence-and-absence across taxa, we constructed the ancestral state of NUMTs given the mammal tree using a coalescent method. We found support on the ancestral node of Fereuungulata within Laurasiatheria, whose subordinal relationships are still controversial. This study broadens our knowledge on NUMT insertion and evolution in mammalian genomes and highlights the merit of NUMTs as alternative genetic markers in phylogenetic inference.

Interrogating 1000 insect genomes for NUMTs: A risk assessment for estimates of species richness

The nuclear genomes of most animal species include NUMTs, segments of the mitogenome incorporated into their chromosomes. Although NUMT counts are known to vary greatly among species, there has been no comprehensive study of their frequency/attributes in the most diverse group of terrestrial organisms, insects. This study examines NUMTs derived from a 658 bp 5' segment of the cytochrome c oxidase I (COI) gene, the barcode region for the animal kingdom. This assessment is important because unrecognized NUMTs can elevate estimates of species richness obtained through DNA barcoding and derived approaches (eDNA, metabarcoding). This investigation detected nearly 10,000 COI NUMTs ≥ 100 bp in the genomes of 1,002 insect species (range = 0-443). Variation in nuclear genome size explained 56% of the mitogenome-wide variation in NUMT counts. Although insect orders with the largest genome sizes possessed the highest NUMT counts, there was considerable variation among their component lineages. Two thirds of COI NUMTs possessed an IPSC (indel and/or premature stop codon) allowing their recognition and exclusion from downstream analyses. The remainder can elevate species richness as they showed 10.1% mean divergence from their mitochondrial homologue. The extent of exposure to "ghost species" is strongly impacted by the target amplicon's length. NUMTs can raise apparent species richness by up to 22% when a 658 bp COI amplicon is examined versus a doubling of apparent richness when 150 bp amplicons are targeted. Given these impacts, metabarcoding and eDNA studies should target the longest possible amplicons while also avoiding use of 12S/16S rDNA as they triple NUMT exposure because IPSC screens cannot be employed.

A Burst of Numt Insertion in the Dasyuridae Family During Marsupial Evolution

Nuclear pseudogenes of mitochondrial origin (numts) are common in all eukaryotes. Our previous scan of numts in sequenced nuclear genomes suggested that the highest numt content currently known in animals is that in the gray short-tailed opossum. The present work sought to determine numt content in marsupials and to compare it to those in placental and monothematic mammals as well as in non-mammalian vertebrates. To achieve this, 70 vertebrate species with available nuclear and mitochondrial genomes were scanned for numt content. An extreme numt content was found in the Dasyuridae, with 3,450 in Sarcophilus harrisii (1,955 kb) and 2,813 in Antechinus flavipes (847 kb). The evolutionarily closest species analyzed, the extinct Thylacinus cynocephalus belonging to the Thylacindae family, had only 435 numts (238 kb). These two Dasyuridae genomes featured the highest numt content identified in animals to date. A phylogenetic analysis of numts longer than 300 bp, using a Diprotodonita mitochondrial tree, indicated a burst of numt insertion that began before the divergence of the Dasyurini and Phascogalini, reaching a peak in the early evolution of the two tribes. No comparable increase was found in the early divergent species T. cynocephalus. Divergence of the Dasyuridae tribes has been previously dated to shortly after the Miocene climate transition, characterized by a rapid temperature decline. Interestingly, deviation from optimal growth temperature is one of the environmental factors reported to increase numt insertions in a laboratory setting.

Using metatranscriptomics to estimate the diversity and composition of zooplankton communities

DNA metabarcoding is a rapid, high‐resolution tool used for biomonitoring complex zooplankton communities. However, diversity estimates derived with this approach can be biased by the co‐detection of sequences from environmental DNA (eDNA), nuclear‐encoded mitochondrial (NUMT) pseudogene contamination, and taxon‐specific PCR primer affinity differences. To avoid these methodological uncertainties, we tested the use of metatranscriptomics as an alternative approach for characterizing zooplankton communities. Specifically, we compared metatranscriptomics with PCR‐based methods using genomic (gDNA) and complementary DNA (cDNA) amplicons, and morphology‐based data for estimating species diversity and composition for both mock communities and field‐collected samples. Mock community analyses showed that the use of gDNA mitochondrial cytochrome c oxidase I (mtCO1) amplicons inflates species richness due to the co‐detection of extra‐organismal eDNA. Significantly more amplicon sequence variants, nucleotide diversity, and indels were observed with gDNA amplicons than with cDNA, indicating the presence of putative NUMT pseudogenes. Moreover, PCR‐based methods failed to detect the most abundant species in mock communities due to priming site mismatch. Overall, metatranscriptomics provided estimates of species richness and composition that closely resembled those derived from morphological data. The use of metatranscriptomics was further tested using field‐collected samples, with the results showing consistent species diversity estimates among biological and technical replicates. Additionally, temporal zooplankton species composition changes could be monitored using different mitochondrial markers. These findings demonstrate the advantages of metatranscriptomics as an effective tool for monitoring diversity in zooplankton research.

Profile hidden Markov model sequence analysis can help remove putative pseudogenes from DNA barcoding and metabarcoding datasets

BACKGROUND

Pseudogenes are non-functional copies of protein coding genes that typically follow a different molecular evolutionary path as compared to functional genes. The inclusion of pseudogene sequences in DNA barcoding and metabarcoding analysis can lead to misleading results. None of the most widely used bioinformatic pipelines used to process marker gene (metabarcode) high throughput sequencing data specifically accounts for the presence of pseudogenes in protein-coding marker genes. The purpose of this study is to develop a method to screen for nuclear mitochondrial DNA segments (nuMTs) in large COI datasets. We do this by: (1) describing gene and nuMT characteristics from an artificial COI barcode dataset, (2) show the impact of two different pseudogene removal methods on perturbed community datasets with simulated nuMTs, and (3) incorporate a pseudogene filtering step in a bioinformatic pipeline that can be used to process Illumina paired-end COI metabarcode sequences. Open reading frame length and sequence bit scores from hidden Markov model (HMM) profile analysis were used to detect pseudogenes.

RESULTS

Our simulations showed that it was more difficult to identify nuMTs from shorter amplicon sequences such as those typically used in metabarcoding compared with full length DNA barcodes that are used in the construction of barcode libraries. It was also more difficult to identify nuMTs in datasets where there is a high percentage of nuMTs. Existing bioinformatic pipelines used to process metabarcode sequences already remove some nuMTs, especially in the rare sequence removal step, but the addition of a pseudogene filtering step can remove up to 5% of sequences even when other filtering steps are in place.

CONCLUSIONS

Open reading frame length filtering alone or combined with hidden Markov model profile analysis can be used to effectively screen out apparent pseudogenes from large datasets. There is more to learn from COI nuMTs such as their frequency in DNA barcoding and metabarcoding studies, their taxonomic distribution, and evolution. Thus, we encourage the submission of verified COI nuMTs to public databases to facilitate future studies.

Concerted and birth-and-death evolution of 26S ribosomal DNA in Camellia L

BACKGROUND AND AIMS

The ribosomal DNA (rDNA) gene family, encoding ribosomal RNA (rRNA), has long been regarded as an archetypal example illustrating the model of concerted evolution. However, controversy is arising, as rDNA in many eukaryotic species has been proved to be polymorphic. Here, a metagenomic strategy was applied to detect the intragenomic polymorphism as well as the evolutionary patterns of 26S rDNA across the genus Camellia.

METHODS

Degenerate primer pairs were designed to amplify the 26S rDNA fragments from different Camellia species. The amplicons were then paired-end sequenced on the Illumina MiSeq platform.

KEY RESULTS

An extremely high level of rDNA polymorphism existed universally in Camellia. However, functional rDNA was still the major component of the family, and was relatively conserved among different Camellia species. Sequence variations mainly came from rRNA pseudogenes and favoured regions that are rich in GC. Specifically, some rRNA pseudogenes have existed in the genome for a long time, and have even experienced several expansion events, which has greatly enriched the abundance of rDNA polymorphism.

CONCLUSIONS

Camellia represents a group in which rDNA is subjected to a mixture of concerted and birth-and-death evolution. Some rRNA pseudogenes may still have potential functions. Conversely, when released from selection constraint, they can evolve in the direction of decreasing GC content and structural stability through a methylation-induced process, and finally be eliminated from the genome.

Mitochondrial heteroplasmy and pseudogenes in the freshwater prawn, Macrobrachium amazonicum (Heller, 1862): DNA barcoding and phylogeographic implications

The mitochondrial cytochrome oxidase c subunit 1 (COI) gene has been widely used in phylogenetic studies of crustaceans and analyses in population genetics. As COI studies have become more popular, there has been an increase in the number of reports of the presence of nuclear insertions of mitochondrial DNA (Numts) and mitochondrial heteroplasmy. Here, we provide evidence of both types of event in the COI sequences of Macrobrachium amazonicum, an economically important freshwater prawn, which is widespread in South America. Heteroplasmy and Numts were confirmed by different methods of DNA extraction (genomic, mitochondrial, and nuclear-enriched DNA), cloning, and sequencing, and were observed in 11 of the 14 populations sampled, primarily in the Amazon region. We discuss how the occurrence of these events affects the interpretation of the genetic relationships among the M. amazonicum populations, and we recommend caution when using COI for genetic inferences in prawns of the genus Macrobrachium, and in particular that any analysis should include nuclear markers.

Tracking the Distribution and Burst of Nuclear Mitochondrial DNA Sequences (NUMTs) in Fig Wasp Genomes

Simple Summary

Nuclear mitochondrial DNA sequences (NUMTs), which result from the insertion of exogenous mtDNA into the nuclear genome, are widely distributed in eukaryotes. However, how NUMTs are inserted into the nuclear genome and their post-insertion fates remain a mystery. Previous studies have suggested that Hymenoptera may be a group rich in NUMTs, which will be helpful to study the biological issues of NUMTs. We here select 11 species of fig wasps (Chalcidoidea, Hymenoptera) to analyze the distribution and evolution of NUMTs at the genomic level. The results show that the distributions of NUMTs are species- or lineage-specific. Furthermore, genomic environmental factors such as genome size, the damage-prone regions, and the mode of TE dynamics can determine the insertion and post-insertion fate of NUMTs. Especially because of TEs, the fragmentation and duplication of NUMTs, and thus their burst, are common. This is a relatively comprehensive investigation of the specific distribution of NUMTs and its influencing factors. Our study will help people to understand the evolution of exogenous fragments in the nuclear genome.

Abstract

Mitochondrial DNA sequences can be transferred into the nuclear genome, giving rise to nuclear mitochondrial DNA sequences (NUMTs). NUMTs have been described in numerous eukaryotes. However, the studies on the distribution of NUMTs and its influencing factors are still inadequate and even controversial. Previous studies have suggested that Hymenoptera may be a group rich in NUMTs, in which we selected 11 species of fig wasps (Chalcidoidea, Hymenoptera) to analyze the distribution and evolution of NUMTs at the genomic level. The results showed that the contents of NUMTs varied greatly in these species, and bursts of NUMTs existed in some species or lineages. Further detailed analyses showed that the large number of NUMTs might be related to the large genomes; NUMTs tended to be inserted into unstable regions of the genomes; and the inserted NUMTs might also be affected by transposable elements (TEs) in the neighbors, leading to fragmentations and duplications, followed by bursts of NUMTs. In summary, our results suggest that a variety of genomic environmental factors can determine the insertion and post-insertion fate of NUMTs, resulting in their species- or lineage-specific distribution patterns, and that studying the evolution of NUMTs can provide good evidence and theoretical basis for exploring the dynamics of exogenous DNA entering into the nuclear genome.

Novel Mitochondrial DNA Lineage Found among Ochlerotatus communis (De Geer, 1776) of the Nordic-Baltic Region

The Ochlerotatus (Oc.) communis complex consist of three Northern American species as well as a common Holarctic mosquito (Diptera: Culicidae) Oc. communis (De Geer, 1776). These sister species exhibit important ecological differences and are capable of transmitting various pathogens, but cannot always be differentiated by morphological traits. To investigate the Oc. communis complex in Europe, we compared three molecular markers (COI, ND5 and ITS2) from 54 Estonian mosquitoes as well as two COI marker sequences from Sweden. These sequences were subjected to phylogenetic analysis and screened for Wolbachia Hertig and Wolbach symbionts. Within and between groups, distances were calculated for each marker to better understand the relationships among individuals. Results demonstrate that a group of samples, extracted from adult female mosquitoes matching the morphology of Oc. communis, show a marked difference from the main species when comparing the mitochondrial markers COI and ND5. However, there is no variance between the same specimens when considering the nuclear ITS2. We conclude that Oc. communis encompasses two distinct mitochondrial DNA lineages in the Nordic-Baltic region. Further research is needed to investigate the origin and extent of these genetic differences.

Mitonuclear genomics and aging

Our cells operate based on two distinct genomes that are enclosed in the nucleus and mitochondria. The mitochondrial genome presumably originates from endosymbiotic bacteria. With time, a large portion of the original genes in the bacterial genome is considered to have been lost or transferred to the nuclear genome, leaving a reduced 16.5 Kb circular mitochondrial DNA (mtDNA). Traditionally only 37 genes, including 13 proteins, were thought to be encoded within mtDNA, its genetic repertoire is expanding with the identification of mitochondrial-derived peptides (MDPs). The biology of aging has been largely unveiled to be regulated by genes that are encoded in the nuclear genome, whereas the mitochondrial genome remained more cryptic. However, recent studies position mitochondria and mtDNA as an important counterpart to the nuclear genome, whereby the two organelles constantly regulate each other. Thus, the genomic network that regulates lifespan and/or healthspan is likely constituted by two unique, yet co-evolved, genomes. Here, we will discuss aspects of mitochondrial biology, especially mitochondrial communication that may add substantial momentum to aging research by accounting for both mitonuclear genomes to more comprehensively and inclusively map the genetic and molecular networks that govern aging and age-related diseases.

Survey of mitochondrial sequences integrated into the bovine nuclear genome

Nuclear copies of the mitochondrial DNA (NUMTs) have already been described in several species. In this context, we identified and analysed 166 bovine NUMT regions with a total length of 430 kbp, representing about 0.02% of the cattle nuclear genome. Copies of all mitochondrial regions were detected in the nuclear genome, with distinct degrees of sequence similarity to the mitogenome. Some NUMT regions include large mitogenome segments and show high similarity to the organelle genome sequence. NUMT regions are frequently modified by insertion of repetitive sequences and by sequence rearrangements. We confirmed the existence of 29 NUMT regions by PCR amplification using DNA from the cow (Dominette) which was used to generate the bovine genome reference sequence, ruling out the possibility that these NUMTs could be artifacts of the genome assembly. As there are NUMT regions with high similarity to the mitogenome, special care is needed when designing primers for mitochondrial DNA amplification. Our results can therefore be used to avoid co-amplification of bovine nuclear sequences similar to mitochondrial DNA.

Nuclear Integrants of Organellar DNA Contribute to Genome Structure and Evolution in Plants

The transfer of genetic material from the mitochondria and plastid to the nucleus gives rise to nuclear integrants of mitochondrial DNA (NUMTs) and nuclear integrants of plastid DNA (NUPTs). This frequently occurring DNA transfer is ongoing and has important evolutionary implications. In this review, based on previous studies and the analysis of NUMT/NUPT insertions of more than 200 sequenced plant genomes, we analyzed and summarized the general features of NUMTs/NUPTs and highlighted the genetic consequence of organellar DNA insertions. The statistics of organellar DNA integrants among various plant genomes revealed that organellar DNA-derived sequence content is positively correlated with the nuclear genome size. After integration, the nuclear organellar DNA could undergo different fates, including elimination, mutation, rearrangement, fragmentation, and proliferation. The integrated organellar DNAs play important roles in increasing genetic diversity, promoting gene and genome evolution, and are involved in sex chromosome evolution in dioecious plants. The integrating mechanisms, involving non-homologous end joining at double-strand breaks were also discussed.

Archaic mitochondrial DNA inserts in modern day nuclear genomes

Background

Traces of interbreeding of Neanderthals and Denisovans with modern humans in the form of archaic DNA have been detected in the genomes of present-day human populations outside sub-Saharan Africa. Up to now, only nuclear archaic DNA has been detected in modern humans; we therefore attempted to identify archaic mitochondrial DNA (mtDNA) residing in modern human nuclear genomes as nuclear inserts of mitochondrial DNA (NUMTs).

Results

We analysed 221 high-coverage genomes from Oceania and Indonesia using an approach which identifies reads that map both to the nuclear and mitochondrial DNA. We then classified reads according to the source of the mtDNA, and found one NUMT of Denisovan mtDNA origin, present in 15 analysed genomes; analysis of the flanking region suggests that this insertion is more likely to have happened in a Denisovan individual and introgressed into modern humans with the Denisovan nuclear DNA, rather than in a descendant of a Denisovan female and a modern human male.

Conclusions

Here we present our pipeline for detecting introgressed NUMTs in next generation sequencing data that can be used on genomes sequenced in the future. Further discovery of such archaic NUMTs in modern humans can be used to detect interbreeding between archaic and modern humans and can reveal new insights into the nature of such interbreeding events.

Conserved numts mask a highly divergent mitochondrial-COI gene in a species complex of Australian stingless bees Tetragonula (Hymenoptera: Apidae)

Tetragonula carbonaria, Tetragonula davenporti, Tetragonula hockingsi and Tetragonula mellipes comprise a species complex of Australian stingless bee species known as the 'Carbonaria' group. The species are difficult to distinguish morphologically and the major species-defining characters relate to comb architecture and nest entrance ornamentation. The taxonomy of the group is further complicated by likely nuclear mitochondrial pseudogenes (numts) and inter-specific hybrids. Here we demonstrate the existence of COI numts and isolate and characterize the 'true' mt-COI gene in T. carbonaria and T. hockingsi. Numts were isolated from enriched-nuclear DNA extraction followed by PCR amplification and Sanger sequencing, and were recognized by the presence of deletions and/or premature stop codons in the translated sequences. The mt-COI sequences were obtained from NGS sequencing using purified mtDNA. In T. carbonaria, two numts (numt1 and numt2) were identified and a third (numt3) was identified in T. hockingsi. Numt2 and numt3 are similar (1.2% sequence divergence), indicating a recent common origin. The genetic distance between the mt-COI of the two Tetragonula species was higher than might be expected for closely related species, 16.5%, corroborating previous studies in which T. carbonaria and T. hockingsi were regarded as separate species. The three numts are more similar to the COI of other stingless bee species, including Australian Austroplebia australis and South American Melipona bicolor (81.7-83.9%) than to the mt-COI of their own species (70-71.4%). This is because the mt-COI of T. carbonaria and T. hockingsi differ greatly from other Meliponinae. Our findings explain some formerly puzzling aspects of Carbonaria biogeography, and misinterpreted amplifications.

Alignment-free approaches for predicting novel Nuclear Mitochondrial Segments (NUMTs) in the human genome

The nuclear human genome harbors sequences of mitochondrial origin, indicating an ancestral transfer of DNA from the mitogenome. Several Nuclear Mitochondrial Segments (NUMTs) have been detected by alignment-based sequence similarity search, as implemented in the Basic Local Alignment Search Tool (BLAST). Identifying NUMTs is important for the comprehensive annotation and understanding of the human genome. Here we explore the possibility of detecting NUMTs in the human genome by alignment-free sequence similarity search, such as k-mers (k-tuples, k-grams, oligos of length k) distributions. We find that when k=6 or larger, the k-mer approach and BLAST search produce almost identical results, e.g., detect the same set of NUMTs longer than 3 kb. However, when k=5 or k=4, certain signals are only detected by the alignment-free approach, and these may indicate yet unrecognized, and potentially more ancestral NUMTs. We introduce a "Manhattan plot" style representation of NUMT predictions across the genome, which are calculated based on the reciprocal of the Jensen-Shannon divergence between the nuclear and mitochondrial k-mer frequencies. The further inspection of the k-mer-based NUMT predictions however shows that most of them contain long-terminal-repeat (LTR) annotations, whereas BLAST-based NUMT predictions do not. Thus, similarity of the mitogenome to LTR sequences is recognized, which we validate by finding the mitochondrial k-mer distribution closer to those for transposable sequences and specifically, close to some types of LTR.

Characterization, phylogenetic distribution and evolutionary trajectories of diverse hydrocarbon degrading microorganisms isolated from refinery sludge

Phylogenic association between bacteria living under harsh conditions can provide important information on adaptive mechanism, survival strategy and their potential application. Indigenous microorganisms isolated from toxic refinery oily sludge with ability to degrade a diverse range of hydrocarbons were identified and characterized. The strains including Pseudomonas aeruginosa RS1, Microbacterium sp. RS2, Bacillus sp. RS3, Acinetobacter baumannii RS4 and Stenotrophomonas sp. RS5 could utilize n-alkanes, cycloalkanes, polynuclear aromatic hydrocarbons (PAHs) with 2-4 rings and also substituted PAHs as sole substrate. The phylogenetic position of Bacillus sp. RS3 and Pseudomonas sp. RS1 was tested by applying the maximum likelihood (ML) method to the aligned 16S rRNA nucleotide sequences of PAH and aliphatic hydrocarbon degrading strains belonging to the corresponding genus. The base substitution matrix created with each set of organisms capable of degrading aromatic and aliphatic hydrocarbons showed significant transitional event with high values of transition: transversion ratio (R) under all conditions. The guanine-cytosine (GC) content of the hydrocarbon degrading test strains was also found to be highest for the clade which harbored them. The test strains consistently occupied a distinct terminal end within the phylogenetic tree constructed by ML analysis. This study reveals that the refinery sludge imposed environmental stress on the bacterial strains which possibly caused significant genetic alteration and phenotypic adaptation. Due to the divergent evolution of the Pseudomonas and Bacillus strains in the sludge, they appeared distinctly different from other hydrocarbon degrading strains of the same genus.

Numtogenesis as a mechanism for development of cancer

Transfer of genetic material from cytoplasmic organelles to the nucleus, an ongoing process, has implications in evolution, aging, and human pathologies such as cancer. The transferred mitochondrial DNA (mtDNA) fragments in the nuclear genome are called nuclear mtDNA or NUMTs. We have named the process numtogenesis, defining the term as the transfer of mtDNA into the nuclear genome, or, less specifically, the transfer of mitochondria or mitochondrial components into the nucleus. There is increasing evidence of the involvement of NUMTs in human biology and pathology. Although information pertaining to NUMTs and numtogenesis is sparse, the role of this aspect of mitochondrial biology to human cancers is apparent. In this review, we present available knowledge about the origin and mechanisms of numtogenesis, with special emphasis on the role of NUMTs in human malignancies. We describe studies undertaken in our laboratory and in others and discuss the influence of NUMTs in tumor initiation and progression and in survival of cancer patients. We describe suppressors of numtogenesis and evolutionary conserved mechanisms underlying numtogenesis in cancer. An understanding the emerging field of numtogenesis should allow comprehension of this process in various malignancies and other diseases and, more generally, in human health.

Analysis of the complete mitochondrial genome and characterization of diverse NUMTs of Macaca leonina

As a non-human primate, the pig-tailed macaque has received wide attention because it can be infected by HIV-1. In this study, we determined the complete mtDNA sequence of the northern pig-tailed macaque (Macaca leonina). Unexpectedly, during the amplification of the mtDNA control region (D-loop region) we observed several D-loop-like sequences, which were NUMTs (nuclear mitochondrial sequences) and a total of 14 D-loop-like NUMT haplotypes were later identified in five individuals. The neighbor-joining tree and estimated divergence time based on these D-loop-like NUMT sequences of M. leonina provide some insights into the understanding of the evolutionary history of NUMTs. D-loop-like haplotypes G and H, which also exist in the nuclear genome of mulatta, appear to have been translocated into the nuclear genome before the divergence of M. mulatta and M. leonina. The other D-loop-like NUMT haplotypes were translocated into the nuclear genome of M. leonina after the divergence of the two species. Later sequence conversion was predicted to occur among these 14 D-loop-like NUMT haplotypes. The overall structure of the mtDNA of M. leonina was found to be similar to that seen in other mammalian mitochondrial genomes. Phylogenetic analysis based on the maximum likelihood method shows M. leonina clustered with Macaca silenus among the analyzed mammalian species.

Rampant nuclear insertion of mtDNA across diverse lineages within Orthoptera (Insecta)

Nuclear mitochondrial pseudogenes (numts) are non-functional fragments of mtDNA inserted into the nuclear genome. Numts are prevalent across eukaryotes and a positive correlation is known to exist between the number of numts and the genome size. Most numt surveys have relied on model organisms with fully sequenced nuclear genomes, but such analyses have limited utilities for making a generalization about the patterns of numt accumulation for any given clade. Among insects, the order Orthoptera is known to have the largest nuclear genome and it is also reported to include several species with a large number of numts. In this study, we use Orthoptera as a case study to document the diversity and abundance of numts by generating numts of three mitochondrial loci across 28 orthopteran families, representing the phylogenetic diversity of the order. We discover that numts are rampant in all lineages, but there is no discernable and consistent pattern of numt accumulation among different lineages. Likewise, we do not find any evidence that a certain mitochondrial gene is more prone to nuclear insertion than others. We also find that numt insertion must have occurred continuously and frequently throughout the diversification of Orthoptera. Although most numts are the result of recent nuclear insertion, we find evidence of very ancient numt insertion shared by highly divergent families dating back to the Jurassic period. Finally, we discuss several factors contributing to the extreme prevalence of numts in Orthoptera and highlight the importance of exploring the utility of numts in evolutionary studies.

Could pseudogenes be widespread in ants? Evidence of numts in the leafcutter ant Acromyrmex striatus (Roger, 1863) (Formicidae: Attini)

The incorporation of fragments of mitochondrial DNA (mtDNA) in the nuclear genome, known as numts (nuclear mitochondrial pseudogenes), undermines general assumptions concerning the use of mtDNA in phylogenetic and phylogeographic studies. Accidental amplifications of these nuclear copies instead of the mitochondrial target can lead to crucial misinterpretations, thus the correct identification of numts and their differentiation from true mitochondrial sequences are important in preventing this kind of error. Our goal was to describe the existence of cytochrome b (cytb) numts in the leafcutter ant Acromyrmex striatus (Roger, 1863). PCR products were directly sequenced using a pair of universal primers designed to amplify the cytb gene of these insects. Other species of leafcutter ants were also sequenced. The sequences were analyzed and the numts were identified by the presence of double peaks, indels and premature stop codons. Only A. striatus clearly showed the presence of numts, while the other species displayed the expected amplification of the mtDNA cytb gene target using the same primer pair. We hope that our report will highlight the benefits and challenges of using mtDNA in the molecular phylogenetic reconstruction and phylogeographic studies of ants, while establishing the importance of numts reports for future studies.

Targeted and robust amplification of mitochondrial DNA in the presence of nuclear-encoded mitochondrial pseudogenes using Φ29 DNA polymerases

Accidental co-amplification of nuclear-encoded mitochondrial pseudogenes (numts) in mitochondrial DNA (mtDNA) analyses is a common problem in biological and medical research alike [Yao et al., J Med Genet 45:769-772, 2008; Song et al., Proc Natl Acad Sci U S A 105:13486-13491, 2008]. A wealth of strategies have been developed to evade the contamination of mtDNA data sets with numts, but time- and cost-effective protocols of general applicability are scarce [Wolff et al., PLoS One 7:e37142, 2012]. The protocol presented here closes this gap by combining serial dilution with standard and rolling circle DNA amplification. This strategy leads to selective enrichment of mtDNA and removal of nuclear DNA (nuDNA) from template solutions, enabling mtDNA analyses in a virtually numt-free environment.

The mitochondrial genome encodes abundant small noncoding RNAs

Small noncoding RNAs identified thus far are all encoded by the nuclear genome. Here, we report that the murine and human mitochondrial genomes encode thousands of small noncoding RNAs, which are predominantly derived from the sense transcripts of the mitochondrial genes (host genes), and we termed these small RNAs mitochondrial genome-encoded small RNAs (mitosRNAs). DICER inactivation affected, but did not completely abolish mitosRNA production. MitosRNAs appear to be products of currently unidentified mitochondrial ribonucleases. Overexpression of mitosRNAs enhanced expression levels of their host genes in vitro, and dysregulated mitosRNA expression was generally associated with aberrant mitochondrial gene expression in vivo. Our data demonstrate that in addition to 37 known mitochondrial genes, the mammalian mitochondrial genome also encodes abundant mitosRNAs, which may play an important regulatory role in the control of mitochondrial gene expression in the cell.

Evaluation of the authenticity of a highly novel environmental sequence from boreal forest soil using ribosomal RNA secondary structure modeling

The number of sequences from both formally described taxa and uncultured environmental DNA deposited in the International Nucleotide Sequence Databases has increased substantially over the last two decades. Although the majority of these sequences represent authentic gene copies, there is evidence of DNA artifacts in these databases as well. These include lab artifacts, such as PCR chimeras, and biological artifacts such as pseudogenes or other paralogous sequences. Sequences that fall in basal positions in phylogenetic trees and appear distant from known sequences are particularly suspect. Phylogenetic analyses suggest that a novel sequence type (NS1) found in two boreal forest soil clone libraries belongs to the fungal kingdom but does not fall unambiguously within any known phylum. We have evaluated this sequence type using an array of secondary-structure analyses. To our knowledge, such analyses have never been used on environmental ribosomal sequences. Ribosomal secondary structure was modeled for four rRNA loci (ITS1, 5.8S, ITS2, 5' LSU). These models were analyzed for the presence of conserved domains, conserved nucleotide motifs, and compensatory base changes. Minimal free energy (MFE) foldings and GC contents of sequences representing the major fungal clades, as well as NS1, were also compared. NS1 displays secondary rRNA structures consistent with other fungi and many, but not all, conserved nucleotide motifs found across eukaryotes. However, our analyses show that many other authentic sequences from basal fungi lack more of these conserved motifs than does NS1. Together our findings suggest that NS1 represents an authentic gene copy. The methods described here can be used on any rRNA-coding sequence, not just environmental fungal sequences. As new-generation sequencing methods that yield shorter sequences become more widely implemented, methods that evaluate sequence authenticity should also be more widely implemented. For fungi, the adjacent 5.8S and ITS2 loci should be prioritized. This region is not only suited to distinguishing between closely related species, but it is also more informative in terms of expected secondary structure.

Factors affecting the relative abundance of nuclear copies of mitochondrial DNA (numts) in hominoids

Although nuclear copies of mitochondrial DNA (numts) can originate from any portion of the mitochondrial genome, evidence from humans suggests that more variable parts of the mitochondrial genome, such as the mitochondrial control region (MCR), are under-represented in the nucleus. This apparent deficit might arise from the erosion of sequence identity in numts originating from rapidly evolving mitochondrial sequences. However, the extent to which mitochondrial sequence properties impacts the number of numts detected in genomic surveys has not been evaluated. In order to address this question, we: (1) conducted exhaustive BLAST searches of MCR numts in three hominoid genomes; (2) assessed numt prevalence across the four MCR sub-domains (HV1, CCD, HV2, and MCR(F)); (3) estimated their insertion rates in great apes (Hominoidea); and (4) examined the relationship between mitochondrial DNA variability and numt prevalence in sequences originating from MCR and coding regions of the mitochondrial genome. Results indicate a marked deficit of numts from HV2 and MCR(F) MCR sub-domains in all three species. These MCR sub-domains exhibited the highest proportion of variable sites and the lowest number of detected numts per mitochondrial site. Variation in MCR insertion rate between lineages was also observed with a pronounced burst in recent integrations within chimpanzees and orangutans. A deficit of numts from HV2/MCR(F) was observed regardless of age, whereas HV1 is under-represented only in older numts (>25 million years). Finally, more variable mitochondrial genes also exhibit a lower identity with nuclear copies and because of this, appear to be under-represented in human numt databases.

Highly conserved d-loop sequences in woolly mouse opossums Marmosa (Micoureus)

This study reports the occurrence of highly conserved d-loop sequences in the mitochondrial genome of the woolly mouse opossum genus Marmosa subgenus Micoureus (Mammalia, Didelphimorphia, Didelphidae). Sixty-six sequences of Marmosa (Micoureus) demerarae, Marmosa (Micoureus) constantiae, and Marmosa (Micoureus) paraguayanus were amplified using universal d-loop primers and virtually no genetic differences were detected within and among species. These sequences matched the control region of the mitochondrial marsupial genome. Analyses of qualitative aspects of these sequences revealed that their structural composition is very similar to the d-loop region of other didelphid species. However, the total lack of variability has not been reported from other closely related species. The data analyzed here support the occurrence of highly conserved d-loop sequences, and we found no support for the hypothesis that these sequences are d-loop-like nuclear pseudogenes. Furthermore, the control and flanking regions obtained with different primers corroborate the lack of variability of the d-loop sequences in the mitochondrial genome of Marmosa (Micoureus).

Role of intercompartmental DNA transfer in producing genetic diversity

In eukaryotic cells, genes are found in three compartments-the nucleus, mitochondria, and plastids-and extensive gene transfer has occurred between them. Most organellar genes in the nucleus migrated there long ago, but transfer is ongoing and ubiquitous. It now generates mostly noncoding nuclear DNA, can also disrupt gene functions, and reshape genes by adding novel exons. Plastid or nuclear sequences have also contributed to the formation of mitochondrial tRNA genes. It is now clear that organelle-to-nucleus DNA transfer involves the escape of DNA molecules from the organelles at times of stress or at certain developmental stages, and their subsequent incorporation at sites of double-stranded breaks in nuclear DNA by nonhomologous recombination. Intercompartmental DNA transfer thus appears to be an inescapable phenomenon that has had a broad impact on eukaryotic evolution, affecting DNA repair, gene and genome evolution, and redirecting proteins to different target compartments.

Molecular poltergeists: mitochondrial DNA copies (numts) in sequenced nuclear genomes

The natural transfer of DNA from mitochondria to the nucleus generates nuclear copies of mitochondrial DNA (numts) and is an ongoing evolutionary process, as genome sequences attest. In humans, five different numts cause genetic disease and a dozen human loci are polymorphic for the presence of numts, underscoring the rapid rate at which mitochondrial sequences reach the nucleus over evolutionary time. In the laboratory and in nature, numts enter the nuclear DNA via non-homolgous end joining (NHEJ) at double-strand breaks (DSBs). The frequency of numt insertions among 85 sequenced eukaryotic genomes reveal that numt content is strongly correlated with genome size, suggesting that the numt insertion rate might be limited by DSB frequency. Polymorphic numts in humans link maternally inherited mitochondrial genotypes to nuclear DNA haplotypes during the past, offering new opportunities to associate nuclear markers with mitochondrial markers back in time.

The mitochondrial genome of the 'twisted-wing parasite' Mengenilla australiensis (Insecta, Strepsiptera): a comparative study

BACKGROUND

Strepsiptera are an unusual group of sexually dimorphic, entomophagous parasitoids whose evolutionary origins remain elusive. The lineage leading to Mengenilla australiensis (Family Mengenillidae) is the sister group to all remaining extant strepsipterans. It is unique in that members of this family have retained a less derived condition, where females are free-living from pupation onwards, and are structurally much less simplified. We sequenced almost the entire mitochondrial genome of M. australiensis as an important comparative data point to the already available genome of its distant relative Xenos vesparum (Family Xenidae). This study represents the first in-depth comparative mitochondrial genomic analysis of Strepsiptera.

RESULTS

The partial genome of M. australiensis is presented as a 13421 bp fragment, across which all 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes and 18 transfer RNA (tRNA) sequences are identified. Two tRNA translocations disrupt an otherwise ancestral insect mitochondrial genome order. A+T content is measured at 84.3%, C-content is also very skewed. Compared with M. australiensis, codon bias in X. vesparum is more balanced. Interestingly, the size of the protein coding genome is truncated in both strepsipterans, especially in X. vesparum which, uniquely, has 4.3% fewer amino acids than the average holometabolan complement. A revised assessment of mitochondrial rRNA secondary structure based on comparative structural considerations is presented for M. australiensis and X. vesparum.

CONCLUSIONS

The mitochondrial genome of X. vesparum has undergone a series of alterations which are probably related to an extremely derived lifestyle. Although M. australiensis shares some of these attributes; it has retained greater signal from the hypothetical most recent common ancestor (MRCA) of Strepsiptera, inviting the possibility that a shift in the mitochondrial selective environment might be related to the specialization accompanying the evolution of a small, morphologically simplified completely host-dependent lifestyle. These results provide useful insights into the nature of the evolutionary transitions that accompanied the emergence of Strepsiptera, but we emphasize the need for adequate sampling across the order in future investigations concerning the extraordinary developmental and evolutionary origins of this group.

Comparative analysis of mitochondrial fragments transferred to the nucleus in vertebrate

Mitochondrial DNA sequences transferred to the nucleus give rise to the socalled nuclear mitochondrial DNA (numt). In the GenBank database, 244 numts have been found in six orders of birds (Anseriformes, Columbiformes, Falconiformes, Charadriiformes, Galliformes and Passeriformes). Sequences alignment (NCBI-BLASTN) was carried out with mitochondrial and corresponding nuclear genome sequences in nine vertebrate species. The sequences with high homology were considered as numts. The number of numts ranged from 15 in chicken to 159 in chimpanzee. The sequences of numts in macaque, chimpanzee, and human spanned 100% of the entire mammalian mitochondrial genome. The reconstructed frequency of the mitochondrial gene transferred to the nucleus demonstrated that the rRNA genes had high frequencies than other mitochondrial genes. Using the RepeatMasker program, the transposable elements were detected in the flanking regions of each numt. The results showed that less than 5% of the flanking sequences were made up of repetitive elements in chicken. The GC content of 5'- and 3'-flanking regions of numts in nine species was less than 44%. The analysis of the flanking sequences provided a valuable understanding for future study on mechanism of mitochondrial gene transfer to the nucleus and the site of numt integration.

DNA transfer from organelles to the nucleus: the idiosyncratic genetics of endosymbiosis

In eukaryotes, DNA is exchanged between endosymbiosis-derived compartments (mitochondria and chloroplasts) and the nucleus. Organelle-to-nucleus DNA transfer involves repair of double-stranded breaks by nonhomologous end-joining, and resulted during early organelle evolution in massive relocation of organelle genes to the nucleus. A large fraction of the products of the nuclear genes so acquired are retargeted to their ancestral compartment; many others now function in new subcellular locations. Almost all present-day nuclear transfers of mitochondrial or plastid DNA give rise to noncoding sequences, dubbed nuclear mitochondrial DNAs (NUMTs) and nuclear plastid DNAs (NUPTs). Some of these sequences were recruited as exons, thus introducing new coding sequences into preexisting nuclear genes by a novel mechanism. In organisms derived from secondary or tertiary endosymbiosis, serial gene transfers involving nucleus-to-nucleus migration of DNA have also occurred. Intercompartmental DNA transfer therefore represents a significant driving force for gene and genome evolution, relocating and refashioning genes and contributing to genetic diversity.

Potential gene flow in natural populations of the Drosophila ananassae species cluster inferred from a nuclear mitochondrial pseudogene

A pseudogene with 94% similarity to mitochondrial cytochrome c oxidase subunit I (COI) was identified and localized to chromosome 4 of Drosophila ananassae. Because this chromosome is believed to have reduced recombination, its history can be traced using the pseudo-COI sequence. Pseudo-COI sequences were obtained from 27 iso-female lines of six taxa belonging to the D. ananassae species cluster in which reproductive isolation is incomplete. The phylogenetic network constructed from seven recognized haplotypes (#0-#6) indicated that different taxa inhabiting the same geographic area share the haplotypes: #1 from Papua New Guinean populations of D. ananassae and pallidosa-like-Wau; #2 from Papua New Guinean populations of D. ananassae, pallidosa-like, and papuensis-like; and #4 from South Pacific populations of D. ananassae and D. pallidosa. Taxon-K has a unique haplotype (#6), and 18 mutation steps separate it from the closest haplotype, #2. We discuss the possibility of chromosome 4 introgression beyond taxon boundaries.

Pseudogenes and DNA-based diet analyses: a cautionary tale from a relatively well sampled predator-prey system

Mitochondrial ribosomal DNA is commonly used in DNA-based dietary analyses. In such studies, these sequences are generally assumed to be the only version present in DNA of the organism of interest. However, nuclear pseudogenes that display variable similarity to the mitochondrial versions are common in many taxa. The presence of nuclear pseudogenes that co-amplify with their mitochondrial paralogues can lead to several possible confounding interpretations when applied to estimating animal diet. Here, we investigate the occurrence of nuclear pseudogenes in fecal samples taken from bottlenose dolphins (Tursiops truncatus) that were assayed for prey DNA with a universal primer technique. We found pseudogenes in 13 of 15 samples and 1-5 pseudogene haplotypes per sample representing 5-100% of all amplicons produced. The proportion of amplicons that were pseudogenes and the diversity of prey DNA recovered per sample were highly variable and appear to be related to PCR cycling characteristics. This is a well-sampled system where we can reliably identify the putative pseudogenes and separate them from their mitochondrial paralogues using a number of recommended means. In many other cases, it would be virtually impossible to determine whether a putative prey sequence is actually a pseudogene derived from either the predator or prey DNA. The implications of this for DNA-based dietary studies, in general, are discussed.

Ancient DNA elucidates the controversy about the flightless island hens (Gallinula sp.) of Tristan da Cunha

A persistent controversy surrounds the flightless island hen of Tristan da Cunha, Gallinula nesiotis. Some believe that it became extinct by the end of the 19th century. Others suppose that it still inhabits Tristan. There is no consensus about Gallinula comeri, the name introduced for the flightless moorhen from the nearby island of Gough. On the basis of DNA sequencing of both recently collected and historical material, we conclude that G. nesiotis and G. comeri are different taxa, that G. nesiotis indeed became extinct, and that G. comeri now inhabits both islands. This study confirms that among gallinules seemingly radical adaptations (such as the loss of flight) can readily evolve in parallel on different islands, while conspicuous changes in other morphological characters fail to occur.

Mitochondrial Introgressions into the Nuclear Genome of the Domestic Cat

Translocation of mtDNA into the nuclear genome, also referred to as numt, was first reported in the domestic cat (Felis catus) by Lopez et al. (1994). The Lopez-numt consisted of a translocation of 7.9 kbp of mtDNA that inserted into the domestic cat chromosome D2 around 1.8 million years ago. More than a decade later, the release of the domestic cat whole-genome shotgun sequences (1.9x coverage) provides the resource to obtain more comprehensive insight into the extent of mtDNA transfer over time in the domestic cat genome. MegaBLAST searches revealed that the cat genome harbors a wide variety of numts (298 320 bp), one-third of which likely correspond to the Lopez-numt tandem repeat, whereas the remaining numts are probably derived from multiple independent insertions, which in some cases were followed by segmental duplication after insertion in the nucleus. Numts were detected across most cat chromosomes, but the number of numts assigned to chromosomes is underestimated due to the relatively high number of numt sequences with insufficient flanking sequence to map. The catalog of cat numts provides a valuable resource for future studies in Felidae species, including its use as a tool to avoid numt contaminations that may confound population genetics and phylogenetic studies.

Transition-transversion bias is not universal: a counter example from grasshopper pseudogenes

Comparisons of the DNA sequences of metazoa show an excess of transitional over transversional substitutions. Part of this bias is due to the relatively high rate of mutation of methylated cytosines to thymine. Postmutation processes also introduce a bias, particularly selection for codon-usage bias in coding regions. It is generally assumed, however, that there is a universal bias in favour of transitions over transversions, possibly as a result of the underlying chemistry of mutation. Surprisingly, this underlying trend has been evaluated only in two types of metazoan, namely Drosophila and the Mammalia. Here, we investigate a third group, and find no such bias. We characterize the point substitution spectrum in Podisma pedestris, a grasshopper species with a very large genome. The accumulation of mutations was surveyed in two pseudogene families, nuclear mitochondrial and ribosomal DNA sequences. The cytosine-guanine (CpG) dinucleotides exhibit the high transition frequencies expected of methylated sites. The transition rate at other cytosine residues is significantly lower. After accounting for this methylation effect, there is no significant difference between transition and transversion rates. These results contrast with reports from other taxa and lead us to reject the hypothesis of a universal transition/transversion bias. Instead we suggest fundamental interspecific differences in point substitution processes.

Some mutations occur more frequently than others. We need to understand these biases if we are to interpret the differences that have accumulated between species and individuals. Applications include estimating the time since evolutionary lineages diverged and detecting the signature of natural selection in DNA sequences. However, mutational biases have been obscured because, since mutations arose, natural selection has eliminated some whilst allowing others to persist to the present. We therefore study mutations that have accumulated in regions of the genome that are free from selection in a grasshopper with a gigantic genome. All other animal studies using this approach find an excess of mutations between DNA bases having similar biochemical properties (transitions rather than transversions). This bias has been widely interpreted as a consequence of the fundamental biochemical basis of mutation. However, once we exclude mutations associated with DNA methylation, we find no evidence of a transition bias, unlike the few comparable animal studies that make the same correction. We propose that this result indicates previously unanticipated differences between species in the selection on or mutation of their DNA.

Molecular analyses of mitochondrial pseudogenes within the nuclear genome of arvicoline rodents

Nuclear sequences of mitochondrial origin (numts) are common among animals and plants. The mechanism(s) by which numts transfer from the mitochondrion to the nucleus is uncertain, but their insertions may be mediated in part by chromosomal repair mechanisms. If so, then lineages where chromosomal rearrangements are common should be good models for the study of numt evolution. Arvicoline rodents are known for their karyotypic plasticity and numt pseudogenes have been discovered in this group. Here, we characterize a 4 kb numt pseudogene in the arvicoline vole Microtus rossiaemeridionalis. This sequence is among the largest numts described for a mammal lacking a completely sequenced genome. It encompasses three protein-coding and six tRNA pseudogenes that span approximately 25% of the entire mammalian mitochondrial genome. It is bordered by a dinucleotide microsatellite repeat and contains four transposable elements within its sequence and flanking regions. To determine the phylogenetic distribution of this numt among the arvicolines, we characterized one of the mitochondrial pseudogenes (cytochrome b) in 21 additional arvicoline species. Average rates of nucleotide substitution in this arvicoline pseudogene are estimated as 2.3 x 10(-8) substitutions/per site/per year. Furthermore, we performed comparative analyses among all species to estimate the age of this mitochondrial transfer at nearly 4 MYA, predating the origin of most arvicolines.

A Comparative Analysis of numt Evolution in Human and Chimpanzee

Mitochondrial DNA sequences are frequently transferred into the nuclear genome, giving rise to numts (nuclear DNA sequences of mitochondrial origin). So far, the evolutionary history of numts has largely been studied by using single genomes. Here, we present the first attempt to study numt evolution in a comparative manner by using a pairwise genomic alignment. The total number of numts was estimated to be 452 in human and 469 in chimpanzee. numts that were found in both genomes at identical loci were deemed to be orthologous; 391 numts (>80%) were classified as such. The preponderance of orthologous numts is due to the very short divergence time between the 2 hominoids. The rest of numts were deemed to be nonorthologous. Nonorthologous numts were subdivided into 1) ancestral numts that have lost an ortholog in one species through deletion (12 in human and 11 in chimpanzee), 2) new numts acquired by the insertion of a mitochondrial sequence after the divergence of the 2 species (34 in human and 46 in chimpanzee), and 3) paralogous numts created by the tandem duplication of a preexisting numt (2 in human). This approach also enabled us to reconstruct the numt repertoire in the common ancestor of humans and chimpanzees (409 numts). Our comparative approach is also useful in identifying the exact boundaries of numts.

Sheep Mitochondrial DNA Variation in European, Caucasian, and Central Asian Areas

Three distinct mitochondrial maternal lineages (haplotype Groups A, B, and C) have been found in the domestic sheep. Group B has been observed primarily in European domestic sheep. The European mouflon carries this haplotype group. This could suggest that European mouflon was independently domesticated in Europe, although archaeological evidence supports sheep domestication in the central part of the Fertile Crescent. To investigate this question, we sequenced a highly variable segment of mitochondrial DNA (mtDNA) in 406 unrelated animals from 48 breeds or local varieties. They originated from a wide area spanning northern Europe and the Balkans to the Altay Mountains in south Siberia. The sample included a representative cross-section of sheep breeds from areas close to the postulated Near Eastern domestication center and breeds from more distant northern areas. Four (A, B, C, and D) highly diverged sheep lineages were observed in Caucasus, 3 (A, B and C) in Central Asia, and 2 (A and B) in the eastern fringe of Europe, which included the area north and west of the Black Sea and the Ural Mountains. Only one example of Group D was detected. The other haplotype groups demonstrated signs of population expansion. Sequence variation within the lineages implied Group A to have expanded first. This group was the most frequent type only in Caucasian and Central Asian breeds. Expansion of Group C appeared most recently. The expansion of Group B involving Caucasian sheep took place at nearly the same time as the expansion of Group A. Group B expansion for the eastern European area started approximately 3,000 years after the earliest inferred expansion. An independent European domestication of sheep is unlikely. The distribution of Group A variation as well as other results are compatible with the Near East being the domestication site. Groups C and D may have been introgressed later into a domestic stock, but larger samples are needed to infer their geographical origin. The results suggest that some mitochondrial lineages arrived in northern Europe from the Near East across Russia.

Evolutionary analysis of a large mtDNA translocation (numt) into the nuclear genome of the Panthera genus species

Translocation of cymtDNA into the nuclear genome, also referred to as numt, has been reported in many species, including several closely related to the domestic cat (Felis catus). We describe the recent transposition of 12,536 bp of the 17 kb mitochondrial genome into the nucleus of the common ancestor of the five Panthera genus species: tiger, P. tigris; snow leopard, P. uncia; jaguar, P. onca; leopard, P. pardus; and lion, P. leo. This nuclear integration, representing 74% of the mitochondrial genome, is one of the largest to be reported in eukaryotes. The Panthera genus numt differs from the numt previously described in the Felis genus in: (1) chromosomal location (F2-telomeric region vs. D2-centromeric region), (2) gene make up (from the ND5 to the ATP8 vs. from the CR to the COII), (3) size (12.5 vs. 7.9 kb), and (4) structure (single monomer vs. tandemly repeated in Felis). These distinctions indicate that the origin of this large numt fragment in the nuclear genome of the Panthera species is an independent insertion from that of the domestic cat lineage, which has been further supported by phylogenetic analyses. The tiger cymtDNA shared around 90% sequence identity with the homologous numt sequence, suggesting an origin for the Panthera numt at around 3.5 million years ago, prior to the radiation of the five extant Panthera species.

Origin, evolution and genetic effects of nuclear insertions of organelle DNA

In eukaryotes, nuclear genomes are subject to an influx of DNA from mitochondria and plastids. The nuclear insertion of organellar sequences can occur during the illegitimate repair of double-stranded breaks. After integration, nuclear organelle DNA is modified by point mutations, and by deletions. Insertion of organelle DNA into nuclear genes is not rare and can potentially have harmful effects. In humans, some insertions of nuclear mitochondrial DNA are associated with heritable diseases. It remains to be determined whether nuclear organelle DNA can contribute beneficially to gene evolution.

Discovery of a large number of previously unrecognized mitochondrial pseudogenes in fish genomes

Nuclear inserted copies of mitochondrial origin (numts) vary widely among eukaryotes, with human and plant genomes harboring the largest repertoires. Numts were previously thought to be absent from fish species, but the recent release of three fish nuclear genome sequences provides the resource to obtain a more comprehensive insight into the extent of mtDNA transfer in fishes. From the sequence analyses of the genomes of Fugu rubripes, Tetraodon nigroviridis, and Danio rerio, we have identified 2, 5, and 10 recent numt integrations, respectively, which integrated into those genomes less than 0.6 million years (Myr) ago. Such results contradict the hypothesis of absence or rarity of numts in fishes, as (i) the ratio of numts to the total size of the nuclear genome in T. nigroviridis was superior to the ratio observed in several higher vertebrate species (e.g., chicken, mouse, and rat), and only surpassed by humans, and (ii) the mtDNA coverage transferred to the nuclear genome of D. rerio is exceeded only by human and mouse, within the whole range of eukaryotic genomes surveyed for numts. Additionally, 335, 336, and 471 old numts (>12.5 Myr) were detected in F. rubripes, T. nigroviridis, and D. rerio, respectively. Surprisingly, old numts are inserted preferentially into known or predicted genes, as inferred for recent numts in human. However, because in fish genomes such integrations are old, they are likely to represent evolutionary successes and they may be considered a potential important evolutionary mechanism for the enhancement of genomic coding regions.