Duration of rupture of membranes and microbiome transmission to the newborn: A prospective study

OBJECTIVE

To assess whether labour variables (i.e. individuals characteristics, labour characteristics and medical interventions) impact maternal and newborn microbiomes.

DESIGN

Prospective monocentric study.

SETTING

Saint-Joseph Hospital tertiary maternity unit, in Paris, France.

POPULATION

All consecutive primiparous women with a physiological pregnancy and term labour from 15 April to 1 June 2017.

METHODS

16S ribosomal RNA gene sequencing of the maternal vaginal, newborn skin and newborn oral microbiomes from 58 mother-baby dyads.

MAIN OUTCOME MEASURES

Analysis of the effects of 19 labour variables on the composition and diversity of these microbiomes.

RESULTS

The 19 labour variables explained a significant part of the variability in the vaginal, newborn oral and skin microbiomes (44%-67%). Strikingly, duration of rupture of membranes was the single factor that explained the greatest variability (adjusted R2 : 7.7%-8.4%, p ≤ 0.002) and conditioned, by itself, the compositions of the three microbiomes under study. Long duration of rupture of membranes was specifically associated with a lower relative abundance of the Lactobacillus genus (1.7-fold to 68-fold reduction, p < 0.0001) as well as an increase in microbiome diversity, including genera implicated in nosocomial infections. The effects of duration of rupture of membranes were also present in newborns delivered by non-elective caesarean section.

CONCLUSIONS

Maternal and newborn microbiomes were greatly affected by labour variables. Duration of rupture of membranes, even in non-elective caesarean sections, should be considered in epidemiological and microbiological studies, as well as in vaginal seeding practices.

Network analyses unveil ageing-associated pathways evolutionarily conserved from fungi to animals

The genetic roots of the diverse paces and shapes of ageing and of the large variations in longevity observed across the tree of life are poorly understood. Indeed, pathways associated with ageing/longevity are incompletely known, both in terms of their constitutive genes/proteins and of their molecular interactions. Moreover, there is limited overlap between the genes constituting these pathways across mammals. Yet, dedicated comparative analyses might still unravel evolutionarily conserved, important pathways associated with longevity or ageing. Here, we used an original strategy with a double evolutionary and systemic focus to analyse protein interactions associated with ageing or longevity during the evolution of five species of Opisthokonta. We ranked these proteins and interactions based on their evolutionary conservation and centrality in past and present protein-protein interaction (PPI) networks, providing a big systemic picture of the evolution of ageing and longevity pathways that identified which pathways emerged in which Opisthokonta lineages, were conserved, and/or central. We confirmed that longevity/ageing-associated proteins (LAPs), be they pro- or anti-longevity, are highly central in extant PPI, consistently with the antagonistic pleiotropy theory of ageing, and identified key antagonistic regulators of ageing/longevity, 52 of which with homologues in humans. While some highly central LAPs were evolutionarily conserved for over a billion years, we report a clear transition in the functionally important components of ageing/longevity within bilaterians. We also predicted 487 novel evolutionarily conserved LAPs in humans, 54% of which are more central than mTOR, and 138 of which are druggable, defining new potential targets for anti-ageing treatments in humans.

Antifungal Potential of the Skin Microbiota of Hibernating Big Brown Bats (Eptesicus fuscus) Infected With the Causal Agent of White-Nose Syndrome

Little is known about skin microbiota in the context of the disease white-nose syndrome (WNS), caused by the fungus Pseudogymnoascus destructans (Pd), that has caused enormous declines of hibernating North American bats over the past decade. Interestingly, some hibernating species, such as the big brown bat (Eptesicus fuscus), appear resistant to the disease and their skin microbiota could play a role. However, a comprehensive analysis of the skin microbiota of E. fuscus in the context of Pd has not been done. In January 2017, we captured hibernating E. fuscus, sampled their skin microbiota, and inoculated them with Pd or sham inoculum. We allowed the bats to hibernate in the lab under controlled conditions for 11 weeks and then sampled their skin microbiota to test the following hypotheses: (1) Pd infection would not disrupt the skin microbiota of Pd-resistant E. fuscus; and (2) microbial taxa with antifungal properties would be abundant both before and after inoculation with Pd. Using high-throughput 16S rRNA gene sequencing, we discovered that beta diversity of Pd-inoculated bats changed more over time than that of sham-inoculated bats. Still, the most abundant taxa in the community were stable throughout the experiment. Among the most abundant taxa, Pseudomonas and Rhodococcus are known for antifungal potential against Pd and other fungi. Thus, in contrast to hypothesis 1, Pd infection destabilized the skin microbiota but consistent with hypothesis 2, bacteria with known antifungal properties remained abundant and stable on the skin. This study is the first to provide a comprehensive survey of skin microbiota of E. fuscus, suggesting potential associations between the bat skin microbiota and resistance to the Pd infection and WNS. These results set the stage for future studies to characterize microbiota gene expression, better understand mechanisms of resistance to WNS, and help develop conservation strategies.

Introducing Trait Networks to Elucidate the Fluidity of Organismal Evolution Using Palaeontological Data

Explaining the evolution of animals requires ecological, developmental, paleontological, and phylogenetic considerations because organismal traits are affected by complex evolutionary processes. Modeling a plurality of processes, operating at distinct time-scales on potentially interdependent traits, can benefit from approaches that are complementary treatments to phylogenetics. Here, we developed an inclusive network approach, implemented in the command line software ComponentGrapher, and analyzed trait co-occurrence of rhinocerotoid mammals. We identified stable, unstable, and pivotal traits, as well as traits contributing to complexes, that may follow to a common developmental regulation, that point to an early implementation of the postcranial Bauplan among rhinocerotoids. Strikingly, most identified traits are highly dissociable, used repeatedly in distinct combinations and in different taxa, which usually do not form clades. Therefore, the genes encoding these traits are likely recruited into novel gene regulation networks during the course of evolution. Our evo-systemic framework, generalizable to other evolved organizations, supports a pluralistic modeling of organismal evolution, including trees and networks.

Mapping responsible conduct in the uncharted field of research-creation: A scoping review

This scoping review addresses the issues of responsible conduct of research (RCR) that can arise in the practice of research-creation (RC), an emergent, interdisciplinary, and heterogeneous field at the interface of academic research and creative activities. Little is yet known about the nature and scope of RCR issues in RC, so our study examined three questions: (1) What are the specific issues in RC in relation to RCR? (2) How does the specificity of RC influence the understanding and practice of RCR? (3) What recommendations could help address the issues highlighted in the literature? To answer these questions, we conducted a scoping review of the academic literature (n = 181 texts) dealing with RCR in RC. We found that researcher-creators faced some very different RCR challenges in comparison with their colleagues in the rest of academia. Addressing these issues is important for both the RCR and RC communities in order to ensure that the rapid development of this field occurs in line with the norms of RCR which, nonetheless, should be adapted to respect the particularities of RC and allow its contributions to the academic world.

CompositeSearch: A Generalized Network Approach for Composite Gene Families Detection

Genes evolve by point mutations, but also by shuffling, fusion, and fission of genetic fragments. Therefore, similarity between two sequences can be due to common ancestry producing homology, and/or partial sharing of component fragments. Disentangling these processes is especially challenging in large molecular data sets, because of computational time. In this article, we present CompositeSearch, a memory-efficient, fast, and scalable method to detect composite gene families in large data sets (typically in the range of several million sequences). CompositeSearch generalizes the use of similarity networks to detect composite and component gene families with a greater recall, accuracy, and precision than recent programs (FusedTriplets and MosaicFinder). Moreover, CompositeSearch provides user-friendly quality descriptions regarding the distribution and primary sequence conservation of these gene families allowing critical biological analyses of these data.

Hundreds of novel composite genes and chimeric genes with bacterial origins contributed to haloarchaeal evolution

BACKGROUND

Haloarchaea, a major group of archaea, are able to metabolize sugars and to live in oxygenated salty environments. Their physiology and lifestyle strongly contrast with that of their archaeal ancestors. Amino acid optimizations, which lowered the isoelectric point of haloarchaeal proteins, and abundant lateral gene transfers from bacteria have been invoked to explain this deep evolutionary transition. We use network analyses to show that the evolution of novel genes exclusive to Haloarchaea also contributed to the evolution of this group.

RESULTS

We report the creation of 320 novel composite genes, both early in the evolution of Haloarchaea during haloarchaeal genesis and later in diverged haloarchaeal groups. One hundred and twenty-six of these novel composite genes derived from genetic material from bacterial genomes. These latter genes, largely involved in metabolic functions but also in oxygenic lifestyle, constitute a different gene pool from the laterally acquired bacterial genes formerly identified. These novel composite genes were likely advantageous for their hosts, since they show significant residence times in haloarchaeal genomes-consistent with a long phylogenetic history involving vertical descent and lateral gene transfer-and encode proteins with optimized isoelectric points.

CONCLUSIONS

Overall, our work encourages a systematic search for composite genes across all archaeal major groups, in order to better understand the origins of novel prokaryotic genes, and in order to test to what extent archaea might have adjusted their lifestyles by incorporating and recycling laterally acquired bacterial genetic fragments into new archaeal genes.

Paternity Analysis of Wood Turtles (Glyptemys insculpta) Reveals Complex Mating Patterns

Mating system characteristics are of great importance as they may influence male and female reproductive success and reproductive isolation. The wood turtle (Glyptemys insculpta) is a terrestrial freshwater species listed as endangered by the International Union for Conservation of Nature. Considering its conservation status and the paucity of information currently available on parentage relationship for the species, we performed a microsatellite analysis to study the mating system of wood turtles in the Shawinigan River (Québec). We sampled 38 clutches over 2 years (14 in 2006 and 24 in 2007), for a total of 248 offspring genotyped with 7 microsatellite loci. The reconstructed genotypes of the fathers revealed that reproductive success in the sampled clutches varied greatly between males and are positively correlated with the number of mates and clutches sired. Frequency of multiple paternity was estimated at 37% through a consensus of 3 different estimation methods. Positive correlation was observed between the genetic diversity of clutches and the number of fathers. Repeat paternity, however, was observed in 88% of the clutches by the same female in successive years, which suggests either a frequent use of sperm storage, or remating with the same partner in successive years.

The Inuit gut microbiome is dynamic over time and shaped by traditional foods

BACKGROUND

The human gut microbiome represents a diverse microbial community that varies across individuals and populations, and is influenced by factors such as host genetics and lifestyle. Diet is a major force shaping the gut microbiome, and the effects of dietary choices on microbiome composition are well documented. However, it remains poorly known how natural temporal variation in diet can affect the microbiome. The traditional Inuit diet is primarily based on animal products, which are thought to vary seasonally according to prey availability. We previously investigated the Inuit gut microbiome sampled at a single time point, and found no detectable differences in overall microbiome community composition attributable to the traditional Inuit diet.

RESULTS

To determine whether seasonal changes in the Inuit diet might induce more pronounced changes in the microbiome, we collected stool and toilet paper samples, and dietary information from Inuit volunteers living in Resolute Bay (Nunavut, Canada), and compared them to individuals of European descent living in Montréal (Québec, Canada) consuming a typical Western diet. We sequenced the V4 region of the 16S rRNA gene to characterize microbiome diversity and composition, and compared samples collected with toilet paper or from stool. Our results show that these sampling methods provide similar, but non-identical portraits of the microbiome. Based on toilet paper samples, we found that much of the variation in microbiome community composition could be explained by individual identity (45-61% of variation explained, depending on the beta diversity metric used), with small but significant variation (3-5%) explained by sex or geography (Nunavut or Montréal). In contrast with our previous study at one time point, sampling over the course of a year revealed that diet explains 11% of variation in community composition across all participants, and 17% of variation specifically among Nunavut participants. However, we observed no clear seasonal shifts in the microbiomes of participants from either Nunavut or Montréal. Within-individual microbial diversity fluctuated more over time in Nunavut than in Montréal, consistent with a more variable and highly individualized diet in Nunavut.

CONCLUSIONS

Together, these results shows that the traditional Inuit diet and lifestyle has an impact on the composition, diversity and stability of the Inuit gut microbiome, even if the seasonality of the diet is less pronounced than expected, perhaps due to an increasingly westernized diet.

Total Evidence, Average Consensus and Matrix Representation with Parsimony: What a Difference Distances Make

Matrix representation with parsimony (MRP) can be used to combine trees in the supertree or the consensus settings. However, despite its popularity, it is still unclear whether MRP is really a consensus method or whether it behaves more like the total evidence approach. Previous simulations have shown that it approximates total evidence trees, whereas other studies have depicted similarities with average consensus trees. In this paper, we assess the hypothesis that MRP is equally related to both approaches. We conducted a simulation study to evaluate the accuracy of total evidence with that or various consensus methods, including MRP. Our results show that the total evidence trees are not significantly more accurate than average consensus trees that accounts for branch lengths, but that both perform better than MRP trees in the consensus setting. The accuracy rate of all methods was similarly affected by the number of taxa, the number of partitions, and the heterogeneity of the data.

A New Effective Method for Estimating Missing Values in the Sequence Data Prior to Phylogenetic Analysis

In this article we address the problem of phylogenetic inference from nucleic acid data containing missing bases. We introduce a new effective approach, called “Probabilistic estimation of missing values” (PEMV), allowing one to estimate unknown nucleotides prior to computing the evolutionary distances between them. We show that the new method improves the accuracy of phylogenetic inference compared to the existing methods “Ignoring Missing Sites” (IMS), “Proportional Distribution of Missing and Ambiguous Bases” (PDMAB) included in the PAUP software [ 26 ]. The proposed strategy for estimating missing nucleotides is based on probabilistic formulae developed in the framework of the Jukes-Cantor [ 10 ] and Kimura 2-parameter [ 11 ] models. The relative performances of the new method were assessed through simulations carried out with the SeqGen program [ 20 ], for data generation, and the BioNJ method [ 7 ], for inferring phylogenies. We also compared the new method to the DNAML program [ 5 ] and “Matrix Representation using Parsimony” (MRP) [ 13 , 19 ] considering an example of 66 eutherian mammals originally analyzed in [ 17 ].

Enrichment of beneficial bacteria in the skin microbiota of bats persisting with white-nose syndrome

BACKGROUND

Infectious diseases of wildlife are increasing worldwide with implications for conservation and human public health. The microbiota (i.e. microbial community living on or in a host) could influence wildlife disease resistance or tolerance. White-nose syndrome (WNS), caused by the fungus Pseudogymnoascus destructans (Pd), has killed millions of hibernating North American bats since 2007. We characterized the skin microbiota of naïve, pre-WNS little brown bats (Myotis lucifugus) from three WNS-negative hibernation sites and persisting, previously exposed bats from three WNS-positive sites to test the hypothesis that the skin microbiota of bats shifts following WNS invasion.

RESULTS

Using high-throughput 16S rRNA gene sequencing on 66 bats and 11 environmental samples, we found that hibernation site strongly influenced the composition and diversity of the skin microbiota. Bats from WNS-positive and WNS-negative sites differed in alpha and beta diversity, as well as in microbiota composition. Alpha diversity was reduced in persisting, WNS-positive bats, and the microbiota profile was enriched with particular taxa such Janthinobacterium, Micrococcaceae, Pseudomonas, Ralstonia, and Rhodococcus. Some of these taxa are recognized for their antifungal activity, and specific strains of Rhodococcus and Pseudomonas are known to inhibit Pd growth. Composition of the microbial community in the hibernaculum environment and the community on bat skin was superficially similar but differed in relative abundance of some bacterial taxa.

CONCLUSIONS

Our results are consistent with the hypothesis that Pd invasion leads to a shift in the skin microbiota of surviving bats and suggest the possibility that the microbiota plays a protective role for bats facing WNS. The detection of what appears to be enrichment of beneficial bacteria in the skin microbiota of persisting bats is a promising discovery for species re-establishment. Our findings highlight not only the potential value of management actions that might encourage transmission, growth, and establishment of beneficial bacteria on bats, and within hibernacula, but also the potential risks of such management actions.

MODELING BRAIN EVOLUTION FROM BEHAVIOR: A PERMUTATIONAL REGRESSION APPROACH

This paper has two complementary purposes: first, to present a method to perform multiple regression on distance matrices, with permutation testing appropriate for path-length matrices representing evolutionary trees, and then, to apply this method to study the joint evolution of brain, behavior and other characteristics in marsupials. To understand the computation method, consider that the dependent matrix is unfolded as a vector y; similarly, consider X to be a table containing the independent matrices, also unfolded as vectors. A multiple regression is computed to express y as a function of X. The parameters of this regression (R² and partial regression coefficients) are tested by permutations, as follows. When the dependent matrix variable y represents a simple distance or similarity matrix, permutations are performed in the same manner as the Mantel permutational test. When it is an ultrametric matrix representing a dendrogram, we use the double-permutation method (Lapointe and Legendre 1990, 1991). When it is a path-length matrix representing an additive tree (cladogram), we use the triple-permutation method (Lapointe and Legendre 1992). The independent matrix variables in X are kept fixed with respect to one another during the permutations. Selection of predictors can be accomplished by forward selection, backward elimination, or a stepwise procedure. A phylogenetic tree, derived from marsupial brain morphology data (28 species), is compared to trees depicting the evolution of diet, sociability, locomotion, and habitat in these animals, as well as their taxonomy and geographical relationships. A model is derived in which brain evolution can be predicted from taxonomy, diet, sociability and locomotion (R² = 0.75). A new tree, derived from the "predicted" data, shows a lot of similarity to the brain evolution tree. The meaning of the taxonomy, diet, sociability, and locomotion predictors are discussed and conclusions are drawn about the evolution of brain and behavior in marsupials.

Environment and host species shape the skin microbiome of captive neotropical bats

BACKGROUND

A wide range of microorganisms inhabit animal skin. This microbial community (microbiome) plays an important role in host defense against pathogens and disease. Bats (Chiroptera: Mammalia) are an ecologically and evolutionarily diversified group with a relatively unexplored skin microbiome. The bat skin microbiome could play a role in disease resistance, for example, to white nose syndrome (WNS), an infection which has been devastating North American bat populations. However, fundamental knowledge of the bat skin microbiome is needed before understanding its role in health and disease resistance. Captive neotropical frugivorous bats Artibeus jamaicensis and Carollia perspicillataprovide a simple controlled system in which to characterize the factors shaping the bat microbiome. Here, we aimed to determine the relative importance of habitat and host species on the bat skin microbiome.

METHODS

We performed high-throughput 16S rRNA gene sequencing of the skin microbiome of two different bat species living in captivity in two different habitats. In the first habitat, A. jamaicensis and C. perspicillata lived together, while the second habitat contained only A. jamaicensis.

RESULTS

We found that both habitat and host species shape the composition and diversity of the skin microbiome, with habitat having the strongest influence. Cohabitating A. jamaicensis and C. perspicillata shared more similar skin microbiomes than members of the same species (A. jamaicensis) across two habitats.

DISCUSSION

These results suggest that in captivity, the skin microbial community is homogenised by the shared environments and individual proximities of bats living together in the same habitat, at the expense of the innate host species factors. The predominant influence of habitat suggests that environmental microorganisms or pathogens might colonize bat skin. We also propose that bat populations could differ in pathogen susceptibility depending on their immediate environment and habitat.

Using hybridization networks to retrace the evolution of Indo-European languages

Background

Curious parallels between the processes of species and language evolution have been observed by many researchers. Retracing the evolution of Indo-European (IE) languages remains one of the most intriguing intellectual challenges in historical linguistics. Most of the IE language studies use the traditional phylogenetic tree model to represent the evolution of natural languages, thus not taking into account reticulate evolutionary events, such as language hybridization and word borrowing which can be associated with species hybridization and horizontal gene transfer, respectively. More recently, implicit evolutionary networks, such as split graphs and minimal lateral networks, have been used to account for reticulate evolution in linguistics.

Results

Striking parallels existing between the evolution of species and natural languages allowed us to apply three computational biology methods for reconstruction of phylogenetic networks to model the evolution of IE languages. We show how the transfer of methods between the two disciplines can be achieved, making necessary methodological adaptations. Considering basic vocabulary data from the well-known Dyen’s lexical database, which contains word forms in 84 IE languages for the meanings of a 200-meaning Swadesh list, we adapt a recently developed computational biology algorithm for building explicit hybridization networks to study the evolution of IE languages and compare our findings to the results provided by the split graph and galled network methods.

Conclusion

We conclude that explicit phylogenetic networks can be successfully used to identify donors and recipients of lexical material as well as the degree of influence of each donor language on the corresponding recipient languages. We show that our algorithm is well suited to detect reticulate relationships among languages, and present some historical and linguistic justification for the results obtained. Our findings could be further refined if relevant syntactic, phonological and morphological data could be analyzed along with the available lexical data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0745-6) contains supplementary material, which is available to authorized users.

Self-tracking the microbiome: where do we go from here?

The quantified self community brings together enthusiasts who are using technological devices to monitor their health and social media to share their personal data with others online. In light of the growing popularity of this movement, self-trackers are challenging the health-care system by raising important questions about data ownership and risk-taking. As we enter a new era of consumer genomics, a significant number of quantified self (QS) individuals are now interested in the monitoring of their microbiome and performing personal interventions. In this paper, we discuss the scientific validity of experiments involving serial observations of a single individual as opposed to randomized clinical trials. We look at self-tracking from an ethical standpoint by questioning the risks and assessing the potential benefits for personalized medicine in general and for microbiome research in particular.

Being human is a gut feeling

Some metagenomic studies have suggested that less than 10% of the cells that comprise our bodies are Homo sapiens cells. The remaining 90% are bacterial cells. The description of this so-called human microbiome is of great interest and importance for several reasons. For one, it helps us redefine what a biological individual is. We suggest that a human individual is now best described as a super-individual in which a large number of different species (including Homo sapiens) coexist. New concepts of biological individuality must extend beyond the traditional limitations of our own skin to include our resident microbes. Besides its important contributions to science, microbiome research raises philosophical questions that strike close to home. What is left of Homo sapiens? If most of our cells are not Homo sapiens cells, what does it mean to be an individual human being? In this paper, we argue that the biological individual is determined by the amount of functional integration among its constitutive parts, a definition that applies perfectly to Homo sapiens and its microbiome.

Genetic structure of the white-footed mouse in the context of the emergence of Lyme disease in southern Québec

The white-footed mouse (Peromyscus leucopus) has expanded its northern limit into southern Québec over the last few decades. P. leucopus is a great disperser and colonizer and is of particular interest because it is considered a primary reservoir for the spirochete bacterium that causes Lyme disease. There is no current information on the gene flow between mouse populations on the mountains and forest fragments found scattered throughout the Montérégie region in southern Québec, and whether various landscape barriers have an effect on their dispersal. We conducted a population genetics analysis on eleven P. leucopus populations using eleven microsatellite markers and showed that isolation by distance was weak, yet barriers were effective. The agricultural matrix had the least effect on gene flow, whereas highways and main rivers were effective barriers. The abundance of ticks collected from mice varied within the study area. Both ticks and mice were screened for the presence of the spirochete bacterium Borrelia burgdorferi, and we predicted areas of greater risk for Lyme disease. Merging our results with ongoing Lyme disease surveillance programs will help determine the future threat of this disease in Québec, and will contribute toward disease prevention and management strategies throughout fragmented landscapes in southern Canada.

Canid hybridization: contemporary evolution in human-modified landscapes

Contemporary evolution through human-induced hybridization occurs throughout the taxonomic range. Formerly allopatric species appear especially susceptible to hybridization. Consequently, hybridization is expected to be more common in regions with recent sympatry owing to human activity than in areas of historical range overlap. Coyotes (Canis latrans) and gray wolves (C. lupus) are historically sympatric in western North America. Following European settlement gray wolf range contracted, whereas coyote range expanded to include eastern North America. Furthermore, wolves with New World (NW) mitochondrial DNA (mtDNA) haplotypes now extend from Manitoba to Québec in Canada and hybridize with gray wolves and coyotes. Using mtDNA and 12 microsatellite markers, we evaluated levels of wolf-coyote hybridization in regions where coyotes were present (the Canadian Prairies, n = 109 samples) and absent historically (Québec, n = 154). Wolves with NW mtDNA extended from central Saskatchewan (51°N, 69°W) to northeastern Québec (54°N, 108°W). On the Prairies, 6.3% of coyotes and 9.2% of wolves had genetic profiles suggesting wolf-coyote hybridization. In contrast, 12.6% of coyotes and 37.4% of wolves in Québec had profiles indicating hybrid origin. Wolves with NW and Old World (C. lupus) mtDNA appear to form integrated populations in both regions. Our results suggest that hybridization is more frequent in historically allopatric populations. Range shifts, now expected across taxa following climate change and other human influence on the environment, might therefore promote contemporary evolution by hybridization.

Of woods and webs: possible alternatives to the tree of life for studying genomic fluidity in E. coli

Background

We introduce several forest-based and network-based methods for exploring microbial evolution, and apply them to the study of thousands of genes from 30 strains of E. coli. This case study illustrates how additional analyses could offer fast heuristic alternatives to standard tree of life (TOL) approaches.

Results

We use gene networks to identify genes with atypical modes of evolution, and genome networks to characterize the evolution of genetic partnerships between E. coli and mobile genetic elements. We develop a novel polychromatic quartet method to capture patterns of recombination within E. coli, to update the clanistic toolkit, and to search for the impact of lateral gene transfer and of pathogenicity on gene evolution in two large forests of trees bearing E. coli. We unravel high rates of lateral gene transfer involving E. coli (about 40% of the trees under study), and show that both core genes and shell genes of E. coli are affected by non-tree-like evolutionary processes. We show that pathogenic lifestyle impacted the structure of 30% of the gene trees, and that pathogenic strains are more likely to transfer genes with one another than with non-pathogenic strains. In addition, we propose five groups of genes as candidate mobile modules of pathogenicity. We also present strong evidence for recent lateral gene transfer between E. coli and mobile genetic elements.

Conclusions

Depending on which evolutionary questions biologists want to address (i.e. the identification of modules, genetic partnerships, recombination, lateral gene transfer, or genes with atypical evolutionary modes, etc.), forest-based and network-based methods are preferable to the reconstruction of a single tree, because they provide insights and produce hypotheses about the dynamics of genome evolution, rather than the relative branching order of species and lineages. Such a methodological pluralism - the use of woods and webs - is to be encouraged to analyse the evolutionary processes at play in microbial evolution.

This manuscript was reviewed by: Ford Doolittle, Tal Pupko, Richard Burian, James McInerney, Didier Raoult, and Yan Boucher

Evaluating phylogenetic congruence in the post-genomic era

Congruence is a broadly applied notion in evolutionary biology used to justify multigene phylogeny or phylogenomics, as well as in studies of coevolution, lateral gene transfer, and as evidence for common descent. Existing methods for identifying incongruence or heterogeneity using character data were designed for data sets that are both small and expected to be rarely incongruent. At the same time, methods that assess incongruence using comparison of trees test a null hypothesis of uncorrelated tree structures, which may be inappropriate for phylogenomic studies. As such, they are ill-suited for the growing number of available genome sequences, most of which are from prokaryotes and viruses, either for phylogenomic analysis or for studies of the evolutionary forces and events that have shaped these genomes. Specifically, many existing methods scale poorly with large numbers of genes, cannot accommodate high levels of incongruence, and do not adequately model patterns of missing taxa for different markers. We propose the development of novel incongruence assessment methods suitable for the analysis of the molecular evolution of the vast majority of life and support the investigation of homogeneity of evolutionary process in cases where markers do not share identical tree structures.

The performance of the Congruence Among Distance Matrices (CADM) test in phylogenetic analysis

BACKGROUND

CADM is a statistical test used to estimate the level of Congruence Among Distance Matrices. It has been shown in previous studies to have a correct rate of type I error and good power when applied to dissimilarity matrices and to ultrametric distance matrices. Contrary to most other tests of incongruence used in phylogenetic analysis, the null hypothesis of the CADM test assumes complete incongruence of the phylogenetic trees instead of congruence. In this study, we performed computer simulations to assess the type I error rate and power of the test. It was applied to additive distance matrices representing phylogenies and to genetic distance matrices obtained from nucleotide sequences of different lengths that were simulated on randomly generated trees of varying sizes, and under different evolutionary conditions.

RESULTS

Our results showed that the test has an accurate type I error rate and good power. As expected, power increased with the number of objects (i.e., taxa), the number of partially or completely congruent matrices and the level of congruence among distance matrices.

CONCLUSIONS

Based on our results, we suggest that CADM is an excellent candidate to test for congruence and, when present, to estimate its level in phylogenomic studies where numerous genes are analysed simultaneously.

Harvesting evolutionary signals in a forest of prokaryotic gene trees

Phylogenomic studies produce increasingly large phylogenetic forests of trees with patchy taxonomical sampling. Typically, prokaryotic data generate thousands of gene trees of all sizes that are difficult, if not impossible, to root. Their topologies do not match the genealogy of lineages, as they are influenced not only by duplication, losses, and vertical descent but also by lateral gene transfer (LGT) and recombination. Because this complexity in part reflects the diversity of evolutionary processes, the study of phylogenetic forests is thus a great opportunity to improve our understanding of prokaryotic evolution. Here, we show how the rich evolutionary content of such novel phylogenetic objects can be exploited through the development of new approaches designed specifically for extracting the multiple evolutionary signals present in the forest of life, that is, by slicing up trees into remarkable bits and pieces: clans, slices, and clips. We harvested a forest of 6,901 unrooted gene trees comprising up to 100 prokaryotic genomes (41 archaea and 59 bacteria) to search for evolutionary events that a species tree would not account for. We identified 1) trees and partitions of trees that reflected the lifestyle of organisms rather than their taxonomy, 2) candidate lifestyle-specific genetic modules, used by distinct unrelated organisms to adapt to the same environment, 3) gene families, nonrandomly distributed in the functional space, that were frequently exchanged between archaea and bacteria, sometimes without major changes in their sequences. Finally, 4) we reconstructed polarized networks of genetic partnerships between archaea and bacteria to describe some of the rules affecting LGT between these two Domains.

Weighted bootstrapping: a correction method for assessing the robustness of phylogenetic trees

Background

Non-parametric bootstrapping is a widely-used statistical procedure for assessing confidence of model parameters based on the empirical distribution of the observed data [1] and, as such, it has become a common method for assessing tree confidence in phylogenetics [2]. Traditional non-parametric bootstrapping does not weigh each tree inferred from resampled (i.e., pseudo-replicated) sequences. Hence, the quality of these trees is not taken into account when computing bootstrap scores associated with the clades of the original phylogeny. As a consequence, traditionally, the trees with different bootstrap support or those providing a different fit to the corresponding pseudo-replicated sequences (the fit quality can be expressed through the LS, ML or parsimony score) contribute in the same way to the computation of the bootstrap support of the original phylogeny.

Results

In this article, we discuss the idea of applying weighted bootstrapping to phylogenetic reconstruction by weighting each phylogeny inferred from resampled sequences. Tree weights can be based either on the least-squares (LS) tree estimate or on the average secondary bootstrap score (SBS) associated with each resampled tree. Secondary bootstrapping consists of the estimation of bootstrap scores of the trees inferred from resampled data. The LS and SBS-based bootstrapping procedures were designed to take into account the quality of each "pseudo-replicated" phylogeny in the final tree estimation. A simulation study was carried out to evaluate the performances of the five weighting strategies which are as follows: LS and SBS-based bootstrapping, LS and SBS-based bootstrapping with data normalization and the traditional unweighted bootstrapping.

Conclusions

The simulations conducted with two real data sets and the five weighting strategies suggest that the SBS-based bootstrapping with the data normalization usually exhibits larger bootstrap scores and a higher robustness compared to the four other competing strategies, including the traditional bootstrapping. The high robustness of the normalized SBS could be particularly useful in situations where observed sequences have been affected by noise or have undergone massive insertion or deletion events. The results provided by the four other strategies were very similar regardless the noise level, thus also demonstrating the stability of the traditional bootstrapping method.

An application of supertree methods to Mammalian mitogenomic sequences

TWO DIFFERENT APPROACHES CAN BE USED IN PHYLOGENOMICS: combined or separate analysis. In the first approach, different datasets are combined in a concatenated supermatrix. In the second, datasets are analyzed separately and the phylogenetic trees are then combined in a supertree. The supertree method is an interesting alternative to avoid missing data, since datasets that are analyzed separately do not need to represent identical taxa. However, the supertree approach and the corresponding consensus methods have been highly criticized for not providing valid phylogenetic hypotheses. In this study, congruence of trees estimated by consensus and supertree approaches were compared to model trees obtained from a combined analysis of complete mitochondrial sequences of 102 species representing 93 mammal families. The consensus methods produced poorly resolved consensus trees and did not perform well, except for the majority rule consensus with compatible groupings. The weighted supertree and matrix representation with parsimony methods performed equally well and were highly congruent with the model trees. The most similar supertree method was the least congruent with the model trees. We conclude that some of the methods tested are worth considering in a phylogenomic context.

Phylogeny of the rock wallabies, Petrogale (Marsupialia: Macropodidae). Part II: Detection of hybridisation among macropodines

Phylogenetic relationships among rock-wallabies, Petrogale (Marsupialia: Macropodidae), have proven difficult to resolve. Given the documented interspecific hybridisation in the wild and the ease with which hybrids can be bred in captivity, introgression and hybrid speciation are likely explanations for these difficulties. In this paper, an attempt is made at using a phylogenetic approach to identify Petrogale hybrids of known origin. The Hybrid Detection Criterion (HDC) test is applied to DNA–DNA hybridisation data for 15 full species, two natural yard-bred hybrids, and two artificial hybrids from the same pairs of parental species. While the yard-bred hybrids elude detection with this technique, the artificial hybrids, consisting of equimolar mixture of parental extracts, are easily identified. Moreover, splitsgraphs constructed from five pairs of natural and artificial hybrids, including those evaluated with HDC, and their parents show that, in all cases but one, these two kinds of hybrids do not group together. Because the HDC assumes an intermediate phylogenetic position of the hybrid between its postulated parents, it is likely that unequal crossing-over, or another recombination event, affects the results of the test. These conclusions cast some doubt on the possibility of accurately detecting Petrogale hybrids with a phylogenetic approach.

Prokaryotic evolution and the tree of life are two different things

BACKGROUND

The concept of a tree of life is prevalent in the evolutionary literature. It stems from attempting to obtain a grand unified natural system that reflects a recurrent process of species and lineage splittings for all forms of life. Traditionally, the discipline of systematics operates in a similar hierarchy of bifurcating (sometimes multifurcating) categories. The assumption of a universal tree of life hinges upon the process of evolution being tree-like throughout all forms of life and all of biological time. In multicellular eukaryotes, the molecular mechanisms and species-level population genetics of variation do indeed mainly cause a tree-like structure over time. In prokaryotes, they do not. Prokaryotic evolution and the tree of life are two different things, and we need to treat them as such, rather than extrapolating from macroscopic life to prokaryotes. In the following we will consider this circumstance from philosophical, scientific, and epistemological perspectives, surmising that phylogeny opted for a single model as a holdover from the Modern Synthesis of evolution.

RESULTS

It was far easier to envision and defend the concept of a universal tree of life before we had data from genomes. But the belief that prokaryotes are related by such a tree has now become stronger than the data to support it. The monistic concept of a single universal tree of life appears, in the face of genome data, increasingly obsolete. This traditional model to describe evolution is no longer the most scientifically productive position to hold, because of the plurality of evolutionary patterns and mechanisms involved. Forcing a single bifurcating scheme onto prokaryotic evolution disregards the non-tree-like nature of natural variation among prokaryotes and accounts for only a minority of observations from genomes.

CONCLUSION

Prokaryotic evolution and the tree of life are two different things. Hence we will briefly set out alternative models to the tree of life to study their evolution. Ultimately, the plurality of evolutionary patterns and mechanisms involved, such as the discontinuity of the process of evolution across the prokaryote-eukaryote divide, summons forth a pluralistic approach to studying evolution.

REVIEWERS

This article was reviewed by Ford Doolittle, John Logsdon and Nicolas Galtier.

Total evidence, average consensus and matrix representation with parsimony: what a difference distances make

Matrix representation with parsimony (MRP) can be used to combine trees in the supertree or the consensus settings. However, despite its popularity, it is still unclear whether MRP is really a consensus method or whether it behaves more like the total evidence approach. Previous simulations have shown that it approximates total evidence trees, whereas other studies have depicted similarities with average consensus trees. In this paper, we assess the hypothesis that MRP is equally related to both approaches. We conducted a simulation study to evaluate the accuracy of total evidence with that or various consensus methods, including MRP. Our results show that the total evidence trees are not significantly more accurate than average consensus trees that accounts for branch lengths, but that both perform better than MRP trees in the consensus setting. The accuracy rate of all methods was similarly affected by the number of taxa, the number of partitions, and the heterogeneity of the data.

A new effective method for estimating missing values in the sequence data prior to phylogenetic analysis

In this article we address the problem of phylogenetic inference from nucleic acid data containing missing bases. We introduce a new effective approach, called "Probabilistic estimation of missing values" (PEMV), allowing one to estimate unknown nucleotides prior to computing the evolutionary distances between them. We show that the new method improves the accuracy of phylogenetic inference compared to the existing methods "Ignoring Missing Sites" (IMS), "Proportional Distribution of Missing and Ambiguous Bases" (PDMAB) included in the PAUP software [26]. The proposed strategy for estimating missing nucleotides is based on probabilistic formulae developed in the framework of the Jukes-Cantor [10] and Kimura 2-parameter [11] models. The relative performances of the new method were assessed through simulations carried out with the SeqGen program [20], for data generation, and the Bio NJ method [7], for inferring phylogenies. We also compared the new method to the DNAML program [5] and "Matrix Representation using Parsimony" (MRP) [13], [19] considering an example of 66 eutherian mammals originally analyzed in [17].

Conservation genetics of the wood turtle (Glyptemys insculpta) in Quebec, Canada

Throughout its range, the long-term persistence of the wood turtle (Glyptemys insculpta Le Conte, 1830) is jeopardized by habitat perturbations and commercial collection. The main objective of this study was to acquire knowledge on the genetic structure of wood turtle populations within Quebec, where the species reaches the northernmost limit of its range, to identify proper conservation units. The six known populations in Quebec were genetically characterized using five microsatellite loci. Genetic analyses revealed high variability among all populations (HOranging from 0.561 to 0.886), suggesting that past population declines have not yet affected their genetic diversity. Since populations are located on both shores of the St. Lawrence River, two colonization hypotheses were tested to determine whether the River acted as a dispersal route or as a barrier. Results of AMOVA tests, phylogenetic analyses, and assignment tests demonstrated the clear distinction between populations from both shores, and three genetically different units were identified: the two populations on the north shore and a single homogenous group on the south shore.

Peromyscuspopulations and theirCuterebraparasites display congruent phylogeographical structure

The relationships between populations of the Deer Mouse (Peromyscus maniculatus) and the White-footed Mouse (P. leucopus) and their respectiveCuterebraparasites were examined. Population genetic structure of hosts and parasites was inferred using cytochrome oxidase mitochondrial sequences of specimens from 7 populations. Genetic analyses revealed that isolation-by-distance applies forP. maniculatusand its associated parasite (C. grisea). A significant correlation was also observed between the genetic distances of these host and parasite species. Furthermore, populations ofP. maniculatusandC. griseafrom the North and South shores of the St Lawrence River were found to be significantly different. This structure may be explained by the St Lawrence River being a dispersal barrier for both species. A robust analysis of the other species pair (P. leucopusandC. fontinella) could not be performed because of limited sample sizes.

A new method to discriminate the deer mouse (Peromyscus maniculatus) from the white-footed mouse (Peromyscus leucopus) using species-specific primers in multiplex PCR

Morphometric measurements or biochemical methods are often required to differentiate deer mice, Peromyscus maniculatus (Wagner, 1845), from white-footed mice, Peromyscus leucopus (Rafinesque, 1818), particularly when they are found in sympatry. However, these approaches cannot easily be applied to juveniles, or to degraded or ancient museum specimens. In this paper, we propose a rapid and non-invasive molecular approach to discriminate these cryptic species from one another. This technique relies on species-specific primers designed in the COIII mitochondrial region to amplify fragments of different lengths in each species. The method developed proved useful for the identification of ethanol-preserved, frozen, degraded, or dry museum specimens.

A weighted least-squares approach for inferring phylogenies from incomplete distance matrices

MOTIVATION

The problem of phylogenetic inference from datasets including incomplete or uncertain entries is among the most relevant issues in systematic biology. In this paper, we propose a new method for reconstructing phylogenetic trees from partial distance matrices. The new method combines the usage of the four-point condition and the ultrametric inequality with a weighted least-squares approximation to solve the problem of missing entries. It can be applied to infer phylogenies from evolutionary data including some missing or uncertain information, for instance, when observed nucleotide or protein sequences contain gaps or missing entries.

RESULTS

In a number of simulations involving incomplete datasets, the proposed method outperformed the well-known Ultrametric and Additive procedures. Generally, the new method also outperformed all the other competing approaches including Triangle and Fitch which is the most popular least-squares method for reconstructing phylogenies. We illustrate the usefulness of the introduced method by analyzing two well-known phylogenies derived from complete mammalian mtDNA sequences. Some interesting theoretical results concerning the NP-hardness of the ordinary and weighted least-squares fitting of a phylogenetic tree to a partial distance matrix are also established.

AVAILABILITY

The T-Rex package including this method is freely available for download at http://www.info.uqam.ca/~makarenv/trex.html