Emmanuel Paradis

APE: Analyses of Phylogenetics and Evolution in R language

UNLABELLEDnAnalysis of Phylogenetics and Evolution (APE) is a package written in the R language for use in molecular evolution and phylogenetics. APE provides both utility functions for reading and writing data and manipulating phylogenetic trees, as well as several advanced methods for phylogenetic and evolutionary analysis (e.g. comparative and population genetic methods). APE takes advantage of the many R functions for statistics and graphics, and also provides a flexible framework for developing and implementing further statistical methods for the analysis of evolutionary processes.nnnAVAILABILITYnThe program is free and available from the official R package archive at http://cran.r-project.org/src/contrib/PACKAGES.html#ape. APE is licensed under the GNU General Public License.

Beyond the Mediterranean peninsulas: evidence of central European glacial refugia for a temperate forest mammal species, the bank vole (Clethrionomys glareolus)

This study details the phylogeographic pattern of the bank vole, Clethrionomys glareolus, a European rodent species strongly associated with forest habitat. We used sequences of 1011 base pairs of the mitochondrial DNA cytochrome b gene from 207 bank voles collected in 62 localities spread throughout its distribution area. Our results reveal the presence of three Mediterranean (Spanish, Italian and Balkan) and three continental (western, eastern and ‘Ural’) phylogroups. The endemic Mediterranean phylogroups did not contribute to the postglacial recolonization of much of the Palaearctic range of species. Instead, the major part of this region was apparently recolonized by bank voles that survived in glacial refugia in central Europe. Moreover, our phylogeographic analyses also reveal differentiated populations of bank voles in the Ural mountains and elsewhere, which carry the mitochondrial DNA of another related vole species, the ruddy vole (Clethrionomys rutilus). In conclusion, this study demonstrates a complex phylogeographic history for a forest species in Europe which is sufficiently adaptable that, facing climate change, survives in relict southern and northern habitats. The high level of genetic diversity characterizing vole populations from parts of central Europe also highlights the importance of such regions as a source of intraspecific genetic biodiversity.

Revealing the demographic histories of species using DNA sequences

Various methodological approaches using molecular sequence data have been developed and applied across several fields, including phylogeography, conservation biology, virology and human evolution. The aim of these approaches is to obtain predictive estimates of population history from DNA sequence data that can then be used for hypothesis testing with empirical data. This recent work provides opportunities to evaluate hypotheses of constant population size through time, of population growth or decline, of the rate of growth or decline, and of migration and growth in subdivided populations. At the core of many of these approaches is the extraction of information from the structure of phylogenetic trees to infer the demographic history of a population, and underlying nearly all methods is coalescent theory. With the increasing availability of DNA sequence data, it is important to review the different ways in which information can be extracted from DNA sequence data to estimate demographic parameters.

Mitochondrial phylogeography of the Woodmouse ( Apodemus sylvaticus ) in the Western Palearctic region

We sequenced 965 base pairs of the mitochondrial DNA cytochrome b from 102 woodmice (Apodemus sylvaticus) collected from 40 European localities. The aims of the study were to answer the following questions. (i) Did the Mediterranean peninsulas play a role as refuge for woodmice? (ii) Is genetic variability of A. sylvaticus higher in the Mediterranean region compared with northern Europe? (iii) Are the patterns of the postglacial colonization of Europe by woodmice similar to those presently recognized for other European species? The results provide a clear picture of the impact of the Quaternary glaciations on the genetic and geographical structure of the woodmouse. Our analyses indicate a higher genetic variability of woodmice in the Mediterranean peninsulas compared to northern Europe, suggesting a role of the former as refuge regions for this small mammal. An original pattern of postglacial colonization is proposed where the Iberian and southern France refuge populations colonized almost all European regions. The Sicilian population appears to be very differentiated and highly variable. This emphasizes the importance of this island as a ‘hot spot’ for the intraspecific genetic diversity of the woodmouse. Finally, woodmice in North Africa originated from southwestern Europe, most probably as a result of a recent anthropogenic introduction.

Independence between developmental stability and canalization in the skull of the house mouse

The relationship between the two components of developmental homeostasis, that is canalization and developmental stability (DS), is currently debated. To appraise this relationship, the levels and morphological patterns of interindividual variation and fluctuating asymmetry were assessed using a geometric morphometric approach applied to the skulls of laboratory samples of the house mouse. These three samples correspond to two random–bred strains of the two European subspecies of the house mouse and their F1 hybrids. The inter– and intraindividual variation levels were found to be smaller in the hybrid group compared to the parental ones, suggesting a common heterotic effect on skull canalization and DS. Both buffering mechanisms might then depend on the same genetic condition, i.e. the level of heterozygosity. However, related morphological patterns did not exhibit any congruence. In contradiction with previous studies on insect wing traits, we therefore suggest that canalization and DS may not act on the same morphological characters. The fact that this discrepancy could be related to the functional importance of the symmetry of the characters under consideration is discussed in the light of our knowledge of the genetic bases of both components of developmental homeostasis.

The evolutionary radiation of Arvicolinae rodents (voles and lemmings): relative contribution of nuclear and mitochondrial DNA phylogenies

BackgroundMitochondrial and nuclear genes have generally been employed for different purposes in molecular systematics, the former to resolve relationships within recently evolved groups and the latter to investigate phylogenies at a deeper level. In the case of rapid and recent evolutionary radiations, mitochondrial genes like cytochrome b (CYB) are often inefficient for resolving phylogenetic relationships. One of the best examples is illustrated by Arvicolinae rodents (Rodentia; Muridae), the most impressive mammalian radiation of the Northern Hemisphere which produced voles, lemmings and muskrats. Here, we compare the relative contribution of a nuclear marker – the exon 10 of the growth hormone receptor (GHR) gene – to the one of the mitochondrial CYB for inferring phylogenetic relationships among the major lineages of arvicoline rodents.ResultsThe analysis of GHR sequences improves the overall resolution of the Arvicolinae phylogeny. Our results show that the Caucasian long-clawed vole (Prometheomys schaposnikowi) is one of the basalmost arvicolines, and confirm that true lemmings (Lemmus) and collared lemmings (Dicrostonyx) are not closely related as suggested by morphology. Red-backed voles (Myodini) are found as the sister-group of a clade encompassing water vole (Arvicola), snow vole (Chionomys), and meadow voles (Microtus and allies). Within the latter, no support is recovered for the generic recognition of Blanfordimys, Lasiopodomys, Neodon, and Phaiomys as suggested by morphology. Comparisons of parameter estimates for branch lengths, base composition, among sites rate heterogeneity, and GTR relative substitution rates indicate that CYB sequences consistently exhibit more heterogeneity among codon positions than GHR. By analyzing the contribution of each codon position to node resolution, we show that the apparent higher efficiency of GHR is due to their third positions. Although we focus on speciation events spanning the last 10 million years (Myr), CYB sequences display highly saturated codon positions contrary to the nuclear exon. Lastly, variable length bootstrap predicts a significant increase in resolution of arvicoline phylogeny through the sequencing of nuclear data in an order of magnitude three to five times greater than the size of GHR exon 10.ConclusionOur survey provides a first resolved gene tree for Arvicolinae. The comparison of CYB and GHR phylogenetic efficiency supports recent assertions that nuclear genes are useful for resolving relationships of recently evolved animals. The superiority of nuclear exons may reside both in (i) less heterogeneity among sites, and (ii) the presence of highly informative sites in third codon positions, that evolve rapidly enough to accumulate synapomorphies, but slow enough to avoid substitutional saturation.

Ecological Correlates and Evolutionary Divergence in the Skull of Turtles: A Geometric Morphometric Assessment

Resource use and phylogeny are often correlated with morphological variation. Moreover, because biological shapes are often complex and evolve depending on several internal constraints, they must be assessed using integrative methods. We analyzed the morphological variation of the turtle skull in the context of an adaptive radiation. Our focus are turtles of the superfamily Testudinoidea, which are remarkably diverse, both in number of species and in ecology. In this study, we depict morphological variation in the turtle skull in three dimensions with respect to diet, phylogeny, and habitat using modern geometric morphometrics. Our study revealed that morphological specialization was related to both diet and habitat. Morphological variation is decomposed in regard of both resource use (habitat and diet) and phylogeny. Feeding mode depending on environment was suggested as a key factor determining morphological evolution and diversification of turtle skulls. Diet (especially durophagy) leads to parallel morphologies in different clades. Phylogeny seemed to constrain only localized features of the skull and remained of minor influence, because overall morphotypes, closely correlated with ecological factors, occurred in both clades. In conclusion, the adaptive radiation of the Testudinoidea is revealed to demonstrate a clear relationship between the skull shape and life style.

Modeling large-scale dispersal distances

We present an approach to analyze dispersal distance data. This approach allows one to take into account accuracy of the recorded dispersal distances. Three distributions were used, all assuming continuous space; a maximum likelihood approach was used for parameter estimation and model selection. Numerical simulations showed that our method is statistically consistent since it selected the correct model with increasing frequency when sample size increased. Ringing data on two species of tits (Parus caeruleus and Parus major) in Britain and Ireland were used to illustrate the potentialities of our method. In both species, adults dispersed significantly further than juveniles. The differences between species within an age-class were not statistically significant. In all species and age-classes, the model finally selected was the one assuming a heavy-tailed half-Cauchy distribution where long-distance dispersers are predicted to be more frequent than in the exponential model. The proposed methodology can potentially be applied to any organisms, and the model selection procedure can be used with any model of the distribution of dispersal distances (DDD). Several extensions are presented in the discussion, such as generalized linear modeling of the dispersal parameters, or interfacing with capture–recapture models.

Phylogeography of Toxoplasma gondii points to a South American origin.

Toxoplasma gondii, a protozoan found ubiquitously in mammals and birds, is the etiologic agent of toxoplasmosis, a disease causing substantial public health burden worldwide, including about 200,000 new cases of congenital toxoplasmosis each year. Clinical severity has been shown to vary across geographical regions, with South America exhibiting the highest burden. Unfortunately, the drivers of these heterogeneities are still poorly understood, and the geographical origin and historical spread of the pathogen worldwide are currently uncertain. A worldwide sample of 168 T. gondii isolates gathered in 13 populations was sequenced for five fragments of genes (140 single nucleotide polymorphisms from 3153bp per isolate). Phylogeny based on Maximum likelihood methods with estimation of the time to the most recent common ancestor (TMRCA) and geostatistical analyses were performed for inferring the putative origin of T. gondii. We show that extant strains of the pathogen likely evolved from a South American ancestor, around 1.5 million years ago, and reconstruct the subsequent spread of the pathogen worldwide. This emergence is much more recent than the appearance of ancestral T. gondii, believed to have taken place about 11 My ago, and follows the arrival of felids in this part of the world. We posit that an ancestral lineage of T. gondii likely arrived in South America with felids and that the evolution of oral infectivity through carnivorism and the radiation of felids in this region enabled a new strain to outcompete the ancestral lineage and undergo a pandemic radiation.

apex: phylogenetics with multiple genes

Genetic sequences of multiple genes are becoming increasingly common for a wide range of organisms including viruses, bacteria and eukaryotes. While such data may sometimes be treated as a single locus, in practice, a number of biological and statistical phenomena can lead to phylogenetic incongruence. In such cases, different loci should, at least as a preliminary step, be examined and analysed separately. The r software has become a popular platform for phylogenetics, with several packages implementing distance‐based, parsimony and likelihood‐based phylogenetic reconstruction, and an even greater number of packages implementing phylogenetic comparative methods. Unfortunately, basic data structures and tools for analysing multiple genes have so far been lacking, thereby limiting potential for investigating phylogenetic incongruence. In this study, we introduce the new r package apex to fill this gap. apex implements new object classes, which extend existing standards for storing DNA and amino acid sequences, and provides a number of convenient tools for handling, visualizing and analysing these data. In this study, we introduce the main features of the package and illustrate its functionalities through the analysis of a simple data set.