Disturbed population genetics: suspected introgressive hybridization between two Mnais damselfly species (Odonata)

Mnais costalis and M. pruinosa are damselflies (Odonata: Calopterygidae) with low dispersal abilities, both during their aquatic stream-living immature stage and their flying adult stage. A previous nuclear DNA (nDNA) sequencing and morphology study showed that these two species are very closely related, and cohabit widely in western Japan. The two species, however, segregate microhabitats along a stream: M. costalis lives in the lower reaches, and M. pruinosa in the upper reaches. In this study, our analyses were based on mitochondrial DNA (mtDNA), which usually mutates faster and is more variable among individuals than nDNA, and which is inherited maternally. We found that most COI haplotypes were shared between the two species, and that for most study sites interspecific riverine genetic structures were not clarified by mtDNA analysis. Incongruent population genetic structures based on nDNA and mtDNA suggested hybridization and introgression of mtDNA between the two species.     

Phylogenetic relationships and biogeography of xantusiid lizards, inferred from mitochondrial DNA sequences.

Portions of two mitochondrial genes (12S ribosomal RNA and cytochrome b) were sequenced in seven species to examine phylogenetic relationships within the lizard family Xantusiidae. Phylogenies derived from these sequences (709 total bp) are concordant and indicate that the Cuban species Cricosaura typica is the sister group to all other xantusiids. The Middle American genus Lepidophyma is the closest relative of Xantusia, and X. riversiana (California Islands) the closest relative of X. vigilis (mainland). These findings are not in agreement either with the results of a recent morphological analysis that united Cricosaura and Lepidophyma as closest relatives or with past studies that have recognized X. riversiana as a separate genus. Levels of sequence divergence, as well as the age and affinities of some mainland fossil taxa, suggest that the origin of Cricosaura was associated with the tectonic evolution of the Greater Antilles in the late Cretaceous. These results further demonstrate that significant resolution of phylogenies can be obtained with relatively short DNA sequences and that these mitochondrial genes are concordant in their estimation of phylogeny.     

Divergent character clines across a recent secondary contact zone in a Hispaniolan lizard

Studies of genetic contact zones provide valuable information regarding the processes of population divergence, adaptation and speciation. In this paper, I examine transitions in morphology, mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) haplotypes across a recent secondary contact zone in a Hispaniolan lizard Ameiva chrysolaema . Maximum likelihood cline fitting analyses suggest non-coincidence of cline centers and that the mtDNA cline is significantly displaced to the west of the remaining clines. nDNA and morphological clines are coincident and tend to be associated with the prevailing environmental gradient. The lack of cytonuclear disequilibrium near the center of the contact zone and the non-coincidence of character clines suggest that this zone does not conform to a tension zone model of hybridization; thus, gene flow across the zone does not seem to be impeded. The extremely narrow width of the dorsal scale size cline and the close association of this cline with the steepness of the environmental (precipitation) gradient suggest that this character may be under environmental selection. Taken together, this contact zone appears to be structured by a combination of mtDNA introgression, possibly associated with eastward movement of the zone, and environmental selection on some characters.     

Hidden support from unpromising data sets strongly unites snakes with anguimorph ‘lizards’

A combined analysis of nuclear, mitochondrial and morphological data robustly resolves snakes as the sister taxon to anguimorph ‘lizards’. Analysed in isolation, nuclear DNA (nDNA) produces a trichotomy between snakes, iguanians and anguimorphs, mitochondrial DNA (mtDNA) is largely uninformative at deeper levels, and morphology tends to nest snakes deep within anguimorphs or with various legless squamate groups. When analysed simultaneously, the nuclear signal is sufficiently strong that mtDNA and morphology are constrained to choose between alternative resolutions of the iguanian–anguimorph–snake trichotomy (generated by the nDNA) – and both support the snake–anguimorph solution. Combined analyses of fast‐evolving or idiosyncratically evolving markers (mtNDA, morphology) with conservative traits (e.g. nuclear genes) might be the best way to resolve ancient, closely spaced divergences. Fast or idiosyncratic markers potentially provide the most information about short, ancient internodes, but can converge on spurious trees if analysed in isolation. However, if constrained to only choosing between plausible trees, such data can contribute unique and valuable phylogenetic signal that resolves such problematic divergences.     

The role of mitochondrial introgression in illuminating the evolutionary history of Nearctic treefrogs

Inferring the evolutionary history of lineages often becomes difficult when gene histories are in conflict with each other. Introgression, for example, can cause DNA sequences from one species to be more similar to sequences of a different species and lead to incongruence amongst gene trees. However, incorporating congruent and incongruent locus‐specific phylogenetic estimates with the geographical distribution of lineages may provide valuable insight into evolutionary processes important to speciation. In this study, we investigated mitochondrial introgression within the Hyla eximia group to better understand its role in illuminating the evolutionary history and phylogeography of these treefrogs. We reconstructed and compared the matrilineal history of the Hyla eximia group with estimates of evolutionary history inferred from nuclear genes. We tested for introgression within the mitochondrial and nuclear genes using a posterior predictive checking approach. Reconstructions of the species tree based on the mitochondrial DNA (mtDNA) and nuclear DNA data were strongly discordant. Introgression between lineages was widespread in the mtDNA data set (145 occurrences amongst 11 of the 16 lineages), but uncommon in the nuclear genes (12 occurrences amongst four of the 16 lineages). Nonetheless, the geographical structuring of mtDNA within species provides valuable information on biogeographical areas, ancient areas of hybridization, and unique histories of lineages within the H. eximia group. These results suggest that the combination of nuclear, mitochondrial, and spatial information can provide a more complete picture of ‘how evolutionary history played out’, particularly in cases where mitochondrial introgression is known to occur. © 2014 The Linnean Society of London     

Introgression and selection shaped the evolutionary history of sympatric sister‐species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister‐species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site‐associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister‐species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174–0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.     

The highly cross-fertile coral species,Acropora hyacinthus and Acropora cytherea,constitute statistically distinguishable lineages

A major challenge for understanding the evolutionary genetics of mass-spawning corals is to explain the maintenance of discrete morphospecies in view of high rates of interspecific fertilization in vitro and nonmonophyletic patterns in molecular phylogenies. In this study, we focused on Acropora cytherea and A. hyacinthus, which have one of the highest potentials for interspecific fertilization. Using sequences of a nuclear intron, we performed phylogenetic and nested clade analyses (NCA). Both species were polyphyletic in molecular phylogenies, but the NCA indicated that they constitute statistically distinguishable lineages. Phylogenetic analysis using an intergenic region of the mitochondrial DNA (mtDNA), was inconclusive because of low levels of variability in this marker. The position of these two species differed between the nuclear DNA (nDNA) and mtDNA phylogenies and was also at odds with a cladistic analysis based on morphology. We conclude that despite the potential for high levels of hybridization and introgression, A. cytherea and A. hyacinthus constitute statistically distinguishable lineages and their taxonomic status is consistent with the cohesion species concept.     

Phylogeography and lineage‐specific patterns of genetic diversity and molecular evolution in a group of North American skinks

Geography influences the evolutionary trajectory of species by mediating opportunities for hybridization, gene flow, demographic shifts and adaptation. We sought to understand how geography and introgression can generate species‐specific patterns of genetic diversity by examining phylogeographical relationships in the North American skink species Plestiodon multivirgatus and P. tetragrammus (Squamata: Scincidae). Using a multilocus dataset (three mitochondrial genes, four nuclear genes; a total of 3455 bp) we discovered mito‐nuclear discordance, consistent with mtDNA introgression. We further tested for evidence of species‐wide mtDNA introgression by using comparisons of genetic diversity, selection tests and extended Bayesian skyline analyses. Our findings suggest that P. multivirgatus acquired its mitochondrial genome from P. tetragrammus after their initial divergence. This putative species‐wide mitochondrial capture was further evidenced by statistically indistinguishable substitution rates between mtDNA and nDNA in P. multivirgatus. This rate discrepancy was observed in P. multivirgatus but not P. tetragrammus, which has important implications for studies that combine mtDNA and nDNA sequences when inferring time since divergence between taxa. Our findings suggest that by facilitating opportunities for interspecific introgression, geography can alter the course of molecular evolution between recently diverged lineages.     

Isolation and introgression in the Intermountain West: contrasting gene genealogies reveal the complex biogeographic history of the American pika (Ochotona princeps)

Aim We studied the history of colonization, diversification and introgression among major phylogroups in the American pika, Ochotona princeps (Lagomorpha), using comparative and statistical phylogeographic methods. Our goal was to understand how Pleistocene climatic fluctuations have shaped the distribution of diversity at mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) loci in this alpine specialist. Location North America’s Intermountain West. Methods We accumulated mtDNA sequence data (c. 560–1700 bp) from 232 pikas representing 64 localities, and sequenced two nuclear introns (mast cell growth factor, c. 550 bp, n = 148; protein kinase C iota, c. 660 bp, n = 139) from a subset of individuals. To determine the distribution of major mtDNA lineages, we conducted a phylogenetic analysis on the mtDNA sequence data, and we calculated divergence times among the lineages using a Bayesian Markov chain Monte Carlo approach. Relationships among nuclear alleles were explored with minimum spanning networks. Finally, we conducted coalescent simulations of alternative models of population history to test for congruence between nDNA and mtDNA responses to Pleistocene glacial cycles. Results We found that: (1) all individuals could be assigned to one of five allopatric mtDNA lineages; (2) lineages are associated with separate mountain provinces; (3) lineages originated from at least two rounds of differentiation; (4) nDNA and mtDNA markers exhibited overall phylogeographic congruence; and (5) introgression among phylogroups has occurred at nuclear loci since their initial isolation. Main conclusions Pika populations associated with different mountain systems have followed separate but not completely independent evolutionary trajectories through multiple glacial cycles. Range expansion associated with climate cooling (i.e. glaciations) promoted genetic admixture among populations within mountain ranges. It also permitted periodic contact and introgression between phylogroups associated with different mountain systems, the record of which is retained at nDNA but not mtDNA loci. Evidence for different histories at nuclear and mtDNA loci (i.e. periodic introgression versus deep isolation, respectively) emphasizes the importance of multilocus perspectives for reconstructing complete population histories.     

Loss of CpG dinucleotides from DNA. VI. Methylation of mitochondrial and chloroplast genes

From nucleotide sequences of mitochondrial and chloroplast genes the probable frequency of the CpG----TpG + CpA substitutions was determined. These substitutions may indicate the level of prior DNA methylation. It was found that the level of this methylation is significantly lower in mitochondrial DNA (mtDNA) and chloroplast DNA (chDNA) than in nuclear DNA (nDNA) of the same species. The species (taxon) specificity of mtDNA and chDNA methylation was revealed. A correlation was found between the level of CpG methylation in nDNA, and mtDNA and chDNA in different organisms. It is shown that cytosine residues in CpG were not subjected to significant methylation in the fungi and invertebrate mtDNA and also in the algae chDNA. In contrast, the vertebrate mtDNA bears the impress of CpG-supression, which is confirmed by direct data on methylation of these DNA. Here the first data on the possible enzymatic methylation of the plant mtDNA and chDNA were obtained. It was shown that the degree of CpG-suppression in the 5S rRNA nuclear genes of lower and higher plants is significantly higher in the chloroplast genes of 4,5S and 5S rRNA. From data on pea chDNA hydrolysis with MspI and HpaII it was established that in CCGG sequences this DNA is not methylated. The role of DNA methylation in increasing the mutation rate and in accelerating the evolutionary rates of vertebrate mtDNA is discussed.     

Sorting through the chaff, nDNA gene trees for phylogenetic inference and hybrid identification of annual sunflowers (Helianthus sect. Helianthus)

The annual sunflowers (Helianthus sect. Helianthus) present a formidable challenge for phylogenetic inference because of ancient hybrid speciation, recent introgression, and suspected issues with deep coalescence. Here we analyze sequence data from 11 nuclear DNA (nDNA) genes for multiple genotypes of species within the section to (1) reconstruct the phylogeny of this group, (2) explore the utility of nDNA gene trees for detecting hybrid speciation and introgression; and (3) test an empirical method of hybrid identification based on the phylogenetic congruence of nDNA gene trees from tightly linked genes. We uncovered considerable topological heterogeneity among gene trees with or without three previously identified hybrid species included in the analyses, as well as a general lack of reciprocal monophyly of species. Nonetheless, partitioned Bayesian analyses provided strong support for the reciprocal monophyly of all species except H. annuus (0.89 PP), the most widespread and abundant annual sunflower. Previous hypotheses of relationships among taxa were generally strongly supported (1.0 PP), except among taxa typically associated with H. annuus, apparently due to the paraphyly of the latter in all gene trees. While the individual nDNA gene trees provided a useful means for detecting recent hybridization, identification of ancient hybridization was problematic for all ancient hybrid species, even when linkage was considered. We discuss biological factors that affect the efficacy of phylogenetic methods for hybrid identification.     

Multiple nuclear gene sequences identify phylogenetic species boundaries in the rapidly radiating clade of Mexican ambystomatid salamanders

Delimiting the boundaries of species involved in radiations is critical to understanding the tempo and mode of lineage formation. Single locus gene trees may or may not reflect the underlying pattern of population divergence and lineage formation, yet they constitute the vast majority of the empirical data in species radiations. In this study we make use of an expressed sequence tag (EST) database to perform nuclear (nDNA) and mitochondrial (mtDNA) genealogical tests of species boundaries in Ambystoma ordinarium, a member of an adaptive radiation of metamorphic and paedomorphic salamanders (the Ambystoma tigrinum complex) that have diversified across terrestrial and aquatic environments. Gene tree comparisons demonstrate extensive nonmonophyly in the mtDNA genealogy of A. ordinarium, while seven of eight independent nuclear loci resolve the species as monophyletic or nearly so, and diagnose it as a well-resolved genealogical species. A differential introgression hypothesis is supported by the observation that western A. ordinarium localities contain mtDNA haplotypes that are identical or minimally diverged from haplotypes sampled from a nearby paedomorphic species, Ambystoma dumerilii, while most nDNA trees place these species in distant phylogenetic positions. These results provide a strong example of how historical introgression can lead to radical differences between gene trees and species histories, even among currently allopatric species with divergent life history adaptations and morphologies. They also demonstrate how EST-based nuclear resources can be used to more fully resolve the phylogenetic history of species radiations.     

Convergent evolution of kin-based sociality in a lizard

Studies of social birds and mammals have produced extensive theory regarding the formation and dynamics of kin-based social groups in vertebrates. However, comparing kin dynamics in birds and mammals to social reptiles provides the opportunity to identify selective factors that promote independent origins of kin sociality across vertebrates. We combined a 5-year mark-recapture study with a DNA microsatellite analysis of relatedness in a social lizard (Xantusia vigilis) to examine the formation and stability of kin groups. We found that these lizards are highly sedentary and that groups often form through the delayed dispersal of offspring. Groups containing juveniles had higher relatedness than adult-only groups, as juveniles were commonly found in aggregations with at least one parent and/or sibling. Groups containing nuclear family members were more stable than groups of less-related lizards, as predicted by social theory. We conclude that X. vigilis aggregations conform to patterns of kin sociality observed in avian and mammalian systems and represent an example of convergent evolution in social systems. We suggest that kin-based sociality in this and other lizards may be a by-product of viviparity, which can promote delayed juvenile dispersal by allowing prolonged interaction between a neonate and its mother.     

Nuclear loci and coalescent methods support ancient hybridization as cause of mitochondrial paraphyly between gadwall and falcated duck (Anas spp.)

Many species have mitochondrial DNA lineages that are phylogenetically intermixed with other species, but studies have rarely tested the cause of such paraphyly. In this study, we tested two hypotheses that could explain mitochondrial paraphyly of Holarctic gadwalls (Anas strepera) with respect to Asian falcated ducks (A. falcata). First, hybridization could have resulted in falcated duck mitochondrial DNA (mtDNA) introgressing into the gadwall gene pool. Second, gadwalls and falcated ducks could have diverged so recently that mtDNA lineages have not sorted to reciprocal monophyly. We used coalescent analyses of three independent loci to distinguish between these two hypotheses. Two lines of evidence support introgression. First, analyses of the three loci combined show that some introgression is necessary to explain current genetic diversity in gadwalls. Second, we generated alternative predictions regarding time since divergence estimated from mtDNA: falcated ducks and gadwalls would have diverged between 65,000 and 700,000 years before present (ybp) under the introgression hypothesis and between 11,000 and 76,000 ybp under the incomplete lineage sorting hypothesis. The two independent nuclear introns indicated that these species diverged between 210,000 and 5,200,000 ybp, which did not overlap the predicted time for incomplete lineage sorting. These analyses also suggested that ancient introgression ( approximately 14,000 ybp) has resulted in the widespread distribution and high frequency of falcated-like mtDNA (5.5% of haplotypes) in North America. This is the first study to use a rigorous quantitative framework to reject incomplete lineage sorting as the cause of mitochondrial paraphyly.     

Structure of animal mitochondrial DNA: nucleotide composition, pyrimidine clusters, and methylation character.

The nucleotide composition, relative concentration of pyrimidine clusters, and the degree of methylation of the mitochondrial and nuclear DNA's of various vertebrates and the protozoan Crithidia oncopelti have been studied. With respect to the relative concentration of GC pairs, the mtDNA of animals (bull, rat) does not differ from the corresponding nDNA. The relative concentration of GC pairs in the mtDNA of certain fish and birds is 1.5-2.5 mole% higher than in the respective nDNA. The kinetoplast DNA of the protozoan C. oncopelti (where the relative concentration of the GC pairs is 42.9 mole %) differs very sharply in composition from the nDNA (where the relative concentration of GC pairs is 51.3 mole %). The mtDNA's and kDNA's studied are distinguished from the respective nDNA'S by a lower degree of clustering of pyrimidine nucleotides. The proportion of mono- and dipyrimidine fragments in the mtDNA and kDNA is 30 mole %, while in the nDNA it does not exceed 23 mole %. The relative concentration of long pyrimidine clusters (hexapyrimidine clusters of larger) in the mtDNA is smaller than in the nDNA by a factor of 2-5. The low degree of clustering of the pyrimidine nucleotides is apparently characteristic of all the known mtDNA's and may support the fact that they have a single type of organization and are of a single origin. All the vertebrate mtDNA's studied contain 5-methylcytosine as a minor base (1.5-3.15 mole %), and their level of methylation is 1.5-2 times greater than that in the respective nDNA's. It has been shown that animals display species specificity with respect to the 5-methylcytosine content in the mtDNA. Its distribution among the pyrimidine clusters in the bovine heart mtDNA differs substantially from that in the nDNA. This suggests that the methylation specificities of nuclear and mitochondrial DNA are different. A DNA methylase, which effects the in vitro methylation of cytosine residues both in the homologous mtDNA and in different heterologous DNA's, has been found in rat liver and bovine heart mitochondria. The specificity of the in vitro methylation of the cytosine residues in the same heterologous Escherichia coli B DNA by the nuclear and mitochondrial enzymes is different: The mitochondrial enzyme methylates predominantly in monopyrimidine fragments, and the nuclear enzyme methylates mostly in di- and tripyrimidine fragments. They, therefore, recognize different nucleotide sequences.     

Oxidative damage to mitochondrial DNA is inversely related to maximum life span in the heart and brain of mammals.

DNA damage is considered of paramount importance in aging. Among causes of this damage, free radical attack, particularly from mitochondrial origin, is receiving special attention. If oxidative damage to DNA is involved in aging, long-lived animals (which age slowly) should show lower levels of markers of this kind of damage than short-lived ones. However, this possibility has not heretofore been investigated. In this study, steady-state levels of 8-oxo-7, 8-dihydro-2'-deoxyguanosine (8-oxodG) referred to deoxyguanosine (dG) were measured by high performance liquid chromatography (HPLC) in the mitochondrial (mtDNA) and nuclear (nDNA) DNA from the heart of eight and the brain of six mammalian species ranging in maximum life span (MLSP) from 3.5 to 46 years. Exactly the same digestion of DNA to deoxynucleosides and HPLC protocols was used for mtDNA and nDNA. Significantly higher (three- to ninefold) 8-oxodG/dG values were found in mtDNA than in nDNA in all the species studied in both tissues. 8-oxodG/dG in nDNA did not correlate with MLSP across species either in the heart (r=-0.68; P<0.06) or brain (r = 0.53; P<0.27). However, 8-oxodG/dG in mtDNA was inversely correlated with MLSP both in heart (r=-0.92; P<0.001) and brain (r=-0.88; P<0.016) tissues following the power function y = a(.)x(b), where y is 8-oxodG/dG and x is the MLSP. This agrees with the consistent observation that mitochondrial free radical generation is also lower in long-lived than in short-lived species. The results obtained agree with the notion that oxygen radicals of mitochondrial origin oxidatively damage mtDNA in a way related to the aging rate of each species.-Barja, G., Herrero, A. Oxidative damage to mitochondrial DNA is inversely related to maximum life span in the heart and brain of mammals.     

Phylogeny and temporal diversification of darters (Percidae: Etheostomatinae)

Discussions aimed at resolution of the Tree of Life are most often focused on the interrelationships of major organismal lineages. In this study, we focus on the resolution of some of the most apical branches in the Tree of Life through exploration of the phylogenetic relationships of darters, a species-rich clade of North American freshwater fishes. With a near-complete taxon sampling of close to 250 species, we aim to investigate strategies for efficient multilocus data sampling and the estimation of divergence times using relaxed-clock methods when a clade lacks a fossil record. Our phylogenetic data set comprises a single mitochondrial DNA (mtDNA) gene and two nuclear genes sampled from 245 of the 248 darter species. This dense sampling allows us to determine if a modest amount of nuclear DNA sequence data can resolve relationships among closely related animal species. Darters lack a fossil record to provide age calibration priors in relaxed-clock analyses. Therefore, we use a near-complete species-sampled phylogeny of the perciform clade Centrarchidae, which has a rich fossil record, to assess two distinct strategies of external calibration in relaxed-clock divergence time estimates of darters: using ages inferred from the fossil record and molecular evolutionary rate estimates. Comparison of Bayesian phylogenies inferred from mtDNA and nuclear genes reveals that heterospecific mtDNA is present in approximately 12.5% of all darter species. We identify three patterns of mtDNA introgression in darters: proximal mtDNA transfer, which involves the transfer of mtDNA among extant and sympatric darter species, indeterminate introgression, which involves the transfer of mtDNA from a lineage that cannot be confidently identified because the introgressed haplotypes are not clearly referable to mtDNA haplotypes in any recognized species, and deep introgression, which is characterized by species diversification within a recipient clade subsequent to the transfer of heterospecific mtDNA. The results of our analyses indicate that DNA sequences sampled from single-copy nuclear genes can provide appreciable phylogenetic resolution for closely related animal species. A well-resolved near-complete species-sampled phylogeny of darters was estimated with Bayesian methods using a concatenated mtDNA and nuclear gene data set with all identified heterospecific mtDNA haplotypes treated as missing data. The relaxed-clock analyses resulted in very similar posterior age estimates across the three sampled genes and methods of calibration and therefore offer a viable strategy for estimating divergence times for clades that lack a fossil record. In addition, an informative rank-free clade-based classification of darters that preserves the rich history of nomenclature in the group and provides formal taxonomic communication of darter clades was constructed using the mtDNA and nuclear gene phylogeny. On the whole, the appeal of mtDNA for phylogeny inference among closely related animal species is diminished by the observations of extensive mtDNA introgression and by finding appreciable phylogenetic signal in a modest sampling of nuclear genes in our phylogenetic analyses of darters.     

Nuclear mitochondrial pseudogenes

As has been demonstrated recently, the transfer of genetic material from mitochondria to the nucleus and its integration into the nuclear genome is a continuous and dynamic process. Fragments of mitochondrial DNA (mtDNA) are incorporated in the nuclear genome as noncoding sequences, which are called nuclear mitochondrial pseudogenes (NUMT pseudogenes or NUMT inserts). In various eukaryotes, NUMT pseudogenes are distributed through different chromosomes to form a “library” of mtDNA fragments, providing important information on genome evolution. The escape of mtDNA from mitochondria is mostly associated with mitochondrial damage and mitophagy. Fragments of mtDNA may be integrated into nuclear DNA (nDNA) during repair of double-strand breaks (DSBs), which are caused by endogenous or exogenous agents. DSB repair of nDNA with a capture of mtDNA fragments may occur via nonhomologous end joining or a similar mechanism that involves microhomologous terminal sequences. An analysis of the available data makes it possible to suppose that the NUMT pseudogene formation rate depends on the DSB rate in nDNA, the activity of the repair systems, and the number of mtDNA fragments leaving organelles and migrating into the nucleus. Such situations are likely after exposure to damaging agents, first and foremost, ionizing radiation. Not only do new NUMT pseudogenes change the genome structure in the regions of their integration, but they may also have a significant impact on the actualization of genetic information. The de novo integration of NUMT pseudogenes in the nuclear genome may play a role in various pathologies and aging. NUMT pseudogenes may cause errors in PCR-based analyses of free mtDNA as a component of total cell DNA because of their coamplification.     

Recurrent introgression of mitochondrial DNA among hares (Lepus spp.) revealed by species-tree inference and coalescent simulations

Understanding recent speciation history requires merging phylogenetic and population genetics approaches, taking into account the persistence of ancestral polymorphism and possible introgression. The emergence of a clear phylogeny of hares (genus Lepus) has been hampered by poor genomic sampling and possible occurrence of mitochondrial DNA (mtDNA) introgression from the arctic/boreal Lepus timidus into several European temperate and possibly American boreal species. However, no formal test of introgression, taking also incomplete lineage sorting into account, has been done. Here, to clarify the yet poorly resolved species phylogeny of hares and test hypotheses of mtDNA introgression, we sequenced 14 nuclear DNA and 2 mtDNA fragments (8205 and 1113 bp, respectively) in 50 specimens from 11 hare species from Eurasia, North America, and Africa. By applying an isolation-with-migration model to the nuclear data on subsets of species, we find evidence for very limited gene flow from L. timidus into most temperate European species, and not into the American boreal ones. Using a multilocus coalescent-based method, we infer the species phylogeny, which we find highly incongruent with mtDNA phylogeny using parametric bootstrap. Simulations of mtDNA evolution under the speciation history inferred from nuclear genes did not support the hypothesis of mtDNA introgression from L. timidus into the American L. townsendii but did suggest introgression from L. timidus into 4 temperate European species. One such event likely resulted in the complete replacement of the aboriginal mtDNA of L. castroviejoi and of its sister species L. corsicanus. It is remarkable that mtDNA introgression in hares is frequent, extensive, and always from the same donor arctic species. We discuss possible explanations for the phenomenon in relation to the dynamics of range expansions and species replacements during the climatic oscillations of the Pleistocene.