Correlating viral phenotypes with phylogeny: Accounting for phylogenetic uncertainty

Viral diversity and co-evolutionary dynamics across the ant phylogeny

Knowledge of viral biodiversity within insects, particularly within ants, is extremely limited with only a few environmental viruses from invasive ant species identified to date. This study documents and explores the viral communities in ants. We comprehensively profile the metagenomes of a phylogenetically broad group of 35 ant species with varied ecological traits and report the discovery of 3710 novel and unique ant-associated viral genomes. These previously unknown viruses discovered within this study constitute over 95% of all currently described ant viruses, significantly increasing our knowledge of the ant virosphere. The identified RNA and DNA viruses fill gaps in insect-associated viral phylogenies and uncover evolutionary histories characterized by both frequent host switching and co-divergence. Many ants also host diverse bacterial communities, and we discovered that approximately one-third of these new ant-associated viruses are bacteriophages. Two ecological categories, bacterial abundance in the host and habitat degradation are both correlated with ant viral diversity and help to structure viral communities within ants. These data demonstrate that the ant virosphere is remarkably diverse phylogenetically and genomically and provide a substantial foundation for studies in virus ecology and evolution within eukaryotes. We highlight the importance of studying insect-associated viruses in natural ecosystems in order to more thoroughly and effectively understand host-microbe evolutionary dynamics.

Bat Rabies in the Americas: Is Myotis the Main Ancestral Spreader?

The rabies virus (RABV) is the exclusive lyssavirus affecting both wild and domestic mammalian hosts in the Americas, including humans. Additionally, the Americas stand out as the sole region where bat rabies occurs. While carnivore rabies is being increasingly managed across the region, bats are emerging as significant reservoirs of RABV infection for humans and domestic animals. Knowledge of the bat species maintaining rabies and comprehending cross-species transmission (CST) and host shift processes are pivotal for directing surveillance as well as ecological research involving wildlife reservoir hosts. Prior research indicates that bat RABV CST is influenced by host genetic similarity and geographic overlap, reflecting host adaptation. In this study, we compiled and analyzed a comprehensive nucleoprotein gene dataset representing bat-borne RABV diversity in Argentina and the broader Americas using Bayesian phylogenetics. We examined the association between host genus and geography, finding both factors shaping the global phylogenetic structure. Utilizing a phylogeographic approach, we inferred CST and identified key bat hosts driving transmission. Consistent with CST determinants, we observed monophyletic/paraphyletic clustering of most bat genera in the RABV phylogeny, with stronger CST evidence between host genera of the same family. We further discuss Myotis as a potential ancestral spreader of much of RABV diversity.

Evolution, recombination and geographic spreading of global Coxsackievirus A6

BACKGROUND

The increasing incidence of hand, foot, and mouth disease (HFMD) associated with Coxsackievirus A6 (CVA6) has become a very significant public health problem. The aim of this study is to investigate the recombination, geographic transmission, and evolutionary characteristics of the global CVA6.

METHODS

From 2019 to 2022, 73 full-length CVA6 sequences were obtained from HFMD patients in China and analyzed in combination with 1032 published whole genome sequences. Based on this dataset, the phylogenetic features, recombinant diversity, Bayesian phylodynamic characteristics, and key amino acid variations in CVA6 were analyzed.

RESULTS

The four genotypes of CVA6, A, D, E, and F, are divided into 24 recombinant forms (RFs, RF-A - RF-X) based on differences in the P3 coding region. The eastern China region plays a key role in the dissemination of CVA6 in China. VP1-137 and VP1-138 are located in the DE loop on the surface of the CVA6 VP1 protein, with the former being a highly variable site and the latter having more non-synonymous substitutions.

CONCLUSIONS

Based on whole genome sequences, this study contributes to the CVA6 monitoring, early warning, and the pathogenic mechanism by studying recombination diversity, geographical transmission characteristics, and the variation of important amino acid sites.

The distinct spatiotemporal evolutionary landscape of HBV and HDV largely determines the unique epidemic features of HDV globally

Chronic infection of hepatitis B virus (HBV) and hepatitis D virus (HDV) causes the most severe form of viral hepatitis. Due to the dependence on HBV, HDV was deemed to co-evolve and co-migrate with HBV. However, we previously found that the naturally occurred HDV/HBV combinations do not always reflect the most efficient virological adaptation (Wang et al., 2021). Moreover, regions with heavy HBV burden do not always correlate with high HDV prevalence (e.g., East Asia), and vice versa (e.g., Central Asia). Herein, we systematically elucidated the spatiotemporal evolutionary landscape of HDV to understand the unique epidemic features of HDV. We found that the MRCA of HDV was from South America around the late 13th century, was globally dispersed mainly via Central Asia, and evolved into eight genotypes from the 19th to 20th century. In contrast, the MRCA of HBV was from Europe ∼23.7 thousand years ago (Kya), globally dispersed mainly via Africa and East Asia, and evolved into eight genotypes ∼1100 years ago. When HDV stepped in, all present-day HBV genotypes had already formed and its global genotypic distribution had stayed stable geographically. Nevertheless, regionalized HDV adapted to local HBV genotypes and human lineages, contributing to the global geographical separation of HDV genotypes. Additionally, a sharp increase in HDV infections was observed after the 20th century. In conclusion, HDV exhibited a distinct spatiotemporal distribution path compared with HBV. This unique evolutionary relationship largely fostered the unique epidemic features we observe nowadays. Moreover, HDV infections may continue to ramp up globally, thus more efforts are urgently needed to combat this disease.

Evolution and host-specific adaptation of Pseudomonas aeruginosa

The major human bacterial pathogen Pseudomonas aeruginosa causes multidrug-resistant infections in people with underlying immunodeficiencies or structural lung diseases such as cystic fibrosis (CF). We show that a few environmental isolates, driven by horizontal gene acquisition, have become dominant epidemic clones that have sequentially emerged and spread through global transmission networks over the past 200 years. These clones demonstrate varying intrinsic propensities for infecting CF or non-CF individuals (linked to specific transcriptional changes enabling survival within macrophages); have undergone multiple rounds of convergent, host-specific adaptation; and have eventually lost their ability to transmit between different patient groups. Our findings thus explain the pathogenic evolution of P. aeruginosa and highlight the importance of global surveillance and cross-infection prevention in averting the emergence of future epidemic clones.

The origin and spread of HIV-1 CRF59_01B epidemic in China: A molecular network and phylogeographic analysis

Human immunodeficiency virus type 1 CRF59_01B, identified in China in 2013, has been detected nationwide, exhibiting notably high prevalence in Guangzhou and its vicinity. This study aimed to unravel its origin and migration. A data set was established, incorporating all available CRF59_01B pol gene sequences and their metadata from Guangzhou and the public database. Bayesian phylogeographic analysis demonstrated that CRF59_01B originated in Shenzhen, the neighboring city of Guangzhou, around 1998 with posterior probability of 0.937. Molecular network analysis detected 1131 transmission links and showed a remarkably high clustering rate (78.9%). Substantial inter-city transmissions (26.5%, 300/1131) were observed between Shenzhen and Guangzhou while inter-region transmissions linked Guangzhou with South (46) and Southwest (64) China. The centre of Guangzhou was the hub of CRF59_01B transmission, including the inflow from Shenzhen (3.57 events/year) and outflow to the outskirts of Guangzhou (>2 events/year). The large-scale analysis revealed significant migration from Shenzhen to Guangzhou (5.08 events/year) and North China (0.59 events/year), and spread from Guangzhou to Central (0.47 events/year), East (0.42 events/year), South (0.76 events/year), Southwest China (0.76 events/year) and Shenzhen (1.89 events/year). Shenzhen and Guangzhou served as the origin and the hub of CRF59_01B circulation, emphasizing inter-city cooperation and data sharing to confine its nationwide diffusion.

Association of disease severity and genetic variation during primary Respiratory Syncytial Virus infections

BACKGROUND

Respiratory Syncytial Virus (RSV) disease in young children ranges from mild cold symptoms to severe symptoms that require hospitalization and sometimes result in death. Studies have shown a statistical association between RSV subtype or phylogenic lineage and RSV disease severity, although these results have been inconsistent. Associations between variation within RSV gene coding regions or residues and RSV disease severity has been largely unexplored.

METHODS

Nasal swabs from children (< 8 months-old) infected with RSV in Rochester, NY between 1977-1998 clinically presenting with either mild or severe disease during their first cold-season were used. Whole-genome RSV sequences were obtained using overlapping PCR and next-generation sequencing. Both whole-genome phylogenetic and non-phylogenetic statistical approaches were performed to associate RSV genotype with disease severity.

RESULTS

The RSVB subtype was statistically associated with disease severity. A significant association between phylogenetic clustering of mild/severe traits and disease severity was also found. GA1 clade sequences were associated with severe disease while GB1 was significantly associated with mild disease. Both G and M2-2 gene variation was significantly associated with disease severity. We identified 16 residues in the G gene and 3 in the M2-2 RSV gene associated with disease severity.

CONCLUSION

These results suggest that phylogenetic lineage and the genetic variability in G or M2-2 genes of RSV may contribute to disease severity in young children undergoing their first infection.

Molecular evolution and phylogeographic analysis of wheat dwarf virus

Wheat dwarf virus (WDV) has caused considerable economic loss in the global production of grain crops. Knowledge of the evolutionary biology and population history of the pathogen remain poorly understood. We performed molecular evolution and worldwide phylodynamic analyses of the virus based on the genes in the protein-coding region of the entire viral genome. Our results showed that host-driven and geography-driven adaptation are major factors that affects the evolution of WDV. Bayesian phylogenetic analysis estimates that the average WDV substitution rate was 4.240 × 10-4 substitutions/site/year (95% credibility interval, 2.828 × 10-4-5.723 × 10-4), and the evolutionary rates of genes encoding proteins with virion-sense transcripts and genes encoding proteins with complementary-sense transcripts were different. The positively selected sites were detected in only two genes encoding proteins with complementary-sense, and WDV-barley are subject to stronger purifying selection than WDV-wheat. The time since the most recent common WDV ancestor was 1746 (95% credibility interval, 1517-1893) CE. Further analyses identified that the WDV-barley population and WDV-wheat population experienced dramatic expansion-decline episodes, and the expansion time of the WDV-barley population was earlier than that of the WDV-wheat population. Our phylogeographic analysis showed that the WDV population originating in Iran was subsequently introduced to Europe, and then spread from Eastern Europe to China.

Investigating the rise of Omicron variant through genomic surveillance of SARS-CoV-2 infections in a highly vaccinated university population

Novel variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continue to emerge as the coronavirus disease 2019 (COVID-19) pandemic extends into its fourth year. Understanding SARS-CoV-2 circulation in university populations is vital for effective interventions in higher education settings and will inform public health policy during pandemics. In this study, we generated 793 whole-genome sequences collected over an entire academic year in a university population in Indiana, USA. We clearly captured the rapidity with which Delta variant was wholly replaced by Omicron variant across the West Lafayette campus over the length of two academic semesters in a community with high vaccination rates. This mirrored the emergence of Omicron throughout the state of Indiana and the USA. Further, phylogenetic analyses demonstrated that there was a more diverse set of potential geographic origins for Omicron viruses introduction into campus when compared to Delta. Lastly, statistics indicated that there was a more significant role for international and out-of-state migration in the establishment of Omicron variants at Purdue. This surveillance workflow, coupled with viral genomic sequencing and phylogeographic analyses, provided critical insights into SARS-CoV-2 transmission dynamics and variant arrival.

Synergetic association between coxsackievirus A16 genotype evolution and recombinant form shifts

Coxsackievirus A16 (CVA16) is a major pathogen that causes hand, foot, and mouth disease (HFMD). The recombination form (RF) shifts and global transmission dynamics of CVA16 remain unknown. In this retrospective study, global sequences of CVA16 were retrieved from the GenBank database and analyzed using comprehensive phylogenetic inference, RF surveys, and population structure. A total of 1,663 sequences were collected, forming a 442-sequences dataset for VP1 coding region analysis and a 345-sequences dataset for RF identification. Based on the VP1 coding region used for serotyping, three genotypes (A, B, and D), two subgenotypes of genotype B (B1 and B2), and three clusters of subgenotype B1 (B1a, B1b, and B1c) were identified. Cluster B1b has dominated the global epidemics, B2 disappeared in 2000, and D is an emerging genotype dating back to August 2002. Globally, four oscillation phases of CVA16 evolution, with a peak in 2013, and three migration pathways were identified. Europe, China, and Japan have served as the seeds for the global transmission of CVA16. Based on the 3D coding region of the RFs, five clusters of RFs (RF-A to -E) were identified. The shift in RFs from RF-B and RF-C to RF-D was accompanied by a change in genotype from B2 to B1a and B1c and then to B1b. In conclusion, the evolution and population dynamics of CVA16, especially the coevolution of 3D and VP1 genes, revealed that genotype evolution and RF replacement were synergistic rather than stochastic.

Reconstructing the phylodynamic history and geographic spread of the CRF01_AE-predominant HIV-1 epidemic in the Philippines from PR/RT sequences sampled from 2008 to 2018

The Philippines has had a rapidly growing human immunodeficiency virus (HIV) epidemic with a shift in the prevalent subtype from B to CRF01_AE. However, the phylodynamic history of CRF01_AE in the Philippines has yet to be reconstructed. We conducted a descriptive retrospective study reconstructing the history of HIV-1 CRF01_AE transmissions in the Philippines through molecular epidemiology. Partial polymerase sequences (n = 1144) collected between 2008 and 2018 from three island groups were collated from the Research Institute for Tropical Medicine drug resistance genotyping database. Estimation of the time to the most recent common ancestor (tMRCA), effective reproductive number (Re), effective viral population size (Ne), relative migration rates, and geographic spread of CRF01_AE was performed with BEAST. Re and Ne were compared between CRF01_AE and B. Most CRF01_AE sequences formed a single clade with a tMRCA of June 1996 [95 per cent highest posterior density (HPD): December 1991, October 1999]. An increasing CRF01_AE Ne was observed from the tMRCA to 2013. The CRF01_AE Re reached peaks of 2.46 [95 per cent HPD: 1.76, 3.27] in 2007 and 2.52 [95 per cent HPD: 1.83, 3.34] in 2015. A decrease of CRF01_AE Re occurred in the intervening years of 2007 to 2011, reaching as low as 1.43 [95 per cent HPD: 1.06, 1.90] in 2011, followed by a rebound. The CRF01_AE epidemic most likely started in Luzon and then spread to the other island groups of the country. Both CRF01_AE and Subtype B exhibited similar patterns of Re fluctuation over time. These results characterize the subtype-specific phylodynamic history of the largest CRF01_AE cluster in the Philippines, which contextualizes and may inform past, present, and future public health measures toward controlling the HIV epidemic in the Philippines.

The episodic resurgence of highly pathogenic avian influenza H5 virus

Highly pathogenic avian influenza (HPAI) H5N1 activity has intensified globally since 2021, increasingly causing mass mortality in wild birds and poultry and incidental infections in mammals1-3. However, the ecological and virological properties that underscore future mitigation strategies still remain unclear. Using epidemiological, spatial and genomic approaches, we demonstrate changes in the origins of resurgent HPAI H5 and reveal significant shifts in virus ecology and evolution. Outbreak data show key resurgent events in 2016-2017 and 2020-2021, contributing to the emergence and panzootic spread of H5N1 in 2021-2022. Genomic analysis reveals that the 2016-2017 epizootics originated in Asia, where HPAI H5 reservoirs are endemic. In 2020-2021, 2.3.4.4b H5N8 viruses emerged in African poultry, featuring mutations altering HA structure and receptor binding. In 2021-2022, a new H5N1 virus evolved through reassortment in wild birds in Europe, undergoing further reassortment with low-pathogenic avian influenza in wild and domestic birds during global dissemination. These results highlight a shift in the HPAI H5 epicentre beyond Asia and indicate that increasing persistence of HPAI H5 in wild birds is facilitating geographic and host range expansion, accelerating dispersion velocity and increasing reassortment potential. As earlier outbreaks of H5N1 and H5N8 were caused by more stable genomic constellations, these recent changes reflect adaptation across the domestic-bird-wild-bird interface. Elimination strategies in domestic birds therefore remain a high priority to limit future epizootics.

Tobamoviruses Show Broad Host Ranges and Little Genetic Diversity Among Four Habitat Types of a Heterogeneous Ecosystem

Host ranges of plant viruses are poorly known, as studies have focused on pathogenic viruses in crops and adjacent wild plants. High-throughput sequencing (HTS) avoids the bias toward plant-virus interactions that result in disease. Here we study the host ranges of tobamoviruses, important pathogens of crops, using HTS analyses of an extensive sample of plant communities in four habitats of a heterogeneous ecosystem. Sequences of 17 virus operational taxonomic units (OTUs) matched references in the Tobamovirus genus, eight had narrow host ranges, and five had wide host ranges. Regardless of host range, the OTU hosts belonged to taxonomically distant families, suggesting no phylogenetic constraints in host use associated with virus adaptation, and that tobamoviruses may be host generalists. The OTUs identified as tobacco mild green mosaic virus (TMGMV), tobacco mosaic virus (TMV), pepper mild mottle virus, and Youcai mosaic virus had the largest realized host ranges that occurred across habitats and exhibited host use unrelated to the degree of human intervention. This result is at odds with assumptions that contact-transmitted viruses would be more abundant in crops than in wild plant communities and could be explained by effective seed-, contact-, or pollinator-mediated transmission or by survival in the soil. TMGMV and TMV had low genetic diversity that was not structured according to habitat or host plant taxonomy, which indicated that phenotypic plasticity allows virus genotypes to infect new hosts with no need for adaptive evolution. Our results underscore the relevance of ecological factors in host range evolution, in addition to the more often studied genetic factors. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Ecological Strategies for Resource Use by Three Bromoviruses in Anthropic and Wild Plant Communities

Ecological strategies for resource utilisation are important features of pathogens, yet have been overshadowed by stronger interest in genetic mechanisms underlying disease emergence. The purpose of this study is to ask whether host range and transmission traits translate into ecological strategies for host-species utilisation in a heterogeneous ecosystem, and whether host utilisation corresponds to genetic differentiation among three bromoviruses. We combine high-throughput sequencing and population genomics with analyses of species co-occurrence to unravel the ecological strategies of the viruses across four habitat types. The results show that the bromoviruses that were more closely related genetically did not share similar ecological strategies, but that the more distantly related pair did. Shared strategies included a broad host range and more frequent co-occurrences, which both were habitat-dependent. Each habitat thus presents as a barrier to gene flow, and each virus has an ecological strategy to navigate limitations to colonising non-natal habitats. Variation in ecological strategies could therefore hold the key to unlocking events that lead to emergence.

Complete genomic characterization of bovine papillomavirus type 1 and 2 strains infers ongoing cross-species transmission between cattle and horses

Infection with bovine papillomavirus (BPV) types 1 and 2 results in the most common skin tumor of horses, termed equine sarcoid. The persistent and recurrent nature of this tumor stands in contrast to the regressive nature of BPV-1/- 2 induced cutaneous papillomas in cattle. The circulation of horse-specific BPV-1/- 2 variants within equine populations has been suggested as a possible explanation for the difference in clinical presentation of BPV-1/- 2 infection between horses and cattle. In order to investigate this hypothesis, we identified 98 complete BPV-1/- 2 genomes using a Nanopore sequencing approach. Separate BPV-1/- 2 alignments were used to infer Bayesian phylogenetic trees. Phylogeny-trait association concerning host species was investigated using Bayesian Tip-association Significance software (BaTS) Overall, 179 unique BPV-1 and 128 BPV-2 substitutions were found. The E2 coding region in the viral genome exhibited an exceptionally high rate of non-synonymous mutations (81 %, n = 13/16). Interestingly, extensive deletions in the L1/L2 region (up to 1.5 kb) were found exclusively in horse-derived samples. Nevertheless, the most frequently detected single nucleotide polymorphisms were shared between equine and bovine hosts, which is in agreement with BaTS results indicating no phylogeny-host correlation. We found indications that horse-specific mutations might exist in subpopulations of equine derived BPV-1/- 2, but these did not result in horse-adapted genetic variants. Based on these observations, cross-species transmission from cattle to horses seems to be an ongoing process, rather than an ancient occurrence that has been followed by the circulation of horse-adapted BPV variants in the horse population..

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

Mining-impacted rice paddies select for Archaeal methylators and reveal a putative (Archaeal) regulator of mercury methylation

Methylmercury (MeHg) is a microbially produced neurotoxin derived from inorganic mercury (Hg), which accumulation in rice represents a major health concern to humans. However, the microbial control of MeHg dynamics in the environment remains elusive. Here, leveraging three rice paddy fields with distinct concentrations of Hg (Total Hg (THg): 0.21-513 mg kg-1 dry wt. soil; MeHg: 1.21-6.82 ng g-1 dry wt. soil), we resorted to metagenomics to determine the microbial determinants involved in MeHg production under contrasted contamination settings. We show that Hg methylating Archaea, along with methane-cycling genes, were enriched in severely contaminated paddy soils. Metagenome-resolved Genomes of novel putative Hg methylators belonging to Nitrospinota (UBA7883), with poorly resolved taxonomy despite high completeness, showed evidence of facultative anaerobic metabolism and adaptations to fluctuating redox potential. Furthermore, we found evidence of environmental filtering effects that influenced the phylogenies of not only hgcA genes under different THg concentrations, but also of two housekeeping genes, rpoB and glnA, highlighting the need for further experimental validation of whether THg drives the evolution of hgcAB. Finally, assessment of the genomic environment surrounding hgcAB suggests that this gene pair may be regulated by an archaeal toxin-antitoxin (TA) system, instead of the more frequently found arsR-like genes in bacterial methylators. This suggests the presence of distinct hgcAB regulation systems in bacteria and archaea. Our results support the emerging role of Archaea in MeHg cycling under mining-impacted environments and shed light on the differential control of the expression of genes involved in MeHg formation between Archaea and Bacteria.

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.

Molecular Characteristics of Enterovirus B83 Strain Isolated from a Patient with Acute Viral Myocarditis and Global Transmission Dynamics

This study determined the global genetic diversity and transmission dynamics of enterovirus B83 (EV-B83) and proposed future disease surveillance directions. Blood samples were collected from a patient with viral myocarditis, and viral isolation was performed. The complete genome sequence of the viral isolate was obtained using Sanger sequencing. A dataset of 15 sequences (from three continents) that had sufficient time signals for Bayesian phylogenetic analysis was set up, and the genetic diversity and transmission dynamics of global EV-B83 were analyzed using bioinformatics methods, including evolutionary dynamics, recombination event analysis, and phylogeographic analysis. Here, we report the complete genome sequence of an EV-B83 strain (S17/YN/CHN/2004) isolated from a patient with acute viral myocarditis in Yunnan Province, China. All 15 EV-B83 strains clustered together in a phylogenetic tree, confirming the classification of these isolates as a single EV type, and the predicted time for the most recent common ancestor appeared in 1998. Recombinant signals were detected in the 5'-untranslated region and 2A-3D coding regions of the S17 genome. The phylogeographic analysis revealed multiple intercontinental transmission routes of EV-B83. This study indicates that EV-B83 is globally distributed. Our findings add to the publicly available EV-B83 genomic sequence data and deepen our understanding of EV-B83 epidemiology.

Emergence, phylogeography, and adaptive evolution of mpox virus

Mpox (Monkeypox) is a zoonotic disease caused by mpox virus (MPXV). A multi-country MPXV outbreak in non-endemic demographics was identified in May 2022. A systematic evaluation of MPXV evolutionary trajectory and genetic diversity could be a timely addition to the MPXV diagnostics and prophylaxis. Herein, we integrated a systematic evolution analysis including phylogenomic and phylogeographic, followed by an in-depth analysis of the adaptive evolution and amino acid variations in type I interferon binding protein (IFNα/βBP). Mutations in IFNα/βBP protein may impair its binding capacity, affecting the MPXV immune evasion strategy. Based on the equilibrated data, we found an evolutionary rate of 7.75 × 10 - 5 substitutions/site/year, and an earlier original time (2021.25) of the clade IIb. We further discovered significant genetic variations in MPXV genomes from different regions and obtained six plausible spread trajectories from its intricate viral flow network, implying that North America might have acted as a bridge for the spread of MPXV from Africa to other continents. We identified two amino acids under positive selection in the Rifampicin resistance protein and extracellular enveloped virus (EEV) type-I membrane glycoprotein, indicating a role in adaptive evolution. Our research sheds light on the emergence, dispersal, and adaptive evolution of MPXV, providing theoretical support for mitigating and containing its expansion.

Influence of Sex on Respiratory Syncytial Virus Genotype Infection Frequency and Nasopharyngeal Microbiome

Respiratory syncytial virus (RSV) has a significant health burden in children, older adults, and the immunocompromised. However, limited effort has been made to identify emergence of new RSV genotypes' frequency of infection and how the combination of nasopharyngeal microbiome and viral genotypes impact RSV disease outcomes. In an observational cohort designed to capture the first infant RSV infection, we employed multi-omics approaches to sequence 349 RSV complete genomes and matched nasopharyngeal microbiomes, during which the 2012/2013 season was dominated by RSV-A, whereas 2013 and 2014 was dominated by RSV-B. We found non-G-72nt-duplicated RSV-A strains were more frequent in male infants (P = 0.02), whereas G-72nt-duplicated genotypes (which is ON1 lineage) were seen equally in both males and females. DESeq2 testing of the nasal microbiome showed Haemophilus was significantly more abundant in infants with RSV-A infection compared to infants with RSV-B infection (adjusted P = 0.002). In addition, the broad microbial clustering of the abundant genera was significantly associated with infant sex (P = 0.03). Overall, we show sex differences in infection by RSV genotype and host nasopharyngeal microbiome, suggesting an interaction between host genetics, virus genotype, and associated nasopharyngeal microbiome. IMPORTANCE Respiratory syncytial virus (RSV) is one of the leading causes of lower respiratory tract infections in young children and is responsible for high hospitalization rates and morbidity in infants and the elderly. To understand how the emergence of RSV viral genotypes and viral-respiratory microbiome interactions contribute to infection frequency and severity, we utilized an observational cohort designed to capture the first infant RSV infection we employed multi-omics approaches to sequence 349 RSV complete genomes and matched nasopharyngeal microbiomes. We found non-G-72nt-duplicated RSV-A genotypes were more frequent in male infants, whereas G-72nt-duplicated RSV-A strains (ON1 lineage) were seen equally in both males and females. Microbiome analysis show Haemophilus was significantly more abundant in infants with RSV-A compared to infants with RSV-B infection and the microbial clustering of the abundant genera was associated with infant sex. Overall, we show sex differences in RSV genotype-nasopharyngeal microbiome, suggesting an interaction host genetics-virus-microbiome interaction.

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.

Population genetics and phylogeography of alfalfa mosaic virus in China and a comparison with other regional epidemics based on the cp gene

Alfalfa mosaic virus (AMV) is the most pervasive epidemic virus affecting alfalfa production. However, detailed investigations on the molecular population genetics and evolutionary dynamics of AMV are scarce. This study aimed to report on a large-scale long-term survey of genetic variability in AMV populations from China and perform a comparative analysis of AMV population genetics in the three most thoroughly studied countries to date: China, Iran, and Spain. The study was based on the analysis of the coat protein gene (cp) using two analytical approaches: an analysis of molecular variance (AMOVA) and a Bayesian Markov Chain Monte Carlo approach that investigates the association between geographical origin and phylogeny. Both analytical approaches found significant genetic differentiation within localities, but not among localities nor among provinces. This observation might result from inappropriate agronomical practices involving extensive exchange of plant materials followed by rapid viral diversification within localities. In the Chinese population, both methods found that genetic diversification in AMV was strongly associated with different bioclimatic zones. Rates of molecular evolution were similar in the three countries. The estimated epidemic exponential population size and growth rate suggest that the epidemics grew faster and with higher incidence in Iran, followed by Spain and China. Estimates of the time to the most recent common ancestors suggest that AMV was first seen in Spain by the beginning of the twentieth century and later on in eastern and central Eurasia. After ruling out the existence of recombination breakpoints within the cp gene, a codon-based selection analysis per population was performed and identified many codons under significant negative selection and a few under significant positive selection; the latter varied among countries, suggesting regional differences in selective pressures.

Metavirome of 31 tick species provides a compendium of 1,801 RNA virus genomes

The increasing prevalence and expanding distribution of tick-borne viruses globally have raised health concerns, but the full repertoire of the tick virome has not been assessed. We sequenced the meta-transcriptomes of 31 different tick species in the Ixodidae and Argasidae families from across mainland China, and identified 724 RNA viruses with distinctive virome compositions among genera. A total of 1,801 assembled and complete or nearly complete viral genomes revealed an extensive diversity of genome architectures of tick-associated viruses, highlighting ticks as a reservoir of RNA viruses. We examined the phylogenies of different virus families to investigate virome evolution and found that the most diverse tick-associated viruses are positive-strand RNA virus families that demonstrate more ancient divergence than other arboviruses. Tick-specific viruses are often associated with only a few tick species, whereas virus clades that can infect vertebrates are found in a wider range of tick species. We hypothesize that tick viruses can exhibit both 'specialist' and 'generalist' evolutionary trends. We hope that our virome dataset will enable much-needed research on vertebrate-pathogenic tick-associated viruses.

Phylogeographic analysis of Tula hantavirus highlights a single introduction to central Europe

Orthohantaviruses are zoonotic pathogens of humans, unique among the bunyaviruses in not being transmitted by an arthropod vector. Tula orthohantavirus (TULV) is an old-world hantavirus, of yet unclear human pathogenicity, with few reported cases of clinically relevant human infection. So far, phylogeographic studies exploring the global pathways of hantaviral migration are scarce and generally do not focus on a specific hantavirus species. The aim of the present study was to reconstruct the dispersal history of TULV lineages across Eurasia based on S segment sequences sampled from different geographic areas. Maximum-likelihood and Bayesian inference methods were used to perform the phylogenetic analysis and phylogeographic reconstructions. Sampling time and trapping localities were obtained for a total of 735 TULV S segment sequences available in public databases at the time of the study. The estimated substitution rate of the analyzed partial S segment alignment was 2.26 × 10-3 substitutions/site/year (95 per cent highest posterior density interval: 1.79 × 10-3 to 2.75 × 10-3). Continuous phylogeography of TULV S segment sequences placed the potential root and origin of TULV spread in the Black Sea region. In our study, we detect a single-lineage introduction of TULV to Europe, followed by local viral circulation further on.

Evolutionary history and global spatiotemporal pattern of alfalfa mosaic virus

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

The Role of Molossidae and Vespertilionidae in Shaping the Diversity of Alphacoronaviruses in the Americas

Bats are reservoirs of diverse coronaviruses (CoVs), including progenitors of severe acute respiratory syndrome CoV (SARS-CoV) and SARS-CoV-2. In the Americas, there is a contrast between alphacoronaviruses (alphaCoVs) and betaCoVs: while cospeciation prevails in the latter, alphaCoV evolution is dominated by deep and recent host switches. AlphaCoV lineages are maintained by two different bat family groups, Phyllostomidae and Vespertilionidae plus Molossidae. In this study, we used a Bayesian framework to analyze the process of diversification of the lineages maintained by Molossidae and Vespertilionidae, adding novel CoV sequences from Argentina. We provide evidence that the observed CoV diversity in these two bat families is shaped by their geographic distribution and that CoVs exhibit clustering at the level of bat genera. We discuss the causes of the cocirculation of two independent clades in Molossus and Tadarida as well as the role of Myotis as the ancestral host and a major evolutionary reservoir of alphaCoVs across the continent. Although more CoV sampling efforts are needed, these findings contribute to a better knowledge of the diversity of alphaCoVs and the links between bat host species. IMPORTANCE Bats harbor the largest diversity of coronaviruses among mammals. In the Americas, seven alphacoronavirus lineages circulate among bats. Three of these lineages are shared by members of two bat families: Vespertilionidae and Molossidae. Uncovering the relationships between these coronaviruses can help us to understand patterns of cross-species transmission and, ultimately, which hosts are more likely to be involved in spillover events. We found that two different lineages cocirculate among the bat genera Molossus and Tadarida, which share roosts and have common viral variants. The bat genus Myotis functions as a reservoir of coronavirus diversity and, as such, is a key host. Although there were some spillovers recorded, there is a strong host association, showing that once a successful host jump takes place, it is transmitted onward to members of the same bat genus.

Genetic Features of Mycobacterium avium subsp. paratuberculosis Strains Circulating in the West of France Deciphered by Whole-Genome Sequencing

Paratuberculosis is a chronic infection of the intestine, mainly the ileum, caused by Mycobacterium avium subsp. paratuberculosis in cattle and other ruminants. This enzootic disease is present worldwide and has a negative impact on the dairy cattle industry. For this subspecies, the current genotyping tools do not provide the needed resolution to investigate the genetic diversity of closely related strains. These limitations can be overcome by the application of whole-genome sequencing (WGS), particularly for clonal populations such as M. avium subsp. paratuberculosis. The purpose of the present study was to undertake a WGS analysis with a panel of 200 animal field M. avium subsp. paratuberculosis strains selected based on a previous large-scale longitudinal study of Prim'Holstein and Normande dairy breeds naturally infected with M. avium subsp. paratuberculosis in the West of France. The pangenome analysis revealed that M. avium subsp. paratuberculosis has a closed pangenome. The phylogeny, based on alignment of 2,786 nonhomoplasic single nucleotide polymorphisms (SNPs), showed that the strain population is structured into three clades independently of the cattle breed or geographic distribution. The increased resolution of phylogeny obtained by WGS confirmed the homoplasic nature of the markers variable-number tandem repeat (VNTR) and short sequence repeat (SSR) used for M. avium subsp. paratuberculosis genotyping. These phylogenetic data also revealed independent introductions of the different genotypes in two main waves since at least 2003. WGS applied to this sampling demonstrated the presence of mixed infections in herds and at the individual animal level. Collectively, the phylogeny results inferred with French isolates compared to M. avium subsp. paratuberculosis isolates from around the world suggest introductions of M. avium subsp. paratuberculosis genotypes through the animal trade. Relationships between genetic traits and epidemiological data can now be investigated to better understand transmission dynamics of the disease. IMPORTANCE Mycobacterium avium subsp. paratuberculosis causes Johne's disease in ruminants, which is present worldwide and has significant negative impacts on the dairy cattle industry and animal welfare. Prevention and control of M. avium subsp. paratuberculosis infection are hampered by knowledge gaps in strain virulence, genotype distribution, and transmission dynamics. This work has revealed new insights into M. avium subsp. paratuberculosis strains currently circulating in western France and how they are related to strains circulating globally. We applied whole-genome sequencing (WGS) to obtain comprehensive information on genome evolution and discrimination of closely related strains. This approach revealed the history of M. avium subsp. paratuberculosis infection in France, refined the pangenomic characteristics of M. avium subsp. paratuberculosis, and demonstrated the existence of mixed infection in animals. Finally, this study identified predominant genotypes, which allow a better understanding of disease transmission dynamics. This information will facilitate tracking of this pathogen on farms and across agricultural regions, thus informing transmission pathways and disease control points.

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.

The diversity and evolutionary relationships of ticks and tick-borne bacteria collected in China

Background

Ticks (order Ixodida) are ectoparasites, vectors and reservoirs of many infectious agents affecting humans and domestic animals. However, the lack of information on tick genomic diversity leaves significant gaps in the understanding of the evolution of ticks and associated bacteria.

Results

We collected > 20,000 contemporary and historical (up to 60 years of preservation) tick samples representing a wide range of tick biodiversity across diverse geographic regions in China. Metagenomic sequencing was performed on individual ticks to obtain the complete or near-complete mitochondrial (mt) genome sequences from 46 tick species, among which mitochondrial genomes of 23 species were recovered for the first time. These new mt genomes data greatly expanded the diversity of many tick groups and revealed five cryptic species. Utilizing the same metagenomic sequence data we identified divergent and abundant bacteria in Haemaphysalis, Ixodes, Dermacentor and Carios ticks, including nine species of pathogenetic bacteria and potentially new species within the genus Borrelia. We also used these data to explore the evolutionary relationship between ticks and their associated bacteria, revealing a pattern of long-term co-divergence relationship between ticks and Rickettsia and Coxiella bacteria.

Conclusions

In sum, our study provides important new information on the genetic diversity of ticks based on an analysis of mitochondrial DNA as well as on the prevalence of tick-borne pathogens in China. It also sheds new light on the long-term evolutionary and ecological relationships between ticks and their associated bacteria.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05485-3.

Lineage BA.2 dominated the Omicron SARS-CoV-2 epidemic wave in the Philippines

The Omicron severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant led to a dramatic global epidemic wave following detection in South Africa in November 2021. The BA.1 Omicron lineage was dominant and responsible for most SARS-CoV-2 outbreaks in countries around the world during December 2021-January 2022, while other Omicron lineages, including BA.2, accounted for the minority of global isolates. Here, we describe the Omicron wave in the Philippines by analysing genomic data. Our results identify the presence of both BA.1 and BA.2 lineages in the Philippines in December 2021, before cases surged in January 2022. We infer that only the BA.2 lineage underwent sustained transmission in the country, with an estimated emergence around 18 November 2021 (95 per cent highest posterior density: 6-28 November), while despite multiple introductions, BA.1 transmission remained limited. These results suggest that the Philippines was one of the earliest areas affected by BA.2 and reiterate the importance of whole genome sequencing for monitoring outbreaks.

Association of Hepatitis C Virus Genotype 2 Spread With Historic Slave Trade and Commerce Routes in Western Africa

The hepatitis C virus genotype 2 (HCV2) is endemic in Western and Central Africa. The HCV2 evolutionary origins remain uncertain due to the paucity of available genomes from African settings. In this study, we investigated the molecular epidemiology of HCV infections in rural Guinea, Western Africa, during 2004 and 2014. Broadly reactive nested RT-PCR-based screening of sera from 1,571 asymptomatic adults resulted in the detection of 25 (1.5%, 95% CI 0.9-2.3) positive samples, with a median viral load of 2.54E+05 IU/mL (IQR 6.72E+05). HCV-infected persons had a median age of 47 years and 62.5% were male and 37.5% female. The full polyprotein-encoding genes were retrieved by a combination of high throughput and Sanger sequencing from 17 samples showing sufficiently high viral loads. Phylogenetic analysis and sequence distances ≥13% averaged over the polyprotein genes compared to other HCV2 subtypes revealed 9 previously unknown HCV2 subtypes. The time to the most recent common ancestor of the Guinean HCV2 strains inferred in a Bayesian framework was 493 years (95% HPD 453-532). Most of the Guinean strains clustered poorly by location both on the level of sampling sites within Guinea and the level of countries in the phylogenetic reconstructions. Ancestral state reconstruction provided decisive support (Bayes factor &gt;100) for an origin of HCV2 in Western Africa. Phylogeographic reconstructions in a Bayesian framework pointed to a radial diffusion of HCV2 from Western African regions encompassing today’s countries like Ghana, Guinea Bissau, or Burkina Faso, to Central and Northern African regions that took place from the 16th century onwards. The spread of HCV2 coincided in time and space with the main historic slave trade and commerce routes, supported by Bayesian tip-association significance testing (p = 0.01). Our study confirms the evolutionary origins of HCV2 in Western Africa and provides a potential link between historic human movements and HCV2 dispersion.

Vaccine Mismatches, Viral Circulation, and Clinical Severity Patterns of Influenza B Victoria and Yamagata Infections in Brazil over the Decade 2010-2020: A Statistical and Phylogeny-Trait Analyses

Worldwide, infections by influenza viruses are considered a major public health challenge. In this study, influenza B vaccine mismatches and clinical aspects of Victoria and Yamagata infections in Brazil were assessed. Clinical samples were collected from patients suspected of influenza infection. In addition, sociodemographic, clinical, and epidemiological information were collected by the epidemiological surveillance teams. Influenza B lineages were determined by real-time RT-PCR and/or Sanger sequencing. In addition, putative phylogeny−trait associations were assessed by using the BaTS program after phylogenetic reconstruction by a Bayesian Markov Chain Monte Carlo method (BEAST software package). Over 2010−2020, B/Victoria and B/Yamagata-like lineages co-circulated in almost all seasonal epidemics, with B/Victoria predominance in most years. Vaccine mismatches between circulating viruses and the trivalent vaccine strains occurred in five of the eleven seasons (45.5%). No significant differences were identified in clinical presentation or disease severity caused by both strains, but subjects infected by B/Victoria-like viruses were significantly younger than their B/Yamagata-like counterparts (16.7 vs. 31.4 years, p < 0.001). This study contributes to a better understanding of the circulation patterns and clinical outcomes of B/Victoria- and B/Yamagata-like lineages in Brazil and advocate for the inclusion of a quadrivalent vaccine in the scope of the Brazilian National Immunization Program.

Occurrence, Genetic Variability of Tomato Yellow Ring Orthotospovirus Population and the Development of Reverse Transcription Loop-Mediated Isothermal Amplification Assay for Its Rapid Detection

Tomato-infecting viruses have been considered as a serious threat to tomato crops in Poland. Therefore, during 2014–2021, 234 tomato samples delivered directly by greenhouse tomato growers to Plant Disease Clinic of IPP-NRI were tested. Eight virus species: pepino mosaic virus (PepMV), tomato yellow ring orthotospovirus (TYRV), tomato spotted wilt orthotospovirus (TSWV), potato virus Y (PVY), cucumber mosaic virus (CMV), tomato black ring virus (TBRV) and tomato mosaic virus (ToMV) were detected in single or mixed infection in 89 samples. The presence of TYRV was established for the first time in Poland in 2014. Since then, its presence has been observed in single and mixed infection with TSWV and CMV. Here, we analysed the genetic variability of TYRV population based on complete nucleocapsid (N) protein gene sequence of 55 TYRV isolates. Maximum-likelihood reconstruction revealed the presence of three distinct, well-supported phylogroups. Moreover, the effect of host species on virus diversity was confirmed. Therefore, RT-LAMP assay was developed for the rapid and efficient detection of TYRV isolates that can be implemented in field and greenhouse conditions.

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.

Genomic surveillance of SARS-CoV-2 in a university community: insights into tracking variants, transmission, and spread of Gamma (P.1) variant

Background

Using a combination of data from routine surveillance, genomic sequencing, and phylogeographic analysis we tracked the spread and introduction events of SARS-CoV-2 variants focusing on a large university community.

Methods

Here, we sequenced and analyzed 677 high-quality SARS-CoV-2 genomes from positive RNA samples collected from Purdue University students, faculty, and staff who tested positive for the virus between January 2021 and May 2021, comprising an average of 32% of weekly cases across the time frame.

Results

Our analysis of circulating SARS-CoV-2 variants over time revealed periods when Variant of Concern (VOC) Alpha (B.1.1.7) and Iota (B.1.526) reached rapid dominance and documented that VOC Gamma (P.1) was increasing in frequency as campus surveillance was ending. Phylodynamic analysis of Gamma genomes from campus alongside a subsampling of &gt;20,000 previously published P.1 genomes revealed ten independent introductions of this variant into the Purdue community, predominantly from elsewhere in the United States, with introductions from within the state of Indiana and from Illinois, and possibly Washington and New York, suggesting a degree of domestic spread.

Conclusions

We conclude that a robust and sustained active and passive surveillance program coupled with genomic sequencing during a pandemic offers important insights into the dynamics of pathogen arrival and spread in a campus community and can help guide mitigation measures.

Comparative Phylodynamics Reveals the Evolutionary History of SARS-CoV-2 Emerging Variants in the Arabian Peninsula

Emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants continue to be responsible for an unprecedented worldwide public health and economic catastrophe. Accurate understanding and comparison of global and regional evolutionary epidemiology of novel SARS-CoV-2 variants are critical to guide current and future interventions. Here, we utilized a Bayesian phylodynamic pipeline to trace and compare the evolutionary dynamics, spatiotemporal origins, and spread of five variants (Alpha, Beta, Delta, Kappa, and Eta) across the Arabian Peninsula. We found variant-specific signatures of evolution and spread that are likely linked to air travel and disease control interventions in the region. Alpha, Beta, and Delta variants went through sequential periods of growth and decline, whereas we inferred inconclusive population growth patterns for the Kappa and Eta variants due to their sporadic introductions in the region. Non-pharmaceutical interventions imposed between mid-2020 and early 2021 likely played a role in reducing the epidemic progression of the Beta and the Alpha variants. In comparison, the combination of the non-pharmaceutical interventions and the rapid rollout of vaccination might have shaped Delta variant dynamics. We found that the Alpha and Beta variants were frequently introduced into the Arab peninsula between mid-2020 and early 2021 from Europe and Africa, respectively, whereas the Delta variant was frequently introduced between early 2021 and mid-2021 from East Asia. For these three variants, we also revealed significant and intense dispersal routes between the Arab region and Africa, Europe, Asia, and Oceania. In contrast, the restricted spread and stable effective population size of the Kappa and the Eta variants suggest that they no longer need to be targeted in genomic surveillance activities in the region. In contrast, the evolutionary characteristics of the Alpha, Beta, and Delta variants confirm the dominance of these variants in the recent outbreaks. Our study highlights the urgent need to establish regional molecular surveillance programs to ensure effective decision making related to the allocation of intervention activities targeted toward the most relevant variants.

Statistical Challenges in Tracking the Evolution of SARS-CoV-2

Genomic surveillance of SARS-CoV-2 has been instrumental in tracking the spread and evolution of the virus during the pandemic. The availability of SARS-CoV-2 molecular sequences isolated from infected individuals, coupled with phylodynamic methods, have provided insights into the origin of the virus, its evolutionary rate, the timing of introductions, the patterns of transmission, and the rise of novel variants that have spread through populations. Despite enormous global efforts of governments, laboratories, and researchers to collect and sequence molecular data, many challenges remain in analyzing and interpreting the data collected. Here, we describe the models and methods currently used to monitor the spread of SARS-CoV-2, discuss long-standing and new statistical challenges, and propose a method for tracking the rise of novel variants during the epidemic.

Phylogenetic analysis of migration, differentiation, and class switching in B cells

B cells undergo rapid mutation and selection for antibody binding affinity when producing antibodies capable of neutralizing pathogens. This evolutionary process can be intermixed with migration between tissues, differentiation between cellular subsets, and switching between functional isotypes. B cell receptor (BCR) sequence data has the potential to elucidate important information about these processes. However, there is currently no robust, generalizable framework for making such inferences from BCR sequence data. To address this, we develop three parsimony-based summary statistics to characterize migration, differentiation, and isotype switching along B cell phylogenetic trees. We use simulations to demonstrate the effectiveness of this approach. We then use this framework to infer patterns of cellular differentiation and isotype switching from high throughput BCR sequence datasets obtained from patients in a study of HIV infection and a study of food allergy. These methods are implemented in the R package dowser, available at https://dowser.readthedocs.io.

B cells produce high affinity antibodies through an evolutionary process of mutation and selection during adaptive immune responses. Migration between tissues, differentiation to cellular subtypes, and switching between different antibody isotypes can be important factors in shaping the role B cells play in response to infection, autoimmune disease, and allergies. B cell receptor (BCR) sequence data has the potential to elucidate important information about these processes. However, there is currently no robust, generalizable framework for making such inferences from BCR sequence data. Here, we develop three parsimony-based summary statistics to characterize migration, differentiation, and isotype switching along B cell phylogenetic trees. We confirm the effectiveness of our approach using simulations and further apply our method to data from patients with HIV and allergy. Our methods are released in the R package dowser: https://dowser.readthedocs.io.

Molecular Epidemiology and Transmission Dynamics of the HIV-1 Epidemic in Ethiopia: Epidemic Decline Coincided With Behavioral Interventions Before ART Scale-Up

Background

Ethiopia is one of the sub-Saharan countries hit hard by the HIV epidemic. Previous studies have shown that subtype C dominates the Ethiopian HIV-1 epidemic, but the evolutionary and temporal dynamics of HIV-1 in Ethiopia have not been closely scrutinized. Understanding the evolutionary and epidemiological pattern of HIV is vital to monitor the spread, evaluate and implement HIV prevention strategies.

Methods

We analyzed 1,276 Ethiopian HIV-1 subtype C polymerase (pol sequences), including 144 newly generated sequences, collected from different parts of the country from 1986 to 2017. We employed state-of-art maximum likelihood and Bayesian phylodynamic analyses to comprehensively describe the evolutionary dynamics of the HIV-1 epidemic in Ethiopia. We used Bayesian phylodynamic models to estimate the dynamics of the effective population size (Ne) and reproductive numbers (Re) through time for the HIV epidemic in Ethiopia.

Results

Our analysis revealed that the Ethiopian HIV-1 epidemic originated from two independent introductions at the beginning of the 1970s and 1980s from eastern and southern African countries, respectively, followed by epidemic growth reaching its maximum in the early 1990s. We identified three large clusters with a majority of Ethiopian sequences. Phylodynamic analyses revealed that all three clusters were characterized by high transmission rates during the early epidemic, followed by a decline in HIV-1 transmissions after 1990. Re was high (4–6) during the earlier time of the epidemic but dropped significantly and remained low (Re &lt; 1) after the mid-1990. Similarly, with an expected shift in time, the effective population size (Ne) steadily increased until the beginning of 2000, followed by a decline and stabilization until recent years. The phylodynamic analyses corroborated the modeled UNAIDS incidence and prevalence estimates.

Conclusion

The rapid decline in the HIV epidemic took place a decade before introducing antiretroviral therapy in Ethiopia and coincided with early behavioral, preventive, and awareness interventions implemented in the country. Our findings highlight the importance of behavioral interventions and antiretroviral therapy scale-up to halt and maintain HIV transmissions at low levels (Re &lt; 1). The phylodynamic analyses provide epidemiological insights not directly available using standard surveillance and may inform the adjustment of public health strategies in HIV prevention in Ethiopia.

Genetic Diversity and Spatiotemporal Dynamics of Chikungunya Infections in Mexico during the Outbreak of 2014-2016

Chikungunya virus (CHIKV) is an alphavirus transmitted by Aedes mosquitoes, which causes Chikungunya fever. Three CHIKV genotypes have been identified: West African, East-Central-South African and Asian. In 2014, CHIKV was detected for the first time in Mexico, accumulating 13,569 confirmed cases in the following three years. Studies on the molecular diversification of CHIKV in Mexico focused on limited geographic regions or investigated only one structural gene of the virus. To describe the dynamics of this outbreak, we analyzed 309 serum samples from CHIKV acute clinical cases from 15 Mexican states. Partial NSP3, E1, and E2 genes were sequenced, mutations were identified, and their genetic variability was estimated. The evolutionary relationship with CHIKV sequences sampled globally were analyzed. Our sequences grouped with the Asian genotype within the Caribbean lineage, suggesting that the Asian was the only circulating genotype during the outbreak. Three non-synonymous mutations (E2 S248F and NSP3 A437T and L451F) were present in our sequences, which were also identified in sequences of the Caribbean lineage and in one Philippine sequence. Based on the phylogeographic analysis, the viral spread was reconstructed, suggesting that after the introduction through the Mexican southern border (Chiapas), CHIKV dispersed to neighboring states before reaching the center and north of the country through the Pacific Ocean states and Quintana Roo. This is the first viral phylogeographic reconstruction in Mexico characterizing the CHIKV outbreak across the country.

Emergence of a Distinct Picobirnavirus Genotype Circulating in Patients Hospitalized with Acute Respiratory Illness

Picobirnaviruses (PBV) are found in a wide range of hosts and typically associated with gastrointestinal infections in immunocompromised individuals. Here, a divergent PBV genome was assembled from a patient hospitalized for acute respiratory illness (ARI) in Colombia. The RdRp protein branched with sequences previously reported in patients with ARI from Cambodia and China. Sputa from hospitalized individuals (n = 130) were screened by RT-qPCR which enabled detection and subsequent metagenomic characterization of 25 additional PBV infections circulating in Colombia and the US. Phylogenetic analysis of RdRp highlighted the emergence of two dominant lineages linked to the index case and Asian strains, which together clustered as a distinct genotype. Bayesian inference further established capsid and RdRp sequences as both significantly associated with ARI. Various respiratory-tropic pathogens were detected in PBV+ patients, yet no specific bacteria was common among them and four individuals lacked co-infections, suggesting PBV may not be a prokaryotic virus nor exclusively opportunistic, respectively. Competing models for the origin and transmission of this PBV genotype are presented that attempt to reconcile vectoring by a bacterial host with human pathogenicity. A high prevalence in patients with ARI, an ability to reassort, and demonstrated global spread indicate PBV warrant greater public health concern.

Spatiotemporal Associations and Molecular Evolution of Highly Pathogenic Avian Influenza A H7N9 Virus in China from 2017 to 2021

Highly pathogenic (HP) H7N9 avian influenza virus (AIV) emerged in China in 2016. HP H7N9 AIV caused at least 33 human infections and has been circulating in poultry farms continuously since wave 5. The genetic divergence, geographic patterns, and hemagglutinin adaptive and parallel molecular evolution of HP H7N9 AIV in China since 2017 are still unclear. Here, 10 new strains of HP H7N9 AIVs from October 2019 to April 2021 were sequenced. We found that HP H7N9 was primarily circulating in Northern China, particularly in the provinces surrounding the Bohai Sea (Liaoning, Hebei, and Shandong) since wave 6. Of note, HP H7N9 AIV phylogenies exhibit a geographical structure compatible with high levels of local transmission after unidirectional rapid geographical expansion towards the north of China in 2017. In addition, we showed that two major subclades were continually expanding with the viral population size undergoing a sharp increase after 2018 with an obvious seasonal tendency. Notably, the hemagglutinin gene showed signs of parallel evolution and positive selection. Our research sheds light on the current epidemiology, evolution, and diversity of HP H7N9 AIV that can help prevent and control the spreading of HP H7N9 AIV.

A Daily-Updated Database and Tools for Comprehensive SARS-CoV-2 Mutation-Annotated Trees

The vast scale of SARS-CoV-2 sequencing data has made it increasingly challenging to comprehensively analyze all available data using existing tools and file formats. To address this, we present a database of SARS-CoV-2 phylogenetic trees inferred with unrestricted public sequences, which we update daily to incorporate new sequences. Our database uses the recently proposed mutation-annotated tree (MAT) format to efficiently encode the tree with branches labeled with parsimony-inferred mutations, as well as Nextstrain clade and Pango lineage labels at clade roots. As of June 9, 2021, our SARS-CoV-2 MAT consists of 834,521 sequences and provides a comprehensive view of the virus' evolutionary history using public data. We also present matUtils-a command-line utility for rapidly querying, interpreting, and manipulating the MATs. Our daily-updated SARS-CoV-2 MAT database and matUtils software are available at http://hgdownload.soe.ucsc.edu/goldenPath/wuhCor1/UShER_SARS-CoV-2/ and https://github.com/yatisht/usher, respectively.

Ongoing Assessment of the Molecular Evolution of Peste Des Petits Ruminants Virus Continues to Question Viral Origins

Understanding the evolution of viral pathogens is critical to being able to define how viruses emerge within different landscapes. Host susceptibility, which is spread between different species and is a contributing factor to the subsequent epidemiology of a disease, is defined by virus detection and subsequent characterization. Peste des petits ruminants virus is a plague of small ruminant species that is a considerable burden to the development of sustainable agriculture across Africa and much of Asia. The virus has also had a significant impact on populations of endangered species in recent years, highlighting its significance as a pathogen of high concern across different regions of the globe. Here, we have re-evaluated the molecular evolution of this virus using novel genetic data to try and further resolve the molecular epidemiology of this disease. Viral isolates are genetically characterized into four lineages (I-IV), and the historic origin of these lineages is of considerable interest to the molecular evolution of the virus. Our re-evaluation of viral emergence using novel genome sequences has demonstrated that lineages I, II and IV likely originated in West Africa, in Senegal (I) and Nigeria (II and IV). Lineage III sequences predicted emergence in either East Africa (Ethiopia) or in the Arabian Peninsula (Oman and/or the United Arab Emirates), with a paucity of data precluding a more refined interpretation. Continual refinements of evolutionary emergence, following the generation of new data, is key to both understanding viral evolution from a historic perspective and informing on the ongoing genetic emergence of this virus.

Geography-Driven Evolution of Potato Virus A Revealed by Genetic Diversity Analysis of the Complete Genome

Potato virus A (PVA), a member of the genus Potyvirus, is an important potato pathogen that causes 30%-40% yield reduction to global potato production. Knowledge on the genetic structure and the evolutionary forces shaping the structure of this pathogen is limited but vital in developing effective management strategies. In this study, we investigated the population structure and molecular evolution of PVA by analyzing novel complete genomic sequences from Chinese isolates combined with available sequences from Europe, South America, Oceania, and North America. High nucleotide diversity was discovered among the populations studied. Pairwise F _ST values between geographical populations of PVA ranged from 0.22 to 0.46, indicating a significant spatial structure for this pathogen. Although purifying selection was detected at the majority of polymorphic sites, significant positive selection was identified in the P1, NIa, and NIb proteins, pointing to adaptive evolution of PVA. Further phylogeny-trait association analysis showed that the clustering of PVA isolates was significantly correlated with geographic regions, suggesting that geography-driven adaptation may be an important determinant of PVA diversification.

Dissemination of Mycobacterium abscessus via global transmission networks

Mycobacterium abscessus, a multidrug-resistant nontuberculous mycobacterium, has emerged as a major pathogen affecting people with cystic fibrosis (CF). Although originally thought to be acquired independently from the environment, most individuals are infected with one of several dominant circulating clones (DCCs), indicating the presence of global transmission networks of M. abscessus. How and when these clones emerged and spread globally is unclear. Here, we use evolutionary analyses of isolates from individuals both with and without CF to reconstruct the population history, spatiotemporal spread and recent transmission networks of the DCCs. We demonstrate synchronous expansion of six unrelated DCCs in the 1960s, a period associated with major changes in CF care and survival. Each of these clones has spread globally as a result of rare intercontinental transmission events. We show that the DCCs, but not environmentally acquired isolates, exhibit a specific smoking-associated mutational signature and that current transmission networks include individuals both with and without CF. We therefore propose that the DCCs initially emerged in non-CF populations but were then amplified and spread through the CF community. While individuals with CF are probably the most permissive host, non-CF individuals continue to play a key role in transmission networks and may facilitate long-distance transmission.

Genomic sequencing of SARS-CoV-2 in Rwanda reveals the importance of incoming travelers on lineage diversity

COVID-19 transmission rates are often linked to locally circulating strains of SARS-CoV-2. Here we describe 203 SARS-CoV-2 whole genome sequences analyzed from strains circulating in Rwanda from May 2020 to February 2021. In particular, we report a shift in variant distribution towards the emerging sub-lineage A.23.1 that is currently dominating. Furthermore, we report the detection of the first Rwandan cases of the B.1.1.7 and B.1.351 variants of concern among incoming travelers tested at Kigali International Airport. To assess the importance of viral introductions from neighboring countries and local transmission, we exploit available individual travel history metadata to inform spatio-temporal phylogeographic inference, enabling us to take into account infections from unsampled locations. We uncover an important role of neighboring countries in seeding introductions into Rwanda, including those from which no genomic sequences were available. Our results highlight the importance of systematic genomic surveillance and regional collaborations for a durable response towards combating COVID-19.

Emergence and adaptive evolution of influenza D virus

As a novel member of the Orthomyxoviridae, influenza D virus (IDV) was firstly isolated from swine. However, cattle were found to serve as its primary reservoir. The study of IDV emergence can shed light into the dynamics of zoonotic infections and interspecies transmission. Although there is an increasing number of strains and sequenced IDV strains, their origin, epidemiology and evolutionary dynamics remain unclear. In this study, we reconstruct the diversity and evolutionary dynamics of IDVs. Molecular detection of swine tissue samples shows that six IDV positive samples were identified in the Eastern China. Phylogenetic analyses suggest three major IDV lineages designated as D/Japan, D/OK and D/660 as well as intermediate lineages. IDVs show strong association with geographical location indicating a high level of local transmission, which suggests IDVs tend to establish a local lineage of in situ evolution. In addition, the D/OK lineage widely circulates in swine in Eastern China, and all of the Chinese virus isolates form a distinct sub-clade (D/China sub-lineage). Furthermore, we identified important amino acids in the HEF gene under positive selection that might affect its receptor binding cavity relevant for its broader cell tropism. The combined results highlight that more attention should be paid to the potential threat of IDV to livestock and farming in China.

The spatial-temporal dynamics of respiratory syncytial virus infections across the east-west coasts of Australia during 2016-17

Respiratory syncytial virus (RSV) is an important human respiratory pathogen. In temperate regions, a distinct seasonality is observed, where peaks of infections typically occur in early winter, often preceding the annual influenza season. Infections are associated with high rates of morbidity and mortality and in some populations exceed that of influenza. Two subtypes, RSV-A and RSV-B, have been described, and molecular epidemiological studies have shown that both viruses mostly co-circulate. This trend also appears to be the case for Australia; however, previous genomic studies have been limited to cases from one Eastern state—New South Wales. As such, the broader spatial patterns and viral traffic networks across the continent are not known. Here, we conducted a whole-genome study of RSV comparing strains across eastern and Western Australia during the period January 2016 to June 2017. In total, 96 new RSV genomes were sequenced, compiled with previously generated data, and examined using a phylodynamic approach. This analysis revealed that both RSV-A and RSV-B strains were circulating, and each subtype was dominated by a single genotype, RSV-A ON1-like and RSV-B BA10-like viruses. Some geographical clustering was evident in strains from both states with multiple distinct sub-lineages observed and relatively low mixing across jurisdictions, suggesting that endemic transmission was likely seeded from imported, unsampled locations. Overall, the RSV phylogenies reflected a complex pattern of interactions across multiple epidemiological scales from fluid virus traffic across global and regional networks to fine-scale local transmission events.