Divergent Influenza-Like Viruses of Amphibians and Fish Support an Ancient Evolutionary Association

Variation of Site-Specific Glycosylation Profiles of Recombinant Influenza Glycoproteins

This work presents a detailed determination of site-specific N-glycan distributions of the recombinant influenza glycoproteins hemagglutinin (HA) and neuraminidase. Variation in glycosylation among recombinant glycoproteins is not predictable and can depend on details of the biomanufacturing process as well as details of protein structure. In this study, recombinant influenza proteins were analyzed from eight strains of four different suppliers. These include five HA and three neuraminidase proteins, each produced from a HEK293 cell line. Digestion was conducted using a series of complex multienzymatic methods designed to isolate glycopeptides containing single N-glycosylated sites. Site-specific glycosylation profiles of intact glycopeptides were produced using a recently developed method and comparisons were made using spectral similarity scores. Variation in glycan abundances and distribution was most pronounced between different strains of virus (similarity score = 383 out of 999), whereas digestion replicates and injection replicates showed relatively little variation (similarity score = 957). Notably, glycan distributions for homologous regions of influenza glycoprotein variants showed low variability. Due to the multiple possible sources of variation and inherent analytical difficulties in site-specific glycan determinations, variations were individually examined for multiple factors, including differences in supplier, production batch, protease digestion, and replicate measurement. After comparing all glycosylation distributions, four distinguishable classes could be identified for the majority of sites. Finally, attempts to identify glycosylation distributions on adjacent potential N-glycosylated sites of one HA variant were made. Only the second site (NnST) was found to be occupied using two rarely used proteases in proteomics, subtilisin and esperase, both of which did selectively cleave these adjacent sites.

Viromes of Antarctic fish resemble the diversity found at lower latitudes

Antarctica harbours some of the most isolated and extreme environments on Earth, concealing a largely unexplored and unique component of the global animal virosphere. To understand the diversity and evolutionary histories of viruses in these polar species, we determined the viromes of gill metatranscriptomes from 11 Antarctic fish species with 248 samples collected from the Ross Sea region spanning the Perciformes, Gadiformes, and Scorpaeniformes orders. The continent's shift southward and cooling temperatures >20 million years ago led to a reduction in biodiversity and subsequent radiation of some marine fauna, such as the notothenioid fishes. Despite decreased host species richness in polar regions, we revealed a surprisingly complex virome diversity in Ross Sea fish, with the types and numbers of viruses per host species and individuals sampled comparable to that of fish in warmer marine environments with higher host community diversity. We also observed a higher number of closely related viruses likely representing instances of recent and historic host-switching events among the Perciformes (all notothenioids) than in the Gadiformes, suggesting that rapid speciation events within this order generated closely related host species with few genetic barriers to cross-species transmission. Additionally, we identified novel genomic variation in an arenavirus with a split nucleoprotein sequence containing a stable helical structure, indicating potential adaptation of viral proteins to extreme temperatures. These findings enhance our understanding of virus evolution and virus-host interactions in response to environmental shifts, especially in less diverse ecosystems that are more vulnerable to the impacts of anthropogenic and climate changes.

A novel tamanavirus (Flaviviridae) of the European common frog (Rana temporaria) from the UK

Flavivirids are small, enveloped, positive-sense RNA viruses from the family Flaviviridae with genomes of ~9-13 kb. Metatranscriptomic analyses of metazoan organisms have revealed a diversity of flavivirus-like or flavivirid viral sequences in fish and marine invertebrate groups. However, no flavivirus-like virus has been identified in amphibians. To remedy this, we investigated the virome of the European common frog (Rana temporaria) in the UK, utilizing high-throughput sequencing at six catch locations. De novo assembly revealed a coding-complete virus contig of a novel flavivirid ~11.2 kb in length. The virus encodes a single ORF of 3456 aa and 5' and 3' untranslated regions (UTRs) of 227 and 666 nt, respectively. We named this virus Rana tamanavirus (RaTV), as BLASTp analysis of the polyprotein showed the closest relationships to Tamana bat virus (TABV) and Cyclopterus lumpus virus from Pteronotus parnellii and Cyclopterus lumpus, respectively. Phylogenetic analysis of the RaTV polyprotein compared to Flavivirus and Flavivirus-like members indicated that RaTV was sufficiently divergent and basal to the vertebrate Tamanavirus clade. In addition to the Mitcham strain, partial but divergent RaTV, sharing 95.64-97.39 % pairwise nucleotide identity, were also obtained from the Poole and Deal samples, indicating that RaTV is widespread in UK frog samples. Bioinformatic analyses of predicted secondary structures in the 3'UTR of RaTV showed the presence of an exoribonuclease-resistant RNA (xrRNA) structure standard in flaviviruses and TABV. To examine this biochemically, we conducted an in vitro Xrn1 digestion assay showing that RaTV probably forms a functional Xrn1-resistant xrRNA.

Evidence for an ancient aquatic origin of the RNA viral order Articulavirales

The emergence of previously unknown disease-causing viruses in mammals is in part the result of a long-term evolutionary process. Reconstructing the deep phylogenetic histories of viruses helps identify major evolutionary transitions and contextualizes the emergence of viruses in new hosts. We used a combination of total RNA sequencing and transcriptome data mining to extend the diversity and evolutionary history of the RNA virus order Articulavirales, which includes the influenza viruses. We identified instances of Articulavirales in the invertebrate phylum Cnidaria (including corals), constituting a novel and divergent family that we provisionally named the "Cnidenomoviridae." We further extended the evolutionary history of the influenza virus lineage by identifying four divergent, fish-associated influenza-like viruses, thereby supporting the hypothesis that fish were among the first hosts of influenza viruses. In addition, we substantially expanded the phylogenetic diversity of quaranjaviruses and proposed that this genus be reclassified as a family-the "Quaranjaviridae." Within this putative family, we identified a novel arachnid-infecting genus, provisionally named "Cheliceravirus." Notably, we observed a close phylogenetic relationship between the Crustacea- and Chelicerata-infecting "Quaranjaviridae" that is inconsistent with virus-host codivergence. Together, these data suggest that the Articulavirales has evolved over at least 600 million years, first emerging in aquatic animals. Importantly, the evolution of the Articulavirales was likely shaped by multiple aquatic-terrestrial transitions and substantial host jumps, some of which are still observable today.

Accumulated metagenomic studies reveal recent migration, whole genome evolution, and undiscovered diversity of orthomyxoviruses

IMPORTANCE

The number of known virus species has increased dramatically through metagenomic studies, which search genetic material sampled from a host for non-host genes. Here, we focus on an important viral family that includes influenza viruses, the Orthomyxoviridae, with over 100 recently discovered viruses infecting hosts from humans to fish. We find that one virus called Wǔhàn mosquito virus 6, discovered in mosquitoes in China, has spread across the globe very recently. Surface proteins used to enter cells show signs of rapid evolution in Wǔhàn mosquito virus 6 and its relatives which suggests an ability to infect vertebrate animals. We compute the rate at which new orthomyxovirus species discovered add evolutionary history to the tree of life, predict that many viruses remain to be discovered, and discuss what appropriately designed future studies can teach us about how diseases cross between continents and species.

450 million years in the making: mapping the evolutionary foundations of germinal centers

Germinal centers (GCs) are distinct microanatomical structures that form in the secondary lymphoid organs of endothermic vertebrates (i.e., mammals and some birds). Within GCs, B cells undergo a Darwinian selection process to identify clones which can respond to pathogen insult as well as affinity mature the B cell repertoire. The GC response ultimately generates memory B cells and bone marrow plasma cells which facilitate humoral immunological memory, the basis for successful vaccination programs. GCs have not been observed in the secondary lymphoid organs of ectothermic jawed vertebrates (i.e., fishes, reptiles, and amphibians). However, abundant research over the past decades has indicated these organisms can produce antigen specific B cell responses and some degree of affinity maturation. This review examines data demonstrating that the fundamentals of B cell selection may be more conserved across vertebrate phylogeny than previously anticipated. Further, research in both conventional mammalian model systems and comparative models raises the question of what evolutionary benefit GCs provide endotherms if they are seemingly unnecessary for generating the basic functional components of jawed vertebrate humoral adaptive immune responses.

Identification of RNA Virus-Derived RdRp Sequences in Publicly Available Transcriptomic Data Sets

RNA viruses are abundant and highly diverse and infect all or most eukaryotic organisms. However, only a tiny fraction of the number and diversity of RNA virus species have been catalogued. To cost-effectively expand the diversity of known RNA virus sequences, we mined publicly available transcriptomic data sets. We developed 77 family-level Hidden Markov Model profiles for the viral RNA-dependent RNA polymerase (RdRp)-the only universal "hallmark" gene of RNA viruses. By using these to search the National Center for Biotechnology Information Transcriptome Shotgun Assembly database, we identified 5,867 contigs encoding RNA virus RdRps or fragments thereof and analyzed their diversity, taxonomic classification, phylogeny, and host associations. Our study expands the known diversity of RNA viruses, and the 77 curated RdRp Profile Hidden Markov Models provide a useful resource for the virus discovery community.

Transcriptome Mining Expands Knowledge of RNA Viruses across the Plant Kingdom

Our current understanding of plant viruses stems largely from those affecting economically important plants. Yet plant species in cultivation represent a small and biased subset of the plant kingdom. Here, we describe virus diversity and abundance in 1,079 transcriptomes from species across the breadth of the plant kingdom (Archaeplastida) by analyzing open-source data from the 1000 Plant Transcriptomes Initiative (1KP). We identified 104 potentially novel viruses, of which 40% were single-stranded positive-sense RNA viruses across eight orders, including members of the Hepelivirales, Tymovirales, Cryppavirales, Martellivirales, and Picornavirales. One-third of the newly described viruses were double-stranded RNA viruses from the orders Durnavirales and Ghabrivirales. The remaining were negative-sense RNA viruses from the Rhabdoviridae, Aspiviridae, Yueviridae, and Phenuiviridae and the newly proposed Viridisbunyaviridae. Our analysis considerably expands the known host range of 13 virus families to include lower plants (e.g., Benyviridae and Secoviridae) and 4 virus families to include alga hosts (e.g., Tymoviridae and Chrysoviridae). More broadly, however, a cophylogeny analysis revealed that the evolutionary history of these families is largely driven by cross-species transmission events. The discovery of the first 30-kDa movement protein in a nonvascular plant suggests that the acquisition of plant virus movement proteins occurred prior to the emergence of the plant vascular system. Together, these data highlight that numerous RNA virus families are associated with older evolutionary plant lineages than previously thought and that the apparent scarcity of RNA viruses found in lower plants likely reflects a lack of investigation rather than their absence. IMPORTANCE Our knowledge of plant viruses is mainly limited to those infecting economically important host species. In particular, we know little about those viruses infecting basal plant lineages such as the ferns, lycophytes, bryophytes, and charophytes. To expand this understanding, we conducted a broad-scale viral survey of species across the breadth of the plant kingdom. We found that basal plants harbor a wide diversity of RNA viruses, including some that are sufficiently divergent to likely compose a new virus family. The basal plant virome revealed offers key insights into the evolutionary history of core plant virus gene modules and genome segments. More broadly, this work emphasizes that the scarcity of viruses found in these species to date most likely reflects the limited research in this area.

Decoding the Role of Temperature in RNA Virus Infections

RNA viruses include respiratory viruses, such as coronaviruses and influenza viruses, as well as vector-borne viruses, like dengue and West Nile virus. RNA viruses like these encounter various environments when they copy themselves and spread from cell to cell or host to host. Ex vivo differences, such as geographical location and humidity, affect their stability and transmission, while in vivo differences, such as pH and host gene expression, impact viral receptor binding, viral replication, and the host immune response against the viral infection. A critical factor affecting RNA viruses both ex vivo and in vivo, and defining the outcome of viral infections and the direction of viral evolution, is temperature. In this minireview, we discuss the impact of temperature on viral replication, stability, transmission, and adaptation, as well as the host innate immune response. Improving our understanding of how RNA viruses function, survive, and spread at different temperatures will improve our models of viral replication and transmission risk analyses.

Structural and inhibitor sensitivity analysis of influenza B-like viral neuraminidases derived from Asiatic toad and spiny eel

Influenza virus neuraminidase (NA) is an important target for antiviral development because it plays a crucial role in releasing newly assembled viruses. Two unique influenza-like virus genomes were recently reported in the Wuhan Asiatic toad and Wuhan spiny eel. Their NA genes appear to be highly divergent from all known influenza NAs, raising key questions as to whether the Asiatic toad influenza-like virus NA (tNA) and spiny eel NA (eNA) have canonical NA activities and structures and whether they show sensitivity to NA inhibitors (NAIs). Here, we found that both tNA and eNA have neuraminidase activities. A detailed structural analysis revealed that tNA and eNA present similar overall structures to currently known NAs, with a conserved calcium binding site. Inhibition assays indicated that tNA is resistant to NAIs, while eNA is still sensitive to NAIs. E119 is conserved in canonical NAs. The P119E substitution in tNA can restore sensitivity to NAIs, and, in contrast, the E119P substitution in eNA decreased its sensitivity to NAIs. The structures of NA-inhibitor complexes further provide a detailed insight into NA-inhibitor interactions at the atomic level. Moreover, tNA and eNA have unique N-glycosylation sites compared with canonical NAs. Collectively, the structural features, NA activities, and sensitivities to NAIs suggest that fish- and amphibian-derived influenza-like viruses may circulate in these vertebrates. More attention should be paid to these influenza-like viruses because their NA molecules may play roles in the emergence of NAI resistance.

Revealing the uncharacterised diversity of amphibian and reptile viruses

Amphibians and non-avian reptiles represent a significant proportion of terrestrial vertebrates, however knowledge of their viruses is not proportional to their abundance. Many amphibians and reptiles have strict habitual environments and localised populations and are vulnerable to viral outbreaks and potential elimination as a result. We sought to identify viruses that were hidden in amphibian and reptile metatranscriptomic data by screening 235 RNA-sequencing datasets from a 122 species covering 25 countries. We identified 26 novel viruses and eight previously characterised viruses from fifteen different viral families. Twenty-five viruses had RNA genomes with identity to Arteriviridae, Tobaniviridae, Hantaviridae, Rhabdoviridae, Astroviridae, Arenaviridae, Hepeviridae, Picornaviridae, Orthomyxoviridae, Reoviridae, Flaviviridae and Caliciviridae. In addition to RNA viruses, we also screened datasets for DNA viral transcripts, which are commonly excluded from transcriptomic analysis. We identified ten DNA viruses with identity to Papillomaviridae, Parvoviridae, Circoviridae and Adomaviridae. With the addition of these viruses, we expand the global amphibian and reptile virome and identify new potentially pathogenic viruses that could challenge populations. We speculate that amphibian viruses often have simpler genomes than those in amniotes, as in the case of the Secondpapillomavirinae and Orthomyxoviridae viruses identified in this study. In addition, we find evidence of inter-family recombination in RNA viruses, and we also identify new members of the recombinant Adomaviridae family. Overall, we provide insights into the uncharacterised diversity of amphibian and reptile viruses with the aim of improving population management, treatment and conservation into the future.

Viruses of Atlantic Bonefish (Albula vulpes) in Florida and the Caribbean show geographic patterns consistent with population declines

Atlantic Bonefish (Albula vulpes) are economically important due to their popularity with recreational anglers. In the State of Florida, USA, bonefish population numbers declined by approximately 60% between the 1990s and 2015. Habitat loss, water quality impairment, chemical inputs, and other anthropogenic factors have been implicated as causes, but the role of pathogens has been largely overlooked, especially with respect to viruses. We used a metagenomic approach to identify and quantify viruses in the blood of 103 A. vulpes sampled throughout their Western Atlantic range, including populations in Florida that have experienced population declines and populations in Belize, Mexico, Puerto Rico, and The Bahamas that have remained apparently stable. We identified four viruses, all of which are members of families known to infect marine fishes (Flaviviridae, Iflaviridae, Narnaviridae, and Nodaviridae), but all of which were previously undescribed. Bonefish from Florida and Mexico had higher viral richness (numbers of distinct viruses per individual fish) than fish sampled from other areas, and bonefish from the Upper Florida Keys had the highest prevalence of viral infection (proportion of positive fish) than fish sampled from any other location. Bonefish from Florida also had markedly higher viral loads than fish sampled from any other area, both for a novel narnavirus and for all viruses combined. Bonefish viruses may be indicators of environmentally driven physiological and immunological compromise, causes of ill health, or both.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10641-022-01306-9.

Structural Analysis of Monomeric RNA-Dependent Polymerases Revisited

In the past few years, our understanding of the RNA virosphere has changed dramatically due to the growth and spurt of metagenomics, exponentially increasing the number of RNA viral sequences, and providing a better understanding of their range of potential hosts. As of today, the only conserved protein among RNA viruses appears to be the monomeric RNA-dependent RNA polymerase. This enzyme belongs to the right-hand DNA-and RNA polymerases, which also includes reverse transcriptases and eukaryotic replicative DNA polymerases. The ubiquity of this protein in RNA viruses makes it a unique evolutionary marker and an appealing broad-spectrum antiviral target. In this work pairwise structural comparisons of viral RdRps and RTs were performed, including tertiary structures that have been obtained in the last few years. The resulting phylogenetic tree shows that the RdRps from (+)ss- and dsRNA viruses might have been recruited several times throughout the evolution of mobile genetic elements. RTs also display multiple evolutionary routes. We have identified a structural core comprising the entire palm, a large moiety of the fingers and the N-terminal helices of the thumb domain, comprising over 300 conserved residues, including two regions that we have named the “knuckles” and the “hypothenar eminence”. The conservation of an helix bundle in the region preceding the polymerase domain confirms that (−)ss and dsRNA Reoviruses’ polymerases share a recent ancestor. Finally, the inclusion of DNA polymerases into our structural analyses suggests that monomeric RNA-dependent polymerases might have diverged from B-family polymerases.

Type B and type A influenza polymerases have evolved distinct binding interfaces to recruit the RNA polymerase II CTD

During annual influenza epidemics, influenza B viruses (IBVs) co-circulate with influenza A viruses (IAVs), can become predominant and cause severe morbidity and mortality. Phylogenetic analyses suggest that IAVs (primarily avian viruses) and IBVs (primarily human viruses) have diverged over long time scales. Identifying their common and distinctive features is an effective approach to increase knowledge about the molecular details of influenza infection. The virus-encoded RNA-dependent RNA polymerases (FluPolB and FluPolA) are PB1-PB2-PA heterotrimers that perform transcription and replication of the viral genome in the nucleus of infected cells. Initiation of viral mRNA synthesis requires a direct association of FluPol with the host RNA polymerase II (RNAP II), in particular the repetitive C-terminal domain (CTD) of the major RNAP II subunit, to enable "cap-snatching" whereby 5'-capped oligomers derived from nascent RNAP II transcripts are pirated to prime viral transcription. Here, we present the first high-resolution co-crystal structure of FluPolB bound to a CTD mimicking peptide at a binding site crossing from PA to PB2. By performing structure-based mutagenesis of FluPolB and FluPolA followed by a systematic investigation of FluPol-CTD binding, FluPol activity and viral phenotype, we demonstrate that IBVs and IAVs have evolved distinct binding interfaces to recruit the RNAP II CTD, despite the CTD sequence being highly conserved across host species. We find that the PB2 627 subdomain, a major determinant of FluPol-host cell interactions and IAV host-range, is involved in CTD-binding for IBVs but not for IAVs, and we show that FluPolB and FluPolA bind to the host RNAP II independently of the CTD. Altogether, our results suggest that the CTD-binding modes of IAV and IBV may represent avian- and human-optimized binding modes, respectively, and that their divergent evolution was shaped by the broader interaction network between the FluPol and the host transcriptional machinery.

Diversity and evolution of the animal virome

The COVID-19 pandemic has given the study of virus evolution and ecology new relevance. Although viruses were first identified more than a century ago, we likely know less about their diversity than that of any other biological entity. Most documented animal viruses have been sampled from just two phyla - the Chordata and the Arthropoda - with a strong bias towards viruses that infect humans or animals of economic and social importance, often in association with strong disease phenotypes. Fortunately, the recent development of unbiased metagenomic next-generation sequencing is providing a richer view of the animal virome and shedding new light on virus evolution. In this Review, we explore our changing understanding of the diversity, composition and evolution of the animal virome. We outline the factors that determine the phylogenetic diversity and genomic structure of animal viruses on evolutionary timescales and show how this impacts assessment of the risk of disease emergence in the short term. We also describe the ongoing challenges in metagenomic analysis and outline key themes for future research. A central question is how major events in the evolutionary history of animals, such as the origin of the vertebrates and periodic mass extinction events, have shaped the diversity and evolution of the viruses they carry.

The science of the host-virus network

Better methods to predict and prevent the emergence of zoonotic viruses could support future efforts to reduce the risk of epidemics. We propose a network science framework for understanding and predicting human and animal susceptibility to viral infections. Related approaches have so far helped to identify basic biological rules that govern cross-species transmission and structure the global virome. We highlight ways to make modelling both accurate and actionable, and discuss the barriers that prevent researchers from translating viral ecology into public health policies that could prevent future pandemics.

Functionality of the putative surface glycoproteins of the Wuhan spiny eel influenza virus

A panel of influenza virus-like sequences were recently documented in fish and amphibians. Of these, the Wuhan spiny eel influenza virus (WSEIV) was found to phylogenetically cluster with influenza B viruses as a sister clade. Influenza B viruses have been documented to circulate only in humans, with certain virus isolates found in harbor seals. It is therefore interesting that a similar virus was potentially found in fish. Here we characterize the putative hemagglutinin (HA) and neuraminidase (NA) surface glycoproteins of the WSEIV. Functionally, we show that the WSEIV NA-like protein has sialidase activity comparable to B/Malaysia/2506/2004 influenza B virus NA, making it a bona fide neuraminidase that is sensitive to NA inhibitors. We tested the functionality of the HA by addressing the receptor specificity, stability, preferential airway protease cleavage, and fusogenicity. We show highly specific binding to monosialic ganglioside 2 (GM2) and fusogenicity at a range of different pH conditions. In addition, we found limited antigenic conservation of the WSEIV HA and NA relative to the B/Malaysia/2506/2004 virus HA and NA. In summary, we perform a functional and antigenic characterization of the glycoproteins of WSEIV to assess if it is indeed a bona fide influenza virus potentially circulating in ray-finned fish.

Influenza B viruses: underestimated and overlooked

Influenza B viruses circulate globally every year causing respiratory disease with significant clinical and socio-economic impacts. IBV are considered exclusive human pathogens with no established animal reservoirs, which suggests with concerted effort it may be possible to eradicate this virus from human circulation. However, this requires a deeper understanding of IBV virology and immunology and the design of vaccines that induce universal immunity to antigenic variants of IBV.