Coordination of MicroRNAs, PhasiRNAs, and NB-LRR Genes in Response to a Plant Pathogen: Insights from Analyses of a Set of Soybean Rps Gene Near-Isogenic Lines

Transposon-Associated Small RNAs Involved in Plant Defense in Poplar

Utilizing high-throughput Illumina sequencing, we examined how small RNA (sRNA) profiles vary in Chinese white poplar (Populus tomentosa) across two pivotal infection stages by the rust fungus Melampsora larici-populina: the biotrophic growth phase (T02; 48 h post infection) and the urediniospore development and dispersal phase (T03; 168 h), both essential for plant colonization and prolonged biotrophic engagement. Far exceeding random expectations, siRNA clusters predominantly arose from transposon regions, with pseudogenes also contributing significantly, and infection-stage-specific variations were notably tied to these transposon-derived siRNAs. As the infection advanced, clusters of 24 nt siRNAs in transposon and intergenic regions exhibited pronounced abundance shifts. An analysis of targets indicated that Populus sRNAs potentially regulate 95% of Melampsora larici-populina genes, with pathogen effector genes showing heightened targeting by sRNAs during the biotrophic and urediniospore phases compared to controls, pointing to selective sRNA-target interactions. In contrast to conserved miRNAs across plant species, Populus-specific miRNAs displayed a markedly greater tendency to target NB-LRR genes. These observations collectively highlight the innovative roles of sRNAs in plant defense, their evolutionary roots, and their dynamic interplay with pathogen coevolution.

RNAi-based biocontrol for crops: a revised expectation for a non-recent technology

MAIN CONCLUSION

The exogenous application of RNAi technology offers new promises for crops improvement. Cell-based or synthetically produced strands are economical, non-transgenic and could induce the same responses. The substantial population growth demands novel strategies to produce crops without further damaging the environment. RNA interference mechanism is one of the promising technologies to biologically control pests and pathogens in crops, suppressing them by cancelling protein synthesis related to parasitism/pathogenesis. The transgenic approach to generate host-induced gene silencing demonstrated high efficacy in controlling pests or pathogens by RNAi mechanism. However, transgenic technology is tightly regulated and still negatively perceived by global consumers. This review presents the basic biology of small RNA, the main actor of the RNAi mechanism, and tested non-transgenic approaches to induce RNAi exogenously. Novel avenues are offered by the discovery of cross-kingdom RNAi, that naturally, plants also deliver small RNA to suppress the growth of their threats. Future applications of non-transgenic RNAi-based biocontrol will involve the production of dsRNA on an industrial scale. Here, the attempts to provide dsRNA for routine application in farms are also discussed, emphasizing that the technology must be accessible by the countries with the greatest population which mostly are poorer ones.

Exogenous dsRNA-Mediated RNAi: Mechanisms, Applications, Delivery Methods and Challenges in the Induction of Viral Disease Resistance in Plants

The increasing challenges posed by plant viral diseases demand innovative and sustainable management strategies to minimize agricultural losses. Exogenous double-stranded RNA (dsRNA)-mediated RNA interference (RNAi) represents a transformative approach to combat plant viral pathogens without the need for genetic transformation. This review explores the mechanisms underlying dsRNA-induced RNAi, highlighting its ability to silence specific viral genes through small interfering RNAs (siRNAs). Key advancements in dsRNA production, including cost-effective microbial synthesis and in vitro methods, are examined alongside delivery techniques such as spray-induced gene silencing (SIGS) and nanocarrier-based systems. Strategies for enhancing dsRNA stability, including the use of nanomaterials like layered double hydroxide nanosheets and carbon dots, are discussed to address environmental degradation challenges. Practical applications of this technology against various plant viruses and its potential to ensure food security are emphasized. The review also delves into regulatory considerations, risk assessments, and the challenges associated with off-target effects and pathogen resistance. By evaluating both opportunities and limitations, this review underscores the role of exogenous dsRNA as a sustainable solution for achieving viral disease resistance 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.

The biogenesis, regulation and functions of transitive siRNA in plants

Small RNA (sRNA)-mediated RNA interference (RNAi) is a sequence-specific gene silencing mechanism that modulates gene expression in eukaryotes. As core molecules of RNAi, various sRNAs are encoded in the plant genome or derived from invading RNA molecules, and their biogenesis depends on distinct genetic pathways. Transitive small interfering RNAs (siRNAs), which are sRNAs produced from double-strand RNA (dsRNA) in a process that depends on RNA-dependent RNA polymerases (RDRs), can amplify and spread silencing signals to additional transcripts, thereby enabling a phenomenon termed "transitive RNAi". Members of this class of siRNAs function in various biological processes ranging from development to stress adaptation. In Arabidopsis thaliana, two RDRs participate in the generation of transitive siRNAs, acting cooperatively with various siRNA generation-related factors, such as the RNA-induced silencing complex (RISC) and aberrant RNAs. Transitive siRNAs are produced in diverse subcellular locations and structures under the control of various mechanisms, highlighting the intricacies of their biogenesis and functions. In this review, we discuss recent advances in understanding the molecular events of transitive siRNA biogenesis and its regulation, with a particular focus on factors involved in RDR recruitment. We aim to provide a comprehensive description of the generalized mechanism governing the biogenesis of transitive siRNAs. Additionally, we present an overview of the diverse biological functions of these siRNAs and raise some pressing questions in this area for further investigation.

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.

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.

Functional analysis of soybean miR156 and miR172 in tobacco highlights their role in plant morphology and floral transition

Soybean (Glycine max), a significant oilseed and protein source for humans and livestock feed, needs short day photoperiod for floral induction. Further, soybean has a paleopolyploid genome with multiple copies of flowering genes adding to the complexity of genetic regulation of flowering, and seed set, especially in investigating the role of the noncoding genome. microRNAs, a class of noncoding RNA, play a regulatory role in plant development. miR156 and miR172 are major components of the essential regulatory hub controlling juvenile and vegetative developments and initiation of reproductive phase change leading to flowering. These microRNAs have been originally isolated and studied from model plant, Arabidopsis. However, a study on soybean microRNAs is lacking. We investigated the temporal expression patterns of gma-miR156a and gma-miR172a and found inversely related - gma-miR156a expression was higher in the vegetative stage, and gma-miR172a expression was elevated under inductive flowering conditions. The functions of gma-miR156a and gma-miR172a were evaluated via heterologous expressions in transgenic tobacco plants (Nicotiana tabacum L.). The analysis of overexpression transgenic lines highlighted that gma-miR156a plays a role in juvenile development via repression of the SPL transcription factor family. In contrast, gma-miR172a plays a pivotal role in the reproductive development phase by down-regulating its target genes, AP2. In addition, ectopic expression of gma-miR156a and gma-miR172a affected plant morphology and physiology during plant growth. Collectively, our results suggest that gma-miR156a and gma-miR172a regulate multiple morpho-physiological traits that could be used to enhance crop yield under changing climate conditions.

RNA-Based Vaccination of Plants for Control of Viruses

Plant viruses cause nearly half of the emerging plant diseases worldwide, contributing to 10-15% of crop yield losses. Control of plant viral diseases is mainly accomplished by extensive chemical applications targeting the vectors (i.e., insects, nematodes, fungi) transmitting these viruses. However, these chemicals have a significant negative effect on human health and the environment. RNA interference is an endogenous, cellular, sequence-specific RNA degradation mechanism in eukaryotes induced by double-stranded RNA molecules that has been exploited as an antiviral strategy through transgenesis. Because genetically modified crop plants are not accepted for cultivation in several countries globally, there is an urgent demand for alternative strategies. This has boosted research on exogenous application of the RNA-based biopesticides that are shown to exhibit significant protective effect against viral infections. Such environment-friendly and efficacious antiviral agents for crop protection will contribute to global food security, without adverse effects on human health.

Role of phasiRNAs in plant-pathogen interactions: molecular perspectives and bioinformatics tools

The genome of an organism is regulated in concert with the organized action of various genetic regulators at different hierarchical levels. Small non-coding RNAs are one of these regulators, among which microRNAs (miRNAs), a distinguished sRNA group with decisive functions in the development, growth and stress-responsive activities of both plants as well as animals, are keenly explored over a good number of years. Recent studies in plants revealed that apart from the silencing activity exhibited by miRNAs on their targets, miRNAs of specific size and structural features can direct the phasing pattern of their target loci to form phased secondary small interfering RNAs (phasiRNAs). These trigger-miRNAs were identified to target both coding and long non-coding RNAs that act as potent phasiRNA precursors or PHAS loci. The phasiRNAs produced thereby exhibit a role in enhancing further downstream regulation either on their own precursors or on those transcripts that are distinct from their genetic source of origin. Hence, these tiny regulators can stimulate an elaborative cascade of interacting RNA networks via cis and trans-regulatory mechanisms. Our review focuses on the comprehensive understanding of phasiRNAs and their trigger miRNAs, by giving much emphasis on their role in the regulation of plant defense responses, together with a summary of the computational tools available for the prediction of the same.

MiR1885 Regulates Disease Tolerance Genes in Brassica rapa during Early Infection with Plasmodiophora brassicae

Clubroot caused by Plasmodiophora brassicae is a severe disease of cruciferous crops that decreases crop quality and productivity. Several clubroot resistance-related quantitative trait loci and candidate genes have been identified. However, the underlying regulatory mechanism, the interrelationships among genes, and how genes are regulated remain unexplored. MicroRNAs (miRNAs) are attracting attention as regulators of gene expression, including during biotic stress responses. The main objective of this study was to understand how miRNAs regulate clubroot resistance-related genes in P. brassicae-infected Brassica rapa. Two Brassica miRNAs, Bra-miR1885a and Bra-miR1885b, were revealed to target TIR-NBS genes. In non-infected plants, both miRNAs were expressed at low levels to maintain the balance between plant development and basal immunity. However, their expression levels increased in P. brassicae-infected plants. Both miRNAs down-regulated the expression of the TIR-NBS genes Bra019412 and Bra019410, which are located at a clubroot resistance-related quantitative trait locus. The Bra-miR1885-mediated down-regulation of both genes was detected for up to 15 days post-inoculation in the clubroot-resistant line CR Shinki and in the clubroot-susceptible line 94SK. A qRT-PCR analysis revealed Bra019412 expression was negatively regulated by miR1885. Both Bra019412 and Bra019410 were more highly expressed in CR Shinki than in 94SK; the same expression pattern was detected in multiple clubroot-resistant and clubroot-susceptible inbred lines. A 5′ rapid amplification of cDNA ends analysis confirmed the cleavage of Bra019412 by Bra-miR1885b. Thus, miR1885s potentially regulate TIR-NBS gene expression during P. brassicae infections of B. rapa.

Roles of small RNAs in crop disease resistance

Small RNAs (sRNAs) are a class of short, non-coding regulatory RNAs that have emerged as critical components of defense regulatory networks across plant kingdoms. Many sRNA-based technologies, such as host-induced gene silencing (HIGS), spray-induced gene silencing (SIGS), virus-induced gene silencing (VIGS), artificial microRNA (amiRNA) and synthetic trans-acting siRNA (syn-tasiRNA)-mediated RNA interference (RNAi), have been developed as disease control strategies in both monocot and dicot plants, particularly in crops. This review aims to highlight our current understanding of the roles of sRNAs including miRNAs, heterochromatic siRNAs (hc-siRNAs), phased, secondary siRNAs (phasiRNAs) and natural antisense siRNAs (nat-siRNAs) in disease resistance, and sRNAs-mediated trade-offs between defense and growth in crops. In particular, we focus on the diverse functions of sRNAs in defense responses to bacterial and fungal pathogens, oomycete and virus in crops. Further, we highlight the application of sRNA-based technologies in protecting crops from pathogens. Further research perspectives are proposed to develop new sRNAs-based efficient strategies to breed non-genetically modified (GMO), disease-tolerant crops for sustainable agriculture.

Genome-Wide Analysis of Coding and Non-coding RNA Reveals a Conserved miR164-NAC-mRNA Regulatory Pathway for Disease Defense in Populus

MicroRNAs (miRNAs) contribute to plant defense responses by increasing the overall genetic diversity; however, their origins and functional importance in plant defense remain unclear. Here, we employed Illumina sequencing technology to assess how miRNA and messenger RNA (mRNA) populations vary in the Chinese white poplar (Populus tomentosa) during a leaf black spot fungus (Marssonina brunnea) infection. We sampled RNAs from infective leaves at conidia germinated stage [12 h post-inoculation (hpi)], infective vesicles stage (24 hpi), and intercellular infective hyphae stage (48 hpi), three essential stages associated with plant colonization and biotrophic growth in M. brunnea fungi. In total, 8,938 conserved miRNA-target gene pairs and 3,901 Populus-specific miRNA-target gene pairs were detected. The result showed that Populus-specific miRNAs (66%) were more involved in the regulation of the disease resistance genes. By contrast, conserved miRNAs (>80%) target more whole-genome duplication (WGD)-derived transcription factors (TFs). Among the 1,023 WGD-derived TF pairs, 44.9% TF pairs had only one paralog being targeted by a miRNA that could be due to either gain or loss of a miRNA binding site after the WGD. A conserved hierarchical regulatory network combining promoter analyses and hierarchical clustering approach uncovered a miR164–NAM, ATAF, and CUC (NAC) transcription factor–mRNA regulatory module that has potential in Marssonina defense responses. Furthermore, analyses of the locations of miRNA precursor sequences reveal that pseudogenes and transposon contributed a certain proportion (∼30%) of the miRNA origin. Together, these observations provide evolutionary insights into the origin and potential roles of miRNAs in plant defense and functional innovation.

Title: Hypermethylation of miRNA Genes During Nodule Development

DNA methylation has recently emerged as a powerful regulatory mechanism controlling the expression of key regulators of various developmental processes, including nodulation. However, the functional role of DNA methylation in regulating the expression of microRNA (miRNA) genes during the formation and development of nitrogen-fixing nodules remains largely unknown. In this study, we profiled DNA methylation patterns of miRNA genes during nodule formation, development, and early senescence stages in soybean (Glycine max) through the analysis of methylC-seq data. Absolute DNA methylation levels in the CG, CHH, and CHH sequence contexts over the promoter and primary transcript regions of miRNA genes were significantly higher in the nodules compared with the corresponding root tissues at these three distinct nodule developmental stages. We identified a total of 82 differentially methylated miRNAs in the nodules compared with roots. Differential DNA methylation of these 82 miRNAs was detected only in the promoter (69), primary transcript region (3), and both in the promoter and primary transcript regions (10). The large majority of these differentially methylated miRNAs were hypermethylated in nodules compared with the corresponding root tissues and were found mainly in the CHH context and showed stage-specific methylation patterns. Differentially methylated regions in the promoters of 25 miRNAs overlapped with transposable elements, a finding that may explain the vulnerability of miRNAs to DNA methylation changes during nodule development. Gene expression analysis of a set of promoter-differentially methylated miRNAs pointed to a negative association between DNA methylation and miRNA expression. Gene Ontology and pathways analyses indicate that changes in DNA methylation of miRNA genes are reprogrammed and contribute to nodule development through indirect regulation of genes involved in cellular processes and pathways with well-established roles in nodulation.

Identification of Differentially Methylated miRNA Genes During Compatible and Incompatible Interactions Between Soybean and Soybean Cyst Nematode

DNA methylation is a widespread epigenetic mark that affects gene expression and transposon mobility during plant development and stress responses. However, the role of DNA methylation in regulating the expression of microRNA (miRNA) genes remains largely unexplored. Here, we analyzed DNA methylation changes of miRNA genes using a pair of soybean (Glycine max) near-isogenic lines (NILs) differing in their response to soybean cyst nematode (SCN; Heterodera glycines). Differences in global DNA methylation levels over miRNA genes in response to SCN infection were observed between the isogenic lines. miRNA genes with significant changes in DNA methylation levels in the promoter and primary transcript-coding regions were detected in both lines. In the susceptible isogenic line (NIL-S), 82 differentially methylated miRNAs were identified in response to SCN infection whereas, in the resistant isogenic line (NIL-R), only 16 differentially methylated miRNAs were identified. Interestingly, gma-miR5032, gma-miR5043, gma-miR1520b, and gma-2107-ch16 showed opposite methylation patterns in the isogenic lines. In addition, the miRNA paralogs gma-miR5770a and gma-miR5770b showed hypermethylation and hypomethylation in NIL-S and NIL-R, respectively. Gene expression quantification of gma-miR5032, gma-miR5043, gma-miR1520b, and gma-miR5770a/b and their confirmed targets indicated a role of DNA methylation in regulating miRNA expression and, thus, their targets upon SCN infection. Furthermore, overexpression of these four miRNAs in NIL-S using transgenic hairy root system enhanced plant resistance to SCN to various degrees with a key role observed for miR5032. Together, our results provide new insights into the role of epigenetic mechanisms in controlling miRNA regulatory function during SCN-soybean interactions.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Bacillus cereus AR156 triggers induced systemic resistance against Pseudomonas syringae pv. tomato DC3000 by suppressing miR472 and activating CNLs-mediated basal immunity in Arabidopsis

Small RNAs play an important role in plant innate immunity. However, their regulatory function in induced systemic resistance (ISR) triggered by plant growth-promoting rhizobacteria remains unclear. Here, using Arabidopsis as a model system, one plant endogenous small RNA, miR472, was identified as an important regulator involved in the process of Bacillus cereus AR156 ISR against Pseudomonas syringae pv. tomato (Pst) DC3000. The results revealed that miR472 was down-regulated with B. cereus AR156 treatment by comparing small RNA profiles and northern blot analysis of Arabidopsis with or without B. cereus AR156 treatment. Plants overexpressing miR472 showed higher susceptibility to Pst DC3000; by contrast, plant lines with miR472 knocked down/out showed the opposite. The transcriptome sequencing revealed thousands of differentially expressed genes in the transgenic plants. Target prediction showed that miR472 targets lots of coiled coil nucleotide-binding site (NBS) and leucine-rich repeat (LRR) type resistance genes and the expression of these targets was negatively correlated with the expression of miR472. In addition, transgenic plants with knocked-out target genes exhibited decreased resistance to Pst DC3000 invasion. Quantitative reverse transcription PCR results indicated that target genes of miR472 were expressed during the process of B. cereus AR156-triggered ISR. Taken together, our results demonstrate that the miR472-mediated silencing pathway is an important regulatory checkpoint occurring via post-transcriptional control of NBS-LRR genes during B. cereus AR156-triggered ISR in Arabidopsis.

Comparative Analysis of miRNA Expression Profiles Between Heat-Tolerant and Heat-Sensitive Genotypes of Flowering Chinese Cabbage Under Heat Stress Using High-Throughput Sequencing

Heat stress disturbs cellular homeostasis, thus usually impairs yield of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). MicroRNAs (miRNAs) play a significant role in plant responses to different stresses by modulating gene expression at the post-transcriptional level. However, the roles that miRNAs and their target genes may play in heat tolerance of flowering Chinese cabbage remain poorly characterized. The current study sequenced six small RNA libraries generated from leaf tissues of flowering Chinese cabbage collected at 0, 6, and 12 h after 38 °C heat treatment, and identified 49 putative novel miRNAs and 43 known miRNAs that differentially expressed between heat-tolerant and heat-sensitive flowering Chinese cabbage. Among them, 14 novel and nine known miRNAs differentially expressed only in the heat-tolerant genotype under heat-stress, therefore, their target genes including disease resistance protein TAO1-like, RPS6, reticuline oxidase-like protein, etc. might play important roles in enhancing heat-tolerance. Gene Ontology (GO) analysis revealed that targets of these differentially expressed miRNAs may play key roles in responses to temperature stimulus, cell part, cellular process, cell, membrane, biological regulation, binding, and catalytic activities. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified their important functions in signal transduction, environmental adaptation, global and overview maps, as well as in stress adaptation and in MAPK signaling pathways such as cell death. These findings provide insight into the functions of the miRNAs in heat stress tolerance of flowering Chinese cabbage.

Stress-responsive miRNAome of Glycine max (L.) Merrill: molecular insights and way forward

MAIN CONCLUSION

Analysis of stress-associated miRNAs of Glycine max (L.) Merrill reveals wider ramifications of small RNA-mediated (conserved and legume-specific miRNAs) gene regulatory foot prints in molecular adaptive responses. MicroRNAs (miRNAs) are indispensable components of gene regulatory mechanism of plants. Soybean is a crop of immense commercial potential grown worldwide for its edible oil and soy meal. Intensive research efforts, using the next generation sequencing and bioinformatics techniques, have led to the identification and characterization of numerous small RNAs, especially microRNAs (miRNAs), in soybean. Furthermore, studies have unequivocally demonstrated the significance of miRNAs during the developmental processes and various stresses in soybean. In this review, we summarize the current state of understanding of miRNA-based abiotic and biotic stress responses in soybean. In addition, the molecular insights gained from the stress-related soybean miRNAs have been compared to the miRNAs of other crops, especially legumes, and the core commonalities have been highlighted, though differences among them were not ignored. Nature of response of soybean-derived conserved miRNAs during various stresses was also analyzed to gain deeper insights regarding sRNAome-based defense responses. This review further provides way forward in legume small RNA transcriptomics based on the adaptive responses of soybean and other legume-derived miRNAs.

Fine Mapping and Identification of a Novel Phytophthora Root Rot Resistance Locus RpsZS18 on Chromosome 2 in Soybean

Phytophthora root rot (PRR) caused by Phytophthora sojae is a major soybean disease that causes severe economic losses worldwide. Using soybean cultivars carrying a Rps resistance gene is the most effective strategy for controlling this disease. We previously detected a novel Phytophthora resistance gene, RpsZS18, on chromosome 2 of the soybean cultivar Zaoshu18. The aim of the present study was to identify and finely map RpsZS18. We used 232 F2:3 families generated from a cross between Zaoshu18 (resistant) and Williams (susceptible) as the mapping population. Simple sequence repeat (SSR) markers distributed on chromosome 2 were used to map RpsZS18. First, 12 SSR markers linked with RpsZS18 were identified by linkage analyses, including two newly developed SSR markers, ZCSSR33 and ZCSSR46, that flanked the gene at distances of 0.9 and 0.5 cM, respectively. Second, PCR-based InDel markers were developed based on sequence differences between the two parents and used to further narrow down the mapping region of RpsZS18 to 71.3 kb. Third, haplotype analyses were carried out in the RpsZS18 region using 14 soybean genotypes with whole-genome resequencing. We detected six genes with unique haplotype sequences in Zaoshu18. Finally, quantitative real-time PCR assays of the six genes revealed an EF-hand calcium-binding domain containing protein encoding gene (Glyma.02g245700), a pfkB carbohydrate kinase encoding gene (Glyma.02g245800), and a gene with no functional annotation (Glyma.02g246300), are putative candidate PRR resistance genes. This study provides useful information for breeding P. sojae-resistant soybean cultivars.

miR3954 is a trigger of phasiRNAs that affects flowering time in citrus

In plant, a few 22-nt miRNAs direct cleavages of their targets and trigger the biogenesis of phased small interfering RNAs (phasiRNAs) in plant. In this study, we characterized a miRNA triggering phasiRNAs generation, miR3954, and explored its downstream target genes and potential function. Our results demonstrated that miR3954 showed specific expression in the flowers of citrus species, and it targeted a NAC transcription factor (Cs7 g22460) and two non-coding RNA transcripts (lncRNAs, Cs1 g09600 and Cs1 g09635). The production of phasiRNAs was detected from transcripts targeted by miR3954, and was further verified in both sequencing data and transient expression experiments. PhasiRNAs derived from the two lncRNAs targeted not only miR3954-targeted NAC gene but also additional NAC homologous genes. No homologous genes of these two lncRNAs were found in plants other than citrus species, implying that this miR3954-lncRNAs-phasiRNAs-NAC pathway is likely citrus-specific. Transgenic analysis indicated that the miR3954-overexpressing lines showed decreased transcripts of lncRNA, elevated abundance of phasiRNAs and reduced expression of NAC genes. Interestingly, the overexpression of miR3954 leads to early flowering in citrus plants. In summary, our results illustrated a model of the regulatory network of miR3954-lncRNA-phasiRNAs-NAC, which may be functionally involved in flowering in citrus.

Small RNAs regulate plant responses to filamentous pathogens

Small RNAs are central players of RNA silencing in eukaryotes. These short RNA molecules (20-25 nucleotides in length) repress target gene expression based on sequence complementarity. While small RNAs are well-known for their essential function in regulating growth and development, recent research has revealed that they also influence plant immunity. Extensive changes in small RNA accumulation have been observed during infection. This review focuses on specific small RNA changes that are involved in plant responses to filamentous eukaryotic pathogens including fungi and oomycetes. We describe how changes in small RNA accumulation influence plant immunity and summarize the cellular processes affected by these small RNAs. In particular, we discuss secondary small interfering RNAs that directly modulate the expression of defense-related genes.

Transcriptome and Small RNAome Dynamics during a Resistant and Susceptible Interaction between Cucumber and Downy Mildew

Cucumber ( L.) downy mildew, caused by the obligate oomycete pathogen (Berk. and Curt.) Rostov., is the primary factor limiting cucumber production. Although sources of resistance have been identified, such as plant introduction line PI 197088, the genes and processes involved in mediating resistance are still unknown. In the current study, we conducted a comprehensive transcriptome and small RNAome analysis of a resistant (PI 197088) and susceptible ('Vlaspik') cucumber during a time course of infection using Illumina sequencing. We identified significantly differentially expressed (DE) genes within and between resistant and susceptible cucumber leaves over a time course of infection. Weighted gene correlation network analyses (WGCNA) created coexpression modules containing genes with unique expression patterns between Vlaspik and PI 197088. Recurring data trends indicated that resistance to cucumber downy mildew is associated with earlier response to the pathogen, hormone signaling, and regulation of nutrient supply. Candidate resistance genes were identified from multiple transcriptome analyses and literature support. Additionally, parallel sequencing of small RNAs (sRNAs) from cucumber and during the infection time course was used to identify and quantify novel and existing microRNA (miRNA) in both species. Predicted miRNA targets of cucumber transcripts suggest a complex interconnectedness of gene expression regulation in this plant-pathogen system. This work bioinformatically uncovered gene expression patterns involved in the mediation of or response to resistance. Herein, we provide the foundation for future work to validate candidate resistance genes and miRNA-based regulation proposed in this study.

Global investigation of the co-evolution of MIRNA genes and microRNA targets during soybean domestication

Although the selection of coding genes during plant domestication has been well studied, the evolution of MIRNA genes (MIRs) and the interaction between microRNAs (miRNAs) and their targets in this process are poorly understood. Here, we present a genome-wide survey of the selection of MIRs and miRNA targets during soybean domestication and improvement. Our results suggest that, overall, MIRs have higher evolutionary rates than miRNA targets. Nonetheless, they do demonstrate certain similar evolutionary patterns during soybean domestication: MIRs and miRNA targets with high expression and duplication status, and with greater numbers of partners, exhibit lower nucleotide divergence than their counterparts without these characteristics, suggesting that expression level, duplication status, and miRNA-target interaction are essential for evolution of MIRs and miRNA targets. Further investigation revealed that miRNA-target pairs that are subjected to strong purifying selection have greater similarities than those that exhibited genetic diversity. Moreover, mediated by domestication and improvement, the similarities of a large number of miRNA-target pairs in cultivated soybean populations were increased compared to those in wild soybeans, whereas a small number of miRNA-target pairs exhibited decreased similarity, which may be associated with the adoption of particular domestication traits. Taken together, our results shed light on the co-evolution of MIRs and miRNA targets during soybean domestication.

Disease Resistance Gene Analogs (RGAs) in Plants

Plants have developed effective mechanisms to recognize and respond to infections caused by pathogens. Plant resistance gene analogs (RGAs), as resistance (R) gene candidates, have conserved domains and motifs that play specific roles in pathogens' resistance. Well-known RGAs are nucleotide binding site leucine rich repeats, receptor like kinases, and receptor like proteins. Others include pentatricopeptide repeats and apoplastic peroxidases. RGAs can be detected using bioinformatics tools based on their conserved structural features. Thousands of RGAs have been identified from sequenced plant genomes. High-density genome-wide RGA genetic maps are useful for designing diagnostic markers and identifying quantitative trait loci (QTL) or markers associated with plant disease resistance. This review focuses on recent advances in structures and mechanisms of RGAs, and their identification from sequenced genomes using bioinformatics tools. Applications in enhancing fine mapping and cloning of plant disease resistance genes are also discussed.