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

Long genetic and social isolation in Neanderthals before their extinction

Neanderthal genomes have been recovered from sites across Eurasia, painting an increasingly complex picture of their populations' structure that mostly indicates that late European Neanderthals belonged to a single metapopulation with no significant evidence of population structure. Here, we report the discovery of a late Neanderthal individual, nicknamed "Thorin," from Grotte Mandrin in Mediterranean France, and his genome. These dentognathic fossils, including a rare example of distomolars, are associated with a rich archeological record of Neanderthal final technological traditions in this region ∼50-42 thousand years ago. Thorin's genome reveals a relatively early divergence of ∼105 ka with other late Neanderthals. Thorin belonged to a population with a small group size that showed no genetic introgression with other known late European Neanderthals, revealing some 50 ka of genetic isolation of his lineage despite them living in neighboring regions. These results have important implications for resolving competing hypotheses about causes of the disappearance of the Neanderthals.

Fast Bayesian Inference of Phylogenies from Multiple Continuous Characters

Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remain challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous morphological characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.

Emergence and spread of Mycobacterium ulcerans at different geographic scales

The classical lineage of Mycobacterium ulcerans is the most prevalent clonal group associated with Buruli ulcer in humans. Its reservoir is strongly associated with the environment. We analyzed together 1,045 isolates collected from 13 countries on two continents to define the evolutionary history and population dynamics of this lineage. We confirm that this lineage spread over 7,000 years from Australia to Africa with the emergence of outbreaks in distinct waves in the 18th and 19th centuries. In sharp contrast with its global spread over the last century, transmission chains are now mostly local, with little or no dissemination between endemic areas. This study provides new insights into the phylogeography and population dynamics of M. ulcerans, highlighting the importance of comparative genomic analyses to improve our understanding of pathogen transmission.

IMPORTANCE

Mycobacterium ulcerans is an environmental mycobacterial pathogen that can cause Buruli ulcer, a severe cutaneous infection, mostly spread in Africa and Australia. We conducted a large genomic study of M. ulcerans, combining genomic and evolutionary approaches to decipher its evolutionary history and pattern of spread at different geographic scales. At the scale of villages in an endemic area of Benin, the circulating genotypes have been introduced in recent decades and are not randomly distributed along the river. On a global scale, M. ulcerans has been spreading for much longer, resulting in distinct and compartmentalized endemic foci across Africa and Australia.

Redefining the treponemal history through pre-Columbian genomes from Brazil

The origins of treponemal diseases have long remained unknown, especially considering the sudden onset of the first syphilis epidemic in the late 15th century in Europe and its hypothesized arrival from the Americas with Columbus' expeditions1,2. Recently, ancient DNA evidence has revealed various treponemal infections circulating in early modern Europe and colonial-era Mexico3-6. However, there has been to our knowledge no genomic evidence of treponematosis recovered from either the Americas or the Old World that can be reliably dated to the time before the first trans-Atlantic contacts. Here, we present treponemal genomes from nearly 2,000-year-old human remains from Brazil. We reconstruct four ancient genomes of a prehistoric treponemal pathogen, most closely related to the bejel-causing agent Treponema pallidum endemicum. Contradicting the modern day geographical niche of bejel in the arid regions of the world, the results call into question the previous palaeopathological characterization of treponeme subspecies and showcase their adaptive potential. A high-coverage genome is used to improve molecular clock date estimations, placing the divergence of modern T. pallidum subspecies firmly in pre-Columbian times. Overall, our study demonstrates the opportunities within archaeogenetics to uncover key events in pathogen evolution and emergence, paving the way to new hypotheses on the origin and spread of treponematoses.

The effect of missing data on evolutionary analysis of sequence capture bycatch, with application to an agricultural pest

Sequence capture is a genomic technique that selectively enriches target sequences before high throughput next-generation sequencing, to generate specific sequences of interest. Off-target or 'bycatch' data are often discarded from capture experiments, but can be leveraged to address evolutionary questions under some circumstances. Here, we investigated the effects of missing data on a variety of evolutionary analyses using bycatch from an exon capture experiment on the global pest moth, Helicoverpa armigera. We added > 200 new samples from across Australia in the form of mitogenomes obtained as bycatch from targeted sequence capture, and combined these into an additional larger dataset to total > 1000 mitochondrial cytochrome c oxidase subunit I (COI) sequences across the species' global distribution. Using discriminant analysis of principal components and Bayesian coalescent analyses, we showed that mitogenomes assembled from bycatch with up to 75% missing data were able to return evolutionary inferences consistent with higher coverage datasets and the broader literature surrounding H. armigera. For example, low-coverage sequences broadly supported the delineation of two H. armigera subspecies and also provided new insights into the potential for geographic turnover among these subspecies. However, we also identified key effects of dataset coverage and composition on our results. Thus, low-coverage bycatch data can offer valuable information for population genetic and phylodynamic analyses, but caution is required to ensure the reduced information does not introduce confounding factors, such as sampling biases, that drive inference. We encourage more researchers to consider maximizing the potential of the targeted sequence approach by examining evolutionary questions with their off-target bycatch where possible-especially in cases where no previous mitochondrial data exists-but recommend stratifying data at different genome coverage thresholds to separate sampling effects from genuine genomic signals, and to understand their implications for evolutionary research.

Ancient chicken remains reveal the origins of virulence in Marek's disease virus

The pronounced growth in livestock populations since the 1950s has altered the epidemiological and evolutionary trajectory of their associated pathogens. For example, Marek's disease virus (MDV), which causes lymphoid tumors in chickens, has experienced a marked increase in virulence over the past century. Today, MDV infections kill >90% of unvaccinated birds, and controlling it costs more than US$1 billion annually. By sequencing MDV genomes derived from archeological chickens, we demonstrate that it has been circulating for at least 1000 years. We functionally tested the Meq oncogene, one of 49 viral genes positively selected in modern strains, demonstrating that ancient MDV was likely incapable of driving tumor formation. Our results demonstrate the power of ancient DNA approaches to trace the molecular basis of virulence in economically relevant pathogens.

Parallel signatures of Mycobacterium tuberculosis and human Y-chromosome phylogeography support the Two Layer model of East Asian population history

The Two Layer hypothesis is fast becoming the favoured narrative describing East Asian population history. Under this model, hunter-gatherer groups who initially peopled East Asia via a route south of the Himalayas were assimilated by agriculturalist migrants who arrived via a northern route across Eurasia. A lack of ancient samples from tropical East Asia limits the resolution of this model. We consider insight afforded by patterns of variation within the human pathogen Mycobacterium tuberculosis (Mtb) by analysing its phylogeographic signatures jointly with the human Y-chromosome. We demonstrate the Y-chromosome lineages enriched in the traditionally hunter-gatherer groups associated with East Asia's first layer of peopling to display deep roots, low long-term effective population size, and diversity patterns consistent with a southern entry route. These characteristics mirror those of the evolutionarily ancient Mtb lineage 1. The remaining East Asian Y-chromosome lineage is almost entirely absent from traditionally hunter-gatherer groups and displays spatial and temporal characteristics which are incompatible with a southern entry route, and which link it to the development of agriculture in modern-day China. These characteristics mirror those of the evolutionarily modern Mtb lineage 2. This model paves the way for novel host-pathogen coevolutionary research hypotheses in East Asia.

An evolutionary epigenetic clock in plants

Molecular clocks are the basis for dating the divergence between lineages over macroevolutionary timescales (~105 to 108 years). However, classical DNA-based clocks tick too slowly to inform us about the recent past. Here, we demonstrate that stochastic DNA methylation changes at a subset of cytosines in plant genomes display a clocklike behavior. This "epimutation clock" is orders of magnitude faster than DNA-based clocks and enables phylogenetic explorations on a scale of years to centuries. We show experimentally that epimutation clocks recapitulate known topologies and branching times of intraspecies phylogenetic trees in the self-fertilizing plant Arabidopsis thaliana and the clonal seagrass Zostera marina, which represent two major modes of plant reproduction. This discovery will open new possibilities for high-resolution temporal studies of plant biodiversity.

A Glimpse into the Past: What Ancient Viral Genomes Reveal About Human History

Humans have battled viruses for millennia. However, directly linking the symptomatology of disease outbreaks to specific viral pathogens was not possible until the twentieth century. With the advent of the genomic era and the development of advanced protocols for isolation, sequencing, and analysis of ancient nucleic acids from diverse human remains, the identification and characterization of ancient viruses became feasible. Recent studies have provided invaluable information about past epidemics and made it possible to examine assumptions and inferences on the origin and evolution of certain viral families. In parallel, the study of ancient viruses also uncovered their importance in the evolution of the human lineage and their key roles in shaping major events in human history. In this review, we describe the strategies used for the study of ancient viruses, along with their limitations, and provide a detailed account of what past viral infections have revealed about human history.

Holocene deglaciation drove rapid genetic diversification of Atlantic walrus

Rapid global warming is severely impacting Arctic ecosystems and is predicted to transform the abundance, distribution and genetic diversity of Arctic species, though these linkages are poorly understood. We address this gap in knowledge using palaeogenomics to examine how earlier periods of global warming influenced the genetic diversity of Atlantic walrus (Odobenus rosmarus rosmarus), a species closely associated with sea ice and shallow-water habitats. We analysed 82 ancient and historical Atlantic walrus mitochondrial genomes (mitogenomes), including now-extinct populations in Iceland and the Canadian Maritimes, to reconstruct the Atlantic walrus' response to Arctic deglaciation. Our results demonstrate that the phylogeography and genetic diversity of Atlantic walrus populations was initially shaped by the last glacial maximum (LGM), surviving in distinct glacial refugia, and subsequently expanding rapidly in multiple migration waves during the late Pleistocene and early Holocene. The timing of diversification and establishment of distinct populations corresponds closely with the chronology of the glacial retreat, pointing to a strong link between walrus phylogeography and sea ice. Our results indicate that accelerated ice loss in the modern Arctic may trigger further dispersal events, likely increasing the connectivity of northern stocks while isolating more southerly stocks putatively caught in small pockets of suitable habitat.

Within-host diversity improves phylogenetic and transmission reconstruction of SARS-CoV-2 outbreaks

Accurate inference of who infected whom in an infectious disease outbreak is critical for the delivery of effective infection prevention and control. The increased resolution of pathogen whole-genome sequencing has significantly improved our ability to infer transmission events. Despite this, transmission inference often remains limited by the lack of genomic variation between the source case and infected contacts. Although within-host genetic diversity is common among a wide variety of pathogens, conventional whole-genome sequencing phylogenetic approaches exclusively use consensus sequences, which consider only the most prevalent nucleotide at each position and therefore fail to capture low-frequency variation within samples. We hypothesized that including within-sample variation in a phylogenetic model would help to identify who infected whom in instances in which this was previously impossible. Using whole-genome sequences from SARS-CoV-2 multi-institutional outbreaks as an example, we show how within-sample diversity is partially maintained among repeated serial samples from the same host, it can transmitted between those cases with known epidemiological links, and how this improves phylogenetic inference and our understanding of who infected whom. Our technique is applicable to other infectious diseases and has immediate clinical utility in infection prevention and control.

Herbarium specimen sequencing allows precise dating of Xanthomonas citri pv. citri diversification history

Herbarium collections are an important source of dated, identified and preserved DNA, whose use in comparative genomics and phylogeography can shed light on the emergence and evolutionary history of plant pathogens. Here, we reconstruct 13 historical genomes of the bacterial crop pathogen Xanthomonas citri pv. citri (Xci) from infected Citrus herbarium specimens. Following authentication based on ancient DNA damage patterns, we compare them with a large set of modern genomes to estimate their phylogenetic relationships, pathogenicity-associated gene content and several evolutionary parameters. Our results indicate that Xci originated in Southern Asia ~11,500 years ago (perhaps in relation to Neolithic climate change and the development of agriculture) and diversified during the beginning of the 13th century, after Citrus diversification and before spreading to the rest of the world (probably via human-driven expansion of citriculture through early East-West trade and colonization).

CNETML: maximum likelihood inference of phylogeny from copy number profiles of multiple samples

Phylogenetic trees based on copy number profiles from multiple samples of a patient are helpful to understand cancer evolution. Here, we develop a new maximum likelihood method, CNETML, to infer phylogenies from such data. CNETML is the first program to jointly infer the tree topology, node ages, and mutation rates from total copy numbers of longitudinal samples. Our extensive simulations suggest CNETML performs well on copy numbers relative to ploidy and under slight violation of model assumptions. The application of CNETML to real data generates results consistent with previous discoveries and provides novel early copy number events for further investigation.

Genomic surveillance uncovers a pandemic clonal lineage of the wheat blast fungus

Wheat, one of the most important food crops, is threatened by a blast disease pandemic. Here, we show that a clonal lineage of the wheat blast fungus recently spread to Asia and Africa following two independent introductions from South America. Through a combination of genome analyses and laboratory experiments, we show that the decade-old blast pandemic lineage can be controlled by the Rmg8 disease resistance gene and is sensitive to strobilurin fungicides. However, we also highlight the potential of the pandemic clone to evolve fungicide-insensitive variants and sexually recombine with African lineages. This underscores the urgent need for genomic surveillance to track and mitigate the spread of wheat blast outside of South America and to guide preemptive wheat breeding for blast resistance.

Genomic surveillance uncovers a pandemic clonal lineage of the wheat blast fungus

Wheat, one of the most important food crops, is threatened by a blast disease pandemic. Here, we show that a clonal lineage of the wheat blast fungus recently spread to Asia and Africa following two independent introductions from South America. Through a combination of genome analyses and laboratory experiments, we show that the decade-old blast pandemic lineage can be controlled by the Rmg8 disease resistance gene and is sensitive to strobilurin fungicides. However, we also highlight the potential of the pandemic clone to evolve fungicide-insensitive variants and sexually recombine with African lineages. This underscores the urgent need for genomic surveillance to track and mitigate the spread of wheat blast outside of South America and to guide preemptive wheat breeding for blast resistance.

Phylostems: a new graphical tool to investigate temporal signal of heterochronous sequences datasets

Motivation

Molecular tip-dating of phylogenetic trees is a growing discipline that uses DNA sequences sampled at different points in time to co-estimate the timing of evolutionary events with rates of molecular evolution. Importantly, such inferences should only be performed on datasets displaying sufficient temporal signal, a feature important to test prior to any tip-dating inference. For this purpose, the most popular method considered to-date has been the 'root-to-tip regression' which consist in fitting a linear regression of the number of substitutions accumulated from the root to the tips of a phylogenetic tree as a function of sampling times. The main limitation of the regression method, in its current implementation, relies in the fact that the temporal signal can only be tested at the whole-tree scale (i.e. its root).

Results

To overcome this limitation we introduce Phylostems, a new graphical user-friendly tool developed to investigate temporal signal within every clade of a phylogenetic tree. We provide a 'how to' guide by running Phylostems on an empirical dataset and supply guidance for results interpretation.

Availability and implementation

Phylostems is freely available at https://pvbmt-apps.cirad.fr/apps/phylostems.

Evolutionary history and global spatiotemporal pattern of alfalfa mosaic virus

Alfalfa mosaic virus (AMV) is an important plant virus causing considerable economic loss to alfalfa production. Knowledge of the evolutionary and demographic history of the pathogen is limited but essential to the development of effective and sustainable pathogen management schemes. In this study, we performed worldwide phylodynamic analyses of AMV based on 154 nucleotide sequences of the coat protein gene, sampled from 1985 to 2020, to understand the epidemiology of this pathogen. Bayesian phylogenetic reconstruction estimates that the crown group of AMV dates back to 1840 (95% credibility interval, 1687-1955). We revealed that AMV continuously evolves at a rate of 4.14 × 10^-4 substitutions/site/year (95% credibility interval, 1.04 × 10^-4 - 6.68 × 10^-4). Our phylogeographic analyses identified multiple migration links between Europe and other regions, implying that Europe played a key role in spreading the virus worldwide. Further analyses showed that the clustering pattern of AMV isolates is significantly correlated to geographic regions, indicating that geography-driven adaptation may be a factor that affects the evolution of AMV. Our findings may be potentially used in the development of effective control strategies for AMV.

Profiling genome-wide recombination in Epstein Barr virus reveals type-specific patterns and associations with endemic-Burkitt lymphoma

BACKGROUND

Endemic Burkitt lymphoma (eBL) is potentiated through the interplay of Epstein Barr virus (EBV) and holoendemic Plasmodium falciparum malaria. To better understand EBV's biology and role in eBL, we characterized genome-wide recombination sites and patterns as a source of genetic diversity in EBV genomes in our well-defined population of eBL cases and controls from Western Kenya.

METHODS

EBV genomes representing 54 eBL cases and 32 healthy children from the same geographic region in Western Kenya that we previously sequenced were analyzed. Whole-genome multiple sequence alignment, recombination analyses, and phylogenetic inference were made using multiple alignment with fast Fourier transform, recombination detection program 4, and molecular evolutionary genetics analysis.

RESULTS

We identified 28 different recombination events and 71 (82.6%) of the 86 EBV genomes analyzed contained evidence of one or more recombinant segments. Associated recombination breakpoints were found to occur in a total of 42 different genes, with only 7 (16.67%) being latent genes. Recombination events were major drivers of clustering within genome-wide phylogenetic trees. The occurrence of recombination segments was comparable between genomes from male and female participants and across age groups. More recombinant segments were found in EBV type 1 genomes (p = 6.4e - 06) and the genomes from the eBLs (p = 0.037). Two recombination events were enriched in the eBLs; event 47 (OR = 4.07, p = 0.038) and event 50 (OR = 14.24, p = 0.012).

CONCLUSIONS

EBV genomes have extensive evidence of recombination likely acquired progressively and cumulatively over time. Recombination patterns display a heterogeneous occurrence rate across the genome with enrichment in lytic genes. Overall, recombination appears to be a major evolutionary force impacting EBV diversity and genome structure with evidence of the association of specific recombinants with eBL.

Phytest: quality control for phylogenetic analyses

MOTIVATION

The ability to automatically conduct quality control checks on phylogenetic analyses is becoming more important with the increase in genetic sequencing and the use of real-time pipelines e.g. in the SARS-CoV-2 era. Implementations of real-time phylogenetic analyses require automated testing to make sure that problems in the data are caught automatically within analysis pipelines and in a timely manner. Here, we present Phytest (version 1.1) a tool for automating quality control checks on sequences, trees and metadata during phylogenetic analyses.

RESULTS

Phytest is a phylogenetic analysis testing program that easily integrates into existing phylogenetic pipelines. We demonstrate the utility of Phytest with real-world examples.

AVAILABILITY AND IMPLEMENTATION

Phytest source code is available on GitHub (https://github.com/phytest-devs/phytest) and can be installed via PyPI with the command 'pip install phytest'. Extensive documentation can be found at https://phytest-devs.github.io/phytest/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Exploring the Emergence and Evolution of Plant Pathogenic Microbes Using Historical and Paleontological Sources

Biotechnological advances now permit broad exploration of past microbial communities preserved in diverse substrates. Despite biomolecular degradation, high-throughput sequencing of preserved materials can yield invaluable genomic and metagenomic data from the past. This line of research has expanded from its initial human- and animal-centric foci to include plant-associated microbes (viruses, archaea, bacteria, fungi, and oomycetes), for which historical, archaeological, and paleontological data illuminate past epidemics and evolutionary history. Genetic mechanisms underlying the acquisition of microbial pathogenicity, including hybridization, polyploidization, and horizontal gene transfer, can now be reconstructed, as can gene-for-gene coevolution with plant hosts. Epidemiological parameters, such as geographic origin and range expansion, can also be assessed. Building on published case studies with individual phytomicrobial taxa, the stage is now set for broader, community-wide studies of preserved plant microbiomes to strengthen mechanistic understanding of microbial interactions and plant disease emergence.

First report of Odontoglossum ringspot virus in Vanilla (Orchidaceae) in Madagascar

Vanilla (Vanilla planifolia, Orchidaceae) is Madagascar's leading agricultural export resource which provides 80% of world's consumption. During a phytosanitary survey conducted from November 2019 to March 2021 in the main vanilla production regions of Madagascar, 250 plots were indexed for cymbidium mosaic virus (CymMV, Potexvirus genus) and odontoglossum ringspot virus (ORSV, Tobamovirus genus) the two most prevalent viruses of cultivated orchids worldwide (Zettler et al., 1990). For each plot, bulk samples (ten leaves taken at random) were assayed using Immunostrips (AGDIA, ISK 13301). A quarter of the plots (63/250) tested positive for CymMV. The highest prevalence of CymMV was observed in the SAVA region (57 out of 153 plots = 37,2%) where the virus has been reported since 1997 (Grisoni et al., 2010). Six plots located in the district of Mahanoro (Atsinanana) tested positive for ORSV. A few plants in these plots showed chlorotic often annular spots on their leaves. They were individually tested positive for ORSV, and negative for CymMV and potyviruses (Immunostrips AGDIA ISK 27200), the other two viruses reported so far in vanilla in Madagascar. To confirm the diagnosis of ORSV, leaf samples from five of the six infected plots were analysed by Tube Capture-RT-PCR (Grisoni et al., 2017) using two pairs of primers flanking the ORSV coat protein (CP) gene: OrCP1 (GGTCGGTAATGGTGTTAG) / OrCP2 (TGCATTATCGTATGCTCC), and CPOR-F(ATGTCTTACACTATTACAGACC) / CPOR-R(TTAGGAAGAGGTCCAAGTAAG). The five samples gave amplicons of the expected size (820 nt and 476 nt, respectively) and were sequenced with Sanger technology (Macrogen, The Netherlands). The ORSV-CP sequences of the Mahanoro isolates showed very close similarity to 198 ORSV-CP sequences from GenBank (95.8% to 99.6% nucleotide and 94.5 to 100% amino-acid identities), and less than 75.4% nucleotide (80.1% amino-acid) identities with Bell pepper mosaic virus (DQ355023), the tobamovirus closest to ORSV. The five ORSV-CP sequences from vanilla were deposited in GenBank under accessions numbers OM847399 to OM847403. These data confirmed that ORSV infects vanilla vines in Madagascar. To our knowledge, this is the first report of this virus in Madagascar and of its ability to infect symptomatically V. planifolia. The five ORSV isolates from vanilla had more than 98.7 % nucleotide identities of CP gene and clustered into a monophyletic group in maximum likelihood phylogenetic tree, suggesting a single origin of these isolates. To further investigate the origin of ORSV in Madagascar, we made use of RNA sequences isolated at different points in time to infer the timing of evolutionary events (Rieux et al., 2016). We estimated the CP gene substitution rate to 4.8E-4 subst/site/year [95%HPD 2.1E-4 - 8.7E-4] which is close to the estimate of He et al. (2019) based on a slightly different sequences set (1.25E-3 subst/site/year). We dated the initial contamination of vanilla plts by ORSV between 2004 and 2013. Both ORSV and CymMV have deleterious effects on many ornamental orchids, and the pathogenicity of CymMV is exacerbated when co-infecting with ORSV (Lee et al., 2021). Therefore, ORSV represents a new threat to the Malagasy vanilla crop, especially in regions where CymMV is already rife. Given the economic importance of vanilla cultivation in the country, the implementation of prophylactic measures aimed at preventing the spread of ORSV, in particular through the sanitary control of cuttings, should be a priority for the vanilla industry.

Variation in the substitution rates among the human mitochondrial haplogroup U sublineages

Resolving the absolute timescale of phylogenetic trees stipulates reliable estimates for the rate of DNA sequence evolution. For this end, various calibration methods have been developed and studied intensively. Intraspecific rate variation among distinct genetic lineages, however, has gained less attention. Here, we have assessed lineage-specific molecular rates of human mitochondrial DNA (mtDNA) by performing tip-calibrated Bayesian phylogenetic analyses. Tip-calibration, as opposed to traditional nodal time stamps from dated fossil evidence or geological events, is based on sample ages and becoming ever more feasible as ancient DNA data from radiocarbon-dated samples accumulate. We focus on subhaplogroups U2, U4, U5a, and U5b, the data including ancient mtDNA genomes from ¹⁴C-dated samples (n = 234), contemporary genomes (n = 301), and two outgroup sequences from haplogroup R. The obtained molecular rates depended on the data sets (with or without contemporary sequences), suggesting time-dependency. More notable was the rate variation between haplogroups: U4 and U5a stand out having a substantially higher rate than U5b. This is also reflected in the divergence times obtained (U5a: 17,700 years and U5b: 29,700 years), a disparity not reported previously. After ruling out various alternative causes (e.g., selection, sampling, and sequence quality), we propose that the substitution rates have been influenced by demographic histories, widely different among populations where U4/U5a or U5b are frequent. As with the Y-chromosomal subhaplogroup R1b, the mitochondrial U4 and U5a have been associated with remarkable range extensions of the Yamnaya culture in the Bronze Age.

Detection of Ancient Viruses and Long-Term Viral Evolution

The COVID-19 outbreak has reminded us of the importance of viral evolutionary studies as regards comprehending complex viral evolution and preventing future pandemics. A unique approach to understanding viral evolution is the use of ancient viral genomes. Ancient viruses are detectable in various archaeological remains, including ancient people's skeletons and mummified tissues. Those specimens have preserved ancient viral DNA and RNA, which have been vigorously analyzed in the last few decades thanks to the development of sequencing technologies. Reconstructed ancient pathogenic viral genomes have been utilized to estimate the past pandemics of pathogenic viruses within the ancient human population and long-term evolutionary events. Recent studies revealed the existence of non-pathogenic viral genomes in ancient people's bodies. These ancient non-pathogenic viruses might be informative for inferring their relationships with ancient people's diets and lifestyles. Here, we reviewed the past and ongoing studies on ancient pathogenic and non-pathogenic viruses and the usage of ancient viral genomes to understand their long-term viral evolution.

Within-host diversity improves phylogenetic and transmission reconstruction of SARS-CoV-2 outbreaks

Accurate inference of who infected whom in an infectious disease outbreak is critical for the delivery of effective infection prevention and control. The increased resolution of pathogen whole-genome sequencing has significantly improved our ability to infer transmission events. Despite this, transmission inference often remains limited by the lack of genomic variation between the source case and infected contacts. Although within-host genetic diversity is common among a wide variety of pathogens, conventional whole-genome sequencing phylogenetic approaches to reconstruct outbreaks exclusively use consensus sequences, which consider only the most prevalent nucleotide at each position and therefore fail to capture low frequency variation within samples. We hypothesized that including within-sample variation in a phylogenetic model would help to identify who infected whom in instances in which this was previously impossible. Using whole-genome sequences from SARS-CoV-2 multi-institutional outbreaks as an example, we show how within-sample diversity is stable among repeated serial samples from the same host, is transmitted between those cases with known epidemiological links, and how this improves phylogenetic inference and our understanding of who infected whom. Our technique is applicable to other infectious diseases and has immediate clinical utility in infection prevention and control.

Dating the Emergence of Human Endemic Coronaviruses

Four endemic coronaviruses infect humans and cause mild symptoms. Because previous analyses were based on a limited number of sequences and did not control for effects that affect molecular dating, we re-assessed the timing of endemic coronavirus emergence. After controlling for recombination, selective pressure, and molecular clock model, we obtained similar tMRCA (time to the most recent common ancestor) estimates for the four coronaviruses, ranging from 72 (HCoV-229E) to 54 (HCoV-NL63) years ago. The split times of HCoV-229E and HCoV-OC43 from camel alphacoronavirus and bovine coronavirus were dated ~268 and ~99 years ago. The split times of HCoV-HKU1 and HCoV-NL63 could not be calculated, as their zoonoticic sources are unknown. To compare the timing of coronavirus emergence to that of another respiratory virus, we recorded the occurrence of influenza pandemics since 1500. Although there is no clear relationship between pandemic occurrence and human population size, the frequency of influenza pandemics seems to intensify starting around 1700, which corresponds with the initial phase of exponential increase of human population and to the emergence of HCoV-229E. The frequency of flu pandemics in the 19th century also suggests that the concurrence of HCoV-OC43 emergence and the Russian flu pandemic may be due to chance.

Role of mobile genetic elements in the global dissemination of the carbapenem resistance gene blaNDM

The mobile resistance gene blaNDM encodes the NDM enzyme which hydrolyses carbapenems, a class of antibiotics used to treat some of the most severe bacterial infections. The blaNDM gene is globally distributed across a variety of Gram-negative bacteria on multiple plasmids, typically located within highly recombining and transposon-rich genomic regions, which leads to the dynamics underlying the global dissemination of blaNDM to remain poorly resolved. Here, we compile a dataset of over 6000 bacterial genomes harbouring the blaNDM gene, including 104 newly generated PacBio hybrid assemblies from clinical and livestock-associated isolates across China. We develop a computational approach to track structural variants surrounding blaNDM, which allows us to identify prevalent genomic contexts, mobile genetic elements, and likely events in the gene's global spread. We estimate that blaNDM emerged on a Tn125 transposon before 1985, but only reached global prevalence around a decade after its first recorded observation in 2005. The Tn125 transposon seems to have played an important role in early plasmid-mediated jumps of blaNDM, but was overtaken in recent years by other elements including IS26-flanked pseudo-composite transposons and Tn3000. We found a strong association between blaNDM-carrying plasmid backbones and the sampling location of isolates. This observation suggests that the global dissemination of the blaNDM gene was primarily driven by successive between-plasmid transposon jumps, with far more restricted subsequent plasmid exchange, possibly due to adaptation of plasmids to their specific bacterial hosts.

Genomic signatures of pre-resistance in Mycobacterium tuberculosis

Recent advances in bacterial whole-genome sequencing have resulted in a comprehensive catalog of antibiotic resistance genomic signatures in Mycobacterium tuberculosis. With a view to pre-empt the emergence of resistance, we hypothesized that pre-existing polymorphisms in susceptible genotypes (pre-resistance mutations) could increase the risk of becoming resistant in the future. We sequenced whole genomes from 3135 isolates sampled over a 17-year period. After reconstructing ancestral genomes on time-calibrated phylogenetic trees, we developed and applied a genome-wide survival analysis to determine the hazard of resistance acquisition. We demonstrate that M. tuberculosis lineage 2 has a higher risk of acquiring resistance than lineage 4, and estimate a higher hazard of rifampicin resistance evolution following isoniazid mono-resistance. Furthermore, we describe loci and genomic polymorphisms associated with a higher risk of resistance acquisition. Identifying markers of future antibiotic resistance could enable targeted therapy to prevent resistance emergence in M. tuberculosis and other pathogens.

Genomic Analysis of 18th-Century Kazakh Individuals and Their Oral Microbiome

The Asian Central Steppe, consisting of current-day Kazakhstan and Russia, has acted as a highway for major migrations throughout history. Therefore, describing the genetic composition of past populations in Central Asia holds value to understanding human mobility in this pivotal region. In this study, we analyse paleogenomic data generated from five humans from Kuygenzhar, Kazakhstan. These individuals date to the early to mid-18th century, shortly after the Kazakh Khanate was founded, a union of nomadic tribes of Mongol Golden Horde and Turkic origins. Genomic analysis identifies that these individuals are admixed with varying proportions of East Asian ancestry, indicating a recent admixture event from East Asia. The high amounts of DNA from the anaerobic Gram-negative bacteria Tannerella forsythia, a periodontal pathogen, recovered from their teeth suggest they may have suffered from periodontitis disease. Genomic analysis of this bacterium identified recently evolved virulence and glycosylation genes including the presence of antibiotic resistance genes predating the antibiotic era. This study provides an integrated analysis of individuals with a diet mostly based on meat (mainly horse and lamb), milk, and dairy products and their oral microbiome.

Genomic analysis of an outbreak of bovine tuberculosis in a man-made multi-host species system: A call for action on wildlife in Brazil

We report on a 15-year-long outbreak of bovine tuberculosis (bTB) in wildlife from a Brazilian safari park. A timeline of diagnostic events and whole-genome sequencing (WGS) of 21 Mycobacterium bovis isolates from deer and llamas were analyzed. Accordingly, from 2003 to 2018, at least 16 animals, from eight species, died due to TB, which is likely an underestimated number. In three occasions since 2013, the deer presented positive tuberculin tests, leading to the park closure and culling of all deer. WGS indicated that multiple M. bovis strains were circulating, with at least three founding introductions since the park inauguration in 1977. Using a previously sequenced dataset of 71 M. bovis genomes from cattle, we found no recent transmission events between nearby farms and the park based on WGS. Lastly, by discussing socio-economic and environmental factors escaping current regulatory gaps that were determinant of this outbreak, we pledge for the development of a plan to report and control bTB in wildlife in Brazil.

Phylodynamic and phylogeographic analysis of the complete genome of the West Nile virus lineage 2 (WNV-2) in the Mediterranean basin

Background

The West Nile virus is a highly contagious agent for a wide range of hosts. Its spread in the Mediterranean region raises several questions about its origin and the risk factors underlying the virus’s dispersal.

Materials and methods

The present study aims to reconstruct the temporal and spatial phylodynamics of West Nile virus lineage 2 in the Mediterranean region using 75 complete genome sequences from different host species retrieved from international databases.

Results

This data set suggests that current strains of WNV-2 began spreading in South Africa or nearby regions in the early twentieth century, and it migrated northwards via at least one route crossing the Mediterranean to reach Hungary in the early 2000s, before spreading throughout Europe. Another introduction event, according to the data set collected and analyses performed, is inferred to have occurred in around 1978. Migratory birds constitute, among others, additional risk factors that enhance the geographical transmission of the infection.

Conclusion

Our data underline the importance of the spatial–temporal tracking of migratory birds and phylodynamic reconstruction in setting up an efficient surveillance system for emerging and reemerging zoonoses in the Mediterranean region.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01902-w.

Quantifying bacterial evolution in the wild: A birthday problem for Campylobacter lineages

Measuring molecular evolution in bacteria typically requires estimation of the rate at which nucleotide changes accumulate in strains sampled at different times that share a common ancestor. This approach has been useful for dating ecological and evolutionary events that coincide with the emergence of important lineages, such as outbreak strains and obligate human pathogens. However, in multi-host (niche) transmission scenarios, where the pathogen is essentially an opportunistic environmental organism, sampling is often sporadic and rarely reflects the overall population, particularly when concentrated on clinical isolates. This means that approaches that assume recent common ancestry are not applicable. Here we present a new approach to estimate the molecular clock rate in Campylobacter that draws on the popular probability conundrum known as the 'birthday problem'. Using large genomic datasets and comparative genomic approaches, we use isolate pairs that share recent common ancestry to estimate the rate of nucleotide change for the population. Identifying synonymous and non-synonymous nucleotide changes, both within and outside of recombined regions of the genome, we quantify clock-like diversification to estimate synonymous rates of nucleotide change for the common pathogenic bacteria Campylobacter coli (2.4 x 10-6 s/s/y) and Campylobacter jejuni (3.4 x 10-6 s/s/y). Finally, using estimated total rates of nucleotide change, we infer the number of effective lineages within the sample time frame-analogous to a shared birthday-and assess the rate of turnover of lineages in our sample set over short evolutionary timescales. This provides a generalizable approach to calibrating rates in populations of environmental bacteria and shows that multiple lineages are maintained, implying that large-scale clonal sweeps may take hundreds of years or more in these species.

Cochlear shape distinguishes southern African early hominin taxa with unique auditory ecologies

Insights into potential differences among the bony labyrinths of Plio-Pleistocene hominins may inform their evolutionary histories and sensory ecologies. We use four recently-discovered bony labyrinths from the site of Kromdraai to significantly expand the sample for Paranthropus robustus. Diffeomorphometry, which provides detailed information about cochlear shape, reveals size-independent differences in cochlear shape between P. robustus and Australopithecus africanus that exceed those among modern humans and the African apes. The cochlea of P. robustus is distinctive and relatively invariant, whereas cochlear shape in A. africanus is more variable, resembles that of early Homo, and shows a degree of morphological polymorphism comparable to that evinced by modern species. The curvature of the P. robustus cochlea is uniquely derived and is consistent with enhanced sensitivity to low-frequency sounds. Combined with evidence for selection, our findings suggest that sound perception shaped distinct ecological adaptations among southern African early hominins.

Estimating the age of the subfamily Orthocoronavirinae using host divergence times as calibration ages at two internal nodes

Viruses of the subfamily Orthocoronavirinae can cause mild to severe disease in people, including COVID-19, MERS and SARS. Their most common natural hosts are bat and bird species, which are mostly split across four virus genera. Molecular clock analyses of orthocoronaviruses suggested the most recent common ancestor of these viruses might have emerged either around 10,000 years ago or, using models accounting for selection, many millions of years. Here, we reassess the evolutionary history of these viruses. We present time-aware phylogenetic analyses of a RNA-dependent RNA polymerase locus from 123 orthocoronaviruses isolated from birds and bats, including those in New Zealand, which were geographically isolated from other bats around 35 million years ago. We used this age, as well as the age of the avian-mammals split, to calibrate the molecular clocks, under the assumption that these ages are applicable to the analyzed viruses. We found that the time to the most recent ancestor common for all orthocoronaviruses is likely 150 or more million years, supporting clock analyses that account for selection.

Ancient viral genomes reveal introduction of human pathogenic viruses into Mexico during the transatlantic slave trade

After the European colonization of the Americas, there was a dramatic population collapse of the Indigenous inhabitants caused in part by the introduction of new pathogens. Although there is much speculation on the etiology of the Colonial epidemics, direct evidence for the presence of specific viruses during the Colonial era is lacking. To uncover the diversity of viral pathogens during this period, we designed an enrichment assay targeting ancient DNA (aDNA) from viruses of clinical importance and applied it to DNA extracts from individuals found in a Colonial hospital and a Colonial chapel (16th-18th century) where records suggest that victims of epidemics were buried during important outbreaks in Mexico City. This allowed us to reconstruct three ancient human parvovirus B19 genomes and one ancient human hepatitis B virus genome from distinct individuals. The viral genomes are similar to African strains, consistent with the inferred morphological and genetic African ancestry of the hosts as well as with the isotopic analysis of the human remains, suggesting an origin on the African continent. This study provides direct molecular evidence of ancient viruses being transported to the Americas during the transatlantic slave trade and their subsequent introduction to New Spain. Altogether, our observations enrich the discussion about the etiology of infectious diseases during the Colonial period in Mexico.

First historical genome of a crop bacterial pathogen from herbarium specimen: Insights into citrus canker emergence

Over the past decade, ancient genomics has been used in the study of various pathogens. In this context, herbarium specimens provide a precious source of dated and preserved DNA material, enabling a better understanding of plant disease emergences and pathogen evolutionary history. We report here the first historical genome of a crop bacterial pathogen, Xanthomonas citri pv. citri (Xci), obtained from an infected herbarium specimen dating back to 1937. Comparing the 1937 genome within a large set of modern genomes, we reconstructed their phylogenetic relationships and estimated evolutionary parameters using Bayesian tip-calibration inferences. The arrival of Xci in the South West Indian Ocean islands was dated to the 19th century, probably linked to human migrations following slavery abolishment. We also assessed the metagenomic community of the herbarium specimen, showed its authenticity using DNA damage patterns, and investigated its genomic features including functional SNPs and gene content, with a focus on virulence factors.

Ancient mitochondrial genomes from the Argentinian Pampas inform the early peopling of the Southern Cone of South America

The Southern Cone of South America (SCSA) is a key region for investigations about the peopling of the Americas. However, little is known about the eastern sector, the Argentinian Pampas. We analyzed 18 mitochondrial genomes-7 of which are novel-from human skeletal remains from 3 Early to Late Holocene archaeological sites. The Pampas present a distinctive genetic makeup compared to other Middle to Late Holocene pre-Columbian SCSA populations. We also report the earliest individuals carrying SCSA-specific mitochondrial haplogroups D1j and D1g from Early and Middle Holocene, respectively. Using these deep calibration time points in Bayesian phylogenetic reconstructions, we suggest that the first settlers of the Pampas were part of a single and rapid dispersal ∼15,600 years ago. Finally, we propose that present-day genetic differences between the Pampas and the rest of the SCSA are due to founder effects, genetic drift, and a partial population replacement ∼9,000 years ago.

Bayesian Evaluation of Temporal Signal in Measurably Evolving Populations

Phylogenetic methods can use the sampling times of molecular sequence data to calibrate the molecular clock, enabling the estimation of evolutionary rates and timescales for rapidly evolving pathogens and data sets containing ancient DNA samples. A key aspect of such calibrations is whether a sufficient amount of molecular evolution has occurred over the sampling time window, that is, whether the data can be treated as having come from a measurably evolving population. Here, we investigate the performance of a fully Bayesian evaluation of temporal signal (BETS) in sequence data. The method involves comparing the fit to the data of two models: a model in which the data are accompanied by the actual (heterochronous) sampling times, and a model in which the samples are constrained to be contemporaneous (isochronous). We conducted simulations under a wide range of conditions to demonstrate that BETS accurately classifies data sets according to whether they contain temporal signal or not, even when there is substantial among-lineage rate variation. We explore the behavior of this classification in analyses of five empirical data sets: modern samples of A/H1N1 influenza virus, the bacterium Bordetella pertussis, coronaviruses from mammalian hosts, ancient DNA from Hepatitis B virus, and mitochondrial genomes of dog species. Our results indicate that BETS is an effective alternative to other tests of temporal signal. In particular, this method has the key advantage of allowing a coherent assessment of the entire model, including the molecular clock and tree prior which are essential aspects of Bayesian phylodynamic analyses.

Development of Phylodynamic Methods for Bacterial Pathogens

Phylodynamic methods have been essential to understand the interplay between the evolution and epidemiology of infectious diseases. To date, the field has centered on viruses. Bacterial pathogens are seldom analyzed under such phylodynamic frameworks, due to their complex genome evolution and, until recently, a paucity of whole-genome sequence data sets with rich associated metadata. We posit that the increasing availability of bacterial genomes and epidemiological data means that the field is now ripe to lay the foundations for applying phylodynamics to bacterial pathogens. The development of new methods that integrate more complex genomic and ecological data will help to inform public heath surveillance and control strategies for bacterial pathogens that represent serious threats to human health.

Cross-Disciplinary Genomics Approaches to Studying Emerging Fungal Infections

Emerging fungal pathogens pose a serious, global and growing threat to food supply systems, wild ecosystems, and human health. However, historic chronic underinvestment in their research has resulted in a limited understanding of their epidemiology relative to bacterial and viral pathogens. Therefore, the untargeted nature of genomics and, more widely, -omics approaches is particularly attractive in addressing the threats posed by and illuminating the biology of these pathogens. Typically, research into plant, human and wildlife mycoses have been largely separated, with limited dialogue between disciplines. However, many serious mycoses facing the world today have common traits irrespective of host species, such as plastic genomes; wide host ranges; large population sizes and an ability to persist outside the host. These commonalities mean that -omics approaches that have been productively applied in one sphere and may also provide important insights in others, where these approaches may have historically been underutilised. In this review, we consider the advances made with genomics approaches in the fields of plant pathology, human medicine and wildlife health and the progress made in linking genomes to other -omics datatypes and sets; we identify the current barriers to linking -omics approaches and how these are being underutilised in each field; and we consider how and which -omics methodologies it is most crucial to build capacity for in the near future.

Introduction of Mycobacterium ulcerans disease in the Bankim Health District of Cameroon follows damming of the Mapé River

Buruli ulcer (BU) is an emerging ulcerative skin disease caused by infection with Mycobacterium ulcerans. Efforts to control its spread have been hampered by our limited understanding of M. ulcerans reservoirs and transmission, and the factors leading to the emergence of BU disease in a particular region. In this report we investigate an anecdotal link between damming the Mapé River in Cameroon and the emergence of BU in the Health Districts bordering Lake Bankim, the impoundment created by the Mapé dam. We used bacterial population genomics and molecular dating to find compelling support for a 2000 M. ulcerans introduction event that followed about 10 years after the filling of the newly created impoundment in 1988. We compared the genomic reconstructions with high-resolution satellite imagery to investigate what major environmental alterations might have driven the emergence of the new focus.

Ancient Bacterial Genomes Reveal a High Diversity of Treponema pallidum Strains in Early Modern Europe

Syphilis is a globally re-emerging disease, which has marked European history with a devastating epidemic at the end of the 15^th century. Together with non-venereal treponemal diseases, like bejel and yaws, which are found today in subtropical and tropical regions, it currently poses a substantial health threat worldwide. The origins and spread of treponemal diseases remain unresolved, including syphilis' potential introduction into Europe from the Americas. Here, we present the first genetic data from archaeological human remains reflecting a high diversity of Treponema pallidum in early modern Europe. Our study demonstrates that a variety of strains related to both venereal syphilis and yaws-causing T. pallidum subspecies were already present in Northern Europe in the early modern period. We also discovered a previously unknown T. pallidum lineage recovered as a sister group to yaws- and bejel-causing lineages. These findings imply a more complex pattern of geographical distribution and etiology of early treponemal epidemics than previously understood.

Genomic epidemiology and evolutionary dynamics of respiratory syncytial virus group B in Kilifi, Kenya, 2015-17

Respiratory syncytial virus (RSV) circulates worldwide, occurring seasonally in communities, and is a leading cause of acute respiratory illness in young children. There is paucity of genomic data from purposively sampled populations by which to investigate evolutionary dynamics and transmission patterns of RSV. Here we present an analysis of 295 RSV group B (RSVB) genomes from Kilifi, coastal Kenya, sampled from individuals seeking outpatient care in nine health facilities across a defined geographical area (∼890 km²), over two RSV epidemics between 2015 and 2017. RSVB diversity was characterized by multiple virus introductions into the area and co-circulation of distinct genetic clusters, which transmitted and diversified locally with varying frequency. Increase in relative genetic diversity paralleled seasonal virus incidence. Importantly, we identified a cluster of viruses that emerged in the 2016/17 epidemic, carrying distinct amino-acid signatures including a novel nonsynonymous change (K68Q) in antigenic site ∅ in the Fusion protein. RSVB diversity was additionally marked by signature nonsynonymous substitutions that were unique to particular genomic clusters, some under diversifying selection. Our findings provide insights into recent evolutionary and epidemiological behaviors of RSVB, and highlight possible emergence of a novel antigenic variant, which has implications on current prophylactic strategies in development.

Genetic Strain Diversity of Multi-Host RNA Viruses that Infect a Wide Range of Pollinators and Associates is Shaped by Geographic Origins

Emerging viruses have caused concerns about pollinator population declines, as multi-host RNA viruses may pose a health threat to pollinators and associated arthropods. In order to understand the ecology and impact these viruses have, we studied their host range and determined to what extent host and spatial variation affect strain diversity. Firstly, we used RT-PCR to screen pollinators and associates, including honey bees (Apis mellifera) and invasive Argentine ants (Linepithema humile), for virus presence and replication. We tested for the black queen cell virus (BQCV), deformed wing virus (DWV), and Kashmir bee virus (KBV) that were initially detected in bees, and the two recently discovered Linepithema humile bunya-like virus 1 (LhuBLV1) and Moku virus (MKV). DWV, KBV, and MKV were detected and replicated in a wide range of hosts and commonly co-infected hymenopterans. Secondly, we placed KBV and DWV in a global phylogeny with sequences from various countries and hosts to determine the association of geographic origin and host with shared ancestry. Both phylogenies showed strong geographic rather than host-specific clustering, suggesting frequent inter-species virus transmission. Transmission routes between hosts are largely unknown. Nonetheless, avoiding the introduction of non-native species and diseased pollinators appears important to limit spill overs and disease emergence.

Plasmodium vivax Malaria Viewed through the Lens of an Eradicated European Strain

The protozoan Plasmodium vivax is responsible for 42% of all cases of malaria outside Africa. The parasite is currently largely restricted to tropical and subtropical latitudes in Asia, Oceania, and the Americas. Though, it was historically present in most of Europe before being finally eradicated during the second half of the 20th century. The lack of genomic information on the extinct European lineage has prevented a clear understanding of historical population structuring and past migrations of P. vivax. We used medical microscope slides prepared in 1944 from malaria-affected patients from the Ebro Delta in Spain, one of the last footholds of malaria in Europe, to generate a genome of a European P. vivax strain. Population genetics and phylogenetic analyses placed this strain basal to a cluster including samples from the Americas. This genome allowed us to calibrate a genomic mutation rate for P. vivax, and to estimate the mean age of the last common ancestor between European and American strains to the 15th century. This date points to an introduction of the parasite during the European colonization of the Americas. In addition, we found that some known variants for resistance to antimalarial drugs, including Chloroquine and Sulfadoxine, were already present in this European strain, predating their use. Our results shed light on the evolution of an important human pathogen and illustrate the value of antique medical collections as a resource for retrieving genomic information on pathogens from the past.

Chytrid fungi and global amphibian declines

Discovering that chytrid fungi cause chytridiomycosis in amphibians represented a paradigm shift in our understanding of how emerging infectious diseases contribute to global patterns of biodiversity loss. In this Review we describe how the use of multidisciplinary biological approaches has been essential to pinpointing the origins of amphibian-parasitizing chytrid fungi, including Batrachochytrium dendrobatidis and Batrachochytrium salamandrivorans, as well as to timing their emergence, tracking their cycles of expansion and identifying the core mechanisms that underpin their pathogenicity. We discuss the development of the experimental methods and bioinformatics toolkits that have provided a fuller understanding of batrachochytrid biology and informed policy and control measures.

A new method for inferring timetrees from temporally sampled molecular sequences

Pathogen timetrees are phylogenies scaled to time. They reveal the temporal history of a pathogen spread through the populations as captured in the evolutionary history of strains. These timetrees are inferred by using molecular sequences of pathogenic strains sampled at different times. That is, temporally sampled sequences enable the inference of sequence divergence times. Here, we present a new approach (RelTime with Dated Tips [RTDT]) to estimating pathogen timetrees based on a relative rate framework underlying the RelTime approach that is algebraic in nature and distinct from all other current methods. RTDT does not require many of the priors demanded by Bayesian approaches, and it has light computing requirements. In analyses of an extensive collection of computer-simulated datasets, we found the accuracy of RTDT time estimates and the coverage probabilities of their confidence intervals (CIs) to be excellent. In analyses of empirical datasets, RTDT produced dates that were similar to those reported in the literature. In comparative benchmarking with Bayesian and non-Bayesian methods (LSD, TreeTime, and treedater), we found that no method performed the best in every scenario. So, we provide a brief guideline for users to select the most appropriate method in empirical data analysis. RTDT is implemented for use via a graphical user interface and in high-throughput settings in the newest release of cross-platform MEGA X software, freely available from http://www.megasoftware.net.

Pathogen timetrees trace the origins and evolutionary histories of strains in populations, hosts, and outbreaks. The tips of these molecular phylogenies often contain sampling time information because the sequences were generally obtained at different times during the disease outbreaks and propagation. We have developed a new method for inferring divergence times and confidence intervals for phylogenies with tip dates. The new Relative Times with Dated Tips (RTDT) methods showed excellent performance in the analysis of computer-simulated datasets, producing similar or better results in several evolutionary scenarios as compared to other fast, non-Bayesian methods. The new method is available in the cross-platform MEGA software package (version 10.1 and higher) that provides a graphical user interface and allows usage via a command line in scripting and high throughput analysis (www.megasoftware.net).

CRISPR elements provide a new framework for the genealogy of the citrus canker pathogen Xanthomonas citri pv. citri

Background

Xanthomonads are an important clade of Gram-negative bacteria infecting a plethora of economically important host plants, including citrus. Knowledge about the pathogen’s diversity and population structure are prerequisite for epidemiological surveillance and efficient disease management. Rapidly evolving genetic loci, such as Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), are of special interest to develop new molecular typing tools.

Results

We analyzed CRISPR loci of 56 Xanthomonas citri pv. citri strains of world-wide origin, a regulated pathogen causing Asiatic citrus canker in several regions of the world. With one exception, 23 unique sequences built up the repertoire of spacers, suggesting that this set of strains originated from a common ancestor that already harbored these 23 spacers. One isolate originating from Pakistan contained a string of 14 additional, probably more recently acquired spacers indicating that this genetic lineage has or had until recently the capacity to acquire new spacers. Comparison of CRISPR arrays with previously obtained molecular typing data, such as amplified fragment length polymorphisms (AFLP), variable-number of tandem-repeats (VNTR) and genome-wide single-nucleotide polymorphisms (SNP), demonstrated that these methods reveal similar evolutionary trajectories. Notably, genome analyses allowed to generate a model for CRISPR array evolution in X. citri pv. citri, which provides a new framework for the genealogy of the citrus canker pathogen.

Conclusions

CRISPR-based typing will further improve the accuracy of the genetic identification of X. citri pv. citri outbreak strains in molecular epidemiology analyses, especially when used concomitantly with another genotyping method.

PhySpeTree: an automated pipeline for reconstructing phylogenetic species trees

Background

Phylogenetic species trees are widely used in inferring evolutionary relationships. Existing software and algorithms mainly focus on phylogenetic inference. However, less attention has been paid to intermediate steps, such as processing extremely large sequences and preparing configure files to connect multiple software. When the species number is large, the intermediate steps become a bottleneck that may seriously affect the efficiency of tree building.

Results

Here, we present an easy-to-use pipeline named PhySpeTree to facilitate the reconstruction of species trees across bacterial, archaeal, and eukaryotic organisms. Users need only to input the abbreviations of species names; PhySpeTree prepares complex configure files for different software, then automatically downloads genomic data, cleans sequences, and builds trees. PhySpeTree allows users to perform critical steps such as sequence alignment and tree construction by adjusting advanced options. PhySpeTree provides two parallel pipelines based on concatenated highly conserved proteins and small subunit ribosomal RNA sequences, respectively. Accessory modules, such as those for inserting new species, generating visualization configurations, and combining trees, are distributed along with PhySpeTree.

Conclusions

Together with accessory modules, PhySpeTree significantly simplifies tree reconstruction. PhySpeTree is implemented in Python running on modern operating systems (Linux, macOS, and Windows). The source code is freely available with detailed documentation (https://github.com/yangfangs/physpetools).