Australian terrestrial environments harbour extensive RNA virus diversity

Australia is home to a diverse range of unique native fauna and flora. To address whether Australian ecosystems also harbour unique viruses, we performed meta-transcriptomic sequencing of 16 farmland and sediment samples taken from the east and west coasts of Australia. We identified 2460 putatively novel RNA viruses across 18 orders, the vast majority of which belonged to the microbe-associated phylum Lenarviricota. In many orders, such as the Nodamuvirales and Ghabrivirales, the novel viruses identified here comprised entirely new clades. Novel viruses also fell between established genera or families, such as in the Cystoviridae and Picornavirales, while highly divergent lineages were identified in the Sobelivirales and Ghabrivirales. Viral read abundance and alpha diversity were influenced by sampling site, soil type and land use, but not by depth from the surface. In sum, Australian soils and sediments are home to remarkable viral diversity, reflecting the biodiversity of local fauna and flora.

Meta-transcriptomic identification of groundnut RNA viruses in western Kenya and the novel detection of groundnut as a host for Cauliflower mosaic virus

BACKGROUND

Groundnut (Arachis hypogaea L.) is the 13th most important global crop grown throughout the tropical and subtropical regions of the world. One of the major constraints to groundnut production is viruses, which are also the most economically important and most abundant pathogens among cultivated legumes. Only a few studies have reported the characterization of RNA viruses in cultivated groundnuts in western Kenya, most of which deployed classical methods of detecting known viruses.

METHODS

We sampled twenty-one symptomatic and three asymptomatic groundnut leaf samples from farmers' fields in western Kenya. Total RNA was extracted from the samples followed by First-strand cDNA synthesis and sequencing on the Illumina HiSeq 2500 platform. After removing host and rRNA sequences, high-quality viral RNA sequences were de novo assembled and viral genomes annotated using the publicly available NCBI virus database. Multiple sequence alignment and phylogenetic analysis were done using MEGA X.

RESULTS

Bioinformatics analyses using as low as ∼3.5 million reads yielded complete and partial genomes for Cauliflower mosaic virus (CaMV), Cowpea polerovirus 2 (CPPV2), Groundnut rosette assistor virus (GRAV), Groundnut rosette virus (GRV), Groundnut rosette virus satellite RNA (satRNA) and Peanut mottle virus (PeMoV) falling within the species demarcation criteria. This is the first report of CaMV and the second report of CPPV2 on groundnut hosts in the world. Confirmation of the detected viruses was further verified through phylogenetic analyses alongside reported publicly available highly similar viruses. PeMoV was the only seed-borne virus reported.

CONCLUSION

Our findings demonstrate the power of Next Generation Sequencing in the discovery and identification of novel viruses in groundnuts. The detection of the new viruses indicates the complexity of virus diseases in groundnuts and would require more focus in future studies to establish the effect of the viruses as sole or mixed infections on the crop. The detection of PeMoV with potential origin from Malawi indicates the importance of seed certification and cross-boundary seed health testing.

The clinical outcome of COVID-19 is strongly associated with microbiome dynamics in the upper respiratory tract

OBJECTIVES

The respiratory tract is the portal of entry for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Although a variety of respiratory pathogens other than SARS-CoV-2 have been associated with severe cases of COVID-19 disease, the dynamics of the upper respiratory microbiota during disease the course of disease, and how they impact disease manifestation, remain uncertain.

METHODS

We collected 349 longitudinal upper respiratory samples from a cohort of 65 COVID-19 patients (cohort 1), 28 samples from 28 recovered COVID-19 patients (cohort 2), and 59 samples from 59 healthy controls (cohort 3). All COVID-19 patients originated from the earliest stage of the epidemic in Wuhan. Based on a modified clinical scale, the disease course was divided into five clinical disease phases (pseudotimes): "Healthy" (pseudotime 0), "Incremental" (pseudotime 1), "Critical" (pseudotime 2), "Complicated" (pseudotime 3), "Convalescent" (pseudotime 4), and "Long-term follow-up" (pseudotime 5). Using meta-transcriptomics, we investigated the features and dynamics of transcriptionally active microbes in the upper respiratory tract (URT) over the course of COVID-19 disease, as well as its association with disease progression and clinical outcomes.

RESULTS

Our results revealed that the URT microbiome exhibits substantial heterogeneity during disease course. Two clusters of microbial communities characterized by low alpha diversity and enrichment for multiple pathogens or potential pathobionts (including Acinetobacter and Candida) were associated with disease progression and a worse clinical outcome. We also identified a series of microbial indicators that classified disease progression into more severe stages. Longitudinal analysis revealed that although the microbiome exhibited complex and changing patterns during COVID-19, a restoration of URT microbiomes from early dysbiosis toward more diverse status in later disease stages was observed in most patients. In addition, a group of potential pathobionts were strongly associated with the concentration of inflammatory indicators and mortality.

CONCLUSION

This study revealed strong links between URT microbiome dynamics and disease progression and clinical outcomes in COVID-19, implying that the treatment of severe disease should consider the full spectrum of microbial pathogens present.

Multi-omics analysis revealed the selective enrichment of partial denitrifying bacteria for the stable coupling of partial-denitrification and anammox process under the influence of low strength magnetic field

The microbial consortium involving anaerobic ammonium oxidation (anammox) and partial denitrification (PD), known as PD-anammox, is an emerging energy-efficient and lower carbon nitrogen removal process from wastewater. However, maintaining a stable PD process by locking nitrate reduction until nitrite was challenging. This study established the first stable connection of anammox with constant nitrite generation by PD bacteria under a low-strength (1.3 mT) magnetic field (MF). When the nitrogen loading rate was 1.81 kg-N/m3/d, the nitrogen removal efficiency of the control reactor (R1) was 75%, lower than that of the experimental reactor (R2), which was 85%. The expression of Thauera and Zoogloea, potential PD bacteria was substantially lower in R1 (5.75% and 1.21%, respectively) than in R2 (10.25 and 6.61%, respectively), according to a meta-transcriptomic analysis. At the same time, the mRNA expression of anammox genera Candidatus Brocadia and Candidatus Kuenenia was 33.53% and 3.83% in R1 and 22.86% and 1.87% in R2. Moreover, carbon and nitrogen metabolism pathways were more abundant under the influence of low-strength MF. The selective enrichment of PD bacteria can be attributed to the increased expression of carbon metabolic pathways like the citrate cycle, glycolysis/gluconeogenesis, and pyruvate metabolism. Interestingly, the control reactor was dominated by a hydroxylamine-dependent anammox process while a low-strength MF-enhanced nitric-oxide-dependent anammox process. For successful anammox-centered nitrogen removal from wastewater, this study demonstrated that low-strength MF is a convenient and applicable technique to lock the nitrate reduction until nitrite.

Nationwide genomic surveillance reveals the prevalence and evolution of honeybee viruses in China

BACKGROUND

The economic and environmental value of honeybees has been severely challenged in recent years by the collapse of their colonies worldwide, often caused by outbreaks of infectious diseases. However, our understanding of the diversity, prevalence, and transmission of honeybee viruses is largely obscure due to a lack of large-scale and longitudinal genomic surveillance on a global scale.

RESULTS

We report the meta-transcriptomic sequencing of nearly 2000 samples of the two most important economic and widely maintained honeybee species, as well as an associated ectoparasite mite, collected across China during 2016-2019. We document the natural diversity and evolution of honeybee viruses in China, providing evidence that multiple viruses commonly co-circulate within individual bee colonies. We also expanded the genomic data for 12 important honeybee viruses and revealed novel genetic variants and lineages associated with China. We identified more than 23 novel viruses from the honeybee and mite viromes, with some exhibiting ongoing replication in their respective hosts. Together, these data provide additional support to the idea that mites are an important reservoir and spill-over host for honeybee viruses.

CONCLUSIONS

Our data show that honeybee viruses are more widespread, prevalent, and genetically diverse than previously realized. The information provided is important in mitigating viral infectious diseases in honeybees, in turn helping to maintain sustainable productive agriculture on a global scale. Video Abstract.

A total infectome approach to understand the etiology of infectious disease in pigs

BACKGROUND

The global pork industry is continuously affected by infectious diseases that can result in large-scale mortality, trade restrictions, and major reductions in production. Nevertheless, the cause of many infectious diseases in pigs remains unclear, largely because commonly used diagnostic tools fail to capture the full diversity of potential pathogens and because pathogen co-infection is common.

RESULTS

We used a meta-transcriptomic approach to systematically characterize the pathogens in 136 clinical cases representing different disease syndromes in pigs, as well as in 12 non-diseased controls. This enabled us to simultaneously determine the diversity, abundance, genomic information, and detailed epidemiological history of a wide range of potential pathogens. We identified 34 species of RNA viruses, nine species of DNA viruses, seven species of bacteria, and three species of fungi, including two novel divergent members of the genus Pneumocystis. While most of these pathogens were only apparent in diseased animals or were at higher abundance in diseased animals than in healthy animals, others were present in healthy controls, suggesting opportunistic infections. Importantly, most of the cases examined here were characterized by co-infection with more than two species of viral, bacterial, or fungal pathogens, some with highly correlated occurrence and abundance levels. Examination of clinical signs and necropsy results in the context of relevant pathogens revealed that a multiple-pathogen model was better associated with the data than a single-pathogen model was.

CONCLUSIONS

Our data demonstrate that most of the pig diseases examined were better explained by the presence of multiple rather than single pathogens and that infection with one pathogen can facilitate infection or increase the prevalence/abundance of another. Consequently, it is generally preferable to consider the cause of a disease based on a panel of co-infecting pathogens rather than on individual infectious agents. Video abstract.

Meta-transcriptomic analysis of the virome and microbiome of the invasive Indian myna (Acridotheres tristis) in Australia

Invasive species exert a serious impact on native fauna and flora and have become the target of eradication and management efforts worldwide. Invasive avian species can also be important pathogen reservoirs, although their viromes and microbiomes have rarely been studied. As one of the top 100 invasive pest species globally, the expansion of Indian mynas (Acridotheres tristis) into peri-urban and rural environments, in conjunction with increasing free-ranging avian agricultural practices, may increase the risk of microbial pathogens jumping species boundaries. Herein, we used a meta-transcriptomic approach to explore the microbes present in brain, liver and large intestine of 16 invasive Indian myna birds in Sydney, Australia. From this, we discovered seven novel viruses from the families Adenoviridae, Caliciviridae, Flaviviridae, Parvoviridae and Picornaviridae. Interestingly, each of the novel viruses identified shared less than 80% genomic similarity with their closest relatives from other avian species, indicative of a lack of detectable virus transmission between invasive mynas to native or domestic species. Of note, we also identified two coccidian protozoa, Isospora superbusi and Isospora greineri, from the liver and gut tissues of mynas. Overall, these data demonstrate that invasive mynas can harbor a diversity of viruses and other microorganisms such that ongoing pathogen surveillance in this species is warranted.

RNA virome abundance and diversity is associated with host age in a bird species

Despite the ongoing interest in virus discovery, little is known about the factors that shape communities of viruses within individual hosts. Here, we address how virus communities might be impacted by the age of the hosts they infect, using total RNA sequencing to reveal the RNA viromes of different age groups of Ruddy Turnstones (Arenaria interpres). From oropharyngeal and cloacal swabs we identified 14 viruses likely infecting birds, 11 of which were novel, including members of the Reoviridae, Astroviridae, and Picornaviridae. Strikingly, 12 viruses identified were from juvenile birds sampled in the first year of their life, compared to only two viruses in adult birds. Both viral abundance and alpha diversity were marginally higher in juvenile than adult birds. As well as informing studies of virus ecology, that host age might be associated with viral composition is an important consideration for the future surveillance of novel and emerging viruses.

Retrospective meta-transcriptomic identification of severe dengue in a traveller returning from Africa to Sweden, 1990

Pathogens associated with haemorrhagic fever commonly have zoonotic origins. The first documented imported case of likely viral severe haemorrhagic fever in Sweden occurred in 1990. Despite extensive study, no aetiological agent was identified. Following retrospective investigation with total RNA-sequencing of samples collected between 7 and 36 days from onset of symptoms we identified dengue virus 3 (DENV-3) and a human pegivirus (HPgV). We conclude that the patient likely suffered from haemorrhagic symptoms due to an atypical severe and undiagnosed dengue infection.

Meta-transcriptomic identification of Trypanosoma spp. in native wildlife species from Australia

Background

Wildlife species carry a remarkable diversity of trypanosomes. The detection of trypanosome infection in native Australian fauna is central to understanding their diversity and host-parasite associations. The implementation of total RNA sequencing (meta-transcriptomics) in trypanosome surveillance and diagnosis provides a powerful methodological approach to better understand the host species distribution of this important group of parasites.

Methods

We implemented a meta-transcriptomic approach to detect trypanosomes in a variety of tissues (brain, liver, lung, skin, gonads) sampled from native Australian wildlife, comprising four marsupials (koala, Phascolarctos cinereus; southern brown bandicoot, Isoodon obesulus; swamp wallaby, Wallabia bicolor; bare-nosed wombat, Vombatus ursinus), one bird (regent honeyeater, Anthochaera phrygia) and one amphibian (eastern dwarf tree frog, Litoria fallax). Samples corresponded to both clinically healthy and diseased individuals. Sequencing reads were de novo assembled into contigs and annotated. The evolutionary relationships among the trypanosomatid sequences identified were determined through phylogenetic analysis of 18S rRNA sequences.

Results

We detected trypanosome sequences in all six species of vertebrates sampled, with positive samples in multiple organs and tissues confirmed by PCR. Phylogenetic analysis indicated that the trypanosomes infecting marsupials were related to those previously detected in placental and marsupial mammals, while the trypanosome in the regent honeyeater grouped with avian trypanosomes. In contrast, we provide the first evidence for a trypanosome in the eastern dwarf tree frog that was phylogenetically distinct from those described in other amphibians.

Conclusions

To our knowledge, this is the first meta-transcriptomic analysis of trypanosomes in native Australian wildlife, expanding the known genetic diversity of these important parasites. We demonstrated that RNA sequencing is sufficiently sensitive to detect low numbers of Trypanosoma transcripts and from diverse hosts and tissues types, thereby representing an effective means to detect trypanosomes that are divergent in genome sequence.

Human pegivirus in brain tissue of a patient with encephalitis

We describe a case of meningoencephalitis in which meta-transcriptomic (RNA) sequencing detected human pegivirus (HPgV) in brain tissue, cerebrospinal fluid, and serum in the absence of other pathogens. This is the first detection of HPgV in antemortem brain tissue, although it is uncertain whether HPgV is responsible for the observed encephalitis.

The fungicide azoxystrobin promotes freshwater cyanobacterial dominance through altering competition

BACKGROUND

Sharp increases in food production worldwide are attributable to agricultural intensification aided by heavy use of agrochemicals. This massive use of pesticides and fertilizers in combination with global climate change has led to collateral damage in freshwater systems, notably an increase in the frequency of harmful cyanobacterial blooms (HCBs). The precise mechanisms and magnitude of effects that pesticides exert on HCBs formation and proliferation have received little research attention and are poorly constrained.

RESULTS

We found that azoxystrobin (AZ), a common strobilurin fungicide, can favor cyanobacterial growth through growth inhibition of eukaryotic competitors (Chlorophyta) and possibly by inhibiting cyanobacterial parasites (fungi) as well as pathogenic bacteria and viruses. Meta-transcriptomic analyses identified AZ-responsive genes and biochemical pathways in eukaryotic plankton and bacteria, potentially explaining the microbial effects of AZ.

CONCLUSIONS

Our study provides novel mechanistic insights into the intertwined effects of a fungicide and eutrophication on microbial planktonic communities and cyanobacterial blooms in a eutrophic freshwater ecosystem. This knowledge may prove useful in mitigating cyanobacteria blooms resulting from agricultural intensification.

Metatranscriptomics as a tool to identify fungal species and subspecies in mixed communities - a proof of concept under laboratory conditions

High-throughput sequencing (HTS) enables the generation of large amounts of genome sequence data at a reasonable cost. Organisms in mixed microbial communities can now be sequenced and identified in a culture-independent way, usually using amplicon sequencing of a DNA barcode. Bulk RNA-seq (metatranscriptomics) has several advantages over DNA-based amplicon sequencing: it is less susceptible to amplification biases, it captures only living organisms, and it enables a larger set of genes to be used for taxonomic identification. Using a model mock community comprising 17 fungal isolates, we evaluated whether metatranscriptomics can accurately identify fungal species and subspecies in mixed communities. Overall, 72.9% of the RNA transcripts were classified, from which the vast majority (99.5%) were correctly identified at the species level. Of the 15 species sequenced, 13 were retrieved and identified correctly. We also detected strain-level variation within the Cryptococcus species complexes: 99.3% of transcripts assigned to Cryptococcus were classified as one of the four strains used in the mock community. Laboratory contaminants and/or misclassifications were diverse, but represented only 0.44% of the transcripts. Hence, these results show that it is possible to obtain accurate species- and strain-level fungal identification from metatranscriptome data as long as taxa identified at low abundance are discarded to avoid false-positives derived from contamination or misclassifications. This study highlights both the advantages and current challenges in the application of metatranscriptomics in clinical mycology and ecological studies.

Meta-transcriptomic identification of hepatitis B virus in cerebrospinal fluid in patients with central nervous system disease

Determining the etiological basis of central nervous system (CNS) infections is inherently challenging, primarily due to the multi-etiological nature. Using RNA sequencing, we aimed to identify microbes present in cerebrospinal fluid (CSF) of two patients suffering CNS infection, previously diagnosed with Cryptococcus sp. and Streptococcus pneumoniae infection, respectively. After meta-transcriptomic analysis, and confirmation with real-time PCR, hepatitis B virus (HBV) was detected in the CSF of two patients diagnosed with CNS syndrome. Phylogenetic analysis of the partial HBV genomes from these patients showed that they belonged to genotypes B and C and clustered with other viruses of Asian origin. In countries with high levels of HBV endemicity, the virus is likely to be found in patients diagnosed with CNS infections, although whether it contributes to symptoms and pathology, or is simply a coincidental infection, is unknown and merits further investigation.

Identification of diverse arthropod associated viruses in native Australian fleas

Fleas are important vectors of zoonotic disease. However, little is known about the natural diversity and abundance of flea viruses, particularly in the absence of disease associations, nor the evolutionary relationships among those viruses found in different parasitic vector species. Herein, we present the first virome scale study of fleas, based on the meta-transcriptomic analysis of 52 fleas collected along the eastern coast of Australia. Our analysis revealed 18 novel RNA viruses belonging to nine viral families with diverse genome organizations, although the majority (72%) possessed single-stranded positive-sense genomes. Notably, a number of the viruses identified belonged to the same phylogenetic groups as those observed in ticks sampled at the same locations, although none were likely associated with mammalian infection. Overall, we identified high levels of genomic diversity and abundance of viruses in the flea species studied, and established that fleas harbor viruses similar to those seen to other vectors.

Meta-transcriptomics reveals a diverse antibiotic resistance gene pool in avian microbiomes

BACKGROUND

Antibiotic resistance is rendering common bacterial infections untreatable. Wildlife can incorporate and disperse antibiotic-resistant bacteria in the environment, such as water systems, which in turn serve as reservoirs of resistance genes for human pathogens. Anthropogenic activity may contribute to the spread of bacterial resistance cycling through natural environments, including through the release of human waste, as sewage treatment only partially removes antibiotic-resistant bacteria. However, empirical data supporting these effects are currently limited. Here we used bulk RNA-sequencing (meta-transcriptomics) to assess the diversity and expression levels of functionally viable resistance genes in the gut microbiome of birds with aquatic habits in diverse locations.

RESULTS

We found antibiotic resistance genes in birds from all localities, from penguins in Antarctica to ducks in a wastewater treatment plant in Australia. Comparative analysis revealed that birds feeding at the wastewater treatment plant carried the greatest resistance gene burden, including genes typically associated with multidrug resistance plasmids as the aac(6)-Ib-cr gene. Differences in resistance gene burden also reflected aspects of bird ecology, taxonomy, and microbial function. Notably, ducks, which feed by dabbling, carried a higher abundance and diversity of resistance genes than turnstones, avocets, and penguins, which usually prey on more pristine waters.

CONCLUSIONS

These transcriptome data suggest that human waste, even if it undergoes treatment, might contribute to the spread of antibiotic resistance genes to the wild. Differences in microbiome functioning across different bird lineages may also play a role in the antibiotic resistance burden carried by wild birds. In summary, we reveal the complex factors explaining the distribution of resistance genes and their exchange routes between humans and wildlife, and show that meta-transcriptomics is a valuable tool to access functional resistance genes in whole microbial communities.

Metagenomic discovery and co-infection of diverse wobbly possum disease viruses and a novel hepacivirus in Australian brushtail possums

BACKGROUND

Australian brushtail possums (Trichosurus vulpecula) are an introduced pest species in New Zealand, but native to Australia where they are protected for biodiversity conservation. Wobbly possum disease (WPD) is a fatal neurological disease of Australian brushtail possums described in New Zealand populations that has been associated with infection by the arterivirus (Arteriviridae) wobbly possum disease virus (WPDV-NZ). Clinically, WPD-infected possums present with chronic meningoencephalitis, choroiditis and multifocal neurological symptoms including ataxia, incoordination, and abnormal gait.

METHODS

We conducted a retrospective investigation to characterise WPD in native Australian brushtail possums, and used a bulk meta-transcriptomic approach (i.e. total RNA-sequencing) to investigate its potential viral aetiology. PCR assays were developed for case diagnosis and full genome recovery in the face of extensive genetic variation.

RESULTS

We identified genetically distinct lineages of arteriviruses from archival tissues of WPD-infected possums in Australia, termed wobbly possum disease virus AU1 and AU2. Phylogenetically, WPDV-AU1 and WPDV-AU2 shared only ~ 70% nucleotide similarity to each other and the WPDV-NZ strain, suggestive of a relatively ancient divergence. Notably, we also identified a novel and divergent hepacivirus (Flaviviridae) - the first in a marsupial - in both WPD-infected and uninfected possums, indicative of virus co-infection.

CONCLUSIONS

We have identified marsupial-specific lineages of arteriviruses in mainland Australia that are genetically distinct from that in New Zealand, in some cases co-infecting animals with a novel hepacivirus. Our study provides new insight into the hidden genetic diversity of arteriviruses, the capacity for virus co-infection, and highlights the utility of meta-transcriptomics for disease investigation in a One Health context.

Microbiome and ecotypic adaption of Holcus lanatus (L.) to extremes of its soil pH range, investigated through transcriptome sequencing

BACKGROUND

Plants can adapt to edaphic stress, such as nutrient deficiency, toxicity and biotic challenges, by controlled transcriptomic responses, including microbiome interactions. Traditionally studied in model plant species with controlled microbiota inoculation treatments, molecular plant-microbiome interactions can be functionally investigated via RNA-Seq. Complex, natural plant-microbiome studies are limited, typically focusing on microbial rRNA and omitting functional microbiome investigations, presenting a fundamental knowledge gap. Here, root and shoot meta-transcriptome analyses, in tandem with shoot elemental content and root staining, were employed to investigate transcriptome responses in the wild grass Holcus lanatus and its associated natural multi-species eukaryotic microbiome. A full factorial reciprocal soil transplant experiment was employed, using plant ecotypes from two widely contrasting natural habitats, acid bog and limestone quarry soil, to investigate naturally occurring, and ecologically meaningful, edaphically driven molecular plant-microbiome interactions.

RESULTS

Arbuscular mycorrhizal (AM) and non-AM fungal colonization was detected in roots in both soils. Staining showed greater levels of non-AM fungi, and transcriptomics indicated a predominance of Ascomycota-annotated genes. Roots in acid bog soil were dominated by Phialocephala-annotated transcripts, a putative growth-promoting endophyte, potentially involved in N nutrition and ion homeostasis. Limestone roots in acid bog soil had greater expression of other Ascomycete genera and Oomycetes and lower expression of Phialocephala-annotated transcripts compared to acid ecotype roots, which corresponded with reduced induction of pathogen defense processes, particularly lignin biosynthesis in limestone ecotypes. Ascomycota dominated in shoots and limestone soil roots, but Phialocephala-annotated transcripts were insignificant, and no single Ascomycete genus dominated. Fusarium-annotated transcripts were the most common genus in shoots, with Colletotrichum and Rhizophagus (AM fungi) most numerous in limestone soil roots. The latter coincided with upregulation of plant genes involved in AM symbiosis initiation and AM-based P acquisition in an environment where P availability is low.

CONCLUSIONS

Meta-transcriptome analyses provided novel insights into H. lanatus transcriptome responses, associated eukaryotic microbiota functions and taxonomic community composition. Significant edaphic and plant ecotype effects were identified, demonstrating that meta-transcriptome-based functional analysis is a powerful tool for the study of natural plant-microbiome interactions.

Meta-transcriptomics and the evolutionary biology of RNA viruses

•

Meta-transcriptomics (bulk RNA-Seq) is a powerful new way to characterise viromes.

•

Meta-transcriptomic data are changing our understanding of virus evolution.

•

Invertebrates harbor an enormous phylogenetic and genomic diversity of RNA viruses.

•

Present sampling schemes have only revealed a miniscule fraction of the virosphere.

•

The new wealth of virus genomic data presents a major challenge to classification.

Metagenomics is transforming the study of virus evolution, allowing the full assemblage of virus genomes within a host sample to be determined rapidly and cheaply. The genomic analysis of complete transcriptomes, so-called meta-transcriptomics, is providing a particularly rich source of data on the global diversity of RNA viruses and their evolutionary history. Herein we review some of the insights that meta-transcriptomics has provided on the fundamental patterns and processes of virus evolution, with a focus on the recent discovery of a multitude of novel invertebrate viruses. In particular, meta-transcriptomics shows that the RNA virus world is more fluid than previously realized, with relatively frequent changes in genome length and structure. As well as having a transformative impact on studies of virus evolution, meta-transcriptomics presents major new challenges for virus classification, with the greater sampling of host taxa now filling many of the gaps on virus phylogenies that were previously used to define taxonomic groups. Given that most viruses in the future will likely be characterized using metagenomics approaches, and that we have evidently only sampled a tiny fraction of the total virosphere, we suggest that proposals for virus classification pay careful attention to the wonders unearthed in this new age of virus discovery.