Relative patterns and rates of evolution in heron nuclear and mitochondrial DNA

Mitochondrial cytochrome b sequence data from 15 species of herons (Aves: Ardeidae), representing 13 genera, were compared with DNA hybridization data of single-copy nuclear DNA (scnDNA) from the same species in a taxonomic congruence assessment of heron phylogeny. The two data sets produced a partially resolved, completely congruent estimate of phylogeny with the following basic structure: (Tigrisoma, Cochlearius, (((Zebrilus, (Ixobrychus, Botaurus)), (((Ardea, Casmerodius), Bubulcus), ((Egretta thula, Egretta caerulea, Egretta tricolor), Syrigma), Butorides, Nycticorax, Nyctanassa)))). Because congruence indicated similar phylogenetic information in the two data sets, we used the relatively unsaturated DNA hybridization distances as surrogates of time to examine graphically the patterns and rates of change in cytochrome b distances. Cytochrome b distances were computed either from whole sequences or from partitioned sequences consisting of transitions, transversions, specific codon site positions, or specific protein-coding regions. These graphical comparisons indicated that unpartitioned cytochrome b has evolved at 5-10 times the rate of scnDNA. Third-position transversions appeared to offer the most useful sequence partition for phylogenetic analysis because of their relatively fast rate of substitution (two times that of scnDNA) and negligible saturation. We also examined lineage-based rates of evolution by comparing branch length patterns between the nuclear and cytochrome b trees. The degree of correlation in corresponding branch lengths between cytochrome b and DNA hybridization trees depended on DNA sequence partitioning. When cytochrome b sequences were not partitioned, branch lengths in the cytochrome b and DNA hybridization trees were not correlated. However, when cytochrome b sequences were reduced to third-position transversions (i.e., unsaturated, relatively fast changing data), branch lengths were correlated. This finding suggests that lineage-based rates of DNA evolution in nuclear and mitochondrial genomes are influenced by common causes.     

Sequence evolution and phylogenetic signal in control-region and cytochrome b sequences of rainbow fishes (Melanotaeniidae)

The nucleotide sequences of segments of the cytochrome b gene (351 bp), the tRNA(Pro) gene (49 bp), and the control region (approximately 313 bp) of mitochondrial DNA were obtained from 26 fish representing different populations and species of Melanotaenia and one species of Glossolepis, freshwater rainbow fishes confined to Australia and New Guinea. The purpose was to investigate relative rates and patterns of sequence evolution. Overall levels of divergence were similar for the cytochrome b and tRNA control-region sequences, both ranging from < 1% within subspecies to 15%-19% between genera. However, the patterns of sequence evolution differed. For the cytochrome b gene, transitions consistently exceeded transversions, the bias ranging from 4.2:1 to 2:1, depending on the level of sequence divergence. However, in the control-region sequence, a bias toward transitions (2:1) was observed only in comparisons between very similar sequences, and transversions outnumbered transitions in comparisons of divergent sequences. Graphic comparisons suggested that the control region was saturated for transitions at relatively low levels of sequence divergence but accumulated transversions at a greater rate than did the cytochrome b sequence. These distinct patterns of base substitution are associated with differences in A+T content, which is 70% for the tRNA control- region segment versus 50% for cytochrome b. A test for skewness in the distribution of lengths of random trees indicated that both segments contained phylogenetic signal. Parsimony analyses of the data from the two regions, with or without weighting schemes appropriate to the respective patterns of sequence evolution, identified the same five groupings of sequences, but the relationships among the groups differed. However, in most cases the branches uniting different combinations of groups were poorly supported, and the differences among topologies were insignificant. Considering the observed patterns of base substitution and the results of the phylogenetic analyses, we deduce that both the control region and cytochrome b are appropriate for population genetic studies but that the control region is less effective than cytochrome b for resolving relationships among divergent lineages of rainbow fishes.      

Molecular phylogenetic relationships based on mitochondrial and nuclear gene sequences for the Todies (Todus, Todidae) of the Caribbean

We used mitochondrial/nuclear gene sequence analyses to determine the historical relationships of the endemic species of Todus (Aves: Todidae) from the Caribbean. We collected 1920-bp of nucleotide sequence data from the mitochondrial genes cytochrome b, ATPase 6, ATPase 8, and 591-bp of the single-copy nuclear gene c-mos for all Todus species and representatives of their outgroup taxa (Hylomanes, Barypthengus, Chloroceryle, Ceryle, and Galbula) to reconstruct the evolutionary history (via parsimony and maximum likelihood) of the five Todus species. The substitution rates among the mitochondrial genes were found to be much higher than the substitution rate for the c-mos gene, consequently resulting in higher substitutional saturation for the mitochondrial genes. When we applied weighting schemes to account for the variance in substitutional heterogeneity among the genes then parsimony and likelihood analyses both demonstrate that the genus Todus is monophyletic and closer to the Hylomanes and Barypthengus genera than the Chloroceryle and Ceryle genera. The mitochondrial-gene trees and nuclear-gene trees both show similar results, thus providing support for the relationships among the taxa from loci within two independently evolving genomes. The nuclear gene c-mos was found, therefore, to be a viable nuclear gene candidate for resolving intermediate and deep divergences.     

Mode and Tempo of Molecular Evolution in the Nematode Caenorhabditis: Cytochrome Oxidase II and Calmodulin Sequences

Through direct sequencing methods, the mitochondrial gene for cytochrome oxidase subunit two (CO II) and the single-copy nuclear gene for calmodulin were compared among strains of Caenorhabidits elegans and two other Caenorhabditis species (C. remanei and C. briggsae). In addition the CO II sequence was determined from a distantly related nematode, Steinernema intermedii. Among the 11 strains of C. elegans tested, there are four types of CO II gene, arising from two major lineages. Levels of intraspecific difference in the CO II gene are low (less than 2.0%) compared to the extraordinary divergence between congeneric species, which is about 50% when corrected for multiple hits. Concordant with the increase in divergence between taxa is a change in the pattern of substitution from a strong transition bias (24 transitions compared to two transversions) within species to a substitution pattern that appears to reflect the base composition of the mitochondrial genome when more divergent nematodes are compared. The base composition of the Caenorhabditis CO II gene is strongly biased toward A + T at all three positions of codons and appears to constrain the amino acid composition of the protein. Both the CO II and calmodulin genes show extreme conservation of amino acid sequences. When the accumulation of changes at silent sites in the two genes is compared among strains, it becomes evident that the mitochondrial gene is changing faster than the nuclear gene.     

Multiple nuclear insertions of mitochondrial cytochrome b sequences in callitrichine primates

We report the presence of four nuclear paralogs of a 380-bp segment of cytochrome b in callitrichine primates (marmosets and tamarins). The mitochondrial cytochrome b sequence and each nuclear paralog were obtained from several species, allowing multiple comparisons of rates and patterns of substitution both between mitochondrial and nuclear sequences and among nuclear sequences. The mitochondrial DNA had high overall rates of molecular evolution and a strong bias toward substitutions at third codon positions. Rates of molecular evolution among the nuclear sequences were low and constant, and there were small differences in substitution patterns among the nuclear clades which were probably attributable to the small number of sites involved. A novel method of phylogenetic reconstruction based on the large difference in rates of evolution at different codon positions among mitochondrial and nuclear clades was used to determine whether different nuclear paralogs represent independent transposition events or duplications following a single insertion. This method is generally applicable in cases where differences in pattern of molecular evolution are known, and it showed that at least three of the four nuclear clades represent independent transposition events. The insertion events giving rise to two of the nuclear clades predate the divergence of the callitrichines, whereas those leading to the other two nuclear clades may have occurred in the common ancestor of marmosets.     

Construction of novel cytochrome b genes in yeast mitochondria by subtraction or addition of introns

The mitochondrial cob-box gene coding for apocytochrome b in yeast has five introns and six exons or two introns and three exons depending on the wild-type strain considered. Some intron mutations in this gene affect not only its expression but also that of another mitochondrial gene: oxi3. To understand better the function of introns in gene expression, we have constructed a series of new strains that differ only by the presence or absence of one of the five wild-type introns in the cytochrome b gene, the rest of the mitochondrial and nuclear genome remaining unchanged. All constructions result from in vivo recombination events between rho- donor and rho+ recipient mtDNA. The following genes have been constructed: [see text]. Interestingly, all the genes lead to the synthesis of cytochrome b, while only the genes having the intron bI4 allow the expression of oxi3. A nuclear gene, when mutated, can compensate for the absence of the intron bI4.     

A re-evaluation of basal phylogenetic relationships within trogons (Aves: Trogonidae) based on nuclear DNA sequences

The avian clade Trogonidae (trogons) consists of approximately 40 species distributed pantropically in the Neotropical, Afrotropical and Indomalayan zoogeographical regions. In this study, we evaluate the basal phylogenetic relationships within the trogons based on DNA sequences from three nuclear introns [myoglobin intron 2, β -fibrinogen intron 7 and glyceraldehydes-3-phosphodehydrogenase (G3PDH) intron 11]. In addition, previously published cytochrome b and 12S sequences were re-analysed and combined with the nuclear data set. The analysis of the three nuclear genes combined suggests a sister group relationship between the Afrotropical ( Apaloderma ) and Indomalayan ( Harpactes ) clades, whereas the Neotropical taxa ( Trogon , Pharomachrus , and Priotelus ) form an unresolved polytomy basal to these two groups. In addition, two of the three individual gene trees also support a sister group relationship between the Afrotropical and Indomalayan trogons. This is at odds with previously published studies based on mitochondrial sequence data and DNA–DNA hybridization. The third nuclear intron (G3PDH), however, suggests that the Afrotropical trogons are basal relative the other trogons. This was also suggested by the mitochondrial data set, as well as the analysis of the combined nuclear and mitochondrial data. Both of these conflicting hypotheses are supported by high posterior probabilities. An insertion in β -fibrinogen further supports a basal position of the Afrotropical clade. Analyses of the myoglobin intron with additional outgroups place the root differently and strongly support monophyly of each of the zoogeographical regions (including the Neotropics), and these three clades form a basal trichotomy. This suggests that that rooting is a serious problem in resolving basal phylogenetic relationships among the trogons.     

Patterns of nucleotide substitution in angiosperm cpDNA trnL (UAA)-trnF (GAA) regions

Patterns of substitution in chloroplast encoded trnL_F regions were compared between species of Actaea (Ranunculales), Digitalis (Scrophulariales), Drosera (Caryophyllales), Panicoideae (Poales), the small chromosome species clade of Pelargonium (Geraniales), each representing a different order of flowering plants, and Huperzia (Lycopodiales). In total, the study included 265 taxa, each with > 900-bp sequences, totaling 0.24 Mb. Both pairwise and phylogeny-based comparisons were used to assess nucleotide substitution patterns. In all six groups, we found that transition/transversion ratios, as estimated by maximum likelihood on most-parsimonious trees, ranged between 0.8 and 1.0 for ingroups. These values occurred both at low sequence divergences, where substitutional saturation, i.e., multiple substitutions having occurred at the same (homologous) nucleotide position, was not expected, and at higher levels of divergence. This suggests that the angiosperm trnL-F regions evolve in a pattern different from that generally observed for nuclear and animal mtDNA (transitional/transversion ratio > or = 2). Transition/transversion ratios in the intron and the spacer region differed in all alignments compared, yet base compositions between the regions were highly similar in all six groups. A>-C transversions were significantly less frequent than the other four substitution types. This correlates with results from studies on fidelity mechanisms in DNA replication that predict A<->T and G<->C transversions to be least likely to occur. It therefore strengthens confidence in the link between mutation bias at the polymerase level and the actual fixation of substitutions as recorded on evolutionary trees, and concomitantly, in the neutrality of nucleotide substitutions as phylogenetic markers.     

Synonymous substitution rates in Drosophila: Mitochondrial versus nuclear genes

Synonymous substitution rates in mitochondrial and nuclear genes of Drosophila were compared. To make accurate comparisons, we considered the following: (1) relative synonymous rates, which do not require divergence time estimates, should be used; (2) methods estimating divergence should take into account base composition; (3) only very closely related species should be used to avoid effects of saturation; (4) the heterogeneity of rates should be examined. We modified the methods estimating synonymous substitution numbers to account for base composition bias. By using these methods, we found that mitochondrial genes have 1.7–3.4 times higher synonymous substitution rates than the fastest nuclear genes or 4.5–9.0 times higher rates than the average nuclear genes. The average rate of synonymous transversions was 2.7 (estimated from the melanogaster species subgroup) or 2.9 (estimated from the obscura group) times higher in mitochondrial genes than in nuclear genes. Synonymous transversions in mitochondrial genes occurred at an approximately equivalent rate to those in the fastest nuclear genes. This last result is not consistent with the hypothesis that the difference in turnover rates between mitochondrial and nuclear genomes is the major factor determining higher synonymous substitution rates in mtDNA. We conclude that the difference in synonymous substitution rates is due to a combination of two factors: a higher transitional mutation rate in mtDNA and constraints on nuclear genes due to selection for codon usage. Received: 27 November 1996 / Accepted: 8 May 1997     

Molecular basis of the 'box effect', A maturase deficiency leading to the absence of splicing of two introns located in two split genes of yeast mitochondrial DNA

In the mitochondrial DNA of Saccharomyces cerevisiae, the genes cob-box and oxi3, coding for apocytochrome b and cytochrome oxidase subunit I respectively, are split. Several mutations located in the introns of the cob-box gene prevent the synthesis of cytochrome b and cytochrome oxidase subunit I (this is known as the 'box effect').-We have elucidated the molecular basis of this phenomenon: these mutants are unable to excise the fourth intron of oxi3 from the cytochrome oxidase subunit I pre-mRNA; the absence of a functional bI4 mRNA maturase, a trans-acting factor encoded by the fourth intron of the cob-box gene explains this phenomenon. This maturase was already known to control the excision of the bI4 intron; consequently we have demonstrated that it is necessary for the processing of two introns located in two different genes. Mutations altering this maturase can be corrected, but only partially, by extragenic suppressors located in the mitochondrial (mim2) or in the nuclear (NAM2) genome. The gene product of these two suppressors should, therefore, control (directly or indirectly) the excision of the two introns as the bI4 mRNA maturase normally does.     

Molecular evolution of cytochrome c oxidase subunit I in primates: is there coevolution between mitochondrial and nuclear genomes?

Phylogenetic analyses carried out on cytochrome c oxidase (COX) subunit I mitochondrial genes from 14 primates representing the major branches of the order and four outgroup nonprimate eutherians revealed that transversions and amino acid replacements (i.e., the more slowly occurring sequence changes) contained lower levels of homoplasy and thus provided more accurate information on cladistic relationships than transitions (i.e., the more rapidly occurring sequence changes). Several amino acids, each with a high likelihood of functionality involving the binding of cytochrome c or interaction with COX VIII, have changed in Anthropoidea, the primate suborder grouping New World monkey, Old World monkey, ape, and human lineages. They are conserved in other mammalian lineages and in nonanthropoid primates. Maximum-likelihood ancestral COX I nucleotide sequences were determined utilizing a near most parsimonious branching arrangement for the primate sequences that was consistent with previously hypothesized primate cladistic relationships based on larger and more diverse data sets. Relative rate tests of COX I mitochondrial sequences showed an elevated nonsynonymous (N) substitution rate for anthropoid-nonanthropoid comparisons. This finding for the largest mitochondrial (mt) DNA-encoded subunit is consistent with previous observations of elevated nonsynonymous substitution/synonymous substitution (S) rates in primates for mt-encoded COX II and for the nuclear-encoded COX IV and COX VIIa-H. Other COX-related proteins, including cytochrome c and cytochrome b, also show elevated amino acid replacement rates or N/S during similar time frames, suggesting that this group of interacting genes is likely to have coevolved during primate evolution.     

Comparative evolution of the mitochondrial cytochrome <Emphasis Type="Italic">b</Emphasis> gene and nuclear <Emphasis Type="Italic">S7</Emphasis> ribosomal protein gene intron 1 in sinipercid fishes and their relatives

The extant sinipercids are a group of freshwater percoid fishes endemic to East Asia. A recent mitochondrial cytochrome b phylogeny of sinipercids has challenged some aspects for their traditional taxonomy and molecular phylogeny, especially for the monophyly of Sinipercidae. In this study, we analyzed mitochondrial cytochrome b and nuclear encoded S7 ribosomal protein gene intron 1 for 10 sinipercid species and 11 related species to compare the phylogenetic signal and nucleotide substitution properties of these two gene sequences. The length of S7 intron 1 ranged from 461 to 719 bp, but alignment was not difficult, and the indels, the proportion of which in the total nucleotides ranged from 3.76 to 45.83%, were phylogenetically informative. Our results indicate that: (1) the relative rate presented by cyt b is five times that of S7 intron 1; (2) the proportion of phylogenetic information is higher in S7 than in cyt b; (3) S7 intron 1 has more base composition bias, but more uniform nucleotide substitution properties; (4) the overall ratio between transitions and transversions in S7 intron 1 is lower than in cyt b. Maximum parsimony and Bayesian analyses of aligned S7 intron 1 and the combined S7 and cyt b dataset resulted in phylogenies that contained the previously identified genera Siniperca and Coreoperca, whereas the monophyly of Coreoperca cannot be corroborated by separate cyt b analysis. The monophyly of Sinipercidae is not supported in separate and combined dataset analyses, although the alternative hypothesis cannot be significantly rejected based on approximately unbiased tests and Shimodaira–Hasegawa tests. Overall, maximum parsimony analyses result in trees with a lack of phylogenetic resolution in deep nodes, and the signal from S7 intron 1 conflicts the cyt b signal in the combined dataset analyses. The reasons for the poor performance of cyt b to S7 intron 1 in the phylogeny are discussed.     

Molecular phylogeny and evolution of primate mitochondrial DNA

We determined nucleotide sequences of homologous 0.9-kb fragments of mitochondrial DNAs (mtDNAs) derived from four species of old-world monkeys, one species of new-world monkeys, and two species of prosimians. With these nucleotide sequences and homologous sequences for five species of hominoids, we constructed a phylogenetic tree for the four groups of primates. The phylogeny obtained is generally consistent with evolutionary trees constructed in previous studies. Our results also suggest that the rate of nucleotide substitution for mtDNAs in hominines (human, chimpanzee, and gorilla) may have slowed down compared with that for old-world monkeys. This evolutionary feature of mitochondrial genes is similar to one found in nuclear genes.      

Incipient mitochondrial evolution in yeasts. II. The complete sequence of the gene coding for cytochrome b in Saccharomyces douglasii reveals the presence of both new and conserved introns and discloses major differences in the fixation of mutations in evolution.

We have determined the complete sequence of the mitochondrial gene coding for cytochrome b in Saccharomyces douglasii. The gene is 6310 base-pairs long and is interrupted by four introns. The first one (1311 base-pairs) belongs to the group ID of secondary structure, contains a fragment open reading frame with a characteristic GIY ... YIG motif, is absent from Saccharomyces cerevisiae and is inserted in the same site in which introns 1 and 2 are inserted in Neurospora crassa and Podospora anserina, respectively. The next three S. douglasii introns are homologous to the first three introns of S. cerevisiae, are inserted at the same positions and display various degrees of similarity ranging from an almost complete identity (intron 2 and 4) to a moderate one (intron 3). We have compared secondary structures of intron RNAs, and nucleotide and amino acid sequences of cytochrome b exons and intron open reading frames in the two Saccharomyces species. The rules that govern fixation of mutations in exon and intron open reading frames are different: the relative proportion of mutations occurring in synonymous codons is low in some introns and high in exons. The overall frequency of mutations in cytochrome b exons is much smaller than in nuclear genes of yeasts, contrary to what has been found in vertebrates, where mitochondrial mutations are more frequent. The divergence of the cytochrome b gene is modular: various parts of the gene have changed with a different mode and tempo of evolution.     

Comparative evolution of the mitochondrial cytochrome b gene and nuclear beta-fibrinogen intron 7 in woodpeckers

Most molecular phylogenetic studies of vertebrates have been based on DNA sequences of mitochondrial-encoded genes. MtDNA evolves rapidly and is thus particularly useful for resolving relationships among recently evolved groups. However, it has the disadvantage that all of the mitochondrial genes are inherited as a single linkage group so that only one independent gene tree can be inferred regardless of the number of genes sequenced. Introns of nuclear genes are attractive candidates for independent sources of rapidly evolving DNA: they are pervasive, most of their nucleotides appear to be unconstrained by selection, and PCR primers can be designed for sequences in adjacent exons where nucleotide sequences are conserved. We sequenced intron 7 of the beta-fibrinogen gene (beta-fibint7) for a diversity of woodpeckers and compared the phylogenetic signal and nucleotide substitution properties of this DNA sequence with that of mitochondrial-encoded cytochrome b (cyt b) from a previous study. A few indels (insertions and deletions) were found in the beta-fibint7 sequences, but alignment was not difficult, and the indels were phylogentically informative. The beta-fibint7 and cyt b gene trees were nearly identical to each other but differed in significant ways from the traditional woodpecker classification. Cyt b evolves 2.8 times as fast as beta-fibint7 (14. 0 times as fast at third codon positions). Despite its relatively slow substitution rate, the phylogenetic signal in beta-fibint7 is comparable to that in cyt b for woodpeckers, because beta-fibint7 has less base composition bias and more uniform nucleotide substitution probabilities. As a consequence, compared with cyt b, beta-fibint7 nucleotide sites are expected to enter more distinct character states over the course of evolution and have fewer multiple substitutions and lower levels of homoplasy. Moreover, in contrast to cyt b, in which nearly two thirds of nucleotide sites rarely vary among closely related taxa, virtually all beta-fibint7 nucleotide sites appear free of selective constraints, which increases informative sites per unit sequenced. However, the estimated gamma distribution used to model rate variation among sites suggests constraints on some beta-fibint7 sites. This study suggests that introns will be useful for phylogenetic studies of recently evolved groups.     

The untranslated leader of nuclear COX4 gene of Saccharomyces cerevisiae contains an intron.

The nuclear gene for subunit IV of cytochrome oxidase (COX4) in Saccharomyces cerevisiae contains a 342 bp intron which is contained entirely within the 5' leader of the message. Splicing of the intron results in removal of several small open reading frames; subsequently, the COX4 AUG becomes the 5' proximal initiation codon. A strain with an rna2- mutation fails to splice mRNA efficiently at restrictive temperature and was used to map the intron splice junctions by RNase protection. Two major mRNA initiation sites were mapped by primer extension of synthetic oligodeoxynucleotides. The splice junctions and internal TACTAAC box conform to consensus sequences previously determined from other yeast introns. One gene for subunit V of cytochrome oxidase (COX5b) has also been shown to contain an intron. The significance of introns in two nuclear genes encoding subunits of cytochrome oxidase is discussed.     

Nuclear and mtDNA phylogenies of the Trimeresurus complex: implications for the gene versus species tree debate

Phylogenies based on mitochondrial DNA (mtDNA) may represent gene trees that may not be congruent with the equivalent species tree. One solution to this problem is to include additional, independent loci from the nuclear genome. Sequence data from the seventh intron of the beta-fibrinogen gene were generated for 25 specimens of vipers, including 8 nominal species of the Trimeresurus complex of Asian pit vipers. Phylogenetic trees were generated using maximum-parsimony and maximum-likelihood methods. The taxonomic level at which the intron provided significant phylogenetic information was examined and the trees were compared to those produced from previously obtained mtDNA cytochrome b sequences. A variety of different approaches (separate analyses, conditional data combination, and consensus) were used in an attempt to provide a sound organismal phylogeny based on both nuclear and mtDNA data sets. We discuss the implications for the gene tree-species tree debate and its particular relevance to medically important organisms.     

Additional support for Afrotheria and Paenungulata,the performance of mitochondrial versus nuclear genes,and the impact of data partitions with heterogeneous base composition

We concatenated sequences for four mitochondrial genes (12S rRNA, tRNA valine, 16S rRNA, cytochrome b) and four nuclear genes [aquaporin, alpha 2B adrenergic receptor (A2AB), interphotoreceptor retinoid-binding protein (IRBP), von Willebrand factor (vWF)] into a multigene data set representing 11 eutherian orders (Artiodactyla, Hyracoidea, Insectivora, Lagomorpha, Macroscelidea, Perissodactyla, Primates, Proboscidea, Rodentia, Sirenia, Tubulidentata). Within this data set, we recognized nine mitochondrial partitions (both stems and loops, for each of 12S rRNA, tRNA valine, and 16S rRNA; and first, second, and third codon positions of cytochrome b) and 12 nuclear partitions (first, second, and third codon positions, respectively, of each of the four nuclear genes). Four of the 21 partitions (third positions of cytochrome b, A2AB, IRBP, and vWF) showed significant heterogeneity in base composition across taxa. Phylogenetic analyses (parsimony, minimum evolution, maximum likelihood) based on sequences for all 21 partitions provide 99-100% bootstrap support for Afrotheria and Paenungulata. With the elimination of the four partitions exhibiting heterogeneity in base composition, there is also high bootstrap support (89-100%) for cow + horse. Statistical tests reject Altungulata, Anagalida, and Ungulata. Data set heterogeneity between mitochondrial and nuclear genes is most evident when all partitions are included in the phylogenetic analyses. Mitochondrial-gene trees associate cow with horse, whereas nuclear-gene trees associate cow with hedgehog and these two with horse. However, after eliminating third positions of A2AB, IRBP, and vWF, nuclear data agree with mitochondrial data in supporting cow + horse. Nuclear genes provide stronger support for both Afrotheria and Paenungulata. Removal of third positions of cytochrome b results in improved performance for the mitochondrial genes in recovering these clades.     

DNA barcoding of African fruit bats (Mammalia, Pteropodidae). The mitochondrial genome does not provide a reliable discrimination between Epomophorus gambianus and Micropteropus pusillus

Sequences of the mitochondrial cytochrome c oxidase subunit I (COI) gene have been shown to be useful for species identification in various groups of animals. However, the DNA barcoding approach has never been tested on African fruit bats of the family Pteropodidae (Mammalia, Chiroptera). In this study, the COI gene was sequenced from 120 bats collected in the Central African Republic and belonging to either Epomophorus?gambianus or Micropteropus?pusillus, two species easily diagnosed on the basis of morphological characters, such as body size, skull shape and palatal ridges. Two additional molecular markers were used for comparisons: the complete mitochondrial cytochrome b gene and the intron 7 of the nuclear β-fibrinogen (FGB) gene. Our results reveal an unexpected discordance between mitochondrial and nuclear genes. The nuclear FGB signal agrees with our morphological identifications, as the three alleles detected for E.?gambianus are divergent from the fourteen alleles found for M.?pusillus. By contrast, this taxonomic distinction is not recovered with the analyses of mitochondrial genes, which support rather a polyphyletic pattern for both species. The conflict between molecular markers is explained by multiple mtDNA introgression events from M.?pusillus into E.?gambianus or, alternatively, by incomplete lineage sorting of mtDNA haplotypes associated with positive selection on FGB alleles of M.?pusillus. Our work shows the failure of DNA barcoding to discriminate between two morphologically distinct fruit bat species and highlights the importance of using both mitochondrial and nuclear markers for taxonomic identification.     

Nuclear gene sequences resolve species phylogeny and mitochondrial introgression in Leptocarabus beetles showing trans-species polymorphisms

We studied the phylogenetic relationships among Japanese Leptocarabus ground beetles, which show extensive trans-species polymorphisms in mitochondrial gene genealogies. Simultaneous analysis of combined nuclear data with partial sequences from the long-wavelength rhodopsin, wingless, phosphoenolpyruvate carboxykinase, and 28S rRNA genes resolved the relationships among the five species, although separate analyses of these genes provided topologies with low resolution. For both the nuclear gene tree resulting from the combined data from four genes and a mitochondrial cytochrome oxidase subunit I (COI) gene tree, we applied a Bayesian divergence time estimation using a common calibration method to identify mitochondrial introgression events that occurred after speciation. Three mitochondrial lineages shared by two or three species were likely subject to introgression due to interspecific hybridization because the coalescent times for these lineages were much shorter than the corresponding speciation times estimated from nuclear gene sequences. We demonstrated that when species phylogeny is fully resolved with nuclear gene sequence data, comparative analysis of nuclear and mitochondrial gene trees can be used to infer introgressive hybridization events that might cause trans-species polymorphisms in mitochondrial gene trees.