Revealing population demographics with environmental RNA

The analyses of environmental DNA (eDNA) and environmental RNA (eRNA) released by organisms into their surrounding environment (water, soil and air) have emerged as powerful tools for monitoring biodiversity. While eDNA has been widely adopted for the non-invasive detection of species and characterization of community composition, the utilization of eRNA is still in its infancy. Due to its functional nature, eRNA holds intriguing potential for biodiversity monitoring offering new avenues of research beyond species detection. For example, conspecifics that are almost genetically identical can exhibit distinct transcriptomic differences depending on their life stage. In this issue of Molecular Ecology Resources, Parsley and Goldberg (2024) demonstrate, through a lab-validated field study, that eRNA can be used to detect distinct life stages of amphibians. This study elegantly demonstrates that eRNA can be used not only to detect invasive or endangered species but also to reveal population demographic information important for guiding effective conservation strategies.

A specific and sensitive droplet digital polymerase chain reaction assay for the detection of tilapia lake virus in fish tissue and environmental samples

Tilapia lake virus (TiLV) causes high mortality in farmed and wild tilapia in various countries. We developed a highly specific and sensitive droplet digital polymerase chain reaction (ddPCR) assay to detect and quantify TiLV. The ddPCR assay could detect the virus at a lower threshold than the reverse transcription-quantitative polymerase reaction (RT-qPCR) method, and the sensitivity of the ddPCR assay was 10-fold higher. The diagnostic sensitivity and specificity of the ddPCR assay were 100% and did not cross-react with tilapia tissues infected with Tilapia parvovirus, Infectious spleen and kidney necrosis virus, Aeromonas hydrophila, Streptococcus agalactiae, S. iniae and Francisella noatunensis. The assay reproducibility was demonstrated by a high correlation coefficient of 0.998, and the inter-assay coefficients of variability indicated that the ddPCR assay exhibited low variability within and between measurements. The detection limit of the TiLV ddPCR assay was 100 fg cDNA, which is equal to 3.3 copies of TiLV. Furthermore, the ddPCR assay could detect TiLV in mucus, water and infected tissue samples and the lowest copy number of TiLV detected in water samples by the ddPCR assay was 7.9 ± 0.99 copies/reaction The results of the clinical samples tested for TiLV revealed that the ddPCR assay had a relatively higher detection rate than the RT-qPCR method. Overall, the ddPCR method offers a highly promising approach for the absolute quantification of TiLV in carrier fish and samples from the environment with low viral concentrations.

Environmental RNA as a Noninvasive Tool for Assessing Toxic Effects in Fish: A Proof-of-concept Study Using Japanese Medaka Exposed to Pyrene

Although environmental RNA (eRNA) is emerging as a noninvasive tool to assess the health status of aquatic macroorganisms, the potential of eRNA in assessing chemical hazards remain largely untested. In this study, we investigated the ability of eRNA to detect changes in gene expression in Japanese medaka fish (Oryzias latipes) in response to sublethal pyrene exposure, as a model toxic chemical. We performed standardized acute toxicity tests and collected eRNA from tank water and RNA from fish tissue after 96 h of exposure. Our results showed that over 1000 genes were detected in eRNA and the sequenced read counts of these genes correlated with those in fish tissue (r = 0.50). Moreover, eRNA detected 86 differentially expressed genes in response to pyrene, some of which were shared by fish RNA, including the suppression of collagen fiber genes. These results suggest that eRNA has the potential to detect changes in gene expression in fish in response to environmental stressors without the need for sacrificing or causing pain to fish. However, we also found that the majority of sequenced reads of eRNA (>99%) were not mapped to the reference medaka genome and they originated from bacteria and fungi, resulting in low sequencing depth. In addition, eRNA, in particular nuclear genes, was highly degraded with a median transcript integrity number (TIN) of <20. These limitations highlight the need for future studies to improve the analytical methods of eRNA application.

Environmental nucleic acids: A field-based comparison for monitoring freshwater habitats using eDNA and eRNA

Nucleic acids released by organisms and isolated from environmental substrates are increasingly being used for molecular biomonitoring. While environmental DNA (eDNA) has received much attention, the potential of environmental RNA as a biomonitoring tool remains under-explored. Several recent studies using paired DNA and RNA metabarcoding of bulk samples suggest that RNA might better reflect "metabolically active" parts of the community. However, such studies mainly capture organismal eDNA and eRNA. For larger eukaryotes, isolation of extra-organismal RNA will be important, but viability needs to be examined in a field-based setting. In this study we evaluate (a) whether extra-organismal eRNA release from macroeukaryotes can be detected given its supposedly rapid degradation, and (b) if the same field collection methods for eDNA can be applied to eRNA. We collected eDNA and eRNA from water in lakes where fish community composition is well documented, enabling a comparison between the two nucleic acids in two different seasons with monitoring using conventional methods. We found that eRNA is released from macroeukaryotes and can be filtered from water and metabarcoded in a similar manner as eDNA to reliably provide species composition information. eRNA had a small but significantly greater true positive rate than eDNA, indicating that it correctly detects more species known to exist in the lakes. Given relatively small differences between the two molecules in describing fish community composition, we conclude that if eRNA provides significant advantages in terms of lability, it is a strong candidate to add to the suite of molecular monitoring tools.

Deciphering the molecular signal from past and alive bacterial communities in aquatic sedimentary archives

Lake sediments accumulate information on biological communities thus acting as natural archives. Traditionally paleolimnology has focussed on fossilized remains of organisms, however, many organisms do not leave fossil evidence, meaning major ecosystem components are missing from environmental reconstructions. Many paleolimnology studies now incorporate molecular methods, including investigating microbial communities using environmental DNA (eDNA), but there is uncertainty about the contribution of living organisms to molecular inventories. In the present study, we obtained DNA and RNA inventories from sediment spanning 700 years to investigate the contribution of past and active communities to the molecular signal from sedimentary archives. Additionally, a droplet digital PCR (ddPCR) targeting the 16S ribosomal RNA (16S rRNA) gene of the photosynthetic cyanobacterial genera Microcystis was used to explore if RNA signals were from legacy RNA. We posit that the RNA signal is a mixture of legacy RNA, dormant cells, living bacteria and modern-day trace level contaminants that were introduced during sampling and preferentially amplified. The presence of legacy RNA was confirmed by the detection of Microcystis in sediments aged to ~200 years ago. Recent comparisons between 16S rRNA gene metabarcoding and traditional paleo proxies showed that past changes in bacterial communities can be reconstructed from sedimentary archives. The recovery of RNA in the present study has provided new insights into the origin of these signals. However, caution is required during analysis and interpretation of 16S rRNA gene metabarcoding data especially in recent sediments were there are potentially active bacteria.

Environmental RNA: A Revolution in Ecological Resolution?

Current advancements in environmental RNA (eRNA) exploit its relatively fast turnover rate relative to environmental DNA (eDNA) to assess 'metabolically active' or temporally/spatially recent community diversity. However, this focus significantly underutilizes the trove of potential ecological information encrypted in eRNA. Here, we argue for pushing beyond current species-level eDNA detection capabilities by using eRNA to detect any organisms with unique eRNA profiles, potentially including different life-history stages, sexes, or even specific phenotypes within a species. We also discuss the future of eRNA as a means of assessing the physiological status of organisms and the ecological health of populations and communities, reflecting ecosystem-level conditions. We posit that eRNA has the potential to significantly improve the resolution of organism detection, biological monitoring, and biomonitoring applications in ecology.

High throughput shotgun sequencing of eRNA reveals taxonomic and derived functional shifts across a benthic productivity gradient

Benthic macrofauna is regularly used in monitoring programmes, however the vast majority of benthic eukaryotic biodiversity lies mostly in microscopic organisms, such as meiofauna (invertebrates < 1 mm) and protists, that rapidly responds to environmental change. These communities have traditionally been hard to sample and handle in the laboratory, but DNA sequencing has made such work less time consuming. While DNA sequencing captures both alive and dead organisms, environmental RNA (eRNA) better targets living organisms or organisms of recent origin in the environment. Here, we assessed the biodiversity of three known bioindicator microeukaryote groups (nematodes, foraminifera, and ciliates) in sediment samples collected at seven coastal sites along an organic carbon (OC) gradient. We aimed to investigate if eRNA shotgun sequencing can be used to simultaneously detect differences in (i) biodiversity of multiple microeukaryotic communities; and (ii) functional feeding traits of nematodes. Results showed that biodiversity was lower for nematodes and foraminifera in high OC (6.2%-6.9%), when compared to low OC sediments (1.2%-2.8%). Dissimilarity in community composition increased for all three groups between Low OC and High OC, as well as the classified feeding type of nematode genera (with more nonselective deposit feeders in high OC sediment). High relative abundant genera included nematode Sabatieria and foraminifera Elphidium in high OC, and Cryptocaryon-like ciliates in low OC sediments. Considering that future sequencing technologies are likely to decrease in cost, the use of eRNA shotgun sequencing to assess biodiversity of benthic microeukaryotes could be a powerful tool in recurring monitoring programmes.

Low Turnover of Soil Bacterial rRNA at Low Temperatures

Ribosomal RNA (rRNA) is used widely to investigate potentially active microorganisms in environmental samples, including soil microorganisms and other microbial communities that are subjected to pronounced seasonal variation in temperature. This raises a question about the turnover of intracellular microbial rRNA at environmentally relevant temperatures. We analyzed the turnover at four temperatures of RNA isolated from soil bacteria amended with ¹⁴C-labeled uridine. We found that the half-life of recently produced RNA increased from 4.0 days at 20°C to 15.8 days at 4°C, and 215 days at −4°C, while no degradation was detected at −18°C during a 1-year period. We discuss the implications of the strong temperature dependency of rRNA turnover for interpretation of microbiome data based on rRNA isolated from environmental samples.

A Transmissible RNA Pathway in Honey Bees

Systemic RNAi, initiated by double-stranded RNA (dsRNA) ingestion, has been reported in diverse invertebrates, including honey bees, demonstrating environmental RNA uptake that undermines homologous gene expression. However, the question why any organism would take up RNA from the environment has remained largely unanswered. Here, we report on horizontal RNA flow among honey bees mediated by secretion and ingestion of worker and royal jelly diets. We demonstrate that transmission of jelly-secreted dsRNA to larvae is biologically active and triggers gene knockdown that lasts into adulthood. Worker and royal jellies harbor differential naturally occurring RNA populations. Jelly RNAs corresponded to honey bee protein-coding genes, transposable elements, and non-coding RNA, as well as bacteria, fungi, and viruses. These results reveal an inherent property of honey bees to share RNA among individuals and generations. Our findings suggest a transmissible RNA pathway, playing a role in social immunity and signaling between members of the hive.

A Secreted RNA Binding Protein Forms RNA-Stabilizing Granules in the Honeybee Royal Jelly

RNA flow between organisms has been documented within and among different kingdoms of life. Recently, we demonstrated horizontal RNA transfer between honeybees involving secretion and ingestion of worker and royal jellies. However, how the jelly facilitates transfer of RNA is still unknown. Here, we show that worker and royal jellies harbor robust RNA-binding activity. We report that a highly abundant jelly component, major royal jelly protein 3 (MRJP-3), acts as an extracellular non-sequence-specific RNA-aggregating factor. Multivalent RNA binding stimulates higher-order assembly of MRJP-3 into extracellular ribonucleoprotein granules that protect RNA from degradation and enhance RNA bioavailability. These findings reveal that honeybees have evolved a secreted dietary RNA-binding factor to concentrate, stabilize, and share RNA among individuals. Our work identifies high-order ribonucleoprotein assemblies with functions outside cells and organisms.