Hop stunt viroid infection induces heterochromatin reorganization.

Viroids are pathogenic noncoding RNAs that completely rely on their host molecular machinery to accomplish their life cycle. Several interactions between viroids and their host molecular machinery have been identified, including interference with epigenetic mechanisms such as DNA methylation. Despite this, whether viroids influence changes in other epigenetic marks such as histone modifications remained unknown. Epigenetic regulation is particularly important during pathogenesis processes because it might be a key regulator of the dynamism of the defense response. Here we have analyzed the changes taking place in Cucumis sativus (cucumber) facultative and constitutive heterochromatin during hop stunt viroid (HSVd) infection using chromatin immunoprecipitation (ChIP) of the two main heterochromatic marks: H3K9me2 and H3K27me3. We find that HSVd infection is associated with changes in both H3K27me3 and H3K9me2, with a tendency to decrease the levels of repressive epigenetic marks through infection progression. These epigenetic changes are connected to the transcriptional regulation of their expected targets, genes, and transposable elements. Indeed, several genes related to the defense response are targets of both epigenetic marks. Our results highlight another host regulatory mechanism affected by viroid infection, providing further information about the complexity of the multiple layers of interactions between pathogens/viroids and hosts/plants.

SARS-CoV-2 remodels the landscape of small non-coding RNAs with infection time and symptom severity.

The COVID-19 pandemic caused by the coronavirus SARS-CoV-2 has significantly impacted global health, stressing the necessity of basic understanding of the host response to this viral infection. In this study, we investigated how SARS-CoV-2 remodels the landscape of small non-coding RNAs (sncRNA) from a large collection of nasopharyngeal swab samples taken at various time points from patients with distinct symptom severity. High-throughput RNA sequencing analysis revealed a global alteration of the sncRNA landscape, with abundance peaks related to species of 21-23 and 32-33 nucleotides. Host-derived sncRNAs, including microRNAs (miRNAs), transfer RNA-derived small RNAs (tsRNAs), and small nucleolar RNA-derived small RNAs (sdRNAs) exhibited significant differential expression in infected patients compared to controls. Importantly, miRNA expression was predominantly down-regulated in response to SARS-CoV-2 infection, especially in patients with severe symptoms. Furthermore, we identified specific tsRNAs derived from Glu- and Gly-tRNAs as major altered elements upon infection, with 5' tRNA halves being the most abundant species and suggesting their potential as biomarkers for viral presence and disease severity prediction. Additionally, down-regulation of C/D-box sdRNAs and altered expression of tinyRNAs (tyRNAs) were observed in infected patients. These findings provide valuable insights into the host sncRNA response to SARS-CoV-2 infection and may contribute to the development of further diagnostic and therapeutic strategies in the clinic.

Integrative time-scale and multi-omics analysis of host responses to viroid infection.

Viroids are circular RNAs of minimal complexity compelled to subvert plant-regulatory networks to accomplish their infectious process. Studies focused on the response to viroid-infection have mostly addressed specific regulatory levels and considered specifics infection-times. Thus, much remains to be done to understand the temporal evolution and complex nature of viroid-host interactions. Here we present an integrative analysis of the temporal evolution of the genome-wide alterations in cucumber plants infected with hop stunt viroid (HSVd) by integrating differential host transcriptome, sRNAnome and methylome. Our results support that HSVd promotes the redesign of the cucumber regulatory-pathways predominantly affecting specific regulatory layers at different infection-phases. The initial response was characterised by a reconfiguration of the host-transcriptome by differential exon-usage, followed by a progressive transcriptional downregulation modulated by epigenetic changes. Regarding endogenous small RNAs, the alterations were limited and mainly occurred at the late stage. Significant host-alterations were predominantly related to the downregulation of transcripts involved in plant-defence mechanisms, the restriction of pathogen-movement and the systemic spreading of defence signals. We expect that these data constituting the first comprehensive temporal-map of the plant-regulatory alterations associated with HSVd infection could contribute to elucidate the molecular basis of the yet poorly known host-response to viroid-induced pathogenesis.

Occurrence of RNA post-transcriptional modifications in plant viruses and viroids and their correlation with structural and functional features.

Post-transcriptional modifications of RNA bases are widespread across all the tree of life and have been linked to RNA maturation, stability, and molecular interactions. RNA modifications have been extensively described in endogenous eukaryotic mRNAs, however, little is known about the presence of RNA modifications in plant viral and subviral RNAs. Here, we used a computational approach to infer RNA modifications in plant-pathogenic viruses and viroids using high-throughput annotation of modified ribonucleotides (HAMR), a software that predicts modified ribonucleotides using high-throughput RNA sequencing data. We analyzed datasets from representative members of different plant viruses and viroids and compared them to plant-endogenous mRNAs. Our approach was able to predict potential RNA chemical modifications (RCMs) in all analyzed pathogens. We found that both DNA and RNA viruses presented a wide range of RCM proportions while viroids had lowest values. Furthermore, we found that for viruses with segmented genomes, some genomic RNAs had a higher proportion of RCM. Interestingly, nuclear-replicating viroids showed most of the predicted modifications located in the pathogenesis region, pointing towards a possible functional role of RCMs in their infectious cycle. Thus, our results strongly suggest that plant viral and subviral RNAs might contain a variety of previously unreported RNA modifications, thus opening a new perspective in the multifaceted process of plant-pathogen interactions.

A circular RNA vector for targeted plant gene silencing based on an asymptomatic viroid.

Gene silencing for functional studies in plants has been largely facilitated by manipulating viral genomes with inserts from host genes to trigger virus-induced gene silencing (VIGS) against the corresponding mRNAs. However, viral genomes encode multiple proteins and can disrupt plant homeostasis by interfering with endogenous cell mechanisms. To try to circumvent this functional limitation, we have developed a silencing method based on the minimal autonomously-infectious nucleic acids currently known: viroids, which lack proven coding capability. The genome of Eggplant latent viroid, an asymptomatic viroid, was manipulated with insertions ranging between 21 and 42 nucleotides. Our results show that, although larger insertions might be tolerated, the maintenance of the secondary structure appears to be critical for viroid genome stability. Remarkably, these modified ELVd molecules are able to induce systemic infection promoting the silencing of target genes in eggplant. Inspired by the design of artificial microRNAs, we have developed a simple and standardized procedure to generate stable insertions into the ELVd genome capable of silencing a specific target gene. Analogously to VIGS, we have termed our approach viroid-induced gene silencing, and demonstrate that it is a promising tool for dissecting gene functions in eggplant.

Plant epigenome alterations: an emergent player in viroid-host interactions.

It is well known that viroids promote significant alterations at diverse host regulatory levels. However the mechanisms by which these tiny RNAs subvert endogenous regulatory networks remain a to a large extent unsolved question. In the last years diverse studies have revealed the existence of a close interplay between viroid infection and host DNA methylation, suggesting that the modulation of the endogenous transcriptional activity by epigenetic alterations of the host genome may emerge as a novel player in plant-viroid interactions. Here, we summarize the more relevant findings related to alteration of the host DNA methylome in response to viroid infection and discus the potential strategies which may be exploited by these non-conventional pathogenic RNAs to hijack and redesign the plant epigenome.

Combined Stress Conditions in Melon Induce Non-additive Effects in the Core miRNA Regulatory Network.

Climate change has been associated with a higher incidence of combined adverse environmental conditions that can promote a significant decrease in crop productivity. However, knowledge on how a combination of stresses might affect plant development is still scarce. MicroRNAs (miRNAs) have been proposed as potential targets for improving crop productivity. Here, we have combined deep-sequencing, computational characterization of responsive miRNAs and validation of their regulatory role in a comprehensive analysis of response of melon to several combinations of four stresses (cold, salinity, short day, and infection with a fungus). Twenty-two miRNA families responding to double and/or triple stresses were identified. The regulatory role of the differentially expressed miRNAs was validated by quantitative measurements of the expression of the corresponding target genes. A high proportion (ca. 60%) of these families (mainly highly conserved miRNAs targeting transcription factors) showed a non-additive response to multiple stresses in comparison with that observed under each one of the stresses individually. Among those miRNAs showing non-additive response to stress combinations, most interactions were negative, suggesting the existence of functional convergence in the miRNA-mediated response to combined stresses. Taken together, our results provide compelling pieces of evidence that the response to combined stresses cannot be easily predicted from the study individual stresses.

Might exogenous circular RNAs act as protein-coding transcripts in plants?

Circular RNAs (circRNAs) are regulatory molecules involved in the modulation of gene expression. Although originally assumed as non-coding RNAs, recent studies have evidenced that animal circRNAs can act as translatable transcripts. The study of plant-circRNAs is incipient, and no autonomous coding plant-circRNA has been described yet. Viroids are the smallest plant-pathogenic circRNAs known to date. Since their discovery 50 years ago, viroids have been considered valuable systems for the study of the structure-function relationships in RNA, essentially because they have not been shown to have coding capacity. We used two pathogenic circRNAs (Hop stunt viroid and Eggplant latent viroid) as experimental tools to explore the coding potential of plant-circRNAs. Our work supports that the analysed viroids contain putative ORFs able to encode peptides carrying subcellular localization signals coincident with the corresponding replication-specific organelle. Bioassays in well-established hosts revealed that mutations in these ORFs diminish their biological efficiency. Interestingly, circular forms of HSVd and ELVd were found to co-sediment with polysomes, revealing their physical interaction with the translational machinery of the plant cell. Based on this evidence we hypothesize about the possibility that plant circRNAs in general, and viroids in particular, can act, under certain cellular conditions, as non-canonical translatable transcripts.

Mapping of Functional Subdomains in the atALKBH9B m6A-Demethylase Required for Its Binding to the Viral RNA and to the Coat Protein of Alfalfa Mosaic Virus.

N 6-methyladenosine (m6A) modification is a dynamically regulated RNA modification that impacts many cellular processes and pathways. This epitranscriptomic methylation relies on the participation of RNA methyltransferases (referred to as "writers") and demethylases (referred to as "erasers"), respectively. We previously demonstrated that the Arabidopsis thaliana protein atALKBH9B showed m6A-demethylase activity and interacted with the coat protein (CP) of alfalfa mosaic virus (AMV), causing a profound impact on the viral infection cycle. To dissect the functional activity of atALKBH9B in AMV infection, we performed a protein-mapping analysis to identify the putative domains required for regulating this process. In this context, the mutational analysis of the protein revealed that the residues between 427 and 467 positions are critical for in vitro binding to the AMV RNA. The atALKBH9B amino acid sequence showed intrinsically disordered regions (IDRs) located at the N-terminal part delimiting the internal AlkB-like domain and at the C-terminal part. We identified an RNA binding domain containing an RGxxxRGG motif that overlaps with the C-terminal IDR. Moreover, bimolecular fluorescent experiments allowed us to determine that residues located between 387 and 427 are critical for the interaction with the AMV CP, which should be critical for modulating the viral infection process. Finally, we observed that atALKBH9B deletions of either N-terminal 20 residues or the C-terminal's last 40 amino acids impede their accumulation in siRNA bodies. The involvement of the regions responsible for RNA and viral CP binding and those required for its localization in stress granules in the viral cycle is discussed.

Symptom Severity, Infection Progression and Plant Responses in Solanum Plants Caused by Three Pospiviroids Vary with the Inoculation Procedure.

Infectious viroid clones consist of dimeric cDNAs used to generate transcripts which mimic the longer-than-unit replication intermediates. These transcripts can be either generated in vitro or produced in vivo by agro-inoculation. We have designed a new plasmid, which allows both inoculation methods, and we have compared them by infecting Solanum lycopersicum and Solanum melongena with clones of Citrus exocortis virod (CEVd), Tomato chlorotic dwarf viroid (TCDVd), and Potato spindle tuber viroid (PSTVd). Our results showed more uniform and severe symptoms in agro-inoculated plants. Viroid accumulation and the proportion of circular and linear forms were different depending on the host and the inoculation method and did not correlate with the symptoms, which correlated with an increase in PR1 induction, accumulation of the defensive signal molecules salicylic (SA) and gentisic (GA) acids, and ribosomal stress in tomato plants. The alteration in ribosome biogenesis was evidenced by both the upregulation of the tomato ribosomal stress marker SlNAC082 and the impairment in 18S rRNA processing, pointing out ribosomal stress as a novel signature of the pathogenesis of nuclear-replicating viroids. In conclusion, this updated binary vector has turned out to be an efficient and reproducible method that will facilitate the studies of viroid-host interactions.

Hop stunt viroid: A polyphagous pathogenic RNA that has shed light on viroid-host interactions.

TAXONOMY: Hop stunt viroid (HSVd) is the type species of the genus Hostuviroid (family Pospiviroidae). The other species of this genus is Dahlia latent viroid, which presents an identical central conserved region (CCR) but lacks other structural hallmarks present in Hop stunt viroid. HSVd replication occurs in the nucleus through an asymmetric rolling-circle model as in the other members of the family Pospiviroidae, which also includes the genera Pospiviroid, Cocadviroid, Apscaviroid, and Coleoviroid. PHYSICAL PROPERTIES: Hop stunt viroid consists of a single-stranded, circular RNA of 295-303 nucleotides depending on isolates and sequence variants. The most stable secondary structure is a rod-like or quasi-rod-like conformation with two characteristic domains: a CCR and a terminal conserved hairpin similar to that of cocadviroids. HSVd lacks a terminal conserved region. HOSTS AND SYMPTOMS: HSVd infects a very broad range of natural hosts and has been reported to be the causal agent of five different diseases (citrus cachexia, cucumber pale fruit, peach and plum apple apricot distortion, and hop stunt). It is distributed worldwide. TRANSMISSION: HSVd is transmitted mechanically and by seed.

HPG-DHunter: an ultrafast, friendly tool for DMR detection and visualization.

BACKGROUND: Software tools for analyzing DNA methylation do not provide graphical results which can be easily identified, but huge text files containing the alignment of the samples and their methylation status at a resolution of base pairs. There have been proposed different tools and methods for finding Differentially Methylated Regions (DMRs) among different samples, but the execution time required by these tools is large, and the visualization of their results is far from being interactive. Additionally, these methods show more accurate results when identifying simulated DM regions that are long and have small within-group variation, but they have low concordance when used with real datasets, probably due to the different approaches they use for DMR identification. Thus, a tool which automatically detects DMRs among different samples and interactively visualizes DMRs at different scales (from a bunch to ten of millions of DNA locations) can be the key for shortening the DNA methylation analysis process in many studies. RESULTS: In this paper, we propose a software tool based on the wavelet transform. This mathematical tool allows the fast automatic DMR detection by simple comparison of different signals at different resolution levels. Also, it allows an interactive visualization of the DMRs found at different resolution levels. The tool is publicly available at https://grev-uv.github.io/ , and it is part of a complete suite of tools which allow to carry out the complete process of DNA alignment and methylation analysis, creation of methylation maps of the whole genome, and the detection and visualization of DMRs between different samples. CONCLUSIONS: The validation of the developed software tool shows similar concordance with other well-known and extended tools when used with real and synthetic data. The batch mode of the tool is capable of automatically detecting the existing DMRs for half (twelve) of the human chromosomes between two sets of six samples (whose.csv files after the alignment and mapping procedures have an aggregated size of 108 Gigabytes) in around three hours and a half. When compared to other well-known tools, HPG-DHunter only requires around 15% of the execution time required by other tools for detecting the DMRs.

Highly efficient construction of infectious viroid-derived clones.

BACKGROUND: Viroid research generally relies on infectious cDNA clones that consist of dimers of the entire viroid sequence. At present, those dimers are generated by self-ligation of monomeric cDNA, a strategy that presents several disadvantages: (i) low efficiency, (ii) it is a non-oriented reaction requiring tedious screenings and (iii) additional steps are required for cloning into a binary vector for agroinfiltration or for in vitro RNA production. RESULTS: We have developed a novel strategy for simultaneous construction of a viroid dimeric cDNA and cloning into a multipurpose binary vector ready for agroinfiltration or in vitro transcription. The assembly is based on IIs restriction enzymes and positive selection and supposes a universal procedure for obtaining infectious clones of a viroid independently of its sequence, with a high efficiency. Thus, infectious clones of one viroid of each family were obtained and its infectivity was analyzed by molecular hybridization. CONCLUSION: This is a zero-background strategy for direct cloning into a binary vector, optimized for the generation of infectious viroids. As a result, this methodology constitutes a powerful tool for viroid research and exemplifies the applicability of type IIs restriction enzymes and the lethal gene ccdB to design efficient and affordable direct cloning approaches of PCR products into binary vectors.

Identification and Characterization of Stress-Responsive TAS3-Derived TasiRNAs in Melon.

Small interfering RNAs (siRNA) are key regulators of gene expression that play essential roles in diverse biological processes. Trans-acting siRNAs (tasiRNAs) are a class of plant-endogenous siRNAs that lead the cleavage of nonidentical transcripts. TasiRNAs are usually involved in fine-tuning development. However, increasing evidence supports that tasiRNAs may be involved in stress response. Melon is a crop of great economic importance extensively cultivated in semiarid regions frequently exposed to changing environmental conditions that limit its productivity. However, knowledge of the precise role of siRNAs in general, and of tasiRNAs in particular, in regulating the response to adverse environmental conditions is limited. Here, we provide the first comprehensive analysis of computationally inferred melon-tasiRNAs responsive to two biotic (viroid-infection) and abiotic (cold treatment) stress conditions. We identify two TAS3-loci encoding to length (TAS3-L) and short (TAS3-S) transcripts. The TAS candidates predicted from small RNA-sequencing data were characterized according to their chromosome localization and expression pattern in response to stress. The functional activity of cmTAS genes was validated by transcript quantification and degradome assays of the tasiRNA precursors and their predicted targets. Finally, the functionality of a representative cmTAS3-derived tasiRNA (TAS3-S) was confirmed by transient assays showing the cleavage of ARF target transcripts.