Complete mitochondrial DNA sequence of the fat dormouse, Glis glis: further evidence of rodent paraphyly

Comparative Phylogenomic Assessment of Mitochondrial Introgression among Several Species of Chipmunks (Tamias)

Many species are not completely reproductively isolated, resulting in hybridization and genetic introgression. Organellar genomes, such as those derived from mitochondria (mtDNA) and chloroplasts, introgress frequently in natural systems; however, the forces shaping patterns of introgression are not always clear. Here, we investigate extensive mtDNA introgression in western chipmunks, focusing on species in the Tamias quadrivittatus group from the central and southern Rocky Mountains. Specifically, we investigate the role of selection in driving patterns of introgression. We sequenced 51 mtDNA genomes from six species and combine these sequences with other published genomic data to yield annotated mitochondrial reference genomes for nine species of chipmunks. Genomic characterization was performed using a series of molecular evolutionary and phylogenetic analyses to test protein-coding genes for positive selection. We fit a series of maximum likelihood models using a model-averaging approach, assessed deviations from neutral expectations, and performed additional tests to search for codons under the influence of selection. We found no evidence for positive selection among these genomes, suggesting that selection has not been the driving force of introgression in these species. Thus, extensive mtDNA introgression among several species of chipmunks likely reflects genetic drift of introgressed alleles in historically fluctuating populations.

Comprehensive species set revealing the phylogeny and biogeography of Feliformia (Mammalia, Carnivora) based on mitochondrial DNA

Extant Feliformia species are one of the most diverse radiations of Carnivora (~123 species). Despite substantial recent interest in their conservation, diversification, and systematic study, no previous phylogeny contains a comprehensive species set, and no biogeography of this group is available. Here, we present a phylogenetic estimate for Feliformia with a comprehensive species set and establish a historical biogeography based on mitochondrial DNA. Both the Bayesian and maximum likelihood phylogeny for Feliformia are elucidated in our analyses and are strongly consistent with many groups recognized in previous studies. The mitochondrial phylogenetic relationships of Felidae were for the first time successfully reconstructed in our analyses with strong supported. When divergence times and dispersal/vicariance histories were compared with historical sea level changes, four dispersal and six vicariance events were identified. These vicariance events were closely related with global sea level changes. The transgression of sea into the lowland plains between Eurasia and Africa may have caused the vicariance in these regions. A fall in the sea level during late Miocene to Pliocene produced the Bering strait land bridge, which assisted the migration of American Feliformia ancestors from Asia to North America. In contrast with the ‘sweepstakes hypothesis’, our results suggest that the climate cooling during 30–27 Ma assisted Feliformia migration from the African mainland to Madagascar by creating a short-lived ice bridge across the Mozambique Channel. Lineages-through-time plots revealed a large increase in lineages since the Mid-Miocene. During the Mid-Miocene Climatic Optimum, the ecosystems and population of Feliformia rapidly expanded. Subsequent climate cooling catalyzed immigration, speciation, and the extinction of Feliformia.

The position of tree shrews in the mammalian tree: Comparing multi-gene analyses with phylogenomic results leaves monophyly of Euarchonta doubtful

The well-accepted Euarchonta grandorder is a pruned version of Archonta nested within the Euarchontoglires (or Supraprimates) clade. At present, it includes tree shrews (Scandentia), flying lemurs (Dermoptera) and primates (Primates). Here, a phylogenomic dataset containing 1912 exons from 22 representative mammals was compiled to investigate the phylogenetic relationships within this group. Phylogenetic analyses and hypothesis testing suggested that tree shrews can be classified as a sister group to Primates or to Glires or even as a basal clade within Euarchontoglires. Further analyses of both modified and original previously published datasets found that the phylogenetic position of tree shrews is unstable. We also found that two of three exonic indels reported as synapomorphies of Euarchonta in a previous study do not unambiguously support the monophyly of such a clade. Therefore, the monophyly of both Euarchonta and Sundatheria (Dermoptera + Scandentia) are suspect. Molecular dating and divergence rate analyses suggested that the ancestor of Euarchontoglires experienced a rapid divergence, which may cause the unresolved position of tree shrews even using the whole genomic data.

Species based synonymous codon usage in fusion protein gene of Newcastle disease virus

Newcastle disease is highly pathogenic to poultry and many other avian species. However, the Newcastle disease virus (NDV) has also been reported from many non-avian species. The NDV fusion protein (F) is a major determinant of its pathogenicity and virulence. The functionalities of F gene have been explored for the development of vaccine and diagnostics against NDV. Although the F protein is well studied but the codon usage and its nucleotide composition from NDV isolated from different species have not yet been explored. In present study, we have analyzed the factors responsible for the determination of codon usage in NDV isolated from four major avian host species. The F gene of NDV is analyzed for its base composition and its correlation with the bias in codon usage. Our result showed that random mutational pressure is responsible for codon usage bias in F protein of NDV isolates. Aromaticity, GC3s, and aliphatic index were not found responsible for species based synonymous codon usage bias in F gene of NDV. Moreover, the low amount of codon usage bias and expression level was further confirmed by a low CAI value. The phylogenetic analysis of isolates was found in corroboration with the relatedness of species based on codon usage bias. The relationship between the host species and the NDV isolates from the host does not represent a significant correlation in our study. The present study provides a basic understanding of the mechanism involved in codon usage among species.

Primate phylogenetic relationships and divergence dates inferred from complete mitochondrial genomes

The origins and the divergence times of the most basal lineages within primates have been difficult to resolve mainly due to the incomplete sampling of early fossil taxa. The main source of contention is related to the discordance between molecular and fossil estimates: while there are no crown primate fossils older than 56Ma, most molecule-based estimates extend the origins of crown primates into the Cretaceous. Here we present a comprehensive mitogenomic study of primates. We assembled 87 mammalian mitochondrial genomes, including 62 primate species representing all the families of the order. We newly sequenced eleven mitochondrial genomes, including eight Old World monkeys and three strepsirrhines. Phylogenetic analyses support a strong topology, confirming the monophyly for all the major primate clades. In contrast to previous mitogenomic studies, the positions of tarsiers and colugos relative to strepsirrhines and anthropoids are well resolved. In order to improve our understanding of how fossil calibrations affect age estimates within primates, we explore the effect of seventeen fossil calibrations across primates and other mammalian groups and we select a subset of calibrations to date our mitogenomic tree. The divergence date estimates of the Strepsirrhine/Haplorhine split support an origin of crown primates in the Late Cretaceous, at around 74Ma. This result supports a short-fuse model of primate origins, whereby relatively little time passed between the origin of the order and the diversification of its major clades. It also suggests that the early primate fossil record is likely poorly sampled.

Jumping and gliding rodents: mitogenomic affinities of Pedetidae and Anomaluridae deduced from an RNA-Seq approach

An RNA-Seq strategy was used to obtain the complete set of protein-coding mitochondrial genes from two rodent taxa. Thanks to the next generation sequencing (NGS) 454 approach, we determined the complete mitochondrial DNA genome from Graphiurus kelleni (Mammalia: Rodentia: Gliridae) and partial mitogenome from Pedetes capensis (Pedetidae), and compared them with published rodent and outgroup mitogenomes. We finished the mitogenome sequencing by a series of amplicons using conserved PCR primers to fill the gaps corresponding to tRNA, rRNA and control regions. Phylogenetic analyses of the mitogenomes suggest a well-supported rodent phylogeny in agreement with nuclear gene trees. Pedetes groups with Anomalurus into the clade Anomaluromorpha, while Graphiurus branches within the squirrel-related clade. Moreover, Pedetes+Anomalurus branch with Castor into the mouse-related clade. Our study demonstrates the utility of NGS for obtaining new mitochondrial genomes as well as the importance of choosing adequate models of sequence evolution to infer the phylogeny of rodents.

Efficient sequencing of Anuran mtDNAs and a mitogenomic exploration of the phylogeny and evolution of frogs

Anura (frogs and toads) constitute over 88% of living amphibian diversity but many important questions about their phylogeny and evolution remain unresolved. For this study, we developed an efficient method for sequencing anuran mitochondrial DNAs (mtDNAs) by amplifying the mitochondrial genome in 12 overlapping fragments using frog-specific universal primer sets. Based on this method, we generated 47 nearly complete, new anuran mitochondrial genomes and discovered nine novel gene arrangements. By combining the new data and published anuran mitochondrial genomes, we assembled a large mitogenomic data set (11,007 nt) including 90 frog species, representing 39 of 53 recognized anuran families, to investigate their phylogenetic relationships and evolutionary history. The resulting tree strongly supported a paraphyletic arrangement of archaeobatrachian (=nonneobatrachian) frogs, with Leiopelmatoidea branching first, followed by Discoglossoidea, Pipoidea, and Pelobatoidea. Within Neobatrachia, the South African Heleophrynidae is the sister-taxon to all other neobatrachian frogs and the Seychelles-endemic Sooglossidae is recovered as the sister-taxon to Ranoidea. These phylogenetic relationships agree with many nuclear gene studies. The chronogram derived from two Bayesian relaxed clock methods (MultiDivTime and BEAST) suggests that modern frogs (Anura) originated in the early Triassic about 244 Ma and the appearance of Neobatrachia took place in the late Jurassic about 163 Ma. The initial diversifications of two species-rich superfamilies Hyloidea and Ranoidea commenced 110 and 133 Ma, respectively. These times are older than some other estimates by approximately 30-40 My. Compared with nuclear data, mtDNA produces compatible time estimates for deep nodes (>150 Ma), but apparently older estimates for more shallow nodes. Our study shows that, although it evolves relatively rapidly and behaves much as a single locus, mtDNA performs well for both phylogenetic and divergence time inferences and will provide important reference hypotheses for the phylogeny and evolution of frogs.

Complete mitochondrial genome of the Eurasian flying squirrel Pteromys volans (Sciuromorpha, Sciuridae) and revision of rodent phylogeny

In this study, the complete mitochondrial genome of the Eurasian flying squirrel Pteromys volans (Rodentia, Sciuromorpha, Sciuridae) was sequenced and characterized in detail. The entire mitochondrial genome of P. volans consisted of 16,513 bp and contained 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and two non-coding regions. Its gene arrangement pattern was consistent with the mammalian ground pattern. The overall base composition and AT contents were similar to those of other rodent mitochondrial genomes. The light-strand origin generally identified between tRNA ( Asn ) and tRNA ( Cys ) consisted of a secondary structure with an 11-bp stem and an 11-bp loop. The large control region was constructed of three characteristic domains, ETAS, CD, and CSB without any repeat sequences. Each domain contained ETAS1, subsequences A, B, and C, and CSB1, respectively. In order to examine phylogenetic contentious issues of the monophyly of rodents and phylogenetic relationships among five rodent suborders, here, phylogenetic analyses based on nucleotide sequence data of the 35 rodent and 3 lagomorph mitochondrial genomes were performed using the Bayesian inference and maximum likelihood method. The result strongly supported the rodent monophyly with high node confidence values (BP 100 % in ML and BPP 1.00 in BI) and also monophylies of four rodent suborders (BP 85-100 % in ML and BPP 1.00 in BI), except for Anomalumorpha in which only one species was examined here. Also, phylogenetic relationships among the five rodent suborders were suggested and discussed in detail.

The rabbit (Oryctolagus cuniculus): a model for mammalian reproduction and early embryology

The rabbit is a valuable animal model for a variety of biomedical research areas including in vitro fertilization, early embryology and organogenesis, neurophysiology, ophthalmology, and cardiovascular research. The rabbit is also used as a model for toxicology studies and analyses of drug effects on embryo and fetal development, as well as for research involving the immune system (e.g., host/pathogen interactions, asthma, arthritis, systemic lupus erythematosus). The rabbit’s most prominent use as a laboratory animal is antibody production; most biomedical research institutions thus provide the infrastructure and veterinary expertise for keeping rabbits. This article discusses current advances in applying modern life science research methods to basic and applied embryological research using the rabbit as a model for early intrauterine development in mammals, including state-of-the-art genomic tools and the sequencing data available (and planned) for the rabbit.

Estimating maximum likelihood phylogenies with PhyML

Our understanding of the origins, the functions and/or the structures of biological sequences strongly depends on our ability to decipher the mechanisms of molecular evolution. These complex processes can be described through the comparison of homologous sequences in a phylogenetic framework. Moreover, phylogenetic inference provides sound statistical tools to exhibit the main features of molecular evolution from the analysis of actual sequences. This chapter focuses on phylogenetic tree estimation under the maximum likelihood (ML) principle. Phylogenies inferred under this probabilistic criterion are usually reliable and important biological hypotheses can be tested through the comparison of different models. Estimating ML phylogenies is computationally demanding, and careful examination of the results is warranted. This chapter focuses on PhyML, a software that implements recent ML phylogenetic methods and algorithms. We illustrate the strengths and pitfalls of this program through the analysis of a real data set. PhyML v3.0 is available from (http://atgc_montpellier.fr/phyml/).

DNA replication and transcription in mammalian mitochondria

The mitochondrion was originally a free-living prokaryotic organism, which explains the presence of a compact mammalian mitochondrial DNA (mtDNA) in contemporary mammalian cells. The genome encodes for key subunits of the electron transport chain and RNA components needed for mitochondrial translation. Nuclear genes encode the enzyme systems responsible for mtDNA replication and transcription. Several of the key components of these systems are related to proteins replicating and transcribing DNA in bacteriophages. This observation has led to the proposition that some genes required for DNA replication and transcription were acquired together from a phage early in the evolution of the eukaryotic cell, already at the time of the mitochondrial endosymbiosis. Recent years have seen a rapid development in our molecular understanding of these machineries, but many aspects still remain unknown.

Automated scanning for phylogenetically informative transposed elements in rodents

Transposed elements constitute an attractive, useful source of phylogenetic markers to elucidate the evolutionary history of their hosts. Frequent and successive amplifications over evolutionary time are important requirements for utilizing their presence or absence as landmarks of evolution. Although transposed elements are well distributed in rodent taxa, the generally high degree of genomic sequence divergence among species complicates our access to presence/absence data. With this in mind we developed a novel, high-throughput computational strategy, called CPAL (Conserved Presence/Absence Locus-finder), to identify genome-wide distributed, phylogenetically informative transposed elements flanked by highly conserved regions. From a total of 232 extracted chromosomal mouse loci we randomly selected 14 of these plus 2 others from previous test screens and attempted to amplify them via PCR in representative rodent species. All loci were amplifiable and ultimately contributed 31 phylogenetically informative markers distributed throughout the major groups of Rodentia.

Site specific rates of mitochondrial genomes and the phylogeny of eutheria

BACKGROUND

Traditionally, most studies employing data from whole mitochondrial genomes to diagnose relationships among the major lineages of mammals have attempted to exclude regions that potentially complicate phylogenetic analysis. Components generally excluded are 3rd codon positions of protein-encoding genes, the control region, rRNAs, tRNAs, and the ND6 gene (encoded on the opposite strand). We present an approach that includes all the data, with the exception of the control region. This approach is based on a site-specific rate model that accommodates excessive homoplasy and that utilizes secondary structure as a reference for proper alignment of rRNAs and tRNAs.

RESULTS

Mitochondrial genomic data for 78 eutherian mammals, 8 metatherians, and 3 monotremes were analyzed with a Bayesian analysis and our site specific rate model. The resultant phylogeny revealed strong support for most nodes and was highly congruent with more recent phylogenies based on nuclear DNA sequences. In addition, many of the conflicting relationships observed by earlier mitochondrial-based analyses were resolved without need for the exclusion of large subsets of the data.

CONCLUSION

Rather than exclusion of data to minimize presumed noise associated with non-protein encoding genes in the mitochondrial genome, our results indicate that selection of an appropriate model that accommodates rate heterogeneity across data partitions and proper treatment of RNA genes can result in a mitochondrial genome-based phylogeny of eutherian mammals that is reasonably congruent with recent phylogenies derived from nuclear genes.

Housekeeping Genes for Phylogenetic Analysis of Eutherian Relationships

The molecular relationship of placental mammals has attracted great interest in recent years. However, 2 crucial and conflicting hypotheses remain, one with respect to the position of the root of the eutherian tree and the other the relationship between the orders Rodentia, Lagomorpha (rabbits, hares), and Primates. Although most mitochondrial (mt) analyses have suggested that rodents have a basal position in the eutherian tree, some nuclear data in combination with mt-rRNA genes have placed the root on the so-called African clade or on a branch that includes this clade and the Xenarthra (e.g., anteater and armadillo). In order to generate a new and independent set of molecular data for phylogenetic analysis, we have established cDNA sequences from different tissues of various mammalian species. With this in mind, we have identified and sequenced 8 housekeeping genes with moderately fast rate of evolution from 22 placental mammals, representing 11 orders. In order to determine the root of the eutherian tree, the same genes were also sequenced for 3 marsupial species, which were used as outgroup. Inconsistent with the analyses of nuclear + mt-rRNA gene data, the current data set did not favor a basal position of the African clade or Xenarthra in the eutherian tree. Similarly, by joining rodents and lagomorphs on the same basal branch (Glires hypothesis), the data set is also inconsistent with the tree commonly favored in mtDNA analyses. The analyses of the currently established sequences have helped examination of problematic parts in the eutherian tree at the same time as they caution against suggestions that have claimed that basal eutherian relationships have been conclusively settled.

LZ complexity distance of DNA sequences and its application in phylogenetic tree reconstruction

DNA sequences can be treated as finite-length symbol strings over a four-letter alphabet (A, C, T, G). As a universal and computable complexity measure, LZ complexity is valid to describe the complexity of DNA sequences. In this study, a concept of conditional LZ complexity between two sequences is proposed according to the principle of LZ complexity measure. An LZ complexity distance metric between two nonnull sequences is defined by utilizing conditional LZ complexity. Based on LZ complexity distance, a phylogenetic tree of 26 species of placental mammals (Eutheria) with three outgroup species was reconstructed from their complete mitochondrial genomes. On the debate that which two of the three main groups of placental mammals, namely Primates, Ferungulates, and Rodents, are more closely related, the phylogenetic tree reconstructed based on LZ complexity distance supports the suggestion that Primates and Ferungulates are more closely related.

Primate sociality in evolutionary context

Much work has been done to further our understanding of the mechanisms that underlie the diversity of primate social organizations, but none has addressed the limits to that diversity or the question of what causes species to either form or not form social networks. The fact that all living primates typically live in social networks makes it highly likely that the last common ancestor of living primates already lived in social networks, and that sociality formed an integral part of the adaptive nature of primate origins. A characterization of primate sociality within the wider mammalian context is therefore essential to further our understanding of the adaptive nature of primate origins. Here we determine correlates of sociality and nonsociality in rodents as a model to infer causes of sociality in primates. We found sociality to be most strongly associated with large-bodied arboreal species that include a significant portion of fruit in their diet. Fruits and other plant products, such as flowers, seeds, and young leaves, are patchily distributed in time and space and are therefore difficult to find. These food resources are, however, predictable and dependable when their location is known. Hence, membership in a social unit can maximize food exploitation if information on feeding sites is shared. Whether sociality evolved in the primate stem lineage or whether it was already present earlier in the evolution of Euarchontoglires remains uncertain, although tentative evidence points to the former scenario. In either case, frugivory is likely to have played an important role in maintaining the presence of a social lifestyle throughout primate evolution.

PCR-based approach for sequencing mitochondrial genomes of decapod crustaceans, with a practical example from kuruma prawn (Marsupenaeus japonicus)

An approach for sequencing the entire mitochondrial genomes (mitogenomes) of decapod crustaceans using 79 newly designed and 7 published polymerase chain reaction (PCR) primers is described. The approach comprises the following steps: (1) the entire mitogenome is amplified in 2 or 3 long PCRs; (2) the 86 primers are used in different combinations to amplify contiguous, overlapping short segments of the entire mitogenome with the diluted long PCR products as templates; (3) direct cycle sequencing is conducted using the short PCR products. This strategy allows a more rapid determination of decapod mitogenomic sequences than a traditional method using cloned mitochondrial DNA and primer walking strategy. As a practical example, the mitogenomic sequence for a kuruma prawn Marsupenaeus japonicus (Crustacea: Decapoda), was determined using the PCR-based approach.

Localization of neuropeptide Y receptor Y5 mRNA in the guinea pig brain by in situ hybridization

Neuropeptide Y (NPY) has prominent stimulatory effects on food intake in virtually all animals that have been studied. In mammals, the effect is primarily mediated by receptors Y1 and Y5, which seem to contribute to different aspects of feeding behavior in guinea pigs and rats/mice. Interestingly, differences in receptor distribution among mammalian species have been reported. To get a broader perspective on the role of Y5, we describe here studies of guinea pig (Cavia porcellus), a species which due to its phylogenetic position in the mammalian radiation is an interesting complement to previous studies in rat and mouse. Guinea pig brain sections were hybridized with two 35S-labeled oligonucleotides complementary to Y5 mRNA. The highest expression levels of Y5 mRNA were observed in the hippocampus and several hypothalamic and brain stem nuclei implicated in the regulation of feeding, such as the paraventricular, arcuate and ventromedial hypothalamic nuclei. This contrasts with autoradiography studies that detected low Y5-like binding in these areas, a discrepancy observed also in rat and human. Y5 mRNA expression was also seen in the striatum, in great contrast to mouse and rat. Taken together, these data show that Y5 mRNA distribution displays some interesting species differences, but that its expression in feeding centers seems to be essentially conserved among mammals, adding further support for an important role in food intake.

A memetic-aided approach to hierarchical clustering from distance matrices: application to gene expression clustering and phylogeny

We propose a heuristic approach to hierarchical clustering from distance matrices based on the use of memetic algorithms (MAs). By using MAs to solve some variants of the Minimum Weight Hamiltonian Path Problem on the input matrix, a sequence of the individual elements to be clustered (referred to as patterns) is first obtained. While this problem is also NP-hard, a probably optimal sequence is easy to find with the current advances for this problem and helps to prune the space of possible solutions and/or to guide the search performed by an actual clustering algorithm. This technique has been successfully applied to both a Branch-and-Bound algorithm, and to evolutionary algorithms and MAs. Experimental results are given in the context of phylogenetic inference and in the hierarchical clustering of gene expression data.

The root of the mammalian tree inferred from whole mitochondrial genomes

Morphological and molecular data are currently contradictory over the position of monotremes with respect to marsupial and placental mammals. As part of a re-evaluation of both forms of data we examine complete mitochondrial genomes in more detail. There is a particularly large discrepancy in the frequencies of thymine and cytosine (T-C) between mitochondrial genomes that appears to affect some deep divergences in the mammalian tree. We report that recoding nucleotides to RY-characters, and partitioning maximum-likelihood analyses among subsets of data reduces such biases, and improves the fit of models to the data, respectively. RY-coding also increases the signal on the internal branches relative to external, and thus increases the phylogenetic signal. In contrast to previous analyses of mitochondrial data, our analyses favor Theria (marsupials plus placentals) over Marsupionta (monotremes plus marsupials). However, a short therian stem lineage is inferred, which is at variance with the traditionally deep placement of monotremes on morphological data.

Revisiting the Glires concept--phylogenetic analysis of nuclear sequences

The so-called Glires hypothesis postulates a sister-group relationship between Rodentia (e.g., rat and mouse) and Lagomorpha (e.g., rabbit). Recent molecular phylogenetic analyses have yielded incongruent results, and either supported or refuted the Glires grouping. In order to study this inconsistency we have reconstructed phylogenetic trees based on data sets of 20 orthologous nuclear protein coding genes (6441 aa, sites) and 12 mitochondrial protein coding genes (3559 aa sites). The size of the nuclear data set is considerably larger than any comparable data set hitherto used to study the Glires concept. Analysis of the nuclear data strongly supported the phylogenetic tree (frog, chicken, ((rat, mouse), (rabbit, (human, (cattle, dog))))), while the mt data could not conclusively resolve the position of rabbit relative to that of human. This result was supported by all methods. Thus, the Glires hypothesis was rejected by this study.

Rodent phylogeny and a timescale for the evolution of Glires: evidence from an extensive taxon sampling using three nuclear genes

Rodentia is the largest order of placental mammals, with approximately 2,050 species divided into 28 families. It is also one of the most controversial with respect to its monophyly, relationships between families, and divergence dates. Here, we have analyzed and compared the performance of three nuclear genes (von Willebrand Factor, interphotoreceptor retinoid-binding protein, and Alpha 2B adrenergic receptor) for a large taxonomic sampling, covering the whole rodent and placental diversity. The phylogenetic results significantly support rodent monophyly, the association of Rodentia with Lagomorpha (the Glires clade), and a Glires + Euarchonta (Primates, Dermoptera, and Scandentia) clade. The resolution of relationships among rodents is also greatly improved. The currently recognized families are divided here into seven well-defined clades (Anomaluromorpha, Castoridae, Ctenohystrica, Geomyoidea, Gliridae, Myodonta, and Sciuroidea) that can be grouped into three major clades: Ctenohystrica, Gliridae + Sciuroidea, and a mouse-related clade (Anomaluromorpha, Castoridae + Geomyoidea, and Myodonta). Molecular datings based on these three genes suggest that the rodent radiation took place at the transition between Paleocene and Eocene. The divergence between rodents and lagomorphs is placed just at the K-T boundary and the first splits among placentals in the Late Cretaceous. Our results thus tend to reconcile molecular and morphological-paleontological insights.

The effects of nucleotide substitution model assumptions on estimates of nonparametric bootstrap support

The use of parameter-rich substitution models in molecular phylogenetics has been criticized on the basis that these models can cause a reduction both in accuracy and in the ability to discriminate among competing topologies. We have explored the relationship between nucleotide substitution model complexity and nonparametric bootstrap support under maximum likelihood (ML) for six data sets for which the true relationships are known with a high degree of certainty. We also performed equally weighted maximum parsimony analyses in order to assess the effects of ignoring branch length information during tree selection. We observed that maximum parsimony gave the lowest mean estimate of bootstrap support for the correct set of nodes relative to the ML models for every data set except one. For several data sets, we established that the exact distribution used to model among-site rate variation was critical for a successful phylogenetic analysis. Site-specific rate models were shown to perform very poorly relative to gamma and invariable sites models for several of the data sets most likely because of the gross underestimation of branch lengths. The invariable sites model also performed poorly for several data sets where this model had a poor fit to the data, suggesting that addition of the gamma distribution can be critical. Estimates of bootstrap support for the correct nodes often increased under gamma and invariable sites models relative to equal rates models. Our observations are contrary to the prediction that such models cause reduced confidence in phylogenetic hypotheses. Our results raise several issues regarding the process of model selection, and we briefly discuss model selection uncertainty and the role of sensitivity analyses in molecular phylogenetics.

Phylogenetic position of Eulipotyphla inferred from the cDNA sequences of pepsinogens A and C

Although to date the phylogenetic position of the provisional order Eulipotyphla has been assessed by various molecular markers, it has not been conclusively clarified due to low statistical supporting values and inconsistent results. To clarify the phylogenetic position of Eulipotyphla, we cloned cDNAs for pepsinogens A and C from five mammalian species belonging to four different orders and determined their nucleotide sequences. Molecular phylogenetic analysis based on the 1st and 2nd codon positions of the protein-coding region of cDNA sequences strongly supported the close relationship between Eulipotyphla and Chiroptera. Carnivora was found to be a sister group to these two orders. The monophyly of the order Rodentia and that of the cohort Glires (Rodentia and Lagomorpha) was also shown by the present phylogenetic trees of pepsinogens.

Sources of incongruence among mammalian mitochondrial sequences: COII, COIII, and ND6 genes are main contributors

To investigate the origins of incongruence among mammalian mitochondrial protein-coding genes, we compiled a matrix that included 13 protein-coding-genes for 41 mammals from 14 different orders. This matrix was examined for congruence using different partitioning strategies. The incongruence length difference test showed significant incongruence among the 13 gene partitions used simultaneously, and the result was not affected by third codon or transversion weighting. In the pair-wise comparisons, significant incongruence was detected between NADH:ubiquinone oxidoreductase subunit 6 gene (ND6), cytochrome oxidase subunit II (COII), or cytochrome oxidase subunit III (COIII) gene partitioned individually against the rest of the genes. Omission of any of the 14 mammalian orders alone or in combinations from the matrix did not result in a statistically significant improvement of congruence, suggesting that taxonomic sampling will not improve congruence among the data sets. However, omission of the ND6, COII, and COIII significantly improved congruence in our data matrix. Possible origins of unusual phylogenetic properties of the three genes are discussed.

Phylogeny of rodentia (Mammalia) inferred from the nuclear-encoded gene IRBP

The order Rodentia includes nearly half of all living mammalian species. Phylogenetic relationships among 22 species of rodents were investigated by use of a 1.2-kb region from exon 1 of the single-copy nuclear gene IRBP. IRBP has been extensively used for study of interordinal phylogeny in mammals, which allowed inclusion of 50 outgroup species, representing every eutherian order plus seven marsupials. Several clades were strongly supported, regardless of analytical method or inclusion/exclusion of data. These include a monophyletic Muroidea, with a clade including Spalax and Rhizomys as the first divergence; a clade uniting Zapus with Dipus, but excluding Sicista; a monophyletic Myodonta (Muroidea plus Dipodidae); and a clade including Aplodontidae as sister to Sciuridae. One bipartition, separating Hystricognathi and Geomyoidea from the remaining rodents, is strongly supported in all analyses that include third-position sites but almost completely absent from analyses that exclude third-position sites. A combination of nonstationary nucleotide composition and branch length effects may be causing all methods examined (including those using the LogDet distance) to support an incorrect conclusion when third-position sites are analyzed together with first- and second-position sites.

Molecular phylogeny and divergence time estimates for major rodent groups: evidence from multiple genes

The order Rodentia contains half of all extant mammal species, and from an evolutionary standpoint, there are persistent controversies surrounding the monophyly of the order, divergence dates for major lineages, and relationships among families. Exons of growth hormone receptor (GHR) and breast cancer susceptibility (BRCA1) genes were sequenced for a wide diversity of rodents and other mammals and combined with sequences of the mitochondrial 12S rRNA gene and previously published sequences of von Willebrand factor (vWF). Rodents exhibit rates of amino acid replacement twice those observed for nonrodents, and this rapid rate of evolution influences estimates of divergence dates. Based on GHR sequences, monophyly is supported, with the estimated divergence between hystricognaths and most sciurognaths dating to about 75 MYA. Most estimated dates of divergence are consistent with the fossil record, including a date of 23 MYA for Mus-Rattus divergence. These dates are considerably later than those derived from some other molecular studies. Among combined and separate analyses of the various gene sequences, moderate to strong support was found for several clades. GHR appears to have greater resolving power than do 12S or vWF. Despite its complete unresponsiveness to growth hormone, Cavia (and other hystricognaths) exhibits a conservative rate of change in the intracellular domain of GHR.

Mitochondrial versus nuclear gene sequences in deep-level mammalian phylogeny reconstruction

Both mitochondrial and nuclear gene sequences have been employed in efforts to reconstruct deep-level phylogenetic relationships. A fundamental question in molecular systematics concerns the efficacy of different types of sequences in recovering clades at different taxonomic levels. We compared the performance of four mitochondrial data sets (cytochrome b, cytochrome oxidase II, NADH dehydrogenase subunit I, 12S rRNA-tRNA-16S rRNA) and eight nuclear data sets (exonic regions of alpha-2B adrenergic receptor, aquaporin, ss-casein, gamma-fibrinogen, interphotoreceptor retinoid binding protein, kappa-casein, protamine, von Willebrand Factor) in recovering deep-level mammalian clades. We employed parsimony and minimum-evolution with a variety of distance corrections for superimposed substitutions. In 32 different pairwise comparisons between these mitochondrial and nuclear data sets, we used the maximum set of overlapping taxa. In each case, the variable-length bootstrap was used to resample at the size of the smaller data set. The nuclear exons consistently performed better than mitochondrial protein and rRNA-tRNA coding genes on a per-residue basis in recovering benchmark clades. We also concatenated nuclear genes for overlapping taxa and made comparisons with concatenated mitochondrial protein-coding genes from complete mitochondrial genomes. The variable-length bootstrap was used to score the recovery of benchmark clades as a function of the number of resampled base pairs. In every case, the nuclear concatenations were more efficient than the mitochondrial concatenations in recovering benchmark clades. Among genes included in our study, the nuclear genes were much less affected by superimposed substitutions. Nuclear genes having appropriate rates of substitution should receive strong consideration in efforts to reconstruct deep-level phylogenetic relationships.

Molecular evidence of an African Phiomorpha-South American Caviomorpha clade and support for Hystricognathi based on the complete mitochondrial genome of the cane rat (Thryonomys swinderianus)

The complete mitochondrial genome of an African cane rat, Thryonomys swinderianus (Rodentia, Hystricognathi), was included in a phylogenetic analysis along with 4 rodents, 14 additional eutherians, and 3 noneutherian outgroups. Monophyly of the suborder Hystricognathi, represented by the cane rat and the South American guinea pig, Cavia porcellus, was strongly supported by maximum-parsimony, neighbor-joining, and maximum-likelihood methods. The molecular-based estimate of the divergence time of Old and New World Hystricognathi (approximately 85 million years before present, MYBP) is consistent with an hypothesis of vicariance divergence due to the rifting of the African and South American continents 86-100 MYBP. Monophyly of Rodentia or the superordinal clade Glires (Rodentia and Lagomorpha) were not supported.

Interordinal relationships and timescale of eutherian evolution as inferred from mitochondrial genome data

Extensive phylogenetic analyses of the updated sequence data of mammalian mitochondrial genomes were carried out using the maximum likelihood method in order to resolve deep branchings in eutherian evolution. The divergence times in the mammalian tree were estimated by a relaxed molecular clock of the mitochondrial proteins calibrated with multiple references. A Chiroptera/Eulipotyphla (i.e. bat/mole) clade and a close relationship of this clade to Fereuungulata (Carnivora+Perissodactyla+Cetartiodactyla) were reconfirmed with high statistical significance. However, a support for a monophyly of Fereuungulata relative to the Chiroptera/Eulipotyphla clade was fragile, and we suggest that the three branchings among Carnivora, Perissodactyla, Cetartiodactyla and Chiroptera/Eulipotyphla occurred successively in a short time period, estimated to be approximately 77Myr BP. The Chiroptera/Eulipotyphla divergence was estimated to roughly coincide with the Cretaceous-Tertiary boundary (65Myr BP). The monophyly of Rodentia, the Lagomorpha/Rodentia clade (traditionally called Glires), and the Afrotheria/Xenarthra clade were preferred over alternative relationships, but the supports of these clades were not strong enough to exclude other possibilities. Although several super-order taxa of eutherians were strongly supported by the analyses of the mitochondrial genome data, the branching order in the deepest part of the eutherian tree remained ambiguous from the data presently available.

Long-branch attraction phenomenon and the impact of among-site rate variation on rodent phylogeny

The phylogenetic relationships among major lineages of rodents is one of the issues most debated by both paleontologists and molecular biologists. In the present study, we have analyzed all complete mammalian mitochondrial genomes available in the databases, including five rodent species (rat, mouse, dormouse, squirrel and guinea-pig). Phylogenetic analyses were performed on H-strand amino acid sequences by means of maximum-likelihood and on H-strand protein-coding and ribosomal genes by means of distance methods. Also, log-likelihood ratio tests were applied to different tree topologies under the assumption of rodent monophyly, paraphyly or polyphyly. The analyses significantly rejected rodent monophyly and showed the existence of two differentiated clades, one containing non-murids (dormouse, squirrel and guinea-pig) and the other containing murids (rat and mouse). Long-branch attraction between murids and the outgroups could not be responsible for the existence of two different rodent clades, as no significant differences in evolutionary rate have been observed, except in the case of the squirrel, which shows a lower rate. The impact of among-site rate variation models on the phylogeny of rodents has been evaluated using the gamma distribution model. Results have shown that relationships among rodents remained unchanged, and the general topology of the tree was not affected, even though some branches were not properly resolved, most likely due to a lack of fit between estimated and real rate heterogeneity parameters.

Diversity and phylogenetic implications of CsCl profiles from rodent DNAs

Buoyant density profiles of high-molecular-weight DNAs sedimented in CsCl gradients, i.e., compositional distributions of 50- to 100-kb genomic fragments, have revealed a clear difference between the murids so far studied and most other mammals, including other rodents. Sequence analyses have revealed other, related, compositional differences between murids and nonmurids. In the present study, we obtained CsCl profiles of 17 rodent species representing 13 families. The modal buoyant densities obtained for rodents span the full range of values observed in other eutherians. More remarkably, the skewness (asymmetry, mean - modal buoyant density) of the rodent profiles extends to values well below those of other eutherians. Scatterplots of these and related CsCl profile parameters show groups of rodent families that agree largely with established rodent taxonomy, in particular with the monophyly of the Geomyoidea superfamily and the position of the Dipodidae family within the Myomorpha. In contrast, while confirming and extending previously reported differences between the profiles of Myomorpha and those of other rodents, the CsCl data question a traditional hypothesis positing Gliridae within Myomorpha, as does the recently sequenced mitochondrial genome of dormouse. Analysis of CsCl profiles is presented here as a rapid, robust method for exploring rodent and other vertebrate systematics.

Lineage-specific evolutionary rate in mammalian mtDNA

The existence of a lineage-specific nucleotide substitution rate in mammalian mtDNA has been investigated by analyzing the mtDNA of all available species, that is, 35 complete mitochondrial genomes from 14 mammalian orders. A detailed study of their evolutionary dynamics has been carried out on both ribosomal RNA and first and second codon positions (P12) of H-strand protein-coding genes by using two different types of relative-rate tests. Results are quite congruent between ribosomal and P12 sites. Significant rate variations have been observed among orders and among species of the same order. However, rate variation does not exceed 1.8-fold between the fastest (Proboscidea and Primates) and the slowest (Perissodactyla) evolving orders. Thus, the observed mitochondrial rate variations among taxa do not invalidate the suitability of mtDNA for drawing mammalian phylogeny. Dependence of evolutionary rate differences on variations in mutation and/or fixation rates was examined. Body size, generation time, and metabolic rate were tested, and no significant correlation was observed between them and the taxon-specific evolutionary rates, most likely because the latter might be influenced by multiple overlapping variable constraints.

Molecular phylogeny of European muroid rodents based on complete cytochrome b sequences

Phylogenetic relationships among 18 species of mainly European muroid rodents that belong to three subfamilies were estimated using complete sequences of the mitochondrial cytochrome b gene. The inferred monophyly of the subfamilies Murinae (mice and rats) and Arvicolinae (voles, lemmings, and muskrats) is in agreement with previous studies. Within the Murinae, the morphology-based division of the genus Apodemus into three subgenera is supported by these DNA sequence data. The relationships among the different genera of the Murinae were generally poorly resolved, and the relationships of Micromys and Acomys to the other murine genera remained unresolved. Within the subfamily Arvicolinae, the relations of the genera Arvicola, Clethrionomys, and Microtus remained tentative with our data. However, within the Microtus group, there is a good molecular support for the phylogenetic relationships. These findings suggest that the origin of the different murine and arvicoline lineages was rapid, indicating an adaptive radiation with fast speciation.

Variance of molecular datings, evolution of rodents and the phylogenetic affinities between Ctenodactylidae and Hystricognathi

The von Willebrand factor (vWF) gene has been used to understand the origin and timing of Rodentia evolution in the context of placental phylogeny vWF exon 28 sequences of 15 rodent families and eight non-rodent eutherian clades are analysed with two different molecular dating methods (uniform clock on a linearized tree; quartet dating). Three main conclusions are drawn from the study of this nuclear exon. First, Ctenodactylidae (gundis) and Hystricognathi (e.g. porcupines, guinea-pigs, chinchillas) robustly cluster together in a newly recognized clade, named 'Ctenohystrica'. The Sciurognathi monophyly is subsequently rejected. Pedetidae (springhares) is an independent and early diverging rodent lineage, suggesting a convergent evolution of the multiserial enamel of rodent incisors. Second, molecular date estimates are here more influenced by accuracy and choice of the palaeontological temporal references used to calibrate the molecular clock than by either characters analysed (nucleotides versus amino acids) or species sampling. The caviomorph radiation at 31 million years (Myr) and the pig porpoise split at 63 Myr appear to be reciprocally compatible dates. Third, during the radiation of Rodentia, at least three lineages (Gliridae, Sciuroidea and Ctenohystrica) emerged close to the Cretaceous-Tertiary boundary, and their common ancestor separated from other placental orders in the Late Cretaceous.

Mitochondrial genes and mammalian phylogenies: increasing the reliability of branch length estimation

Since branch lengths provide important information about the timing and the extent of evolutionary divergence among taxa, accurate resolution of evolutionary history depends as much on branch length estimates as on recovery of the correct topology. However, the empirical relationship between the choice of genes to sequence and the quality of branch length estimation remains ill defined. To address this issue, we evaluated the accuracy of branch lengths estimated from subsets of the mitochondrial genome for a mammalian phylogeny with known subordinal relationships. Using maximum-likelihood methods, we estimated branch lengths from an 11-kb sequence of all 13 protein-coding genes and compared them with estimates from single genes (0.2-1.8 kb) and from 7 different combinations of genes (2-3.5 kb). For each sequence, we separated the component of the log-likelihood deviation due to branch length differences associated with alternative topologies from that due to those that are independent of the topology. Even among the sequences that recovered the same tree topology, some produced significantly better branch length estimates than others did. The combination of correct topology and significantly better branch length estimation suggests that these gene combinations may prove useful in estimating phylogenetic relationships for mammalian divergences below the ordinal level. Thus, the proper choice of genes to sequence is a critical factor for reliable estimation of evolutionary history from molecular data.

Model dependence of the phylogenetic inference: relationship among carnivores, Perissodactyls and cetartiodactyls as inferred from mitochondrial genome sequences

Some previous analysis of mitochondrial proteins strongly support the Carnivora/Perissodactyla grouping excluding Cetartiodactyla (Artiodactyla + Cetacea) as an outgroup, but the support of the hypothesis remains equivocal from the analysis of several nuclear-encoded proteins. In order to evaluate the strength of the support by mitochondrial proteins, phylogenetic relationship among Carnivora, Perissodactyla, and Cetartiodactyla was estimated with the ML method by using the updated data set of the 12 mitochondrial proteins with several alternative models. The analyses demonstrate that the phylogenetic inference depends on the model used in the ML analysis; i.e., whether the site-heterogeneity is taken into account and whether the rate parameters are estimated for each individual proteins or for the concatenated sequences. Although the analysis of concatenated sequences strongly supports the Carnivora/Perissodactyla grouping, the total evaluation of the separate analyses of individual proteins, which approximates the data better than the concatenated analysis, gives only ambiguous results, and therefore it is concluded that more data are needed to resolve this trichotomy.

Evolutionary genomics in Metazoa: the mitochondrial DNA as a model system

One of the most important aspects of mitochondrial (mt) genome evolution in Metazoa is constancy of size and gene content of mtDNA, whose plasticity is maintained through a great variety of gene rearrangements probably mediated by tRNA genes. The trend of mtDNA to maintain the same genetic structure within a phylum (e.g., Chordata) is generally accepted, although more recent reports show that a considerable number of transpositions are observed also between closely related organisms. Base composition of mtDNA is extremely variable. Genome GC content is often low and, when it increases, the two complementary bases distribute asymmetrically, creating, particularly in vertebrates, a negative GC-skew. In mammals, we have found coding strand base composition and average degree of gene conservation to be related to the asymmetric replication mechanism of mtDNA. A quantitative measurement of mtDNA evolutionary rate has revealed that each of the various components has a different evolutionary rate. Non-synonymous rates are gene specific and fall in a range comparable to that of nuclear genes, whereas synonymous rates are about 22-fold higher in mt than in nuclear genes. tRNA genes are among the most conserved but, when compared to their nuclear counterparts, they evolve 100 times faster. Finally, we describe some molecular phylogenetic reconstructions which have produced unexpected outcomes, and might change our vision of the classification of living organisms.

Animal mitochondrial genomes

Animal mitochondrial DNA is a small, extrachromosomal genome, typically approximately 16 kb in size. With few exceptions, all animal mitochondrial genomes contain the same 37 genes: two for rRNAs, 13 for proteins and 22 for tRNAs. The products of these genes, along with RNAs and proteins imported from the cytoplasm, endow mitochondria with their own systems for DNA replication, transcription, mRNA processing and translation of proteins. The study of these genomes as they function in mitochondrial systems-'mitochondrial genomics'-serves as a model for genome evolution. Furthermore, the comparison of animal mitochondrial gene arrangements has become a very powerful means for inferring ancient evolutionary relationships, since rearrangements appear to be unique, generally rare events that are unlikely to arise independently in separate evolutionary lineages. Complete mitochondrial gene arrangements have been published for 58 chordate species and 29 non-chordate species, and partial arrangements for hundreds of other taxa. This review compares and summarizes these gene arrangements and points out some of the questions that may be addressed by comparing mitochondrial systems.

How good are deep phylogenetic trees?

Great interest is given to species emerging early in phylogenetic reconstruction because they are often assumed to represent an ancestor. Recent studies indicate, however, that species branching deep in molecular trees are often fast-evolving ones, misplaced because of the long-branch artefact. The detection of genuinely deep-branching organisms remains an elusive task.