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

Mitochondrial DNA intra-individual variation in a bumblebee species: A challenge for evolutionary studies and molecular identification

Mitochondrial DNA (mtDNA) regions have been widely used as molecular markers in evolutionary studies and species identification. However, the presence of heteroplasmy and NUMTs may represent obstacles. Heteroplasmy is a state where an organism has different mitochondrial haplotypes. NUMTs are nuclear pseudogenes originating from mtDNA sequences transferred to nuclear DNA. Evidences of heteroplasmy were already verified in the bumblebee Bombus morio in an earlier study. The present work investigated in more detail the presence of intra-individual haplotypes variation in this species. Heteroplasmy was detected in individuals from all the ten sampled locations, with an average of six heteroplasmic haplotypes per individual. In addition, some of these heteroplasmic haplotypes were shared among individuals from different locations, suggesting the existence of stable heteroplasmy in B. morio. These results demonstrated that heteroplasmy is likely to affect inferences based on mtDNA analysis, especially in phylogenetic, phylogeographic and population genetics studies. In addition, NUMTs were also detected. These sequences showed divergence of 2.7% to 12% in relation to the mitochondrial haplotypes. These levels of divergence could mislead conclusions in evolutionary studies and affect species identification through DNA barcoding.

Chimeric mitochondrial peptides from contiguous regular and swinger RNA

Previous mass spectrometry analyses described human mitochondrial peptides entirely translated from swinger RNAs, RNAs where polymerization systematically exchanged nucleotides. Exchanges follow one among 23 bijective transformation rules, nine symmetric exchanges (X ↔ Y, e.g. A ↔ C) and fourteen asymmetric exchanges (X → Y → Z → X, e.g. A → C → G → A), multiplying by 24 DNA's protein coding potential. Abrupt switches from regular to swinger polymerization produce chimeric RNAs. Here, human mitochondrial proteomic analyses assuming abrupt switches between regular and swinger transcriptions, detect chimeric peptides, encoded by part regular, part swinger RNA. Contiguous regular- and swinger-encoded residues within single peptides are stronger evidence for translation of swinger RNA than previously detected, entirely swinger-encoded peptides: regular parts are positive controls matched with contiguous swinger parts, increasing confidence in results. Chimeric peptides are 200 × rarer than swinger peptides (3/100,000 versus 6/1000). Among 186 peptides with > 8 residues for each regular and swinger parts, regular parts of eleven chimeric peptides correspond to six among the thirteen recognized, mitochondrial protein-coding genes. Chimeric peptides matching partly regular proteins are rarer and less expressed than chimeric peptides matching non-coding sequences, suggesting targeted degradation of misfolded proteins. Present results strengthen hypotheses that the short mitogenome encodes far more proteins than hitherto assumed. Entirely swinger-encoded proteins could exist.

•

Chimeric peptides are translated from contiguous regular and swinger RNA

•

They are 200x rarer than mitochondrial swinger peptides

•

Chimeric peptides integrated in regular mitochondrial proteins are downregulated

•

Contiguous regular parts are matched positive controls for swinger parts

•

The last point validates results beyond other statistical tests for robustness

Odintifier--A computational method for identifying insertions of organellar origin from modern and ancient high-throughput sequencing data based on haplotype phasing

Background

Cellular organelles with genomes of their own (e.g. plastids and mitochondria) can pass genetic sequences to other organellar genomes within the cell in many species across the eukaryote phylogeny. The extent of the occurrence of these organellar-derived inserted sequences (odins) is still unknown, but if not accounted for in genomic and phylogenetic studies, they can be a source of error. However, if correctly identified, these inserted sequences can be used for evolutionary and comparative genomic studies. Although such insertions can be detected using various laboratory and bioinformatic strategies, there is currently no straightforward way to apply them as a standard organellar genome assembly on next-generation sequencing data. Furthermore, most current methods for identification of such insertions are unsuitable for use on non-model organisms or ancient DNA datasets.

Results

We present a bioinformatic method that uses phasing algorithms to reconstruct both source and inserted organelle sequences. The method was tested in different shotgun and organellar-enriched DNA high-throughput sequencing (HTS) datasets from ancient and modern samples. Specifically, we used datasets from lions (Panthera leo ssp. and Panthera leo leo) to characterize insertions from mitochondrial origin, and from common grapevine (Vitis vinifera) and bugle (Ajuga reptans) to characterize insertions derived from plastid genomes. Comparison of the results against other available organelle genome assembly methods demonstrated that our new method provides an improvement in the sequence assembly.

Conclusion

Using datasets from a wide range of species and different levels of complexity we showed that our novel bioinformatic method based on phasing algorithms can be used to achieve the next two goals: i) reference-guided assembly of chloroplast/mitochondrial genomes from HTS data and ii) identification and simultaneous assembly of odins. This method represents the first application of haplotype phasing for automatic detection of odins and reference-based organellar genome assembly.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0682-1) contains supplementary material, which is available to authorized users.

Pleistocene aridification cycles shaped the contemporary genetic architecture of Southern African baboons

Plio-Pleistocene environmental change influenced the evolutionary history of many animal lineages in Africa, highlighting key roles for both climate and tectonics in the evolution of Africa’s faunal diversity. Here, we explore diversification in the southern African chacma baboon Papio ursinus sensu lato and reveal a dominant role for increasingly arid landscapes during past glacial cycles in shaping contemporary genetic structure. Recent work on baboons (Papio spp.) supports complex lineage structuring with a dominant pulse of diversification occurring 1-2Ma, and yet the link to palaeoenvironmental change remains largely untested. Phylogeographic reconstruction based on mitochondrial DNA sequence data supports a scenario where chacma baboon populations were likely restricted to refugia during periods of regional cooling and drying through the Late Pleistocene. The two lineages of chacma baboon, ursinus and griseipes, are strongly geographically structured, and demographic reconstruction together with spatial analysis of genetic variation point to possible climate-driven isolating events where baboons may have retreated to more optimum conditions during cooler, drier periods. Our analysis highlights a period of continuous population growth beginning in the Middle to Late Pleistocene in both the ursinus and the PG2 griseipes lineages. All three clades identified in the study then enter a state of declining population size (Nef) through to the Holocene; this is particularly marked in the last 20,000 years, most likely coincident with the Last Glacial Maximum. The pattern recovered here conforms to expectations based on the dynamic regional climate trends in southern Africa through the Pleistocene and provides further support for complex patterns of diversification in the region’s biodiversity.

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.

Complete chloroplast genome sequence of holoparasite Cistanche deserticola (Orobanchaceae) reveals gene loss and horizontal gene transfer from its host Haloxylon ammodendron (Chenopodiaceae)

BACKGROUND

The central function of chloroplasts is to carry out photosynthesis, and its gene content and structure are highly conserved across land plants. Parasitic plants, which have reduced photosynthetic ability, suffer gene losses from the chloroplast (cp) genome accompanied by the relaxation of selective constraints. Compared with the rapid rise in the number of cp genome sequences of photosynthetic organisms, there are limited data sets from parasitic plants. PRINCIPAL FINDINGS/SIGNIFICANCE: Here we report the complete sequence of the cp genome of Cistanche deserticola, a holoparasitic desert species belonging to the family Orobanchaceae. The cp genome of C. deserticola is greatly reduced both in size (102,657 bp) and in gene content, indicating that all genes required for photosynthesis suffer from gene loss and pseudogenization, except for psbM. The striking difference from other holoparasitic plants is that it retains almost a full set of tRNA genes, and it has lower dN/dS for most genes than another close holoparasitic plant, E. virginiana, suggesting that Cistanche deserticola has undergone fewer losses, either due to a reduced level of holoparasitism, or to a recent switch to this life history. We also found that the rpoC2 gene was present in two copies within C. deserticola. Its own copy has much shortened and turned out to be a pseudogene. Another copy, which was not located in its cp genome, was a homolog of the host plant, Haloxylon ammodendron (Chenopodiaceae), suggesting that it was acquired from its host via a horizontal gene transfer.

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.

Understanding phylogenetic incongruence: lessons from phyllostomid bats

All characters and trait systems in an organism share a common evolutionary history that can be estimated using phylogenetic methods. However, differential rates of change and the evolutionary mechanisms driving those rates result in pervasive phylogenetic conflict. These drivers need to be uncovered because mismatches between evolutionary processes and phylogenetic models can lead to high confidence in incorrect hypotheses. Incongruence between phylogenies derived from morphological versus molecular analyses, and between trees based on different subsets of molecular sequences has become pervasive as datasets have expanded rapidly in both characters and species. For more than a decade, evolutionary relationships among members of the New World bat family Phyllostomidae inferred from morphological and molecular data have been in conflict. Here, we develop and apply methods to minimize systematic biases, uncover the biological mechanisms underlying phylogenetic conflict, and outline data requirements for future phylogenomic and morphological data collection. We introduce new morphological data for phyllostomids and outgroups and expand previous molecular analyses to eliminate methodological sources of phylogenetic conflict such as taxonomic sampling, sparse character sampling, or use of different algorithms to estimate the phylogeny. We also evaluate the impact of biological sources of conflict: saturation in morphological changes and molecular substitutions, and other processes that result in incongruent trees, including convergent morphological and molecular evolution. Methodological sources of incongruence play some role in generating phylogenetic conflict, and are relatively easy to eliminate by matching taxa, collecting more characters, and applying the same algorithms to optimize phylogeny. The evolutionary patterns uncovered are consistent with multiple biological sources of conflict, including saturation in morphological and molecular changes, adaptive morphological convergence among nectar-feeding lineages, and incongruent gene trees. Applying methods to account for nucleotide sequence saturation reduces, but does not completely eliminate, phylogenetic conflict. We ruled out paralogy, lateral gene transfer, and poor taxon sampling and outgroup choices among the processes leading to incongruent gene trees in phyllostomid bats. Uncovering and countering the possible effects of introgression and lineage sorting of ancestral polymorphism on gene trees will require great leaps in genomic and allelic sequencing in this species-rich mammalian family. We also found evidence for adaptive molecular evolution leading to convergence in mitochondrial proteins among nectar-feeding lineages. In conclusion, the biological processes that generate phylogenetic conflict are ubiquitous, and overcoming incongruence requires better models and more data than have been collected even in well-studied organisms such as phyllostomid bats.

Primate numts and reticulate evolution of capped and golden leaf monkeys (Primates: Colobinae)

A recent phylogenetic study of langurs and leaf monkeys of South Asia suggested a reticulate evolution of capped and golden leaf monkeys through ancient hybridization between Semnopithecus and Trachypithecus .To test this hybridization scenario, I analysed nuclear copies of the mitochondrial cytochrome b gene (numts) from capped,golden and Phayre's leaf monkeys. These numts were aligned with mitochondrial cytochrome b sequences of various species belonging to the genera Semnopithecus and Trachypithecus .In the phylogenetic tree derived from this alignment,the numts fell into three distinct clades (A,B and C) suggesting three independent integration events.Clade A was basal to Semnopithecus, and clades B and C were basal to Trachypithecus. Among the numts in clades A and C were sequences derived from species not represented in their respective sister mitochondrial groups.This unusual placement of certain numts is taken as additional support for the hybridization scenario. Based on the molecular dating of these integration events,hybridization is estimated to have occurred around 7.1 to 3.4 million years ago.Capped and golden leaf monkeys might have to be assigned to a new genus to reconcile their unique evolutionary history. Additionally, northeast India appears to be a 'hot spot' for lineages that might have evolved through reticulate evolution.

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.

Molecular phylogeny and biogeography of langurs and leaf monkeys of South Asia (Primates: Colobinae)

The two recently proposed taxonomies of the langurs and leaf monkeys (Subfamily Colobinae) provide different implications to our understanding of the evolution of Nilgiri and purple-faced langurs. Groves (2001) [Groves, C.P., 2001. Primate Taxonomy. Smithsonian Institute Press, Washington], placed Nilgiri and purple-faced langurs in the genus Trachypithecus, thereby suggesting disjunct distribution of the genus Trachypithecus. [Brandon-Jones, D., Eudey, A.A., Geissmann, T., Groves, C.P., Melnick, D.J., Morales, J.C., Shekelle, M., Stewart, C.-B., 2003. Asian primate classification. Int. J. Primatol. 25, 97-162] placed these langurs in the genus Semnopithecus, which suggests convergence of morphological characters in Nilgiri and purple-faced langurs with Trachypithecus. To test these scenarios, we sequenced and analyzed the mitochondrial cytochrome b gene and two nuclear DNA-encoded genes, lysozyme and protamine P1, from a variety of colobine species. All three markers support the clustering of Nilgiri and purple-faced langurs with Hanuman langur (Semnopithecus), while leaf monkeys of Southeast Asian (Trachypithecus) form a distinct clade. The phylogenetic position of capped and golden leaf monkeys is still unresolved. It is likely that this species group might have evolved due to past hybridization between Semnopithecus and Trachypithecus clades.

Detection of mitochondrial insertions in the nucleus (NuMts) of Pleistocene and modern muskoxen

Background

Nuclear insertions of mitochondrial sequences (NuMts) have been identified in a wide variety of organisms. Trafficking of genetic material from the mitochondria to the nucleus has occurred frequently during mammalian evolution and can lead to the production of a large pool of sequences with varying degrees of homology to organellar mitochondrial DNA (mtDNA) sequences. This presents both opportunities and challenges for forensics, population genetics, evolutionary genetics, conservation biology and the study of DNA from ancient samples. Here we present a case in which difficulties in ascertaining the organellar mtDNA sequence from modern samples hindered their comparison to ancient DNA sequences.

Results

We obtained mitochondrial hypervariable region (HVR) sequences from six ancient samples of tundra muskox (Ovibos moschatus) that were reproducible but distinct from modern muskox sequences reported previously. Using the same PCR primers applied to the ancient specimens and the primers used to generate the modern muskox DNA sequences in a previous study, we failed to definitively identify the organellar sequence from the two modern muskox samples tested. Instead of anticipated sequence homogeneity, we obtained multiple unique sequences from both hair and blood of one modern specimen. Sequencing individual clones of a >1 kb PCR fragment from modern samples did not alleviate the problem as there was not a consistent match across the entire length of the sequences to Ovibos when compared to sequences in GenBank.

Conclusion

In specific taxa, due to nuclear insertions some regions of the mitochondrial genome may not be useful for the characterization of modern or ancient DNA.

Extensive variation in nuclear mitochondrial DNA content between the genomes of Phytophthora sojae and Phytophthora ramorum

Fragments of mitochondrial DNA (mtDNA) transferred to the nuclear genome are called nuclear mitochondrial DNAs (NUMTs). We report here a comparison of NUMT content between genomes from two species of the same genus. Analysis of the genomes of Phytophthora sojae and P. ramorum revealed large differences in the NUMT content of the two genomes: 16.27 x 10(-3) and 2.28 x 10(-3)% of each genome, respectively. Substantial differences also exist between the two species in the sizes of the NUMTs found in each genome, with ranges of 20 to 405 bp for P. sojae and 19 to 137 bp for P. ramorum. Furthermore, in P. sojae, fragments from the mitochondrial genes rns, rnl, coxl, and nad (various subunits) are found most frequently, whereas P. ramorum NUMTs most often originate from the cox3, rpsl4, nad4, and nad5 genes. The large differences in the presumptive mtDNA insertions suggest that the insertions occurred subsequent to the divergence of the two species, and this is supported by sequence comparisons among the NUMTs and the mtDNA sequences of the two species. P. sojae mtDNA sequences inserted in the nuclear genome appear to have been altered as a result of insertions, deletions, inversions, and translocations and provide insights into active mechanisms of sequence divergence in this plant pathogen. No clear examples were found of NUMTs forming functional nuclear genes or of NUMTs inserted into exons or introns of any nuclear gene.

Multiple origins of the Juan Fernández kelpfish fauna and evidence for frequent and unidirectional dispersal of cirrhitoid fishes across the South Pacific

Phylogenetic relationships were reconstructed among chironemid fishes based on morphological and molecular (lrRNA, NADH4, S7 ribosomal protein) characters. Two sympatric species from Juan Fernández in the southeast Pacific are not sister taxa, but rather exhibit independent relationships to Australian/New Zealand chironemids. The most plausible explanation for these relationships and contemporary distributions is an Australian/New Zealand origin of the family, followed by two trans-Pacific dispersal and colonization events, facilitated by larval entrapment within the West Wind Drift. This study demonstrates that the diversity of taxa on an island can reflect multiple colonizations, rather than in situ diversification, even in the case of very small, isolated, and geologically recent islands. When taken in conjunction with studies of related taxa, our results indicate that transoceanic dispersal of temperate cirrhitoid fishes in the South Pacific has been frequent and unidirectional. Molecular estimates of divergence time between southeast Pacific chironemids and their western relatives predate the emergence of Juan Fernández, consistent with hypotheses that much of the marine nearshore faunas of young southeast Pacific islands may be the product of successive transfer from older, now submerged islands.

Horizontal gene transfer of a plastid gene in the non-photosynthetic flowering plants Orobanche and Phelipanche (Orobanchaceae)

Plastid sequences are among the most widely used in phylogenetic and phylogeographic studies in flowering plants, where they are usually assumed to evolve like non-recombining, uniparentally transmitted, single-copy genes. Among others, this assumption can be violated by intracellular gene transfer (IGT) within cells or by the exchange of genes across mating barriers (horizontal gene transfer, HGT). We report on HGT of a plastid region including rps2, trnL-F, and rbcL in a group of non-photosynthetic flowering plants. Species of the parasitic broomrape genus Phelipanche harbor two copies of rps2, a plastid ribosomal gene, one corresponding to the phylogenetic position of the respective species, the other being horizontally acquired from the related broomrape genus Orobanche. While the vertically transmitted copies probably reside within the plastid genome, the localization of the horizontally acquired copies is not known. With both donor and recipient being parasitic plants, a possible pathway for the exchange of genetic material is via a commonly attacked host.

MITOCHONDRIAL DNA AND HUMAN EVOLUTION

Several unique properties of human mitochondrial DNA (mtDNA), including its high copy number, maternal inheritance, lack of recombination, and high mutation rate, have made it the molecule of choice for studies of human population history and evolution. Here we review the current state of knowledge concerning these properties, how mtDNA variation is studied, what we have learned, and what the future likely holds. We conclude that increasingly, mtDNA studies are (and should be) supplemented with analyses of the Y-chromosome and other nuclear DNA variation. Some serious issues need to be addressed concerning nuclear inserts, database quality, and the possible influence of selection on mtDNA variation. Nonetheless, mtDNA studies will continue to play an important role in such areas as examining socio-cultural influences on human genetic variation, ancient DNA, certain forensic DNA applications, and in tracing personal genetic history.

Intraordinal phylogenetics of treeshrews (Mammalia: Scandentia) based on evidence from the mitochondrial 12S rRNA gene

Despite their traditional and continuing prominence in studies of interordinal mammalian phylogenetics, treeshrews (order Scandentia) remain relatively unstudied with respect to their intraordinal relationships. At the same time, significant morphological variation among living treeshrews has been shown to have direct relevance to higher-level interpretations of character state change as reconstructed in traditional interordinal studies, which have often included only a single species of treeshrew. Therefore, the importance of resolving relationships among treeshrews extends well beyond a better understanding of patterns of diversification within the order. A recent review highlighted several shortcomings in published studies of treeshrew phylogenetics based on morphology. Here we present the first investigation of treeshrew phylogenetics based on DNA sequences, utilizing previously published sequences from the mitochondrial 12S rRNA gene and combining them with newly generated sequence data from 15 species. Parsimony, likelihood, and Bayesian analyses all strongly support a sister relationship between Ptilocercus and the remaining species, further substantiating its recent elevation to familial status. Dendrogale is consistently recovered as the next taxon to diverge, but relationships among the remaining taxa are poorly supported by these data. We provide evidence for a relatively rapid radiation within the genera Tupaia and Urogale, but limited resolution precludes more than a cursory interpretation of biogeographic patterns.

Unreliable mtDNA data due to nuclear insertions: a cautionary tale from analysis of humans and other great apes

Analysis of mitochondrial DNA sequence variation has been used extensively to study the evolutionary relationships of individuals and populations, both within and across species. So ubiquitous and easily acquired are mtDNA data that it has been suggested that such data could serve as a taxonomic 'barcode' for an objective species classification scheme. However, there are technical pitfalls associated with the acquisition of mtDNA data. One problem is the presence of translocated pieces of mtDNA in the nuclear genome of many taxa that may be mistaken for authentic organellar mtDNA. We assessed the extent to which such 'numt' sequences may pose an overlooked problem in analyses of mtDNA from humans and apes. Using long-range polymerase chain reaction (PCR), we generated necessarily authentic mtDNA sequences for comparison with sequences obtained using typical methods for a segment of the mtDNA control region in humans, chimpanzees, bonobos, gorillas and orangutans. Results revealed that gorillas are notable for having such a variety of numt sequences bearing high similarity to authentic mtDNA that any analysis of mtDNA using standard approaches is rendered impossible. Studies on humans, chimpanzees, bonobos or orangutans are apparently less problematic. One implication is that explicit measures need to be taken to authenticate mtDNA sequences in newly studied taxa or when any irregularities arise. Furthermore, some taxa may not be amenable to analysis of mtDNA variation at all.

Molecular systematics of armadillos (Xenarthra, Dasypodidae): contribution of maximum likelihood and Bayesian analyses of mitochondrial and nuclear genes

The 30 living species of armadillos, anteaters, and sloths (Mammalia: Xenarthra) represent one of the three major clades of placentals. Armadillos (Cingulata: Dasypodidae) are the earliest and most speciose xenarthran lineage with 21 described species. The question of their tricky phylogeny was here studied by adding two mitochondrial genes (NADH dehydrogenase subunit 1 [ND1] and 12S ribosomal RNA [12S rRNA]) to the three protein-coding nuclear genes (alpha2B adrenergic receptor [ADRA2B], breast cancer susceptibility exon 11 [BRCA1], and von Willebrand factor exon 28 [VWF]) yielding a total of 6869 aligned nucleotide sites for thirteen xenarthran species. The two mitochondrial genes were characterized by marked excesses of transitions over transversions-with a strong bias toward CT transitions for the 12S rRNA-and exhibited two- to fivefold faster evolutionary rates than the fastest nuclear gene (ADRA2B). Maximum likelihood and Bayesian phylogenetic analyses supported the monophyly of Dasypodinae, Tolypeutinae, and Euphractinae, with the latter two armadillo subfamilies strongly clustering together. Conflicting branching points between individual genes involved relationships within the subfamilies Tolypeutinae and Euphractinae. Owing to a greater number of informative sites, the overall concatenation favored the mitochondrial topology with the classical grouping of Cabassous and Priodontes within Tolypeutinae, and a close relationship between Euphractus and Chaetophractus within Euphractinae. However, low statistical support values associated with almost equal distributions of apomorphies among alternatives suggested that two parallel events of rapid speciation occurred within these two armadillo subfamilies.