Ecological divergence of wild birds drives avian influenza spillover and global spread.

The diversity of influenza A viruses (IAV) is primarily hosted by two highly divergent avian orders: Anseriformes (ducks, swans and geese) and Charadriiformes (gulls, terns and shorebirds). Studies of IAV have historically focused on Anseriformes, specifically dabbling ducks, overlooking the diversity of hosts in nature, including gull and goose species that have successfully adapted to human habitats. This study sought to address this imbalance by characterizing spillover dynamics and global transmission patterns of IAV over 10 years at greater taxonomic resolution than previously considered. Furthermore, the circulation of viral subtypes in birds that are either host-adapted (low pathogenic H13, H16) or host-generalist (highly pathogenic avian influenza-HPAI H5) provided a unique opportunity to test and extend models of viral evolution. Using Bayesian phylodynamic modelling we uncovered a complex transmission network that relied on ecologically divergent bird hosts. The generalist subtype, HPAI H5 was driven largely by wild geese and swans that acted as a source for wild ducks, gulls, land birds, and domestic geese. Gulls were responsible for moving HPAI H5 more rapidly than any other host, a finding that may reflect their long-distance, pelagic movements and their immuno-naïve status against this subtype. Wild ducks, long viewed as primary hosts for spillover, occupied an optimal space for viral transmission, contributing to geographic expansion and rapid dispersal of HPAI H5. Evidence of inter-hemispheric dispersal via both the Pacific and Atlantic Rims was detected, supporting surveillance at high latitudes along continental margins to achieve early detection. Both neutral (geographic expansion) and non-neutral (antigenic selection) evolutionary processes were found to shape subtype evolution which manifested as unique geographic hotspots for each subtype at the global scale. This study reveals how a diversity of avian hosts contribute to viral spread and spillover with the potential to improve surveillance in an era of rapid global change.

Accelerating Bayesian inference of dependency between mixed-type biological traits.

Inferring dependencies between mixed-type biological traits while accounting for evolutionary relationships between specimens is of great scientific interest yet remains infeasible when trait and specimen counts grow large. The state-of-the-art approach uses a phylogenetic multivariate probit model to accommodate binary and continuous traits via a latent variable framework, and utilizes an efficient bouncy particle sampler (BPS) to tackle the computational bottleneck-integrating many latent variables from a high-dimensional truncated normal distribution. This approach breaks down as the number of specimens grows and fails to reliably characterize conditional dependencies between traits. Here, we propose an inference pipeline for phylogenetic probit models that greatly outperforms BPS. The novelty lies in 1) a combination of the recent Zigzag Hamiltonian Monte Carlo (Zigzag-HMC) with linear-time gradient evaluations and 2) a joint sampling scheme for highly correlated latent variables and correlation matrix elements. In an application exploring HIV-1 evolution from 535 viruses, the inference requires joint sampling from an 11,235-dimensional truncated normal and a 24-dimensional covariance matrix. Our method yields a 5-fold speedup compared to BPS and makes it possible to learn partial correlations between candidate viral mutations and virulence. Computational speedup now enables us to tackle even larger problems: we study the evolution of influenza H1N1 glycosylations on around 900 viruses. For broader applicability, we extend the phylogenetic probit model to incorporate categorical traits, and demonstrate its use to study Aquilegia flower and pollinator co-evolution.

Sequential viral introductions and spread of BA.1 across Pakistan provinces during the Omicron wave.

BACKGROUND: COVID-19 waves caused by specific SARS-CoV-2 variants have occurred globally at different times. We focused on Omicron variants to understand the genomic diversity and phylogenetic relatedness of SARS-CoV-2 strains in various regions of Pakistan. METHODS: We studied 276,525 COVID-19 cases and 1,031 genomes sequenced from December 2021 to August 2022. Sequences were analyzed and visualized using phylogenetic trees. RESULTS: The highest case numbers and deaths were recorded in Sindh and Punjab, the most populous provinces in Pakistan. Omicron variants comprised 93% of all genomes, with BA.2 (32.6%) and BA.5 (38.4%) predominating. The first Omicron wave was associated with the sequential identification of BA.1 in Sindh, then Islamabad Capital Territory, Punjab, Khyber Pakhtunkhwa (KP), Azad Jammu Kashmir (AJK), Gilgit-Baltistan (GB) and Balochistan. Phylogenetic analysis revealed Sindh to be the source of BA.1 and BA.2 introductions into Punjab and Balochistan during early 2022. BA.4 was first introduced in AJK and BA.5 in Punjab. Most recent common ancestor (MRCA) analysis revealed relatedness between the earliest BA.1 genome from Sindh with Balochistan, AJK, Punjab and ICT, and that of first BA.1 from Punjab with strains from KPK and GB. CONCLUSIONS: Phylogenetic analysis provides insights into the introduction and transmission dynamics of the Omicron variant in Pakistan, identifying Sindh as a hotspot for viral dissemination. Such data linked with public health efforts can help limit surges of new infections.

Global dissemination of influenza A virus is driven by wild bird migration through arctic and subarctic zones.

Influenza A viruses (IAV) circulate endemically among many wild aquatic bird populations that seasonally migrate between wintering grounds in southern latitudes to breeding ranges along the perimeter of the circumpolar arctic. Arctic and subarctic zones are hypothesized to serve as ecologic drivers of the intercontinental movement and reassortment of IAVs due to high densities of disparate populations of long distance migratory and native bird species present during breeding seasons. Iceland is a staging ground that connects the East Atlantic and North Atlantic American flyways, providing a unique study system for characterizing viral flow between eastern and western hemispheres. Using Bayesian phylodynamic analyses, we sought to evaluate the viral connectivity of Iceland to proximal regions and how inter-species transmission and reassortment dynamics in this region influence the geographic spread of low and highly pathogenic IAVs. Findings demonstrate that IAV movement in the arctic and subarctic reflects wild bird migration around the perimeter of the circumpolar north, favouring short-distance flights between proximal regions rather than long distance flights over the polar interior. Iceland connects virus movement between mainland Europe and North America, consistent with the westward migration of wild birds from mainland Europe to Northeastern Canada and Greenland. Though virus diffusion rates were similar among avian taxonomic groups in Iceland, gulls play an outsized role as sinks of IAVs from other avian hosts prior to onward migration. These data identify patterns of virus movement in northern latitudes and inform future surveillance strategies related to seasonal and emergent IAVs with potential public health concern.

Human-Origin Influenza A(H3N2) Reassortant Viruses in Swine, Southeast Mexico.

The genetic diversity of influenza A viruses circulating in swine in Mexico complicates control efforts in animals and presents a threat to humans, as shown by influenza A(H1N1)pdm09 virus. To describe evolution of swine influenza A viruses in Mexico and evaluate strains for vaccine development, we sequenced the genomes of 59 viruses and performed antigenic cartography on strains from 5 regions. We found that genetic and antigenic diversity were particularly high in southeast Mexico because of repeated introductions of viruses from humans and swine in other regions in Mexico. We identified novel reassortant H3N2 viruses with genome segments derived from 2 different viruses that were independently introduced from humans into swine: pandemic H1N1 viruses and seasonal H3N2 viruses. The Mexico swine viruses are antigenically distinct from US swine lineages. Protection against these viruses is unlikely to be afforded by US virus vaccines and would require development of new vaccines specifically targeting these diverse strains.

The emergence and evolution of influenza A (H1α) viruses in swine in Canada and the United States.

Swine are a key reservoir host for influenza A viruses (IAVs), with the potential to cause global pandemics in humans. Gaps in surveillance in many of the world's largest swine populations impede our understanding of how novel viruses emerge and expand their spatial range in pigs. Although US swine are intensively sampled, little is known about IAV diversity in Canada's population of ~12 million pigs. By sequencing 168 viruses from multiple regions of Canada, our study reveals that IAV diversity has been underestimated in Canadian pigs for many years. Critically, a new H1 clade has emerged in Canada (H1α-3), with a two-amino acid deletion at H1 positions 146-147, that experienced rapid growth in Manitoba's swine herds during 2014-2015. H1α-3 viruses also exhibit a higher capacity to invade US swine herds, resulting in multiple recent introductions of the virus into the US Heartland following large-scale movements of pigs in this direction. From the Heartland, H1α-3 viruses have disseminated onward to both the east and west coasts of the United States, and may become established in Appalachia. These findings demonstrate how long-distance trading of live pigs facilitates the spread of IAVs, increasing viral genetic diversity and complicating pathogen control. The proliferation of novel H1α-3 viruses also highlights the need for expanded surveillance in a Canadian swine population that has long been overlooked, and may have implications for vaccine design.

Introduction, Evolution, and Dissemination of Influenza A Viruses in Exhibition Swine in the United States during 2009 to 2013.

The swine-human interface created at agricultural fairs, along with the generation of and maintenance of influenza A virus diversity in exhibition swine, presents an ongoing threat to public health. Nucleotide sequences of influenza A virus isolates collected from exhibition swine in Ohio (n = 262) and Indiana (n = 103) during 2009 to 2013 were used to investigate viral evolution and movement within this niche sector of the swine industry. Phylogenetic and Bayesian analyses were employed to identify introductions of influenza A virus to exhibition swine and study viral population dynamics. In 2013 alone, we identified 10 independent introductions of influenza A virus into Ohio and/or Indiana exhibition swine. Frequently, viruses from the same introduction were identified at multiple fairs within the region, providing evidence of rapid and widespread viral movement within the exhibition swine populations of the two states. While pigs moving from fair to fair to fair is possible in some locations, the concurrent detection of nearly identical strains at several fairs indicates that a common viral source was more likely. Importantly, we detected an association between the high number of human variant H3N2 (H3N2v) virus infections in 2012 and the widespread circulation of influenza A viruses of the same genotype in exhibition swine in Ohio fairs sampled that year. The extent of viral diversity observed in exhibition swine and the rapidity with which it disseminated across long distances indicate that novel strains of influenza A virus will continue to emerge and spread within exhibition swine populations, presenting an ongoing threat to humans. IMPORTANCE: Understanding the underlying population dynamics of influenza A viruses in commercial and exhibition swine is central to assessing the risk for human infections with variant viruses, including H3N2v. We used viral genomic sequences from isolates collected from exhibition swine during 2009 to 2013 to understand how the peak of H3N2v cases in 2012 relates to long-term trends in the population dynamics of pandemic viruses recently introduced into commercial and exhibition swine in the United States. The results of our spatial analysis underscore the key role of rapid viral dispersal in spreading multiple genetic lineages throughout a multistate network of agricultural fairs, providing opportunities for divergent lineages to coinfect, reassort, and generate new viral genotypes. The higher genetic diversity of genotypes cocirculating in exhibition swine since 2013 could facilitate the evolution of new reassortants, potentially with even greater ability to cause severe infections in humans or cause human-to-human transmission, highlighting the need for continued vigilance.

The evolutionary and transmission dynamics of HIV-1 CRF08_BC.

HIV-1 CRF08_BC is a significant subtype in China, though its origin and spread remain incompletely understood. Previous studies using partial genomic data have provided insights but lack comprehensive analysis. Here, we investigate the early evolutionary and spatiotemporal dynamics of HIV-1 CRF08_BC in China and Myanmar using near-complete genome sequences. We analyzed 28 near-complete HIV-1 CRF08_BC genomes from China and Myanmar (1997-2013). Phylogenetic, molecular clock, and Bayesian discrete trait analyses were performed to infer the virus's origin, spread, and associated risk groups. Based on Bayesian time-scaled inference with the best-fitting combination of models determined by marginal likelihood estimation (MLE), we inferred the time to the most recent common ancestor (TMRCA) and evolutionary rate of HIV-1 CRF08_BC to be at 3 October 1991 (95% HPD: 22 February1989-27 November 1993) and 2.30 × 10-3 substitutions per site per year (95% HPD: 1.96 × 10-3-2.63 × 10-3), respectively. Our analysis suggests that HIV-1 CRF08_BC originated in Yunnan Province, China, among injecting drug users, and subsequently spread to other regions. This study provides valuable insights into the early dynamics of HIV-1 CRF08_BC through combined genomic and epidemiological data, which may inform effective prevention and mitigation efforts. However, the limited genomic data influenced the extent of our findings, and challenges in collecting accurate risk group information during surveillance were evident.

The emergence and circulation of human immunodeficiency virus (HIV)-1 subtype C.

Introduction. Human immunodeficiency virus (HIV)-1 subtype C is the most prevalent globally and is thought to have originated in non-human primates in the Democratic Republic of Congo.Hypothesis/Gap Statement. Although the global dominance of HIV-1 subtype C is well established, a thorough understanding of its evolutionary history and transmission dynamics across various risk populations remains elusive. The current knowledge is insufficient to fully capture the global diversification and dissemination of this subtype.Aim. We for the first time sought to investigate the global evolutionary history and spatiotemporal dynamics of HIV-1 subtype C using a selection of maximum-likelihood-based phylodynamic approaches on a total of 1210 near full-length genomic sequences sampled from 32 countries, collected in 4 continents, with sampling dates between 1986-2019 among various risk groups were analysed.Methodology. We subsampled the HIV-1 subtype C genomic datasets based on continent and risk group traits, and performed nucleotide substitution model selection analysis, maximum likelihood (ML) phylogenetic reconstruction, phylogenetic tree topology similarity analysis, temporal signal analysis and traced the timings of viral spread both geographically and by risk group.Results. Based on the phylodynamic analyses of four datasets (full1210, locrisk626, loc562 and risk393), we inferred the time to the most recent common ancestor (TMRCA) in the 1930s and an evolutionary rate of 0.0023 substitutions per site per year. The total number of introduction events of HIV-1 subtype C between continents and between risk groups is estimated to be 71 and 115, respectively. The largest number of introductions occurred from Africa to Europe (n=32), from not-recorded to heterosexual (n=40) and from heterosexual to not-recorded (n=51) risk groups.Conclusion. Our results emphasize that HIV subtype C has mainly spread from Africa to Europe, likely through heterosexual transmission.

Comparative evolution of influenza A virus H1 and H3 head and stalk domains across host species.

IMPORTANCE: For decades, researchers have studied the rapid evolution of influenza A viruses for vaccine design and as a useful model system for the study of host/parasite evolution. By performing an exhaustive analysis of hemagglutinin protein (HA) sequences from 49 lineages independently evolving in birds, swine, canines, equines, and humans over the last century, our work uncovers surprising features of HA evolution. In particular, the canine H3 stalk, unlike human H3 and H1 stalk domains, is not evolving slowly, suggesting that evolution in the stalk domain is not universally constrained across all host species. Therefore, a broader multi-host perspective on HA evolution may be useful during the evaluation and design of stalk-targeted vaccine candidates.

A comprehensive SARS-CoV-2 and COVID-19 review, Part 2: host extracellular to systemic effects of SARS-CoV-2 infection.

COVID-19, the disease caused by SARS-CoV-2, has caused significant morbidity and mortality worldwide. The betacoronavirus continues to evolve with global health implications as we race to learn more to curb its transmission, evolution, and sequelae. The focus of this review, the second of a three-part series, is on the biological effects of the SARS-CoV-2 virus on post-acute disease in the context of tissue and organ adaptations and damage. We highlight the current knowledge and describe how virological, animal, and clinical studies have shed light on the mechanisms driving the varied clinical diagnoses and observations of COVID-19 patients. Moreover, we describe how investigations into SARS-CoV-2 effects have informed the understanding of viral pathogenesis and provide innovative pathways for future research on the mechanisms of viral diseases.

Space radiation damage rescued by inhibition of key spaceflight associated miRNAs.

Our previous research revealed a key microRNA signature that is associated with spaceflight that can be used as a biomarker and to develop countermeasure treatments to mitigate the damage caused by space radiation. Here, we expand on this work to determine the biological factors rescued by the countermeasure treatment. We performed RNA-sequencing and transcriptomic analysis on 3D microvessel cell cultures exposed to simulated deep space radiation (0.5 Gy of Galactic Cosmic Radiation) with and without the antagonists to three microRNAs: miR-16-5p, miR-125b-5p, and let-7a-5p (i.e., antagomirs). Significant reduction of inflammation and DNA double strand breaks (DSBs) activity and rescue of mitochondria functions are observed after antagomir treatment. Using data from astronaut participants in the NASA Twin Study, Inspiration4, and JAXA missions, we reveal the genes and pathways implicated in the action of these antagomirs are altered in humans. Our findings indicate a countermeasure strategy that can potentially be utilized by astronauts in spaceflight missions to mitigate space radiation damage.

Optimizing ancestral trait reconstruction of large HIV Subtype C datasets through multiple-trait subsampling.

Large datasets along with sampling bias represent a challenge for phylodynamic reconstructions, particularly when the study data are obtained from various heterogeneous sources and/or through convenience sampling. In this study, we evaluate the presence of unbalanced sampled distribution by collection date, location, and risk group of human immunodeficiency virus Type 1 Subtype C using a comprehensive subsampling strategy and assess their impact on the reconstruction of the viral spatial and risk group dynamics using phylogenetic comparative methods. Our study shows that a most suitable dataset for ancestral trait reconstruction can be obtained through subsampling by all available traits, particularly using multigene datasets. We also demonstrate that sampling bias is inflated when considerable information for a given trait is unavailable or of poor quality, as we observed for the trait risk group. In conclusion, we suggest that, even if traits are not well recorded, including them deliberately optimizes the representativeness of the original dataset rather than completely excluding them. Therefore, we advise the inclusion of as many traits as possible with the aid of subsampling approaches in order to optimize the dataset for phylodynamic analysis while reducing the computational burden. This will benefit research communities investigating the evolutionary and spatio-temporal patterns of infectious diseases.

Evolutionary and spatiotemporal analyses reveal multiple introductions and cryptic transmission of SARS-CoV-2 VOC/VOI in Malta.

IMPORTANCE: Our study provides insights into the evolution of the coronavirus disease 2019 (COVID-19) pandemic in Malta, a highly connected and understudied country. We combined epidemiological and phylodynamic analyses to analyze trends in the number of new cases, deaths, tests, positivity rates, and evolutionary and dispersal patterns from August 2020 to January 2022. Our reconstructions inferred 173 independent severe acute respiratory syndrome coronavirus 2 introductions into Malta from various global regions. Our study demonstrates that characterizing epidemiological trends coupled with phylodynamic modeling can inform the implementation of public health interventions to help control COVID-19 transmission in the community.

SARS-CoV-2 Orphan Gene ORF10 Contributes to More Severe COVID-19 Disease.

The orphan gene of SARS-CoV-2, ORF10, is the least studied gene in the virus responsible for the COVID-19 pandemic. Recent experimentation indicated ORF10 expression moderates innate immunity in vitro. However, whether ORF10 affects COVID-19 in humans remained unknown. We determine that the ORF10 sequence is identical to the Wuhan-Hu-1 ancestral haplotype in 95% of genomes across five variants of concern (VOC). Four ORF10 variants are associated with less virulent clinical outcomes in the human host: three of these affect ORF10 protein structure, one affects ORF10 RNA structural dynamics. RNA-Seq data from 2070 samples from diverse human cells and tissues reveals ORF10 accumulation is conditionally discordant from that of other SARS-CoV-2 transcripts. Expression of ORF10 in A549 and HEK293 cells perturbs immune-related gene expression networks, alters expression of the majority of mitochondrially-encoded genes of oxidative respiration, and leads to large shifts in levels of 14 newly-identified transcripts. We conclude ORF10 contributes to more severe COVID-19 clinical outcomes in the human host.

Core mitochondrial genes are down-regulated during SARS-CoV-2 infection of rodent and human hosts.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral proteins bind to host mitochondrial proteins, likely inhibiting oxidative phosphorylation (OXPHOS) and stimulating glycolysis. We analyzed mitochondrial gene expression in nasopharyngeal and autopsy tissues from patients with coronavirus disease 2019 (COVID-19). In nasopharyngeal samples with declining viral titers, the virus blocked the transcription of a subset of nuclear DNA (nDNA)-encoded mitochondrial OXPHOS genes, induced the expression of microRNA 2392, activated HIF-1α to induce glycolysis, and activated host immune defenses including the integrated stress response. In autopsy tissues from patients with COVID-19, SARS-CoV-2 was no longer present, and mitochondrial gene transcription had recovered in the lungs. However, nDNA mitochondrial gene expression remained suppressed in autopsy tissue from the heart and, to a lesser extent, kidney, and liver, whereas mitochondrial DNA transcription was induced and host-immune defense pathways were activated. During early SARS-CoV-2 infection of hamsters with peak lung viral load, mitochondrial gene expression in the lung was minimally perturbed but was down-regulated in the cerebellum and up-regulated in the striatum even though no SARS-CoV-2 was detected in the brain. During the mid-phase SARS-CoV-2 infection of mice, mitochondrial gene expression was starting to recover in mouse lungs. These data suggest that when the viral titer first peaks, there is a systemic host response followed by viral suppression of mitochondrial gene transcription and induction of glycolysis leading to the deployment of antiviral immune defenses. Even when the virus was cleared and lung mitochondrial function had recovered, mitochondrial function in the heart, kidney, liver, and lymph nodes remained impaired, potentially leading to severe COVID-19 pathology.

Readdressing the genetic diversity and taxonomy of the Mesoniviridae family, as well as its relationships with other nidoviruses and putative mesonivirus-like viral sequences.

Research on the recently established Mesoniviridae family (Order Nidovirales), RNA genome insect-specific viruses, has been steadily growing in the last decade. However, after the last detailed phylogenetic characterization of mesoniviruses in 2014, numerous new sequences, even in organisms other than mosquitos, have been identified and characterized. In this study, we analyzed nucleotide and protein sequences of mesoniviruses with a wide range of molecular tools including genetic distance, Shannon entropy, selective pressure analysis, polymorphism identification, principal coordinate analysis, likelihood mapping and phylodynamic reconstruction. We also sought to revaluate new mesoniviruses sequence positions within the family, proposing a taxonomic revision. The different sub-lineages of mosquito mesoniviruses sequences presented low sequence diversity and entropy, with incongruences to the existing taxonomy being found after an extensive phylogenetic characterization. High sequence discrepancy and differences in genome organization were found between mosquito mesoniviruses and other mesoniviruses, so their future classification, as other meso-like viruses that are found in other organisms, should be approached with caution. No evidence of frequent recombination was found, and mesonivirus genomes seem to evolve under strong purifying selection. Insufficient data by root-to-tip analysis did not yet allow for an adequate phylogeographic reconstruction.

Insect-specific viruses in the Parvoviridae family: Genetic lineage characterization and spatiotemporal dynamics of the recently established Brevihamaparvovirus genus.

The analysis of the viruses allocated to the recently established Brevihamaparvovirus genus (Parvoviridae family), which includes all previously known brevidensoviruses, has not yet been carried out on an extensive basis. As a result, no detailed genetic lineage characterization has ever been performed for this group of insect-specific viruses. Using a wide range of molecular tools, we have explored this taxon by calculating Shannon entropy values, intra- and inter-taxon genetic distances, analysed sequence polymorphisms, and evaluated selective pressures acting on the viral genome. While the calculated Brevihamaparvovirus mutation rates were within the range of those of other parvoviruses, their genomes look to be under strong purifying selection, and are also characterized by low diversity and entropy. Furthermore, even though recombination events are quite common among parvoviruses, no evidence of recombination (either intra or intergenic) was found in the Brevihamaparvoviruses sequences analyzed. An extended taxonomic analysis and reevaluation of existing Brevihamaparvoviruses sequences, many still unclassified, was performed using cut-off values defining NS1 identity between viral sequences from the Parvovirus family. Two existing genetic lineages, Dipteran Brevihamaparvovirus 1 and Dipteran Brevihamaparvovirus 2, were rearranged and the creation of a new one, Dipteran Brevihamaparvovirus 3, was suggested. Finally, despite the uncertainties associated with both the time estimates of the most recent common ancestors, which could span from twenty thousand years before the current era to way earlier (in the last century), and the dispersal routes proposed for Brevihamaparvoviruses sequences by phylodynamic reconstruction, the analyses here presented could help define how future studies should be conducted as more isolates continue to be identified in the future, and contribute to eliminating possible analytical biases.