Inferring epidemic contact structure from phylogenetic trees.

Contact structure is believed to have a large impact on epidemic spreading and consequently using networks to model such contact structure continues to gain interest in epidemiology. However, detailed knowledge of the exact contact structure underlying real epidemics is limited. Here we address the question whether the structure of the contact network leaves a detectable genetic fingerprint in the pathogen population. To this end we compare phylogenies generated by disease outbreaks in simulated populations with different types of contact networks. We find that the shape of these phylogenies strongly depends on contact structure. In particular, measures of tree imbalance allow us to quantify to what extent the contact structure underlying an epidemic deviates from a null model contact network and illustrate this in the case of random mixing. Using a phylogeny from the Swiss HIV epidemic, we show that this epidemic has a significantly more unbalanced tree than would be expected from random mixing.

Phylogenetic relationships and divergence times among dormice (Rodentia, Gliridae) based on three nuclear genes

We examined phylogenetic relationships among six species representing three subfamilies, Glirinae, Graphiurinae and Leithiinae with sequences from three nuclear protein-coding genes (apolipoprotein B, APOB; interphotoreceptor retinoid-binding protein, IRBP; recombination-activating gene 1, RAG1). Phylogenetic trees reconstructed from maximum-parsimony (MP), maximum-likelihood (ML) and Bayesian-inference (BI) analyses showed the monophyly of Glirinae (Glis and Glirulus) and Leithiinae (Dryomys, Eliomys and Muscardinus) with strong support, although the branch length maintaining this relationship was very short, implying rapid diversification among the three subfamilies. Divergence time estimates were calculated from ML (local clock model) and Bayesian-dating method using a calibration point of 25 Myr (million years) ago for the divergence between Glis and Glirulus, and 55 Myr ago for the split between lineages of Gliridae and Sciuridae on the basis of fossil records. The results showed that each lineage of Graphiuros, Glis, Glirulus and Muscardinus dates from the Late Oligocene to the Early Miocene period, which is mostly in agreement with fossil records. Taking into account that warm climate harbouring a glirid-favoured forest dominated from Europe to Asia during this period, it is considered that this warm environment triggered the prosperity of the glirid species through the rapid diversification. Glirulus japonicas is suggested to be a relict of this ancient diversification during the warm period.

Improved multiple-locus variable-number tandem-repeat assay for Staphylococcus aureus genotyping, providing a highly informative technique together with strong phylogenetic value.

We describe an improved multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) scheme for genotyping Staphylococcus aureus. We compare its performance to those of multilocus sequence typing (MLST) and spa typing in a survey of 309 strains. This collection includes 87 epidemic methicillin-resistant S. aureus (MRSA) strains of the Harmony collection, 75 clinical strains representing the major MLST clonal complexes (CCs) (50 methicillin-sensitive S. aureus [MSSA] and 25 MRSA), 135 nasal carriage strains (133 MSSA and 2 MRSA), and 13 published S. aureus genome sequences. The results show excellent concordance between the techniques' results and demonstrate that the discriminatory power of MLVA is higher than those of both MLST and spa typing. Two hundred forty-two genotypes are discriminated with 14 VNTR loci (diversity index, 0.9965; 95% confidence interval, 0.9947 to 0.9984). Using a cutoff value of 45%, 21 clusters are observed, corresponding to the CCs previously defined by MLST. The variability of the different tandem repeats allows epidemiological studies, as well as follow-up of the evolution of CCs and the identification of potential ancestors. The 14 loci can conveniently be analyzed in two steps, based upon a first-line simplified assay comprising a subset of 10 loci (panel 1) and a second subset of 4 loci (panel 2) that provides higher resolution when needed. In conclusion, the MLVA scheme proposed here, in combination with available on-line genotyping databases (including http://mlva.u-psud.fr/), multiplexing, and automatic sizing, can provide a basis for almost-real-time large-scale population monitoring of S. aureus.

The immunological, environmental, and phylogenetic perpetrators of metastatic leishmaniasis.

Cutaneous leishmaniases have persisted for centuries as chronically disfiguring parasitic infections affecting millions of people across the subtropics. Symptoms range from the more prevalent single, self-healing cutaneous lesion to a persistent, metastatic disease, where ulcerations and granulomatous nodules can affect multiple secondary sites of the skin and delicate facial mucosa, even sometimes diffusing throughout the cutaneous system as a papular rash. The basis for such diverse pathologies is multifactorial, ranging from parasite phylogeny to host immunocompetence and various environmental factors. Although complex, these pathologies often prey on weaknesses in the innate immune system and its pattern recognition receptors. This review explores the observed and potential associations among the multifactorial perpetrators of infectious metastasis and components of the innate immune system.

Arbuscular mycorrhizal fungi mediate below-ground plant-herbivore interactions: a phylogenetic study

Ecological interactions are complex networks, but have typically been studied in a pairwise fashion. Examining how third-party species can modify the outcome of pairwise interactions may allow us to better predict their outcomes in realistic systems. For instance, arbuscular mycorrhizal fungi (AMF) can affect plant interactions with other organisms, including below-ground herbivores, but the mechanisms underlying these effects remain unclear. Here, we use a comparative, phylogenetically controlled approach to test the relative importance of mycorrhizal colonization and plant chemical defences (cardenolides) in predicting plant survival and the abundance of a generalist below-ground herbivore across 14 species of milkweeds (Asclepias spp.). Plants were inoculated with a mixture of four generalist AMF species or left uninoculated. After 1month, larvae of Bradysia sp. (Diptera: Sciaridae), a generalist below-ground herbivore, colonized plant roots. We performed phylogenetically controlled analyses to assess the influence of AMF colonization and toxic cardenolides on plant growth, mortality and infestation by fungus gnats. Overall, plants inoculated with AMF exhibited greater survival than did uninoculated plants. Additionally, surviving inoculated plants had lower numbers of larvae in their roots and fewer non-AM fungi than surviving uninoculated plants. In phylogenetic controlled regressions, gnat density in roots was better predicted by the extent of root colonized by AMF than by root cardenolide concentration. Taken as a whole, AMF modify the effect of below-ground herbivores on plants in a species-specific manner, independent of changes in chemical defence. This study adds to the growing body of literature demonstrating that mycorrhizal fungi may improve plant fitness by conferring protection against antagonists, rather than growth benefits. In addition, we advocate using comparative analyses to disentangle the roles of shared history and ecology in shaping trait expression and to better predict the outcomes of complex multitrophic interactions.

Evolution of specialization: a phylogenetic study of host range in the red milkweed beetle (Tetraopes tetraophthalmus).

Specialization is common in most lineages of insect herbivores, one of the most diverse groups of organisms on earth. To address how and why specialization is maintained over evolutionary time, we hypothesized that plant defense and other ecological attributes of potential host plants would predict the performance of a specialist root-feeding herbivore (the red milkweed beetle, Tetraopes tetraophthalmus). Using a comparative phylogenetic and functional trait approach, we assessed the determinants of insect host range across 18 species of Asclepias. Larval survivorship decreased with increasing phylogenetic distance from the true host, Asclepias syriaca, suggesting that adaptation to plant traits drives specialization. Among several root traits measured, only cardenolides (toxic defense chemicals) correlated with larval survival, and cardenolides also explained the phylogenetic distance effect in phylogenetically controlled multiple regression analyses. Additionally, milkweed species having a known association with other Tetraopes beetles were better hosts than species lacking Tetraopes herbivores, and milkweeds with specific leaf area values (a trait related to leaf function and habitat affiliation) similar to those of A. syriaca were better hosts than species having divergent values. We thus conclude that phylogenetic distance is an integrated measure of phenotypic and ecological attributes of Asclepias species, especially defensive cardenolides, which can be used to explain specialization and constraints on host shifts over evolutionary time.

Optimization strategies for fast detection of positive selection on phylogenetic trees.

MOTIVATION: The detection of positive selection is widely used to study gene and genome evolution, but its application remains limited by the high computational cost of existing implementations. We present a series of computational optimizations for more efficient estimation of the likelihood function on large-scale phylogenetic problems. We illustrate our approach using the branch-site model of codon evolution. RESULTS: We introduce novel optimization techniques that substantially outperform both CodeML from the PAML package and our previously optimized sequential version SlimCodeML. These techniques can also be applied to other likelihood-based phylogeny software. Our implementation scales well for large numbers of codons and/or species. It can therefore analyse substantially larger datasets than CodeML. We evaluated FastCodeML on different platforms and measured average sequential speedups of FastCodeML (single-threaded) versus CodeML of up to 5.8, average speedups of FastCodeML (multi-threaded) versus CodeML on a single node (shared memory) of up to 36.9 for 12 CPU cores, and average speedups of the distributed FastCodeML versus CodeML of up to 170.9 on eight nodes (96 CPU cores in total).Availability and implementation: ftp://ftp.vital-it.ch/tools/FastCodeML/. CONTACT: selectome@unil.ch or nicolas.salamin@unil.ch.

Evidence for adaptive radiation from a phylogenetic study of plant defenses.

One signature of adaptive radiation is a high level of trait change early during the diversification process and a plateau toward the end of the radiation. Although the study of the tempo of evolution has historically been the domain of paleontologists, recently developed phylogenetic tools allow for the rigorous examination of trait evolution in a tremendous diversity of organisms. Enemy-driven adaptive radiation was a key prediction of Ehrlich and Raven's coevolutionary hypothesis [Ehrlich PR, Raven PH (1964) Evolution 18:586-608], yet has remained largely untested. Here we examine patterns of trait evolution in 51 North American milkweed species (Asclepias), using maximum likelihood methods. We study 7 traits of the milkweeds, ranging from seed size and foliar physiological traits to defense traits (cardenolides, latex, and trichomes) previously shown to impact herbivores, including the monarch butterfly. We compare the fit of simple random-walk models of trait evolution to models that incorporate stabilizing selection (Ornstein-Ulenbeck process), as well as time-varying rates of trait evolution. Early bursts of trait evolution were implicated for 2 traits, while stabilizing selection was implicated for several others. We further modeled the relationship between trait change and species diversification while allowing rates of trait evolution to vary during the radiation. Species-rich lineages underwent a proportionately greater decline in latex and cardenolides relative to species-poor lineages, and the rate of trait change was most rapid early in the radiation. An interpretation of this result is that reduced investment in defensive traits accelerated diversification, and disproportionately so, early in the adaptive radiation of milkweeds.

Phylogenetic analysis of plastid and nuclear DNA sequences reveals a rapid late Miocene radiation of the temperate bamboo tribe Arundinarieae (Poaceae, Bambusoideae)

Background: Arundinarieae are a large tribe of temperate woody bamboos for which phylogenetics are poorly understood because of limited taxon sampling and lack of informative characters. Aims: This study assessed phylogenetic relationships, origins and classification of Arundinarieae. Methods: DNA sequences (plastid trnL-F; nuclear ITS) were used for parsimony and Bayesian inference including 41 woody bamboo taxa. Divergence dates were estimated using a relaxed Bayesian clock. Results: Arundinarieae were monophyletic but their molecular divergence was low compared to the tropical Bambuseae. Ancestors of the Arundinarieae lineage were estimated to have diverged from the other bamboos 23 (15-30) million years ago (Mya). However, the Arundinarieae radiation occurred 10 (6-16) Mya compared to 18 (11-25) Mya for the tropical Bambuseae. Some groups could be defined within Arundinarieae, but they do not correspond to recognised subtribes such as Arundinariinae or Shibataeinae. Conclusions: Arundinarieae are a relatively ancient bambusoid lineage that underwent a rapid radiation in the late Miocene. The radiation coincides with the continental collision of the Indo-Australian and Eurasian Plates. Arundinarieae are distributed primarily in East Asia and the Himalayas to northern Southeast Asia. It is unknown whether they were present in Asia long before their radiation, but we believe recent dispersal is a more likely scenario. Keywords: Arundinarieae; Bambuseae; internal transcribed spacer (ITS); molecular clock; phylogenetics; radiation; temperate bamboos; Thamnocalaminae; trnL-F

Methicillin-resistant Staphylococcus aureus: phylogenetic relatedness between European epidemic clones and Swiss sporadic strains.

We have compared the phylogenetic diversity of methicillin-resistant Staphylococcus aureus (MRSA) strains from Switzerland and their phylogenetic relationships with European epidemic clones, using multiprimer random amplification polymorphic DNA (RAPD). Strains included 24 European epidemic clones (59 strains), 66 sporadic strains isolated in Switzerland in 1996-1997, and 15 reference strains of five other Staphylococcus species. Similarity and clustering analysis with the Jaccard's coefficient showed that the maximum genetic distance between MRSA strains was 0.43, whereas the minimum genetic distance between the six Staphylococcus species was 0.97, indicating that the method permits phylogenetic hierarchization. The 24 MRSA clones reported to be epidemic in European countries during the 1990s were distributed into seven different genetic clusters with a maximum distance of 0.29 among them. This clustering pattern was confirmed by the analysis of a subset of MRSA strains by multilocus enzyme electrophoresis at 12 loci. Most of the sporadic Swiss strains were distributed into these seven different genetic clusters, together with the epidemic MRSA clones. This suggests that there is no phylogenetic cluster specific to epidemic clones of MRSA.

Phylogenetic position of Fit insect toxin producing pseudomonads

Fit produced by Pseudomonas fluorescens CHA0 is a novel insect toxin in root colonizing pseudomonads, of which a homologue is described in Photorhabdus species.However, occurrence and abundance of insect pathogenicity in plant-associated pseudomonads is still unclear. An extensive screening outside the P. fluorescens complex identified strains of Pseudomonas chlororaphis as further Fit toxin producing candidates. Sequences of five different P. chlororaphis strains generated in this study were used to reconstruct the evolutionary history of the Fit toxin gene and to analyse its mode of evolution. We found that P. chlororaphis is closely associated with a small subgroup of 2,4-diacetylphloroglucinol and pyoluteorin- producing pseudomonads, both when analyzing four housekeeping genes and the nucleotide sequences for the Fit toxin gene. Additionally, we identified purifying selection to be the predominant mode of Fit toxin evolution.

Phylogenetic study of the Grimmia Hedw. (Grimmiales, Bryopsida) based on a combination of morphological and molecular characters

Résumé Les mousses sont la plus ancienne lignée de plantes terrestres et leur longue évolution a été accompagnée par des tendances à la simplification des caractères morphologiques. Ce phénomène a quelque peu compliqué les reconstructions phylogénétiques basées sur la morphologie. Les analyses génétiques ont permis de donner de nouvelles informations dans le cadre des analyses phylogénétiques et une réévaluation de certains caractères morphologiques. La plupart des études combinant les données morphologiques et moléculaires ne concernent que des niveaux systématiques élevés comme l'ordre ou la famille et très peu considèrent le niveau du genre. La présente étude tend à tester les relations phylogénétiques du genre Grimmia à l'aide d'une combinaison de caractères morphologiques et moléculaires. Les 40 espèces de Grimmia utilisées dans la première partie de cette étude représentent la majorité des espèces trouvées en Eurasie, un des centres de diversification du genre. Lors de l'analyse morphologique, 52 caractères morphologiques/anatomiques (33 du gamétophyte et 19 du sporophyte) ont été numérisés. Malgré le peu de support statistique des arbres, la topologie des arbres est stable. Les Grimmia, comme décrit précédemment, sont paraphylétiques. Trois clades, correspondant respectivement aux sous-genres Rhabdogrimmia Limpr, Litoneuron I.Hagen et Gasterogrimmia Schimp. sont présents, tandis que le restant des taxons appartenant aux Grimmia forment un groupe non-résolu et indistinct des autres Grimmiaceae. Les séquences chloroplastiques trnL-trnF et rps4 combinés à la morphologie ont été ensuite utilisés pour reconstruire la phylogénie des Grimmia. Les arbres obtenus soutiennent la monophylie des Grimmiaceae tandis que les Grimmia, sont paraphylétiques. Deux clades principaux correspondant aux "Rhabdogrimmia" et aux "Grimmia" se détachent. Seules les espèces de "Rhabdogrimmia" produisent des gemmules foliaires (reproduction asexuée). Dans une étude considèrant 91 séquences trrIL-trnF les espèces appartenant aux "Rhabdogrimrnia" (reproduction asexuée essentiellement) ont des variabilités intraspécifique très faible et interspécifique relativement élevée tandis que les "Grimmia" possèdent la tendance inverse (plus de reproduction sexuée). Summary The mosses are a very old land plant lineage and their long evolutionary history has been accompanied by a trend of morphological character simplifications. This phenomenon has somewhat complicated morphological based phylogenetic reconstructions. Genetic analyses have provided new insights for phylogenetic studies, and have allowed morphological data to be re¬evaluated. Most of the studies combining morphological and molecular data have concerned the higher systematic levels of order and family and only have few considered the genus. The present study aims to test the phylogenetic relationships of the genus Grimmia using a combination of morphological and molecular characters. The 40 chosen Grimmia species represent the majority of those found in Eurasia, one diversification centers of the genus. For the morphological analysis, 52 morphological/anatomical characters (33 gametophyte and 19 sporophyte characters) were numerized. Although the internal statistical support was relatively low, the tree topologies were stable. Grimmia as currently defined was found to be paraphyletic. Three subclades, corresponding to the subgenera Rhabdogrimmia Limpr., Litoneuron I.Hagen, and Gasterogrimmia Schimp. were observed in the trees, while the reminder of the Grimmia species formed an unresolved group indistinct from other Grimmiaceae. Chloroplast (trnL-trnF and rps4) DNA sequences combined with morphology were used to reconstruct the phylogeny of Grimmia. The resulting trees supported the monophyly of Grimmiaceae and that the genus Grimmia, as currently defined, as paraphyletic. Two main clades were resolved corresponding to "Rhabdogrimmia" and "Grimmia". The species belonging to "Rhabdogrimmia" produce foliar-gemmae (asexual reproduction). In a study using 91 sequences of trnL-trnF,"Rhabdogrimmia" species (mainly asexual reproduction) have very low intraspecific variability and high interspecific variability whereas the "Grimmia" species possess the inverse tendency.

Phylogenetic relationships of Apodemus Kaup, 1829 (Rodentia: Muridae) species in the Eastern Mediterranean inferred from mitochondrial DNA, with emphasis on iranian species

Wood mice of the genus Apodemus are widely distributed in Eurasia, with the Eastern Mediterranean being considered as a hotspot. Indeed, numerous species have been documented in Iran, including A. witherbyi, A. hyrcanicus, A. uralensis, A. avicennicus, A. hermonensis, and A. arianus. In this study, 129 specimens were collected from different Iranian localities and two specimens from Afghanistan. The animals were identified taxonomically and their phylogenetic relationships were investigated using cytochrome b mitochondrial DNA sequences. Five species of the genus Apodemus were identified in Iran, including A. hyrcanicus, A. witherbyi, A. cf. ponticus, A. uralensis, and A. mystacinus, beside, A. pallipes from Afghanistan. This study found no evidence of A. flavicollis or A. sylvaticus in Iran, despite their occurrence in Turkey, shedding doubt on the status of A. flavicollis in Iran, Asia Minor, and the Levant. Phylogenetic analyses imply that A. witherbyi has priority over A. avicennicus, A. hermonensis, and A. iconicus. Estimation of the divergence time for these taxa suggests a separation at around 7.2 Ma for the subgenera Karstomys (including A. mystacinus and A. epimelas) and Sylvaemus (including A. flavicollis, A. sylvaticus, A. uralensis, A. pallipes, A. hyrcanicus, A. witherbyi, and A. cf. ponticus). Within the subgenus Karstomys, the divergence times for A. mystacinus and A. epimelas were between 3.0 and 6.1 Ma, and divergence times within the subgenus Sylvaemus were between 5.2 and 6.9 Ma for A. witherbyi and other species in this subgenus. It is postulated that vicariance events including the uplifting of the Zagros Mountains and Anatolian Plateau in the middle Miocene and climate oscillations during the Messinian Salinity Crisis besides formation of the Hyrcanian tertiary forests during the Neogene probably played substantial roles in the radiation and distribution of the genus Apodemus in the Eastern Mediterranean.

Coev-web: a web platform designed to simulate and evaluate coevolving positions along a phylogenetic tree.

BACKGROUND: Available methods to simulate nucleotide or amino acid data typically use Markov models to simulate each position independently. These approaches are not appropriate to assess the performance of combinatorial and probabilistic methods that look for coevolving positions in nucleotide or amino acid sequences. RESULTS: We have developed a web-based platform that gives a user-friendly access to two phylogenetic-based methods implementing the Coev model: the evaluation of coevolving scores and the simulation of coevolving positions. We have also extended the capabilities of the Coev model to allow for the generalization of the alphabet used in the Markov model, which can now analyse both nucleotide and amino acid data sets. The simulation of coevolving positions is novel and builds upon the developments of the Coev model. It allows user to simulate pairs of dependent nucleotide or amino acid positions. CONCLUSIONS: The main focus of our paper is the new simulation method we present for coevolving positions. The implementation of this method is embedded within the web platform Coev-web that is freely accessible at http://coev.vital-it.ch/, and was tested in most modern web browsers.