Multiple nuclear pseudogenes of mitochondrial cytochrome b in Ctenomys (Caviomorpha, rodentia) with either great similarity to or high divergence from the true mitochondrial sequence

Numerous insertions of mitochondrial DNA in the genome of the northern mole vole, Ellobius talpinus

BACKGROUND

Ellobius talpinus is a subterranean rodent representing an attractive model in population ecology studies due to its highly special lifestyle and sociality. In such studies, mitochondrial DNA (mtDNA) is widely used. However, if nuclear copies of mtDNA, aka NUMTs, are present, they may co-amplify with the target mtDNA fragment, generating misleading results. The aim of this study was to determine whether NUMTs are present in E. talpinus.

METHODS AND RESULTS

PCR amplification of the putative mtDNA CytB-D-loop fragment using 'universal' primers from 56 E. talpinus samples produced multiple double peaks in 90% of the sequencing chromatograms. To reveal NUMTs, molecular cloning and sequencing of PCR products of three specimens was conducted, followed by phylogenetic analysis. The pseudogene nature of three out of the seven detected haplotypes was confirmed by their basal positions in relation to other Ellobius haplotypes in the phylogenetic tree. Additionally, 'haplotype B' was basal in relation to other E. talpinus haplotypes and found present in very distant sampling sites. BLASTN search revealed 195 NUMTs in the E. talpinus nuclear genome, including fragments of all four PCR amplified pseudogenes. Although the majority of the NUMTs studied were short, the entire mtDNA had copies in the nuclear genome. The most numerous NUMTs were found for rrnL, COXI, and D-loop.

CONCLUSIONS

Numerous NUMTs are present in E. talpinus and can be difficult to discriminate against mtDNA sequences. Thus, in future population or phylogenetic studies in E. talpinus, the possibility of cryptic NUMTs amplification should always be taken into account.

Whole mitochondrial genome of long-clawed mole vole (Prometheomys schaposchnikowi) from Turkey, with its phylogenetic relationships

The mitogenome of Prometheomys schaposchnikowi was characterized for the first time as a circular DNA molecule (16.284 bp), containing 37 coding and 2 non-coding regions. In the mitogenome, ND6 and 8 tRNA genes were encoded on the light chain, while 12 PCGs, 14 tRNAs, 2 rRNAs, D-loop and O_L were encoded on the heavy chain. The most common initiation codon in PCGs was ATG. As in many mammals, incomplete stop codons in P. schaposchnikowi were in the COX3, ND1 and ND4. Phylogenetic relationships were revealed using Bayesian method and the 13 PCGs. Seven genera (Arvicola, Dicrostonyx, Lasiopodomys, Myodes, Ondatra, Proedromys and Prometheomys) formed a monophyletic group, while Eothenomys, Microtus and Neodon were paraphyletic. P. schaposchnikowi constituted the most basal group within Arvicolinae. Divergence time estimation suggested that P. schaposchnikowi diversified during the Miocene (16.28 Mya). Further molecular studies are needed to test the distinctiveness and diversity of the genus Prometheomys.

The Anatolian glacial refugium and human-mediated colonization: a phylogeographical study of the stone marten (Martes foina) in Turkey

The Anatolian Peninsula, comprising most of modern Turkey, is situated at the junction of Europe, the Middle East and Asia. Together with its complex geomorphological and climatic history, this has given rise to a rich fauna and flora, which exhibits a wide range of historical biogeographical patterns. The stone marten (Martes foina) is a small carnivore that is found across the temperate Palaearctic region including Anatolia, where it is often associated with habitats modified by humans, but few genetic data exist for this species. We sequenced a 1840-bp region of the mitochondrial genome from 97 martens sampled across the peninsula and intron 7 of the nuclear β-fibrinogen gene from 53 of these. Two mitochondrial lineages were recovered, with overlapping eastern and western distributions, but there was no geographical structure for the autosomal marker. Coalescent analyses indicated that both of the lineages originated during the Last Glacial Maximum, one of them within an eastern Anatolian refugium and the other in a western Anatolian or Balkan refugium. The western lineage colonized most of Europe in the Holocene, while the eastern lineage may be endemic to Anatolia, from where it colonized the Iberian Peninsula via human translocation. The presence of at least one refugial stone marten population highlights the importance of Anatolia to the preservation of genetic variation and biodiversity.

DNA barcoding of the rodent genus Oligoryzomys (Cricetidae: Sigmodontinae): mitogenomic-anchored database and identification of nuclear mitochondrial translocations (Numts)

DNA barcoding has become a standard method for species identification in taxonomically complex groups. An important step of the barcoding process is the construction of a library of voucher-based material that was properly identified by independent methods, free of inaccurate identification, and paralogs. We provide here a cytochrome oxidase I (mt-Co1) DNA barcode database for species of the genus Oligoryzomys, based on type material and karyotyped specimens, and anchored on the mitochondrial genome of one species of Oligoryzomys, O. stramineus. To evaluate the taxonomic determination of new COI sequences, we assessed species intra/interspecific genetic distances (barcode gap), performed the General Mixed Yule Coalescent method (GMYC) for lineages' delimitation, and identified diagnostic nucleotides for each species of Oligoryzomys. Phylogenetic analyses of Oligoryzomys were performed on 2 datasets including 14 of the 23 recognized species of this genus: a mt-Co1 only matrix, and a concatenated matrix including mt-Co1, cytochrome b (mt-Cytb), and intron 7 of the nuclear fibrinogen beta chain gene (i7Fgb). We recovered nuclear-mitochondrial translocated (Numts) pseudogenes on our samples and identified several published sequences that are cases of Numts. We analyzed the rate of non-synonymous and synonymous substitution, which were higher in Numts in comparison to mtDNA sequences. GMYC delimitations and DNA barcode gap results highlight the need for further work that integrate molecular, karyotypic, and morphological analyses, as well as additional sampling, to tackle persistent problems in the taxonomy of Oligoryzomys.

Mitogenomic phylogenetics of the bank vole Clethrionomys glareolus, a model system for studying end-glacial colonization of Europe

We have revisited the mtDNA phylogeny of the bank vole Clethrionomys glareolus based on Sanger and next-generation Illumina sequencing of 32 complete mitochondrial genomes. The bank vole is a key study species for understanding the response of European fauna to the climate change following the Last Glacial Maximum (LGM) and one of the most convincing examples of a woodland mammal surviving in cryptic northern glacial refugia in Europe. The genomes sequenced included multiple representatives of each of the eight bank vole clades previously described based on cytochrome b (cob) sequences. All clades with the exception of the Basque - likely a misidentified pseudogene clade - were highly supported in all phylogenetic analyses and the relationships between the clades were resolved with high confidence. Our data extend the distribution of the Carpathian clade, the marker of a northern glacial refugium in the Carpathian Mountains, to include Britain and Fennoscandia (but not adjacent areas of continental Europe). The Carpathian sub-clade that colonized Britain and Fennoscandia had a somewhat different history from the sub-clade currently found in or close to the Carpathians and may have derived from a more north-westerly refugial area. The two bank vole populations that colonized Britain at the end of the last glaciation are for the first time linked with particular continental clades, the first colonists with the Carpathian clade and the second colonists with the western clade originating in a more southerly refugium in the vicinity of the Alps. We however found no evidence that a functional divergence of proteins encoded in the mitochondrial genome promoted the partial genetic replacement of the first colonists by the second colonists detected previously in southern Britain. We did identify one codon site that changed more often and more radically in the tree than expected and where the observed amino acid change may affect the reductase activity of the cytochrome bc1 complex, but the change was not specific to a particular clade. We also found an excess of radical changes to the primary protein structure for geographically restricted clades from southern Italy and Norway, respectively, possibly related to stronger selective pressure at the latitudinal extremes of the bank vole distribution. However, overall, we find little evidence of pervasive effects of deviation from neutrality on bank vole mtDNA phylogeography.

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.

Genetic identification of Iberian rodent species using both mitochondrial and nuclear loci: application to noninvasive sampling

Species identification through noninvasive sampling is increasingly used in animal conservation genetics, given that it obviates the need to handle free-living individuals. Noninvasive sampling is particularly valuable for elusive and small species such as rodents. Although rodents are not usually assumed to be the most obvious target for conservation, of the 21 species or near-species present in Iberia, three are considered endangered and declining, while several others are poorly studied. Here, we develop a genetic tool for identifying all rodent species in Iberia by noninvasive genetic sampling. To achieve this purpose, we selected one mitochondrial gene [cytochrome b (cyt-b)] and one nuclear gene [interphotoreceptor retinoid-binding protein (IRBP)], which we first sequenced using tissue samples. Both genes allow for the phylogenetic distinction of all species except the sibling species Microtus lusitanicus and Microtus duodecimcostatus. Overall, cyt-b showed higher resolution than IRBP, revealing a clear barcoding gap. To allow these markers to be applied to noninvasive samples, we selected a short highly diagnostic fragment from each gene, which we used to obtain sequences from faeces and bones from owl pellets. Amplification success for the cyt-b and IRBP fragment was 85% and 43% in faecal and 88% and 64% in owl-pellet DNA extractions, respectively. The method allows the unambiguous identification of the great majority of Iberian rodent species from noninvasive samples, with application in studies of distribution, spatial ecology and population dynamics, and for conservation.

Frequency matrix approach demonstrates high sequence quality in avian BARCODEs and highlights cryptic pseudogenes

The accuracy of DNA barcode databases is critical for research and practical applications. Here we apply a frequency matrix to assess sequencing errors in a very large set of avian BARCODEs. Using 11,000 sequences from 2,700 bird species, we show most avian cytochrome c oxidase I (COI) nucleotide and amino acid sequences vary within a narrow range. Except for third codon positions, nearly all (96%) sites were highly conserved or limited to two nucleotides or two amino acids. A large number of positions had very low frequency variants present in single individuals of a species; these were strongly concentrated at the ends of the barcode segment, consistent with sequencing error. In addition, a small fraction (0.1%) of BARCODEs had multiple very low frequency variants shared among individuals of a species; these were found to represent overlooked cryptic pseudogenes lacking stop codons. The calculated upper limit of sequencing error was 8 × 10(-5) errors/nucleotide, which was relatively high for direct Sanger sequencing of amplified DNA, but unlikely to compromise species identification. Our results confirm the high quality of the avian BARCODE database and demonstrate significant quality improvement in avian COI records deposited in GenBank over the past decade. This approach has potential application for genetic database quality control, discovery of cryptic pseudogenes, and studies of low-level genetic variation.

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.

Nuclear copies of mitochondrial genes: another problem for ancient DNA

The application of ancient DNA techniques is subject to many problems caused primarily by low quality and by low quantity of DNA. For these reasons most studies employing ancient DNA rely on the characterization of mitochondrial DNA, which is present in many more copies per cell than nuclear DNA and hence more copies are likely to survive. We used universal and taxon specific mitochondrial primers to amplify DNA from museum specimens, and found many instances where the amplification of nuclear copies of the mitochondrial gene (numts) instead of the targeted mitochondrial fragment had occurred. Furthermore, the likelihood of amplifying numts increased dramatically when universal primers were utilized. Here we suggest that ancient DNA practitioners must consider the possibility that numts can be amplified at higher rates than previously thought. This is another complication for ancient DNA studies, but it also suggests that more extensive inclusion of nuclear markers in ancient DNA studies should be feasible.

Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified

Nuclear mitochondrial pseudogenes (numts) are nonfunctional copies of mtDNA in the nucleus that have been found in major clades of eukaryotic organisms. They can be easily coamplified with orthologous mtDNA by using conserved universal primers; however, this is especially problematic for DNA barcoding, which attempts to characterize all living organisms by using a short fragment of the mitochondrial cytochrome c oxidase I (COI) gene. Here, we study the effect of numts on DNA barcoding based on phylogenetic and barcoding analyses of numt and mtDNA sequences in two divergent lineages of arthropods: grasshoppers and crayfish. Single individuals from both organisms have numts of the COI gene, many of which are highly divergent from orthologous mtDNA sequences, and DNA barcoding analysis incorrectly overestimates the number of unique species based on the standard metric of 3% sequence divergence. Removal of numts based on a careful examination of sequence characteristics, including indels, in-frame stop codons, and nucleotide composition, drastically reduces the incorrect inferences of the number of unique species, but even such rigorous quality control measures fail to identify certain numts. We also show that the distribution of numts is lineage-specific and the presence of numts cannot be known a priori. Whereas DNA barcoding strives for rapid and inexpensive generation of molecular species tags, we demonstrate that the presence of COI numts makes this goal difficult to achieve when numts are prevalent and can introduce serious ambiguity into DNA barcoding.

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.

Co-amplification of mitochondrial pseudogenes inCalomys musculinus(Rodentia, Cricetidae): a source of error in phylogeographic studies

In a previous phylogeographic study of the rodent Calomys musculinus , 24 haplotypes of the mitochondrial DNA D-loop region were detected using the restriction fragment length polymorphism technique (PCR-RFLP). Seven percent of the individuals showed patterns in which the sum of the sizes of the restriction fragments exceeded the size of the original PCR product. In the present paper we analyze possible causes of these atypical haplotypes. PCR products were cloned, and two or three different clones from a single individual were detected by their RFLP patterns. Nine clones with different restriction patterns were selected for sequence analyses. A maximum parsimony phylogenetic analysis revealed two well-supported paraphyletic groups. One group comprised sequences showing low nucleotide divergence compared with the most common haplotypes detected in the phylogeographic study. The other group was basal to the three species of Calomys other than C. musculinus included in the study; the mutations in the short portion of the cytochrome b gene amplified corresponded to 12 amino acid substitutions. The results suggest that two independent insertions of mtDNA sequences into the nucleus occurred; these sequences would co-amplify in the PCR procedure. Identification of pseudogenes is crucial to obtain reliable reconstruction of the intraspecific genealogy in phylogeographic studies.

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

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

Multigenic and morphometric differentiation of ground squirrels (Spermophilus, Scuiridae, Rodentia) in Turkey, with a description of a new species

This study reports on the molecular phylogenetics of ground squirrels, genus Spermophilus, in Turkey using cytochrome b (1140bp), part of the D-loop and flanking tRNAs (572bp), X chromosome (867-1051bp) and Y chromosome (983-989bp) DNA sequences. Individuals also were characterized by karyotype and with geometric morphometric analyses of mandibles and skulls. Two hundred fourteen individuals from 91 localities were studied. All the data support the recognition of a new species in SW Anatolia: the Taurus ground squirrel Spermophilus taurensis sp. nov. The new species has a small distribution in the Taurus Mountains in an area that is a hotspot for biodiversity. Molecular clock analysis suggests that the new species diverged from the European ground squirrel, Spermophilus citellus, about 2.5 million years ago and that the ancestor of these two species diverged from the widespread Anatolian ground squirrel, Spermophilus xanthoprymnus, about 5 million years ago. Morphometric differentiation in skull and mandible shape among the three species is incomplete, but statistically significant. S. xanthoprymnus is subdivided into five cytochrome b phylogroups and we use these data to infer the location of glacial refugia where the species lived during the last glacial maximum. This study illustrates the potential of combined molecular and morphometric studies to uncover new Anatolian species and to reconstruct their phylogeographic history. The new species is important for squirrel taxonomy and for understanding Eurasian mammal evolution.

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.

An unusual phylogeography in the bushcricket Ephippiger ephippiger from Southern France

Pleistocene glaciations have played a major role in species divergence. The bushcricket Ephippiger ephippiger shows unusual patterns of intraspecific variation in multiple traits across Southern Europe. This is centred in Southern France, and evidence implies that it results from secondary contact after differentiation in Pleistocene refugia. However, the possible time scales involved, locations of the refugia and patterns of expansion remain obscure. This study sequenced the COII (507 BP) and cyt b (428 BP) mitochondrial genes to examine the intraspecific phylogeography of Western European samples of E. ephippiger. A minimum evolution tree revealed little resolution between described subspecies of E. ephippiger. Strikingly, populations from the Pyrenees and Mediterranean coastal region contained a complex genetic structure corresponding to major river valleys, independent of the traditional taxonomy. Samples of the subspecies E. e. vitium formed a distinct clade, perhaps supporting their taxonomic status. However, other forms (cruciger and cunii) were not genetically distinct, which is surprising given differences in their morphology and behaviour. The extent of the genetic divergence between Pyreneen valleys is unexpectedly deep, with average Tamura-Nei distances of around 14% (net distances of 11%) separating the main clades of coding COII sequences. Cyt b showed a similar pattern, but was confounded by some non-coding probable pseudogenes. If a conventional insect molecular clock is applied, these cryptic clades must pre-date the Pleistocene, and hypotheses for their history are discussed. However, mtDNA divergence in Ephippiger is not evolving in a clock-like manner, because a likelihood ratio test rejects clock assumptions for the COII sequences.

A simple method for domestic animal identification in Argentina using PCR-RFLP analysis of cytochrome b gene

We developed a simple, quick assay in order to discriminate forensic samples among human, and common domestic and livestock species of the Pampean region, Argentina. A mitochondrial cytochrome b fragment amplified with universal primers was separately digested with three restriction enzymes (AluI, HaeIII, and HinfI) and the resulting fragments were resolved through electrophoresis in polyacrylamide gels. This PCR-RFLP method allowed us to identify the target species and worked on degraded samples. The assay was successfully applied in livestock robbery cases in Argentine, and may be useful when attempting a first assessment as to the specific status of a forensic evidence.

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

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

Molecular phylogeny of the speciose vole genus Microtus (Arvicolinae, Rodentia) inferred from mitochondrial DNA sequences

Voles of the genus Microtus represent one of the most speciose mammalian genera in the Holarctic. We established a molecular phylogeny for Microtus to resolve contentious issues of systematic relationships and evolutionary history in this genus. A total of 81 specimens representing ten Microtus species endemic to Europe as well as eight Eurasian, six Asian and one Holarctic species were sequenced for the entire cytochrome b gene (1140 bp). A further 25 sequences were retrieved from GenBank, providing data on an additional 23, mainly Nearctic, Microtus species. Phylogenetic analysis of these 48 species generated four well-supported monophyletic lineages. The genus Chionomys, snow voles, formed a distinct and well-supported lineage separate from the genus Microtus. The subgenus Microtus formed the strongest supported lineage with two sublineages displaying a close relationship between the arvalis species group (common voles) and the socialis species group (social voles). Monophyly of the Palearctic pitymyid voles, subgenus Terricola, was supported, and this subgenus was also subdivided into two monophyletic species groups. Together, these groupings clarify long-standing taxonomic uncertainties in Microtus. In addition, the "Asian" and the Nearctic lineages reported previously were identified although the latter group was not supported. However, relationships among the main Microtus branches were not resolved, suggesting a rapid and potentially simultaneous radiation of a widespread ancestor early in the history of the genus. This and subsequent radiations discernible in the cytochrome b phylogeny, show the considerable potential of Microtus for analysis of historical and ecological determinants of speciation in small mammals. It is evident that speciation is an ongoing process in the genus and that the molecular data provides a vital insight into current species limits as well as cladogenic events of the past.

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.

A highly divergent mitochondrial DNA lineage of Microtus agrestis in southern Europe

The Mediterranean peninsulas constitute important areas for endemism and intraspecific variation, and are likely places for cryptic biodiversity. We assessed the phylogeographic pattern of field voles (Microtus agrestis) in southern and central Europe by sequence analysis of a 385-bp fragment of the mitochondrial cytochrome b gene in 74 specimens from 44 localities. The majority of samples consisted of skulls collected from owl pellets. The data revealed a highly distinct cytochrome b lineage in an area ranging from Portugal to Hungary. This southern field vole phylogroup differed by a sequence divergence of 5.6-7.1% from the remaining haplotypes, a level of divergence comparable to that found between known Microtus sibling species. However, this ancient phylogeographic break that dates back many glacial cycles has not been recognised previously by either morphology or karyotype. The southern cytochrome b lineage was further divided into two well-defined sublineages that appear to have derived from different glacial refugia in the Iberian Peninsula.

How did pygmy shrews colonize Ireland? Clues from a phylogenetic analysis of mitochondrial cytochrome b sequences

There is a long-standing debate as to how Ireland attained its present fauna; we help to inform this debate with a molecular study of one species. A 1110 base pair fragment of the mitochondrial cytochrome b gene was sequenced in 74 specimens of the pygmy shrew, Sorex minutus, collected from throughout its western Palaearctic range. Phylogenetic analysis of these sequences revealed several well-supported lineages. Most of the 65 haplotypes belonged to a northern lineage, which ranged from Britain in the west to Lake Baikal in the east. The other lineages were largely limited to Iberia, Italy and the Balkans. One exception, however, was a lineage found in both Ireland and Andorra. This affinity, and the large difference between the mitochondrial sequences of Irish and British individuals, suggest that pygmy shrews did not colonize Ireland via a land connection from Britain, as has been previously supposed, but instead were introduced by boat from southwest continental Europe. All the Irish pygmy shrews analysed were identical or very similar in cytochrome b sequence, suggesting an extreme founding event.

Out of Anatolia: longitudinal gradients in genetic diversity support an eastern origin for a circum-Mediterranean oak gallwasp Andricus quercustozae

Many studies have addressed the latitudinal gradients in intraspecific genetic diversity of European taxa generated during postglacial range expansion from southern refugia. Although Asia Minor is known to be a centre of diversity for many taxa, relatively few studies have considered its potential role as a Pleistocene refugium or a potential source for more ancient westward range expansion into Europe. Here we address these issues for an oak gallwasp, Andricus quercustozae (Hymenoptera: Cynipidae), whose distribution extends from Morocco along the northern coast of the Mediterranean through Turkey to Iran. We use sequence data for a fragment of the mitochondrial gene cytochrome b and allele frequency data for 12 polymorphic allozyme loci to answer the following questions: (1). which regions represent current centres of genetic diversity for A. quercustozae? Do eastern populations represent one refuge or several discrete glacial refugia? (2). Can we infer the timescale and sequence of the colonization processes linking current centres of diversity? Our results suggest that A. quercustozae was present in five distinct refugia (Iberia, Italy, the Balkans, southwestern Turkey and northeastern Turkey) with recent genetic exchange between Italy and Hungary. Genetic diversity is greatest in the Turkish refugia, suggesting that European populations are either (a). derived from Asia Minor, or (b). subject to more frequent population bottlenecks. Although Iberian populations show the lowest diversity for putatively selectively neutral markers, they have colonized a new oak host and represent a genetically and biologically discrete entity within the species.

Lifecycle closure, lineage sorting, and hybridization revealed in a phylogenetic analysis of European oak gallwasps (Hymenoptera: Cynipidae: Cynipini) using mitochondrial sequence data

Oak gallwasps are cyclically parthenogenetic insects that induce a wide diversity of highly complex species- and generation-specific galls on oaks and other Fagaceae. Phylogenetic relationships within oak gallwasps remain to be established, while sexual and parthenogenetic generations of many species remain unpaired. Previous work on oak gallwasps has revealed substantial intra-specific variation, particularly between regions known to represent discrete Pleistocene glacial refuges. Here we use statistical phylogenetic inference methods on sequence data for a fragment of the mitochondrial cytochrome b gene to reconstruct the relationships among 62 oak gallwasp species. For 16 of these we also include 23 additional cytochrome b haplotype sequences from different Pleistocene refuge areas to test the effect of intra-specific variation on inter-specific phylogeny reconstruction. The reconstructed phylogenies show good intra-generic resolution and identify several conserved clades, but fail to reconstruct either very recent or very ancient divergences. Nine of the 16 species represented by multiple haplotypes are not monophyletic. The apparent discordance between the recovered gene tree and the current taxonomic classification can be explained through: (a) collapsing of some species currently known only from either a sexual or a parthenogenetic generation into a single cyclically parthenogenetic entity; (b) sorting of ancestral polymorphism in diverging lineages, and (c) horizontal transfer of haplotypes, perhaps due to hybridization within glacial refuges. Our conclusions emphasise the need for careful intra-specific sampling when reconstructing phylogenies for radiations of closely related species and imply that for certain taxonomic groups full phylogenetic resolution (using molecular markers) may not be attainable.

Genetic diversity and population history of two related seabird species based on mitochondrial DNA control region sequences

Geographical variation in two related seabird species, the razorbill (Alca torda) and common guillemot (Uria aalge), was investigated using sequence analysis of mitochondrial DNA (mtDNA) control regions. We determined the nucleotide sequence of the variable 5' segment of the control region in razorbills and common guillemots from breeding colonies across the Atlantic Ocean. The ecology and life history characteristics of razorbill and common guillemot are in many respects similar. They are both considered highly philopatric and have largely overlapping distributions in temperate and subarctic regions of the North Atlantic, yet the species were found to differ widely in the extent and spatial distribution of mtDNA variation. Moreover, the differences in genetic differentiation and diversity were in the opposite direction to that expected from a consideration of traditional classifications and current population sizes. Indices of genetic diversity were highest in razorbill and varied among colonies, as did genotype frequencies, suggestive of restrictions to gene flow. The distribution of genetic variation suggests that razorbills originated from a refugial population in the south-western Atlantic Ocean through sequential founder events and subsequent expansion in the east and north. In common guillemots, genetic diversity was low and there was a lack of geographical structure, consistent with a recent population bottleneck, expansion and gene flow. We suggest that the reduced level of genetic diversity and differentiation in the common guillemot is caused by an inherent propensity for repeated population bottlenecks and concomitantly unstable population structure related to their specialized feeding ecology.

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.