Mitochondrial pseudogenes: evolution's misplaced witnesses

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

青蟹线粒体COI假基因的分离和特征分析

线粒体DNA标记在遗传结构和系统进化研究中得到广泛应用,然而核假基因的存在对此有很大威胁。本文以中国东南沿海的青蟹(Scylla paramamosain)为研究对象,利用线粒体COI基因的通用引物和特异性引物进行扩增,分别得到34个假基因（nuclear mitochondrial pseudogenes, Numts）和5个线粒体COI基因序列。在所获得的34个假基因中共定义了29种单倍型,根据序列的相似度,这些假基因可以分为2类,每类假基因都有各自保守的核苷酸序列。第Ⅰ类假基因存在2处插入序列和1处8 bp的缺失序列,这些位点导致了整个阅读框的移位;在第Ⅱ类假基因和5个线粒体COI序列中只有碱基替换,未发现插入和缺失序列。实验结果分析表明,这两类假基因分别代表了2次核整合事件,即核转移事件的最低值。研究结果提示了     

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.     

Structure and chromosomal distribution of human mitochondrial pseudogenes

Nuclear mitochondrial pseudogenes (Numts) have been found in the genome of many eukaryote species, including humans. Using a BLAST approach, we found 1105 DNA sequences homologous to mitochondrial DNA (mtDNA) in the August 2001 Goldenpath human genome database. We assembled these sequences manually into 286 pseudogenes on the basis of single insertion events and constructed a chromosomal map of these Numts. Some pseudogenes appeared highly modified, containing inversions, deletions, duplications, and displaced sequences. In the case of four randomly selected Numts, we used PCR tests on cells lacking mtDNA to ensure that our technique was free from genome-sequencing artifacts. Furthermore, phylogenetic investigation suggested that one Numt, apparently inserted into the nuclear genome 25-30 million years ago, had been duplicated at least 10 times in various chromosomes during the course of evolution. Thus, these pseudogenes should be very useful in the study of ancient mtDNA and nuclear genome evolution.     

Mitochondrial DNA Restriction Map and Cytochrome c Oxidase Subunits I and II Sequence Divergence of Corn Stalk Borer Sesamia nonagrioides (Lepidoptera: Noctuidae)

Corn stalk borer Sesamia nonagrioides (Lefebvre) (Lepidoptera: Noctuidae) is among the most important insect pests of corn in the Mediterranean basin. The mitochondrial DNA of this insect was purified and a restriction map was constructed. The size of the mtDNA genome is 16.3 kb. Genetic analysis of four corn stalk borer populations, collected from Greece (three populations) and Spain (one population), was undertaken using DNA sequences of the mtDNA cytochrome oxidase (CO) I and II genes. Sequencing of a 2079 bp region of these genes revealed 25 polymorphic sites among the populations. Five molecular RFLP markers, located in the mtDNA COI and COII genes, were surveyed, and two different haplotypes were detected. Phylogenetic analysis based on COI/COII nucleotide sequences revealed genetic differentiation between samples, and the results are discussed in relation to the geographic distribution of the corn stalk borer in two Mediterranean countries.     

Phylogenetic discordance at the species boundary: comparative gene genealogies among rapidly radiating Heliconius butterflies

Recent adaptive radiations provide excellent model systems for understanding speciation, but rapid diversification can cause problems for phylogenetic inference. Here we use gene genealogies to investigate the phylogeny of recent speciation in the heliconiine butterflies. We sequenced three gene regions, intron 3 ( approximately 550 bp) of sex-linked triose-phosphate isomerase (Tpi), intron 3 ( approximately 450 bp) of autosomal mannose-phosphate isomerase (Mpi), and 1,603 bp of mitochondrial cytochrome oxidase subunits I and II (COI and COII), for 37 individuals from 25 species of Heliconius and related genera. The nuclear intron sequences evolved at rates similar to those of mitochondrial coding sequences, but the phylogenetic utility of introns was restricted to closely related geographic populations and species due to high levels of indel variation. For two sister species pairs, Heliconius erato-Heliconius himera and Heliconius melpomene-Heliconius cydno, there was highly significant discordance between the three genes. At mtDNA and Tpi, the hypotheses of reciprocal monophyly and paraphyly of at least one species with respect to its sister could not be distinguished. In contrast alleles sampled from the third locus, Mpi, showed polyphyletic relationships between both species pairs. In all cases, recent coalescence of mtDNA lineages within species suggests that polyphyly of nuclear genes is not unexpected. In addition, very similar alleles were shared between melpomene and cydno, implying recent gene flow. Our finding of discordant genealogies between genes is consistent with models of adaptive speciation with ongoing gene flow and highlights the need for multiple locus comparisons to resolve phylogeny among closely related species.     

Mitochondrial DNA restriction map and cytochrome c oxidase subunits I and II sequence divergence of corn stalk borer Sesamia nonagrioides (Lepidoptera: Noctuidae)

Corn stalk borer Sesamia nonagrioides (Lefebvre) (Lepidoptera: Noctuidae) is among the most important insect pests of corn in the Mediterranean basin. The mitochondrial DNA of this insect was purified and a restriction map was constructed. The size of the mtDNA genome is 16.3 kb. Genetic analysis of four corn stalk borer populations, collected from Greece (three populations) and Spain (one population), was undertaken using DNA sequences of the mtDNA cytochrome oxidase (CO) I and II genes. Sequencing of a 2079 bp region of these genes revealed 25 polymorphic sites among the populations. Five molecular RFLP markers, located in the mtDNA COI and COII genes, were surveyed, and two different haplotypes were detected. Phylogenetic analysis based on COI/COII nucleotide sequences revealed genetic differentiation between samples, and the results are discussed in relation to the geographic distribution of the corn stalk borer in two Mediterranean countries.     

Frequent assimilation of mitochondrial DNA by grasshopper nuclear genomes

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

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

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

Highly conserved D-loop-like nuclear mitochondrial sequences (Numts) in tiger (Panthera tigris)

Using oligonucleotide primers designed to match hypervariable segments I (HVS-1) ofPanthera tigris mitochondrial DNA (mtDNA), we amplified two different PCR products (500 bp and 287 bp) in the tiger (Panthera tigris), but got only one PCR product (287 bp) in the leopard (Panthera pardus). Sequence analyses indicated that the sequence of 287 bp was a D-loop-like nuclear mitochondrial sequence (Numts), indicating a nuclear transfer that occurred approximately 4.8–17 million years ago in the tiger and 4.6–16 million years ago in the leopard. Although the mtDNA D-loop sequence has a rapid rate of evolution, the 287-bp Numts are highly conserved; they are nearly identical in tiger subspecies and only 1.742% different between tiger and leopard. Thus, such sequences represent molecular ‘fossils’ that can shed light on evolution of the mitochondrial genome and may be the most appropriate outgroup for phylogenetic analysis. This is also proved by comparing the phylogenetic trees reconstructed using the D-loop sequence of snow leopard and the 287-bp Numts as outgroup.     

A cryptic, intergeneric cytochrome c oxidase I pseudogene in tyrant flycatchers (family: Tyrannidae)

Nuclear mitochondrial pseudogenes, or "numts", are nonfunctional copies of mitochondrial genes that have been translocated to the nuclear genome. Numts have been used to study differences in mutation rates between the nuclear and mitochondrial genomes, but have also been implicated as troublesome for phylogenetic studies and DNA-based species identification (i.e., DNA barcoding). In this study, a suspected numt discovered during a study of mitochondrial cytochrome c oxidase I (COI) diversity in North American birds was targeted and sequenced from tyrant flycatchers (family: Tyrannidae). In total, the numt was found in five taxa representing two genera. Substitution rates were compared between COI and numt sequences. None of the numt sequences harboured stop codons nor frameshift mutations, but phylogenetic analysis revealed they had accumulated more amino acid substitutions than the mitochondrial COI sequences. Mitochondrial COI appeared to be preferentially amplified in most cases, but methods for numt detection are discussed for cases like this where sequences lack obvious features for identification. Because of its persistence across a broad taxonomic lineage, this numt could form a valuable model system for studying evolution in numts. The full size of the numt and its location within the nuclear genome are yet to be determined.     

家马核基因组中线粒体核内插入序列分析

在各种真核生物核基因组中,存在一些由线粒体基因组转移进入核基因组中的DNA片段,这些被认为是分子化石的片段叫做线粒体核内插入序列(Numt)。由于Numt与真实的线粒体序列高度相似,因此它的存在必然会成为PCR扩增线粒体DNA的不利因素。利用已经公布的家马(Equus caballus)基因组序列(2007年9月公布,GenBank登录号为NC_009144-NC_009175)对家马Numt进行了深入分析,共发现200个可能的Numt,长度范围为29到3727bp,其中有10个的长度大于800bp。分析结果显示由于不存在线粒体控制区域的疑似Numt,因此对基于此区域的群体遗传学研究不会产生影响。本研究还发现在家马进化过程中,第1号和27号染色体更倾向于接受线粒体序列的转移。以上结果将为今后马科动物的研究提供重要的参考信息,有助于避免在线粒体DNA研究中由于Numt污染的存在而得出错误的实验结果。     

Mitochondrial DNA-like sequences in the nuclear genome of the opossum genus Didelphis (Marsupialia: Didelphidae).

This study reports the occurrence of mtDNA-like sequences in the nuclear genome of the opossum genus Didelphis (Didelphidae, Marsupialia). A specific primer pair designed to amplify a region encompassing a 3' terminal 118 bp region of the cytochrome b gene, the Thr and Pro tRNA genes, and a 489 bp region of the D-loop of the D. virginiana mtDNA, was used in highly stringent PCR reactions. These PCR reactions resulted in several fragments per individual varying in size from 259 bp to 1 kb. The sequencing of some of these fragments showed the occurrence of paralogous mtDNA-like sequences among the PCR amplified fragments. Analyses of qualitative aspects of these sequences, their transition/transversion ratios, and phylogenetic relationships were conclusive in showing the occurrence of mtDNA-like sequences in the nuclear genome of the genus Didelphis. Comparisons and phylogenetic analysis of orthologous mtDNA from the four Didelphis species and paralogous nuclear sequences suggested that mtDNA migration to the nuclear genome occurred more than once in Didelphis evolution.     

Estimates of linkage disequilibrium and effective population size in rainbow trout

Background

Mitochondrial DNA (mtDNA) is widely used in population genetic and phylogenetic studies in animals. However, such studies can generate misleading results if the species concerned contain nuclear copies of mtDNA (Numts) as these may amplify in addition to, or even instead of, the authentic target mtDNA. The aim of this study was to determine if Numts are present in Aedes aegypti mosquitoes, to characterise any Numts detected, and to assess the utility of using mtDNA for population genetics studies in this species.

Results

BLAST searches revealed large numbers of Numts in the Ae. aegypti nuclear genome on 146 supercontigs. Although the majority are short (80% < 300 bp), some Numts are almost full length mtDNA copies. These long Numts are not due to misassembly of the nuclear genome sequence as the Numt-nuclear genome junctions could be recovered by amplification and sequencing. Numt evolution appears to be a complex process in Ae. aegypti with ongoing genomic integration, fragmentation and mutation and the secondary movement of Numts within the nuclear genome. The PCR amplification of the putative mtDNA nicotinamide adenine dinucleotide dehydrogenase subunit 4 (ND4) gene from 166 Southeast Asian Ae. aegypti mosquitoes generated a network with two highly divergent lineages (clade 1 and clade 2). Approximately 15% of the ND4 sequences were a composite of those from each clade indicating Numt amplification in addition to, or instead of, mtDNA. Clade 1 was shown to be composed at least partially of Numts by the removal of clade 1-specific bases from composite sequences following enrichment of the mtDNA. It is possible that all the clade 1 sequences in the network were Numts since the clade 2 sequences correspond to the known mitochondrial genome sequence and since all the individuals that produced clade 1 sequences were also found to contain clade 2 mtDNA-like sequences using clade 2-specific primers. However, either or both sets of clade sequences could have Numts since the BLAST searches revealed two long Numts that match clade 2 and one long Numt that matches clade 1. The substantial numbers of mutations in cloned ND4 PCR products also suggest there are both recently-derived clade 1 and clade 2 Numt sequences.

Conclusion

We conclude that Numts are prevalent in Ae. aegypti and that it is difficult to distinguish mtDNA sequences due to the presence of recently formed Numts. Given this, future population genetic or phylogenetic studies in Ae. aegypti should use nuclear, rather than mtDNA, markers.     

单个拟穴青蟹核内多拷贝线粒体COI假基因序列的研究(简报)

线粒体DNA(mtDNA)因它的高拷贝数、易扩增、高突变率、中性、低重组和母系遗传等特征,已广泛应用在系统进化和群体遗传研究方面[1,2].但是,由于mtDNA本身的特征,如异质性.     

Fugu genome does not contain mitochondrial pseudogenes

Contrary to previous observations that fish genomes are devoid of nuclear mitochondrial pseudogenes, a genome-wide survey identified a large number of "recent" and "ancient" nuclear mitochondrial DNA fragments (Numts) in the whole-genome sequences of the fugu (Takifugu rubripes), Tetraodon nigroviridis, and zebrafish (Danio rerio). We have analyzed the latest assembly (v4.0) of the fugu genome and show that, like the Anopheles genome, the fugu nuclear genome does not contain mitochondrial pseudogenes. Fugu assembly v4.0 contains a single scaffold representing the near complete sequence of the fugu mitochondria. The "recent" Numts identified by the previous study in fugu assembly v2.0 are in fact shotgun sequences of mitochondrial DNA that were misassembled with the nuclear sequences, whereas the "ancient" Numts appear to be the result of spurious matches. It is likely that the Numts identified in the genomes of Tetraodon and zebrafish are also similar artifacts. Shotgun sequences of whole genomes often include some mitochondrial sequences. Therefore, any Numts identified in shotgun-sequence assemblies should be verified by Southern hybridization or PCR amplification.     

Phylogeny of Anophelinae (Diptera: Culicidae) based on nuclear ribosomal and mitochondrial DNA sequences

Abstract Phylogenetic relationships among thirty-two species of mosquitoes in subfamily Anophelinae are inferred from portions of the mitochondrial genes COI and COII, the nuclear 18S small subunit rRNA gene and the expansion D2 region of the nuclear large subunit 28S rRNA gene. Sequences were obtained from the genera Anopheles , Bironella and Chagasia . Representatives of all six subgenera of Anopheles were included: Anopheles , Cellia , Kerteszia , Lophopodomyia , Nyssorhynchus and Stethomyia. Using parsimony and maximum likelihood methods, various combinations of these DNA sequence data were analysed separately: 18S, 28S, combined 18S and 28S, combined COI and COII, and combined 18S, 28S, COI and COII ('total evidence'). The combined rDNA data contain strong phylogenetic signal, moderately to strongly supporting most clades in MP and ML analyses; however, the mtDNA data (analysed as either nucleotide or amino acid sequences) contain little phylogenetic signal, except for relationships of very recently derived groups of species and, at the deepest level, for the monophyly of Anophelinae. The paraphyly of Anopheles relative to Bironella is confirmed by most analyses and statistical tests. Support for the monophyly of subgenera Anopheles , Cellia , Kerteszia and Nyssorhynchus is indicated by most analyses. Subgenus Lophopodomyia is reconstructed as the sister to Bironella , nested within a clade also containing Nyssorhynchus and Kerteszia . The most basal relationships within genus Anopheles are not well resolved by any of the data partitions, although the results of statistical analyses of the rDNA data (S-H-tests, likelihood ratio tests for monophyly and Bayesian MCMC analyses) suggest that the clade consisting of Bironella , Lophopodomyia , Nyssorhynchus and Kerteszia is the sister to the clade containing Cellia and Anopheles .     

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.     

Phylogenetic relationships of the New World Troidini swallowtails (Lepidoptera: Papilionidae) based on COI, COII, and EF-1alpha genes

A phylogeny of the Neotropical members of the Tribe Troidini (Lepidoptera: Papilionidae) was obtained with sequences of three protein-coding genes: two mitochondrial (COI and COII), and one nuclear (EF-1alpha). Parsimony and Bayesian analyses of 33 taxa resulted in very similar trees regardless of method used with the 27 troidines always forming a monophyletic clade. Within Troidini, the genus Battus is sister group to the remaining troidines, followed by a clade formed by the Paleotropical taxa (here represented by three exemplars). The genus Euryades is the next branch, and sister group of Parides. The genus Parides is monophyletic, and is divided into four main groups by Maximum Parsimony analysis, with the most basal group composed of tailed species restricted to SE Brazil. Character optimization of ecological and morphological traits over the phylogeny proposed for troidines indicated that the use of several species of Aristolochia is ancestral over the use of few or a single host-plant. For the other three characters, the ancestral states were the absence of long tails, forest as the primary habitat and oviposition solitary or in loose group of several eggs.     

云南边境地区果实蝇属种类DNA条形码鉴定

【目的】果实蝇属Bactrocera中有国际上重要的检疫性害虫, 基于形态的物种鉴定有一定的局限性。另一方面, 云南边境地区为东南亚地区实蝇入侵我国的重要通道。因此, 对该地区实蝇分子鉴定方法的研究对于该属物种的快速准确鉴定具有重要意义。本研究旨在探讨DNA条形码技术在果实蝇属物种鉴定中的有效性。【方法】使用线粒体基因COI和COII序列的通用引物对果实蝇属20个物种60份样品进行PCR扩增、测序和序列分析; 采取距离方法和建树方法评价2种序列的鉴别能力。【结果】COI和COII序列平均长度分别为682 bp和339 bp, 种内和种间遗传差异较大, 有较明显的遗传距离间隔（barcoding gap）, 鉴定成功率分别为91.2%和90.7%。另外, 分子系统树表明华实蝇亚属Sinodacus不是单系群。【结论】COI和COII序列均能够将绝大多数果实蝇属物种进行准确鉴别, 应用COI或COII序列进行果实蝇属物种鉴定具有一定的可行性。