Multiple Regression Analysis Reveals MicroRNA Regulatory Networks in Oryza sativa under Drought Stress

MicroRNAs in Plant Genetic Regulation of Drought Tolerance and Their Function in Enhancing Stress Adaptation

Adverse environmental conditions, including drought stress, pose a significant threat to plant survival and agricultural productivity, necessitating innovative and efficient approaches to enhance their resilience. MicroRNAs (miRNAs) are recognized as key elements in regulating plant adaptation to drought stress, with a notable ability to modulate various physiological and molecular mechanisms. This review provides an in-depth analysis of the role of miRNAs in drought response mechanisms, including abscisic acid (ABA) signaling, reactive oxygen species (ROS) detoxification, and the optimization of root system architecture. Additionally, it examines the effectiveness of bioinformatics tools, such as those employed in in silico analyses, for studying miRNA-mRNA interactions, as well as the potential for their integration with experimental methods. Advanced methods such as microarray analysis, high-throughput sequencing (HTS), and RACE-PCR are discussed for their contributions to miRNA target identification and validation. Moreover, new data and perspectives are presented on the role of miRNAs in plant responses to abiotic stresses, particularly drought adaptation. This review aims to deepen the understanding of genetic regulatory mechanisms in plants and to establish a robust scientific foundation for the development of drought-tolerant crop varieties.

Identification of hub genes and key pathways in arsenic-treated rice (Oryza sativa L.) based on 9 topological analysis methods of CytoHubba

BACKGROUND

Arsenic is a toxic metalloid that can cause acute and chronic adverse health problems. Unfortunately, rice, the primary staple food for more than half of the world's population, is generally regarded as a typical arsenic-accumulating crop plant. Evidence indicates that arsenic stress can influence the growth and development of the rice plant, and lead to high concentrations of arsenic in rice grain. But the underlying mechanisms remain unclear.

METHODS

In the present research, the possible molecules and pathways involved in rice roots in response to arsenic stress were explored using bioinformatics methods. Datasets that involving arsenic-treated rice root and the "study type" that was restricted to "Expression profiling by array" were selected and downloaded from Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between the arsenic-treated group and the control group were obtained using the online web tool GEO2R. Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed to investigate the functions of DEGs. The protein-protein interactions (PPI) network and the molecular complex detection algorithm (MCODE) of DEGs were analyzed using STRING and Cystoscope, respectively. Important nodes and hub genes in the PPI network were predicted and explored using the Cytoscape-cytoHubba plug-in.

RESULTS

Two datasets, GSE25206 and GSE71492, were downloaded from Gene Expression Omnibus (GEO) database. Eighty common DEGs from the two datasets, including sixty-three up-regulated and seventeen down-regulated genes, were then selected. After functional enrichment analysis, these common DEGs were enriched mainly in 10 GO items, including glutathione transferase activity, glutathione metabolic process, toxin catabolic process, and 7 KEGG pathways related to metabolism. After PPI network and MCODE analysis, 49 nodes from the DEGs PPI network were identified, filtering two significant modules. Next, the Cytoscape-cytoHubba plug-in was used to predict important nodes and hub genes. Finally, five genes [Os01g0644000, PRDX6 (Os07g0638400), PRX112 (Os07g0677300), ENO1(Os06g0136600), LOGL9 (Os09g0547500)] were verified and could serve as the best candidates associated with rice root in response to arsenic stress.

CONCLUSIONS

In summary, we elucidated the potential pathways and genes in rice root in response to arsenic stress through a comprehensive bioinformatics analysis.

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.

Plant responses to drought stress: microRNAs in action

Drought is common in most regions of the world, and it has a significant impact on plant growth and development. Plants, on the other hand, have evolved their own defense systems to deal with the extreme weather. The reprogramming of gene expression by microRNAs (miRNAs) is one of these defense mechanisms. miRNAs are short noncoding RNAs that have emerged as key post-transcriptional gene regulators in a variety of species. Drought stress modulates the expression of certain miRNAs that are functionally conserved across plant species. These characteristics imply that miRNA-based genetic changes might improve drought resistance in plants. This study highlights current knowledge of plant miRNA biogenesis, regulatory mechanisms and their role in drought stress responses. miRNAs functions and their adaptations by plants during drought stress has also been explained that can be exploited to promote drought-resistance among economically important crops.

Regulatory role of microRNAs (miRNAs) in the recent development of abiotic stress tolerance of plants

MicroRNAs (miRNAs) are a distinct groups of single-stranded non-coding, tiny regulatory RNAs approximately 20-24 nucleotides in length. miRNAs negatively influence gene expression at the post-transcriptional level and have evolved considerably in the development of abiotic stress tolerance in a number of model plants and economically important crop species. The present review aims to deliver the information on miRNA-mediated regulation of the expression of major genes or Transcription Factors (TFs), as well as genetic and regulatory pathways. Also, the information on adaptive mechanisms involved in plant abiotic stress responses, prediction, and validation of targets, computational tools, and databases available for plant miRNAs, specifically focus on their exploration for engineering abiotic stress tolerance in plants. The regulatory function of miRNAs in plant growth, development, and abiotic stresses consider in this review, which uses high-throughput sequencing (HTS) technologies to generate large-scale libraries of small RNAs (sRNAs) for conventional screening of known and novel abiotic stress-responsive miRNAs adds complexity to regulatory networks in plants. The discoveries of miRNA-mediated tolerance to multiple abiotic stresses, including salinity, drought, cold, heat stress, nutritional deficiency, UV-radiation, oxidative stress, hypoxia, and heavy metal toxicity, are highlighted and discussed in this review.

Cerebrospinal fluid proteome shows disrupted neuronal development in multiple sclerosis

Despite intensive research, the aetiology of multiple sclerosis (MS) remains unknown. Cerebrospinal fluid proteomics has the potential to reveal mechanisms of MS pathogenesis, but analyses must account for disease heterogeneity. We previously reported explorative multivariate analysis by hierarchical clustering of proteomics data of MS patients and controls, which resulted in two groups of individuals. Grouping reflected increased levels of intrathecal inflammatory response proteins and decreased levels of proteins involved in neural development in one group relative to the other group. MS patients and controls were present in both groups. Here we reanalysed these data and we also reanalysed data from an independent cohort of patients diagnosed with clinically isolated syndrome (CIS), who have symptoms of MS without evidence of dissemination in space and/or time. Some, but not all, CIS patients had intrathecal inflammation. The analyses reported here identified a common protein signature of MS/CIS that was not linked to elevated intrathecal inflammation. The signature included low levels of complement proteins, semaphorin-7A, reelin, neural cell adhesion molecules, inter-alpha-trypsin inhibitor heavy chain H2, transforming growth factor beta 1, follistatin-related protein 1, malate dehydrogenase 1 cytoplasmic, plasma retinol-binding protein, biotinidase, and transferrin, all known to play roles in neural development. Low levels of these proteins suggest that MS/CIS patients suffer from abnormally low oxidative capacity that results in disrupted neural development from an early stage of the disease.

Utilization of Transcriptome, Small RNA, and Degradome Sequencing to Provide Insights Into Drought Stress and Rewatering Treatment in Medicago ruthenica

Drought is a major limiting factor in foraging grass yield and quality. Medicago ruthenica (M. ruthenica) is a high-quality forage legume with drought resistance, cold tolerance, and strong adaptability. In this study, we integrated transcriptome, small RNA, and degradome sequencing in identifying drought response genes, microRNAs (miRNAs), and key miRNA-target pairs in M. ruthenica under drought and rewatering treatment conditions. A total of 3,905 genes and 50 miRNAs (45 conserved and 5 novel miRNAs) were significantly differentially expressed in three test conditions (CK: control, DS: plants under drought stress, and RW: plants rewatering after drought stress). The degradome sequencing (AllenScore < 4) analysis revealed that 104 miRNAs (11 novel and 93 conserved miRNAs) were identified with 263 target transcripts, forming 296 miRNA-target pairs in three libraries. There were 38 differentially expressed targets from 16 miRNAs in DS vs. CK, 31 from 11 miRNAs in DS vs. RW, and 6 from 3 miRNAs in RW vs. CK; 21, 18, and 3 miRNA-target gene pairs showed reverse expression patterns in DS vs. CK, DS vs. RW, and RW vs. CK comparison groups, respectively. These findings provide valuable information for further functional characterization of genes and miRNAs in response to abiotic stress, in general, and drought stress in M. ruthenica, and potentially contribute to drought resistance breeding of forage in the future.

Exploring Heat-Response Mechanisms of MicroRNAs Based on Microarray Data of Rice Post-meiosis Panicle

To explore heat response mechanisms of mircoRNAs (miRNAs) in rice post-meiosis panicle, microarray analysis was performed on RNA isolated from rice post-meiosis panicles which were treated at 40°C for 0 min, 10 min, 20 min, 60 min, and 2 h. By integrating paired differentially expressed (DE) miRNAs and mRNA expression profiles, we found that the expression levels of 29 DE-miRNA families were negatively correlated to their 178 DE-target genes. Further analysis showed that the majority of miRNAs in 29 DE-miRNA families resisted the heat stress by downregulating their target genes and a time lag existed between expression of miRNAs and their target genes. Then, GO-Slim classification and functional identification of these 178 target genes showed that (1) miRNAs were mainly involved in a series of basic biological processes even under heat conditions; (2) some miRNAs might play important roles in the heat resistance (such as osa-miR164, osa-miR166, osa-miR169, osa-miR319, osa-miR390, osa-miR395, and osa-miR399); (3) osa-miR172 might play important roles in protecting the rice panicle under the heat stress, but osa-miR437, osa-miR418, osa-miR164, miR156, and miR529 might negatively affect rice fertility and panicle flower; and (4) osa-miR414 might inhibit the flowering gene expression by downregulation of LOC_Os 05g51830 to delay the heading of rice. Finally, a heat-induced miRNA-PPI (protein-protein interaction) network was constructed, and three miRNA coregulatory modules were discovered.

Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

Background

Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant.

Results

Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes.

Conclusions

Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

In silico prediction of human genes as potential targets for rice miRNAs

BACKGROUND

Exogenous microRNAs (miRNAs) enter the human body through food, and their effects on metabolic processes can be considerable. It is important to determine which miRNAs from plants affect the expression of human genes and the extent of their influence.

METHOD

The binding sites of 738Oryza sativa miRNAs (osa-miRNAs) that interact with 17 508 mRNAs of human genes were determined using the MirTarget program.

RESULT

The characteristics of the binding of 46 single osa-miRNAs to 86 mRNAs of human genes with a value of free energy (ΔG) interaction equal 94%-100% from maximum ΔG were established. The findings showed that osa-miR2102-5p, osa-miR5075-3p, osa-miR2097-5p, osa-miR2919 targeted the largest number of genes at 38, 36, 23, 19 sites, respectively. mRNAs of 86 human genes were identified as targets for 93 osa-miRNAs of all family osa-miRNAs with ΔG values equal 94%-98% from maximum ΔG. Each miRNA of the osa-miR156-5p, osa-miR164-5p, osa-miR168-5p, osa-miR395-3p, osa-miR396-3p, osa-miR396-5p, osa-miR444-3p, osa-miR529-3p, osa-miR1846-3p, osa-miR2907-3p families had binding sites in mRNAs of several human target genes. The binding sites of osa-miRNAs in mRNAs of the target genes for each family of osa-miRNAs were conserved when compared to flanking nucleotide sequences.

CONCLUSION

Target mRNA human genes of osa-miRNAs are also candidate genes of cancer, cardiovascular and neurodegenerative diseases.

Genome-Wide Identification and Characterization of Drought Stress Responsive microRNAs in Tibetan Wild Barley

Drought stress is a major obstacle to agricultural production. Tibetan wild barley with rich genetic diversity is useful for drought-tolerant improvement of cereals. MicroRNAs (miRNAs) play critical roles in controlling gene expression in response to various environment perturbations in plants. However, the genome-wide expression profiles of miRNAs and their targets in response to drought stress are largely unknown in wild barley. In this study, a polyethylene glycol (PEG) induced drought stress hydroponic experiment was performed, and the expression profiles of miRNAs from the roots of two contrasting Tibetan wild barley genotypes XZ5 (drought-tolerant) and XZ54 (drought-sensitive), and one cultivated barley Tadmor (drought-tolerant) generated by high-throughput sequencing were compared. There were 69 conserved miRNAs and 1574 novel miRNAs in the dataset of three genotypes under control and drought conditions. Among them, seven conserved miRNAs and 36 novel miRNAs showed significantly genotype-specific expression patterns in response to drought stress. And 12 miRNAs were further regarded as drought tolerant associated miRNAs in XZ5, which mostly participate in gene expression, metabolism, signaling and transportation, suggesting that they and their target genes play important roles in plant drought tolerance. This is the first comparation study on the miRNA transcriptome in the roots of two Tibetan wild barley genotypes differing in drought tolerance and one drought tolerant cultivar in response to PEG treatment. Further results revealed the candidate drought tolerant miRNAs and target genes in the miRNA regulation mechanism in wild barley under drought stress. Our findings provide valuable understandings for the functional characterization of miRNAs in drought tolerance.

Identification of Exogenous Nitric Oxide-Responsive miRNAs from Alfalfa (Medicago sativa L.) under Drought Stress by High-Throughput Sequencing

Alfalfa (Medicago sativa L.) is a high quality leguminous forage. Drought stress is one of the main factors that restrict the development of the alfalfa industry. High-throughput sequencing was used to analyze the microRNA (miRNA) profiles of alfalfa plants treated with CK (normal water), PEG (polyethylene glycol-6000; drought stress), and PEG + SNP (sodium nitroprusside; nitric oxide (NO) sprayed externally under drought stress). We identified 90 known miRNAs belonging to 46 families and predicted 177 new miRNAs. Real-time quantitative fluorescent PCR (qRT-PCR) was used to validate high-throughput expression analysis data. A total of 32 (14 known miRNAs and 18 new miRNAs) and 55 (24 known miRNAs and 31 new miRNAs) differentially expressed miRNAs were identified in PEG and PEG + SNP samples. This suggested that exogenous NO can induce more new miRNAs. The differentially expressed miRNA maturation sequences in the two treatment groups were targeted by 86 and 157 potential target genes, separately. The function of target genes was annotated by gene ontology (GO) enrichment and kyoto encyclopedia of genes and genomes (KEGG) analysis. The expression profiles of nine selected miRNAs and their target genes verified that their expression patterns were opposite. This study has documented that analysis of miRNA under PEG and PEG + SNP conditions provides important insights into the improvement of drought resistance of alfalfa by exogenous NO at the molecular level. This has important scientific value and practical significance for the improvement of plant drought resistance by exogenous NO.