tipdatingbeast: an r package to assist the implementation of phylogenetic tip-dating tests using beast

Genomic epidemiology of Vibrio cholerae reveals the regional and global spread of two epidemic non-toxigenic lineages

Non-toxigenic Vibrio cholerae isolates have been found associated with diarrheal disease globally, however, the global picture of non-toxigenic infections is largely unknown. Among non-toxigenic V. cholerae, ctxAB negative, tcpA positive (CNTP) isolates have the highest risk of disease. From 2001 to 2012, 71 infectious diarrhea cases were reported in Hangzhou, China, caused by CNTP serogroup O1 isolates. We sequenced 119 V. cholerae genomes isolated from patients, carriers and the environment in Hangzhou between 2001 and 2012, and compared them with 850 publicly available global isolates. We found that CNTP isolates from Hangzhou belonged to two distinctive lineages, named L3b and L9. Both lineages caused disease over a long time period with usually mild or moderate clinical symptoms. Within Hangzhou, the spread route of the L3b lineage was apparently from rural to urban areas, with aquatic food products being the most likely medium. Both lineages had been previously reported as causing local endemic disease in Latin America, but here we show that global spread of them has occurred, with the most likely origin of L3b lineage being in Central Asia. The L3b lineage has spread to China on at least three occasions. Other spread events, including from China to Thailand and to Latin America were also observed. We fill the missing links in the global spread of the two non-toxigenic serogroup O1 V. cholerae lineages that can cause human infection. The results are important for the design of future disease control strategies: surveillance of V. cholerae should not be limited to ctxAB positive strains.

Phyloanatomic characterization of the distinct T cell and monocyte contributions to the peripheral blood HIV population within the host

Human immunodeficiency virus (HIV) is a rapidly evolving virus, allowing its genetic sequence to act as a fingerprint for epidemiological processes among, as well as within, individual infected hosts. Though primarily infecting the CD4+ T-cell population, HIV can also be found in monocytes, an immune cell population that differs in several aspects from the canonical T-cell viral target. Using single genome viral sequencing and statistical phylogenetic inference, we investigated the viral RNA diversity and relative contribution of each of these immune cell types to the viral population within the peripheral blood. Results provide evidence of an increased prevalence of circulating monocytes harboring virus in individuals with high viral load in the absence of suppressive antiretroviral therapy. Bayesian phyloanatomic analysis of three of these individuals demonstrated a measurable role for these cells, but not the circulating T-cell population, as a source of cell-free virus in the plasma, supporting the hypothesis that these cells can act as an additional conduit of virus spread.

The Effect of Sample Bias and Experimental Artefacts on the Statistical Phylogenetic Analysis of Picornaviruses

Statistical phylogenetic methods are a powerful tool for inferring the evolutionary history of viruses through time and space. The selection of mathematical models and analysis parameters has a major impact on the outcome, and has been relatively well-described in the literature. The preparation of a sequence dataset is less formalized, but its impact can be even more profound. This article used simulated datasets of enterovirus sequences to evaluate the effect of sample bias on picornavirus phylogenetic studies. Possible approaches to the reduction of large datasets and their potential for introducing additional artefacts were demonstrated. The most consistent results were obtained using “smart sampling”, which reduced sequence subsets from large studies more than those from smaller ones in order to preserve the rare sequences in a dataset. The effect of sequences with technical or annotation errors in the Bayesian framework was also analyzed. Sequences with about 0.5% sequencing errors or incorrect isolation dates altered by just 5 years could be detected by various approaches, but the efficiency of identification depended upon sequence position in a phylogenetic tree. Even a single erroneous sequence could profoundly destabilize the whole analysis by increasing the variance of the inferred evolutionary parameters.

Genomic epidemiology of syphilis reveals independent emergence of macrolide resistance across multiple circulating lineages

Syphilis is a sexually transmitted infection caused by Treponema pallidum subspecies pallidum and may lead to severe complications. Recent years have seen striking increases in syphilis in many countries. Previous analyses have suggested one lineage of syphilis, SS14, may have expanded recently, indicating emergence of a single pandemic azithromycin-resistant cluster. Here we use direct sequencing of T. pallidum combined with phylogenomic analyses to show that both SS14- and Nichols-lineages are simultaneously circulating in clinically relevant populations in multiple countries. We correlate the appearance of genotypic macrolide resistance with multiple independently evolved SS14 sub-lineages and show that genotypically resistant and sensitive sub-lineages are spreading contemporaneously. These findings inform our understanding of the current syphilis epidemic by demonstrating how macrolide resistance evolves in Treponema subspecies and provide a warning on broader issues of antimicrobial resistance.

Two Coxsackievirus B3 outbreaks associated with hand, foot, and mouth disease in China and the evolutionary history worldwide

BACKGROUND

Coxsackievirus B3 (CV-B3) is usually associated with aseptic meningitis and myocarditis; however, the association between CV-B3 and hand, foot, and mouth disease (HFMD) has not been clearly demonstrated, and the phylogenetic dynamics and transmission history of CV-B3 have not been well summarized.

METHOD

Two HFMD outbreaks caused by CV-B3 were described in Hebei Province in 2012 and in Shandong Province in 2016 in China. To analyze the epidemiological features of two CV-B3 outbreaks, a retrospective analysis was conducted. All clinical specimens from CV-B3 outbreaks were collected and disposed according to the standard procedures supported by the WHO Global Poliovirus Specialized Laboratory. EV genotyping and phylogenetic analysis were performed to illustrate the genetic characteristics of CV-B3 in China and worldwide.

RESULTS

Two transmissible lineages (lineage 2 and 3) were observed in Northern China, which acted as an important "reservoir" for the transmission of CV-B3. Sporadic exporting and importing of cases were observed in almost all regions. In addition, the global sequences of CV-B3 showed a tendency of geographic-specific clustering, indicating that geographic-driven adaptation plays a major role in the diversification and evolution of CV-B3.

CONCLUSIONS

Overall, our study indicated that CV-B3 is a causative agent of HFMD outbreak and revealed the phylogenetic dynamics of CV-B3 worldwide, as well as provided an insight on CV-B3 outbreaks for effective intervention and countermeasures.

Evolutionary history and palaeoecology of brown bear in North-East Siberia re-examined using ancient DNA and stable isotopes from skeletal remains

Over 60% of the modern distribution range of brown bears falls within Russia, yet palaeoecological data from the region remain scarce. Complete modern Russian brown bear mitogenomes are abundant in the published literature, yet examples of their ancient counterparts are absent. Similarly, there is only limited stable isotopic data of prehistoric brown bears from the region. We used ancient DNA and stable carbon (δ¹³C) and nitrogen (δ¹⁵N) isotopes retrieved from five Pleistocene Yakutian brown bears (one Middle Pleistocene and four Late Pleistocene), to elucidate the evolutionary history and palaeoecology of the species in the region. We were able to reconstruct the complete mitogenome of one of the Late Pleistocene specimens, but we were unable to assign it to any of the previously published brown bear mitogenome clades. A subsequent analysis of published mtDNA control region sequences, which included sequences of extinct clades from other geographic regions, assigned the ancient Yakutian bear to the extinct clade 3c; a clade previously identified from Late Quaternary specimens from Eastern Beringia and Northern Spain. Our analyses of stable isotopes showed relatively high δ¹⁵N values in the Pleistocene Yakutian brown bears, suggesting a more carnivorous diet than contemporary brown bears from Eastern Beringia.

Bayesian phylodynamic analysis reveals the dispersal patterns of tobacco mosaic virus in China

Tobacco mosaic virus (TMV) is widespread in China and causes considerable economic losses to tobacco production. The molecular epidemiology of this virus is, however, poorly understood. In this study, we sequenced the genomes of 51 TMV isolates from five tobacco-producing regions in China and investigated the dispersal patterns of this virus. Our phylogenetic analysis showed that TMV might have been introduced to China in the early 1900s, probably first to southwest China. However, TMV then moved to the north of the country, where it expanded. The north became the main seeding region for the subsequent movements of the virus within China. The north-to-south movement of TMV coincides with a shift of major tobacco-producing areas from north to south in this century, suggesting a link between human activities and the dispersal of TMV in China.

Genomic insights into the 2016-2017 cholera epidemic in Yemen

Yemen is currently experiencing, to our knowledge, the largest cholera epidemic in recent history. The first cases were declared in September 2016, and over 1.1 million cases and 2,300 deaths have since been reported1. Here we investigate the phylogenetic relationships, pathogenesis and determinants of antimicrobial resistance by sequencing the genomes of Vibrio cholerae isolates from the epidemic in Yemen and recent isolates from neighbouring regions. These 116 genomic sequences were placed within the phylogenetic context of a global collection of 1,087 isolates of the seventh pandemic V. cholerae serogroups O1 and O139 biotype El Tor2-4. We show that the isolates from Yemen that were collected during the two epidemiological waves of the epidemic1-the first between 28 September 2016 and 23 April 2017 (25,839 suspected cases) and the second beginning on 24 April 2017 (more than 1 million suspected cases)-are V. cholerae serotype Ogawa isolates from a single sublineage of the seventh pandemic V. cholerae O1 El Tor (7PET) lineage. Using genomic approaches, we link the epidemic in Yemen to global radiations of pandemic V. cholerae and show that this sublineage originated from South Asia and that it caused outbreaks in East Africa before appearing in Yemen. Furthermore, we show that the isolates from Yemen are susceptible to several antibiotics that are commonly used to treat cholera and to polymyxin B, resistance to which is used as a marker of the El Tor biotype.

Phylogenomic Analysis of Extraintestinal Pathogenic Escherichia coli Sequence Type 1193, an Emerging Multidrug-Resistant Clonal Group

The fluoroquinolone-resistant sequence type 1193 (ST1193) of Escherichia coli, from the ST14 clonal complex (STc14) within phylogenetic group B2, has appeared recently as an important cause of extraintestinal disease in humans. Although this emerging lineage has been characterized to some extent using conventional methods, it has not been studied extensively at the genomic level. Here, we used whole-genome sequence analysis to compare 355 ST1193 isolates with 72 isolates from other STs within STc14. Using core genome phylogeny, the ST1193 isolates formed a tightly clustered clade with many genotypic similarities, unlike ST14 isolates. All ST1193 isolates possessed the same set of three chromosomal mutations conferring fluoroquinolone resistance, carried the fimH64 allele, and were lactose non-fermenting. Analysis revealed an evolutionary progression from K1 to K5 capsular types and acquisition of an F-type virulence plasmid, followed by changes in plasmid structure congruent with genome phylogeny. In contrast, the numerous identified antimicrobial resistance genes were distributed incongruently with the underlying phylogeny, suggesting frequent gain or loss of the corresponding resistance gene cassettes despite retention of the presumed carrier plasmids. Pangenome analysis revealed gains and losses of genetic loci occurring during the transition from ST14 to ST1193 and from the K1 to K5 capsular types. Using time-scaled phylogenetic analysis, we estimated that current ST1193 clades first emerged approximately 25 years ago. Overall, ST1193 appears to be a recently emerged clone in which both stepwise and mosaic evolution have contributed to epidemiologic success.

Integrative approach using Yersinia pestis genomes to revisit the historical landscape of plague during the Medieval Period

Over the last few years, genomic studies on Yersinia pestis, the causative agent of all known plague epidemics, have considerably increased in numbers, spanning a period of about 5,000 y. Nonetheless, questions concerning historical reservoirs and routes of transmission remain open. Here, we present and describe five genomes from the second half of the 14th century and reconstruct the evolutionary history of Y. pestis by reanalyzing previously published genomes and by building a comprehensive phylogeny focused on strains attributed to the Second Plague Pandemic (14th to 18th century). Corroborated by historical and ecological evidence, the presented phylogeny, which includes our Y. pestis genomes, could support the hypothesis of an entry of plague into Western European ports through distinct waves of introduction during the Medieval Period, possibly by means of fur trade routes, as well as the recirculation of plague within the human population via trade routes and human movement.

Europe was a hub for the global spread of potato virus S in the 19th century

Potato virus S (PVS) is a major plant pathogen that causes considerable losses in global potato production. Knowledge of the evolutionary history and spatio-temporal dynamics of PVS is vital for developing sustainable management schemes. In this study, we investigated the phylodynamics of the virus by analysing 103 nucleotide sequences of the coat protein gene, sampled between 1985 and 2014. Our Bayesian phylogenetic analyses showed that PVS has been evolving at a rate of 3.32 × 10^-4 substitutions/site/year (95% credibility interval 1.33 × 10^-4-5.58 × 10^-4). We dated the crown group to the year 1325 CE (95% credibility interval 762-1743 CE). Our phylogeographic analyses pointed to viral origins in South America and identified multiple migration pathways between Europe and other regions, suggesting that Europe has been a major hub for PVS transmission. The results of our study have potential implications for developing effective strategies for the control of this pathogen.

Global expansion of Mycobacterium tuberculosis lineage 4 shaped by colonial migration and local adaptation

On the basis of population genomic and phylogeographic analyses of 1669 Mycobacterium tuberculosis lineage 4 (L4) genomes, we find that dispersal of L4 has been completely dominated by historical migrations out of Europe. We demonstrate an intimate temporal relationship between European colonial expansion into Africa and the Americas and the spread of L4 tuberculosis (TB). Markedly, in the age of antibiotics, mutations conferring antimicrobial resistance overwhelmingly emerged locally (at the level of nations), with minimal cross-border transmission of resistance. The latter finding was found to reflect the relatively recent emergence of these mutations, as a similar degree of local restriction was observed for susceptible variants emerging on comparable time scales. The restricted international transmission of drug-resistant TB suggests that containment efforts at the level of individual countries could be successful.

Adaptive evolution and demographic history contribute to the divergent population genetic structure of Potato virus Y between China and Japan

Potato virus Y (PVY) is an important plant pathogen causing considerable economic loss to potato production. Knowledge of the population genetic structure and evolutionary biology of the pathogen, particularly at a transnational scale, is limited but vital in developing sustainable management schemes. In this study, the population genetic structure and molecular evolution of PVY were studied using 127 first protein (P1) and 137 coat protein (CP) sequences generated from isolates collected from potato in China and Japan. High genetic differentiation was found between the populations from the two countries, with higher nucleotide diversity in Japan than China in both genes and a KST value of .216 in the concatenated sequences of the two genes. Sequences from the two countries clustered together according to their geographic origin. Further analyses showed that spatial genetic structure in the PVY populations was likely caused by demographic dynamics of the pathogen and natural selection generated by habitat heterogeneity. Purifying selection was detected at the majority of polymorphic sites although some clade-specific codons were under positive selection. In past decades, PVY has undergone a population expansion in China, whereas in Japan, the population size of the pathogen has remained relatively constant.

DNA analysis of a 30,000-year-old Urocitellus glacialis from northeastern Siberia reveals phylogenetic relationships between ancient and present-day arctic ground squirrels

In contrast to the abundant fossil record of arctic ground squirrels, Urocitellus parryii, from eastern Beringia, only a limited number of fossils is known from its western part. In 1946, unnamed GULAG prisoners discovered a nest with three mummified carcasses of arctic ground squirrels in the permafrost sediments of the El’ga river, Yakutia, Russia, that were later attributed to a new species, Citellus (Urocitellus) glacialis Vinogr. To verify this assignment and to explore phylogenetic relationships between ancient and present-day arctic ground squirrels, we performed ¹⁴C dating and ancient DNA analyses of one of the El’ga mummies and four contemporaneous fossils from Duvanny Yar, northeastern Yakutia. Phylogenetic reconstructions, based on complete cytochrome b gene sequences of five Late Pleistocene arctic ground squirrels and those of modern U. parryii from 21 locations across western Beringia, provided no support for earlier proposals that ancient arctic ground squirrels from Siberia constitute a distinct species. In fact, we observed genetic continuity of the glacialis mitochondrial DNA lineage in modern U. parryii of the Kamchatka peninsula. When viewed in a broader geographic perspective, our findings provide new insights into the genetic history of U. parryii in Late Pleistocene Beringia.