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.     

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.     

Phylogeography of the banded killifish (Fundulus diaphanus): glacial races and secondary contact

Spatial patterns of morphological and genetic polymorphisms between the two banded killifish sub-species ( Fundulus diaphanus diaphanus and F. d. menona ) provide evidence that these sub-species represent distinct evolutionary lineages corresponding to glacial races of Atlantic and Mississippian origins. Individuals with intermediate phenotypes in Lake Ontario and the upper St Lawrence River probably result from secondary contacts between these races. The hybrid zone unexpectedly extends to Lake Erie and involves the introgression of the F. d. diaphanus mitochondrial genome into fish with morphology like that of F. d. menona .     

Comparative analysis of mitochondrial fragments transferred to the nucleus in vertebrate

Mitochondfial DNA sequences transferred to the nucleus give rise to the so-called 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.     

Mitochondrial DNA insertions in the nuclear horse genome

The insertion of mitochondrial DNA in the nuclear genome generates numts, nuclear sequences of mitochondrial origin. In the horse reference genome, we identified 82 numts and showed that the entire horse mitochondrial DNA is represented as numts without gross bias. Numts were inserted in the horse nuclear genome at random sites and were probably generated during the repair of DNA double-strand breaks. We then analysed 12 numt loci in 20 unrelated horses and found that null alleles, lacking the mitochondrial DNA insertion, were present at six of these loci. At some loci, the null allele is prevalent in the sample analysed, suggesting that, in the horse population, the number of numt loci may be higher than 82 present in the reference genome. Contrary to humans, the insertion polymorphism of numts is extremely frequent in the horse population, supporting the hypothesis that the genome of this species is in a rapidly evolving state.     

Rates of DNA Duplication and Mitochondrial DNA Insertion in the Human Genome

The hundreds of mitochondrial pseudogenes in the human nuclear genome sequence (numts) constitute an excellent system for studying and dating DNA duplications and insertions. These pseudogenes are associated with many complete mitochondrial genome sequences and through those with a good fossil record. By comparing individual numts with primate and other mammalian mitochondrial genome sequences, we estimate that these numts arose continuously over the last 58 million years. Our pairwise comparisons between numts suggest that most human numts arose from different mitochondrial insertion events and not by DNA duplication within the nuclear genome. The nuclear genome appears to accumulate mtDNA insertions at a rate high enough to predict within-population polymorphism for the presence/absence of many recent mtDNA insertions. Pairwise analysis of numts and their flanking DNA produces an estimate for the DNA duplication rate in humans of 2.2 × 10^–9 per numt per year. Thus, a nucleotide site is about as likely to be involved in a duplication event as it is to change by point substitution. This estimate of the rate of DNA duplication of noncoding DNA is based on sequences that are not in duplication hotspots, and is close to the rate reported for functional genes in other species.     

Genetic subdivision, glacial refugia and postglacial recolonization in the golden-striped salamander, Chioglossa lusitanica (Amphibia: urodela)

The golden-striped salamander (Chioglossa lusitanica) is an ecologically specialized species, endemic to north-western Iberia. Patterns of genetic variation were assessed at seven polymorphic enzyme loci and one mitochondrial DNA (mtDNA) marker (cytochrome b) in 17 populations across its range. Estimates of enzyme genetic diversity revealed a high degree of genetic subdivision (FST = 0.68), mainly attributable to the existence of two groups of populations. The groups were located, respectively, north and south of the Mondego River, indicating that this river coincided with a major historical barrier to gene flow. A significant decrease in genetic variability from the Mondego northwards was associated with the Douro and Minho rivers. mtDNA sequence variation revealed a congruent pattern of two haplotype groups (d = 2.2%), with a geographical distribution resembling that of allozymes. The pattern and depth of genetic variation is consistent with the following hypotheses: (i) subdivision of an ancestral range of the species prior to the middle Pleistocene; (ii) secondary contact between populations representing historical refugia; (iii) relatively recent range expansion giving rise to the northern part of the species range; and (iv) loss of genetic variation through founder effects during range expansion across major rivers.     

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 haplotype variability of brown trout populations from Northwestern Iberian Peninsula,a secondary contact area between lineages

Previous studies using mitochondrial markers and microsatellites revealed a new highly differentiated Atlantic lineage of brown trout in the Iberian Peninsula mainly confined to the Upper Duero Basin. Data suggest that Duero lineage (DU) hybridizes with previously described Atlantic (AT) lineages in that basin. However, nuclear and cytogenetic data suggest another hybrid zone between divergent lineages might be present in the Miño Basin, another Atlantic drainage located to the north. Here, we extend those studies by analyzing the mitochondrial haplotype variability of brown trout populations from northwestern Spain, a poorly studied area, to obtain an accurate picture of the distribution of DU. Results show that the DU haplotype is fixed in several isolated populations from the upper course of the Miño River and its main tributary the Sil River. DU haplotype co-occurs with two related AT haplotypes in the Miño Basin revealing a new secondary contact zone. The identification of both longtime stable populations of DU and several other populations showing introgression with widespread AT lineages will be useful for the management and conservation of such an endemic lineage.     

Assessing the effects of primer specificity on eliminating numt coamplification in DNA barcoding: a case study from Orthoptera (Arthropoda: Insecta)

DNA barcoding is a diagnostic method of species identification based on sequencing a short mitochondrial DNA fragment of cytochrome oxidase I (COI), but its ability to correctly diagnose species is limited by the presence of nuclear mitochondrial pseudogenes (numts). Numts can be coamplified with the mitochondrial orthologue when using universal primers, which can lead to incorrect species identification and an overestimation of the number of species. Some researchers have proposed that using more specific primers may help eliminate numt coamplification, but the efficacy of this method has not been thoroughly tested. In this study, we investigate the taxonomic distribution of numts in 11 lineages within the insect order Orthoptera, by analysing cloned COI sequences and further test the effects of primer specificity on eliminating numt coamplification in four lineages. We find that numts are coamplified in all 11 taxa using universal (barcoding) primers, which suggests that numts may be widespread in other taxonomic groups as well. Increased primer specificity is only effective at reducing numt coamplification in some species tested, and only eliminates it in one species tested. Furthermore, we find that a number of numts do not have stop codons or indels, making it difficult to distinguish them from mitochondrial orthologues, thus putting the efficacy of barcoding quality control measures under question. Our findings suggest that numt coamplification is a serious problem for DNA barcoding and more quality control measures should be implemented to identify and eliminate numts prior to using mitochondrial barcodes for species diagnoses.     

Phylogeography of the rare Balkan endemic Martino's vole, Dinaromys bogdanovi, reveals strong differentiation within the western Balkan Peninsula

The spatial genetic structure of Martino's vole, a rare palaeoendemic species of the western Balkans, was investigated using DNA isolated from archived museum samples. The study was based on partial sequencing (555 bp) of the mitochondrial cytochrome b gene for 63 specimens from 20 different localities throughout the species' range. Three highly divergent allopatric phylogenetic lineages (Northwestern, Central and Southeastern) were recognized among 47 haplotypes, suggesting three independent glacial differentiation centres within the western Balkans. The Northwestern lineage, which showed the highest divergence from all other samples (mean sequence divergence of 6.64% +/- 1.10), comprised samples collected from northwest of the Neretva River in Croatia, western Bosnia and Herzegovina. The Central and Southeastern lineages were separated by an average sequence divergence of 2.95% +/- 0.66 and were geographically divided by the Drim River (the Kosovo Basin in Serbia). Overall, haplotype diversity decreased from the Northwestern lineage to the Central and subsequently the Southeastern lineage, in a geographical pattern consistent with a stepping stone colonization. The observed distribution indicates a gradual southerly expansion with subsequent allopatry across the Neretva River and Drim River approximately 1 and 0.3 million years ago, respectively. Such a scenario is concordant with palaeontological evidence. Several highly divergent sublineages within the Northwestern and Central lineages showed no significant geographical structuring, suggesting secondary contact of allopatrically evolved lineages. We hypothesize that the topographical complexity of the Balkans promoted allopatry and isolation on a small geographical scale during interglacial periods, with secondary contact during glacial maxima. Furthermore, the three main lineages should be regarded as evolutionary significant units with important implications for conservation. Ecological data show that the Northwestern lineage in particular fulfils all criteria for a highly endangered, evolutionarily significant unit.     

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 ∼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 × 10⁻⁸ 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. All sequences generated in this study have been deposited within the GenBank database.     

Few mitochondrial DNA sequences are inserted into the turkey (Meleagris gallopavo) nuclear genome: evolutionary analyses and informativity in the domestic lineage

Mitochondrial DNA (mtDNA) insertions have been detected in the nuclear genome of many eukaryotes. These sequences are pseudogenes originated by horizontal transfer of mtDNA fragments into the nuclear genome, producing nuclear DNA sequences of mitochondrial origin (numt). In this study we determined the frequency and distribution of mtDNA‐originated pseudogenes in the turkey (Meleagris gallopavo) nuclear genome. The turkey reference genome (Turkey_2.01) was aligned with the reference linearized mtDNA sequence using last . A total of 32 numt sequences (corresponding to 18 numt regions derived by unique insertional events) were identified in the turkey nuclear genome (size ranging from 66 to 1415 bp; identity against the modern turkey mtDNA corresponding region ranging from 62% to 100%). Numts were distributed in nine chromosomes and in one scaffold. They derived from parts of 10 mtDNA protein‐coding genes, ribosomal genes, the control region and 10 tRNA genes. Seven numt regions reported in the turkey genome were identified in orthologues positions in the Gallus gallus genome and therefore were present in the ancestral genome that in the Cretaceous originated the lineages of the modern crown Galliformes. Five recently integrated turkey numts were validated by PCR in 168 turkeys of six different domestic populations. None of the analysed numts were polymorphic (i.e. absence of the inserted sequence, as reported in numts of recent integration in other species), suggesting that the reticulate speciation model is not useful for explaining the origin of the domesticated turkey lineage.     

Reduced introgression of the Y chromosome between subspecies of the European rabbit (Oryctolagus cuniculus) in the Iberian Peninsula

The role of the Y chromosome in speciation is unclear. Hybrid zones provide natural arenas for studying speciation, as differential introgression of markers may reveal selection acting against incompatibilities. Two subspecies of the European rabbit (Oryctolagus cuniculus) form a hybrid zone in the Iberian Peninsula. Previous work on mitochondrial DNA (mtDNA), Y- and X-linked loci revealed the existence of two divergent lineages in the rabbit genome and that these lineages are largely subspecies-specific for mtDNA and two X-linked loci. Here we investigated the geographic distribution of the two Y chromosome lineages by genotyping two diagnostic single nucleotide polymorphisms in a sample of 353 male rabbits representing both subspecies, and found that Y chromosome lineages are also largely subspecies-specific. We then sequenced three autosomal loci and discovered considerable variation in levels of differentiation at these loci. Finally, we compared estimates of population differentiation between rabbit subspecies at 26 markers and found a surprising bimodal distribution of F(ST)values. The vast majority of loci showed little or no differentiation between rabbit subspecies while a few loci, including the SRY gene, showed little or no introgression across the hybrid zone. Estimates of population differentiation for the Y chromosome were surprisingly high given that there is male-biased dispersal in rabbits. Taken together, these data indicate that there is a clear dichotomy in the rabbit genome and that some loci remain highly differentiated despite extensive gene flow following secondary contact.     

Analysis of nuclear copies of mitochondrial sequences in honeybee (Apis mellifera) genome

At least 0.08% of the Apis mellifera nuclear genome contains sequences that originated from mitochondria. These nuclear copies of mitochondrial sequences (numts) are scattered all over the honeybee chromosomes and have originated by multiple independent insertions of mitochondrial DNA (mtDNA) as evident by phylogenetic analysis. Apart from original insertions, moderate duplications of numts also contributed to the present pattern and distribution of mitochondrial sequences in honeybee chromosomes. Assimilation of mitochondrial genes in the nuclear genome is mediated by extensive fragmentations of the original inserts. Replication slippage seems to be a major mechanism by which small sequences are inserted or deleted from mtDNA destined to nucleus. Most of the honeybee numts (84%) are located in the nongenic regions. The majority (94%) of the numts that are located in predicted nuclear genes have originated from mitochondrial genes coding for cytochrome oxidase and NADH dehydrogenase subunits. On the other hand, the mitochondrial rRNA or tRNA gene sequences are predominantly (88%) located in nongenic regions of the genome. Evidences also support for exertion of purifying selection on numts located in specific genes. Comparative analysis of numts of European, African, and Africanized honeybees suggests that numt evolution in A. mellifera is probably not demarked by speciation time frame but may be a continuous and dynamic process.     

Multilocus phylogeography of a widespread savanna–woodland‐adapted rodent reveals the influence of Pleistocene geomorphology and climate change in Africa's Zambezi region

Understanding historical influences of climate and physiographic barriers in shaping patterns of biodiversity remains limited for many regions of the world. For mammals of continental Africa, phylogeographic studies, particularly for West African lineages, implicate both geographic barriers and climate oscillations in shaping small mammal diversity. In contrast, studies for southern African species have revealed conflicting phylogenetic patterns for how mammalian lineages respond to both climate change and geologic events such as river formation, especially during the Pleistocene. However, these studies were often biased by limited geographic sampling or exclusively focused on large‐bodied taxa. We exploited the broad southern African distribution of a savanna–woodland‐adapted African rodent, Gerbilliscus leucogaster (bushveld gerbil) and generated mitochondrial, autosomal and sex chromosome data to quantify regional signatures of climatic and vicariant biogeographic phenomena. Results indicate the most recent common ancestor for all G. leucogaster lineages occurred during the early Pleistocene. We documented six divergent mitochondrial lineages that diverged ~0.270–0.100 mya, each of which was geographically isolated during periods characterized by alterations to the course of the Zambezi River and its tributaries as well as regional ‘megadroughts’. Results demonstrate the presence of a widespread lineage exhibiting demographic expansion ~0.065–0.035 mya, a time that coincides with savanna–woodland expansion across southern Africa. A multilocus autosomal perspective revealed the influence of the Kafue River as a current barrier to gene flow and regions of secondary contact among divergent mitochondrial lineages. Our results demonstrate the importance of both climatic fluctuations and physiographic vicariance in shaping the distribution of southern African biodiversity.     

Introgression at differentially aged hybrid zones in red-tailed chipmunks

Hybrid zones allow us to investigate the maintenance and the break down of reproductive isolation; they are a window into the speciation process. Tamias ruficaudus (red-tailed chipmunk) has a roughly ring-like distribution in the Inland Northwest and includes two morphologically well-differentiated subspecies, T. r. ruficaudus (in the eastern portion of its range) and T. r. simulans (in the western portion). These taxa meet at two contact zones: the Lochsa River in Idaho and 200 km to the north, near Whitefish, Montana. The Lochsa Zone is encompassed within the Clearwater River Drainage, which has been proposed as a glacial refugium for many taxa throughout the Pleistocene, whereas the Whitefish Zone was under the Cordilleran ice sheet during the most recent glacial maxima approximately 10,000 years ago. Mitochondrial DNA introgression has been documented at both contact zones, yet the subspecies remain significantly distinct with respect to bacular morphology and no intermediate morphologies have ever been reported. Here, we elucidate differentiation and introgression of the nuclear genome using ten microsatellite loci and compare findings to previously described mitochondrial DNA haplotype distribution and introgression. We found significant substructure in the nuclear data; each subspecies is divided into at least two genetically distinct demes. At the Lochsa contact zone, individuals restricted to the mtDNA zone of introgression form a distinct deme at microsatellite loci whereas in the younger, Whitefish contact zone, there is no hybrid-zone specific group. The genetic distances of the demes within these two subspecies indicate recent northward expansion.     

Preventing the pollution of mitochondrial datasets with nuclear mitochondrial paralogs (numts)

Molecular tools have become prominent in ecology and evolution. A target of choice for molecular ecologists and evolutionists is mitochondrial DNA (mtDNA), whose many advantages have also convinced broad-scale, pragmatic programmes such as barcode initiatives. Of course, mtDNA is also of interest to human geneticists investigating mitochondrial diseases. Studies using mtDNA are however put at great risk by the inadvertent co-amplification or preferred amplification of nuclear pseudogenes (numts). A posteriori analysis of putative mtDNA sequences can help in removing numts but faces severe limitations (e.g. recently translocated numts will most of the time go unnoticed). Counter-measures taken a priori, i.e. explicitly designed for avoiding numt co-amplification or preferred amplification, are appealing but have never been properly assessed. Here we investigate the efficiency of four such measures (mtDNA enrichment, cDNA amplification, long-range amplification and pre-PCR dilution) on a common set of numt cases, showing that mtDNA enrichment is the worst performer while the use of pre-PCR dilution is a simple, yet robust method to prevent the pollution of putative mtDNA datasets with numts. Therefore, straightforward recommendations can be made that, if followed, will considerably increase the confidence in the mitochondrial origin of any mtDNA-like sequence.     

Mito‐nuclear discordance across a recent contact zone for California voles

To examine the processes that maintain genetic diversity among closely related taxa, we investigated the dynamics of introgression across a contact zone between two lineages of California voles (Microtus californicus). We tested the prediction that introgression of nuclear loci would be greater than that for mitochondrial loci, assuming ongoing gene flow across the contact zone. We also predicted that genomic markers would show a mosaic pattern of differentiation across this zone, consistent with genomes that are semi‐permeable. Using mitochondrial cytochrome b sequences and genome‐wide loci developed via ddRAD‐seq, we analyzed genetic variation for 10 vole populations distributed along the central California coast; this transect included populations from within the distributions of both parental lineages as well as the putative contact zone. Our analyses revealed that (1) the two lineages examined are relatively young, having diverged ca. 8.5–54 kya, (2) voles from the contact zone in Santa Barbara County did not include F1 or early generation backcrossed individuals, and (3) there appeared to be little to no recurrent gene flow across the contact zone. Introgression patterns for mitochondrial and nuclear markers were not concordant; only mitochondrial markers revealed evidence of introgression, putatively due to historical hybridization. These differences in genetic signatures are intriguing given that the contact zone occurs in a region of continuous vole habitat, with no evidence of past or present physical barriers. Future studies that examine specific isolating mechanisms, such as microhabitat use and mate choice, will facilitate our understanding of how genetic boundaries are maintained in this system.     

Nuclear insertions help and hinder inference of the evolutionary history of gorilla mtDNA

Numts are fragments of mitochondrial DNA (mtDNA) that have been translocated to the nucleus, where they can persist while their mitochondrial counterparts continue to rapidly evolve. Thus, numts represent 'molecular fossils' useful for comparison with mitochondrial variation, and are particularly suited for studies of the fast-evolving hypervariable segment of the mitochondrial control region (HV1). Here we used information from numts found in western gorillas (Gorilla gorilla) and eastern gorillas (Gorilla beringei) to estimate that these two species diverged about 1.3 million years ago (Ma), an estimate similar to recent calculations for the divergence of chimpanzee and bonobo. We also describe the sequence of a gorilla numt still possessing a segment lost from all contemporary gorilla mtDNAs. In contrast to that sequence, many numts of the HV1 are highly similar to authentic mitochondrial organellar sequences, making it difficult to determine whether purported mitochondrial sequences truly derive from that genome. We used all available organellar HV1 and corresponding numt sequences from gorillas in a phylogenetic analysis aimed at distinguishing these two types of sequences. Numts were found in several clades in the tree. This, in combination with the fact that only a limited amount of the extant variation in gorillas has been sampled, suggests that categorization of new sequences by the indirect means of phylogenetic comparison would be prone to uncertainty. We conclude that for taxa such as gorillas that contain numerous numts, direct approaches to the authentication of HV1 sequences, such as amplification strategies relying upon the circularity of the mtDNA molecule, remain necessary.