miR396-OsGRFs Module Balances Growth and Rice Blast Disease-Resistance

Refinement of rice blast disease resistance QTLs and gene networks through meta-QTL analysis

Rice blast disease is the most devastating disease constraining crop productivity. Vertical resistance to blast disease is widely studied despite its instability. Clusters of genes or QTLs conferring blast resistance that offer durable horizontal resistance are important in resistance breeding. In this study, we aimed to refine the reported QTLs and identify stable meta-QTLs (MQTLs) associated with rice blast resistance. A total of 435 QTLs were used to project 71 MQTLs across all the rice chromosomes. As many as 199 putative rice blast resistance genes were identified within 53 MQTL regions. The genes included 48 characterized resistance gene analogs and related proteins, such as NBS-LRR type, LRR receptor-like kinase, NB-ARC domain, pathogenesis-related TF/ERF domain, elicitor-induced defense and proteins involved in defense signaling. MQTL regions with clusters of RGA were also identified. Fifteen highly significant MQTLs included 29 candidate genes and genes characterized for blast resistance, such as Piz, Nbs-Pi9, pi55-1, pi55-2, Pi3/Pi5-1, Pi3/Pi5-2, Pikh, Pi54, Pik/Pikm/Pikp, Pb1 and Pb2. Furthermore, the candidate genes (42) were associated with differential expression (in silico) in compatible and incompatible reactions upon disease infection. Moreover, nearly half of the genes within the MQTL regions were orthologous to those in O. sativa indica, Z. mays and A. thaliana, which confirmed their significance. The peak markers within three significant MQTLs differentiated blast-resistant and susceptible lines and serve as potential surrogates for the selection of blast-resistant lines. These MQTLs are potential candidates for durable and broad-spectrum rice blast resistance and could be utilized in blast resistance breeding.

Rice breeding for low input agriculture

A low-input-based farming system can reduce the adverse effects of modern agriculture through proper utilization of natural resources. Modern varieties often need to improve in low-input settings since they are not adapted to these systems. In addition, rice is one of the most widely cultivated crops worldwide. Enhancing rice performance under a low input system will significantly reduce the environmental concerns related to rice cultivation. Traits that help rice to maintain yield performance under minimum inputs like seedling vigor, appropriate root architecture for nutrient use efficiency should be incorporated into varieties for low input systems through integrated breeding approaches. Genes or QTLs controlling nutrient uptake, nutrient assimilation, nutrient remobilization, and root morphology need to be properly incorporated into the rice breeding pipeline. Also, genes/QTLs controlling suitable rice cultivars for sustainable farming. Since several variables influence performance under low input conditions, conventional breeding techniques make it challenging to work on many traits. However, recent advances in omics technologies have created enormous opportunities for rapidly improving multiple characteristics. This review highlights current research on features pertinent to low-input agriculture and provides an overview of alternative genomics-based breeding strategies for enhancing genetic gain in rice suitable for low-input farming practices.

Genome-wide characterization and function analysis of ginger (Zingiber officinale Roscoe) ZoGRFs in responding to adverse stresses

Growth-regulating factors (GRFs) play crucial roles in plant growth, development, hormone signaling, and stress response. Despite their significance, the roles of GRFs in ginger remain largely unknown. Herein, 31 ginger ZoGRFs were identified and designated as ZoGRF1-ZoGRF31 according to their phylogenetic relationships. All ZoGRFs were characterized as unstable, hydrophilic proteins, with 29 predicted to be located in the nucleus. Functional cis-elements related to growth and development were enriched in ZoGRF's promoter regions. RNA-seq and RT-qPCR analysis revealed that ZoGRF12, ZoGRF24, and ZoGRF28 were highly induced in various growth and development stages, displaying differential regulation under waterlogging, chilling, drought, and salt stresses, indicating diverse expression patterns of ZoGRFs. Transient expression analysis in Nicotiana benthamiana indicated that overexpressing ZoGRF28 regulated the transcription levels of salicylic acid, jasmonic acid, and pattern-triggered immunity-related genes, increased chlorophyll content and contributed to reduced disease lesions and an increased net photosynthetic rate. This research lays the foundation for further understanding the biological roles of ZoGRFs.

microRNAs: Key Regulators in Plant Responses to Abiotic and Biotic Stresses via Endogenous and Cross-Kingdom Mechanisms

Dramatic shifts in global climate have intensified abiotic and biotic stress faced by plants. Plant microRNAs (miRNAs)-20-24 nucleotide non-coding RNA molecules-form a key regulatory system of plant gene expression; playing crucial roles in plant growth; development; and defense against abiotic and biotic stress. Moreover, they participate in cross-kingdom communication. This communication encompasses interactions with other plants, microorganisms, and insect species, collectively exerting a profound influence on the agronomic traits of crops. This article comprehensively reviews the biosynthesis of plant miRNAs and explores their impact on plant growth, development, and stress resistance through endogenous, non-transboundary mechanisms. Furthermore, this review delves into the cross-kingdom regulatory effects of plant miRNAs on plants, microorganisms, and pests. It proceeds to specifically discuss the design and modification strategies for artificial miRNAs (amiRNAs), as well as the protection and transport of miRNAs by exosome-like nanovesicles (ELNVs), expanding the potential applications of plant miRNAs in crop breeding. Finally, the current limitations associated with harnessing plant miRNAs are addressed, and the utilization of synthetic biology is proposed to facilitate the heterologous expression and large-scale production of miRNAs. This novel approach suggests a plant-based solution to address future biosafety concerns in agriculture.

Unveiling the Roles of LncRNA MOIRAs in Rice Blast Disease Resistance

Rice blast disease, caused by the fungal pathogen Magnaporthe oryzae, is a major threat to rice production worldwide. This study investigates the role of long non-coding RNAs (lncRNAs) in rice's response to this destructive disease, with a focus on their impacts on disease resistance and yield traits. Three specific lncRNAs coded by M. oryzae infection-responsive lncRNAs (MOIRAs), MOIRA1, MOIRA2, and MOIRA3, were identified as key regulators of rice's response to M. oryzae infection. Strikingly, when MOIRA1 and MOIRA2 were overexpressed, they exhibited a dual function: they increased rice's susceptibility to blast fungus, indicating a negative role in disease resistance, while simultaneously enhancing tiller numbers and single-plant yield, with no adverse effects on other yield-related traits. This unexpected improvement in productivity suggests the possibility of overcoming the traditional trade-off between disease resistance and crop yield. These findings provide a novel perspective on crop enhancement, offering a promising solution to global food security challenges by developing rice varieties that effectively balance disease resistance and increased productivity.

Growing Conditions and Varietal Ecologies Differently Regulates the Growth-regulating-factor (GRFs) Gene Family in Rice

Background

Growth-regulating factors (GRFs) are crucial in rice for controlling plant growth and development. Among the rice cultivation practices, aerobic methods are water efficient but result in significant yield reduction relative to non-aerobic cultivation. Therefore, mechanistic insights into aerobic rice cultivation are important for improving the aerobic performance of rice.

Objectives

This study aimed to examine the evolution of GRFs in different rice species, analyse the phenotypic differences between aerobic and non-aerobic conditions in three rice varieties, and assess the expression of GRFs in these varieties under both aerobic and non-aerobic conditions.

Materials and Methods

This study comprehensively examined the GRFs gene family in 11 rice species (Oryza barthii, Oryza brachyantha, Oryza glaberrima, Oryza glumipatula, Oryza sativa subsp. indica, Oryza longistaminata, Oryza meridionalis, Oryza nivara, Oryza punctata, Oryza rufipogon, Oryza sativa subsp. japonica) focusing on phylogenetic analysis. Additionally, the expression patterns of 12 GRFs were investigated in three distinct genotypes of O. sativa subsp. indica rice, under both non-aerobic and aerobic conditions.

Results

Three major phylogenetic clades were formed based on conserved motifs in the 123 GRFs proteins in eleven rice species. Further, novel motifs were identified especially in O. longistaminata indicative of the species level evolutionary differences in rice. Among the trait performance, the number of tillers was reduced by ~ 36% under aerobic conditions, but the reduction was found to be less in CR Dhan 201, an aerobic variety. Besides, three GRFs namely GRF3, GRF4, and GRF7 were found to be distinct in expression between aerobic and non-aerobic conditions.

Conclusion

Three GRF genes namely GRF3, GRF4, and GRF7 could be associated with the aerobic adaptation in rice.

Bacillus subtilis 26D Triggers Induced Systemic Resistance against Rhopalosiphum padi L. by Regulating the Expression of Genes AGO, DCL and microRNA in Bread Spring Wheat

Bacillus subtilis 26D is a plant growth-promoting endophytic bacteria capable of inducing systemic resistance through the priming mechanism, which includes plant genome reprogramming and the phenomenon of RNA interference (RNAi) and microRNA (miRNAs). The phloem-feeding insect bird cherry-oat aphid Rhopalosiphum padi L. is a serious pest that causes significant damage to crops throughout the world. However, the function of plant miRNAs in the response to aphid infestation remains unclear. The results of this work showed that B. subtilis 26D stimulated aphid resistance in wheat plants, inducing the expression of genes of hormonal signaling pathways ICS, WRKY13, PR1, ACS, EIN3, PR3, and ABI5. In addition, B. subtilis 26D activated the RNAi mechanism and regulated the expression of nine conserved miRNAs through activation of the ethylene, salicylic acid (SA), and abscisic acid (ABA) signaling pathways, which was demonstrated by using treatments with phytohormones. Treatment of plants with SA, ethylene, and ABA acted in a similar manner to B. subtilis 26D on induction of the expression of the AGO4, AGO5 and DCL2, DCL4 genes, as well as the expression of nine conserved miRNAs. Different patterns of miRNA expression were found in aphid-infested plants and in plants treated with B. subtilis 26D or SA, ethylene, and ABA and infested by aphids, suggesting that miRNAs play multiple roles in the plant response to phloem-feeding insects, associated with effects on hormonal signaling pathways, redox metabolism, and the synthesis of secondary metabolites. Our study provides new data to further elucidate the fine mechanisms of bacterial-induced priming. However, further extensive work is needed to fully unravel these mechanisms.

Genome-wide identification, characterization, and expression analysis unveil the roles of pseudouridine synthase (PUS) family proteins in rice development and stress response

Pseudouridylation, the conversion of uridine (U) to pseudouridine (Ѱ), is one of the most prevalent and evolutionary conserved RNA modifications, which is catalyzed by pseudouridine synthase (PUS) enzymes. Ѱs play a crucial epitranscriptomic role by regulating attributes of cellular RNAs across diverse organisms. However, the precise biological functions of PUSs in plants remain largely elusive. In this study, we identified and characterized 21 members in the rice PUS family which were categorized into six distinct subfamilies, with RluA and TruA emerging as the most extensive. A comprehensive analysis of domain structures, motifs, and homology modeling revealed that OsPUSs possess all canonical features of true PUS proteins, essential for substrate recognition and catalysis. The exploration of OsPUS promoters revealed presence of cis-acting regulatory elements associated with hormone and abiotic stress responses. Expression analysis of OsPUS genes showed differential expression at developmental stages and under stress conditions. Notably, OsTruB3 displayed high expression in salt, heat, and drought stresses. Several OsRluA members showed induction in heat stress, while a significant decline in expression was observed for various OsTruA members in drought and salinity. Furthermore, miRNAs predicted to target OsPUSs were themselves responsive to variable stresses, adding an additional layer of regulatory complexity of OsPUSs. Study of protein-protein interaction networks provided substantial support for the potential regulatory role of OsPUSs in numerous cellular and stress response pathways. Conclusively, our study provides functional insights into the OsPUS family, contributing to a better understanding of their crucial roles in shaping the development and stress adaptation in rice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01396-4.

Unraveling the transcriptional network regulated by miRNAs in blast-resistant and blast-susceptible rice genotypes during Magnaporthe oryzae interaction

The plant pathogen Magnaporthe oryzae poses a significant threat to global food security, and its management through the cultivation of resistant varieties and crop husbandry practices, including fungicidal sprays, has proven to be inadequate. To address this issue, we conducted small-RNA sequencing to identify the roles of miRNAs and their target genes in both resistant (PB1637) and susceptible (PB1) rice genotypes. We confirmed the expression of differentially expressed miRNAs using stem-loop qRT-PCR analysis and correlated them with rice patho-phenotypic and physio-biochemical responses. Our findings revealed several noteworthy differences between the resistant and susceptible genotypes. The resistant genotype exhibited reduced levels of total chlorophyll and carotenoids compared to the susceptible genotype. However, it showed increased levels of total protein, callose, H2O2, antioxidants, flavonoids, and total polyphenols. Additionally, among the defense-associated enzymes, guaiacol peroxidase and polyphenol oxidase responses were higher in the susceptible genotypes. In our comparative analysis, we identified 27 up-regulated and 43 down-regulated miRNAs in the resistant genotype, while the susceptible genotype exhibited 44 up-regulated and 62 down-regulated miRNAs. Furthermore, we discovered eight up-regulated and five down-regulated miRNAs shared between the resistant and susceptible genotypes. Notably, we also identified six novel miRNAs in the resistant genotype and eight novel miRNAs in the susceptible genotype. These novel miRNAs, namely Chr8_26996, Chr12_40110, and Chr12_41899, were found to negatively correlate with the expression of predicted target genes, including Cyt-P450 monooxygenase, serine carboxypeptidase, and zinc finger A20 domain-containing stress-associated protein, respectively. The results of our study on miRNA and transcriptional responses provide valuable insights for the development of future rice lines that are resistant to blast disease. By understanding the roles of specific miRNAs and their target genes in conferring resistance, we can enhance breeding strategies and improve crop management practices to ensure global food security.

OsmiRNA5488 Regulates the Development of Embryo Sacs and Targets OsARF25 in Rice (Oryza sativa L.)

Small RNAs are a class of non-coding RNAs that typically range from 20 to 24 nucleotides in length. Among them, microRNAs (miRNAs) are particularly important regulators for plant development. The biological function of the conserved miRNAs has been studied extensively in plants, while that of the species-specific miRNAs has been studied in-depth. In this study, the regulatory role of a rice-specific OsmiRNA5488 (OsmiR5488) was characterized with the miR5488-overexpressed line (miR5488-OE) and miR5488-silenced line (STTM-5488). The seed-setting rate was notably reduced in miR5488-OE lines, but not in STTM-5488 lines. Cytological observation demonstrated the different types of abnormal mature embryo sacs, including the degeneration of embryo sacs and other variant types, in miR5488-OE lines. The percentage of the abnormal mature embryo sacs accounted for the reduced value of the seed-setting rate. Furthermore, OsARF25 was identified as a target of OsmiR5488 via RNA ligase-mediated 3'-amplifification of cDNA ends, dual luciferase assays, and transient expression assays. The primary root length was decreased with the increases in auxin concentrations in miR5488-OE lines compared to wild-type rice. Summarily, our results suggested that OsmiR5488 regulates the seed-setting rate and down-regulates the targeted gene OsARF25.

Osa-miR11117 Targets OsPAO4 to Regulate Rice Immunity against the Blast Fungus Magnaporthe oryzae

The intricate regulatory process governing rice immunity against the blast fungus Magnaporthe oryzae remains a central focus in plant-pathogen interactions. In this study, we investigated the important role of Osa-miR11117, an intergenic microRNA, in regulating rice defense mechanisms. Stem-loop qRT-PCR analysis showed that Osa-miR11117 is responsive to M. oryzae infection, and overexpression of Osa-miR11117 compromises blast resistance. Green fluorescent protein (GFP)-based reporter assay indicated OsPAO4 is one direct target of Osa-miR11117. Furthermore, qRT-PCR analysis showed that OsPAO4 reacts to M. oryzae infection and polyamine (PA) treatment. In addition, OsPAO4 regulates rice resistance to M. oryzae through the regulation of PA accumulation and the expression of the ethylene (ETH) signaling genes. Taken together, these results suggest that Osa-miR11117 is targeting OsPAO4 to regulate blast resistance by adjusting PA metabolism and ETH signaling pathways.

Development and Genome-Wide Analysis of a Blast-Resistant japonica Rice Variety

Rice is one of the most important crops in the world, and its production is severely affected by the rice blast disease caused by the fungus Magnaporthe oryzae. Several major blast resistance genes and QTLs associated with blast resistance have been described and mostly identified in indica rice varieties. In this work, we report the obtention of a blast-resistant rice breeding line derived from crosses between the resistant indica variety CT13432 and the japonica elite cultivar JSendra (highly susceptible to blast). The breeding line, named COPSEMAR9, was found to exhibit resistance to leaf blast and panicle blast, as demonstrated by disease assays under controlled and field conditions. Furthermore, a high-quality genome sequence of the blast-resistant breeding line was obtained using a strategy that combines short-read sequencing (Illumina sequencing) and long-read sequencing (Pacbio sequencing). The use of a whole-genome approach allowed the fine mapping of DNA regions of indica and japonica origin present in the COPSEMAR9 genome and the identification of parental gene regions potentially contributing to blast resistance in the breeding line. Rice blast resistance genes (including Pi33 derived from the resistant parent) and defense-related genes in the genome of COPSEMAR9 were identified. Whole-genome analyses also revealed the presence of microRNAs (miRNAs) with a known function in the rice response to M. oryzae infection in COPSEMAR9, which might also contribute to its phenotype of blast resistance. From this study, the genomic information and analysis methods provide valuable knowledge that will be useful in breeding programs for blast resistance in japonica rice cultivars.

Cytological observation and RNA-seq analysis reveal novel miRNAs high expression associated with the pollen fertility of neo-tetraploid rice

BACKGROUND

Neo-tetraploid rice lines exhibit high fertility and strong heterosis and harbor novel specific alleles, which are useful germplasm for polyploid rice breeding. However, the mechanism of the fertility associated with miRNAs remains unknown. In this study, a neo-tetraploid rice line, termed Huaduo21 (H21), was used. Cytological observation and RNA-sequencing were employed to identify the fertility-related miRNAs in neo-tetraploid rice.

RESULTS

H21 showed high pollen fertility (88.08%), a lower percentage of the pollen mother cell (PMC) abnormalities, and lower abnormalities during double fertilization and embryogenesis compared with autotetraploid rice. A total of 166 non-additive miRNAs and 3108 non-additive genes were detected between H21 and its parents. GO and KEGG analysis of non-additive genes revealed significant enrichments in the DNA replication, Chromosome and associated proteins, and Replication and repair pathways. Comprehensive multi-omics analysis identified 32 pairs of miRNA/target that were associated with the fertility in H21. Of these, osa-miR408-3p and osa-miR528-5p displayed high expression patterns, targeted the phytocyanin genes, and were associated with high pollen fertility. Suppression of osa-miR528-5p in Huaduo1 resulted in a low seed set and a decrease in the number of grains. Moreover, transgenic analysis implied that osa-MIR397b-p3, osa-miR5492, and osa-MIR5495-p5 might participate in the fertility of H21.

CONCLUSION

Taken together, the regulation network of fertility-related miRNAs-targets pairs might contribute to the high seed setting in neo-tetraploid rice. These findings enhance our understanding of the regulatory mechanisms of pollen fertility associated with miRNAs in neo-tetraploid rice.

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.

Combined analysis of microRNA and mRNA profiles provides insights into the pathogenic resistant mechanisms of the lotus rhizome rot

Lotus rhizome rot caused by Fusarium oxysporum is a common vascular fungal disease in plants that significantly impacts the yield. However, only a few studies have studied the mechanism of Nelumbo nucifera responding to lotus rhizome rot. Here, we investigated the pathogenic genes and miRNAs in lotus rhizome rot to uncover the pathogenic resistant mechanisms by transcriptome and small RNA sequencing of lotus roots after inoculation with Fusarium oxysporum. GO and KEGG functional enrichment analysis showed that differential miRNAs were mostly enriched in starch and sucrose metabolism, biosynthesis of secondary metabolites, glutathione metabolism, brassinosteroid biosynthesis and flavonoid biosynthesis pathways. Twenty-seven upregulated miRNAs, 19 downregulated miRNAs and their target genes were identified. Correlation analysis found that miRNAs negatively regulate target genes, which were also enriched in starch and sucrose metabolism and glutathione metabolism pathways. Their expression was measured by reverse transcription quantitative PCR (qRT-PCR), and the results were consistent with the transcriptome analysis, thus verifying the reliability of transcriptome data. We selected three miRNAs (miRNA858-y, miRNA171-z and a novel miRNA novel-m0005-5p) to test the relationship between miRNAs and their target genes. The activity of the GUS testing assay indicated that miRNA could decrease the GUS activity by inhibiting the expression of their target genes. Collectively, this study provides a comprehensive analysis of transcriptome and small RNA sequencing of lotus root after inoculation with Fusarium oxysporum, and we identified candidate miRNAs and their target genes for breeding strategies of Nelumbo nucifera.

Comprehensive identification and characterization of lncRNAs and circRNAs reveal potential brown planthopper-responsive ceRNA networks in rice

Brown planthopper (Nilaparvata lugens Stål, BPH) is one of the most destructive pests of rice. Non-coding RNA plays an important regulatory role in various biological processes. However, comprehensive identification and characterization of long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) in BPH-infested rice have not been performed. Here, we performed a genome-wide analysis of lncRNAs and circRNAs in BPH6-transgenic (resistant, BPH6G) and Nipponbare (susceptible, NIP) rice plants before and after BPH feeding (early and late stage) via deep RNA-sequencing. A total of 310 lncRNAs and 129 circRNAs were found to be differentially expressed. To reveal the different responses of resistant and susceptible rice to BPH herbivory, the potential functions of these lncRNAs and circRNAs as competitive endogenous RNAs (ceRNAs) were predicted and investigated using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses. Dual-luciferase reporter assays revealed that miR1846c and miR530 were targeted by the lncRNAs XLOC_042442 and XLOC_028297, respectively. In responsive to BPH infestation, 39 lncRNAs and 21 circRNAs were predicted to combine with 133 common miRNAs and compete for miRNA binding sites with 834 mRNAs. These mRNAs predictably participated in cell wall organization or biogenesis, developmental growth, single-organism cellular process, and the response to stress. This study comprehensively identified and characterized lncRNAs and circRNAs, and integrated their potential ceRNA functions, to reveal the rice BPH-resistance network. These results lay a foundation for further study on the functions of lncRNAs and circRNAs in the rice-BPH interaction, and enriched our understanding of the BPH-resistance response in rice.

Osa-miR162a Enhances the Resistance to the Brown Planthopper via α-Linolenic Acid Metabolism in Rice (Oryza sativa)

The brown planthopper (BPH) is the most serious pest causing yield losses in rice. MicroRNAs (miRNAs) are emerging as key modulators of plant-pest interactions. In the study, we found that osa-miR162a is induced in response to BPH attack in the seedling stage and tunes rice resistance to the BPH via the α-linolenic acid metabolism pathway as indicated by gas chromatography/liquid chromatography-mass spectrometry analysis. Overexpression of osa-miR162a inhibited the development and growth of the BPH and simultaneously reduced the release of 3-hexenal and 3-hexen-1-ol to block host recognition in the BPH. Moreover, knockdown of OsDCL1, which is targeted by osa-miR162a, inhibited α-linolenic acid metabolism to enhance the resistance to the BPH, which was similar to that in miR162a-overexpressing plants. Our study revealed a novel defense mechanism mediated by plant miRNAs developed during the long-term evolution of plant-host interaction, provided new ideas for the identification of rice resistance resources, and promoted a better understanding of pest control.

Recent perspective of non-coding RNAs at the nexus of plant-pathogen interaction

In natural habitats, plants are exploited by pathogens in biotrophic or necrotrophic ways. Concurrently, plants have evolved their defense systems for rapid perception of pathogenic effectors and begin concerted cellular reprogramming pathways to confine the pathogens at the entry sites. During the reorganization of cellular signaling mechanisms following pathogen attack, non-coding RNAs serves an indispensable role either as a source of resistance or susceptibility. Besides the well-studied functions of non-coding RNAs related to plant development and abiotic stress responses, previous and recent discoveries have established that non-coding RNAs like miRNAs, siRNAs, lncRNAs and phasi-RNAs can fine tune plant defense responses by targeting various signaling pathways. In this review, recapitulation of previous reports associated with non-coding RNAs as a defense responder against virus, bacteria and fungus attacks and insightful discussion will lead us to conceive innovative ideas to fight against approaching threats of resistant breaking pathogens.

MicroRNAs Involved in Nutritional Regulation During Plant-Microbe Symbiotic and Pathogenic Interactions with Rice as a Model

Plants are constantly challenged with numerous adverse environmental conditions, including biotic and abiotic stresses. Coordinated regulation of plant responses requires crosstalk between regulatory pathways initiated by different external cues. Stress induced by excessiveness or deficiency of nutrients has been shown to positively or negatively interact with pathogen-induced immune responses. Also, colonization by arbuscular mycorrhizal (AM) fungi can improve plant nutrition, mainly phosphorus and resistance to pathogen infection. The proposed review addresses these issues about a new question that integrates adaptation to nutrient stress and disease resistance. The main goal of the current review is to provide insights into the interconnected regulation between nutrient signaling and immune signaling pathways in rice, focusing on phosphate, potassium and iron signaling. The underpinnings of plant/pathogen/AM fungus interaction concerning rice/M. oryzae/R. irregularis is highlighted. The role of microRNAs (miRNAs) involved in Pi (miR399, miR827) and Fe (miR7695) homeostasis in pathogenic/symbiotic interactions in rice is discussed. The intracellular dynamics of membrane proteins that function in nutrient transport transgenic rice lines expressing fluorescent protein fusion genes are outlined. Integrating functional genomic, nutritional and metal content, molecular and cell biology approaches to understand how disease resistance is regulated by nutrient status leading to novel concepts in fundamental processes underlying plant disease resistance will help to devise novel strategies for crop protection with less input of pesticides and fertilizers.

Rice microRNA156/529-SQUAMOSA PROMOTER BINDING PROTEIN-LIKE7/14/17 modules regulate defenses against bacteria

Rice (Oryza sativa L.) microRNA156/529-SQUAMOSA PROMOTER BINDING PROTEIN-LIKE7/14/17 (miR156/529-SPL7/14/17) modules have pleiotropic effects on many biological pathways. OsSPL7/14 can interact with DELLA protein SLENDER RICE1 (SLR1) to modulate gibberellin acid (GA) signal transduction against the bacterial pathogen Xanthomonas oryzae pv. oryzae. However, whether the miR156/529-OsSPL7/14/17 modules also regulate resistance against other pathogens is unclear. Notably, OsSPL7/14/17 functioning as transcriptional activators, their target genes, and the corresponding downstream signaling pathways remain largely unexplored. Here, we demonstrate that miR156/529 play negative roles in plant immunity and that miR156/529-regulated OsSPL7/14/17 confer broad-spectrum resistance against 2 devastating bacterial pathogens. Three OsSPL7/14/17 proteins directly bind to the promoters of rice Allene Oxide Synthase 2 (OsAOS2) and NONEXPRESSOR OF PATHOGENESIS-RELATED GENES1 (OsNPR1) and activate their transcription, regulating jasmonic acid (JA) accumulation and the salicylic acid (SA) signaling pathway, respectively. Overexpression of OsAOS2 or OsNPR1 impairs the susceptibility of the osspl7/14/17 triple mutant. Exogenous application of JA enhances resistance of the osspl7/14/17 triple mutant and the miR156 overexpressing plants. In addition, genetic evidence confirms that bacterial pathogen-activated miR156/529 negatively regulate pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) responses, such as pattern recognition receptor Xa3/Xa26-initiated PTI. Our findings demonstrate that bacterial pathogens modulate miR156/529-OsSPL7/14/17 modules to suppress OsAOS2-catalyzed JA accumulation and the OsNPR1-promoted SA signaling pathway, facilitating pathogen infection. The uncovered miR156/529-OsSPL7/14/17-OsAOS2/OsNPR1 regulatory network provides a potential strategy to genetically improve rice disease resistance.

OsUGE1 is directly targeted by OsGRF6 to regulate root hair length in rice

KEY MESSAGE

Root hairs are required for water and nutrient acquisition in plants. Here, we report a novel mechanism that OsUGE1 is negatively controlled by OsGRF6 to regulate root hair elongation in rice. Root hairs are tubular outgrowths generated by the root epidermal cells. They effectively enlarge the soil-root contact area and play essential roles for nutrient and water absorption. Here, in this study, we demonstrated that the Oryza sativa UDP-glucose 4-epimerase 1-like (OsUGE1) negatively regulated root hair elongation and was directly targeted by Oryza sativa growth regulating factor 6 (OsGRF6). Knockout mutants of OsUGE1 using CRISPR-Cas9 technology showed longer root hairs than those of wild type. In contrast, overexpression lines of OsUGE1 displayed shorter root hair compared with those of wild type. GUS staining showed that it could specifically express in root hair. Subcellular localization analysis indicates that OsUGE1 is located in endoplasmic reticulum, nucleus and plasma membrane. More importantly, ChIP-qPCR, Yeast-one-hybrid and BiFC experiments revealed that OsGRF6 could bind to the promoter of OsUGE1. Furthermore, knockout mutants of OsGRF6 showed shorter root hair than those of wild type, and OsGRF6 dominantly expressed in root. In addition, the expression level of OsUGE1 is significantly downregulated in Osgrf6 mutant. Taken together, our study reveals a novel pathway that OsUGE1 is negatively controlled by OsGRF6 to regulate root hair elongation in rice.

An OsPRMT5-OsAGO2/miR1875-OsHXK1 module regulates rice immunity to blast disease

Rice ARGONAUTE2 (OsAGO2) is a core component of the rice RNA-induced silencing complex (RISC), which is repressed by Magnaporthe oryzae (M. oryzae) infection. Whether and how OsAGO2-mediated gene silencing plays a role in rice blast resistance and which sRNAs participate in this process are unknown. Our results indicate that OsAGO2 is a key immune player that manipulates rice defense responses against blast disease. OsAGO2 associates with the 24-nt miR1875 and binds to the promoter region of HEXOKINASE1 (OsHXK1), which causes DNA methylation and leads to gene silencing. Our multiple genetic evidence showed that, without M. oryzae infection, OsAGO2/miR1875 RISC promoted OsHXK1 promoter DNA methylation and OsHXK1 silencing; after M. oryzae infection, the reduced OsAGO2/miR1875 led to a relatively activated OsHXK1 expression. OsHXK1 acts as a positive regulator of blast disease resistance that OsHXK1-OE rice exhibited enhanced resistance, whereas Cas9-Oshxk1 rice showed reduced resistance against M. oryzae infection. OsHXK1 may function through its sugar sensor activity as glucose induced defense-related gene expression and reactive oxygen species (ROS) accumulation in Nipponbare and OsHXK1-OE but not in Cas9-Oshxk1 rice. OsAGO2 itself is delicately regulated by OsPRMT5, which senses M. oryzae infection and attenuates OsAGO2-mediated gene silencing through OsAGO2 arginine methylation. Our study reveals an OsPRMT5-OsAGO2/miR1875-OsHXK1 regulatory module that fine tunes the rice defense response to blast disease.

Enhancing alfalfa resistance to Spodoptera herbivory by sequestering microRNA396 expression

KEY MESSAGE

Sequestering microRNA396 by overexpression of MIM396 enhanced alfalfa resistance to Spodoptera litura larvae, which may be due to increased lignin content and enhanced low-molecular weight flavonoids and glucosinolates biosynthesis. Alfalfa (Medicago sativa), the most important leguminous forage crop, suffers from the outbreak of defoliator insects, especially Spodoptera litura, resulting in heavy losses in yield and forage quality. Here, we found that the expression of alfalfa microRNA396 (miR396) precursor genes and mature miR396 was significantly up-regulated in wounding treatment that simulates feeding injury by defoliator insects. To verify the function of miR396 in alfalfa resistance to insect, we generated MIM396 transgenic alfalfa plants with significantly down-regulated miR396 expression by Agrobacterium-mediated genetic transformation. The MIM396 transgenic alfalfa plants exhibited improved resistance to Spodoptera litura larvae with increased lignin content but decreased JA accumulation. Most of the miR396 putative target GRF genes were up-regulated in MIM396 transgenic lines, and responded to the wounding treatment. By RNA sequencing analysis, we found that the differentially expressed genes related to insect resistance between WT and MIM396 transgenic plants mainly clustered in biosynthesis pathways in lignin, flavonoids and glucosinolates. In addition to the phenotype of enhanced insect resistance, MIM396 transgenic plants also displayed reduced biomass yield and forage quality. Our results broaden the function of miR396 in alfalfa and provide genetic resources for studying alfalfa insect resistance.

A conserved SNP variation in the pre-miR396c flanking region in Oryza sativa indica landraces correlates with mature miRNA abundance

Plant precursor miRNAs (pre-miRNA) have conserved evolutionary footprints that correlate with mode of miRNA biogenesis. In plants, base to loop and loop to base modes of biogenesis have been reported. Conserved structural element(s) in pre-miRNA play a major role in turn over and abundance of mature miRNA. Pre-miR396c sequences and secondary structural characteristics across Oryza species are presented. Based on secondary structure, twelve Oryza pre-miR396c sequences are divided into three groups, with the precursor from halophytic Oryza coarctata forming a distinct group. The miRNA-miRNA* duplex region is completely conserved across eleven Oryza species as are other structural elements in the pre-miRNA, suggestive of an evolutionarily conserved base-to-loop mode of miRNA biogenesis. SNPs within O. coarctata mature miR396c sequence and miRNA* region have the potential to alter target specificity and association with the RNA-induced silencing complex. A conserved SNP variation, rs10234287911 (G/A), identified in O. sativa pre-miR396c sequences alters base pairing above the miRNA-miRNA* duplex. The more stable structure conferred by the 'A_10234287911' allele may promote better processing vis-à-vis the structure conferred by 'G_10234287911' allele. We also examine pri- and pre-miR396c expression in cultivated rice under heat and salinity and their correlation with miR396c expression.

Genome-wide identification and validation of tomato-encoded sRNA as the cross-species antifungal factors targeting the virulence genes of Botrytis cinerea

Recent evidence shows that small RNAs are transferred from a species to another through cross-species transmission and exhibit biological activities in the receptor. In this study, we focused on tomato-derived sRNAs play a role of defense against Botrytis cinerea. Bioinformatics method was firstly employed to identify tomato-encoded sRNAs as the cross-species antifungal factors targeting B. cinerea genes. Then the expression levels of some identifed sRNAs were checked in B. cinerea-infected plant using qRT-PCR method. Exogenic RNA-induced gene silences analysis were performed to investigate the antifungal roles of the sRNAs, and the target genes in B. cinerea of antifungal sRNAs would be confirmed by using co-expression analysis. Results showed that a total of 21 B.cinerea-induced sRNAs with high abundance were identified as the cross-kingdom regulator candidates. Among them, three sRNAs containing a miRNA (miR396a-5p) and two siRNA (siR3 and siR14) were selected for experimental validation and bioassay analysis. qRT-PCR confirmed that all of these 3 sRNAs were induced in tomato leaves by B. cinerea infection. Correspondingly, 4 virulence genes of B. cinerea respectively targeted by these 3 sRNAs were down-regulated. Bioassay revealed that all of these 3 cross-species sRNAs could inhibit the virulence and spore gemination of B. cinerea. Correspondingly, the coding genes of B. cinerea targeted by these sRNAs were also down-regulated. Moreover, the virulence inhibition by double strand sRNA was more effective than that by single strand sRNA. The inhibition efficiency of sRNA against B. cinerea increased with the increase of its concentration. Our findings provide new evidence into the coevolution of pathogens and host plants, as well as new directions for the use of plant-derived sRNAs to control pathogens.

Genome-wide analysis of growth-regulating factor genes in grape (Vitis vinifera L.): identification, characterization and their responsive expression to osmotic stress

Identification, characterization and osmotic stress responsive expression of growth-regulating factor genes in grape. The growth and fruit production of grape vine are severely affected by adverse environmental conditions. Growth-regulating factors (GRFs) play a vital role in the regulation of plant growth, reproduction and stress tolerance. However, their biological functions in fruit vine crops are still largely unknown. In the present study, a total number of nine VvGRFs were identified in the grape genome. Phylogenetic and collinear relationship analysis revealed that they formed seven subfamilies, and have gone through three segmental duplication events. All VvGRFs were predicted to be nucleic localized and contained both the conserved QLQ and WRC domains at their N-terminals, one of the typical structural features of GRF proteins. Quantitative real-time PCR analyses demonstrated that all VvGRFs, with a predominant expression of VvGRF7, were constitutively expressed in roots, leaves and stems of grape plants, and showed responsive expression to osmotic stress. Further growth phenotypic analysis demonstrated that ectopic expression of VvGRF7 promoted the growth and sensitivity of transgenic Arabidopsis plants to osmotic stress. Our findings provide important information for the future study of VvGRF gene functions, and potential gene resources for the genetic breeding of new fruit vine varieties with improved fruit yield and stress tolerance.

Genome‑wide identification and characterization of miR396 family members and their target genes GRF in sorghum (Sorghum bicolor (L.) moench)

MicroRNAs (miRNAs) widely participate in plant growth and development. The miR396 family, one of the most conserved miRNA families, remains poorly understood in sorghum. To reveal the evolution and expression pattern of Sbi-miR396 gene family in sorghum, bioinformatics analysis and target gene prediction were performed on the sequences of the Sbi-miR396 gene family members. The results showed that five Sbi-miR396 members, located on chromosomes 4, 6, and 10, were identified at the whole-genome level. The secondary structure analysis showed that the precursor sequences of all five Sbi-miR396 potentially form a stable secondary stem-loop structure, and the mature miRNA sequences were generated on the 5' arm of the precursors. Sequence analysis identified the mature sequences of the five sbi-miR396 genes were high identity, with differences only at the 1st, 9th and 21st bases at the 5' end. Phylogenetic analysis revealed that Sbi-miR396a, Sbi-miR396b, and Sbi-miR396c were clustered into Group I, and Sbi-miR396d and Sbi-miR396e were clustered into Group II, and all five sbi-miR396 genes were closely related to those of maize and foxtail millet. Expression analysis of different tissue found that Sbi-miR396d/e and Sbi-miR396a/b/c were preferentially and barely expressed, respectively, in leaves, flowers, and panicles. Target gene prediction indicates that the growth-regulating factor family members (SbiGRF1/2/3/4/5/6/7/8/10) were target genes of Sbi-miR396d/e. Thus, Sbi-miR396d/e may affect the growth and development of sorghum by targeting SbiGRFs. In addition, expression analysis of different tissues and developmental stages found that all Sbi-miR396 target genes, SbiGRFs, were barely expressed in leaves, root and shoot, but were predominantly expressed in inflorescence and seed development stage, especially SbiGRF1/5/8. Therefore, inhibition the expression of sbi-miR396d/e may increase the expression of SbiGRF1/5/8, thereby affecting floral organ and seed development in sorghum. These findings provide the basis for studying the expression of the Sbi-mir396 family members and the function of their target genes.

Relationship between reduction in rice (Nipponbare) leaf blade size under elevated CO2 and miR396-GRF module

Elevated CO₂ (eCO₂; 1000 ppm) influences developing rice leaf formation, reducing leaf blade length and width as compared to rice grown under ambient CO₂ (aCO₂; 400 ppm). Since micro RNAs (miRNAs) are known to play multiple roles in plant development, we hypothesized that miRNAs might be involved in modulating leaf size under eCO₂ conditions. To identify miRNAs responding to eCO₂, we profiled miRNA levels in developing rice leaves (P4; plastochron number of the fourth-youngest leaf) under eCO₂ using small RNA-seq. We detected 18 mature miRNA sequences for which expression levels varied more than two-fold between the eCO₂ and aCO₂ conditions. Among them, only miR396e and miR396f significantly differed between the two conditions. Additionally, the expression of growth-regulating factors (GRFs), potential target mRNA of miR396s, were repressed under the eCO₂ condition. We used an antisense oligonucleotide approach to confirm that single-strand DNA corresponding to the miR396e sequence effectively downregulated GRF expression in developing leaves, reducing the leaf blade length, such as for rice grown under eCO₂. These results suggest that the miR396-GRF module is crucially relevant to controlling rice leaf blade length in eCO₂ environments.

Leaf rust responsive miRNA and their target genes in wheat

Small RNA sequencing (sRNA-seq) and degradome analysis were used for the identification of miRNAs and their target host genes in a pair of near-isogenic lines (NILs), which differed for the presence of leaf rust resistance gene Lr28. The study led to identification of (i) 506 known and 346 novel miRNAs; and (ii) 5054 target genes including 4557 in silico predicted and 497 degradome-based genes using 105 differentially expressed (DE) miRNAs. A subset of 128 targets (67 in silico + 61 degradome-based) was differentially expressed in RNA-seq data that was generated by us earlier using the same pair of NILs; among these 128 targets, 58 target genes exhibited an inverse relationship with the DE miRNAs (expression of miRNAs and activation/suppression of target genes). Eight miRNAs which belonged to the conserved miRNA families and were known to be induced in response to fungal diseases in plants included the following: miR156, miR158, miR159, miR168, miR169, miR172, miR319, miR396. The target genes belonged to the following classes of genes known to be involved in downstream disease resistance pathways; peroxidases, sugar transporters, auxin response signaling, oxidation-reduction, etc. It was also noticed that although a majority of miRNAs and target genes followed the above classical inverse relationship, there were also examples, where no such relationship was observed. Among the target genes, there were also 51 genes that were not only regulated by miRNAs, but were also differentially methylated at sequences including the following segments: promotors, introns, TSS, exons. The results of the present study suggest a complex interplay among miRNA genes, target genes, and various epigenetic controls, which regulate the expression of genes involved in downstream pathways for disease resistance.

Fungal Secondary Metabolites and Small RNAs Enhance Pathogenicity during Plant-Fungal Pathogen Interactions

Fungal plant pathogens use proteinaceous effectors as well as newly identified secondary metabolites (SMs) and small non-coding RNA (sRNA) effectors to manipulate the host plant's defense system via diverse plant cell compartments, distinct organelles, and many host genes. However, most molecular studies of plant-fungal interactions have focused on secreted effector proteins without exploring the possibly equivalent functions performed by fungal (SMs) and sRNAs, which are collectively known as "non-proteinaceous effectors". Fungal SMs have been shown to be generated throughout the plant colonization process, particularly in the early biotrophic stages of infection. The fungal repertoire of non-proteinaceous effectors has been broadened by the discovery of fungal sRNAs that specifically target plant genes involved in resistance and defense responses. Many RNAs, particularly sRNAs involved in gene silencing, have been shown to transmit bidirectionally between fungal pathogens and their hosts. However, there are no clear functional approaches to study the role of these SM and sRNA effectors. Undoubtedly, fungal SM and sRNA effectors are now a treasured land to seek. Therefore, understanding the role of fungal SM and sRNA effectors may provide insights into the infection process and identification of the interacting host genes that are targeted by these effectors. This review discusses the role of fungal SMs and sRNAs during plant-fungal interactions. It will also focus on the translocation of sRNA effectors across kingdoms, the application of cross-kingdom RNA interference in managing plant diseases and the tools that can be used to predict and study these non-proteinaceous effectors.

Plant Disease Resistance-Related Signaling Pathways: Recent Progress and Future Prospects

Plant-pathogen interactions induce a signal transmission series that stimulates the plant's host defense system against pathogens and this, in turn, leads to disease resistance responses. Plant innate immunity mainly includes two lines of the defense system, called pathogen-associated molecular pattern-triggered immunity (PTI) and effector-triggered immunity (ETI). There is extensive signal exchange and recognition in the process of triggering the plant immune signaling network. Plant messenger signaling molecules, such as calcium ions, reactive oxygen species, and nitric oxide, and plant hormone signaling molecules, such as salicylic acid, jasmonic acid, and ethylene, play key roles in inducing plant defense responses. In addition, heterotrimeric G proteins, the mitogen-activated protein kinase cascade, and non-coding RNAs (ncRNAs) play important roles in regulating disease resistance and the defense signal transduction network. This paper summarizes the status and progress in plant disease resistance and disease resistance signal transduction pathway research in recent years; discusses the complexities of, and interactions among, defense signal pathways; and forecasts future research prospects to provide new ideas for the prevention and control of plant diseases.

Evolutionary characterization of miR396s in Poaceae exemplified by their genetic effects in wheat and maize

MiR396s play important roles in regulating plant growth and stress response, and great potential for crop yield promotion was anticipated. For more comprehensive and precise understanding of miR396s in Poaceae, we analyzed the phylogenetic linkage, gene expression, and chromosomal distribution of miR396s in this study. Although the mature miR396s' sequences were mostly conserved, differential expression patterns and chromosomal distribution were found among Poaceae species including the major cereal crops rice, wheat, and maize. Consistently, in comparison with rice, wheat and maize plants transformed with the target mimicry construct of miR396 (MIM396) exhibited differential effects on grain size and disease resistance. While the TaMIM396 plants showed increased grain size, panicle length and sensitivity to B. graminis, the ZmMIM396 plants didn't show obvious changes in grain size and disease resistance. In Addition, several GROWTH-REGULATING FACTOR (GRF) genes in wheat and maize were repressed by miR396s, which could be reversed by MIM396, confirming the conserved regulatory roles of miR396 on GRFs. While providing new solution to enhance grain yield in wheat and revealing potential regulatory variations of miR396s in controlling grain size and disease resistance in different crops, this study gives clues to further explore miR396s' functions in other Poaceae species.

Osa-miR160a confers broad-spectrum resistance to fungal and bacterial pathogens in rice

Rice production is threatened by multiple pathogens. Breeding cultivars with broad-spectrum disease resistance is necessary to maintain and improve crop production. Previously we found that overexpression of miR160a enhanced rice blast disease resistance. However, it is unclear whether miR160a also regulates resistance against other pathogens, and what the downstream signaling pathways are. Here, we demonstrate that miR160a positively regulates broad-spectrum resistance against the causative agents of blast, leaf blight and sheath blight in rice. Mutations of miR160a-targeted Auxin Response Factors result in different alteration of resistance conferred by miR160a. miR160a enhances disease resistance partially by suppressing ARF8, as mutation of ARF8 in MIM160 background partially restores the compromised resistance resulting from MIM160. ARF8 protein binds directly to the promoter and suppresses the expression of WRKY45, which acts as a positive regulator of rice immunity. Mutation of WRKY45 compromises the enhanced blast resistance and bacterial leaf blight resistance conferred by arf8 mutant. Overall, our results reveal that a microRNA coordinates rice broad-spectrum disease resistance by suppressing multiple target genes that play different roles in disease resistance, and uncover a new regulatory pathway mediated by the miR160a-ARF8 module. These findings provide new resources to potentially improve disease resistance for breeding in rice.

Analysis of the MIR396 gene family and the role of MIR396b in regulating fiber length in cotton

Cotton fiber is one of the most important natural raw materials in the world textile industry. Improving fiber yield and quality has always been the main goal. MicroRNAs, as typical small noncoding RNAs, could affect fiber length during different stages of fiber development. Based on differentially expressed microRNA in the two interspecific backcross inbred lines (BILs) with a significant difference in fiber length, we identified the miR396 gene family in the two tetraploid cotton genomes and found MIR396b_D13 as the functional precursor to produce mature miR396 during the fiber elongation stage. Among 46 target genes regulated by miR396b, the GROWTH-REGULATING FACTOR 5 gene (GRF5, Gh_A10G0492) had a differential expression level in the two BILs during fiber elongation stage. The expression patterns indicated that the miR396b-GRF5 regulatory module has a critical role in fiber development. Furthermore, virus-induced gene silencing (VIGS) of miR396b significantly produced longer fiber than the wild type, and the expression level of GRF5 showed the reverse trends of the miR396b expression level. The analysis of co-expression network for the GRF5 gene suggested that a cytochrome P450 gene functions as an allene oxide synthase (Gh_D06G0089, AOS), which plays a critical role in jasmonate biosynthetic pathway. In conclusion, our results revealed that the miR396b-GRF5 module has a critical role in fiber development. These findings provide a molecular foundation for fiber quality improvement in the future.

Systematic analysis of the CsmiR396-CsGRFs/CsGIFs module and the opposite role of CsGRF3 and CsGRF5 in regulating cell proliferation in cucumber

Growth-regulating factors (GRFs) are plant-specific transcription factors, and their activities are regulated by miR396 and the GRF-GIF interaction. The miR396-GRFs/GIFs module determines organ size by regulating cell proliferation. However, it is largely unknown in cucumber. In this study, the CsmiR396-CsGRFs/CsGIFs module was investigated in cucumber. Five CsMIR396 loci (CsMIR396A-E), eight CsGRFs and two CsGIFs were identified. CsMIR396A-E was distributed within two clusters and coded three different mature CsmiR396, and all CsGRFs acted as the target of CsmiR396. Bioinformatic analyses showed that miR396s were classified into five types, while GRFs were classified into six groups in plants. The GRFs from group Ⅰ exhibited high diversity and harbored specific characteristics (truncated C-terminus or two WRC domains). qRT-PCR results showed that CsMIR396s (CsMIR396A, CsMIR396B and CsMIR396D) and mature CsmiR396 increased, whereas CsGRFs declined as leaf age increased. In contrast, CsMIR396E was highly expressed in young leaves and shoot tissue, and it was expressed in an age-independent pattern. Yeast two-hybrid assays showed that CsGRF3 strongly interacted with CsGIFs, while CsGRF5 weakly interacted with CsGIFs. Overexpression of CsGRF3 resulted in an enlarged organ size; in contrast, overexpression of CsGRF5, which belonged to group Ⅰ and harbored two WRC domains, resulted in a reduced organ size in Arabidopsis. Section analysis showed that cell proliferation was increased in CsGRF3OE plants, whereas it was decreased in CsGRF5OE plants. In summary, our results reveal the diversity of the CsmiR396-CsGRFs/CsGIFs module in cucumber, and that CsGRF3 and CsGRF5 play an opposite role in regulating cell proliferation.

R gene triplication confers European fodder turnip with improved clubroot resistance

Clubroot is one of the most important diseases for many important cruciferous vegetables and oilseed crops worldwide. Different clubroot resistance (CR) loci have been identified from only limited species in Brassica, making it difficult to compare and utilize these loci. European fodder turnip ECD04 is considered one of the most valuable resources for CR breeding. To explore the genetic and evolutionary basis of CR in ECD04, we sequenced the genome of ECD04 using de novo assembly and identified 978 candidate R genes. Subsequently, the 28 published CR loci were physically mapped to 15 loci in the ECD04 genome, including 62 candidate CR genes. Among them, two CR genes, CRA3.7.1 and CRA8.2.4, were functionally validated. Phylogenetic analysis revealed that CRA3.7.1 and CRA8.2.4 originated from a common ancestor before the whole-genome triplication (WGT) event. In clubroot susceptible Brassica species, CR-gene homologues were affected by transposable element (TE) insertion, resulting in the loss of CR function. It can be concluded that the current functional CR genes in Brassica rapa and non-functional CR genes in other Brassica species were derived from a common ancestral gene before WGT. Finally, a hypothesis for CR gene evolution is proposed for further discussion.

MicroRNA-mediated host defense mechanisms against pathogens and herbivores in rice: balancing gains from genetic resistance with trade-offs to productivity potential

BACKGROUND

Rice (Oryza sativa L.) is the major source of daily caloric intake for more than 30% of the human population. However, the sustained productivity of this staple food crop is continuously threatened by various pathogens and herbivores. Breeding has been successful in utilizing various mechanisms of defense by gene pyramiding in elite cultivars, but the continuous resurgence of highly resistant races of pathogens and herbivores often overcomes the inherent capacity of host plant immunity. MicroRNAs (miRNAs) are endogenous, short, single-stranded, non-coding RNA molecules that regulate gene expression by sequence-specific cleavage of target mRNA or suppressing target mRNA translation. While miRNAs function as upstream regulators of plant growth, development, and host immunity, their direct effects on growth and development in the context of balancing defenses with agronomic potential have not been extensively discussed and explored as a more viable strategy in breeding for disease and pest resistant cultivars of rice with optimal agronomic potentials.

RESULTS

Using the available knowledge in rice and other model plants, this review examines the important roles of miRNAs in regulating host responses to various fungal, bacterial, and viral pathogens, and insect pests, in the context of gains and trade-offs to crop yield. Gains from R-gene-mediated resistance deployed in modern rice cultivars are often undermined by the rapid breakdown of resistance, negative pleiotropic effects, and linkage drags with undesirable traits. In stark contrast, several classes of miRNAs are known to efficiently balance the positive gains from host immunity without significant costs in terms of losses in agronomic potentials (i.e., yield penalty) in rice. Defense-related miRNAs such as Osa-miR156, Osa-miR159, Osa-miR162, Osa-miR396, Osa-530, Osa-miR1432, Osa-miR1871, and Osa-miR1873 are critical in fine-tuning and integrating immune responses with physiological processes that are necessary to the maintenance of grain yield. Recent research has shown that many defense-related miRNAs regulate complex and agronomically important traits.

CONCLUSIONS

Identification of novel immune-responsive miRNAs that orchestrate physiological processes critical to the full expression of agronomic potential will facilitate the stacking of optimal combinations of miRNA-encoding genes to develop high-yielding cultivars with durable resistance to disease and insect pests with minimal penalties to yield.

NAC transcription factors ATAF1 and ANAC055 affect the heat stress response in Arabidopsis

Pre-exposing (priming) plants to mild, non-lethal elevated temperature improves their tolerance to a later higher-temperature stress (triggering stimulus), which is of great ecological importance. ‘Thermomemory’ is maintaining this tolerance for an extended period of time. NAM/ATAF1/2/CUC2 (NAC) proteins are plant-specific transcription factors (TFs) that modulate responses to abiotic stresses, including heat stress (HS). Here, we investigated the potential role of NACs for thermomemory. We determined the expression of 104 Arabidopsis NAC genes after priming and triggering heat stimuli, and found ATAF1 expression is strongly induced right after priming and declines below control levels thereafter during thermorecovery. Knockout mutants of ATAF1 show better thermomemory than wild type, revealing a negative regulatory role. Differential expression analyses of RNA-seq data from ATAF1 overexpressor, ataf1 mutant and wild-type plants after heat priming revealed five genes that might be priming-associated direct targets of ATAF1: AT2G31260 (ATG9), AT2G41640 (GT61), AT3G44990 (XTH31), AT4G27720 and AT3G23540. Based on co-expression analyses applied to the aforementioned RNA-seq profiles, we identified ANAC055 to be transcriptionally co-regulated with ATAF1. Like ataf1, anac055 mutants show improved thermomemory, revealing a potential co-control of both NAC TFs over thermomemory. Our data reveals a core importance of two NAC transcription factors, ATAF1 and ANAC055, for thermomemory.

Genome-Wide Identification of Powdery Mildew Responsive Long Non-Coding RNAs in Cucurbita pepo

Cucurbita pepo L. is an essential economic vegetable crop worldwide, and its production is severely affected by powdery mildew (PM). However, our understanding of the molecular mechanism of PM resistance in C. pepo is very limited. Long non-coding RNAs (lncRNAs) play an important role in regulating plant responses to biotic stress. Here, we systematically identified 2,363 reliably expressed lncRNAs from the leaves of PM-susceptible (PS) and PM-resistant (PR) C. pepo. The C. pepo lncRNAs are shorter in length and expressed at a lower level than the protein-coding transcripts. Among the 2,363 lncRNAs, a total of 113 and 146 PM-responsive lncRNAs were identified in PS and PR, respectively. Six PM-responsive lncRNAs were predicted as potential precursors of microRNAs (miRNAs). In addition, 58 PM-responsive lncRNAs were predicted as targets of miRNAs and one PM-responsive lncRNA was predicted as an endogenous target mimic (eTM). Furthermore, a total of 5,200 potential cis target genes and 5,625 potential trans target genes were predicted for PM-responsive lncRNAs. Functional enrichment analysis showed that these potential target genes are involved in different biological processes, such as the plant-pathogen interaction pathway, MAPK signaling pathway, and plant hormone signal transduction pathway. Taken together, this study provides a comprehensive view of C. pepo lncRNAs and explores the putative functions of PM-responsive lncRNAs, thus laying the foundation for further study of the regulatory mechanisms of lncRNAs responding to PM.

Systematical characterization of GRF gene family in sorghum, and their potential functions in aphid resistance

Sorghum (Sorghum bicolor) is the fifth important cereal and an industrial energy crop in the world. Growth Regulation Factors (GRFs) play an important role in response to environmental stress, however, the knowledge of GRFs relating to the pest resistance is lacking. Here, we identified 8 GRF genes harboring the typical QLQ (glutamine, leucine, glutamine) and WRC (tryptophan, arginine, cysteine) domains in Sorghum, which could be classified into 4 clades through phylogenetic analysis. The SbGRF genes express in most tissues, while more than half of them express at the highest level in inflorescence. To further investigate their possible role in stress response, we analyzed the transcriptomics data. The results showed that SbGRFs could respond to the abiotic stresses including heat, salt and drought stress. Furthermore, combined the data with qRT-PCR, SbGRF1, 2, 4 and 7 were identified as dominant genes response to the aphid-induced stress. SSR markers close to these genes were also searched. Above all, we summarized the SbGRFs and provided their potential roles in aphid response.

Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.)

The NBS-LRR (NLR) gene family plays a pivotal role in regulating disease defense response in plants. Cucumber is one of the most important vegetable crops in the world, and various plant diseases, including powdery mildew (PM), cause severe losses in both cucumber productivity and quality annually. To characterize and understand the role of the CC-NBS-LRR(CNL) family of genes in disease defense response in cucumber plants, we performed bioinformatical analysis to characterize these genes systematically. We identified 33 members of the CNL gene family in cucumber plants, and they are distributed on each chromosome with chromosome 4 harboring the largest cluster of five different genes. The corresponding CNL family member varies in the number of amino acids and exons, molecular weight, theoretical isoelectric point (pI) and subcellular localization. Cis-acting element analysis of the CNL genes reveals the presence of multiple phytohormone, abiotic and biotic responsive elements in their promoters, suggesting that these genes might be responsive to plant hormones and stress. Phylogenetic and synteny analysis indicated that the CNL proteins are conserved evolutionarily in different plant species, and they can be divided into four subfamilies based on their conserved domains. MEME analysis and multiple sequence alignment showed that conserved motifs exist in the sequence of CNLs. Further DNA sequence analysis suggests that CsCNL genes might be subject to the regulation of different miRNAs upon PM infection. By mining available RNA-seq data followed by real-time quantitative PCR (qRT-PCR) analysis, we characterized expression patterns of the CNL genes, and found that those genes exhibit a temporospatial expression pattern, and their expression is also responsive to PM infection, ethylene, salicylic acid, and methyl jasmonate treatment in cucumber plants. Finally, the CNL genes targeted by miRNAs were predicted in cucumber plants. Our results in this study provided some basic information for further study of the functions of the CNL gene family in cucumber plants.

OsAT1, an anion transporter, negatively regulates grain size and yield in rice

Improving the grain yield of rice is a central goal of basic and applied scientific research. Here, we identified an anion transporter, OsAT1, localized in the endoplasmic reticulum and Golgi. OsAT1 is highly expressed in flag, stem, and sheath as monitored using qRT-PCR and pOsAT1::GUS. Thousand-grain weight, grain weight per plant, and content of starch were significantly increased in OsAT1 knock-down mutants (OsAT1-Ri) but significantly decreased in OsAT1 overexpressed lines (OsAT1-OE). In addition, the grain weight per plant increased by 6.17% to 6.78% in OsAT1-RNAi lines, whereas it decreased by 45.93% to 46.76% in OsAT1-OE lines, compared to wild-type. Moreover, the copper content was noticeably reduced in flag leaf of OsAT1-Ri lines and increased in OsAT1-OE lines. RNA-sequencing analysis of OsAT1-OE lines revealed that the genes related to starch biosynthesis and metabolism pathway were enriched in the down-regulated category. Thus, our results suggest that knock-down of OsAT1 in rice possibly reduces copper accumulation and improves the accumulation of storage starch, hence, increasing the grain size and weight. OsAT1 may be a useful gene to consider for cereal breeding programs.

Recent advances and potential applications of cross-kingdom movement of miRNAs in modulating plant's disease response

In the recent past, cross-kingdom movement of miRNAs, small (20–25 bases), and endogenous regulatory RNA molecules has emerged as one of the major research areas to understand the potential implications in modulating the plant’s biotic stress response. The current review discussed the recent developments in the mechanism of cross-kingdom movement (long and short distance) and critical cross-talk between host’s miRNAs in regulating gene function in bacteria, fungi, viruses, insects, and nematodes, and vice-versa during host-pathogen interaction and their potential implications in crop protection. Moreover, cross-kingdom movement during symbiotic interaction, the emerging role of plant’s miRNAs in modulating animal’s gene function, and feasibility of spray-induced gene silencing (SIGS) in combating biotic stresses in plants are also critically evaluated. The current review article analysed the horizontal transfer of miRNAs among plants, animals, and microbes that regulates gene expression in the host or pathogenic organisms, contributing to crop protection. Further, it highlighted the challenges and opportunities to harness the full potential of this emerging approach to mitigate biotic stress efficiently.

Osa-miR535 targets SQUAMOSA promoter binding protein-like 4 to regulate blast disease resistance in rice

Many rice microRNAs have been identified as fine-tuning factors in the regulation of agronomic traits and immunity. Among them, Osa-miR535 targets SQUAMOSA promoter binding protein-like 14 (OsSPL14) to positively regulate tillers but negatively regulate yield and immunity. Here, we uncovered that Osa-miR535 targets another SPL gene, OsSPL4, to suppress rice immunity against Magnaporthe oryzae. Overexpression of Osa-miR535 significantly decreased the accumulation of the fusion protein SPL4_TBS -YFP that contains the target site of Osa-miR535 in OsSPL4. Consistently, Osa-miR535 mediated the cleavage of OsSPL4 mRNA between the 10th and 11th base pair of the predicted binding site at the 3' untranslated region. Transgenic rice lines overexpressing OsSPL4 (OXSPL4) displayed enhanced blast disease resistance accompanied by enhanced immune responses, including increased expression of defense-relative genes and up-accumulated H₂ O₂ . By contrast, the knockout mutant osspl4 exhibited susceptibility. Moreover, OsSPL4 binds to the promoter of GH3.2, an indole-3-acetic acid-amido synthetase, and promotes its expression. Together, these data indicate that Os-miR535 targets OsSPL4 and OsSPL4-GH3.2, which may parallel the OsSPL14-WRKY45 module in rice blast disease resistance.

Blocking Osa-miR1871 enhances rice resistance against Magnaporthe oryzae and yield

MicroRNAs (miRNAs) play vital roles in plant development and defence responses against various stresses. Here, we show that blocking miR1871 improves rice resistance against Magnaporthe oryzae and enhances grain yield simultaneously. The transgenic lines overexpressing miR1871 (OX1871) exhibit compromised resistance, suppressed defence responses and reduced panicle number resulting in slightly decreased yield. In contrast, the transgenic lines blocking miR1871 (MIM1871) show improved resistance, enhanced defence responses and significantly increased panicle number leading to enhanced yield per plant. The RNA-seq assay and defence response assays reveal that blocking miR1871 resulted in the enhancement of PAMP-triggered immunity (PTI). Intriguingly, miR1871 suppresses the expression of LOC_Os06g22850, which encodes a microfibrillar-associated protein (MFAP1) locating nearby the cell wall and positively regulating PTI responses. The mutants of MFAP1 resemble the phenotype of OX1871. Conversely, the transgenic lines overexpressing MFAP1 (OXMFAP1) or overexpressing both MFAP1 and miR1871 (OXMFAP1/OX1871) resemble the resistance of MIM1871. The time-course experiment data reveal that the expression of miR1871 and MFAP1 in rice leaves, panicles and basal internode is dynamic during the whole growth period to manipulate the resistance and yield traits. Our results suggest that miR1871 regulates rice yield and immunity via MFAP1, and the miR8171-MFAP1 module could be used in rice breeding to improve both immunity and yield.

microRNAs and Their Roles in Plant Development

Small RNAs are short non-coding RNAs with a length ranging between 20 and 24 nucleotides. Of these, microRNAs (miRNAs) play a distinct role in plant development. miRNAs control target gene expression at the post-transcriptional level, either through direct cleavage or inhibition of translation. miRNAs participate in nearly all the developmental processes in plants, such as juvenile-to-adult transition, shoot apical meristem development, leaf morphogenesis, floral organ formation, and flowering time determination. This review summarizes the research progress in miRNA-mediated gene regulation and its role in plant development, to provide the basis for further in-depth exploration regarding the function of miRNAs and the elucidation of the molecular mechanism underlying the interaction of miRNAs and other pathways.

Identification and functional deciphering suggested the regulatory roles of long intergenic ncRNAs (lincRNAs) in increasing grafting pepper resistance to Phytophthora capsici

BACKGROUND

As a popular and valuable technique, grafting is widely used to protect against soil-borne diseases and nematodes in vegetable production. Growing evidences have revealed that long intergenic ncRNAs (lincRNAs) are strictly regulated and play essential roles in plants development and stress responses. Nevertheless, genome-wide identification and function deciphering of pepper lincRNAs, especially for their roles in improving grafting pepper resistance to Phytophthora capsici is largely unknown.

RESULTS

In this study, RNA-seq data of grafting and control pepper plants with or without P. capsici inoculation were used to identify lincRNAs. In total, 2,388 reliable lincRNAs were identified. They were relatively longer and contained few exons than protein-coding genes. Similar to coding genes, lincRNAs had higher densities in euchromatin regions; and longer chromosome transcribed more lincRNAs. Expression pattern profiling suggested that lincRNAs commonly had lower expression than mRNAs. Totally, 607 differentially expressed lincRNAs (DE-lincRANs) were identified, of which 172 were found between P. capsici resistance grafting pepper sample GR and susceptible sample LDS. The neighboring genes of DE-lincRNAs and miRNAs competitively sponged by DE-lincRNAs were identified. Subsequently, the expression level of DE-lincRNAs was further confirmed by qRT-PCR and regulation patterns between DE-lincRNAs and neighboring mRNAs were also validated. Function annotation revealed that DE-lincRNAs increased the resistance of grafting prepper to P. capsici by modulating the expression of disease-defense related genes through cis-regulating and/or lincRNA-miRNA-mRNA interaction networks.

CONCLUSIONS

This study identified pepper lincRNAs and suggested their potential roles in increasing the resistance level of grafting pepper to P. capsici.

Identification of miRNA-Target Gene Pairs Responsive to Fusarium Wilt of Cucumber via an Integrated Analysis of miRNA and Transcriptome Profiles

Fusarium wilt (FW) of cucumber (Cucumis sativus L.) caused by Fusarium oxysporum f. sp. cucumerinum (Foc) is a destructive soil-borne disease that severely decreases cucumber yield and quality worldwide. MicroRNAs (miRNAs) are small non-coding RNAs (sRNAs) that are important for regulating host immunity because they affect target gene expression. However, the specific miRNAs and the miRNA/target gene crosstalk involved in cucumber resistance to FW remain unknown. In this study, we compared sRNA-seq and RNA-seq data for cucumber cultivar 'Rijiecheng', which is resistant to FW. The integrated analysis identified FW-responsive miRNAs and their target genes. On the basis of verified expression levels, we detected two highly expressed miRNAs with down-regulated expression in response to Foc. Moreover, an analysis of 21 target genes in cucumber inoculated with Foc indicated that JRL3 (Csa2G362470), which is targeted by miR319a, and BEE1 (Csa1G024150), DAHP1 (Csa2G369040), and PERK2 (Csa4G642480), which are targeted by miR6300, are expressed at high levels, but their expression is further up-regulated after Foc inoculation. These results imply that miR319a-JRL3, miR6300-BEE1, miR6300-DAHP1 and miR6300-PERK2 regulate cucumber defenses against FW, and provide the gene resources that may be useful for breeding programs focused on developing new cucumber varieties with enhanced resistance to FW.

Rice miR1432 Fine-Tunes the Balance of Yield and Blast Disease Resistance via Different Modules

microRNAs act as fine-tuners in the regulation of plant growth and resistance against biotic and abiotic stress. Here we demonstrate that rice miR1432 fine-tunes yield and blast disease resistance via different modules. Overexpression of miR1432 leads to compromised resistance and decreased yield, whereas blocking miR1432 using a target mimic of miR1432 results in enhanced resistance and yield. miR1432 suppresses the expression of LOC_Os03g59790, which encodes an EF-hand family protein 1 (OsEFH1). Overexpression of OsEFH1 leads to enhanced rice resistance but decreased grain yield. Further study revealed that miR1432 and OsEFH1 are differentially responsive to chitin, a fungus-derived pathogen/microbe-associated molecular pattern (PAMP/MAMP). Consistently, blocking miR1432 or overexpression of OsEFH1 improves chitin-triggered immunity responses. In contrast, overexpression of ACOT, another target gene regulating rice yield traits, has no significant effects on rice blast disease resistance. Altogether, these results indicate that miR1432 balances yield and resistance via different target genes, and blocking miR1432 can simultaneously improve yield and resistance.

Identification of miRNAs Contributing to the Broad-Spectrum and Durable Blast Resistance in the Yunnan Local Rice Germplasm

MicroRNAs are 20-24 nucleotide non-coding RNAs and play important roles in plant-environment interactions. In recent years, many microRNAs (miRNAs) have been found to regulate rice immunity against rice blast fungus. However, there are limited studies about miRNAs that directly target resistance (R) genes to regulate rice immunity. In this study, by deep sequencing, small RNA libraries were constructed from four-leaf stage seedlings of the resistant variety Ziyu44 and susceptible variety Jiangnanxiangnuo (JNXN) upon Magnaporthe oryzae infection, we found that much more miRNAs were significantly differentially expressed in Ziyu44 than in JNXN. Among these miRNAs, we focused on miR9664, a newly identified rice miRNA in our sequencing, which was upregulated lightly in Ziyu44 and drastically in JNXN at 24-48 h post-inoculation (hpi). The transgenic plants overexpressing miR9664 (miR9664-oe) displayed reduced defense responses to M. oryzae, while those knocking down miR9664 (miR9664-m) displayed enhanced defense responses to M. oryzae. Most of the detected miR9664 predicted target genes were reduced in the miR9664-oe lines while increased in the miR9664-m lines. The cleavage site of LOC_Os08g07774 was confirmed by RLM-RACE. Meanwhile, after being inoculated with M. oryzae, the genes were expressed differently between Ziyu44 and JNXN. The results suggest that miR9664-mediated R gene turnover contributes to Ziyu44 broad-spectrum resistance to rice blast fungus. Taken together, our research identified a new rice miRNA that directly targets R genes to regulate rice immunity against rice blast fungus, adding significant information to the study of rice-M. oryzae interaction.