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

Antimicrobial resistance and population genomics of emerging multidrug-resistant Salmonella 4,[5],12:i:- in Guangdong, China

Salmonella 4,[5],12:i:-, a monophasic variant of Salmonella Typhimurium, has emerged as a global cause of multidrug-resistant salmonellosis and has become endemic in many developing and developed countries, especially in China. Here, we have sequenced 352 clinical isolates in Guangdong, China, during 2009-2019 and performed a large-scale collection of Salmonella 4,[5],12:i:- with whole genome sequencing (WGS) data across the globe, to better understand the population structure, antimicrobial resistance (AMR) genomic characterization, and transmission routes of Salmonella 4,[5],12:i:- across Guangdong. Salmonella 4,[5],12:i:- strains showed broad genetic diversity; Guangdong isolates were found to be widely distributed among the global lineages. Of note, we identified the formation of a novel Guangdong clade (Bayesian analysis of population structure lineage 1 [BAPS1]) genetically diversified from the global isolates and likely emerged around 1990s. BAPS1 exhibits unique genomic features, including large pan-genome, decreased ciprofloxacin susceptibility due to mutation in gyrA and carriage of plasmid-mediated quinolone resistance (PMQR) genes, and the multidrug-resistant IncHI2 plasmid. Furthermore, high genetic similarity was found between strains collected from Guangdong, Europe, and North America, indicating the association with multiple introductions from overseas. These results suggested that global dissemination and local clonal expansion simultaneously occurred in Guangdong, China, and horizontally acquired resistance to first-line and last-line antimicrobials at local level, underlying emergences of extensive drug and pan-drug resistance. Our findings have increased the knowledge of global and local epidemics of Salmonella 4,[5],12:i:- in Guangdong, China, and provided a comprehensive baseline data set essential for future molecular surveillance.IMPORTANCESalmonella 4,[5],12:i:- has been regarded as the predominant pandemic serotype causing diarrheal diseases globally, while multidrug resistance (MDR) constitutes great public health concerns. This study provided a detailed and comprehensive genome-scale analysis of this important Salmonella serovar in the past decade in Guangdong, China. Our results revealed the complexity of two distinct transmission modes, namely global transmission and local expansion, circulating in Guangdong over a decade. Using phylogeography models, the origin of Salmonella 4,[5],12:i:- was predicted from two aspects, year and country, that is, Salmonella 4,[5],12:i:- emerged in 1983, and was introduced from the UK, and subsequently differentiated into the local endemic lineage circa 1991. Additionally, based on the pan-genome analysis, it was found that the gene accumulation rate in local endemic BAPS 1 lineage was higher than in other lineages, and the horizontal transmission of MDR IncHI2 plasmid associated with high resistance played a major role, which showed the potential threat to public health.

Whole-Genome sequencing in routine Mycobacterium bovis epidemiology - scoping the potential

Mycobacterium bovis the main agent of bovine tuberculosis (bTB), presents as a series of spatially-localised micro-epidemics across landscapes. Classical molecular typing methods applied to these micro-epidemics, based on genotyping a few variable loci, have significantly improved our understanding of potential epidemiological links between outbreaks. However, they have limited utility owing to low resolution. Conversely, whole-genome sequencing (WGS) provides the highest resolution data available for molecular epidemiology, producing richer outbreak tracing, insights into phylogeography and epidemic evolutionary history. We illustrate these advantages by focusing on a common single lineage of M. bovis (1.140) from Northern Ireland. Specifically, we investigate the spatial sub-structure of 20 years of herd-level multi locus VNTR analysis (MLVA) surveillance data and WGS data from a down sampled subset of isolates of this MLVA type over the same time frame. We mapped 2108 isolate locations of MLVA type 1.140 over the years 2000-2022. We also mapped the locations of 148 contemporary WGS isolates from this lineage, over a similar geographic range, stratifying by single nucleotide polymorphism (SNP) relatedness cut-offs of 15 SNPs. We determined a putative core range for the 1.140 MLVA type and SNP-defined sequence clusters using a 50 % kernel density estimate, using cattle movement data to inform on likely sources of WGS isolates found outside of core ranges. Finally, we applied Bayesian phylogenetic methods to investigate past population history and reproductive number of the 1.140 M. bovis lineage. We demonstrate that WGS SNP-defined clusters exhibit smaller core ranges than the established MLVA type - facilitating superior disease tracing. We also demonstrate the superior functionality of WGS data in determining how this lineage was disseminated across the landscape, likely via cattle movement and to infer how its effective population size and reproductive number has been in flux since its emergence. These initial findings highlight the potential of WGS data for routine monitoring of bTB outbreaks.

Genomic analysis reveals the population structure and antimicrobial resistance of avian Pasteurella multocida in China

OBJECTIVES

To investigate the population structure and antimicrobial resistance (AMR) of avian Pasteurella multocida in China.

METHODS

Utilizing WGS analysis, we explored the phylogeny using a dataset of 546 genomes, comprising avian P. multocida isolates from China (n = 121), the USA (n = 165), Australia(n = 153), Bangladesh (n = 3) and isolates of other hosts from China (n = 104). We examined the integrative and conjugative element (ICE) structures and the distribution of their components carrying resistance genes, and reconstructed the evolutionary history of A:L1:ST129 (n = 110).

RESULTS

The population structure of avian P. multocida in China was dominated by the A:L1:ST129 clone with limited genetic diversity. A:L1:ST129 isolates possessed a broader spectrum of resistance genes at comparatively higher frequencies than those from other hosts and countries. The novel putative ICEs harboured complex resistant clusters that were prevalent in A:L1:ST129. Bayesian analysis predicted that the A:L1:ST129 clone emerged around 1923, and evolved slowly.

CONCLUSIONS

A:L1:ST129 appears to possess a host predilection towards avian species in China, posing a potential health threat to other animals. The complex AMR determinants coupled with high frequencies may strengthen the population dominance of A:L1:ST129. The extensive antimicrobial utilization in poultry farming and the mixed rearing practices could have accelerated AMR accumulation in A:L1:ST129. ICEs, together with their resistant clusters, significantly contribute to resistance gene transfer and facilitate the adaptation of A:L1:ST129 to ecological niches. Despite the genetic stability and slow evolution rate, A:L1:ST129 deserves continued monitoring due to its propensity to retain resistance genes, warranting global attention to preclude substantial economic losses.

Probability of outbreaks and cross-border dissemination of the emerging pathogen: a genomic survey of Elizabethkingia meningoseptica

IMPORTANCE

Elizabethkingia meningoseptica is an emerging infectious agent associated with life-threatening infections in immunocompromised individuals. However, there are limited data available on the genomic features of E. meningoseptica. This study aims to characterize the geographical distribution, phylogenetic evolution, pathogenesis, and transmission of this bacterium. A systematic analysis of the E. meningoseptica genome revealed that a common ancestor of this bacterium existed 90 years ago. The evolutionary history showed no significant relationship with the sample source, origin, or region, despite the presence of genetic diversity. Whole genome sequencing data also demonstrated that E. meningoseptica bacteria possess inherent resistance and pathogenicity, enabling them to spread within the same hospital and even across borders. This study highlights the potential for E. meningoseptica to cause severe nosocomial outbreaks and horizontal transmission between countries worldwide. The available evidence is crucial for the development of evidence-based public health policies to prevent global outbreaks caused by emerging pathogens.

An integrated nationwide genomics study reveals transmission modes of typhoid fever in China

IMPORTANCE

Typhoid fever is a life-threatening disease caused by Salmonella enterica serovar Typhi, resulting in a significant disease burden across developing countries. Historically, China was very much close to the global epicenter of typhoid, but the role of typhoid transmission within China and among epicenter remains overlooked in previous investigations. By using newly produced genomics on a national scale, we clarify the complex local and global transmission history of such a notorious disease agent in China spanning the most recent five decades, which largely undermines the global public health network.

Host-specific adaptation drove the coevolution of leek yellow stripe virus and Allium plants

Host adaptation plays a crucial role in virus evolution and is a consequence of long-term interactions between virus and host in a complex arms race between host RNA silencing and viral RNA silencing suppressor (RSS) as counterdefense. Leek yellow stripe virus (LYSV), a potyvirus causing yield loss of garlic, infects several species of Allium plants. The unexpected discovery of an interspecific hybrid of garlic, leek, and great-headed (GH) garlic motivated us to explore the host-adaptive evolution of LYSV. Here, using Bayesian phylogenetic comparative methods and a functional assay of viral RSS activity, we show that the evolutionary context of LYSV has been shaped by the host adaptation of the virus during its coevolution with Allium plants. Our phylogenetic analysis revealed that LYSV isolates from leek and their taxonomic relatives (Allium ampeloprasum complex; AAC) formed a distinct monophyletic clade separate from garlic isolates and are likely to be uniquely adapted to AAC. Our comparative studies on viral accumulation indicated that LYSV accumulated at a low level in leek, whereas LYSVs were abundant in other Allium species such as garlic and its relatives. When RSS activity of the viral P1 and HC-Pro of leek LYSV isolate was analyzed, significant synergism in RSS activity between the two proteins was observed in leek but not in other species, suggesting that viral RSS activity may be important for the viral host-specific adaptation. We thus consider that LYSV may have undergone host-specific evolution at least in leek, which must be driven by speciation of its Allium hosts. IMPORTANCE Potyviruses are the most abundant plant RNA viruses and are extremely diversified in terms of their wide host range. Due to frequent host switching during their evolution, host-specific adaptation of potyviruses may have been shaped by numerous host factors. However, any critical determinants for viral host range remain largely unknown, possibly because of the repeated gain and loss of virus infectivity of plants. Leek yellow stripe virus (LYSV) is a species of the genus Potyvirus, which has a relatively narrow host range, generally limited to hosts in the genus Allium. Our investigations on leek and leek relatives (Allium ampeloprasum complex), which must have been generated through interspecies hybridization, revealed that LYSV accumulation remained low in leek as a result of viral host adaptation in competition with host resistance such as RNA silencing. This study presents LYSV as an ideal model to study the process of host-adaptive evolution and virus-host coevolution.

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).

Global molecular evolution and phylogeographic analysis of barley yellow dwarf virus based on the cp and mp genes

Barley yellow dwarf virus (BYDV) has caused considerable losses in the global production of grain crops such as wheat, barley and maize. We investigated the phylodynamics of the virus by analysing 379 and 485 nucleotide sequences of the genes encoding the coat protein and movement protein, respectively. The maximum clade credibility tree indicated that BYDV-GAV and BYDV-MAV, BYDV-PAV and BYDV-PAS share the same evolutionary lineage, respectively. The diversification of BYDV arises from its adaptability to vector insects and geography. Bayesian phylogenetic analyses showed that the mean substitution rates of the coat and movement proteins of BYDV ranged from 8.327 × 10^- 4 (4.700 × 10^- 4-1.228 × 10^- 3) and 8.671 × 10^- 4 (6.143 × 10^- 4-1.130 × 10^- 3) substitutions/site/year, respectively. The time since the most recent common BYDV ancestor was 1434 (1040-1766) CE (Common Era). The Bayesian skyline plot (BSP) showed that the BYDV population experienced dramatic expansions approximately 8 years into the 21st century, followed by a dramatic decline in less than 15 years. Our phylogeographic analysis showed that the BYDV population originating in the United States was subsequently introduced to Europe, South America, Australia and Asia. The migration pathways of BYDV suggest that the global spread of BYDV is associated with human activities.

Mitogenome evidence shows two radiation events and dispersals of matrilineal ancestry from northern coastal China to the Americas and Japan

Although it is widely recognized that the ancestors of Native Americans (NAs) primarily came from Siberia, the link between mitochondrial DNA (mtDNA) lineage D4h3a (typical of NAs) and D4h3b (found so far only in East China and Thailand) raises the possibility that the ancestral sources for early NAs were more variegated than hypothesized. Here, we analyze 216 contemporary (including 106 newly sequenced) D4h mitogenomes and 39 previously reported ancient D4h data. The results reveal two radiation events of D4h in northern coastal China, one during the Last Glacial Maximum and the other within the last deglaciation, which facilitated the dispersals of D4h sub-branches to different areas including the Americas and the Japanese archipelago. The coastal distributions of the NA (D4h3a) and Japanese lineages (D4h1a and D4h2), in combination with the Paleolithic archaeological similarities among Northern China, the Americas, and Japan, lend support to the coastal dispersal scenario of early NAs.

Genomic Epidemiology and Transmission Dynamics of Global Coxsackievirus B4

The aim of this study was to determine the global genetic diversity and transmission dynamics of coxsackievirus B4 (CVB4) and to propose future directions for disease surveillance. Next-generation sequencing was performed to obtain the complete genome sequence of CVB4, and the genetic diversity and transmission dynamics of CVB4 worldwide were analyzed using bioinformatics methods such as phylogenetic analysis, evolutionary dynamics, and phylogeographic analysis. Forty complete genomes of CVB4 were identified from asymptomatic infected individuals and hand, foot, and mouth disease (HFMD) patients. Frequent recombination between CVB4 and EV-B multiple serotypes in the 3D^pol region was found and formed 12 recombinant patterns (A-L). Among these, the CVB4 isolated from asymptomatic infected persons and HFMD patients belonged to lineages H and I, respectively. Transmission dynamics analysis based on the VP1 region revealed that CVB4 epidemics in countries outside China were dominated by the D genotype, whereas the E genotype was dominant in China, and both genotypes evolved at a rate of > 6.50 × 10^-3 substitutions/site/year. CVB4 spreads through the population unseen, with the risk of disease outbreaks persisting as susceptible individuals accumulate. Our findings add to publicly available CVB4 genomic sequence data and deepen our understanding of CVB4 molecular epidemiology.

Suspicions of two bridgehead invasions of Xylella fastidiosa subsp. multiplex in France

Of American origin, a wide diversity of Xylella fastidiosa strains belonging to different subspecies have been reported in Europe since 2013 and its discovery in Italian olive groves. Strains from the subspecies multiplex (ST6 and ST7) were first identified in France in 2015 in urban and natural areas. To trace back the most probable scenario of introduction in France, the molecular evolution rate of this subspecies was estimated at 3.2165 × 10^-7 substitutions per site per year, based on heterochronous genome sequences collected worldwide. This rate allowed the dating of the divergence between French and American strains in 1987 for ST6 and in 1971 for ST7. The development of a new VNTR-13 scheme allowed tracing the spread of the bacterium in France, hypothesizing an American origin. Our results suggest that both sequence types were initially introduced and spread in Provence-Alpes-Côte d'Azur (PACA); then they were introduced in Corsica in two waves from the PACA bridgehead populations.

Habitat connectivity and host relatedness influence virus spread across an urbanising landscape in a fragmentation-sensitive carnivore

Spatially heterogeneous landscape factors such as urbanisation can have substantial effects on the severity and spread of wildlife diseases. However, research linking patterns of pathogen transmission to landscape features remains rare. Using a combination of phylogeographic and machine learning approaches, we tested the influence of landscape and host factors on feline immunodeficiency virus (FIV_Lru) genetic variation and spread among bobcats (Lynx rufus) sampled from coastal southern California. We found evidence for increased rates of FIV_Lru lineage spread through areas of higher vegetation density. Furthermore, single-nucleotide polymorphism (SNP) variation among FIV_Lru sequences was associated with host genetic distances and geographic location, with FIV_Lru genetic discontinuities precisely correlating with known urban barriers to host dispersal. An effect of forest land cover on FIV_Lru SNP variation was likely attributable to host population structure and differences in forest land cover between different populations. Taken together, these results suggest that the spread of FIV_Lru is constrained by large-scale urban barriers to host movement. Although urbanisation at fine spatial scales did not appear to directly influence virus transmission or spread, we found evidence that viruses transmit and spread more quickly through areas containing higher proportions of natural habitat. These multiple lines of evidence demonstrate how urbanisation can change patterns of contact-dependent pathogen transmission and provide insights into how continued urban development may influence the incidence and management of wildlife disease.

Salmonella enterica serovar Cerro displays a phylogenetic structure and genomic features consistent with virulence attenuation and adaptation to cattle

Salmonella enterica subsp. enterica (S.) serovar Cerro is rarely isolated from human clinical cases of salmonellosis but represents the most common serovar isolated from cattle without clinical signs of illness in the United States. In this study, using a large, diverse set of 316 isolates, we utilized genomic methods to further elucidate the evolutionary history of S. Cerro and to identify genomic features associated with its apparent virulence attenuation in humans. Phylogenetic analyses showed that within this polyphyletic serovar, 98.4% of isolates (311/316) represent a monophyletic clade within section Typhi and the remaining 1.6% of isolates (5/316) form a monophyletic clade within subspecies enterica Clade A1. Of the section Typhi S. Cerro isolates, 93.2% of isolates (290/311) clustered into a large clonal clade comprised of predominantly sequence type (ST) 367 cattle and environmental isolates, while the remaining 6.8% of isolates (21/311), primarily from human clinical sources, clustered outside of this clonal clade. A tip-dated phylogeny of S. Cerro ST367 identified two major clades (I and II), one of which overwhelmingly consisted of cattle isolates that share a most recent common ancestor that existed circa 1975. Gene presence/absence and rarefaction curve analyses suggested that the pangenome of section Typhi S. Cerro is open, potentially reflecting the gain/loss of prophage; human isolates contained the most open pangenome, while cattle isolates had the least open pangenome. Hypothetically disrupted coding sequences (HDCs) displayed clade-specific losses of intact speC and sopA virulence genes within the large clonal S. Cerro clade, while loss of intact vgrG, araH, and vapC occurred in all section Typhi S. Cerro isolates. Further phenotypic analysis suggested that the presence of a premature stop codon in speC does not abolish ornithine decarboxylase activity in S. Cerro, likely due to the activity of the second ornithine decarboxylase encoded by speF, which remained intact in all isolates. Overall, our study identifies specific genomic features associated with S. Cerro's infrequent isolation from humans and its apparent adaptation to cattle, which has broader implications for informing our understanding of the evolutionary events facilitating host adaptation in Salmonella.

Genomic Epidemiology and Phylodynamic Analysis of Enterovirus A71 Reveal Its Transmission Dynamics in Asia

Enterovirus A71 (EV-A71) is one of the main pathogens causing hand, foot, and mouth disease (HFMD) outbreaks in Asian children under 5 years of age. In severe cases, it can cause neurological complications and be life-threatening. In this study, 200 newly sequenced EV-A71 whole-genome sequences were combined with 772 EV-A71 sequences from GenBank for large-scale analysis to investigate global EV-A71 epidemiology, phylogeny, and Bayesian phylodynamic characteristics. Based on the phylogenetic analysis of the EV-A71 3D^pol region, six new evolutionary lineages (lineages B, J, K, O, P, and Q) were found in this study, and the number of evolutionary lineages was expanded from 11 to 17. Temporal dynamics and recombination breakpoint analyses based on genotype C revealed that recombination of nonstructural protein-coding regions, including 3D^pol, is an important reason for the emergence of new lineages. The EV-A71 epidemic in the Asia-Pacific region is complex, and phylogeographic analysis found that Vietnam played a key role in the spread of subgenotypes B5 and C4. The origin of EV-A71 subgenotype C4 in China is East China, which is closely related to the prevalence of subgenotype C4 in the south and throughout China. Selection pressure analysis revealed that, in addition to VP1 amino acid residues VP1-98 and VP1-145, which are associated with EV-A71 pathogenicity, amino acid residues VP1-184 and VP1-249 were also positively selected, and their functions still need to be determined by biology and immunology. This study aimed to provide a solid theoretical basis for EV-A71-related disease surveillance and prevention, antiviral research, and vaccine development through a comprehensive analysis. IMPORTANCE EV-A71 is one of the most important pathogens causing HFMD outbreaks; however, large-scale studies of EV-A71 genomic epidemiology are currently lacking. In this study, 200 new EV-A71 whole-genome sequences were determined. Combining these with 772 EV-A71 whole-genome sequences in the GenBank database, the evolutionary and transmission characteristics of global and Asian EV-A71 were analyzed. Six new evolutionary lineages were identified in this study. We also found that recombination in nonstructural protein-coding regions, including 3D^pol, is an important cause for the emergence of new lineages. The results provided a solid theoretical basis for EV-A71-related disease surveillance and prevention, antiviral research, and vaccine development.

Population genomics of emerging Elizabethkingia anophelis pathogens reveals potential outbreak and rapid global dissemination

Elizabethkingia anophelis is an emerging species and have increasingly been reported to cause life-threatening infections and even outbreaks in humans. Nevertheless, there is little data regarding the E. anophelis geographical distribution, phylogenetic structure, and transmission across the globe, especially in Asia. We utilize whole genome sequencing (WGS) data to define a global population framework, phylogenetic structure, geographical distribution, and transmission evaluation of E. anophelis pathogens. The geographical distribution diagram revealed the emerging pathogenic bacteria already distributed in various countries worldwide, especially in the USA and China. Strikingly, phylogenetic analysis showed a part of our China original E. anophelis shared the same ancestor with the USA outbreak strain, which implies the possibility of localized outbreaks and global spread. These closer related strains also contained ICEEaI, which might insert into a disrupted DNA repair mutY gene and made the strain more liable to mutation and outbreak infection. BEAST analysis showed that the most recent common ancestor for ICEEaI E. anophelis was dated twelve years ago, and China might be the most likely recent source of this bacteria. Our study sheds light on the potential possibility of E. anophelis causing the large-scale outbreak and rapid global dissemination. Continued genomic surveillance of the dynamics of E. anophelis populations will generate further knowledge for optimizing future prevent global outbreak infections.

Transcontinental spread and evolution of Mycobacterium tuberculosis W148 European/Russian clade toward extensively drug resistant tuberculosis

Transmission-driven multi-/extensively drug resistant (M/XDR) tuberculosis (TB) is the largest single contributor to human mortality due to antimicrobial resistance. A few major clades of the Mycobacterium tuberculosis complex belonging to lineage 2, responsible for high prevalence of MDR-TB in Eurasia, show outstanding transnational distributions. Here, we determined factors underlying the emergence and epidemic spread of the W148 clade by genome sequencing and Bayesian demogenetic analyses of 720 isolates from 23 countries. We dated a common ancestor around 1963 and identified two successive epidemic expansions in the late 1980s and late 1990s, coinciding with major socio-economic changes in the post-Soviet Era. These population expansions favored accumulation of resistance mutations to up to 11 anti-TB drugs, with MDR evolving toward additional resistances to fluoroquinolones and second-line injectable drugs within 20 years on average. Timescaled haplotypic density analysis revealed that widespread acquisition of compensatory mutations was associated with transmission success of XDR strains. Virtually all W148 strains harbored a hypervirulence-associated ppe38 gene locus, and incipient recurrent emergence of prpR mutation-mediated drug tolerance was detected. The outstanding genetic arsenal of this geographically widespread M/XDR strain clade represents a “perfect storm” that jeopardizes the successful introduction of new anti-M/XDR-TB antibiotic regimens.

An outbreak of the multidrug-resistant Mycobacterium tuberculosis lineage W148 has spread widely across Russia, Central Asia and Europe. Here, the authors use whole genome sequences of ~700 isolates of this lineage collected over ~20 years to analyze its spread, evolution of drug resistance, and impact of compensatory mutations.

Phylogenomic analysis of Salmonella Indiana ST17, an emerging MDR clonal group in China

OBJECTIVES

To reconstruct the genomic epidemiology and evolution of MDR Salmonella Indiana in China.

METHODS

A total of 108 Salmonella Indiana strains were collected from humans and livestock in China. All isolates were subjected to WGS and antimicrobial susceptibility testing. Phylogenetic relationships and evolutionary analyses were conducted using WGS data from this study and the NCBI database.

RESULTS

Almost all 108 Salmonella Indiana strains displayed the MDR phenotype. Importantly, 84 isolates possessed concurrent resistance to ciprofloxacin and cefotaxime. WGS analysis revealed that class 1 integrons on the chromosome and IncHI2 plasmids were the key vectors responsible for multiple antibiotic resistance gene (ARG) [including ESBL and plasmid-mediated quinolone resistance (PMQR) genes] transmission among Salmonella Indiana. The 108 Salmonella Indiana dataset displayed a relatively large core genome and ST17 was the predominant ST. Moreover, the global ST17 Salmonella Indiana strains could be divided into five distinct lineages, each of which was significantly associated with a geographical distribution. Genomic analysis revealed multiple antimicrobial resistance determinants and QRDR mutations in Chinese lineages, which almost did not occur in other global lineages. Using molecular clock analysis, we hypothesized that ST17 isolates have existed since 1956 and underwent a major population expansion from the 1980s to the 2000s and the genetic diversity started to decrease around 2011, probably due to geographical barriers, antimicrobial selective pressure and MDR, favouring the establishment of this prevalent multiple antibiotic-resistant lineage and local epidemics.

CONCLUSIONS

This study revealed that adaptation to antimicrobial pressure was possibly pivotal in the recent evolutionary trajectory for the clonal spread of ST17 Salmonella Indiana in China.

Transmission history of SARS-CoV-2 in humans and white-tailed deer

The emergence of a novel pathogen in a susceptible population can cause rapid spread of infection. High prevalence of SARS-CoV-2 infection in white-tailed deer (Odocoileus virginianus) has been reported in multiple locations, likely resulting from several human-to-deer spillover events followed by deer-to-deer transmission. Knowledge of the risk and direction of SARS-CoV-2 transmission between humans and potential reservoir hosts is essential for effective disease control and prioritisation of interventions. Using genomic data, we reconstruct the transmission history of SARS-CoV-2 in humans and deer, estimate the case finding rate and attempt to infer relative rates of transmission between species. We found no evidence of direct or indirect transmission from deer to human. However, with an estimated case finding rate of only 4.2%, spillback to humans cannot be ruled out. The extensive transmission of SARS-CoV-2 within deer populations and the large number of unsampled cases highlights the need for active surveillance at the human–animal interface.

The Diversity and Spatiotemporally Evolutionary Dynamic of Atypical Porcine Pestivirus in China

The presence of congenital tremor (CT) type A-II in newborn piglets, caused by atypical porcine pestivirus (APPV), has been a focus since 2016. However, the source, evolutionary history, and transmission pattern of APPV in China remain poorly understood. In this study, we undertook phylogenetic analyses based on available complete E2 gene sequences along with 98 newly sequenced E2 genes between 2016 and 2020 in China within the context of global genetic diversity. The phylogenies revealed four distinct lineages of APPV, and interestingly, all lineages could be detected in China with the greatest diversity. Bayesian phylogenetic analyses showed that the E2 gene evolves at a mean rate of 1.22 × 10⁻³ (8.54 × 10⁻⁴-1.60 × 10⁻³) substitutions/site/year. The most recent common ancestor for APPVs is dated to 1886 (1837–1924) CE, somewhat earlier than the documented emergence of CT (1922 CE). Our phylogeographic analyses suggested that the APPV population possibly originated in the Netherlands, a country with developed livestock husbandry, and was introduced into China during the period 1837–2010. Guangdong, as a primary seeding population together with Central and Southwest China as epidemic linkers, was responsible for the dispersal of APPVs in China. The transmission pattern of “China lineages” (lineage 3 and lineage 4) presented a “south to north” movement tendency, which was likely associated with the implementation of strict environmental policy in China since 2000. Reconstruction of demographic history showed that APPV population size experienced multiple changes, which correlated well with the dynamic of the number of pigs in the past decades in China. Besides, positively selected pressure and geography-driven adaptation were supposed to be key factors for the diversification of APPV lineages. Our findings provide comprehensive insights into the diversity and spatiotemporal dynamic of APPV in China.

Divergent Viruses Discovered in Swine Alter the Understanding of Evolutionary History and Genetic Diversity of the Respirovirus Genus and Related Porcine Parainfluenza Viruses

Paramyxoviridae is a rapidly growing family of viruses, whose potential for cross-species transmission makes it difficult to predict the harm of newly emerging viruses to humans and animals. To better understand their diversity, evolutionary history, and co-evolution with their hosts, we analyzed a collection of porcine parainfluenza virus (PPIV) genomes to reconstruct the species classification basis and evolutionary history of the Respirovirus genus. We sequenced 17 complete genomes of porcine respirovirus 1 (also known as porcine parainfluenza virus 1; PPIV-1), thereby nearly tripling the number of currently available PPIV-1 genomes. We found that PPIV-1 was widely prevalent in China with two divergent lineages, PPIV-1a and PPIV-1b. We further provided evidence that a new species, porcine parainfluenza virus 2 (PPIV-2), had recently emerged in China. Our results pointed to a need for revising the current species demarcation criteria of the Respirovirus genus. In addition, we used PPIV-1 as an example to explore recombination and diversity of the Respirovirus genus. Interestingly, we only detected heterosubtypic recombination events between PPIV-1a and PPIV-1b with no intrasubtypic recombination events. The recombination hotspots highlighted a diverse geography-dependent genome structure of paramyxovirus infecting swine in China. Furthermore, we found no evidence of co-evolution between respirovirus and its host, indicating frequent cross-species transmission. In summary, our analyses showed that swine can be infected with a broad range of respiroviruses and recombination may serve as an important evolutionary mechanism for the Respirovirus genus’ greater diversity in genome structure than previously anticipated.

IMPORTANCE Livestock have emerged as critically underrecognized sources of paramyxovirus diversity, including pigs serving as the source of Nipah virus (NiV) and swine parainfluenza virus type 3, and goats and bovines harboring highly divergent viral lineages. Here, we identified a new species of Respirovirus genus named PPIV-2 in swine and proposed to revise the species demarcation criteria of the Respirovirus genus. We found heterosubtypic recombination events and high genetic diversity in PPIV-1. Further, we showed that genetic recombination may have occurred in the Respirovirus genus which may be associated with host range expansion. The continued expansion of Respirovirus genus diversity in livestock with relatively high human contact rates requires enhanced surveillance and ongoing evaluation of emerging cross-species transmission threats.

Genetic evidence for post-glacial expansion from a southern refugium in the eastern moa (Emeus crassus)

Cycles of glacial expansion and contraction throughout the Pleistocene drove increases and decreases, respectively, in the geographical range and population size of many animal species. Genetic data have revealed that during glacial maxima the distribution of many Eurasian animals was restricted to small refugial areas, from which species expanded to reoccupy parts of their former range as the climate warmed. It has been suggested that the extinct eastern moa (Emeus crassus)-a large, flightless bird from New Zealand-behaved analogously during glacial maxima, possibly surviving only in a restricted area of lowland habitat in the southern South Island of New Zealand during the Last Glacial Maximum (LGM). However, previous studies have lacked the power and geographical sampling to explicitly test this hypothesis using genetic data. Here we analyse 46 ancient mitochondrial genomes from Late Pleistocene and Holocene bones of the eastern moa from across their post-LGM distribution. Our results are consistent with a post-LGM increase in the population size and genetic diversity of eastern moa. We also demonstrate that genetic diversity was higher in eastern moa from the southern extent of their range, supporting the hypothesis that they expanded from a single glacial refugium following the LGM.

Bayesian Phylodynamic Analysis Reveals the Dispersal Patterns of African Swine Fever Virus

The evolutionary and demographic history of African swine fever virus (ASFV) is potentially quite valuable for developing efficient and sustainable management strategies. In this study, we performed phylogenetic, phylodynamic, and phylogeographic analyses of worldwide ASFV based on complete ASFV genomes, B646L gene, and E183L gene sequences obtained from NCBI to understand the epidemiology of ASFV. Bayesian phylodynamic analysis and phylogenetic analysis showed highly similar results of group clustering between E183L and the complete genome. The evidence of migration and the demographic history of ASFV were also revealed by the Bayesian phylodynamic analysis. The evolutionary rate was estimated to be 1.14 × 10⁻⁵ substitution/site/year. The large out-migration from the viral population in South Africa played a crucial role in spreading the virus worldwide. Our study not only provides resources for the better utilization of genomic data but also reveals the comprehensive worldwide evolutionary history of ASFV with a broad sampling window across ~70 years. The characteristics of the virus spatiotemporal transmission are also elucidated, which could be of great importance for devising strategies to control the virus.

Recent Zoonotic Spillover and Tropism Shift of a Canine Coronavirus Is Associated with Relaxed Selection and Putative Loss of Function in NTD Subdomain of Spike Protein

A canine coronavirus (CCoV) has now been reported from two independent human samples from Malaysia (respiratory, collected in 2017–2018; CCoV-HuPn-2018) and Haiti (urine, collected in 2017); these two viruses were nearly genetically identical. In an effort to identify any novel adaptations associated with this apparent shift in tropism we carried out detailed evolutionary analyses of the spike gene of this virus in the context of related Alphacoronavirus 1 species. The spike 0-domain retains homology to CCoV2b (enteric infections) and Transmissible Gastroenteritis Virus (TGEV; enteric and respiratory). This domain is subject to relaxed selection pressure and an increased rate of molecular evolution. It contains unique amino acid substitutions, including within a region important for sialic acid binding and pathogenesis in TGEV. Overall, the spike gene is extensively recombinant, with a feline coronavirus type II strain serving a prominent role in the recombinant history of the virus. Molecular divergence time for a segment of the gene where temporal signal could be determined, was estimated at around 60 years ago. We hypothesize that the virus had an enteric origin, but that it may be losing that particular tropism, possibly because of mutations in the sialic acid binding region of the spike 0-domain.

Global phylogeny of Treponema pallidum lineages reveals recent expansion and spread of contemporary syphilis

Syphilis, which is caused by the sexually transmitted bacterium Treponema pallidum subsp. pallidum, has an estimated 6.3 million cases worldwide per annum. In the past ten years, the incidence of syphilis has increased by more than 150% in some high-income countries, but the evolution and epidemiology of the epidemic are poorly understood. To characterize the global population structure of T. pallidum, we assembled a geographically and temporally diverse collection of 726 genomes from 626 clinical and 100 laboratory samples collected in 23 countries. We applied phylogenetic analyses and clustering, and found that the global syphilis population comprises just two deeply branching lineages, Nichols and SS14. Both lineages are currently circulating in 12 of the 23 countries sampled. We subdivided T. p. pallidum into 17 distinct sublineages to provide further phylodynamic resolution. Importantly, two Nichols sublineages have expanded clonally across 9 countries contemporaneously with SS14. Moreover, pairwise genome analyses revealed examples of isolates collected within the last 20 years from 14 different countries that had genetically identical core genomes, which might indicate frequent exchange through international transmission. It is striking that most samples collected before 1983 are phylogenetically distinct from more recently isolated sublineages. Using Bayesian temporal analysis, we detected a population bottleneck occurring during the late 1990s, followed by rapid population expansion in the 2000s that was driven by the dominant T. pallidum sublineages circulating today. This expansion may be linked to changing epidemiology, immune evasion or fitness under antimicrobial selection pressure, since many of the contemporary syphilis lineages we have characterized are resistant to macrolides.

Inferring Mycobacterium bovis transmission between cattle and badgers using isolates from the Randomised Badger Culling Trial

Mycobacterium bovis (M. bovis) is a causative agent of bovine tuberculosis, a significant source of morbidity and mortality in the global cattle industry. The Randomised Badger Culling Trial was a field experiment carried out between 1998 and 2005 in the South West of England. As part of this trial, M. bovis isolates were collected from contemporaneous and overlapping populations of badgers and cattle within ten defined trial areas. We combined whole genome sequences from 1,442 isolates with location and cattle movement data, identifying transmission clusters and inferred rates and routes of transmission of M. bovis. Most trial areas contained a single transmission cluster that had been established shortly before sampling, often contemporaneous with the expansion of bovine tuberculosis in the 1980s. The estimated rate of transmission from badger to cattle was approximately two times higher than from cattle to badger, and the rate of within-species transmission considerably exceeded these for both species. We identified long distance transmission events linked to cattle movement, recurrence of herd breakdown by infection within the same transmission clusters and superspreader events driven by cattle but not badgers. Overall, our data suggests that the transmission clusters in different parts of South West England that are still evident today were established by long-distance seeding events involving cattle movement, not by recrudescence from a long-established wildlife reservoir. Clusters are maintained primarily by within-species transmission, with less frequent spill-over both from badger to cattle and cattle to badger.

Contribution of historical herbarium small RNAs to the reconstruction of a cassava mosaic geminivirus evolutionary history

Emerging viral diseases of plants are recognised as a growing threat to global food security. However, little is known about the evolutionary processes and ecological factors underlying the emergence and success of viruses that have caused past epidemics. With technological advances in the field of ancient genomics, it is now possible to sequence historical genomes to provide a better understanding of viral plant disease emergence and pathogen evolutionary history. In this context, herbarium specimens represent a valuable source of dated and preserved material. We report here the first historical genome of a crop pathogen DNA virus, a 90-year-old African cassava mosaic virus (ACMV), reconstructed from small RNA sequences bearing hallmarks of small interfering RNAs. Relative to tip-calibrated dating inferences using only modern data, those performed with the historical genome yielded both molecular evolution rate estimates that were significantly lower, and lineage divergence times that were significantly older. Crucially, divergence times estimated without the historical genome appeared in discordance with both historical disease reports and the existence of the historical genome itself. In conclusion, our study reports an updated time-frame for the history and evolution of ACMV and illustrates how the study of crop viral diseases could benefit from natural history collections.

Genomic analysis of the brassica pathogen turnip mosaic potyvirus reveals its spread along the former trade routes of the Silk Road

Our article presents the most comprehensive reconstruction of the evolutionary and phylogeographic history of a major plant pathogen of brassica vegetables in Eurasia. Sampling across such a large landmass poses considerable challenges, and our study attempts to describe the spatial and temporal patterns of migration for a plant pathogen on a large scale. Our phylogeographic and molecular clock analyses show that the migration pathways of turnip mosaic potyvirus retrace some of the historical trade arteries of the Silk Road. This study demonstrates how a comprehensive genetic analysis can provide a large-scale view of the epidemiology and human-mediated spread of a plant pathogen across centuries of evolutionary history.

Plant pathogens have agricultural impacts on a global scale and resolving the timing and route of their spread can aid crop protection and inform control strategies. However, the evolutionary and phylogeographic history of plant pathogens in Eurasia remains largely unknown because of the difficulties in sampling across such a large landmass. Here, we show that turnip mosaic potyvirus (TuMV), a significant pathogen of brassica crops, spread from west to east across Eurasia from about the 17th century CE. We used a Bayesian phylogenetic approach to analyze 579 whole genome sequences and up to 713 partial sequences of TuMV, including 122 previously unknown genome sequences from isolates that we collected over the past five decades. Our phylogeographic and molecular clock analyses showed that TuMV isolates of the Asian-Brassica/Raphanus (BR) and basal-BR groups and world-Brassica3 (B3) subgroup spread from the center of emergence to the rest of Eurasia in relation to the host plants grown in each country. The migration pathways of TuMV have retraced some of the major historical trade arteries in Eurasia, a network that formed the Silk Road, and the regional variation of the virus is partly characterized by different type patterns of recombinants. Our study presents a complex and detailed picture of the timescale and major transmission routes of an important plant pathogen.

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.

Evolution and genomic insight into methicillin-resistant Staphylococcus aureus ST9 in China

OBJECTIVES

To reconstruct the evolutionary history and genomic epidemiology of Staphylococcus aureus ST9 in China.

METHODS

Using WGS analysis, we described the phylogeny of 131 S. aureus ST9 isolates collected between 2002 and 2016 from 11 provinces in China, including six clinical samples from Taiwan. We also investigated the complex structure and distribution of the lsa(E)-carrying multiresistance gene cluster, and genotyped prophages in the genomes of the ST9 isolates.

RESULTS

ST9 was subdivided into one major (n = 122) and one minor (n = 9) clade. Bayesian phylogeny predicted the divergence of ST9 isolates in pig farming in China as early as 1987, which then evolved rapidly in the following three decades. ST9 isolates shared similar multiresistance properties, which were likely acquired before the ST9 emergence in China. The accessory genome is highly conserved, and ST9 harboured similar sets of phages, but lacked certain virulence genes.

CONCLUSIONS

Host exchange and regional transmission of ST9 have occurred between pigs and humans. Pig rearing and trading might have favoured gene exchanges between ST9 isolates. Resistance genes, obtained from the environment and other isolates, were stably integrated into the chromosomal DNA. The abundance of resistance genes among ST9 is likely attributed to the extensive use of antimicrobial agents in livestock. Phages are present in the genomes of ST9 and may play a role in the rapid evolution of this ST. Although human ST9 infections are rare, ST9 isolates may constitute a potential risk to public health as a repository of antimicrobial resistance genes.

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.

Genomic and temporal analyses of Mycobacterium bovis in southern Brazil

Mycobacterium bovis is a causal agent of bovine tuberculosis (bTB), one of the most important diseases currently facing the cattle industry worldwide. Tracing the source of M. bovis infections of livestock is an important tool for understanding the epidemiology of bTB and defining control/eradication strategies. In this study, whole genome sequencing (WGS) of 74 M. bovis isolates sourced from naturally infected cattle in the State of Rio Grande do Sul (RS), southern Brazil, was used to evaluate the population structure of M. bovis in the region, identify potential transmission events and date the introduction of clonal complex (CC) European 2 (Eu2). In silico spoligotyping identified 11 distinct patterns including four new profiles and two CCs, European 1 (Eu1) and Eu2. The analyses revealed a high level of genetic diversity in the majority of herds and identified putative transmission clusters that suggested that within- and between-herd transmission is occurring in RS. In addition, a comparison with other published M. bovis isolates from Argentina, Brazil, Paraguay and Uruguay demonstrated some evidence for a possible cross-border transmission of CC Eu1 into RS from Uruguay or Argentina. An estimated date for the introduction of CC Eu2 into RS in the middle of the 19th century correlated with the historical introduction of cattle into RS to improve existing local breeds. These findings contribute to the understanding of the population structure of M. bovis in southern Brazil and highlight the potential of WGS in surveillance and helping to identify bTB transmission.

Genetic and evolutionary analysis of emerging HoBi-like pestivirus

HoBi-like pestivirus, an emerging species within the Pestivirus genus, is an important pathogen associated with a variety of clinical manifestations of ruminants, especially cattle. HoBi-like pestiviruses were identified in several countries and from different hosts, and raised concerns with regard to their acute and persistent infections, which is implicated in economic losses for cattle farmers. However, the transmission path, codon usage bias, and host adaptation of the virus has not been studied. Hence, we performed the analysis the spatio-temporal transmission based on the available 5'-UTR sequences of HoBi-like pestivirus, and then conducted codon analysis of the complete coding sequence of the virus. The results show the virus appeared in 1952 (95% HPD: 1905-1985) and may have originated in India. In addition, Italy is the hub for the spread of the virus. Moreover, six potential recombination events and two complex recombination events were discovered. Analysis of codon usage patterns revealed that the effective number of codon (ENC) values with an average of 50.85, and the codon usage bias is greatly affected by natural selection, which is different from the previous BVDV-1, 2. Finally, codon adaptation index (CAI) analysis shows that pigs may be the potential origin species of the HoBi-like pestivirus. These findings will contribute to more effective control of the spread of the virus, extend the knowledge about the genetic and evolutionary features of HoBi-like viruses and provide some information for vaccine research.

Exploring the phylogeography and population dynamics of the giant deer (Megaloceros giganteus) using Late Quaternary mitogenomes

Late Quaternary climatic fluctuations in the Northern Hemisphere had drastic effects on large mammal species, leading to the extinction of a substantial number of them. The giant deer (Megaloceros giganteus) was one of the species that became extinct in the Holocene, around 7660 calendar years before present. In the Late Pleistocene, the species ranged from western Europe to central Asia. However, during the Holocene, its range contracted to eastern Europe and western Siberia, where the last populations of the species occurred. Here, we generated 35 Late Pleistocene and Holocene giant deer mitogenomes to explore the genetics of the demise of this iconic species. Bayesian phylogenetic analyses of the mitogenomes suggested five main clades for the species: three pre-Last Glacial Maximum clades that did not appear in the post-Last Glacial Maximum genetic pool, and two clades that showed continuity into the Holocene. Our study also identified a decrease in genetic diversity starting in Marine Isotope Stage 3 and accelerating during the Last Glacial Maximum. This reduction in genetic diversity during the Last Glacial Maximum, coupled with a major contraction of fossil occurrences, suggests that climate was a major driver in the dynamics of the giant deer.

Population genomics provides insights into the evolution and adaptation to humans of the waterborne pathogen Mycobacterium kansasii

Mycobacterium kansasii can cause serious pulmonary disease. It belongs to a group of closely-related species of non-tuberculous mycobacteria known as the M. kansasii complex (MKC). Here, we report a population genomics analysis of 358 MKC isolates from worldwide water and clinical sources. We find that recombination, likely mediated by distributive conjugative transfer, has contributed to speciation and on-going diversification of the MKC. Our analyses support municipal water as a main source of MKC infections. Furthermore, nearly 80% of the MKC infections are due to closely-related M. kansasii strains, forming a main cluster that apparently originated in the 1900s and subsequently expanded globally. Bioinformatic analyses indicate that several genes involved in metabolism (e.g., maintenance of the methylcitrate cycle), ESX-I secretion, metal ion homeostasis and cell surface remodelling may have contributed to M. kansasii's success and its ongoing adaptation to the human host.

Population structure and transmission of Mycobacterium bovis in Ethiopia

Bovine tuberculosis (bTB) is endemic in cattle in Ethiopia, a country that hosts the largest national cattle herd in Africa. The intensive dairy sector, most of which is peri-urban, has the highest prevalence of disease. Previous studies in Ethiopia have demonstrated that the main cause is Mycobacterium bovis, which has been investigated using conventional molecular tools including deletion typing, spoligotyping and Mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR). Here we use whole-genome sequencing to examine the population structure of M. bovis in Ethiopia. A total of 134 M. bovis isolates were sequenced including 128 genomes from 85 mainly dairy cattle and six genomes isolated from humans, originating from 12 study sites across Ethiopia. These genomes provided a good representation of the previously described population structure of M. bovis, based on spoligotyping and demonstrated that the population is dominated by the clonal complexes African 2 (Af2) and European 3 (Eu3). A range of within-host diversity was observed amongst the isolates and evidence was found for both short- and long-distance transmission. Detailed analysis of available genomes from the Eu3 clonal complex combined with previously published genomes revealed two distinct introductions of this clonal complex into Ethiopia between 1950 and 1987, likely from Europe. This work is important to help better understand bTB transmission in cattle in Ethiopia and can potentially inform national strategies for bTB control in Ethiopia and beyond.

Molecular Evolution of Infectious Pancreatic Necrosis Virus in China

Passive virus surveillance was performed in twenty-nine salmon and trout farms from seven provinces and districts in China during the period 2017–2020. A total of 25 infectious pancreatic necrosis virus (IPNV) isolates were obtained, mainly from rainbow trout (Oncorhynchus mykiss). The molecular evolution of these Chinese IPNV isolates and the previously reported Chinese IPNV strains ChRtm213 and WZ2016 was analyzed, based on their VP2 gene coding region sequences (CDS). All 27 Chinese IPNV isolates clustered within genogroups I and V, with 24 of the IPNV isolates belonging to genogroup I (including ChRtm213 and WZ2016), and only three isolates clustering in genogroup V. The Chinese genogroup I IPNV isolates lacked diversity, composing six haplotypes with 41 polymorphic sites, and the identity of nucleotide and amino acid sequences among the entire VP2 gene CDS from these isolates was 97.44%–100% and 98.19%–100%, respectively. Divergence time analyses revealed that the Chinese genogroup I IPNV isolates likely diverged from Japanese IPNV isolates in 1985 (95% highest posterior density (HPD), 1965–1997), and diverged again in 2006 (95% HPD, 1996–2013) in China. Each of the three Chinese genogroup V IPNV isolates has a unique VP2 gene CDS, with a total of 21 polymorphic sites; the identity of nucleotide and amino acid sequences among all VP2 gene CDS from these isolates was 98.5%–99.5% and 98.6%–99.0%, respectively. The data demonstrate that genogroups I and V are more likely the currently prevalent Chinese IPNV genotypes.

Evolutionary analysis of SARS-CoV-2 spike protein for its different clades

The spike protein of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) has become the main target for antiviral and vaccine development. Despite its relevance, e information is scarse about its evolutionary traces. The aim of this study was to investigate the diversification patterns of the spike for each clade of SARS‐CoV‐2 through different approaches. Two thousand and one hundred sequences representing the seven clades of the SARS‐CoV‐2 were included. Patterns of genetic diversifications and nucleotide evolutionary rate were estimated for the spike genomic region. The haplotype networks showed a star shape, where multiple haplotypes with few nucleotide differences diverge from a common ancestor. Four hundred seventy‐nine different haplotypes were defined in the seven analyzed clades. The main haplotype, named Hap‐1, was the most frequent for clades G (54%), GH (54%), and GR (56%) and a different haplotype (named Hap‐252) was the most important for clades L (63.3%), O (39.7%), S (51.7%), and V (70%). The evolutionary rate for the spike protein was estimated as 1.08 × 10⁻³ nucleotide substitutions/site/year. Moreover, the nucleotide evolutionary rate after nine months of the pandemic was similar for each clade. In conclusion, the present evolutionary analysis is relevant as the spike protein of SARS‐CoV‐2 is the target for most therapeutic candidates; besides, changes in this protein could have consequences on viral transmission, response to antivirals and efficacy of vaccines. Moreover, the evolutionary characterization of clades improves knowledge of SARS‐CoV‐2 and deserves to be assessed in more detail as re‐infection by different phylogenetic clades has been reported.

Phylogenetic and genetic characterization of Treponema pallidum strains from syphilis patients in Japan by whole-genome sequence analysis from global perspectives

Japan has had a substantial increase in syphilis cases since 2013. However, research on the genomic features of the Treponema pallidum subspecies pallidum (TPA) strains from these cases has been limited. Here, we elucidated the genetic variations and relationships between TPA strains in Japan (detected between 2014 and 2018) and other countries by whole-genome sequencing and phylogenetic analyses, including syphilis epidemiological surveillance data and information on patient sexual orientation. Seventeen of the 20 strains in Japan were SS14- and the remaining 3 were Nichols-lineage. Sixteen of the 17 SS14-lineage strains were classified into previously reported Sub-lineage 1B. Sub-lineage 1B strains in Japan have formed distinct sub-clusters of strains from heterosexuals and strains from men who have sex with men. These strains were closely related to reported TPA strains in China, forming an East-Asian cluster. However, those strains in these countries evolved independently after diverging from their most recent common ancestor and expanded their genetic diversity during the time of syphilis outbreak in each country. The genetic difference between the TPA strains in these countries was characterized by single-nucleotide-polymorphism analyses of their penicillin binding protein genes. Taken together, our results elucidated the detailed phylogenetic features and transmission networks of syphilis.

The evolutionary dynamics of endemic human coronaviruses

Community protective immunity can affect RNA virus evolution by selecting for new antigenic variants on the scale of years, exemplified by the need of annual evaluation of influenza vaccines. The extent to which this process termed antigenic drift affects coronaviruses remains unknown. Alike the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), seasonal human coronaviruses (HCoV) likely emerged from animal reservoirs as new human pathogens in the past. We therefore analyzed the long-term evolutionary dynamics of the ubiquitous HCoV-229E and HCoV-OC43 in comparison with human influenza A virus (IAV) subtype H3N2. We focus on viral glycoprotein genes that mediate viral entry into cells and are major targets of host neutralizing antibody responses. Maximum likelihood and Bayesian phylogenies of publicly available gene datasets representing about three decades of HCoV and IAV evolution showed that all viruses had similar ladder-like tree shapes compatible with antigenic drift, supported by different tree shape statistics. Evolutionary rates inferred in a Bayesian framework were 6.5 × 10-4 (95% highest posterior density (HPD), 5.4-7.5 × 10-4) substitutions per site per year (s/s/y) for HCoV-229E spike (S) genes and 5.7 × 10-4 (95% HPD, 5-6.5 × 10-4) s/s/y for HCoV-OC43 S genes, which were about fourfold lower than the 2.5 × 10-3 (95% HPD, 2.3-2.7 × 10-3) s/s/y rate for IAV hemagglutinin (HA) genes. Coronavirus S genes accumulated about threefold less (P < 0.001) non-synonymous mutations (dN) over time than IAV HA genes. In both IAV and HCoV, the average rate of dN within the receptor binding domains (RBD) was about fivefold higher (P < 0.0001) than in other glycoprotein gene regions. Similarly, most sites showing evidence for positive selection occurred within the RBD (HCoV-229E, 6/14 sites, P < 0.05; HCoV-OC43, 23/38 sites, P < 0.01; IAV, 13/15 sites, P = 0.08). In sum, the evolutionary dynamics of HCoV and IAV showed several similarities, yet amino acid changes potentially representing antigenic drift occurred on a lower scale in endemic HCoV compared to IAV. It seems likely that pandemic SARS-CoV-2 evolution will bear similarities with IAV evolution including accumulation of adaptive changes in the RBD, requiring vaccines to be updated regularly, whereas higher SARS-CoV-2 evolutionary stability resembling endemic HCoV can be expected in the post-pandemic stage.

Time-calibrated genomic evolution of a monomorphic bacterium during its establishment as an endemic crop pathogen

Horizontal gene transfer is of major evolutionary importance as it allows for the redistribution of phenotypically important genes among lineages. Such genes with essential functions include those involved in resistance to antimicrobial compounds and virulence factors in pathogenic bacteria. Understanding gene turnover at microevolutionary scales is critical to assess the pace of this evolutionary process. Here, we characterized and quantified gene turnover for the epidemic lineage of a bacterial plant pathogen of major agricultural importance worldwide. Relying on a dense geographic sampling spanning 39 years of evolution, we estimated both the dynamics of single nucleotide polymorphism accumulation and gene content turnover. We identified extensive gene content variation among lineages even at the smallest phylogenetic and geographic scales. Gene turnover rate exceeded nucleotide substitution rate by three orders of magnitude. Accessory genes were found preferentially located on plasmids, but we identified a highly plastic chromosomal region hosting ecologically important genes such as transcription activator-like effectors. Whereas most changes in the gene content are probably transient, the rapid spread of a mobile element conferring resistance to copper compounds widely used for the management of plant bacterial pathogens illustrates how some accessory genes can become ubiquitous within a population over short timeframes.

Genomic epidemiology of coxsackievirus A16 in mainland of China, 2000-18

Hand, foot, and mouth disease (HFMD), which is a frequently reported and concerning disease worldwide, is a severe burden on societies globally, especially in the countries of East and Southeast Asia. Coxsackievirus A16 (CV-A16) is one of the most important causes of HFMD and a severe threat to human health, especially in children under 5 years of age. To investigate the epidemiological characteristics, spread dynamics, recombinant forms (RFs), and other features of CV-A16, we leveraged the continuous surveillance data of CV-A16-related HFMD cases collected over an 18-year period. With the advent of the EV-A71 vaccine since 2016, which targeted the EV-A71-related HFMD cases, EV-A71-related HFMD cases decreased dramatically, whereas the CV-A16-related HFMD cases showed an upward trend from 2017 to October 2019. The CV-A16 strains observed in this study were genetically related and widely distributed in the mainland of China. Our results show that three clusters (B1a-B1c) existed in the mainland of China and that the cluster of B1b dominates the diffusion of CV-A16 in China. We found that eastern China played a decisive role in seeding the diffusion of CV-A16 in China, with a more complex and variant transmission trend. Although EV-A71 vaccine was launched in China in 2016, it did not affect the genetic diversity of CV-A16, and its genetic diversity did not decline, which confirmed the epidemiological surveillance trend of CV-A16. Two discontinuous clusters (2000-13 and 2014-18) were observed in the full-length genome and arranged along the time gradient, which revealed the reason why the relative genetic diversity of CV-A16 increased and experienced more complex fluctuation model after 2014. In addition, the switch from RFs B (RF-B) and RF-C co-circulation to RF-D contributes to the prevalence of B1b cluster in China after 2008. The correlation between genotype and RFs partially explained the current prevalence of B1b. This study provides unprecedented full-length genomic sequences of CV-A16 in China, with a wider geographic distribution and a long-term time scale. The study presents valuable information about CV-A16, aimed at developing effective control strategies, as well as a call for a more robust surveillance system, especially in the Asia-Pacific region.

Phylogenetic analysis of SARS-CoV-2 in the first few months since its emergence

During the first few months of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) evolution in a new host, contrasting hypotheses have been proposed about the way the virus has evolved and diversified worldwide. The aim of this study was to perform a comprehensive evolutionary analysis to describe the human outbreak and the evolutionary rate of different genomic regions of SARS-CoV-2. The molecular evolution in nine genomic regions of SARS-CoV-2 was analyzed using three different approaches: phylogenetic signal assessment, emergence of amino acid substitutions, and Bayesian evolutionary rate estimation in eight successive fortnights since the virus emergence. All observed phylogenetic signals were very low and tree topologies were in agreement with those signals. However, after 4 months of evolution, it was possible to identify regions revealing an incipient viral lineage formation, despite the low phylogenetic signal since fortnight 3. Finally, the SARS-CoV-2 evolutionary rate for regions nsp3 and S, the ones presenting greater variability, was estimated as 1.37 × 10^-3 and 2.19 × 10^-3 substitution/site/year, respectively. In conclusion, results from this study about the variable diversity of crucial viral regions and determination of the evolutionary rate are consequently decisive to understand essential features of viral emergence. In turn, findings may allow the first-time characterization of the evolutionary rate of S protein, crucial for vaccine development.

Does the virus cross the road? Viral phylogeographic patterns among bobcat populations reflect a history of urban development

Urban development has major impacts on connectivity among wildlife populations and is thus likely an important factor shaping pathogen transmission in wildlife. However, most investigations of wildlife diseases in urban areas focus on prevalence and infection risk rather than potential effects of urbanization on transmission itself. Feline immunodeficiency virus (FIV) is a directly transmitted retrovirus that infects many felid species and can be used as a model for studying pathogen transmission at landscape scales. We investigated phylogenetic relationships among FIV isolates sampled from five bobcat (Lynx rufus) populations in coastal southern California that appear isolated due to major highways and dense urban development. Divergence dates among FIV phylogenetic lineages in several cases reflected historical urban growth and construction of major highways. We found strong FIV phylogeographic structure among three host populations north-west of Los Angeles, largely coincident with host genetic structure. In contrast, relatively little FIV phylogeographic structure existed among two genetically distinct host populations south-east of Los Angeles. Rates of FIV transfer among host populations did not vary significantly, with the lack of phylogenetic structure south-east of Los Angeles unlikely to reflect frequent contemporary transmission among populations. Our results indicate that major barriers to host gene flow can also act as barriers to pathogen spread, suggesting potentially reduced susceptibility of fragmented populations to novel directly transmitted pathogens. Infrequent exchange of FIV among host populations suggests that populations would best be managed as distinct units in the event of a severe disease outbreak. Phylogeographic inference of pathogen transmission is useful for estimating the ability of geographic barriers to constrain disease spread and can provide insights into contemporary and historical drivers of host population connectivity.

Emergence and evolution of highly pathogenic porcine epidemic diarrhea virus by natural recombination of a low pathogenic vaccine isolate and a highly pathogenic strain in the spike gene

Outbreaks of a new variant of porcine epidemic diarrhea virus (PEDV) at the end of 2010 have raised interest in the mutation and recombination of PEDV. A PEDV strain (CN/Liaoning25/2018) isolated from a clinical outbreak of piglet diarrhea contained a 49-bp deletion in the ORF3 gene. This deletion is considered a genetic characteristic of low pathogenic attenuated vaccine strains. However, CN/Liaoning25/2018 was highly pathogenic. Complete genome sequencing, identity analysis, phylogenetic tree construction, and recombination analysis showed that this virus was a recombinant strain containing the Spike (S) gene from the highly pathogenic CN/GDZQ/2014 strain and the remaining genomic regions from the low pathogenic vaccine isolate SQ2014. Histopathology and immunohistochemistry results confirmed that this strain was highly pathogenic and indicated that intestinal epithelial cell vacuolation was positively correlated with the intensity and density of PEDV antigens. A new natural recombination model for PEDV was identified. Our results suggest that new highly pathogenic recombinant strains in the field may be generated by recombination between low pathogenic attenuated live PEDV vaccines and pathogenic circulating PEDV strains. Our findings also highlight that the 49-bp deletion of the ORF3 gene in low pathogenic attenuated vaccine strains will no longer be a reliable standard to differentiate the classical vaccine attenuated from the field strains.

A Large-Scale Outbreak of Echovirus 30 in Gansu Province of China in 2015 and Its Phylodynamic Characterization

Background

Echovirus 30 (E-30) has been investigated and reported worldwide and is closely associated with several infectious diseases, including encephalitis; myocarditis; and hand, foot, and mouth disease. Although many E-30 outbreaks associated with encephalitis have been reported around the world, it was not reported in northwest China until 2015.

Methods

The clinical samples, including the feces, serum, throat swabs, and cerebrospinal fluid, were collected for this study and were analyzed for diagnosis. E-30 was isolated and processed according to the standard procedures. The epidemiological and phylogenetic analysis were performed to indicate the characteristics of E-30 outbreaks and phylodynamics of E-30 in China.

Results

The E-30 outbreaks affected nine towns of Gansu Province in 2015, starting at a school of Nancha town and spreading to other towns within 1 month. The epidemiological features showed that children aged 6–15 years were more susceptible to E-30 infection. The genotypes B and C cocirculated in the world, whereas the latter dominated the circulation of E-30 in China. The genome sequences of this outbreak present 99.3–100% similarity among these strains, indicating a genetic-linked aggregate outbreak of E-30 in this study. Two larger genetic diversity expansions and three small fluctuations of E-30 were observed from 1987 to 2016 in China, which revealed the oscillating patterns of E-30 in China. In addition, the coastal provinces of China, such as Zhejiang, Fujian, and Shandong, were initially infected, followed by other parts of the country. The E-30 strains isolated from mainland of China may have originated from Taiwan of China in the last century.

Conclusion

The highly similar E-30 genomes in this outbreak showed an aggregate outbreak of E-30, with nine towns affected. Our results suggested that, although the genetic diversity of E-30 oscillates, the dominant lineages of E-30 in China has complex genetic transmission. The coastal provinces played an important role in E-30 spread, which implied further development of effective countermeasures. This study provides a further insight into the E-30 outbreak and transmission and illustrates the importance of valuable surveillance in the future.

The Gene Flow Direction of Geographically Distinct Phytophthora infestans Populations in China Corresponds With the Route of Seed Potato Exchange

Phytophthora infestans is a widespread destructive plant pathogen that causes economic losses worldwide to potato production. In this study, we sequenced four mitochondrial DNA gene sequences of 101 P. infestans isolates from five potato-growing regions in China to investigate the population structure and dispersal pattern of this pathogen. The concatenated mtDNA sequences in the populations showed high haplotype diversity, but low nucleotide diversity. Although there was a degree of spatial structure, our phylogeographic analyses support frequent gene flow between populations and the direction of gene flow, primarily from north to south, corresponds to the route of seed potato transportation, suggesting a role of human activities in the dispersal of P. infestans in China.

Emergence of a Novel Salmonella enterica Serotype Reading Clonal Group Is Linked to Its Expansion in Commercial Turkey Production, Resulting in Unanticipated Human Illness in North America

Two separate human outbreaks of Salmonella enterica serotype Reading occurred between 2017 and 2019 in the United States and Canada, and both outbreaks were linked to the consumption of raw turkey products. In this study, a comprehensive genomic investigation was conducted to reconstruct the evolutionary history of S. Reading from turkeys and to determine the genomic context of outbreaks involving this infrequently isolated Salmonella serotype. A total of 988 isolates of U.S. origin were examined using whole-genome-based approaches, including current and historical isolates from humans, meat, and live food animals. Broadly, isolates clustered into three major clades, with one apparently highly adapted turkey clade. Within the turkey clade, isolates clustered into three subclades, including an "emergent" clade that contained only isolates dated 2016 or later, with many of the isolates from these outbreaks. Genomic differences were identified between emergent and other turkey subclades, suggesting that the apparent success of currently circulating subclades is, in part, attributable to plasmid acquisitions conferring antimicrobial resistance, gain of phage-like sequences with cargo virulence factors, and mutations in systems that may be involved in beta-glucuronidase activity and resistance towards colicins. U.S. and Canadian outbreak isolates were found interspersed throughout the emergent subclade and the other circulating subclade. The emergence of a novel S Reading turkey subclade, coinciding temporally with expansion in commercial turkey production and with U.S. and Canadian human outbreaks, indicates that emergent strains with higher potential for niche success were likely vertically transferred and rapidly disseminated from a common source.IMPORTANCE Increasingly, outbreak investigations involving foodborne pathogens are difficult due to the interconnectedness of food animal production and distribution, and homogeneous nature of industry integration, necessitating high-resolution genomic investigations to determine their basis. Fortunately, surveillance and whole-genome sequencing, combined with the public availability of these data, enable comprehensive queries to determine underlying causes of such outbreaks. Utilizing this pipeline, it was determined that a novel clone of Salmonella Reading has emerged that coincided with increased abundance in raw turkey products and two outbreaks of human illness in North America. The rapid dissemination of this highly adapted and conserved clone indicates that it was likely obtained from a common source and rapidly disseminated across turkey production. Key genomic changes may have contributed to its apparent continued success in commercial turkeys and ability to cause illness in humans.

Phylodynamics of sunflower chlorotic mottle virus, an emerging pathosystem

Distribution and epidemiological patterns of sunflower chlorotic mottle virus (SCMoV) in sunflower (Helianthus annuus L.) growing areas in Argentina were studied from 2006 to 2017. The virus was detected exclusively in the Pampas region (Entre Ríos, Santa Fe, Córdoba, La Pampa and Buenos Aires provinces). Phylodynamic analyses performed using the coat protein gene of SCMoV isolates from sunflower and weeds dated the most recent common ancestor (MRCA) back to 1887 (HPD95% = 1572-1971), which coincides with the dates of sunflower introduction in Argentina. The MRCA was located in the south of Buenos Aires province and was associated with sunflower host (posterior probability for the ancestral host, ppah = 0.98). The Bayesian phylodynamic analyses revealed the dispersal patterns of SCMoV, suggesting a link between natural host diversity, crop displacement by human activities and virus spread.

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.