The Functional RNA Identification (FRID) Pipeline: Identification of Potential Pseudoknot-Containing RNA Elements as Therapeutic Targets for SARS-CoV-2

The COVID-19 pandemic persists despite the development of effective vaccines. As such, it remains crucial to identify new targets for antiviral therapies. The causative virus of COVID-19, SARS-CoV-2, is a positive-sense RNA virus with RNA structures that could serve as therapeutic targets. One such RNA with established function is the frameshift stimulatory element (FSE), which promotes programmed ribosomal frameshifting. To accelerate identification of additional functional RNA elements, we introduce a novel computational approach termed the Functional RNA Identification (FRID) pipeline. The guiding principle of our pipeline, which uses established component programs as well as customized component programs, is that functional RNA elements have conserved secondary and pseudoknot structures that facilitate function. To assess the presence and conservation of putative functional RNA elements in SARS-CoV-2, we compared over 6,000 SARS-CoV-2 genomic isolates. We identified 22 functional RNA elements from the SARS-CoV-2 genome, 14 of which have conserved pseudoknots and serve as potential targets for small molecule or antisense oligonucleotide therapeutics. The FRID pipeline is general and can be applied to identify pseudoknotted RNAs for targeted therapeutics in genomes or transcriptomes from any virus or organism.

Experimental demonstration and pan-structurome prediction of climate-associated riboSNitches in Arabidopsis

Background

Genome-wide association studies (GWAS) aim to correlate phenotypic changes with genotypic variation. Upon transcription, single nucleotide variants (SNVs) may alter mRNA structure, with potential impacts on transcript stability, macromolecular interactions, and translation. However, plant genomes have not been assessed for the presence of these structure-altering polymorphisms or “riboSNitches.”

Results

We experimentally demonstrate the presence of riboSNitches in transcripts of two Arabidopsis genes, ZINC RIBBON 3 (ZR3) and COTTON GOLGI-RELATED 3 (CGR3), which are associated with continentality and temperature variation in the natural environment. These riboSNitches are also associated with differences in the abundance of their respective transcripts, implying a role in regulating the gene's expression in adaptation to local climate conditions. We then computationally predict riboSNitches transcriptome-wide in mRNAs of 879 naturally inbred Arabidopsis accessions. We characterize correlations between SNPs/riboSNitches in these accessions and 434 climate descriptors of their local environments, suggesting a role of these variants in local adaptation. We integrate this information in CLIMtools V2.0 and provide a new web resource, T-CLIM, that reveals associations between transcript abundance variation and local environmental variation.

Conclusion

We functionally validate two plant riboSNitches and, for the first time, demonstrate riboSNitch conditionality dependent on temperature, coining the term “conditional riboSNitch.” We provide the first pan-genome-wide prediction of riboSNitches in plants. We expand our previous CLIMtools web resource with riboSNitch information and with 1868 additional Arabidopsis genomes and 269 additional climate conditions, which will greatly facilitate in silico studies of natural genetic variation, its phenotypic consequences, and its role in local adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02656-4.