Inferences from tip-calibrated phylogenies: a review and a practical guide

Molecular dating of phylogenetic trees is a growing discipline using sequence data to co-estimate the timing of evolutionary events and rates of molecular evolution. All molecular-dating methods require converting genetic divergence between sequences into absolute time. Historically, this could only be achieved by associating externally derived dates obtained from fossil or biogeographical evidence to internal nodes of the tree. In some cases, notably for fast-evolving genomes such as viruses and some bacteria, the time span over which samples were collected may cover a significant proportion of the time since they last shared a common ancestor. This situation allows phylogenetic trees to be calibrated by associating sampling dates directly to the sequences representing the tips (terminal nodes) of the tree. The increasing availability of genomic data from ancient DNA extends the applicability of such tip-based calibration to a variety of taxa including humans, extinct megafauna and various microorganisms which typically have a scarce fossil record. The development of statistical models accounting for heterogeneity in different aspects of the evolutionary process while accommodating very large data sets (e.g. whole genomes) has allowed using tip-dating methods to reach inferences on divergence times, substitution rates, past demography or the age of specific mutations on a variety of spatiotemporal scales. In this review, we summarize the current state of the art of tip dating, discuss some recent applications, highlight common pitfalls and provide a 'how to' guide to thoroughly perform such analyses

Four decades of transmission of a multidrug-resistant Mycobacterium tuberculosis outbreak strain

The rise of drug-resistant strains is a major challenge to containing the tuberculosis (TB) pandemic. Yet, little is known about the extent of resistance in early years of chemotherapy and when transmission of resistant strains on a larger scale became a major public health issue. Here we reconstruct the timeline of the acquisition of antimicrobial resistance during a major ongoing outbreak of multidrug-resistant TB in Argentina. We estimate that the progenitor of the outbreak strain acquired resistance to isoniazid, streptomycin and rifampicin by around 1973, indicating continuous circulation of a multidrug-resistant TB strain for four decades. By around 1979 the strain had acquired additional resistance to three more drugs. Our results indicate that Mycobacterium tuberculosis (Mtb) with extensive resistance profiles circulated 15 years before the outbreak was detected, and about one decade before the earliest documented transmission of Mtb strains with such extensive resistance profiles globally.Fil: Eldholm, Vegard. Norwegian Institute of Public Health; NoruegaFil: Monteserin, Johana. Dirección Nacional de Institutos de Investigación. Administración Nacional de Laboratorios e Institutos de Salud. Instituto Nacional de Enfermedades Infecciosas; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Rieux, Adrien. Colegio Universitario de Londres; Reino UnidoFil: Lopez, Beatriz. Dirección Nacional de Institutos de Investigación. Administración Nacional de Laboratorios e Institutos de Salud. Instituto Nacional de Enfermedades Infecciosas; ArgentinaFil: Sobkowiak, Benjamin. Colegio Universitario de Londres; Reino UnidoFil: Ritacco, Gloria Viviana. Dirección Nacional de Institutos de Investigación. Administración Nacional de Laboratorios e Institutos de Salud. Instituto Nacional de Enfermedades Infecciosas; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Balloux, Francois. Colegio Universitario de Londres; Reino Unid

Antimicrobial Resistance in Mycobacterium tuberculosis: The Odd One Out

Antimicrobial resistance (AMR) threats are typically represented by bacteria capable of extensive horizontal gene transfer (HGT). One clear exception is Mycobacterium tuberculosis (Mtb). It is an obligate human pathogen with limited genetic diversity and a low mutation rate which lacks any evidence for HGT. Such features should, in principle, reduce its ability to rapidly evolve AMR. We identify key features in its biology and epidemiology that allow it to overcome its low adaptive potential. We focus in particular on its innate resistance to drugs, its unusual life cycle, including an often extensive latent phase, and its ability to shelter from exposure to antimicrobial drugs within cavities it induces in the lungs

Spatiotemporal dynamics in the early stages of the 2009 A/H1N1 influenza pandemic.

Epidemiology and public health planning will increasingly rely on the analysis of genetic sequence data. The ongoing influenza A/H1N1 pandemic may represent a tipping point in this trend, with A/H1N1 being the first human pathogen routinely genotyped from the beginning of its spread. To take full advantage of this genetic information, we introduce a novel method to reconstruct the spatiotemporal dynamics of outbreaks from sequence data. The approach is based on a new paradigm were ancestries are inferred directly rather than through the reconstruction of most recent common ancestors (MRCAs) as in phylogenetics. Using 279 A/H1N1 hemagglutinin (HA) sequences, we confirm the emergence of the 2009 flu pandemic in Mexico. The virus initially spread to the US, and then to the rest of the world with both Mexico and the US acting as the main sources. While compatible with current epidemiological understanding of the 2009 H1N1 pandemic, our results provide a much finer picture of the spatiotemporal dynamics. The results also highlight how much additional epidemiological information can be gathered from genetic monitoring of a disease outbreak

The population genetics of clonal and partially clonal diploids.

The consequences of variable rates of clonal reproduction on the population genetics of neutral markers are explored in diploid organisms within a subdivided population (island model). We use both analytical and stochastic simulation approaches. High rates of clonal reproduction will positively affect heterozygosity. As a consequence, nearly twice as many alleles per locus can be maintained and population differentiation estimated as F(ST) value is strongly decreased in purely clonal populations as compared to purely sexual ones. With increasing clonal reproduction, effective population size first slowly increases and then points toward extreme values when the reproductive system tends toward strict clonality. This reflects the fact that polymorphism is protected within individuals due to fixed heterozygosity. Contrarily, genotypic diversity smoothly decreases with increasing rates of clonal reproduction. Asexual populations thus maintain higher genetic diversity at each single locus but a lower number of different genotypes. Mixed clonal/sexual reproduction is nearly indistinguishable from strict sexual reproduction as long as the proportion of clonal reproduction is not strongly predominant for all quantities investigated, except for genotypic diversities (both at individual loci and over multiple loci)

Do riverine barriers, history or introgression shape the genetic structuring of a common shrew (Sorex araneus) population?

The common shrew (Sorex araneus) is subdivided into numerous chromosome races. The Valais and Cordon chromosome races meet and hybridize at a mountain river in Les Houches (French Alps). Significant genetic structuring was recently reported among populations found on the Valais side of this hybrid zone. In this paper, a phylogenetic analysis and partial Mantel tests are used to investigate the patterns and causes of this structuring. A total of 185 shrews were trapped at 12 localities. All individuals were typed for nine microsatellite loci. Although several mountain rivers are found in the study area, riverine barriers do not have a significant influence on gene flow. Partial Mantel tests show that our result is caused by the influence of the hybrid zone with the Cordon race. The geographical patterns of this structuring are discussed in the context of the contact zone, which appears to extend up to a group of two rivers. The glacier they originate from is known to have cut the Arve valley as recently as 1818. The recent history of this glacier, its moraine and possibly rivers, may therefore be linked to the history of this hybrid zone

Islands of linkage in an ocean of pervasive recombination reveals two-speed evolution of human cytomegalovirus genomes

Human cytomegalovirus (HCMV) infects most of the population worldwide, persisting throughout the host's life in a latent state with periodic episodes of reactivation. While typically asymptomatic, HCMV can cause fatal disease among congenitally infected infants and immunocompromised patients. These clinical issues are compounded by the emergence of antiviral resistance and the absence of an effective vaccine, the development of which is likely complicated by the numerous immune evasins encoded by HCMV to counter the host's adaptive immune responses, a feature that facilitates frequent super-infections. Understanding the evolutionary dynamics of HCMV is essential for the development of effective new drugs and vaccines. By comparing viral genomes from uncultivated or low-passaged clinical samples of diverse origins, we observe evidence of frequent homologous recombination events, both recent and ancient, and no structure of HCMV genetic diversity at the whole-genome scale. Analysis of individual gene-scale loci reveals a striking dichotomy: while most of the genome is highly conserved, recombines essentially freely and has evolved under purifying selection, 21 genes display extreme diversity, structured into distinct genotypes that do not recombine with each other. Most of these hyper-variable genes encode glycoproteins involved in cell entry or escape of host immunity. Evidence that half of them have diverged through episodes of intense positive selection suggests that rapid evolution of hyper-variable loci is likely driven by interactions with host immunity. It appears that this process is enabled by recombination unlinking hyper-variable loci from strongly constrained neighboring sites. It is conceivable that viral mechanisms facilitating super-infection have evolved to promote recombination between diverged genotypes, allowing the virus to continuously diversify at key loci to escape immune detection, while maintaining a genome optimally adapted to its asymptomatic infectious lifecycle

Population genetics of trypanosoma brucei rhodesiense: clonality and diversity within and between foci

African trypanosomes are unusual among pathogenic protozoa in that they can undergo their complete morphological life cycle in the tsetse fly vector with mating as a non-obligatory part of this development. Trypanosoma brucei rhodesiense, which infects humans and livestock in East and Southern Africa, has classically been described as a host-range variant of the non-human infective Trypanosoma brucei that occurs as stable clonal lineages. We have examined T. b. rhodesiense populations from East (Uganda) and Southern (Malawi) Africa using a panel of microsatellite markers, incorporating both spatial and temporal analyses. Our data demonstrate that Ugandan T. b. rhodesiense existed as clonal populations, with a small number of highly related genotypes and substantial linkage disequilibrium between pairs of loci. However, these populations were not stable as the dominant genotypes changed and the genetic diversity also reduced over time. Thus these populations do not conform to one of the criteria for strict clonality, namely stability of predominant genotypes over time, and our results show that, in a period in the mid 1990s, the previously predominant genotypes were not detected but were replaced by a novel clonal population with limited genetic relationship to the original population present between 1970 and 1990. In contrast, the Malawi T. b. rhodesiense population demonstrated significantly greater diversity and evidence for frequent genetic exchange. Therefore, the population genetics of T. b. rhodesiense is more complex than previously described. This has important implications for the spread of the single copy T. b. rhodesiense gene that allows human infectivity, and therefore the epidemiology of the human disease, as well as suggesting that these parasites represent an important organism to study the influence of optional recombination upon population genetic dynamics

Discovery of mating in the major African livestock pathogen Trypanosoma congolense

The protozoan parasite, Trypanosoma congolense, is one of the most economically important pathogens of livestock in Africa and, through its impact on cattle health and productivity, has a significant effect on human health and well being. Despite the importance of this parasite our knowledge of some of the fundamental biological processes is limited. For example, it is unknown whether mating takes place. In this paper we have taken a population genetics based approach to address this question. The availability of genome sequence of the parasite allowed us to identify polymorphic microsatellite markers, which were used to genotype T. congolense isolates from livestock in a discrete geographical area of The Gambia. The data showed a high level of diversity with a large number of distinct genotypes, but a deficit in heterozygotes. Further analysis identified cryptic genetic subdivision into four sub-populations. In one of these, parasite genotypic diversity could only be explained by the occurrence of frequent mating in T. congolense. These data are completely inconsistent with previous suggestions that the parasite expands asexually in the absence of mating. The discovery of mating in this species of trypanosome has significant consequences for the spread of critical traits, such as drug resistance, as well as for fundamental aspects of the biology and epidemiology of this neglected but economically important pathogen

Islands of linkage in an ocean of pervasive recombination reveals two-speed evolution of human cytomegalovirus genomes

Human cytomegalovirus (HCMV) infects most of the population worldwide, persisting throughout the host's life in a latent state with periodic episodes of reactivation. While typically asymptomatic, HCMV can cause fatal disease among congenitally infected infants and immunocompromised patients. These clinical issues are compounded by the emergence of antiviral resistance and the absence of an effective vaccine, the development of which is likely complicated by the numerous immune evasins encoded by HCMV to counter the host's adaptive immune responses, a feature that facilitates frequent super-infections. Understanding the evolutionary dynamics of HCMV is essential for the development of effective new drugs and vaccines. By comparing viral genomes from uncultivated or low-passaged clinical samples of diverse origins, we observe evidence of frequent homologous recombination events, both recent and ancient, and no structure of HCMV genetic diversity at the whole-genome scale. Analysis of individual gene-scale loci reveals a striking dichotomy: while most of the genome is highly conserved, recombines essentially freely and has evolved under purifying selection, 21 genes display extreme diversity, structured into distinct genotypes that do not recombine with each other. Most of these hyper-variable genes encode glycoproteins involved in cell entry or escape of host immunity. Evidence that half of them have diverged through episodes of intense positive selection suggests that rapid evolution of hyper-variable loci is likely driven by interactions with host immunity. It appears that this process is enabled by recombination unlinking hyper-variable loci from strongly constrained neighboring sites. It is conceivable that viral mechanisms facilitating super-infection have evolved to promote recombination between diverged genotypes, allowing the virus to continuously diversify at key loci to escape immune detection, while maintaining a genome optimally adapted to its asymptomatic infectious lifecycle