Phased, secondary, small interfering RNAs in posttranscriptional regulatory networks

syn-tasiRnas targeting the coat protein of potato virus Y confer antiviral resistance in Nicotiana benthamiana

Trans-acting small interfering RNAs (tasiRNAs) are 21-nt phased (phased siRNAs) resulting from successive DCL-catalyzed processing from the end of a double-stranded RNA substrate originating from the RDR of an AGO-catalyzed cleaved RNA at a micro RNA target site. Plant tasiRNAs have been synthesized to produce synthetic tasiRNAs (syn-tasiRNAs) targeting viral RNAs that confer viral resistance. In this study, we engineered syn-tasiRNAs to target potato virus Y (PVY) infection by replacing five native siRNAs of TAS1c with 210-bp fragments from the coat protein (CP) region of the PVY genome. The results showed that the transient expression of syn-tasiR-CPpvy2 in Nicotiana benthamiana (N. benthamiana) plants conferred antiviral resistance, supported by the absence of PVY infection symptoms and viral accumulation. This indicated that syn-tasiR-CPpvy2 successfully targeted and silenced the PVY CP gene, effectively inhibiting viral infection. syn-tasiR-CPpvy1 displayed attenuated symptoms and decreased viral accumulation in these plants However, severe symptoms of PVY infection and a similar amount of viral accumulation as the control were observed in plants expressing syn-tasiR-CPpvy3. syn-tasiR-CPpvy/pvx, which targets both PVY and potato virus X (PVX), was engineered using a single precursor. After the transient expression of syn-tasiR-CPpvy/pvx3 and syn-tasiR-CPpvy/pvx5 in N. benthamiana, the plants were resistant to both PVY and PVX. These results suggested that engineered syn-tasiRNAs could not only specifically induce antiviral resistance against one target virus but could also be designed for multi-targeted silencing of different viruses, thereby preventing complex virus infection in plants.

DNA methylation analysis reveals local changes in resistant and susceptible soybean lines in response to Phytophthora sansomeana

Phytophthora sansomeana is an emerging oomycete pathogen causing root rot in many agricultural species including soybean. However, as of now, only one potential resistance gene has been identified in soybean, and our understanding of how genetic and epigenetic regulation in soybean contributes to responses against this pathogen remains largely unknown. In this study, we performed whole genome bisulfite sequencing (WGBS) on two soybean lines, Colfax (resistant) and Williams 82 (susceptible), in response to P. sansomeana at two time points: 4 and 16 hours post-inoculation to compare their methylation changes. Our findings revealed that there were no significant changes in genome-wide CG, CHG (H = A, T, or C), and CHH methylation. However, we observed local methylation changes, specially an increase in CHH methylation around genes and transposable elements (TEs) after inoculation, which occurred earlier in the susceptible line and later in the resistant line. After inoculation, we identified differentially methylated regions (DMRs) in both Colfax and Williams 82, with a predominant presence in TEs. Notably, our data also indicated that more TEs exhibited changes in their methylomes in the susceptible line compared to the resistant line. Furthermore, we discovered 837 DMRs within or flanking 772 differentially expressed genes (DEGs) in Colfax and 166 DMRs within or flanking 138 DEGs in Williams 82. These DEGs had diverse functions, with Colfax primarily showing involvement in metabolic process, defense response, plant and pathogen interaction, anion and nucleotide binding, and catalytic activity, while Williams 82 exhibited a significant association with photosynthesis. These findings suggest distinct molecular responses to P. sansomeana infection in the resistant and susceptible soybean lines.

Structural determinants in the miRNA/miRNA* duplex and the DCL1 PAZ domain for precise and efficient plant miRNA processing

The accessory proteins Hyponastic-like 1 (HYL1) and Serrated (SE) enhance the precise and efficient processing of miRNAs by Dicer-like 1 (DCL1), which is important for proper miRNA function. However, other factors determining the precision and efficiency of miRNA biogenesis are not well-known. Here, we found that an asymmetric bulge (AB) at the 3' end of miR-5p (produced from the 5' arm of the pre-miRNA) reduced the precision of the second cleavage, whereas an AB at other sites of miR-5p mainly affected the accumulation level of miR-5p in transient expression in Nicotiana benthamiana. In contrast, many ABs in miR-3p (produced from the 3' arm of the pre-miRNA) impose strong negative impact on the processing precision and the accumulation level of miR-5p in N. benthamiana. Arabidopsis DCL1/SE/HYL1 complex-mediated miRNA processing was reconstituted in Saccharomyces cerevisiae to further investigate AB-mediated interference with DCL1 processing. With this system, the positional effect of AB on miRNA processing was tested. The results showed that ABs on the middle of miR-5p have less of an impact on DCL1 cleavage efficiency and precision, whereas those on miR-3p or near the ends of miR-5p strongly reduce DCL1 cleavage activity, precision or both. Studies using the yeast miRNA processing system and transgenic Arabidopsis also revealed the importance of the interaction between the 2-nt 3' overhang of pre-miRNA and the 3' overhang binding pocket (3'BP) on the precision of the second cleavage reaction for many endogenous miRNAs. These findings provide new insights into the mechanism of miRNA biogenesis.

A phased small interfering RNA-derived pathway mediates lead stress tolerance in maize

Phased small interfering RNAs (phasiRNAs) are a distinct class of endogenous small interfering RNAs, which regulate plant growth, development, and environmental stress response. To determine the effect of phasiRNAs on maize (Zea mays L.) tolerance to lead (Pb) stress, the roots of 305 maize lines under Pb treatment were subjected to generation of individual databases of small RNAs. We identified 55 high-confidence phasiRNAs derived from 13 PHAS genes (genes producing phasiRNAs) in this maize panel, of which 41 derived from 9 PHAS loci were negatively correlated with Pb content in the roots. The potential targets of the 41 phasiRNAs were enriched in ion transport and import. Only the expression of PHAS_1 (ZmTAS3j, Trans-Acting Short Interference RNA3) was regulated by its cis-expression quantitative trait locus and thus affected the Pb content in the roots. Using the Nicotiana benthamiana transient expression system, 5'-rapid amplification of cDNA ends, and Arabidopsis heterologously expressed, we verified that ZmTAS3j was cleaved by zma-miR390 and thus generated tasiRNA targeting ARF genes (tasiARFs), and that the 5' and 3' zma-miR390 target sites of ZmTAS3j were both necessary for efficient biosynthesis and functional integrity of tasiARFs. We validated the involvement of the zma-miR390-ZmTAS3j-tasiARF-ZmARF3-ZmHMA3 pathway in Pb accumulation in maize seedlings using genetic, molecular, and cytological methods. Moreover, the increased Pb tolerance in ZmTAS3j-overexpressed lines was likely attributed to the zma-miR390-ZmTAS3j-tasiARF-ZmARF3-SAURs pathway, which elevated indole acetic acid levels and thus reactive oxygen species-scavenging capacity in maize roots. Our study reveals the importance of the TAS3-derived tasiRNA pathway in plant adaptation to Pb stress.

Long noncoding RNAs in ubiquitination, protein degradation, and human diseases

Protein stability and turnover is critical in normal cellular and physiological process and their misregulation may contribute to accumulation of unwanted proteins causing cellular malfunction, neurodegeneration, mitochondrial malfunction, and disrupted metabolism. Signaling mechanism associated with protein degradation is complex and is extensively studied. Many protein and enzyme machineries have been implicated in regulation of protein degradation. Despite these insights, our understanding of protein degradation mechanisms remains limited. Emerging studies suggest that long non-coding RNAs (lncRNAs) play critical roles in various cellular and physiological processes including metabolism, cellular homeostasis, and protein turnover. LncRNAs, being large nucleic acids (>200 nt long) can interact with various proteins and other nucleic acids and modulate protein structure and function leading to regulation of cell signaling processes. LncRNAs are widely distributed across cell types and may exhibit tissue specific expression. They are detected in body fluids including blood and urine. Their expressions are also altered in various human diseases including cancer, neurological disorders, immune disorder, and others. LncRNAs are being recognized as novel biomarkers and therapeutic targets. This review article focuses on the emerging role of noncoding RNAs (ncRNAs), particularly long noncoding RNAs (lncRNAs), in the regulation of protein polyubiquitination and proteasomal degradation.

Integrating multiregulatory analysis reveals the negative regulatory function of miR482a in the response of poplar to canker pathogen infection

Canker disease caused by the bacterium Lonsdalea populi is one of the most destructive diseases affecting poplar stems. However, the detailed stress response mechanisms of poplar have not been widely characterized. To explore the diverse regulatory RNA landscape and the function of key regulators in poplar subjected to L. populi stress, we integrated time-course experiment with mock-inoculation (CK) and inoculation (IN) with L. populi at the first, third, and sixth day (IN1, IN3, IN6) on Populus × euramericana cv. '74/76' (107), small RNA-seq, whole transcriptome-wide analysis, degradome analysis and transgenic experiments. A total of 98 differentially expressed (DE) miRNA, 17 974 DEmRNA, and 807 DElncRNA were identified in poplar infected by L. populi, presenting dynamic changes over the infection course. Regulatory networks among RNAs were further constructed. Notably, a network centered on ptc-miR482a in CK-vs-IN3 contained most DEGs. We show that miR482a and miR1448 are located in one transcript as a polycistron. Overexpression of pre-miR482a-miR1448 (OX482-1448) and pre-miR482a (OX482) increased poplar susceptibility to canker pathogen with reduced accumulation of reactive oxygen species, while the suppression of miR482a (STTM482) conferred poplar disease resistance. PHA7 was validated as the target of miR482a with degradome sequencing and tobacco transient co-transformation, its expression being downregulated in OX482-1448 and OX482 lines. Additionally, a series of phasiRNAs were triggered by miR482a targeting PHA7, forming regulatory cascades with more RLP, NBS-LRR, and PK genes, further verifying the defense function of miR482a. These findings provide insights for understanding the roles of ncRNAs and regulatory networks involved in poplar immunity.

Analysis of the miR482 Gene Family in Plants

MicroRNA482 (miR482) is a conserved microRNA family in plants, playing critical regulatory roles in different biological activities. Though the members of the miR482 gene family have been identified in plants, a systematic study has not been reported yet. In the present research, 140 mature sequences generated by 106 precursors were used for molecular characterization, phylogenetic analysis, and target gene prediction, and the competing endogenous RNA (ceRNA) network mediated by miR482 was summarized. The length of mature sequences ranged from 17 nt to 25 nt, with 22 nt being the most abundant, and the start and end of the mature sequences had a preference for uracil (U). By sequence multiplex comparison, it was found that the mature sequences of 5p were clustered into one group, and others were clustered into the other group. Phylogenetic analysis revealed that the 140 mature sequences were categorized into six groups. Meanwhile, all the precursor sequences formed a stable hairpin structure, and the 106 precursors were divided into five groups. However, the expression of miR482 varied significantly between different species and tissues. In total, 149 target genes were predicted and their functions focused on single-organism process, cellular process, and cell and cell part. The ceRNA network of miR482 in tomato, cotton, and peanut was summarized based on related publications. In conclusion, this research will provide a foundation for further understanding of the miR482 gene family.

The role of microRNAs in NBS-LRR gene expression and its implications for plant immunity and crop development

Plants evolved, over millions of years, complex defense systems against pathogens. Once infected, the interaction between pathogen effector molecules and host receptors triggers plant immune responses, which include apoptosis, systemic immune response, among others. An important protein family responsible for pathogen effector recognition is the nucleotide binding site-leucine repeat rich (NBS-LRR) proteins. The NBS-LRR gene family is the largest disease resistance gene class in plants. These proteins are widely distributed in vascular plants and have a complex multigenic cluster distribution in plant genomes. To counteract the genetic load of such a large gene family on fitness cost, plants evolved a mechanism using post transcriptional gene silencing induced by small RNAs, particularly microRNAs. For the NBS-LRR gene family, the small RNAs involved in this silencing mechanism are mainly the microRNA482/2118 superfamily. This suppression mechanism is relieved upon pathogen infection, thus allowing increased NBS-LRR expression and triggering plant immunity. In this review, we will discuss the biogenesis of microRNAs and secondary RNAs involved in this silencing mechanism, biochemical and structural features of NBS-LRR proteins in response to pathogen effectors and the evolution of microRNA-based silencing mechanism with a focus on the miR482/2118 family. Furthermore, the biotechnological manipulation of microRNA expression, using both transgenic or genome editing approaches to improve cultivated plants will be discussed, with a focus on the miR482/2118 family in soybean.

The miR172/TOE3 module regulates resistance to tobacco mosaic virus in tobacco

The outcome of certain plant-virus interaction is symptom recovery, which is accompanied with the emergence of asymptomatic tissues in which the virus accumulation decreased dramatically. This phenomenon shows the potential to reveal critical molecular factors for controlling viral disease. MicroRNAs act as master regulators in plant growth, development, and immunity. However, the mechanism by which miRNA participates in regulating symptom recovery remains largely unknown. Here, we reported that miR172 was scavenged in the recovered tissue of tobacco mosaic virus (TMV)-infected Nicotiana tabacum plants. Overexpression of miR172 promoted TMV infection, whereas silencing of miR172 inhibited TMV infection. Then, TARGET OF EAT3 (TOE3), an APETALA2 transcription factor, was identified as a downstream target of miR172. Overexpression of NtTOE3 significantly improved plant resistance to TMV infection, while knockout of NtTOE3 facilitated virus infection. Furthermore, transcriptome analysis indicated that TOE3 promoted the expression of defense-related genes, such as KL1 and MLP43. Overexpression of these genes conferred resistance of plant against TMV infection. Importantly, results of dual-luciferase assay, chromatin immunoprecipitation-quantitative PCR, and electrophoretic mobility shift assay proved that TOE3 activated the transcription of KL1 and MLP43 by binding their promoters. Moreover, overexpression of rTOE3 (the miR172-resistant form of TOE3) significantly reduced TMV accumulation compared to the overexpression of TOE3 (the normal form of TOE3) in miR172 overexpressing Nicotiana benthamiana plants. Taken together, our study reveals the pivotal role of miR172/TOE3 module in regulating plant immunity and in the establishment of recovery in virus-infected tobacco plants, elucidating a regulatory mechanism integrating plant growth, development, and immune response.

The microRNA miR482 regulates NBS-LRR genes in response to ALT1 infection in apple

Plants are frequently attacked by a variety of pathogens and thus have evolved a series of defense mechanisms, one important mechanism is resistance gene (R gene)-mediated disease resistance, but its expression is tightly regulated. NBS-LRR genes are the largest gene family of R genes. microRNAs (miRNAs) target to a number of NBS-LRR genes and trigger the production of phased small interfering RNAs (phasiRNAs) from these transcripts. phasiRNAs cis or trans regulate NBS-LRR genes, which can result in the repression of R gene expression. In this study, we screened for upregulated miR482 in the susceptible apple cultivar 'Golden Delicious' (GD) after inoculation with the fungal pathogen Alternaria alternata f. sp. mali (ALT1). Additionally, through combined degradome sequencing, we identified a gene targeted by miR482, named MdTNL1, a gene encoding a TIR-NBS-LRR (Toll/interleukin1 receptor-nucleotide binding site-leucine-rich repeat) protein. This gene exhibited a significant down-regulation post ALT1 inoculation, suggesting an impact on gene expression mediated by miRNA regulation. miR482 could cleave MdTNL1 and generate phasiRNAs at the cleavage site. We found that overexpression of miR482 inhibited the expression of MdTNL1 and thus reduced the disease resistance of GD, while silencing of miR482 increased the expression of MdTNL1 and thus improved the disease resistance of GD. This work elucidates key mechanisms underlying the immune response to Alternaria infection in apple. Identification of the resistance genes involved will enable molecular breeding for prevention and control of Alternaria leaf spot disease in this important fruit crop.

After silencing suppression: miRNA targets strike back

Within the continuous tug-of-war between plants and microbes, RNA silencing stands out as a key battleground. Pathogens, in their quest to colonize host plants, have evolved a diverse arsenal of silencing suppressors as a common strategy to undermine the host's RNA silencing-based defenses. When RNA silencing malfunctions in the host, genes that are usually targeted and silenced by microRNAs (miRNAs) become active and can contribute to the reprogramming of host cells, providing an additional defense mechanism. A growing body of evidence suggests that miRNAs may act as intracellular sensors to enable a rapid response to pathogen threats. Herein we review how plant miRNA targets play a crucial role in immune responses against different pathogens.

Genome-wide identification of phasiRNAs in Arabidopsis thaliana, and insights into biogenesis, temperature sensitivity, and organ specificity

The knowledge of biogenesis and target regulation of the phased small interfering RNAs (phasiRNAs) needs continuous update, since the phasiRNA loci are dynamically evolved in plants. Here, hundreds of phasiRNA loci of Arabidopsis thaliana were identified in distinct tissues and under different temperature. In flowers, most of the 24-nt loci are RNA-dependent RNA polymerase 2 (RDR2)-dependent, while the 21-nt loci are RDR6-dependent. Among the RDR-dependent loci, a significant portion is Dicer-like 1-dependent, indicating the involvement of microRNAs in their expression. Besides, two TAS candidates were discovered. Some interesting features of the phasiRNA loci were observed, such as the strong strand bias of phasiRNA generation, and the capacity of one locus for producing phasiRNAs by different increments. Both organ specificity and temperature sensitivity were observed for phasiRNA expression. In leaves, the TAS genes are highly activated under low temperature. Several trans-acting siRNA-target pairs are also temperature-sensitive. In many cases, the phasiRNA expression patterns correlate well with those of the processing signals. Analysis of the rRNA-depleted degradome uncovered several phasiRNA loci to be RNA polymerase II-independent. Our results should advance the understanding on phasiRNA biogenesis and regulation in plants.

Comparative analysis, diversification, and functional validation of plant nucleotide-binding site domain genes

Nucleotide-binding site (NBS) domain genes are one of the superfamily of resistance genes involved in plant responses to pathogens. The current study identified 12,820 NBS-domain-containing genes across 34 species covering from mosses to monocots and dicots. These identified genes are classified into 168 classes with several novel domain architecture patterns encompassing significant diversity among plant species. Several classical (NBS, NBS-LRR, TIR-NBS, TIR-NBS-LRR, etc.) and species-specific structural patterns (TIR-NBS-TIR-Cupin_1-Cupin_1, TIR-NBS-Prenyltransf, Sugar_tr-NBS etc.) were discovered. We observed 603 orthogroups (OGs) with some core (most common orthogroups; OG0, OG1, OG2, etc.) and unique (highly specific to species; OG80, OG82, etc.) OGs with tandem duplications. The expression profiling presented the putative upregulation of OG2, OG6, and OG15 in different tissues under various biotic and abiotic stresses in susceptible and tolerant plants to cotton leaf curl disease (CLCuD). The genetic variation between susceptible (Coker 312) and tolerant (Mac7) Gossypium hirsutum accessions identified several unique variants in NBS genes of Mac7 (6583 variants) and Coker312 (5173 variants). The protein-ligand and proteins-protein interaction showed a strong interaction of some putative NBS proteins with ADP/ATP and different core proteins of the cotton leaf curl disease virus. The silencing of GaNBS (OG2) in resistant cotton through virus-induced gene silencing (VIGS) demonstrated its putative role in virus tittering. The presented study will be further helpful in understanding the plant adaptation mechanism.

Two adjacent NLR genes conferring quantitative resistance to clubroot disease in Arabidopsis are regulated by a stably inherited epiallelic variation

Clubroot caused by the protist Plasmodiophora brassicae is a major disease affecting cultivated Brassicaceae. Using a combination of quantitative trait locus (QTL) fine mapping, CRISPR-Cas9 validation, and extensive analyses of DNA sequence and methylation patterns, we revealed that the two adjacent neighboring NLR (nucleotide-binding and leucine-rich repeat) genes AT5G47260 and AT5G47280 cooperate in controlling broad-spectrum quantitative partial resistance to the root pathogen P. brassicae in Arabidopsis and that they are epigenetically regulated. The variation in DNA methylation is not associated with any nucleotide variation or any transposable element presence/absence variants and is stably inherited. Variations in DNA methylation at the Pb-At5.2 QTL are widespread across Arabidopsis accessions and correlate negatively with variations in expression of the two genes. Our study demonstrates that natural, stable, and transgenerationally inherited epigenetic variations can play an important role in shaping resistance to plant pathogens by modulating the expression of immune receptors.

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

The unusual predominance of maintenance DNA methylation in Spirodela polyrhiza

Duckweeds are among the fastest reproducing plants, able to clonally divide at exponential rates. However, the genetic and epigenetic impact of clonality on plant genomes is poorly understood. 5-methylcytosine (5mC) is a modified base often described as necessary for the proper regulation of certain genes and transposons and for the maintenance of genome integrity in plants. However, the extent of this dogma is limited by the current phylogenetic sampling of land plant species diversity. Here we analyzed DNA methylomes, small RNAs, mRNA-seq, and H3K9me2 histone modification for Spirodela polyrhiza. S. polyrhiza has lost highly conserved genes involved in de novo methylation of DNA at sites often associated with repetitive DNA, and within genes, however, symmetrical DNA methylation and heterochromatin are maintained during cell division at certain transposons and repeats. Consequently, small RNAs that normally guide methylation to silence repetitive DNA like retrotransposons are diminished. Despite the loss of a highly conserved methylation pathway, and the reduction of small RNAs that normally target repetitive DNA, transposons have not proliferated in the genome, perhaps due in part to the rapid, clonal growth lifestyle of duckweeds.

Review: Plant microRNAs in pathogen defense: A panacea or a piece of the puzzle?

Plant microRNAs (miRNAs) control key agronomic traits that are associated with their conserved role(s) in development. However, despite a multitude of studies, the utility of miRNAs in plant-pathogen resistance remains less certain. Reviewing the literature identifies three general classes of miRNAs regarding plant pathogen defense. Firstly, a number of evolutionary dynamic 22 nucleotide miRNA families that repress large numbers of plant immunity genes, either directly, or through triggering the biogenesis of secondary siRNAs. However, understanding of their role in defense and of their manipulation to enhance pathogen resistance are still lacking. Secondly, highly conserved miRNAs that indirectly impact disease resistance through their targets that are primarily regulating development or hormone signaling. Any alteration of these miRNAs usually results in pleiotropic impacts, which may alter disease resistance in some plant species, and against some pathogens. Thirdly, are the comparatively diverse and evolutionary dynamic set of non-conserved miRNAs, some of which contribute to pathogen resistance, but whose narrow evolutionary presence will likely restrict their utility. Therefore, reflecting the diverse and evolving nature of plant-pathogen interactions, a complex interplay of plant miRNAs with pathogen responses exists. Any miRNA-based solution for pathogen resistance will likely be highly specific, rather than a general panacea.

sRNAminer: A multifunctional toolkit for next-generation sequencing small RNA data mining in plants

Small RNAs (sRNAs), found extensively in plants, play an essential role in plant growth and development. Although various sRNA analysis tools have been developed for plants, the use of most of them depends on programming and command-line environments, which is a challenge for many wet-lab biologists. Furthermore, current sRNA analysis tools mostly focus on the analysis of certain type of sRNAs and are resource-intensive, normally demanding an immense amount of time and effort to learn the use of numerous tools or scripts and assemble them into a workable pipeline to get the final results. Here, we present sRNAminer, a powerful stand-alone toolkit with a user-friendly interface that integrates all common functions for the analysis of three major types of plant sRNAs: microRNAs (miRNAs), phased small interfering RNAs (phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). We constructed a curated or "golden" set of MIRNA and PHAS loci, which was used to assess the performance of sRNAminer in comparison to other existing tools. The results showed that sRNAminer outperformed these tools in multiple aspects, highlighting its functionality. In addition, to enable an efficient evaluation of sRNA annotation results, we developed Integrative Genomics Viewer (IGV)-sRNA, a modified genome browser optimized from IGV and we incorporated it as a functional module in sRNAminer. IGV-sRNA can display a wealth of sRNA-specific features, enabling a more comprehensive understanding of sRNA data. sRNAminer and IGV-sRNA are both platform-independent software that can be run under all operating systems. They are now freely available at https://github.com/kli28/sRNAminer and https://gitee.com/CJchen/IGV-sRNA.

mRNA cleavage by 21-nucleotide phasiRNAs determines temperature-sensitive male sterility in rice

Temperature-sensitive male sterility is one of the core components for hybrid rice (Oryza sativa) breeding based on the 2-line system. We previously found that knockout of ARGONAUTE 1d (AGO1d) causes temperature-sensitive male sterility in rice by influencing phased small interfering RNA (phasiRNA) biogenesis and function. However, the specific phasiRNAs and their targets underlying the temperature-sensitive male sterility in the ago1d mutant remain unknown. Here, we demonstrate that the ago1d mutant displays normal female fertility but complete male sterility at low temperature. Through a multiomics analysis of small RNA (sRNA), degradome, and transcriptome, we found that 21-nt phasiRNAs account for the greatest proportion of the 21-nt sRNA species in rice anthers and are sensitive to low temperature and markedly downregulated in the ago1d mutant. Moreover, we found that 21-nt phasiRNAs are essential for the mRNA cleavage of a set of fertility- and cold tolerance-associated genes, such as Earlier Degraded Tapetum 1 (EDT1), Tapetum Degeneration Retardation (TDR), OsPCF5, and OsTCP21, directly or indirectly determined by AGO1d-mediated gene silencing. The loss of function of 21-nt phasiRNAs can result in upregulation of their targets and causes varying degrees of defects in male fertility and grain setting. Our results highlight the essential functions of 21-nt phasiRNAs in temperature-sensitive male sterility in rice and suggest their promising application in 2-line hybrid rice breeding in the future.

Assessing small RNA profiles in potato diploid hybrid and its resynthesized allopolyploid reveals conserved abundance with distinct genomic distribution

KEY MESSAGE

The shock produced by the allopolyploidization process on a potato interspecific diploid hybrid displays a non-random remobilization of the small RNAs profile on a variety of genomic features. Allopolyploidy, a complex process involving interspecific hybridization and whole genome duplication, significantly impacts plant evolution, leading to the emergence of novel phenotypes. Polyploids often present phenotypic nuances that enhance adaptability, enabling them to compete better and occasionally to colonize new habitats. Whole-genome duplication represents a genomic "shock" that can trigger genetic and epigenetic changes that yield novel expression patterns. In this work, we investigate the polyploidization effect on a diploid interspecific hybrid obtained through the cross between the cultivated potato Solanum tuberosum and the wild potato Solanum kurtzianum, by assessing the small RNAs (sRNAs) profile of the parental diploid hybrid and its derived allopolyploid. Small RNAs are key components of the epigenetic mechanisms involved in silencing by RNA-directed DNA Methylation (RdDM). A sRNA sequencing (sRNA-Seq) analysis was performed to individually profile the 21 to 22 nucleotide (21 to 22-nt) and 24-nt sRNA size classes due to their unique mechanism of biogenesis and mode of function. The composition and distribution of different genomic features and differentially accumulated (DA) sRNAs were evaluated throughout the potato genome. We selected a subset of genes associated with DA sRNAs for messenger RNA (mRNA) expression analysis to assess potential impacts on the transcriptome. Interestingly, we noted that 24-nt DA sRNAs that exclusively mapped to exons were correlated with differentially expressed mRNAs between genotypes, while this behavior was not observed when 24-nt DA sRNAs were mapped to intronic regions. These findings collectively emphasize the nonstochastic nature of sRNA remobilization in response to the genomic shock induced by allopolyploidization.

Epigenetic weapons of plants against fungal pathogens

In the natural environment, plants face constant exposure to biotic stress caused by fungal attacks. The plant's response to various biotic stresses relies heavily on its ability to rapidly adjust the transcriptome. External signals are transmitted to the nucleus, leading to activation of transcription factors that subsequently enhance the expression of specific defense-related genes. Epigenetic mechanisms, including histone modifications and DNA methylation, which are closely linked to chromatin states, regulate gene expression associated with defense against biotic stress. Additionally, chromatin remodelers and non-coding RNA play a significant role in plant defense against stressors. These molecular modifications enable plants to exhibit enhanced resistance and productivity under diverse environmental conditions. Epigenetic mechanisms also contribute to stress-induced environmental epigenetic memory and priming in plants, enabling them to recall past molecular experiences and utilize this stored information for adaptation to new conditions. In the arms race between fungi and plants, a significant aspect is the cross-kingdom RNAi mechanism, whereby sRNAs can traverse organismal boundaries. Fungi utilize sRNA as an effector molecule to silence plant resistance genes, while plants transport sRNA, primarily through extracellular vesicles, to pathogens in order to suppress virulence-related genes. In this review, we summarize contemporary knowledge on epigenetic mechanisms of plant defense against attack by pathogenic fungi. The role of epigenetic mechanisms during plant-fungus symbiotic interactions is also considered.

Improved super-resolution ribosome profiling reveals prevalent translation of upstream ORFs and small ORFs in Arabidopsis

A crucial step in functional genomics is identifying actively translated ORFs and linking them to biological functions. The challenge lies in identifying short ORFs, as their identification is greatly influenced by data quality and depth. Here, we improved the coverage of super-resolution Ribo-seq in Arabidopsis (Arabidopsis thaliana), revealing uncharacterized translation events for nuclear, chloroplastic, and mitochondrial genes. Assisted by a transcriptome assembly, we identified 7,751 unconventional translation events, comprising 6,996 upstream ORFs (uORFs) and 209 downstream ORFs on annotated protein-coding genes, as well as 546 ORFs in presumed noncoding RNAs. Proteomic data confirmed the production of stable proteins from some of these unannotated translation events. We present evidence of active translation from primary transcripts of trans-acting small interfering RNAs (TAS1-4) and microRNAs (pri-MIR163 and pri-MIR169) and periodic ribosome stalling supporting cotranslational decay. Additionally, we developed a method for identifying extremely short uORFs, including 370 minimum uORFs (AUG-stop), and 2,921 tiny uORFs (2 to 10 amino acids) and 681 uORFs that overlap with each other. Remarkably, these short uORFs exhibit strong translational repression as do longer uORFs. We also systematically discovered 594 uORFs regulated by alternative splicing, suggesting widespread isoform-specific translational control. Finally, these prevalent uORFs are associated with numerous important pathways. In summary, our improved Arabidopsis translational landscape provides valuable resources to study gene expression regulation.

The P2 nucleic acid binding protein of Sugarcane bacilliform virus is a viral pathogenic factor

Background

Saccharum spp. is the primary source of sugar and plays a significant role in global renewable bioenergy. Sugarcane bacilliform virus (SCBV) is one of the most important viruses infecting sugarcane, causing severe yield losses and quality degradation. It is of great significance to reveal the pathogenesis of SCBV and resistance breeding. However, little is known about the viral virulence factors or RNA silencing suppressors and the molecular mechanism of pathogenesis.

Methods

To systematically investigate the functions of the unknown protein P2 encoded by SCBV ORF2. Phylogenetic analysis was implemented to infer the evolutionary relationship between the P2 of SCBV and other badnaviruses. The precise subcellular localization of P2 was verified in the transient infiltrated Nicotiana benthamiana epidermal mesophyll cells and protoplasts using the Laser scanning confocal microscope (LSCM). The post-transcriptional gene silencing (PTGS) and transcriptional gene silencing (TGS) RNA silencing suppressor activity of P2 was analyzed, respectively. Furthermore, restriction digestion and RT-qPCR assays were conducted to verify the probable mechanism of P2 on repressing DNA methylation. To explore the pathogenicity of P2, a potato virus X-based viral vector was used to heterologously express SCBV P2 and the consequent H2O2 accumulation was detected by the 3,3'-diaminobenzidine (DAB) staining method.

Results

Phylogenetic analysis shows that SCBV has no obvious sequence similarity and low genetic relatedness to Badnavirus and Tungrovirus representatives. LSCM studies show that P2 is localized in both the cytoplasm and nucleus. Moreover, P2 is shown to be a suppressor of PTGS and TGS, which can not only repress ssRNA-induced gene silencing but also disrupt the host RNA-directed DNA methylation (RdDM) pathway. In addition, P2 can trigger an oxidative burst and cause typical hypersensitive-like response (HLR) necrosis in systemic leaves of N. benthamiana when expressed by PVX. Overall, our results laid a foundation for deciphering the molecular mechanism of SCBV pathogenesis and made progress for resistance breeding.

The trans-kingdom communication of noncoding RNAs in plant-environment interactions

As conserved regulatory agents, noncoding RNAs (ncRNAs) have an important impact on many aspects of plant life, including growth, development, and environmental response. Noncoding RNAs can travel through not only plasmodesma and phloem but also intercellular barriers to regulate distinct processes. Increasing evidence shows that the intercellular trans-kingdom transmission of ncRNAs is able to modulate many important interactions between plants and other organisms, such as plant response to pathogen attack, the symbiosis between legume plants and rhizobia and the interactions with parasitic plants. In these interactions, plant ncRNAs are believed to be sorted into extracellular vesicles (EVs) or other nonvesicular vehicles to pass through cell barriers and trigger trans-kingdom RNA interference (RNAi) in recipient cells from different species. There is evidence that the features of extracellular RNAs and associated RNA-binding proteins (RBPs) play a role in defining the RNAs to retain in cell or secrete outside cells. Despite the few reports about RNA secretion pathway in plants, the export of extracellular ncRNAs is orchestrated by a series of pathways in plants. The identification and functional analysis of mobile small RNAs (sRNAs) are attracting increasing attention in recent years. In this review, we discuss recent advances in our understanding of the function, sorting, transport, and regulation of plant extracellular ncRNAs.

Phased secondary small interfering RNAs in Camellia sinensis var. assamica

Phased secondary small interfering RNAs (phasiRNAs) in plants play important roles in regulating genome stability, plant development and stress adaption. Camellia sinensis var. assamica has immense economic, medicinal and cultural significance. However, there are still no studies of phasiRNAs and their putative functions in this valuable plant. We identified 476 and 43 PHAS loci which generated 4290 twenty one nucleotide (nt) and 264 twenty four nt phasiRNAs, respectively. Moreover, the analysis of degradome revealed more than 35000 potential targets for these phasiRNAs. We identified several conserved 21 nt phasiRNA generation pathways in tea plant, including miR390 → TAS3, miR482/miR2118 → NB-LRR, miR393 → F-box, miR828 → MYB/TAS4, and miR7122 → PPR in this study. Furthermore, we found that some transposase and plant mobile domain genes could generate phasiRNAs. Our results show that phasiRNAs target genes in the same family in cis- or trans-manners, and different members of the same gene family may generate the same phasiRNAs. The phasiRNAs, generated by transposase and plant mobile domain genes, and their targets, suggest that phasiRNAs may be involved in the inhibition of transposable elements in tea plant. To summarize, these results provide a comprehensive view of phasiRNAs in Camellia sinensis var. assamica.

The P1 protein of wheat yellow mosaic virus exerts RNA silencing suppression activity to facilitate virus infection in wheat plants

Wheat yellow mosaic virus (WYMV) causes severe wheat viral disease in Asia. However, the viral suppressor of RNA silencing (VSR) encoded by WYMV has not been identified. Here, the P1 protein encoded by WYMV RNA2 was shown to suppress RNA silencing in Nicotiana benthamiana. Mutagenesis assays revealed that the alanine substitution mutant G175A of P1 abolished VSR activity and mutant Y10A VSR activity remained only in younger leaves. P1, but not G175A, interacted with gene silencing-related protein, N. benthamiana calmodulin-like protein (NbCaM), and calmodulin-binding transcription activator 3 (NbCAMTA3), and Y10A interacted with NbCAMTA3 only. Competitive Bimolecular fluorescence complementation and co-immunoprecipitation assays showed that the ability of P1 disturbing the interaction between NbCaM and NbCAMTA3 was stronger than Y10A, Y10A was stronger than G175A. In vitro transcript inoculation of infectious WYMV clones further demonstrated that VSR-defective mutants G175A and Y10A reduced WYMV infection in wheat (Triticum aestivum L.), G175A had a more significant effect on virus accumulation in upper leaves of wheat than Y10A. Moreover, RNA silencing, temperature, and autophagy have significant effects on the accumulation of P1 in N. benthamiana. Taken together, WYMV P1 acts as VSR by interfering with calmodulin-associated antiviral RNAi defense to facilitate virus infection in wheat, which has provided clear insights into the function of P1 in the process of WYMV infection.

Genome-Wide Identification and Characterization of Major RNAi Genes Highlighting Their Associated Factors in Cowpea (Vigna unguiculata (L.) Walp.)

In different regions of the world, cowpea (Vigna unguiculata (L.) Walp.) is an important vegetable and an excellent source of protein. It lessens the malnutrition of the underprivileged in developing nations and has some positive effects on health, such as a reduction in the prevalence of cancer and cardiovascular disease. However, occasionally, certain biotic and abiotic stresses caused a sharp fall in cowpea yield. Major RNA interference (RNAi) genes like Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) are essential for the synthesis of their associated factors like domain, small RNAs (sRNAs), transcription factors, micro-RNAs, and cis-acting factors that shield plants from biotic and abiotic stresses. In this study, applying BLASTP search and phylogenetic tree analysis with reference to the Arabidopsis RNAi (AtRNAi) genes, we discovered 28 VuRNAi genes, including 7 VuDCL, 14 VuAGO, and 7 VuRDR genes in cowpea. We looked at the domains, motifs, gene structures, chromosomal locations, subcellular locations, gene ontology (GO) terms, and regulatory factors (transcription factors, micro-RNAs, and cis-acting elements (CAEs)) to characterize the VuRNAi genes and proteins in cowpea in response to stresses. Predicted VuDCL1, VuDCL2(a, b), VuAGO7, VuAGO10, and VuRDR6 genes might have an impact on cowpea growth, development of the vegetative and flowering stages, and antiviral defense. The VuRNAi gene regulatory features miR395 and miR396 might contribute to grain quality improvement, immunity boosting, and pathogen infection resistance under salinity and drought conditions. Predicted CAEs from the VuRNAi genes might play a role in plant growth and development, improving grain quality and production and protecting plants from biotic and abiotic stresses. Therefore, our study provides crucial information about the functional roles of VuRNAi genes and their associated components, which would aid in the development of future cowpeas that are more resilient to biotic and abiotic stress. The manuscript is available as a preprint at this link: doi:10.1101/2023.02.15.528631v1.

The Pivotal Role of Noncoding RNAs in Flowering Time Regulation

Noncoding RNAs constitute a substantial portion of the transcriptome and play pivotal roles in plant growth and development. Among these processes, flowering stands out as a crucial trait, ensuring reproductive success and seed set, and is meticulously controlled by genetic and environmental factors. With remarkable advancements in the identification and characterization of noncoding RNAs in plants, it has become evident that noncoding RNAs are intricately linked to the regulation of flowering time. In this article, we present an overview of the classification of plant noncoding RNAs and delve into their functions in the regulation of flowering time. Furthermore, we review their molecular mechanisms and their involvement in flowering pathways. Our comprehensive review enhances the understanding of how noncoding RNAs contribute to the regulation of flowering time and sheds light on their potential implications in crop breeding.

The phased short-interfering RNA siRD29(-) regulates GIBBERELLIN 3-OXIDASE 3 during stolon-to-tuber transitions in potato

Phased short-interfering RNAs (phasiRNAs) fine tune various stages of growth, development, and stress responses in plants. Potato (Solanum tuberosum) tuberization is a complex process, wherein a belowground modified stem (stolon) passes through developmental stages like swollen stolon and minituber before it matures to a potato. Previously, we identified several phasiRNA-producing loci (PHAS) from stolon-to-tuber transition stages. However, whether phasiRNAs mediate tuber development remains unknown. Here, we show that a gene encoding NB-ARC DOMAIN-CONTAINING DISEASE RESISTANCE PROTEIN (StRGA4; a PHAS locus) is targeted by Stu-microRNA482c to generate phasiRNAs. Interestingly, we observed that one of the phasiRNAs, referred as short-interfering RNA D29(-), i.e. siRD29(-), targets the gibberellin (GA) biosynthesis gene GIBBERELLIN 3-OXIDASE 3 (StGA3ox3). Since regulation of bioactive GA levels in stolons controls tuber development, we hypothesized that a gene regulatory module, Stu-miR482c-StRGA4-siRD29(-)-StGA3ox3, could govern tuber development. Through transient expression assays and small RNA sequencing, generation of siRD29(-) and its phase was confirmed in planta. Notably, the expression of StGA3ox3 was higher in swollen stolon compared to stolon, whereas siRD29(-) showed a negative association with StGA3ox3 expression. Antisense (AS) lines of StGA3ox3 produced more tubers compared to wild type. As expected, StRGA4 overexpression (OE) lines had high levels of siRD29(-) and mimicked the phenotypes of StGA3ox3-AS lines, indicating the functionality of this module in potato. In vitro tuberization assays (with or without a GA inhibitor) using StGA3ox3 antisense lines and overexpression lines of StGA3ox3 or StRGA4 revealed that StGA3ox3 controls the tuber stalk development. Taken together, our findings suggest that a phasiRNA, siRD29(-), mediates the regulation of StGA3ox3 during stolon-to-tuber transitions in potato.

Multiomics analysis reveals a link between Brassica-specific miR1885 and rapeseed tolerance to low temperature

Brassica crops include various edible vegetable and plant oil crops, and their production is limited by low temperature beyond their tolerant capability. The key regulators of low-temperature resistance in Brassica remain largely unexplored. To identify posttranscriptional regulators of plant response to low temperature, we performed small RNA profiling, and found that 16 known miRNAs responded to cold treatment in Brassica rapa. The cold response of seven of those miRNAs were further confirmed by qRT-PCR and/or northern blot analyses. In parallel, a genome-wide association study of 220 accessions of Brassica napus identified four candidate MIRNA genes, all of which were cold-responsive, at the loci associated with low-temperature resistance. Specifically, these large-scale data analyses revealed a link between miR1885 and the plant response to low temperature in both B. rapa and B. napus. Using 5' rapid amplification of cDNA ends approach, we validated that miR1885 can cleave its putative target gene transcripts, Bn.TIR.A09 and Bn.TNL.A03, in B. napus. Furthermore, overexpression of miR1885 in Semiwinter type B. napus decreased the mRNA abundance of Bn.TIR.A09 and Bn.TNL.A03 and resulted in increased sensitivity to low temperature. Knocking down of miR1885 in Spring type B. napus led to increased mRNA abundance of its targets and improved rapeseed tolerance to low temperature. Together, our results suggested that the loci of miR1885 and its targets could be potential candidates for the molecular breeding of low temperature-tolerant Spring type Brassica crops.

Regulation of plant immunity via small RNA-mediated control of NLR expression

Plants use different receptors to detect potential pathogens: membrane-anchored pattern recognition receptors (PRRs) activated upon perception of pathogen-associated molecular patterns (PAMPs) that elicit pattern-triggered immunity (PTI); and intracellular nucleotide-binding leucine-rich repeat proteins (NLRs) activated by detection of pathogen-derived effectors, activating effector-triggered immunity (ETI). The interconnections between PTI and ETI responses have been increasingly reported. Elevated NLR levels may cause autoimmunity, with symptoms ranging from fitness cost to developmental arrest, sometimes combined with run-away cell death, making accurate control of NLR dosage key for plant survival. Small RNA-mediated gene regulation has emerged as a major mechanism of control of NLR dosage. Twenty-two nucleotide miRNAs with the unique ability to trigger secondary siRNA production from target transcripts are particularly prevalent in NLR regulation. They enhance repression of the primary NLR target, but also bring about repression of NLRs only complementary to secondary siRNAs. We summarize current knowledge on miRNAs and siRNAs in the regulation of NLR expression with an emphasis on 22 nt miRNAs and propose that miRNA and siRNA regulation of NLR levels provides additional links between PTI and NLR defense pathways to increase plant responsiveness against a broad spectrum of pathogens and control an efficient deployment of defenses.

Genome-wide identification of Sclerotinia sclerotiorum small RNAs and their endogenous targets

BACKGROUND

Several phytopathogens produce small non-coding RNAs of approximately 18-30 nucleotides (nt) which post-transcriptionally regulate gene expression. Commonly called small RNAs (sRNAs), these small molecules were also reported to be present in the necrotrophic pathogen Sclerotinia sclerotiorum. S. sclerotiorum causes diseases in more than 400 plant species, including the important oilseed crop Brassica napus. sRNAs can further be classified as microRNAs (miRNAs) and short interfering RNAs (siRNAs). Certain miRNAs can activate loci that produce further sRNAs; these secondary sRNA-producing loci are called 'phased siRNA' (PHAS) loci and have only been described in plants. To date, very few studies have characterized sRNAs and their endogenous targets in S. sclerotiorum.

RESULTS

We used Illumina sequencing to characterize sRNAs from fungal mycelial mats of S. sclerotiorum spread over B. napus leaves. In total, eight sRNA libraries were prepared from in vitro, 12 h post-inoculation (HPI), and 24 HPI mycelial mat samples. Cluster analysis identified 354 abundant sRNA clusters with reads of more than 100 Reads Per Million (RPM). Differential expression analysis revealed upregulation of 34 and 57 loci at 12 and 24 HPI, respectively, in comparison to in vitro samples. Among these, 25 loci were commonly upregulated. Altogether, 343 endogenous targets were identified from the major RNAs of 25 loci. Almost 88% of these targets were annotated as repeat element genes, while the remaining targets were non-repeat element genes. Fungal degradome reads confirmed cleavage of two transposable elements by one upregulated sRNA. Altogether, 24 milRNA loci were predicted with both mature and milRNA* (star) sequences; these are both criteria associated previously with experimentally verified miRNAs. Degradome sequencing data confirmed the cleavage of 14 targets. These targets were related to repeat element genes, phosphate acetyltransferases, RNA-binding factor, and exchange factor. A PHAS gene prediction tool identified 26 possible phased interfering loci with 147 phasiRNAs from the S. sclerotiorum genome, suggesting this pathogen might produce sRNAs that function similarly to miRNAs in higher eukaryotes.

CONCLUSIONS

Our results provide new insights into sRNA populations and add a new resource for the study of sRNAs in S. sclerotiorum.

Plastid dsRNA transgenes trigger phased small RNA-based gene silencing of nuclear-encoded genes

Plastid transformation technology has been widely used to express traits of potential commercial importance, though the technology has been limited to traits that function while sequestered in the organelle. Prior research indicates that plastid contents can escape from the organelle, suggesting a possible mechanism for engineering plastid transgenes to function in other cellular locations. To test this hypothesis, we created tobacco (Nicotiana tabacum cv. Petit Havana) plastid transformants that express a fragment of the nuclear-encoded Phytoene desaturase (PDS) gene capable of catalyzing post-transcriptional gene silencing if RNA escapes into the cytoplasm. We found multiple lines of direct evidence that plastid-encoded PDS transgenes affect nuclear PDS gene silencing: knockdown of the nuclear-encoded PDS mRNA and/or its apparent translational inhibition, biogenesis of 21-nucleotide (nt) phased small interfering RNAs (phasiRNAs), and pigment-deficient plants. Furthermore, plastid-expressed dsRNA with no cognate nuclear-encoded pairing partner also produced abundant 21-nt phasiRNAs in the cytoplasm, demonstrating that a nuclear-encoded template is not required for siRNA biogenesis. Our results indicate that RNA escape from plastids to the cytoplasm occurs generally, with functional consequences that include entry into the gene silencing pathway. Furthermore, we uncover a method to produce plastid-encoded traits with functions outside of the organelle and open additional fields of study in plastid development, compartmentalization, and small RNA biogenesis.

Revitalizing miRNAs mediated agronomical advantageous traits improvement in rice

One of the key enigmas in conventional and modern crop improvement programmes is how to introduce beneficial traits without any penalty impairment. Rice (Oryza sativa L.), among the essential staple food crops grown and utilized worldwide, needs to improve genotypes in multifaceted ways. With the global view to feed ten billion under the climatic perturbation, only a potent functional master regulator can withstand with hope for the next green revolution and food security. miRNAs are such, miniature, fine tuners for crop improvement and provide a value addition in emerging technologies, namely large-scale genotyping, phenotyping, genome editing, marker-assisted selection, and genomic selection, to make rice production feasible. There has been surplus research output generated since the last decade on miRNAs in rice, however, recent functional knowledge is limited to reaping the benefits for conventional and modern improvements in rice to avoid ambiguity and redundancy in the generated data. Here, we present the latest functional understanding of miRNAs in rice. In addition, their biogenesis, intra- and inter-kingdom signaling and communication, implication of amiRNAs, and consequences upon integration with CRISPR-Cas9. Further, highlights refer to the application of miRNAs for rice agronomical trait improvements, broadly classified into three functional domains. The majority of functionally established miRNAs are responsible for growth and development, followed by biotic and abiotic stresses. Tabular cataloguing reveals and highlights two multifaceted modules that were extensively studied. These belong to miRNA families 156 and 396, orchestrate multifarious aspects of advantageous agronomical traits. Moreover, updated and exhaustive functional aspects of different supplemental miRNA modules that would strengthen rice improvement are also being discussed.

Noncoding RNAs and their roles in regulating the agronomic traits of crops

Molecular breeding is one of the most effective methods for improving the performance of crops. Understanding the genome features of crops, especially the physiological functions of individual genes, is of great importance to molecular breeding. Evidence has shown that genomes of both animals and plants transcribe numerous non-coding RNAs, which are involved in almost every aspect of development. In crops, an increasing number of studies have proven that non-coding RNAs are new genetic resources for regulating crop traits. In this review, we summarize the current knowledge of non-coding RNAs, which are potential crop trait regulators, and focus on the functions of long non-coding RNAs (lncRNAs) in determining crop grain yield, phased small-interfering RNAs (phasiRNAs) in regulating fertility, small interfering RNAs (siRNAs) and microRNAs (miRNAs) in facilitating plant immune response and disease resistance, and miRNAs mediating nutrient and metal stress. Finally, we also discuss the next-generation method for ncRNA application in crop domestication and breeding.

The role of microRNAs in responses to drought and heat stress in peanut (Arachis hypogaea)

MicroRNAs (miRNAs) are 21-24 nt small RNAs (sRNAs) that negatively regulate protein-coding genes and/or trigger phased small-interfering RNA (phasiRNA) production. Two thousand nine hundred miRNA families, of which ∼40 are deeply conserved, have been identified in ∼80 different plant species genomes. miRNA functions in response to abiotic stresses is less understood than their roles in development. Only seven peanut MIRNA families are documented in miRBase, yet a reference genome assembly is now published and over 480 plant-like MIRNA loci were predicted in the diploid peanut progenitor Arachis duranensis genome. We explored by computational analysis of a leaf sRNA library and publicly available sRNA, degradome, and transcriptome datasets the miRNA and phasiRNA space associated with drought and heat stresses in peanut. We characterized 33 novel candidate and 33 ancient conserved families of MIRNAs and present degradome evidence for their cleavage activities on mRNA targets, including several noncanonical targets and novel phasiRNA-producing noncoding and mRNA loci with validated novel targets such as miR1509 targeting serine/threonine-protein phosphatase7 and miRc20 and ahy-miR3514 targeting penta-tricopeptide repeats (PPRs), in contradistinction to other claims of miR1509/173/7122 superfamily miRNAs indirectly targeting PPRs via TAS-like noncoding RNA loci. We characterized the inverse correlations of significantly differentially expressed drought- and heat-regulated miRNAs, assayed by sRNA blots or transcriptome datasets, with target mRNA expressions in the same datasets. Meta-analysis of an expression atlas and over representation of miRNA target genes in co-expression networks suggest that miRNAs have functions in unique aspects of peanut gynophore development. Genome-wide MIRNA annotation of the published allopolyploid peanut genome can facilitate molecular breeding of value-added traits.

RNA-Based Control of Fungal Pathogens in Plants

Our duty to conserve global natural ecosystems is increasingly in conflict with our need to feed an expanding population. The use of conventional pesticides not only damages the environment and vulnerable biodiversity but can also still fail to prevent crop losses of 20-40% due to pests and pathogens. There is a growing call for more ecologically sustainable pathogen control measures. RNA-based biopesticides offer an eco-friendly alternative to the use of conventional fungicides for crop protection. The genetic modification (GM) of crops remains controversial in many countries, though expression of transgenes inducing pathogen-specific RNA interference (RNAi) has been proven effective against many agronomically important fungal pathogens. The topical application of pathogen-specific RNAi-inducing sprays is a more responsive, GM-free approach to conventional RNAi transgene-based crop protection. The specific targeting of essential pathogen genes, the development of RNAi-nanoparticle carrier spray formulations, and the possible structural modifications to the RNA molecules themselves are crucial to the success of this novel technology. Here, we outline the current understanding of gene silencing pathways in plants and fungi and summarize the pioneering and recent work exploring RNA-based biopesticides for crop protection against fungal pathogens, with a focus on spray-induced gene silencing (SIGS). Further, we discuss factors that could affect the success of RNA-based control strategies, including RNA uptake, stability, amplification, and movement within and between the plant host and pathogen, as well as the cost and design of RNA pesticides.

The OsSGS3-tasiRNA-OsARF3 module orchestrates abiotic-biotic stress response trade-off in rice

Recurrent heat stress and pathogen invasion seriously threaten crop production, and abiotic stress often antagonizes biotic stress response against pathogens. However, the molecular mechanisms of trade-offs between thermotolerance and defense remain obscure. Here, we identify a rice thermo-sensitive mutant that displays a defect in floret development under high temperature with a mutation in SUPPRESSOR OF GENE SILENCING 3a (OsSGS3a). OsSGS3a interacts with its homolog OsSGS3b and modulates the biogenesis of trans-acting small interfering RNA (tasiRNA) targeting AUXIN RESPONSE FACTORS (ARFs). We find that OsSGS3a/b positively, while OsARF3a/b and OsARF3la/lb negatively modulate thermotolerance. Moreover, OsSGS3a negatively, while OsARF3a/b and OsARF3la/lb positively regulate disease resistance to the bacterial pathogen Xanthomonas oryzae pv. oryzae (Xoo) and the fungal pathogen Magnaporthe oryzae (M. oryzae). Taken together, our study uncovers a previously unknown trade-off mechanism that regulates distinct immunity and thermotolerance through the OsSGS3-tasiRNA-OsARF3 module, highlighting the regulation of abiotic-biotic stress response trade-off in plants.

Genome-wide identification of DCL, AGO and RDR gene families and their associated functional regulatory element analyses in sunflower (Helianthus annuus)

RNA interference (RNAi) regulates a variety of eukaryotic gene expressions that are engaged in response to stress, growth, and the conservation of genomic stability during developmental phases. It is also intimately connected to the post-transcriptional gene silencing (PTGS) process and chromatin modification levels. The entire process of RNA interference (RNAi) pathway gene families mediates RNA silencing. The main factors of RNA silencing are the Dicer-Like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) gene families. To the best of our knowledge, genome-wide identification of RNAi gene families like DCL, AGO, and RDR in sunflower (Helianthus annuus) has not yet been studied despite being discovered in some species. So, the goal of this study is to find the RNAi gene families like DCL, AGO, and RDR in sunflower based on bioinformatics approaches. Therefore, we accomplished an inclusive in silico investigation for genome-wide identification of RNAi pathway gene families DCL, AGO, and RDR through bioinformatics approaches such as (sequence homogeneity, phylogenetic relationship, gene structure, chromosomal localization, PPIs, GO, sub-cellular localization). In this study, we have identified five DCL (HaDCLs), fifteen AGO (HaAGOs), and ten RDR (HaRDRs) in the sunflower genome database corresponding to the RNAi genes of model plant Arabidopsis thaliana based on genome-wide analysis and a phylogenetic method. The analysis of the gene structure that contains exon-intron numbers, conserved domain, and motif composition analyses for all HaDCL, HaAGO, and HaRDR gene families indicated almost homogeneity among the same gene family. The protein-protein interaction (PPI) network analysis illustrated that there exists interconnection among identified three gene families. The analysis of the Gene Ontology (GO) enrichment showed that the detected genes directly contribute to the RNA gene-silencing and were involved in crucial pathways. It was observed that the cis-acting regulatory components connected to the identified genes were shown to be responsive to hormone, light, stress, and other functions. That was found in HaDCL, HaAGO, and HaRDR genes associated with the development and growth of plants. Finally, we are able to provide some essential information about the components of sunflower RNA silencing through our genome-wide comparison and integrated bioinformatics analysis, which open the door for further research into the functional mechanisms of the identified genes and their regulatory elements.

Large Artificial microRNA Cluster Genes Confer Effective Resistance against Multiple Tomato Yellow Leaf Curl Viruses in Transgenic Tomato

Tomato yellow leaf curl disease (TYLCD) has become the key limiting factor for the production of tomato in many areas because of the continuous infection and recombination of several tomato yellow leaf curl virus (TYLCV)-like species (TYLCLV) which produce novel and destructive viruses. Artificial microRNA (AMIR) is a recent and effective technology used to create viral resistance in major crops. This study applies AMIR technology in two ways, i.e., amiRNA in introns (AMINs) and amiRNA in exons (AMIEs), to express 14 amiRNAs targeting conserved regions in seven TYLCLV genes and their satellite DNA. The resulting pAMIN14 and pAMIE14 vectors can encode large AMIR clusters and their function in silencing reporter genes was validated with transient assays and stable transgenic N. tabacum plants. To assess the efficacy of conferring resistance against TYLCLV, pAMIE14 and pAMIN14 were transformed into tomato cultivar A57 and the resulting transgenic tomato plants were evaluated for their level of resistance to mixed TYLCLV infection. The results suggest that pAMIN14 transgenic lines have a more effective resistance than pAMIE14 transgenic lines, reaching a resistance level comparable to plants carrying the TY1 resistance gene.

Phylogenetic analyses of seven protein families refine the evolution of small RNA pathways in green plants

Several protein families participate in the biogenesis and function of small RNAs (sRNAs) in plants. Those with primary roles include Dicer-like (DCL), RNA-dependent RNA polymerase (RDR), and Argonaute (AGO) proteins. Protein families such as double-stranded RNA-binding (DRB), SERRATE (SE), and SUPPRESSION OF SILENCING 3 (SGS3) act as partners of DCL or RDR proteins. Here, we present curated annotations and phylogenetic analyses of seven sRNA pathway protein families performed on 196 species in the Viridiplantae (aka green plants) lineage. Our results suggest that the RDR3 proteins emerged earlier than RDR1/2/6. RDR6 is found in filamentous green algae and all land plants, suggesting that the evolution of RDR6 proteins coincides with the evolution of phased small interfering RNAs (siRNAs). We traced the origin of the 24-nt reproductive phased siRNA-associated DCL5 protein back to the American sweet flag (Acorus americanus), the earliest diverged, extant monocot species. Our analyses of AGOs identified multiple duplication events of AGO genes that were lost, retained, or further duplicated in subgroups, indicating that the evolution of AGOs is complex in monocots. The results also refine the evolution of several clades of AGO proteins, such as AGO4, AGO6, AGO17, and AGO18. Analyses of nuclear localization signal sequences and catalytic triads of AGO proteins shed light on the regulatory roles of diverse AGOs. Collectively, this work generates a curated and evolutionarily coherent annotation for gene families involved in plant sRNA biogenesis/function and provides insights into the evolution of major sRNA pathways.

Trans-species microRNA loci in the parasitic plant Cuscuta campestris have a U6-like snRNA promoter

Small regulatory RNAs can move between organisms and regulate gene expression in the recipient. Whether the trans-species small RNAs being exported are distinguished from the normal endogenous small RNAs of the source organism is not known. The parasitic plant Cuscuta campestris (dodder) produces many microRNAs that specifically accumulate at the host-parasite interface, several of which have trans-species activity. We found that induction of C. campestris interface-induced microRNAs is similar regardless of host species and occurs in C. campestris haustoria produced in the absence of any host. The loci-encoding C. campestris interface-induced microRNAs are distinguished by a common cis-regulatory element. This element is identical to a conserved upstream sequence element (USE) used by plant small nuclear RNA loci. The properties of the interface-induced microRNA primary transcripts strongly suggest that they are produced via U6-like transcription by RNA polymerase III. The USE promotes accumulation of interface-induced miRNAs (IIMs) in a heterologous system. This promoter element distinguishes C. campestris IIM loci from other plant small RNAs. Our data suggest that C. campestris IIMs are produced in a manner distinct from canonical miRNAs. All confirmed C. campestris microRNAs with documented trans-species activity are interface-induced and possess these features. We speculate that RNA polymerase III transcription of IIMs may allow these miRNAs to be exported to hosts.

Indian Himalayan natural Arabidopsis thaliana accessions with abolished miR158 levels exhibit robust miR173-initiated trans-acting cascade silencing

Small RNAs (sRNAs) such as microRNAs (miRNAs) and small interfering RNAs (siRNAs) are short 20-24-nucleotide non-coding RNAs. They are key regulators of gene expression in plants and other organisms. Several 22-nucleotide miRNAs trigger biogenesis cascades of trans-acting secondary siRNAs, which are involved in various developmental and stress responses. Here we show that Himalayan Arabidopsis thaliana accessions having natural mutations in the miR158 locus exhibit robust cascade silencing of the pentatricopeptide repeat (PPR)-like locus. Furthermore, we show that these cascade sRNAs trigger tertiary silencing of a gene involved in transpiration and stomatal opening. The natural deletions or insertions in MIR158 led to improper processing of miR158 precursors, thereby blocking synthesis of mature miR158. Reduced miR158 levels led to increased levels of its target, a pseudo-PPR gene that is targeted by tasiRNAs generated by the miR173 cascade in other accessions. Using sRNA datasets derived from Indian Himalayan accessions, as well as overexpression and knockout lines of miR158, we show that absence of miR158 led to buildup of pseudo-PPR-derived tertiary sRNAs. These tertiary sRNAs mediated robust silencing of a gene involved in stomatal closure in Himalayan accessions lacking miR158 expression. We functionally validated the tertiary phasiRNA that targets NHX2, which encodes a Na⁺ -K⁺ /H⁺ antiporter protein, thereby regulating transpiration and stomatal conductance. Overall, we report the role of the miRNA-TAS-siRNA-pseudogene-tertiary phasiRNA-NHX2 pathway in plant adaptation.

Strategies for Engineering Virus Resistance in Potato

Potato (Solanum tuberosum L.) is an important vegetable crop that plays a pivotal role in the world, especially given its potential to feed the world population and to act as the major staple food in many developing countries. Every year, significant crop loss is caused by viral diseases due to a lack of effective agrochemical treatments, since only transmission by insect vectors can be combated with the use of insecticides, and this has been an important factor hindering potato production. With the rapid development of molecular biology and plant genetic engineering technology, transgenic approaches and non-transgenic techniques (RNA interference and CRISPR-cas9) have been effectively employed to improve potato protection against devastating viruses. Moreover, the availability of viral sequences, potato genome sequences, and host immune mechanisms has remarkably facilitated potato genetic engineering. In this study, we summarize the progress of antiviral strategies applied in potato through engineering either virus-derived or plant-derived genes. These recent molecular insights into engineering approaches provide the necessary framework to develop viral resistance in potato in order to provide durable and broad-spectrum protection against important viral diseases of solanaceous crops.

Plant non-coding RNAs function in pollen development and male sterility

Male sterility is classified as either cytoplasmic male sterility (CMS) or genic male sterility (GMS). Generally, CMS involves mitochondrial genomes interacting with the nuclear genome, while GMS is caused by nuclear genes alone. Male sterility is regulated by multilevel mechanisms in which non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and phased small interfering RNAs (phasiRNAs), which have been proven to be critical elements. The development of high-throughput sequencing technology offers new opportunities to evaluate the genetic mechanism of ncRNAs in plant male sterility. In this review, we summarize the critical ncRNAs that regulate gene expression in ways dependent on or independent of hormones, which involve the differentiation of the stamen primordia, degradation of the tapetum, formation of microspores, and the release of pollen. In addition, the key mechanisms of the miRNA-lncRNA-mRNA interaction networks mediating male sterility in plants are elaborated. We present a different perspective on exploring the ncRNA-mediated regulatory pathways that control CMS in plants and create male-sterile lines through hormones or genome editing. A refined understanding of the ncRNA regulatory mechanisms in plant male sterility for the development of new sterile lines would be conducive to improve hybridization breeding.

Taxon-specific, phased siRNAs underlie a speciation locus in monkeyflowers

Taxon-specific small RNA loci are widespread in eukaryotic genomes, yet their role in lineage-specific adaptation, phenotypic diversification, and speciation is poorly understood. Here, we report that a speciation locus in monkeyflowers (Mimulus), YELLOW UPPER (YUP), contains an inverted repeat region that produces small interfering RNAs (siRNAs) in a phased pattern. Although the inverted repeat is derived from a partial duplication of a protein-coding gene that is not involved in flower pigmentation, one of the siRNAs targets and represses a master regulator of floral carotenoid pigmentation. YUP emerged with two protein-coding genes that control other aspects of flower coloration as a "superlocus" in a subclade of Mimulus and has contributed to subsequent phenotypic diversification and pollinator-mediated speciation in the descendant species.

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.

Tweaking the Small Non-Coding RNAs to Improve Desirable Traits in Plant

Plant transcriptome contains an enormous amount of non-coding RNAs (ncRNAs) that do not code for proteins but take part in regulating gene expression. Since their discovery in the early 1990s, much research has been conducted to elucidate their function in the gene regulatory network and their involvement in plants' response to biotic/abiotic stresses. Typically, 20-30 nucleotide-long small ncRNAs are a potential target for plant molecular breeders because of their agricultural importance. This review summarizes the current understanding of three major classes of small ncRNAs: short-interfering RNAs (siRNAs), microRNA (miRNA), and transacting siRNAs (tasiRNAs). Furthermore, their biogenesis, mode of action, and how they have been utilized to improve crop productivity and disease resistance are discussed here.

Temperature-sensitive male sterility in rice determined by the roles of AGO1d in reproductive phasiRNA biogenesis and function

Phased secondary siRNAs (phasiRNAs) are broadly present in the reproductive tissues of flowering plants, with spatial-temporal specificity. However, the ARGONAUTE (AGO) proteins associated with phasiRNAs and their miRNA triggers remain elusive. Here, through histological and high-throughput sequencing analyses, we show that rice AGO1d, which is specifically expressed in anther wall cells before and during meiosis, associates with both miR2118 and miR2275 to mediate phasiRNA biogenesis. AGO1d preferentially binds to miR2118-triggered 21-nucleotide (nt) phasiRNAs with a 5'-terminal uridine, suggesting a dual role in phasiRNA biogenesis and function. Depletion of AGO1d causes a reduction of 21- and 24-nt phasiRNAs and temperature-sensitive male sterility. At lower temperatures, anthers of the ago1d mutant predominantly show excessive tapetal cells with little starch accumulation during pollen formation, possibly caused by the dysregulation of cell metabolism. These results uncover an essential role of AGO1d in rice anther development at lower temperatures and demonstrate coordinative roles of AGO proteins during reproductive phasiRNA biogenesis and function.