Virus-derived siRNAs and piRNAs in immunity and pathogenesis

VmRDR2 of Valsa mali mediates the generation of VmR2-siR1 that suppresses apple resistance by RNA interference

Cross-kingdom RNA interference (RNAi) is a crucial mechanism in host-pathogen interactions, with RNA-dependent RNA polymerase (RdRP) playing a vital role in signal amplification during RNAi. However, the role of pathogenic fungal RdRP in siRNAs generation and the regulation of plant-pathogen interactions remains elusive. Using deep sequencing, molecular, genetic, and biochemical approaches, this study revealed that VmRDR2 of Valsa mali regulates VmR2-siR1 to suppress the disease resistance-related gene MdLRP14 in apple. Both VmRDR1 and VmRDR2 are essential for the pathogenicity of V. mali in apple, with VmRDR2 mediating the generation of endogenous siRNAs, including an infection-related siRNA, VmR2-siR1. This siRNA specifically degrades the apple intracellular LRR-RI protein gene MdLRP14 in a sequence-specific manner, and overexpression of MdLRP14 enhances apple resistance against V. mali, which can be suppressed by VmR2-siR1. Conversely, MdLRP14 knockdown reduces resistance. In summary, this study demonstrates that VmRDR2 contributes to the generation of VmR2-siR1, which silences the host's intracellular LRR protein gene, thereby inhibiting host resistance. These findings offer novel insights into the fungi-mediated pathogenicity mechanism through RNAi.

Differential expression of microRNAs in response to Papaya ringspot virus infection in differentially responding genotypes of papaya (Carica papaya L.) and its wild relative

Papaya ringspot virus (PRSV) is one of the most devastating viruses of papaya that has significantly hampered papaya production across the globe. Although PRSV resistance is known in some of its wild relatives, such as Vasconcellea cauliflora and in some of the improved papaya genotypes, the molecular basis of this resistance mechanism has not been studied and understood. Plant microRNAs are an important class of small RNAs that regulate the gene expression in several plant species against the invading plant pathogens. These miRNAs are known to manifest the expression of genes involved in resistance against plant pathogens, through modulation of the plant's biochemistry and physiology. In this study we made an attempt to study the overall expression pattern of small RNAs and more specifically the miRNAs in different papaya genotypes from India, that exhibit varying levels of tolerance or resistance to PRSV. Our study found that the expression of some of the miRNAs was differentially regulated in these papaya genotypes and they had entirely different miRNA expression profile in healthy and PRSV infected symptomatic plants. This data may help in improvement of papaya cultivars for resistance against PRSV through new breeding initiatives or biotechnological approaches such as genome editing.

RNA virus diversity in three parasitoid wasps of tephritid flies: insights from novel and known species

IMPORTANCE

Parasitoid wasp populations have developed persistent beneficial symbiotic relationships with several viruses through repeated evolution. However, there have been limited reports on RNA viruses in parasitoid wasps of tephritid flies, a significant pest group affecting fruits and vegetables. This study explores the diversity of RNA viruses in three parasitoid wasps of tephritid flies and highlights the potential biological significance of specific viruses in Diachasmimorpha longicaudata. These findings have important implications for the development of sustainable pest management strategies and the enhancement of artificial rearing techniques for parasitoid wasps.

Infection Defects of RNA and DNA Viruses Induced by Antiviral RNA Interference

Immune recognition of viral genome-derived double-stranded RNA (dsRNA) molecules and their subsequent processing into small interfering RNAs (siRNAs) in plants, invertebrates, and mammals trigger specific antiviral immunity known as antiviral RNA interference (RNAi). Immune sensing of viral dsRNA is sequence-independent, and most regions of viral RNAs are targeted by virus-derived siRNAs which extensively overlap in sequence. Thus, the high mutation rates of viruses do not drive immune escape from antiviral RNAi, in contrast to other mechanisms involving specific virus recognition by host immune proteins such as antibodies and resistance (R) proteins in mammals and plants, respectively. Instead, viruses actively suppress antiviral RNAi at various key steps with a group of proteins known as viral suppressors of RNAi (VSRs). Some VSRs are so effective in virus counter-defense that potent inhibition of virus infection by antiviral RNAi is undetectable unless the cognate VSR is rendered nonexpressing or nonfunctional. Since viral proteins are often multifunctional, resistance phenotypes of antiviral RNAi are accurately defined by those infection defects of VSR-deletion mutant viruses that are efficiently rescued by host deficiency in antiviral RNAi. Here, we review and discuss in vivo infection defects of VSR-deficient RNA and DNA viruses resulting from the actions of host antiviral RNAi in model systems.

Translation Arrest: A Key Player in Plant Antiviral Response

Plants evolved several mechanisms to protect themselves against viruses. Besides recessive resistance, where compatible host factors required for viral proliferation are absent or incompatible, there are (at least) two types of inducible antiviral immunity: RNA silencing (RNAi) and immune responses mounted upon activation of nucleotide-binding domain leucine-rich repeat (NLR) receptors. RNAi is associated with viral symptom recovery through translational repression and transcript degradation following recognition of viral double-stranded RNA produced during infection. NLR-mediated immunity is induced upon (in)direct recognition of a viral protein by an NLR receptor, triggering either a hypersensitive response (HR) or an extreme resistance response (ER). During ER, host cell death is not apparent, and it has been proposed that this resistance is mediated by a translational arrest (TA) of viral transcripts. Recent research indicates that translational repression plays a crucial role in plant antiviral resistance. This paper reviews current knowledge on viral translational repression during viral recovery and NLR-mediated immunity. Our findings are summarized in a model detailing the pathways and processes leading to translational arrest of plant viruses. This model can serve as a framework to formulate hypotheses on how TA halts viral replication, inspiring new leads for the development of antiviral resistance in crops.

Genes Involved in miRNA Biogenesis Are Not Downregulated in SARS-CoV-2 Infection

miRNAs, small non-coding RNAs that regulate gene expression, are involved in various pathological processes, including viral infections. Virus infections may interfere with the miRNA pathway through the inhibition of genes involved in miRNA biogenesis. A reduction in the number and the levels of miRNAs expressed in nasopharyngeal swabs of patients with severe COVID-19 was lately observed by us, pointing towards the potential of miRNAs as possible diagnostic or prognostic biomarkers for predicting outcomes among patients with severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection. The objective of the present study was to investigate whether SARS-CoV-2 infection influences the expression levels of messenger RNAs (mRNAs) of key genes involved in miRNA biogenesis. mRNA levels of AGO2, DICER1, DGCR8, DROSHA, and Exportin-5 (XPO5) were measured by quantitative reverse-transcription polymerase chain reaction (RT-qPCR) in nasopharyngeal swab specimens from patients with COVID-19 and controls, as well as in cells infected with SARS-CoV-2 in vitro. Our data showed that the mRNA expression levels of AGO2, DICER1, DGCR8, DROSHA, and XPO5 were not significantly different in patients with severe COVID-19 when compared to patients with non-severe COVID-19 and controls. Similarly, the mRNA expression of these genes was not affected by SARS-CoV-2 infection in NHBE and Calu-3 cells. However, in Vero E6 cells, AGO2, DICER1, DGCR8, and XPO5 mRNA levels were slightly upregulated 24 h after infection with SARS-CoV-2. In conclusion, we did not find evidence for downregulation of mRNA levels of miRNA biogenesis genes during SARS-CoV-2 infection, neither ex vivo nor in vitro.

An efficient deep learning based predictor for identifying miRNA-triggered phasiRNA loci in plant

Phasic small interfering RNAs are plant secondary small interference RNAs that typically generated by the convergence of miRNAs and polyadenylated mRNAs. A growing number of studies have shown that miRNA-initiated phasiRNA plays crucial roles in regulating plant growth and stress responses. Experimental verification of miRNA-initiated phasiRNA loci may take considerable time, energy and labor. Therefore, computational methods capable of processing high throughput data have been proposed one by one. In this work, we proposed a predictor (DIGITAL) for identifying miRNA-initiated phasiRNAs in plant, which combined a multi-scale residual network with a bi-directional long-short term memory network. The negative dataset was constructed based on positive data, through replacing 60% of nucleotides randomly in each positive sample. Our predictor achieved the accuracy of 98.48% and 94.02% respectively on two independent test datasets with different sequence length. These independent testing results indicate the effectiveness of our model. Furthermore, DIGITAL is of robustness and generalization ability, and thus can be easily extended and applied for miRNA target recognition of other species. We provide the source code of DIGITAL, which is freely available at https://github.com/yuanyuanbu/DIGITAL.

Micro"bee"ota: Honey Bee Normal Microbiota as a Part of Superorganism

Honey bees are model organisms for microbiota research. Gut microbiomes are very interesting for surveys due to their simple structure and relationship with hive production. Long-term studies reveal the gut microbiota patterns of various hive members, as well as the functions, sources, and interactions of the majority of its bacteria. But the fungal non-pathogenic part of gut microbiota is almost unexplored, likewise some other related microbiota. Honey bees, as superorganisms, interact with their own microorganisms, the microbial communities of food stores, hive surfaces, and other environments. Understanding microbiota diversity, its transition ways, and hive niche colonization control are necessary for understanding any separate microbiota niche because of their interplay. The long coevolution of bees with the microorganisms populating these niches makes these systems co-dependent, integrated, and stable. Interaction with the environment, hive, and other bees determines caste lifestyle as well as individual microbiota. In this article, we bring together studies on the microbiota of the western honey bee. We show a possible relationship between caste determination and microbiota composition. And what is primary: caste differentiation or microbiota composition?

The Diversity of Viral Community in Invasive Fruit Flies (Bactrocera and Zeugodacus) Revealed by Meta-transcriptomics

RNA viruses are extremely diverse and rapidly evolving in various organisms. Our knowledge on viral evolution with interacted hosts in the manner of ecology is still limited. In the agricultural ecosystem, invasive insect species are posing a great threat to sustainable crop production. Among them, fruit flies (Diptera: Tephritidae Bactrocera and Zeugodacus) are destructive to fruits and vegetables, which are also closely related and often share similar ecological niches. Thus, they are ideal models for investigating RNA virome dynamics in host species. Using meta-transcriptomics, we found 39 viral sequences in samples from 12 fly species. These viral species represented the diversity of the viromes including Dicistroviridae, negev-like virus clades, Thika virus clades, Solemoviridae, Narnaviridae, Nodaviridae, Iflaviridae, Orthomyxoviridae, Bunyavirales, Partitiviridae, and Reoviridae. In particular, dicistrovirus, negev-like virus, orthomyxovirus, and orbivirus were common in over four of the fly species, which suggests a positive interaction between fly viromes that exist under the same ecological conditions. For most of the viruses, the virus-derived small RNAs displayed significantly high peaks in 21 nt and were symmetrically distributed throughout the viral genome. These results suggest that infection by these viruses can activate the host's RNAi immunity. Our study provides RNA virome diversity and evidence on their infection activity in ecologically associated invasive fruit fly species, which could help our understanding of interactions between complex species and viruses.

Virus and endogenous viral element-derived small non-coding RNAs and their roles in insect-virus interaction

RNA interference pathways mediated by different types of small non-coding RNAs (siRNAs, miRNAs and piRNAs) are conserved biological responses to exotic stresses, including viral infection. Aside from the well-established siRNA pathway, the miRNA pathway and the piRNA pathway process viral sequences, exogenously or endogenously, into miRNAs and piRNAs, respectively. During the host-virus interaction, viral sequences, including both coding and non-coding sequences, can be integrated as endogenous viral elements (EVEs) and thereby become present within the germline of a non-viral organism. In recent years, significant progress has been made in characterizing the biogenesis and function of viruses and EVEs associated with snRNAs. Overall, the siRNA pathway acts as the primarily antiviral defense against a wide range of exogenous viruses; the miRNA pathways associated with viruses or EVEs function in antiviral response and host gene regulation; EVE derived piRNAs with a ping-pong signature have the potential to limit cognate viral infection.

Next-generation sequencing identification and multiplex RT-PCR detection for viruses infecting cigar and flue-cured tobacco

BACKGROUND

Early, precise and simultaneous identification of plant viruses is of great significance for preventing virus spread and reducing losses in agricultural yields.

METHODS AND RESULTS

In this study, the identification of plant viruses from symptomatic samples collected from a cigar tobacco planting area in Deyang and a flue-cured tobacco planting area in Luzhou city, Sichuan Province, China, was conducted by deep sequencing of small RNAs (sRNAs) through an Illumina sequencing platform, and plant virus-specific contigs were generated based on virus-derived siRNA sequences. Additionally, sequence alignment and phylogenetic analysis were performed to determine the species or strains of these viruses. A total of 27930450, 21537662 and 28194021 clean reads were generated from three pooled samples, with a total of 105 contigs mapped to the closest plant viruses with lengths ranging from 34 ~ 1720 nt. The results indicated that the major viruses were potato virus Y, Chilli veinal mottle virus, tobacco vein banding mosaic virus, tobacco mosaic virus and cucumber mosaic virus. Subsequently, a fast and sensitive multiplex reverse transcription polymerase chain reaction assay was developed for the simultaneous detection of the most frequent RNA viruses infecting cigar and flue-cured tobacco in Sichuan.

CONCLUSIONS

These results provide a theoretical basis and convenient methods for the rapid detection and control of viruses in cigar- and flue-cured tobacco.

Use of deep sequencing to profile small RNAs derived from tomato spotted wilt orthotospovirus and hippeastrum chlorotic ringspot orthotospovirus in infected Capsicum annuum

Virus-derived small RNAs are one of the key factors of RNA silencing in plant defence against viruses. We obtained virus-derived small interfering RNA profiles from Tomato spotted wilt orthotospovirus and Hippeastrum chlorotic ringspot orthotospovirus infected Capsicum annuum XX19 and XY11 by deep sequencing one day after inoculation. The vsiRNAs data were mapped to the TSWV and HCRV genomes, and the results showed that the vsiRNAs measured 19-24 nucleotides in length. Most of the vsiRNAs were mapped to the S segment of the viral genome. For XX19 and XY11 infected with HCRV, the distribution range of vsiRNAs in S RNA was 52.06-55.20%, while for XX19 and XY11 infected with TSWV, the distribution range of vsiRNAs in S RNA was 87.76-89.07%. The first base at the 5' end of the siRNA from TSWV and HCRV was primarily biased towards A, U, or C. Compared with mock-inoculated XX19 and XY11, the expression level of CaRDR1 was upregulated in TSWV- and HCRV-inoculated XX19 and XY11. CaAGO2 and CaAGO5 were upregulated in XY11 against HCRV infection, and CaRDR2 was downregulated in TSWV-infected XY11 and XX19. The profile of HCRV and TSWV vsiRNA verified in this study could be useful for selecting key vsiRNA such as those in disease-resistant varieties by artificially synthesizing amiRNA.

Non-Coding RNAs in Response to Drought Stress

Drought stress causes changes in the morphological, physiological, biochemical and molecular characteristics of plants. The response to drought in different plants may vary from avoidance, tolerance and escape to recovery from stress. This response is genetically programmed and regulated in a very complex yet synchronized manner. The crucial genetic regulations mediated by non-coding RNAs (ncRNAs) have emerged as game-changers in modulating the plant responses to drought and other abiotic stresses. The ncRNAs interact with their targets to form potentially subtle regulatory networks that control multiple genes to determine the overall response of plants. Many long and small drought-responsive ncRNAs have been identified and characterized in different plant varieties. The miRNA-based research is better documented, while lncRNA and transposon-derived RNAs are relatively new, and their cellular role is beginning to be understood. In this review, we have compiled the information on the categorization of non-coding RNAs based on their biogenesis and function. We also discuss the available literature on the role of long and small non-coding RNAs in mitigating drought stress in plants.

Ferredoxin 1 is downregulated by the accumulation of abscisic acid in an ABI5-dependent manner to facilitate rice stripe virus infection in Nicotiana benthamiana and rice

Ferredoxin 1 (FD1) accepts and distributes electrons in the electron transfer chain of plants. Its expression is universally downregulated by viruses and its roles in plant immunity have been brought into focus over the past decade. However, the mechanism by which viruses regulate FD1 remains to be defined. In a previous report, we found that the expression of Nicotiana benthamiana FD1 (NbFD1) was downregulated following infection with potato virus X (PVX) and that NbFD1 regulates callose deposition at plasmodesmata to play a role in defense against PVX infection. We now report that NbFD1 is downregulated by rice stripe virus (RSV) infection and that silencing of NbFD1 also facilitates RSV infection, while viral infection was inhibited in a transgenic line overexpressing NbFD1, indicating that NbFD1 also functions in defense against RSV infection. Next, a RSV-derived small interfering RNA was identified that contributes to the downregulation of FD1 transcripts. Further analysis showed that the abscisic acid (ABA) which accumulates in RSV-infected plants also represses NbFD1 transcription. It does this by stimulating expression of ABA insensitive 5 (ABI5), which binds the ABA response element motifs in the NbFD1 promoter, resulting in negative regulation. Regulation of FD1 by ABA was also confirmed in RSV-infected plants of the natural host rice. The results therefore suggest a mechanism by which virus regulates chloroplast-related genes to suppress their defense roles.

Circular RNA as an Additional Player in the Conflicts Between the Host and the Virus

Circular RNA (circRNA), a relatively new member of the non-coding RNA family, has spurred great interest among researchers following its discovery as a ubiquitous class within the RNA world. Rapid progress in circRNA biology has coincided with its identification in a plethora of diverse roles including regulation of gene expression and probable coding potential, as well as competing interactions with proteins and microRNAs in various pathological conditions. Emerging evidence suggests that circRNAs also function in viral infections. The deregulation of circRNAs during viral infection has prompted investigations into the possibilities of circRNA as a competing endogenous RNA (ceRNA) that modulates response to infection. Recently, viruses have been shown to encode circRNAs with proviral functions, providing a strong impetus for focused efforts to elucidate the networks coaxed by circRNAs during infection. This review elaborates on recent insights gained on the roles of circRNAs during virus infection and immunity.

High Throughput Sequencing-Aided Survey Reveals Widespread Mixed Infections of Whitefly-Transmitted Viruses in Cucurbits in Georgia, USA

Viruses transmitted by the sweet potato whitefly (Bemisia tabaci) have been detrimental to the sustainable production of cucurbits in the southeastern USA. Surveys were conducted in the fall of 2019 and 2020 in Georgia, a major cucurbit-producing state of the USA, to identify the viruses infecting cucurbits and their distribution. Symptomatic samples were collected and small RNA libraries were prepared and sequenced from three cantaloupes, four cucumbers, and two yellow squash samples. An analysis of the sequences revealed the presence of the criniviruses cucurbit chlorotic yellows virus (CCYV), cucurbit yellow stunting disorder virus (CYSDV), and the begomovirus cucurbit leaf crumple virus (CuLCrV). CuLCrV was detected in 76%, CCYV in 60%, and CYSDV in 43% of the total samples (n = 820) tested. The level of mixed infections was high in all the cucurbits, with most plants tested being infected with at least two of these viruses. Near-complete genome sequences of two criniviruses, CCYV and CYSDV, were assembled from the small RNA sequences. An analysis of the coding regions showed low genetic variability among isolates from different hosts. In phylogenetic analysis, the CCYV isolates from Georgia clustered with Asian isolates, while CYSDV isolates clustered with European and USA isolates. This work enhances our understanding of the distribution of viruses on cucurbits in South Georgia and will be useful to develop strategies for managing the complex of whitefly-transmitted viruses in the region.

From genomes to molecular dynamics - A bottom up approach in extrication of SARS CoV-2 main protease inhibitors

The recent pandemic Coronavirus disease-19 outbreak had traumatized global countries since its origin in late December 2019. Though the virus originated in China, it has spread rapidly across the world due its firmly established community transmission. To successfully tackle the spread and further infection, there needs a clear multidimensional understanding of the molecular mechanisms. Henceforth, 942 viral genome sequences were analysed to predict the core genomes crucial in virus life cycle. Additionally, 35 small interfering RNA transcripts were predicted that can target specifically the viral core proteins and reduce pathogenesis. The crystal structure of Covid-19 main protease-6LU7 was chosen as an attractive target due to the factors that there were fewer mutations and whose structure had significant identity to the annotated protein sequence of the core genome. Drug repurposing of both recruiting and non recruiting drugs was carried out through molecular docking procedures to recognize bitolterol as a good inhibitor of Covid-19 protease. The study was extended further to screen antiviral phytocompounds through quantitative structure activity relationship and molecular docking to identify davidigenin, from licorice as the best novel lead with good interactions and binding energy. The docking of the best compounds in all three categories was validated with molecular dynamics simulations which implied stable binding of the drug and lead molecule. Though the studies need clinical evaluations, the results are suggestive of curbing the pandemic.

Characterization of Pathogenicity-Associated V2 Protein of Tobacco Curly Shoot Virus

V2 proteins encoded by some whitefly-transmitted geminiviruses were reported to be functionally important proteins. However, the functions of the V2 protein of tobacco curly shoot virus (TbCSV), a monopartite begomovirus that causes leaf curl disease on tomato and tobacco in China, remains to be characterized. In our report, an Agrobacterium infiltration-mediated transient expression assay indicated that TbCSV V2 can suppress local and systemic RNA silencing and the deletion analyses demonstrated that the amino acid region 1-92 of V2, including the five predicted α-helices, are required for local RNA silencing suppression. Site-directed substitutions showed that the conserved basic and ring-structured amino acids in TbCSV V2 are critical for its suppressor activity. Potato virus X-mediated heteroexpression of TbCSV V2 in Nicotiana benthamiana induced hypersensitive response-like (HR-like) cell death and systemic necrosis in a manner independent of V2's suppressor activity. Furthermore, TbCSV infectious clone mutant with untranslated V2 protein (TbCSV∆V2) could not induce visual symptoms, and coinfection with betasatellite (TbCSB) could obviously elevate the viral accumulation and symptom development. Interestingly, symptom recovery occurred at 15 days postinoculation (dpi) and onward in TbCSV∆V2/TbCSB-inoculated plants. The presented work contributes to understanding the RNA silencing suppression activity of TbCSV V2 and extends our knowledge of the multifunctional role of begomovirus-encoded V2 proteins during viral infections.

Viromes of Ten Alfalfa Plants in Australia Reveal Diverse Known Viruses and a Novel RNA Virus

Alfalfa plants in the field can display a range of virus-like symptoms, especially when grown over many years for seed production. Most known alfalfa viruses have RNA genomes, some of which can be detected using diagnostic assays, but many viruses of alfalfa are not well characterized. This study aims to identify the RNA and DNA virus complexes associated with alfalfa plants in Australia. To maximize the detection of RNA viruses, we purified double-stranded RNA (dsRNA) for high throughput sequencing and characterized the viromes of ten alfalfa samples that showed diverse virus-like symptoms. Using Illumina sequencing of tagged cDNA libraries from immune-captured dsRNA, we identified sequences of the single-stranded RNA viruses, alfalfa mosaic virus (AMV), bean leafroll virus, a new emaravirus tentatively named alfalfa ringspot-associated virus, and persistent dsRNA viruses belonging to the families Amalgaviridae and Partitiviridae. Furthermore, rolling circle amplification and restriction enzyme digestion revealed the complete genome of chickpea chlorosis Australia virus, a mastrevirus (family Geminiviridae) previously reported only from chickpea and French bean that was 97% identical to the chickpea isolate. The sequence data also enabled the assembly of the first complete genome (RNAs 1–3) of an Australian AMV isolate from alfalfa.

Identification of microRNAs regulated by tobacco curly shoot virus co-infection with its betasatellite in Nicotiana benthamiana

BACKGROUND

MicroRNAs (miRNAs) are a class of 21-24 nucleotide endogenous non-coding small RNAs that play important roles in plant development and defense responses to biotic and abiotic stresses. Tobacco curly shoot virus (TbCSV) is a monopartite begomovirus, cause leaf curling and plant stunting symptoms in many Solanaceae plants. The betasatellite of TbCSV (TbCSB) induces more severe symptoms and enhances virus accumulation when co-infect the plants with TbCSV.

METHODS

In this study, miRNAs regulated by TbCSV and TbCSB co-infection in Nicotiana benthamiana were characterized using high-throughput sequencing technology.

RESULTS

Small RNA sequencing analysis revealed that a total of 13 known miRNAs and 42 novel miRNAs were differentially expressed in TbCSV and TbCSB co-infected N. benthamiana plants. Several potential miRNA-targeted genes were identified through data mining and were involved in both catalytic and metabolic processes, in addition to plant defense mechanisms against virus infections according to Gene Ontology (GO) analyses. In addition, the expressions of several differentially expressed miRNAs and their miRNA-targeted gene were validated through quantitative real time polymerase chain reaction (qRT-PCR) approach.

CONCLUSIONS

A large number of miRNAs are identified, and their target genes, functional annotations also have been explored. Our results provide the information on N. benthamiana miRNAs and would be useful to further understand miRNA regulatory mechanisms after TbCSV and TbCSB co-infection.

Widespread selection for extremely high and low levels of secondary structure in coding sequences across all domains of life

Codon composition, GC content and local RNA secondary structures can have a profound effect on gene expression, and mutations affecting these parameters, even though they do not alter the protein sequence, are not neutral in terms of selection. Although evidence exists that, in some cases, selection favours more stable RNA secondary structures, we currently lack a concrete idea of how many genes are affected within a species, and whether this is a universal phenomenon in nature. We searched for signs of structural selection in a global manner, analysing a set of 1 million coding sequences from 73 species representing all domains of life, as well as viruses, by means of our newly developed software PACKEIS. We show that codon composition and amino acid identity are main determinants of RNA secondary structure. In addition, we show that the arrangement of synonymous codons within coding sequences is non-random, yielding extremely high, but also extremely low, RNA structuredness significantly more often than expected by chance. Taken together, we demonstrate that selection for high and low levels of secondary structure is a widespread phenomenon. Our results provide another line of evidence that synonymous mutations are less neutral than commonly thought, which is of importance for many evolutionary models.

Tomato yellow leaf curl virus intergenic siRNAs target a host long noncoding RNA to modulate disease symptoms

Tomato yellow leaf curl virus (TYLCV) and its related begomoviruses cause fast-spreading diseases in tomato worldwide. How this virus induces diseases remains largely unclear. Here we report a noncoding RNA-mediated model to elucidate the molecular mechanisms of TYLCV-tomato interaction and disease development. The circular ssDNA genome of TYLCV contains a noncoding intergenic region (IR), which is known to mediate viral DNA replication and transcription in host cells, but has not been reported to contribute directly to viral disease development. We demonstrate that the IR is transcribed in dual orientations during plant infection and confers abnormal phenotypes in tomato independently of protein-coding regions of the viral genome. We show that the IR sequence has a 25-nt segment that is almost perfectly complementary to a long noncoding RNA (lncRNA, designated as SlLNR1) in TYLCV-susceptible tomato cultivars but not in resistant cultivars which contains a 14-nt deletion in the 25-nt region. Consequently, we show that viral small-interfering RNAs (vsRNAs) derived from the 25-nt IR sequence induces silencing of SlLNR1 in susceptible tomato plants but not resistant plants, and this SlLNR1 downregulation is associated with stunted and curled leaf phenotypes reminiscent of TYLCV symptoms. These results suggest that the lncRNA interacts with the IR-derived vsRNAs to control disease development during TYLCV infection. Consistent with its possible function in virus disease development, over-expression of SlLNR1 in tomato reduces the accumulation of TYLCV. Furthermore, gene silencing of the SlLNR1 in the tomato plants induced TYLCV-like leaf phenotypes without viral infection. Our results uncover a previously unknown interaction between vsRNAs and host lncRNA, and provide a plausible model for TYLCV-induced diseases and host antiviral immunity, which would help to develop effective strategies for the control of this important viral pathogen.

Impact of Two Reoviruses and Their Coinfection on the Rice RNAi System and vsiRNA Production

Both Southern rice black-streaked dwarf virus (SRBSDV) and Rice ragged stunt virus (RRSV) belong to the family Reoviridae, and synergistic infection of these two viruses commonly occurs in the field. This study revealed that both SRBSDV and RRSV affect the RNA interference (RNAi) pathway and form different virus-derived interfering RNA (vsiRNA) profiles in rice. Co-infection of rice by SRBSDV and RRSV up-regulated the expression of rice DICER-like (DCL) proteins but down-regulated the expression of rice RNA-dependent RNA polymerases (RDRs), and the accumulation of vsiRNAs of either RBSDV or RRSV was decreased compared with that in singly infected plants. The majority of SRBSDV vsiRNAs were 21 nt or 22 nt in length, whether plants were singly infected with SRBSDV or co-infected with RRSV. On the other hand, the majority of RRSV vsiRNAs were 20 nt, 21 nt, or 22 nt in length, among which those 20 nt in length accounted for the largest proportion; co-infection with SRBSDV further increased the proportion of 20 nt vsiRNAs and decreased the proportion of 21 nt vsiRNAs. Co-infection had no effects on the strand favoritism and hot spots of the vsiRNAs, but changed the bias of the 5′ terminal nucleotide significantly. This study provides a reference for further study on the pathogenesis and synergistic mechanism of SRBSDV and RRSV.

Transgene-Assisted Genetic Screen Identifies rsd-6 and Novel Genes as Key Components of Antiviral RNA Interference in Caenorhabditis elegans

RNA interference (RNAi) is a widespread antiviral mechanism triggered by virus-produced double-stranded RNAs (dsRNAs). In Caenorhabditis elegans, antiviral RNAi involves a RIG-I-like RNA helicase, termed DRH-1 (dicer related RNA helicase 1), that is not required for classical RNAi triggered by artificial dsRNA. Currently, whether antiviral RNAi in C. elegans involves novel factors that are dispensable for classical RNAi remains an open question. To address this question, we designed and carried out a genetic screen that aims to identify novel genes involved in worm antiviral RNAi. By introducing extra copies of known antiviral RNAi genes into the reporter worms, we managed to reject alleles derived from 4 known antiviral RNAi genes, including the DRH-1 coding gene, during the screen. Our genetic screen altogether identified 25 alleles, which were assigned to 11 candidate genes and 2 known antiviral RNAi genes through genetic complementation tests. Using a mapping-by-sequencing strategy, we identified one of the candidate genes as rsd-6, a gene that helps maintain genome integrity through an endogenous gene-silencing pathway but was not known to be required for antiviral RNAi. More importantly, we found that two of the candidate genes are required for antiviral RNAi targeting Orsay virus, a natural viral pathogen of C. elegans, but dispensable for classical RNAi. Since drh-1 is so far the only antiviral RNAi gene not required for classical RNAi, we believe that our genetic screen led to identification of novel worm genes that may target virus-specific features to function in RNAi.IMPORTANCE In nematode worms, drh-1 detects virus-produced double-stranded RNA (dsRNA), thereby specifically contributing to antiviral RNA silencing. To identify drh-1-like genes with dedicated function in antiviral RNAi, we recently carried out a genetic screen that was designed to automatically reject all alleles derived from 4 known antiviral silencing genes, including drh-1 Of the 11 candidate genes identified, we found two of them to be required for antiviral silencing targeting a natural viral pathogen of C. elegans but not for classical RNA silencing triggered by artificial dsRNA. We believe that these two genes are novel components of worm antiviral RNAi, considering the fact that drh-1 is the only known antiviral RNAi gene that is dispensable for classical RNAi. This genetic screen also identified rsd-6, a gene that maintains genome integrity under unfavorable conditions, as a key regulator of worm antiviral silencing, demonstrating an interplay between antiviral immunity and genome integrity maintenance.

MicroRNA and cellular targets profiling reveal miR-217 and miR-576-3p as proviral factors during Oropouche infection

Oropouche Virus is the etiological agent of an arbovirus febrile disease that affects thousands of people and is widespread throughout Central and South American countries. Although isolated in 1950’s, still there is scarce information regarding the virus biology and its prevalence is likely underestimated. In order to identify and elucidate interactions with host cells factors and increase the understanding about the Oropouche Virus biology, we performed microRNA (miRNA) and target genes screening in human hepatocarcinoma cell line HuH-7. Cellular miRNAs are short non-coding RNAs that regulates gene expression post-transcriptionally and play key roles in several steps of viral infections. The large scale RT-qPCR based screening found 13 differentially expressed miRNAs in Oropouche infected cells. Further validation confirmed that miR-217 and miR-576-3p were 5.5 fold up-regulated at early stages of virus infection (6 hours post-infection). Using bioinformatics and pathway enrichment analysis, we predicted the cellular targets genes for miR-217 and miR-576-3p. Differential expression analysis of RNA from 95 selected targets revealed genes involved in innate immunity modulation, viral release and neurological disorder outcomes. Further analysis revealed the gene of decapping protein 2 (DCP2), a previous known restriction factor for bunyaviruses transcription, as a miR-217 candidate target that is progressively down-regulated during Oropouche infection. Our analysis also showed that activators genes involved in innate immune response through IFN-β pathway, as STING (Stimulator of Interferon Genes) and TRAF3 (TNF-Receptor Associated Factor 3), were down-regulated as the infection progress. Inhibition of miR-217 or miR-576-3p restricts OROV replication, decreasing viral RNA (up to 8.3 fold) and virus titer (3 fold). Finally, we showed that virus escape IFN-β mediated immune response increasing the levels of cellular miR-576-3p resulting in a decreasing of its partners STING and TRAF3. We concluded stating that the present study, the first for a Peribunyaviridae member, gives insights in its prospective pathways that could help to understand virus biology, interactions with host cells and pathogenesis, suggesting that the virus escapes the antiviral cellular pathways increasing the expression of cognates miRNAs.

Oropouche Virus causes typical arboviral febrile illness and is widely distributed in tropical region of Americas, mainly Amazon region, associated with cases of encephalitis. 500,000 people are estimated to be infected with Oropouche worldwide and some states in Brazil detected higher number of cases among other arboviruses such as Dengue and Chikungunya. As much as climate change, human migration and vector and host availability might increase the risk of virus transmission. Despite its estimated high prevalence in Central and South America populations, the literature concerning the main aspects of viral biology remain scarce and began to be investigated only in the last two decades. Nonetheless, little is known about virus-host cell interactions and pathogenesis. Virus infection regulates cellular pathways either promoting its replication or escaping from immune response through microRNAs. Knowing which microRNAs and target genes are modulated in infection could give us new insights to understand multiple aspects of infection. Here, we depicted candidate miRNAs, genes and pathways affected by Oropouche Virus infection in hepatocyte cells. We hope this work serve as guideline for prospective studies in order to assess the complexity regarding the orthobunyaviruses infections.

The protective function of noncoding DNA in genome defense of eukaryotic male germ cells

Peripheral and abundant noncoding DNA has been hypothesized to protect the genome and the central protein-coding sequences against DNA damage in somatic genome. In the cytosol, invading exogenous nucleic acids may first be deactivated by small RNAs encoded by noncoding DNA via mechanisms similar to the prokaryotic CRISPR-Cas system. In the nucleus, the radicals generated by radiation in the cytosol, radiation energy and invading exogenous nucleic acids are absorbed, blocked and/or reduced by peripheral heterochromatin, and damaged DNA in heterochromatin is removed and excluded from the nucleus to the cytoplasm through nuclear pore complexes. To further strengthen the hypothesis, this review summarizes the experimental evidence supporting the protective function of noncoding DNA in the genome of male germ cells. Based on these data, this review provides evidence supporting the protective role of noncoding DNA in the genome defense of sperm genome through similar mechanisms to those of the somatic genome.

Differential miRNA expression in B cells is associated with inter-individual differences in humoral immune response to measles vaccination

Background

MicroRNAs are important mediators of post-transcriptional regulation of gene expression through RNA degradation and translational repression, and are emerging biomarkers of immune system activation/response after vaccination.

Methods

We performed Next Generation Sequencing (mRNA-Seq) of intracellular miRNAs in measles virus-stimulated B and CD4⁺ T cells from high and low antibody responders to measles vaccine. Negative binomial generalized estimating equation (GEE) models were used for miRNA assessment and the DIANA tool was used for gene/target prediction and pathway enrichment analysis.

Results

We identified a set of B cell-specific miRNAs (e.g., miR-151a-5p, miR-223, miR-29, miR-15a-5p, miR-199a-3p, miR-103a, and miR-15a/16 cluster) and biological processes/pathways, including regulation of adherens junction proteins, Fc-receptor signaling pathway, phosphatidylinositol-mediated signaling pathway, growth factor signaling pathway/pathways, transcriptional regulation, apoptosis and virus-related processes, significantly associated with neutralizing antibody titers after measles vaccination. No CD4⁺ T cell-specific miRNA expression differences between high and low antibody responders were found.

Conclusion

Our study demonstrates that miRNA expression directly or indirectly influences humoral immunity to measles vaccination and suggests that B cell-specific miRNAs may serve as useful predictive biomarkers of vaccine humoral immune response.

Dicer-like Proteins Regulate the Growth, Conidiation, and Pathogenicity of Colletotrichum gloeosporioides from Hevea brasiliensis

Colletotrichum gloeosporioides from Hevea brasiliensis is the hemibiotrophic fungi which could cause anthracnose in rubber trees. Dicer like proteins (DCL) were the core enzymes for generation of small RNAs. In the present study, the knocking-out mutants of two dicer like proteins encoding genes of C. gloeosporioides were constructed; and functions of two proteins were investigated. The results showed that DCL play important roles in regulating the growth, conidiation and pathogenicity of C. gloeosporioides; and there is a functional redundancy between DCL1 and DCL2. Microscopy analysis and DAB staining revealed that loss of penetration ability into the host cells, instead of the decreased growth rate, was the main cause for the impaired pathogenicity of the ΔDcl1ΔDcl2 double mutant. Proteomics analysis suggested that DCL proteins affected the expression of functional proteins to regulating multiple biological processes of C. gloeosporioides. These data lead to a better understanding of the functions of DCL proteins in regulating the development and pathogenesis of C. gloeosporioides.

PVsiRNAdb: a database for plant exclusive virus-derived small interfering RNAs

Ribonucleic acids (RNA) interference mechanism has been proved to be an important regulator of both transcriptional and post-transcription controls of gene expression during biotic and abiotic stresses in plants. Virus-derived small interfering RNAs (vsiRNAs) are established components of the RNA silencing mechanism for incurring anti-viral resistance in plants. Some databases like siRNAdb, HIVsirDB and VIRsiRNAdb are available online pertaining to siRNAs as well as vsiRNAs generated during viral infection in humans; however, currently there is a lack of repository for plant exclusive vsiRNAs. We have developed `PVsiRNAdb (http://www.nipgr.res.in/PVsiRNAdb)', a manually curated plant-exclusive database harboring information related to vsiRNAs found in different virus-infected plants collected by exhaustive data mining of published literature so far. This database contains a total of 322 214 entries and 282 549 unique sequences of vsiRNAs. In PVsiRNAdb, detailed and comprehensive information is available for each vsiRNA sequence. Apart from the core information consisting of plant, tissue, virus name and vsiRNA sequence, additional information of each vsiRNAs (map position, length, coordinates, strand information and predicted structure) may be of high utility to the user. Different types of search and browse modules with three different tools namely BLAST, Smith-Waterman Align and Mapping are provided at PVsiRNAdb. Thus, this database being one of its kind will surely be of much use to molecular biologists for exploring the complex viral genetics and genomics, viral-host interactions and beneficial to the scientific community and can prove to be very advantageous in the field of agriculture for producing viral resistance transgenic crops.

Deep sequencing of mycovirus-derived small RNAs from Botrytis species

RNA silencing is an ancient regulatory mechanism operating in all eukaryotic cells. In fungi, it was first discovered in Neurospora crassa, although its potential as a defence mechanism against mycoviruses was first reported in Cryphonectria parasitica and, later, in several fungal species. There is little evidence of the antiviral potential of RNA silencing in the phytopathogenic species of the fungal genus Botrytis. Moreover, little is known about the RNA silencing components in these fungi, although the analysis of public genome databases identified two Dicer-like genes in B. cinerea, as in most of the ascomycetes sequenced to date. In this work, we used deep sequencing to study the virus-derived small RNA (vsiRNA) populations from different mycoviruses infecting field isolates of Botrytis spp. The mycoviruses under study belong to different genera and species, and have different types of genome [double-stranded RNA (dsRNA), (+)single-stranded RNA (ssRNA) and (-)ssRNA]. In general, vsiRNAs derived from mycoviruses are mostly of 21, 20 and 22 nucleotides in length, possess sense or antisense orientation, either in a similar ratio or with a predominance of sense polarity depending on the virus species, have predominantly U at their 5' end, and are unevenly distributed along the viral genome, showing conspicuous hotspots of vsiRNA accumulation. These characteristics reveal striking similarities with vsiRNAs produced by plant viruses, suggesting similar pathways of viral targeting in plants and fungi. We have shown that the fungal RNA silencing machinery acts against the mycoviruses used in this work in a similar manner independent of their viral or fungal origin.

The Shift of the Intestinal Microbiome in the Innate Immunity-Deficient Mutant rde-1 Strain of C. elegans upon Orsay Virus Infection

The status of intestinal microbiota is a determinant of host health. However, the alteration of the gut microbiota caused by the innate immune response to virus infection is unclear. Caenorhabditis elegans and its natural virus Orsay provide an excellent model of host–virus interactions. We evaluated the intestinal microbial community complexity of the wild-type N2 and the innate immunity-deficient mutant rde-1 (ne219) strains of C. elegans upon Orsay virus infection. The gut microbiota diversity was decreased in rde-1 (ne219) mutant animals, and a large number of genes were associated with the difference between infected and uninfected rde-1 (ne219) mutant animals. Therefore, this study provides the first evaluation of the alterations caused by Orsay virus on intestinal microbiota in wildtype and innate immunity-deficient animals using C. elegans as the model species. Our findings indicate that virus infection may alters the microbiome in animals with defective immune response.

Synergistic infection of BrYV and PEMV 2 increases the accumulations of both BrYV and BrYV-derived siRNAs in Nicotiana benthamiana

Viral synergism is caused by co-infection of two unrelated viruses, leading to more severe symptoms or increased titres of one or both viruses. Synergistic infection of phloem-restricted poleroviruses and umbraviruses has destructive effects on crop plants. The mechanism underlying this synergy remains elusive. In our study, synergism was observed in co-infections of a polerovirus Brassica yellows virus (BrYV) and an umbravirus Pea enation mosaic virus 2 (PEMV 2) on Nicotiana benthamiana, which led to (1) increased titres of BrYV, (2) appearance of severe symptoms, (3) gain of mechanical transmission capacity of BrYV, (4) broader distribution of BrYV to non-vascular tissues. Besides, profiles of virus-derived small interfering RNAs (vsiRNAs) from BrYV and PEMV 2 in singly and doubly infected plants were obtained by small RNA deep sequencing. Our results showed that accumulation of BrYV vsiRNAs increased tremendously and ratio of positive to negative strand BrYV vsiRNAs differed between singly infected and co-infected plants. Positions to which the BrYV vsiRNAs mapped to the viral genome varied considerably during synergistic infection. Moreover, target genes of vsiRNAs were predicted and annotated. Our results revealed the synergistic characteristics during co-infection of BrYV and PEMV 2, and implied possible effects of synergism have on vsiRNAs.

Caenorhabditis elegans RIG-I Homolog Mediates Antiviral RNA Interference Downstream of Dicer-Dependent Biogenesis of Viral Small Interfering RNAs

Dicer enzymes process virus-specific double-stranded RNA (dsRNA) into small interfering RNAs (siRNAs) to initiate specific antiviral defense by related RNA interference (RNAi) pathways in plants, insects, nematodes, and mammals. Antiviral RNAi in Caenorhabditis elegans requires Dicer-related helicase 1 (DRH-1), not found in plants and insects but highly homologous to mammalian retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs), intracellular viral RNA sensors that trigger innate immunity against RNA virus infection. However, it remains unclear if DRH-1 acts analogously to initiate antiviral RNAi in C. elegans. Here, we performed a forward genetic screen to characterize antiviral RNAi in C. elegans. Using a mapping-by-sequencing strategy, we uncovered four loss-of-function alleles of drh-1, three of which caused mutations in the helicase and C-terminal domains conserved in RLRs. Deep sequencing of small RNAs revealed an abundant population of Dicer-dependent virus-derived small interfering RNAs (vsiRNAs) in drh-1 single and double mutant animals after infection with Orsay virus, a positive-strand RNA virus. These findings provide further genetic evidence for the antiviral function of DRH-1 and illustrate that DRH-1 is not essential for the sensing and Dicer-mediated processing of the viral dsRNA replicative intermediates. Interestingly, vsiRNAs produced by drh-1 mutants were mapped overwhelmingly to the terminal regions of the viral genomic RNAs, in contrast to random distribution of vsiRNA hot spots when DRH-1 is functional. As RIG-I translocates on long dsRNA and DRH-1 exists in a complex with Dicer, we propose that DRH-1 facilitates the biogenesis of vsiRNAs in nematodes by catalyzing translocation of the Dicer complex on the viral long dsRNA precursors.

The helicase and C-terminal domains of mammalian RLRs sense intracellular viral RNAs to initiate the interferon-regulated innate immunity against RNA virus infection. Both of the domains from human RIG-I can substitute for the corresponding domains of DRH-1 to mediate antiviral RNAi in C. elegans, suggesting an analogous role for DRH-1 as an intracellular dsRNA sensor to initiate antiviral RNAi. Here, we developed a forward genetic screen for the identification of host factors required for antiviral RNAi in C. elegans. Characterization of four distinct drh-1 mutants obtained from the screen revealed that DRH-1 did not function to initiate antiviral RNAi. We show that DRH-1 acted in a downstream step to enhance Dicer-dependent biogenesis of viral siRNAs in C. elegans. As mammals produce Dicer-dependent viral siRNAs to target RNA viruses, our findings suggest a possible role for mammalian RLRs and interferon signaling in the biogenesis of viral siRNAs.

Characteristics of siRNAs derived from Southern rice black-streaked dwarf virus in infected rice and their potential role in host gene regulation

BACKGROUND

Virus-derived siRNAs (vsiRNAs)-mediated RNA silencing plays important roles in interaction between plant viruses and their hosts. Southern rice black-streaked dwarf virus (SRBSDV) is a newly emerged devastating rice reovirus with ten dsRNA genomic segments. The characteristics of SRBSDV-derived siRNAs and their biological implications in SRBSDV-rice interaction remain unexplored.

METHODS

VsiRNAs profiling from SRBSDV-infected rice samples was done via small RNA deep sequencing. The putative rice targets of abundantly expressed vsiRNAs were bioinformatically predicted and subjected to functional annotation. Differential expression analysis of rice targets and RNA silencing components between infected and healthy samples was done using RT-qPCR.

RESULTS

The vsiRNA was barely detectable at 14 days post infection (dpi) but abundantly present along with elevated expression level of the viral genome at 28 dpi. From the 28-dpi sample, 70,878 reads of 18 ~ 30-nt vsiRNAs were recognized (which mostly were 21-nt and 22-nt), covering 75 ~ 91% of the length of the ten genomic segments respectively. 86% of the vsiRNAs had a <50% GC content and 79% of them were 5'-uridylated or adenylated. The production of vsiRNAs had no strand polarity but varied among segment origins. Each segment had a few hotspot regions where vsiRNAs of high abundance were produced. 151 most abundant vsiRNAs were predicted to target 844 rice genes, including several types of host resistance or pathogenesis related genes encoding F-box/LRR proteins, receptor-like protein kinases, universal stress proteins, tobamovirus multiplication proteins, and RNA silencing components OsDCL2a and OsAGO17 respectively, some of which showed down regulation in infected plants in RT-qPCR. GO and KEGG classification showed that a majority of the predicted targets were related to cell parts and cellular processes and involved in carbohydrate metabolism, translation, and signal transduction. The silencing component genes OsDCL2a, OsDCL2b, OsDCL4, and OsAGO18 were down regulated, while OsAGO1d, OsAGO2, OsRDR1 and OsRDR6 were up regulated, significantly, upon SRBSDV infection.

CONCLUSIONS

SRBSDV can regulate the expression of rice RNA silencing pathway components and the virus might compromise host defense and influence host pathogenesis via siRNA pathways.

Characteristics of siRNAs derived from Southern rice black-streaked dwarf virus in infected rice and their potential role in host gene regulation

Background

Virus-derived siRNAs (vsiRNAs)-mediated RNA silencing plays important roles in interaction between plant viruses and their hosts. Southern rice black-streaked dwarf virus (SRBSDV) is a newly emerged devastating rice reovirus with ten dsRNA genomic segments. The characteristics of SRBSDV-derived siRNAs and their biological implications in SRBSDV-rice interaction remain unexplored.

Methods

VsiRNAs profiling from SRBSDV-infected rice samples was done via small RNA deep sequencing. The putative rice targets of abundantly expressed vsiRNAs were bioinformatically predicted and subjected to functional annotation. Differential expression analysis of rice targets and RNA silencing components between infected and healthy samples was done using RT-qPCR.

Results

The vsiRNA was barely detectable at 14 days post infection (dpi) but abundantly present along with elevated expression level of the viral genome at 28 dpi. From the 28-dpi sample, 70,878 reads of 18 ~ 30-nt vsiRNAs were recognized (which mostly were 21-nt and 22-nt), covering 75 ~ 91% of the length of the ten genomic segments respectively. 86% of the vsiRNAs had a <50% GC content and 79% of them were 5’-uridylated or adenylated. The production of vsiRNAs had no strand polarity but varied among segment origins. Each segment had a few hotspot regions where vsiRNAs of high abundance were produced. 151 most abundant vsiRNAs were predicted to target 844 rice genes, including several types of host resistance or pathogenesis related genes encoding F-box/LRR proteins, receptor-like protein kinases, universal stress proteins, tobamovirus multiplication proteins, and RNA silencing components OsDCL2a and OsAGO17 respectively, some of which showed down regulation in infected plants in RT-qPCR. GO and KEGG classification showed that a majority of the predicted targets were related to cell parts and cellular processes and involved in carbohydrate metabolism, translation, and signal transduction. The silencing component genes OsDCL2a, OsDCL2b, OsDCL4, and OsAGO18 were down regulated, while OsAGO1d, OsAGO2, OsRDR1 and OsRDR6 were up regulated, significantly, upon SRBSDV infection.

Conclusions

SRBSDV can regulate the expression of rice RNA silencing pathway components and the virus might compromise host defense and influence host pathogenesis via siRNA pathways.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-017-0699-3) contains supplementary material, which is available to authorized users.

Small RNA Based Genetic Engineering for Plant Viral Resistance: Application in Crop Protection

Small RNAs regulate a large set of gene expression in all plants and constitute a natural immunity against viruses. Small RNA based genetic engineering (SRGE) technology had been explored for crop protection against viruses for nearly 30 years. Viral resistance has been developed in diverse crops with SRGE technology and a few viral resistant crops have been approved for commercial release. In this review we summarized the efforts generating viral resistance with SRGE in different crops, analyzed the evolution of the technology, its efficacy in different crops for different viruses and its application status in different crops. The challenge and potential solution for application of SRGE in crop protection are also discussed.

CaRDR1, an RNA-Dependent RNA Polymerase Plays a Positive Role in Pepper Resistance against TMV

RNA silencing functions as a major natural antiviral defense mechanism in plants. RNA-dependent RNA polymerases (RDRs) that catalyze the synthesis of double-stranded RNAs, are considered as a fundamental element in RNA silencing pathways. In Arabidopsis thaliana, RDR1, 2 and 6 play important roles in anti-viral RNA silencing. Expression of RDR1 can be elevated following plant treatment with defense hormones and virus infection. RDR1 has been studied in several crop species, but not in pepper (Capsicum annuum L.). Here, a RDR1 gene was isolated from Capsicum annuum L., designated as CaRDR1. The full-length cDNA of CaRDR1 was 3,351 bp, encoding a 1,116-amino acid protein, which contains conserved regions, such as the most remarkable motif DLDGD. The transcripts of CaRDR1 could be induced by salicylic acid (SA), abscisic acid (ABA), H₂O₂, and tobacco mosaic virus (TMV). Silencing of CaRDR1 in pepper resulted in increased susceptibility to TMV as evident by severe symptom, increased of TMV-CP transcript, higher malondialdehyde (MDA) content and lower antioxidant enzymes activities compared with that of control plants. CaRDR1-overexpressing in Nicotiana benthamiana showed mild disease symptom and reduced TMV-CP transcripts than that of empty vector (EV) following TMV inoculation. The RNA silencing related genes, including NbAGO2, NbDCL2, NbDCL3, and NbDCL4 elevated expression in overexpressed plants. Alternative oxidase (AOX), the terminal oxidase of the cyanide (CN)-resistant alternative respiratory pathway, catalyze oxygen-dependent oxidation of ubiquinol in plants. It has an important function in plant defense against TMV. In addition, CaRDR1 overexpression promoted the expression of NbAOX1a and NbAOX1b. In conclusion, these results suggest that CaRDR1 plays a positive role in TMV resistance by regulating antioxidant enzymes activities and RNA silencing-related genes expression to suppress the replication and movement of TMV.

Incomplete DRB4-dependence of the DCL4-mediated antiviral defense

The double-stranded RNA-binding protein DRB4 of Arabidopsis was shown previously to contribute to the DICER-LIKE 4 (DCL4)-mediated biogenesis of viral small interfering RNAs (vsiRNAs) of 21 nucleotides (nt) in size. However, it is unclear whether all 21-nt vsiRNAs are dependent on this DRB4-DCL4 partnership. To resolve this question, we generated dcl2drb4 and dcl4drb4 double knockout mutants, and subjected them to infections with CPB-CC-PDS, a turnip crinkle virus mutant capable of inducing silencing of the PHYTOENE DESATURASE gene. The dcl2drb4 double knockouts caused a far smaller loss of antiviral silencing than dcl2dcl4. In addition, although both drb4 and dcl4 single mutants permitted a consistent (but small) increase in viral RNA levels, the drb4 mutant correlated with a less pronounced reduction of 21-nt vsiRNAs. Therefore, a substantial subset of DCL4 antiviral activity is DRB4-independent, and may involve other DRB proteins that compensate for loss of DRB4.

Novel Functional Genomics Approaches: A Promising Future in the Combat Against Plant Viruses

Advances in functional genomics and genome editing approaches have provided new opportunities and potential to accelerate plant virus control efforts through modification of host and viral genomes in a precise and predictable manner. Here, we discuss application of RNA-based technologies, including artificial micro RNA, transacting small interfering RNA, and Cas9 (clustered regularly interspaced short palindromic repeat-associated protein 9), which are currently being successfully deployed in generating virus-resistant plants. We further discuss the reverse genetics approach, targeting induced local lesions in genomes (TILLING) and its variant, known as EcoTILLING, that are used in the identification of plant virus recessive resistance gene alleles. In addition to describing specific applications of these technologies in plant virus control, this review discusses their advantages and limitations.

Arabidopsis RNA Polymerases IV and V Are Required To Establish H3K9 Methylation, but Not Cytosine Methylation, on Geminivirus Chromatin

In plants, RNA-directed DNA methylation (RdDM) employs small RNAs to target enzymes that methylate cytosine residues. Cytosine methylation and dimethylation of histone 3 lysine 9 (H3K9me2) are often linked. Together they condition an epigenetic defense that results in chromatin compaction and transcriptional silencing of transposons and viral chromatin. Canonical RdDM (Pol IV-RdDM), involving RNA polymerases IV and V (Pol IV and Pol V), was believed to be necessary to establish cytosine methylation, which in turn could recruit H3K9 methyltransferases. However, recent studies have revealed that a pathway involving Pol II and RNA-dependent RNA polymerase 6 (RDR6) (RDR6-RdDM) is likely responsible for establishing cytosine methylation at naive loci, while Pol IV-RdDM acts to reinforce and maintain it. We used the geminivirus Beet curly top virus (BCTV) as a model to examine the roles of Pol IV and Pol V in establishing repressive viral chromatin methylation. As geminivirus chromatin is formed de novo in infected cells, these viruses are unique models for processes involved in the establishment of epigenetic marks. We confirm that Pol IV and Pol V are not needed to establish viral DNA methylation but are essential for its amplification. Remarkably, however, both Pol IV and Pol V are required for deposition of H3K9me2 on viral chromatin. Our findings suggest that cytosine methylation alone is not sufficient to trigger de novo deposition of H3K9me2 and further that Pol IV-RdDM is responsible for recruiting H3K9 methyltransferases to viral chromatin.

IMPORTANCE

In plants, RNA-directed DNA methylation (RdDM) uses small RNAs to target cytosine methylation, which is often linked to H3K9me2. These epigenetic marks silence transposable elements and DNA virus genomes, but how they are established is not well understood. Canonical RdDM, involving Pol IV and Pol V, was thought to establish cytosine methylation that in turn could recruit H3K9 methyltransferases, but recent studies compel a reevaluation of this view. We used BCTV to investigate the roles of Pol IV and Pol V in chromatin methylation. We found that both are needed to amplify, but not to establish, DNA methylation. However, both are required for deposition of H3K9me2. Our findings suggest that cytosine methylation is not sufficient to recruit H3K9 methyltransferases to naive viral chromatin and further that Pol IV-RdDM is responsible.

Putative RNA-directed adaptive mutations in cancer evolution

Understanding the molecular mechanisms behind the capacity of cancer cells to adapt to the tumor microenvironment and to anticancer therapies is a major challenge. In this context, cancer is believed to be an evolutionary process where random mutations and the selection process shape the mutational pattern and phenotype of cancer cells. This article challenges the notion of randomness of some cancer-associated mutations by describing molecular mechanisms involving stress-mediated biogenesis of mRNA-derived small RNAs able to target and increase the local mutation rate of the genomic loci they originate from. It is proposed that the probability of some mutations at specific loci could be increased in a stress-specific and RNA-depending manner. This would increase the probability of generating mutations that could alleviate stress situations, such as those triggered by anticancer drugs. Such a mechanism is made possible because tumor- and anticancer drug-associated stress situations trigger both cellular reprogramming and inflammation, which leads cancer cells to express molecular tools allowing them to “attack” and mutate their own genome in an RNA-directed manner.

Arabidopsis RNA Polymerases IV and V Are Required To Establish H3K9 Methylation, but Not Cytosine Methylation, on Geminivirus Chromatin

In plants, RNA-directed DNA methylation (RdDM) employs small RNAs to target enzymes that methylate cytosine residues. Cytosine methylation and dimethylation of histone 3 lysine 9 (H3K9me2) are often linked. Together they condition an epigenetic defense that results in chromatin compaction and transcriptional silencing of transposons and viral chromatin. Canonical RdDM (Pol IV-RdDM), involving RNA polymerases IV and V (Pol IV and Pol V), was believed to be necessary to establish cytosine methylation, which in turn could recruit H3K9 methyltransferases. However, recent studies have revealed that a pathway involving Pol II and RNA-dependent RNA polymerase 6 (RDR6) (RDR6-RdDM) is likely responsible for establishing cytosine methylation at naive loci, while Pol IV-RdDM acts to reinforce and maintain it. We used the geminivirus Beet curly top virus (BCTV) as a model to examine the roles of Pol IV and Pol V in establishing repressive viral chromatin methylation. As geminivirus chromatin is formed de novo in infected cells, these viruses are unique models for processes involved in the establishment of epigenetic marks. We confirm that Pol IV and Pol V are not needed to establish viral DNA methylation but are essential for its amplification. Remarkably, however, both Pol IV and Pol V are required for deposition of H3K9me2 on viral chromatin. Our findings suggest that cytosine methylation alone is not sufficient to trigger de novo deposition of H3K9me2 and further that Pol IV-RdDM is responsible for recruiting H3K9 methyltransferases to viral chromatin.

IMPORTANCE

In plants, RNA-directed DNA methylation (RdDM) uses small RNAs to target cytosine methylation, which is often linked to H3K9me2. These epigenetic marks silence transposable elements and DNA virus genomes, but how they are established is not well understood. Canonical RdDM, involving Pol IV and Pol V, was thought to establish cytosine methylation that in turn could recruit H3K9 methyltransferases, but recent studies compel a reevaluation of this view. We used BCTV to investigate the roles of Pol IV and Pol V in chromatin methylation. We found that both are needed to amplify, but not to establish, DNA methylation. However, both are required for deposition of H3K9me2. Our findings suggest that cytosine methylation is not sufficient to recruit H3K9 methyltransferases to naive viral chromatin and further that Pol IV-RdDM is responsible.

Virus-derived small RNAs in the penaeid shrimp Fenneropenaeus chinensis during acute infection of the DNA virus WSSV

Small interfering RNAs (siRNAs) and microRNAs (miRNAs) are two classes of small RNAs (sRNAs) that are critical for virus-host interplay via the RNA interference (RNAi) pathway. One virus-derived siRNA and numerous miRNAs has been reported for the double-stranded DNA virus white spot syndrome virus (WSSV), however, the expression profiles of these different types of sRNAs have not been assessed. Here, by sequencing the sRNAs and mRNAs of WSSV-infected Chinese shrimp (Fenneropenaeus chinensis), we found that the viral transcripts were universally targeted by WSSV-derived siRNAs, supporting a pivotal role for RNAi in the anti-viral immunity of shrimp. The genesis of WSSV-derived siRNAs was associated with long RNA structures. Moreover, by separating miRNAs from siRNAs, 12 WSSV miRNAs were identified. Investigation of conserved viral miRNA targets in different host species indicated the involvement of viral miRNAs in host immune responses. Collectively, our data provide new insights into the role of the RNAi pathway in the interplay between DNA viruses and crustaceans.

Properties of satellite tobacco mosaic virus phenotypes expressed in the presence and absence of helper virus

In this study, we assembled an Agrobacterium-based transient expression system for the ectopic expression of Satellite tobacco mosaic virus (STMV) (+) or (-) transcripts and their biological activity was confirmed when Nicotiana benthamiana plants were co-expressed with helper Tobacco mosaic virus replicase. Characterization of STMV in the presence and absence of its HV revealed: (i) HV-dependent expression of STMV (+) in N. benthamiana, but not in N. tabacum, generated a replication-deficient but translation and encapsidation competent variant lacking the highly conserved 3' 150 nucleotides (nt) (STMVΔ150); (ii) mutational analysis demonstrated that a conserved 3' stem-loop structure in wild type and STMVΔ150 located between nt 874 and 897 is essential for translation of CP; (iii) helper virus-independent expression of CP from wt STMV was competent for the assembly of empty aberrant virion-like particles; whereas, CP translated from STMVΔ150 resulted in disorganized CP aggregates suggesting a role for the 3'tRNA-like structure in STMV assembly.

MicroRNA-Dependent Transcriptional Silencing of Transposable Elements in Drosophila Follicle Cells

RNA interference-related silencing mechanisms concern very diverse and distinct biological processes, from gene regulation (via the microRNA pathway) to defense against molecular parasites (through the small interfering RNA and the Piwi-interacting RNA pathways). Small non-coding RNAs serve as specificity factors that guide effector proteins to ribonucleic acid targets via base-pairing interactions, to achieve transcriptional or post-transcriptional regulation. Because of the small sequence complementarity required for microRNA-dependent post-transcriptional regulation, thousands of microRNA (miRNA) putative targets have been annotated in Drosophila. In Drosophila somatic ovarian cells, genomic parasites, such as transposable elements (TEs), are transcriptionally repressed by chromatin changes induced by Piwi-interacting RNAs (piRNAs) that prevent them from invading the germinal genome. Here we show, for the first time, that a functional miRNA pathway is required for the piRNA-mediated transcriptional silencing of TEs in this tissue. Global miRNA depletion, caused by tissue- and stage-specific knock down of drosha (involved in miRNA biogenesis), AGO1 or gawky (both responsible for miRNA activity), resulted in loss of TE-derived piRNAs and chromatin-mediated transcriptional de-silencing of TEs. This specific TE de-repression was also observed upon individual titration (by expression of the complementary miRNA sponge) of two miRNAs (miR-14 and miR-34) as well as in a miR-14 loss-of-function mutant background. Interestingly, the miRNA defects differentially affected TE- and 3' UTR-derived piRNAs. To our knowledge, this is the first indication of possible differences in the biogenesis or stability of TE- and 3' UTR-derived piRNAs. This work is one of the examples of detectable phenotypes caused by loss of individual miRNAs in Drosophila and the first genetic evidence that miRNAs have a role in the maintenance of genome stability via piRNA-mediated TE repression.