Identification and integrated analysis of glyphosate stress-responsive microRNAs, lncRNAs, and mRNAs in rice using genome-wide high-throughput sequencing

Switching action modes of miR408-5p mediates auxin signaling in rice

MicroRNAs (miRNAs) play fundamental roles in many developmental and physiological processes in eukaryotes. MiRNAs in plants generally regulate their targets via either mRNA cleavage or translation repression; however, which approach plays a major role and whether these two function modes can shift remains elusive. Here, we identify a miRNA, miR408-5p that regulates AUXIN/INDOLE ACETIC ACID 30 (IAA30), a critical repressor in the auxin pathway via switching action modes in rice. We find that miR408-5p usually inhibits IAA30 protein translation, but in a high auxin environment, it promotes the decay of IAA30 mRNA when it is overproduced. We further demonstrate that IDEAL PLANT ARCHITECTURE1 (IPA1), an SPL transcription factor regulated by miR156, mediates leaf inclination through association with miR408-5p precursor promoter. We finally show that the miR156-IPA1-miR408-5p-IAA30 module could be controlled by miR393, which silences auxin receptors. Together, our results define an alternative auxin transduction signaling pathway in rice that involves the switching of function modes by miR408-5p, which contributes to a better understanding of the action machinery as well as the cooperative network of miRNAs in plants.

Construction and analysis of the tapping panel dryness-related lncRNA/circRNA-miRNA-mRNA ceRNA network in latex of Hevea brasiliensis

Tapping panel dryness (TPD) results in a severe reduction in latex yield in Hevea brasiliensis. However, the molecular regulatory mechanisms of TPD occurrence are still largely unclear. In this study, whole-transcriptome sequencing was carried out on latex from TPD and healthy trees. In total, 7078 long noncoding RNAs (lncRNAs), 3077 circular RNAs (circRNAs), 4956 miRNAs, and 25041 mRNAs were identified in latex, among which 435 lncRNAs, 68 circRNAs, 320 miRNAs, and 1574 mRNAs were differentially expressed in the latex of TPD trees. GO and KEGG analyses indicated that plant hormone signal transduction, MAPK signaling pathway, and ubiquitin-mediated proteolysis were the key pathways associated with TPD onset. Phytohormone profiling revealed significant changes in the contents of 28 hormonal compounds, among which ACC, ABA, IAA, GA, and JA contents were increased, while SA content was reduced in TPD latex, suggesting that hormone homeostasis is disrupted in TPD trees. Furthermore, we constructed a TPD-related competitive endogenous RNA (ceRNA) regulatory network of lncRNA/circRNA-miRNA-mRNA with 561 edges and 434 nodes (188 lncRNAs, 5 circRNAs, 191 miRNAs, and 50 mRNAs) and identified two hub lncRNAs (MSTRG.11908.1 and MSTRG.8791.1) and four hub miRNAs (hbr-miR156, miR156-x, miRf10477-y, and novel-m0452-3p). Notably, the lncRNA-miR156/157-SPL module containing three hubs probably plays a crucial role in TPD onset. The expression of network hubs and the lncRNA-miR156/157-SPL module were further validated by qRT-PCR. Our results reveal the TPD-associated ceRNA regulatory network of lncRNA/circRNA-miRNA-mRNA in latex and lay a foundation for further investigation of molecular regulatory mechanisms for TPD onset in H. brasiliensis.

The roles of miR156 in abiotic and biotic stresses in plants

MicroRNAs (miRNAs), known as a kind of non-coding RNA, can negatively regulate its target genes. To date, the roles of various miRNAs in plant development and resistance to abiotic and biotic stresses have been widely explored. The present review summarized and discussed the functions of miR156 or miR156-SPL module in abiotic and biotic stresses, such as drought, salt, heat, cold stress, UV-B radiation, heavy mental hazards, nutritional starvation, as well as plant viruses, plant diseases, etc. Based on this, the regulation of miR156-involved stress tolerance was better understood, thus, it would be much easier for plant biologists to carry out suitable strategies to help plants suffer from unfavorable living environments.

Genome-wide identification and in-silico expression analysis of carotenoid cleavage oxygenases gene family in Oryza sativa (rice) in response to abiotic stress

Rice constitutes a foundational cereal and plays a vital role in the culinary sector. However, the detriments of abiotic stress on rice quality and productivity are noteworthy. Carotenoid cleavage oxygenases (CCO) hold vital importance as they enable the particular breakdown of carotenoids and significantly contribute towards the growth and response to abiotic stress in rice. Due to the insufficient information regarding rice CCOs and their potential role in abiotic stress, their utilization in stress-resistant genetic breeding remains limited. The current research identified 16 CCO genes within the Oryza sativa japonica group. These OsCCO genes can be bifurcated into three categories based on their conserved sequences: NCEDs (9-Cis-epoxycarotenoid dioxygenases), CCDs (Carotenoid cleavage dioxygenases) and CCD-like (Carotenoid cleavage dioxygenases-like). Conserved motifs were found in the OsCCO gene sequence via MEME analysis and multiple sequence alignment. Stress-related cis-elements were detected in the promoter regions of OsCCOs genes, indicating their involvement in stress response. Additionally, the promoters of these genes had various components related to plant light, development, and hormone responsiveness, suggesting they may be responsive to plant hormones and involved in developmental processes. MicroRNAs play a pivotal role in the regulation of these 16 genes, underscoring their significance in rice gene regulation. Transcriptome data analysis suggests a tissue-specific expression pattern for rice CCOs. Only OsNCED6 and OsNCED10 significantly up-regulated during salt stress, as per RNA seq analyses. CCD7 and CCD8 levels were also higher in the CCD group during the inflorescence growth stage. This provides insight into the function of rice CCOs in abiotic stress response and identifies possible genes that could be beneficial for stress-resistant breeding.

Role of non-coding RNAs against salinity stress in Oryza species: Strategies and challenges in analyzing miRNAs, tRFs and circRNAs

Salinity is an imbalanced concentration of mineral salts in the soil or water that causes yield loss in salt-sensitive crops. Rice plant is vulnerable to soil salinity stress at seedling and reproductive stages. Different non-coding RNAs (ncRNAs) post-transcriptionally regulate different sets of genes during different developmental stages under varying salinity tolerance levels. While microRNAs (miRNAs) are well known small endogenous ncRNAs, tRNA-derived RNA fragments (tRFs) are an emerging class of small ncRNAs derived from tRNA genes with a demonstrated regulatory role, like miRNAs, in humans but unexplored in plants. Circular RNA (circRNA), another ncRNA produced by back-splicing events, acts as target mimics by preventing miRNAs from binding with their target mRNAs, thereby reducing the miRNA's action upon its target. Same may hold true between circRNAs and tRFs. Hence, the work done on these ncRNAs was reviewed and no reports were found for circRNAs and tRFs under salinity stress in rice, either at seedling or reproductive stages. Even the reports on miRNAs are restricted to seedling stage only, in spite of severe effects on rice crop production due to salt stress during reproductive stage. Moreover, this review sheds light on strategies to predict and analyze these ncRNAs in an effective manner.

Roles of long non-coding RNAs in plant immunity

Robust plant immune systems are fine-tuned by both protein-coding genes and non-coding RNAs. Long non-coding RNAs (lncRNAs) refer to RNAs with a length of more than 200 nt and usually do not have protein-coding function and do not belong to any other well-known non-coding RNA types. The non-protein-coding, low expression, and non-conservative characteristics of lncRNAs restrict their recognition. Although studies of lncRNAs in plants are in the early stage, emerging studies have shown that plants employ lncRNAs to regulate plant immunity. Moreover, in response to stresses, numerous lncRNAs are differentially expressed, which manifests the actions of low-expressed lncRNAs and makes plant-microbe/insect interactions a convenient system to study the functions of lncRNAs. Here, we summarize the current advances in plant lncRNAs, discuss their regulatory effects in different stages of plant immunity, and highlight their roles in diverse plant-microbe/insect interactions. These insights will not only strengthen our understanding of the roles and actions of lncRNAs in plant-microbe/insect interactions but also provide novel insight into plant immune responses and a basis for further research in this field.

Multi-omics analyses reveal the crosstalk between the circadian clock and the response to herbicide application in Oryza sativa

Plants have evolved circadian clock systems that enable biological processes to occur in tandem with periodic changes in the environment. However, it is largely unknown whether crosstalk occurs between the circadian clock and the response to herbicide in rice. We identified 19 conserved rhythmic metabolites which were response to pesticide application and their metabolic abundance peaked mainly at ZT2 or ZT14-ZT18. We found a series of glyphosate, s-Metolachlor, fenclorim, metcamifen and GA3 response genes were expressed following stable circadian rhythms. In order to determine the patterns of their temporal expression, co-expression network analysis was done on 10,467 genes that were periodically expressed throughout a 24-hour period. Next, we identified 4,031 potential direct target genes of OsCCA1 in using DAP-seq data for OsCCA1. Of these, 339, 22, 53, 53 and 63 genes showed a response to glyphosate, s-Metolachlor, fenclorim, metcamifen and GA3 application, respectively. And they were mainly phased from dusk to midnight. Interestingly, we identified significant OsCCA1 binding peaks in the promoter regions of four herbicide resistance genes, including OsCYP81A12, OsCYP81E22, OsCYP76C2, and OsCYP76C4. Finally, we found that herbicide application could affects the expression of some of the central oscillator genes of the rice circadian clock. Here, we used multi-omics data to reveal the crosstalk between the circadian clock and herbicide response processes at the epigenomics, transcriptome, and metabolome levels in rice. This work will serve as a theoretical guide for identifying rhythmic herbicide targets, leading to the creation of new herbicides or the breeding of crops resistant to herbicides.

Non-coding RNA expression analysis revealed the molecular mechanism of flag leaf heterosis in inter-subspecific hybrid rice

Heterosis has been used widespread in agriculture, but its molecular mechanism is inadequately understood. Plants have a large number of non-coding RNAs (ncRNAs), among them, functional ncRNAs that have been studied widely containing long non-coding RNA (lncRNA) and circular RNA (circRNA) that play a role in varied biological processes, as well as microRNA (miRNA), which can not only regulate the post-transcriptional expression of target genes, but also target lncRNA and circRNA then participate the competing endogenous RNA (ceRNA) regulatory network. However, the influence of these three ncRNAs and their regulatory relationships on heterosis is unknown in rice. In this study, the expression profile of ncRNAs and the ncRNA regulatory network related to heterosis were comprehensively analyzed in inter-subspecific hybrid rice. A total of 867 miRNAs, 3,278 lncRNAs and 2,521 circRNAs were identified in the hybrid and its parents. Analysis of the global profiles of these three types of ncRNAs indicated that significant differences existed in the distribution and sequence characteristics of the corresponding genes. The numbers of miRNA and lncRNA in hybrid were higher than those in its parents. A total of 784 ncRNAs (169 miRNAs, 573 lncRNAs and 42 circRNAs) showed differentially expressed in the hybrid, and their target/host genes were vital in stress tolerance, growth and development in rice. These discoveries suggested that the expression plasticity of ncRNA has an important role of inter-subspecific hybrid rice heterosis. It is worth mentioning that miRNAs exhibited substantially more variations between hybrid and parents compared with observed variation for lncRNA and circRNA. Non-additive expression ncRNAs and allele-specific expression genes-related ncRNAs in hybrid were provided in this study, and multiple sets of ncRNA regulatory networks closely related to heterosis were obtained. Meanwhile, heterosis-related regulatory networks of ceRNA (lncRNA and circRNA) and miRNA were also demonstrated.

Prediction of banana maturity based on the sweetness and color values of different segments during ripening

To predict the maturity of bananas, the present study used non-destructive methods to analyze changes in the sweetness and color of the stalks, middles, and tips of bananas during ripening. The results indicated that the respective maturation of these three segments did not occur simultaneously, as indicated by the differential enzyme activity and gene expression levels recorded in these segments. A principal component analysis and cluster plots were used to review the classification of banana maturity, highlighting that banana maturation can be divided into six stages. Two distinct maturity prediction algorithms were established using random forest, artificial neural network, and support vector machines, and they also indicated that dividing the maturity of bananas into six stages was adequate. These findings contribute to the development of quality evaluation and of a rapid grading system for processing, which improves the quality and sale of banana fruits and the related processed products.

•

Sweetness and color during ripening were assessed along banana fingers.

•

A new maturity prediction model was established for bananas.

•

Banana maturity was divided in six stages.

•

The theoretical basis for developing a maturity grading detection device was set.

An Integrated Regulatory Network of mRNAs, microRNAs, and lncRNAs Involved in Nitrogen Metabolism of Moso Bamboo

Nitrogen is a key macronutrient essential for plant growth and development, and its availability has a strong influence on biological processes. Nitrogen fertilizer has been widely applied in bamboo forests in recent decades; however, the mechanism of nitrogen metabolism in bamboo is not fully elucidated. Here, we characterized the morphological, physiological, and transcriptome changes of moso bamboo in response to different schemes for nitrogen addition to illuminate the regulation mechanism of nitrogen metabolism. The appropriate addition of nitrogen improved the chlorophyll content and Pn (net photosynthetic rate) of leaves, the nitrogen and ammonium contents of the seedling roots, the biomass of the whole seedling, the number of lateral roots, and the activity of enzymes involved in nitrogen metabolism in the roots. Based on the whole transcriptome data of the roots, a total of 8,632 differentially expressed mRNAs (DEGs) were identified under different nitrogen additions, such as 52 nitrate transporter genes, 6 nitrate reductase genes, 2 nitrite reductase genes, 2 glutamine synthase genes, 2 glutamate synthase genes (GOGAT), 3 glutamate dehydrogenase genes, and 431 TFs belonging to 23 families. Meanwhile, 123 differentially expressed miRNAs (DEMs) and 396 differentially expressed lncRNAs (DELs) were characterized as nitrogen responsive, respectively. Furthermore, 94 DEM-DEG pairs and 23 DEL-DEG pairs involved in nitrogen metabolism were identified. Finally, a predicted regulatory network of nitrogen metabolism was initially constructed, which included 17 nitrogen metabolic pathway genes, 15 TFs, 4 miRNAs, and 10 lncRNAs by conjoint analysis of DEGs, DEMs, and DELs and their regulatory relationships, which was supported by RNA-seq data and qPCR results. The lncRNA-miRNA-mRNA network provides new insights into the regulation mechanism of nitrogen metabolism in bamboo, which facilitates further genetic improvement for bamboo to adapt to the fluctuating nitrogen environment.

Whole-Transcriptome Analysis Reveals Autophagy Is Involved in Early Senescence of zj-es Mutant Rice

Senescence is a necessary stage of plant growth and development, and the early senescence of rice will lead to yield reduction and quality decline. However, the mechanisms of rice senescence remain obscure. In this study, we characterized an early-senescence rice mutant, designated zj-es (ZheJing-early senescence), which was derived from the japonica rice cultivar Zhejing22. The mutant zj-es exhibited obvious early-senescence phenotype, such as collapsed chloroplast, lesions in leaves, declined fertility, plant dwarf, and decreased agronomic traits. The ZJ-ES gene was mapped in a 458 kb-interval between the molecular markers RM5992 and RM5813 on Chromosome 3, and analysis suggested that ZJ-ES is a novel gene controlling rice early senescence. Subsequently, whole-transcriptome RNA sequencing was performed on zj-es and its wild-type rice to dissect the underlying molecular mechanism for early senescence. Totally, 10,085 differentially expressed mRNAs (DEmRNAs), 1,253 differentially expressed lncRNAs (DElncRNAs), and 614 differentially expressed miRNAs (DEmiRNAs) were identified, respectively, in different comparison groups. Based on the weighted gene co-expression network analysis (WGCNA), the co-expression turquoise module was found to be the key for the occurrence of rice early senescence. Furthermore, analysis on the competing endogenous RNA (CeRNA) network revealed that 14 lncRNAs possibly regulated 16 co-expressed mRNAs through 8 miRNAs, and enrichment analysis showed that most of the DEmRNAs and the targets of DElncRNAs and DEmiRNAs were involved in reactive oxygen species (ROS)-triggered autophagy-related pathways. Further analysis showed that, in zj-es, ROS-related enzyme activities were markedly changed, ROS were largely accumulated, autophagosomes were obviously observed, cell death was significantly detected, and lesions were notably appeared in leaves. Totally, combining our results here and the remaining research, we infer that ROS-triggered autophagy induces the programmed cell death (PCD) and its coupled early senescence in zj-es mutant rice.

miRNAs: A potentially valuable tool in pesticide toxicology assessment-current experimental and epidemiological data review

miRNAs are responsible for the regulation of many cellular processes such as development, cell differentiation, proliferation, apoptosis, and tumor growth. Several studies showed that they can also serve as specific, stable, and sensitive markers of chemical exposure. In this review, current experimental and epidemiological data evidencing deregulation in miRNA expression in response to fungicides, insecticides or herbicides were analyzed. As shown by Venn's diagrams, miR-363 and miR-9 deregulation is associated with fungicide exposure in vitro and in vivo, while let-7, miR-155, miR-181 and miR-21 were found to be commonly deregulated by at least three different insecticides. Furthermore, let-7, miR-30, miR-126, miR-181 and miR-320 were commonly deregulated by 3 different herbicides. Notably, these 5 miRNAs were also found to be deregulated by one or more insecticides, suggesting their participation in the cellular response to pesticides, regardless of their chemical structure. All these miRNAs have been proposed as potential biomarkers for fungicide, insecticide, or herbicide exposure. These results allow us to improve our understanding of the molecular mechanisms of toxicity upon pesticide exposure, although further studies are needed to confirm these miRNAs as definitive (not potential) biomarkers of pesticide exposure.

Integrated analysis of the lncRNA/circRNA-miRNA-mRNA expression profiles reveals novel insights into potential mechanisms in response to root-knot nematodes in peanut

Background

Peanut is the most essential oil and food crop globally due to its high oil and protein content. Root-knot nematode infects peanut roots, causing poor development and severely limiting peanut yields worldwide. The discovery of peanut genome identified a considerable number of genetic loci controlling the peanut root-knot nematode; however, the molecular mechanism of root-knot nematode remains unknown.

Results

The heterogeneous response to root-knot nematode stress in peanut roots was identified using whole-transcriptome RNA-seq. A total of 430 mRNAs, 111 miRNAs, 4453 lncRNAs, and 123 circRNAs were found to have differential expression between infected and non-infected peanuts. The expression profiles of the lncRNA/circRNA-miRNA-mRNA network were developed to understand the potential pathways that lead to root-knot nematodes in peanut roots. During root-knot nematodes stress, a total of 10 lncRNAs, 4 circRNAs, 5 miRNAs, and 13 mRNAs can create competing endogenous RNA and participate in the oxidation–reduction process as well as other biological metabolism processes in peanuts. The findings will highlight the role of peanut ceRNAs in response to root-knot nematodes.

Conclusion

The GO classification and KEGG pathway enrichment study of core regulatory networks revealed that ceRNAs are involved in oxidation–reduction, peroxidase activity, lignin synthesis in the xylem, and flavonoid synthesis. Overall, these findings may help researchers better understand the role of non-coding RNAs in response to root-knot nematodes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08470-3.

Development of 'multiresistance rice' by an assembly of herbicide, insect and disease resistance genes with a transgene stacking system

BACKGROUND

Weeds, diseases and pests pose serious threats to rice production and cause significant economic losses. Cultivation of rice varieties with resistance to herbicides, diseases and pests is believed to be the most economical and environmentally friendly method to deal with these problems.

RESULTS

In this study, a highly efficient transgene stacking system was used to assembly the synthetic glyphosate-tolerance gene (I. variabilis-EPSPS*), lepidopteran pest resistance gene (Cry1C*), brown planthopper resistance genes (Bph14* and OsLecRK1*), bacterial blight resistance gene (Xa23*) and rice blast resistance gene (Pi9*) onto a transformable artificial chromosome vector. The construct was transferred into ZH11 (a widely used japonica rice cultivar Zhonghua 11) via Agrobacterium-mediated transformation and 'multiresistance rice' (MRR) with desirable agronomic traits was obtained. The results showed that MRR had significantly improved resistance to glyphosate, borers, brown planthopper, bacterial blight and rice blast relative to the recipient cultivar ZH11. Besides, under the natural occurrence of pests and diseases in the field, the yield of MRR was significantly higher than that of ZH11.

CONCLUSION

A multigene transformation strategy was employed to successfully develop rice lines with multiresistance to glyphosate, borers, brown planthopper, bacterial blight and rice blast, and the obtained MRR is expected to have great application potential. © 2020 Society of Chemical Industry.

Integrated Analysis of MicroRNA and Target Genes in Brachypodium distachyon Infected by Magnaporthe oryzae by Small RNA and Degradome Sequencing

Rice blast caused by Magnaporthe oryzae is one of the most important diseases that seriously threaten rice production. Brachypodium distachyon is a grass species closely related to grain crops, such as rice, barley, and wheat, and has become a new model plant of Gramineae. In this study, 15 small RNA samples were sequenced to examine the dynamic changes in microRNA (miRNA) expression in B. distachyon infected by M. oryzae at 0, 24, and 48 h after inoculation. We identified 432 conserved miRNAs and 288 predicted candidate miRNAs in B. distachyon. Additionally, there were 7 and 19 differentially expressed miRNAs at 24 and 48 h post-inoculation, respectively. Furthermore, using degradome sequencing, we identified 2,126 genes as targets for 308 miRNAs; using quantitative real-time PCR (qRT-PCR), we validated five miRNA/target regulatory units involved in B. distachyon-M. oryzae interactions. Moreover, using co-transformation technology, we demonstrated that BdNAC21 was negatively regulated by miR164c. This study provides a new approach for identifying resistance genes in B. distachyon by mining the miRNA regulatory network of host-pathogen interactions.