Deep sequencing of wheat sRNA transcriptome reveals distinct temporal expression pattern of miRNAs in response to heat, light and UV

The multifaceted role of RNA-based regulation in plant stress memory

Plants have evolved interconnected regulatory pathways which enable them to respond and adapt to their environments. In plants, stress memory enhances stress tolerance through the molecular retention of prior stressful experiences, fostering rapid and robust responses to subsequent challenges. Mounting evidence suggests a close link between the formation of stress memories and effective future stress responses. However, the mechanism by which environmental stressors trigger stress memory formation is poorly understood. Here, we review the current state of knowledge regarding the RNA-based regulation on stress memory formation in plants and discuss research challenges and future directions. Specifically, we focus on the involvement of microRNAs (miRNAs), small interfering RNAs (siRNAs), long non-coding RNAs (lncRNAs), and alternative splicing (AS) in stress memory formation. miRNAs regulate target genes via post-transcriptional silencing, while siRNAs trigger stress memory formation through RNA-directed DNA methylation (RdDM). lncRNAs guide protein complexes for epigenetic regulation, and AS of pre-mRNAs is crucial to plant stress memory. Unraveling the mechanisms underpinning RNA-mediated stress memory formation not only advances our knowledge of plant biology but also aids in the development of improved stress tolerance in crops, enhancing crop performance and global food security.

Genome-wide expression analysis of novel heat-responsive microRNAs and their targets in contrasting wheat genotypes at reproductive stage under terminal heat stress

Introduction

Heat stress at terminal stage of wheat is critical and leads to huge yield losses worldwide. microRNAs (miRNAs) play significant regulatory roles in gene expression associated with abiotic and biotic stress at the post-transcriptional level.

Methods

In the present study, we carried out a comparative analysis of miRNAs and their targets in flag leaves as well as developing seeds of heat tolerant (RAJ3765) and heat susceptible (HUW510) wheat genotypes under heat stress and normal conditions using small RNA and degradome sequencing.

Results and discussion

A total of 84 conserved miRNAs belonging to 35 miRNA families and 93 novel miRNAs were identified in the 8 libraries. Tae-miR9672a-3p, tae-miR9774, tae-miR9669-5p, and tae-miR5048-5p showed the highest expression under heat stress. Tae-miR9775, tae-miR9662b-3p, tae-miR1120a, tae-miR5084, tae-miR1122a, tae-miR5085, tae-miR1118, tae-miR1130a, tae-miR9678-3p, tae-miR7757-5p, tae-miR9668-5p, tae-miR5050, tae-miR9652-5p, and tae-miR9679-5p were expressed only in the tolerant genotype, indicating their role in heat tolerance. Comparison between heat-treated and control groups revealed that 146 known and 57 novel miRNAs were differentially expressed in the various tissues. Eight degradome libraries sequence identified 457 targets of the differentially expressed miRNAs. Functional analysis of the targets indicated their involvement in photosynthesis, spliceosome, biosynthesis of nucleotide sugars and protein processing in the endoplasmic reticulum, arginine and proline metabolism and endocytosis.

Conclusion

This study increases the number of identified and novel miRNAs along with their roles involved in heat stress response in contrasting genotypes at two developing stages of wheat.

Non-coding RNAs (ncRNAs) in plant: Master regulators for adapting to extreme temperature conditions

Unusual daily temperature fluctuations caused by climate change and climate variability adversely impact agricultural crop production. Since plants are immobile and constantly receive external environmental signals, such as extreme high (heat) and low (cold) temperatures, they have developed complex molecular regulatory mechanisms to cope with stressful situations to sustain their natural growth and development. Among these mechanisms, non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), small-interfering RNAs (siRNAs), and long-non-coding RNAs (lncRNAs), play a significant role in enhancing heat and cold stress tolerance. This review explores the pivotal findings related to miRNAs, siRNAs, and lncRNAs, elucidating how they functionally regulate plant adaptation to extreme temperatures. In addition, this review addresses the challenges associated with uncovering these non-coding RNAs and understanding their roles in orchestrating heat and cold tolerance in plants.

Heat-responsive microRNAs participate in regulating the pollen fertility stability of CMS-D2 restorer line under high-temperature stress

Anther development and pollen fertility of cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) restorer lines are susceptible to continuous high-temperature (HT) stress in summer, which seriously hinders the large-scale application of "three-line" hybrids in production. Here, integrated small RNA, transcriptome, degradome, and hormone profiling was performed to explore the roles of microRNAs (miRNAs) in regulating fertility stability in mature pollens of isonuclear alloplasmic near-isogenic restorer lines NH and SH under HT stress at two environments. A total of 211 known and 248 novel miRNAs were identified, of which 159 were differentially expressed miRNAs (DEMs). Additionally, 45 DEMs in 39 miRNA clusters (PmCs) were also identified, and most highly expressed miRNAs were significantly induced in SH under extreme HT, especially four MIR482 and six MIR6300 family miRNAs. PmC28 was located in the fine-mapped interval of the Rf1 gene and contained two DEMs, gra-miR482_L-2R + 2 and gma-miR2118a-3p_R + 1_1ss18TG. Transcriptome sequencing identified 6281 differentially expressed genes, of which heat shock protein (HSP)-related genes, such as HSP70, HSP22, HSP18.5-C, HSP18.2 and HSP17.3-B, presented significantly reduced expression levels in SH under HT stress. Through integrating multi-omics data, we constructed a comprehensive molecular network of miRNA-mRNA-gene-KEGG containing 35 pairs of miRNA/target genes involved in regulating the pollen development in response to HT, among which the mtr-miR167a_R + 1, tcc-miR167c and ghr-miR390a, tcc-miR396c_L-1 and ghr-MIR169b-p3_1ss6AG regulated the pollen fertility by influencing ARF8 responsible for the auxin signal transduction, ascorbate and aldarate metabolism, and the sugar and lipid metabolism and transport pathways, respectively. Further combination with hormone analysis revealed that HT-induced jasmonic acid signaling could activate the expression of downstream auxin synthesis-related genes and cause excessive auxin accumulation, followed by a cascade of auxin signal transduction, ultimately resulting in pollen abortion. The results provide a new understanding of how heat-responsive miRNAs regulate the stability of fertility restoration for CMS-D2 cotton under heat stress.

High Daytime Temperature Responsive MicroRNA Profiles in Developing Grains of Rice Varieties with Contrasting Chalkiness

High temperature impairs starch biosynthesis in developing rice grains and thereby increases chalkiness, affecting the grain quality. Genome encoded microRNAs (miRNAs) fine-tune target transcript abundances in a spatio-temporal specific manner, and this mode of gene regulation is critical for a myriad of developmental processes as well as stress responses. However, the role of miRNAs in maintaining rice grain quality/chalkiness during high daytime temperature (HDT) stress is relatively unknown. To uncover the role of miRNAs in this process, we used five contrasting rice genotypes (low chalky lines Cyp, Ben, and KB and high chalky lines LaGrue and NB) and compared the miRNA profiles in the R6 stage caryopsis samples from plants subjected to prolonged HDT (from the onset of fertilization through R6 stage of caryopsis development). Our small RNA analysis has identified approximately 744 miRNAs that can be grouped into 291 families. Of these, 186 miRNAs belonging to 103 families are differentially regulated under HDT. Only two miRNAs, Osa-miR444f and Osa-miR1866-5p, were upregulated in all genotypes, implying that the regulations greatly varied between the genotypes. Furthermore, not even a single miRNA was commonly up/down regulated specifically in the three tolerant genotypes. However, three miRNAs (Osa-miR1866-3p, Osa-miR5150-3p and canH-miR9774a,b-3p) were commonly upregulated and onemiRNA (Osa-miR393b-5p) was commonly downregulated specifically in the sensitive genotypes (LaGrue and NB). These observations suggest that few similarities exist within the low chalky or high chalky genotypes, possibly due to high genetic variation. Among the five genotypes used, Cypress and LaGrue are genetically closely related, but exhibit contrasting chalkiness under HDT, and thus, a comparison between them is most relevant. This comparison revealed a general tendency for Cypress to display miRNA regulations that could decrease chalkiness under HDT compared with LaGrue. This study suggests that miRNAs could play an important role in maintaining grain quality in HDT-stressed rice.

Site-specific analysis reveals candidate cross-kingdom small RNAs, tRNA and rRNA fragments, and signs of fungal RNA phasing in the barley-powdery mildew interaction

The establishment of host-microbe interactions requires molecular communication between both partners, which may involve the mutual transfer of noncoding small RNAs. Previous evidence suggests that this is also true for powdery mildew disease in barley, which is caused by the fungal pathogen Blumeria hordei. However, previous studies lacked spatial resolution regarding the accumulation of small RNAs upon host infection by B. hordei. Here, we analysed site-specific small RNA repertoires in the context of the barley-B. hordei interaction. To this end, we dissected infected leaves into separate fractions representing different sites that are key to the pathogenic process: epiphytic fungal mycelium, infected plant epidermis, isolated haustoria, a vesicle-enriched fraction from infected epidermis, and extracellular vesicles. Unexpectedly, we discovered enrichment of specific 31-33-base 5'-terminal fragments of barley 5.8S ribosomal RNA in extracellular vesicles and infected epidermis, as well as particular B. hordei transfer RNA fragments in haustoria. We describe canonical small RNAs from both the plant host and the fungal pathogen that may confer cross-kingdom RNA interference activity. Interestingly, we found first evidence of phased small interfering RNAs in B. hordei, a feature usually attributed to plants, which may be associated with the posttranscriptional control of fungal coding genes, pseudogenes, and transposable elements. Our data suggest a key and possibly site-specific role for cross-kingdom RNA interference and noncoding RNA fragments in the host-pathogen communication between B. hordei and its host barley.

Role of plant microRNAs and their corresponding pathways in fluctuating light conditions

In recent years, it has been established that microRNAs (miRNAs) are critical for various plant physiological regulations in numerous species. Next-generation sequencing technologies have aided to our understandings related to the critical role of miRNAs during environmental stress conditions and plant development. Light influences not just miRNA accumulation but also their biological activities via regulating miRNA gene transcription, biosynthesis, and RNA-induced silencing complex (RISC) activity. Light-regulated routes, processes, and activities can all be affected by miRNAs. Here, we will explore how light affects miRNA gene expression and how conserved and novel miRNAs exhibit altered expression across different plant species in response to variable light quality. Here, we will mainly discuss recent advances in understanding how miRNAs are involved in photomorphogenesis, and photoperiod-dependent plant biological processes such as cell proliferation, metabolism, chlorophyll pigment synthesis and axillary bud growth. The review concludes by presenting future prospects via hoping that light-responsive miRNAs can be exploited in a better way to engineer economically important crops to ensure future food security.

Wheat Long Noncoding RNAs from Organelle and Nuclear Genomes Carry Conserved microRNA Precursors Which May Together Comprise Intricate Networks in Insect Responses

Long noncoding RNAs (lncRNAs) are a diverse class of noncoding RNAs that are typically longer than 200 nucleotides but lack coding potentials. Advances in deep sequencing technologies enabled a better exploration of this type of noncoding transcripts. The poor sequence conservation, however, complicates the identification and annotation of lncRNAs at a large scale. Wheat is among the leading food staples worldwide whose production is threatened by both biotic and abiotic stressors. Here, we identified putative lncRNAs from durum wheat varieties that differ in stem solidness, a major source of defense against wheat stem sawfly, a devastating insect pest. We also analyzed and annotated lncRNAs from two bread wheat varieties, resistant and susceptible to another destructive pest, orange wheat blossom midge, with and without infestation. Several putative lncRNAs contained potential precursor sequences and/or target regions for microRNAs, another type of regulatory noncoding RNAs, which may indicate functional networks. Interestingly, in contrast to lncRNAs themselves, microRNAs with potential precursors within the lncRNA sequences appeared to be highly conserved at the sequence and family levels. We also observed a few putative lncRNAs that have perfect to near-perfect matches to organellar genomes, supporting the recent observations that organellar genomes may contribute to the noncoding transcript pool of the cell.

Comprehensive, integrative genomic analysis of microRNA expression profiles in different tissues of two wheat cultivars with different traits

Wheat is one of the most important food sources on Earth. MicroRNAs (miRNA) play important roles in wheat productivity. To identify wheat miRNAs, we constructed and sequenced sRNA libraries from leaves and roots of two wheat cultivars (RAC875 and Kukri) with many different traits. Given that available miRNA wheat complement in the plant-specific database PmiREN ( https://pmiren.com ) does not include root tissues and root-associated miRNAs might thus be missing, we performed first the prediction of novel miRNAs using the sRNAbench tool. We found a total of 150 putatively novel miRNA genes with expression of both arms from 289 unique mature sequences and nearly 30% of all miRNA reads in roots corresponded to novel miRNAs. In contrast, this figure in leaves dropped to under 3%, confirming the undersampling of roots in the complement of known miRNAs. By using 120 publicly available wheat datasets, 598 Zea mays small RNA libraries, 64 plant species genomes, wheat degradome library, and functional enrichment analysis, a subset of novel miRNAs were confirmed as bona-fide miRNAs. Of the total 605 miRNAs identified in this study inclusive of 316 known miRNAs, 528 miRNAs were shared by both cultivars, 429 miRNAs were shared by both root tissues and 329 miRNAs were shared by both leaf tissues. In addition, 32 miRNAs were specific to Kukri while 45 miRNAs were specific to RAC875. These miRNAs had diverse functions, such as regulation of gene transcription, protein translation, energy metabolism, and cell cycle progression. Our data provide a genome-wide miRNA expression profile in these two wheat cultivars and help functional studies of wheat genomics.

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.

Roles of microRNAs in abiotic stress response and characteristics regulation of plant

MicroRNAs (miRNAs) are a class of non-coding endogenous small RNAs (long 20-24 nucleotides) that negatively regulate eukaryotes gene expression at post-transcriptional level via cleavage or/and translational inhibition of targeting mRNA. Based on the diverse roles of miRNA in regulating eukaryotes gene expression, research on the identification of miRNA target genes has been carried out, and a growing body of research has demonstrated that miRNAs act on target genes and are involved in various biological functions of plants. It has an important influence on plant growth and development, morphogenesis, and stress response. Recent case studies indicate that miRNA-mediated regulation pattern may improve agronomic properties and confer abiotic stress resistance of plants, so as to ensure sustainable agricultural production. In this regard, we focus on the recent updates on miRNAs and their targets involved in responding to abiotic stress including low temperature, high temperature, drought, soil salinity, and heavy metals, as well as plant-growing development. In particular, this review highlights the diverse functions of miRNAs on achieving the desirable agronomic traits in important crops. Herein, the main research strategies of miRNAs involved in abiotic stress resistance and crop traits improvement were summarized. Furthermore, the miRNA-related challenges and future perspectives of plants have been discussed. miRNA-based research lays the foundation for exploring miRNA regulatory mechanism, which aims to provide insights into a potential form of crop improvement and stress resistance breeding.

Arbuscular mycorrhizal fungi improve the growth and drought tolerance of Cinnamomum migao by enhancing physio-biochemical responses

Drought is the main limiting factor for plant growth in karst areas with a fragile ecological environment. Cinnamomum migao H.W. Li is an endemic medicinal woody plant present in the karst areas of southwestern China, and it is endangered due to poor drought tolerance. Arbuscular mycorrhizal fungi (AMF) are known to enhance the drought tolerance of plants. However, few studies have examined the contribution of AMF in improving the drought tolerance of C. migao seedlings. Therefore, we conducted a series of experiments to determine whether a single inoculation and coinoculation of AMF (Claroideoglomus lamellosum and Claroideoglomus etunicatum) enhanced the drought tolerance of C. migao. Furthermore, we compared the effects of single inoculation and coinoculation with different inoculum sizes (20, 40, 60, and 100 g; four replicates per treatment) on mycorrhizal colonization rate, plant growth, photosynthetic parameters, antioxidant enzyme activity, and malondialdehyde (MDA) and osmoregulatory substance contents. The results showed that compared with nonmycorrhizal plants, AMF colonization significantly improved plant growing status; net photosynthetic rate; superoxide dismutase, catalase, and peroxidase activities; and soluble sugar, soluble protein, and proline contents. Furthermore, AMF colonization increased relative water content and reduced MDA content in cells. These combined cumulative effects of AMF symbiosis ultimately enhanced the drought tolerance of seedlings and were closely related to the inoculum size. With an increase in inoculum size, the growth rate and drought tolerance of plants first increased and then decreased. The damage caused by drought stress could be reduced by inoculating 40-60 g of AMF, and the effect of coinoculation was significantly better than that of single inoculation at 60 g of AMF, while the effect was opposite at 40 g of AMF. Additionally, the interaction between AMF and inoculum sizes had a significant effect on drought tolerance. In conclusion, the inoculation of the AMF (Cl. lamellosum and Cl. etunicatum) improved photosynthesis, activated antioxidant enzymes, regulated cell osmotic state, and enhanced the drought tolerance of C. migao, enabling its growth in fragile ecological environments.

The genetic and molecular basis for improving heat stress tolerance in wheat

Wheat production requires at least ~ 2.4% increase per year rate by 2050 globally to meet food demands. However, heat stress results in serious yield loss of wheat worldwide. Correspondingly, wheat has evolved genetic basis and molecular mechanisms to protect themselves from heat-induced damage. Thus, it is very urgent to understand the underlying genetic basis and molecular mechanisms responsive to elevated temperatures to provide important strategies for heat-tolerant varieties breeding. In this review, we focused on the impact of heat stress on morphology variation at adult stage in wheat breeding programs. We also summarize the recent studies of genetic and molecular factors regulating heat tolerance, including identification of heat stress tolerance related QTLs/genes, and the regulation pathway in response to heat stress. In addition, we discuss the potential ways to improve heat tolerance by developing new technologies such as genome editing. This review of wheat responses to heat stress may shed light on the understanding heat-responsive mechanisms, although the regulatory network of heat tolerance is still ambiguous in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-021-00064-z.

Physiological and molecular insights on wheat responses to heat stress

Increasing temperature is a key component of global climate change, affecting crop growth and productivity worldwide. Wheat is a major cereal crop grown in various parts of the globe, which is affected severely by heat stress. The morphological parameters affected include germination, seedling establishment, source-sink activity, leaf area, shoot and root growth. The physiological parameters such as photosynthesis, respiration, leaf senescence, water and nutrient relation are also affected by heat. At the cellular level, heat stress leads to the generation of reactive oxygen species that disrupt the membrane system of thylakoid, chloroplast and plasma membrane. The deactivation of the photosystem, reduction in photosynthesis and inactivation of rubisco affect the production of photoassimilates and their allocation. This ultimately affects anthesis, grain filling, size, number and maturity of wheat grains, which hamper crop productivity. The interplay of various systems comprising antioxidants and hormones plays a crucial role in imparting heat stress tolerance in wheat. Thus, implementation of various omics technologies could foster in-depth insights on heat stress effects, eventually devising heat stress mitigation strategies by conventional and modern breeding to develop heat-tolerant wheat varieties. This review provides an integrative view of heat stress responses in wheat and also discusses approaches to develop heat-tolerant wheat varieties.

Small RNAs: The Essential Regulators in Plant Thermotolerance

Small RNAs (sRNAs) are a class of non-coding RNAs that consist of 21-24 nucleotides. They have been extensively investigated as critical regulators in a variety of biological processes in plants. sRNAs include two major classes: microRNAs (miRNAs) and small interfering RNAs (siRNAs), which differ in their biogenesis and functional pathways. Due to global warming, high-temperature stress has become one of the primary causes for crop loss worldwide. Recent studies have shown that sRNAs are involved in heat stress responses in plants and play essential roles in high-temperature acclimation. Genome-wide studies for heat-responsive sRNAs have been conducted in many plant species using high-throughput sequencing. The roles for these sRNAs in heat stress response were also unraveled subsequently in model plants and crops. Exploring how sRNAs regulate gene expression and their regulatory mechanisms will broaden our understanding of sRNAs in thermal stress responses of plant. Here, we highlight the roles of currently known miRNAs and siRNAs in heat stress responses and acclimation of plants. We also discuss the regulatory mechanisms of sRNAs and their targets that are responsive to heat stress, which will provide powerful molecular biological resources for engineering crops with improved thermotolerance.

Discovery of miRNAs and Development of Heat-Responsive miRNA-SSR Markers for Characterization of Wheat Germplasm for Terminal Heat Tolerance Breeding

A large proportion of the Asian population fulfills their energy requirements from wheat (Triticum aestivum L.). Wheat quality and yield are critically affected by the terminal heat stress across the globe. It affects approximately 40% of the wheat-cultivating regions of the world. Therefore, there is a critical need to develop improved terminal heat-tolerant wheat varieties. Marker-assisted breeding with genic simple sequence repeats (SSR) markers have been used for developing terminal heat-tolerant wheat varieties; however, only few studies involved the use of microRNA (miRNA)-based SSR markers (miRNA-SSRs) in wheat, which were found as key players in various abiotic stresses. In the present study, we identified 104 heat-stress-responsive miRNAs reported in various crops. Out of these, 70 miRNA-SSR markers have been validated on a set of 20 terminal heat-tolerant and heat-susceptible wheat genotypes. Among these, only 19 miRNA-SSR markers were found to be polymorphic, which were further used to study the genetic diversity and population structure. The polymorphic miRNA-SSRs amplified 61 SSR loci with an average of 2.9 alleles per locus. The polymorphic information content (PIC) value of polymorphic miRNA-SSRs ranged from 0.10 to 0.87 with a mean value of 0.48. The dendrogram constructed using unweighted neighbor-joining method and population structure analysis clustered these 20 wheat genotypes into 3 clusters. The target genes of these miRNAs are involved either directly or indirectly in providing tolerance to heat stress. Furthermore, two polymorphic markers miR159c and miR165b were declared as very promising diagnostic markers, since these markers showed specific alleles and discriminated terminal heat-tolerant genotypes from the susceptible genotypes. Thus, these identified miRNA-SSR markers will prove useful in the characterization of wheat germplasm through the study of genetic diversity and population structural analysis and in wheat molecular breeding programs aimed at terminal heat tolerance of wheat varieties.

Small RNAs, Degradome, and Transcriptome Sequencing Provide Insights into Papaya Fruit Ripening Regulated by 1-MCP

As an inhibitor of ethylene receptors, 1-methylcyclopropene (1-MCP) can delay the ripening of papaya. However, improper 1-MCP treatment will cause a rubbery texture in papaya. Understanding of the underlying mechanism is still lacking. In the present work, a comparative sRNA analysis was conducted after different 1-MCP treatments and identified a total of 213 miRNAs, of which 44 were known miRNAs and 169 were novel miRNAs in papaya. Comprehensive functional enrichment analysis indicated that plant hormone signal pathways play an important role in fruit ripening. Through the comparative analysis of sRNAs and transcriptome sequencing, a total of 11 miRNAs and 12 target genes were associated with the ethylene and auxin signaling pathways. A total of 1741 target genes of miRNAs were identified by degradome sequencing, and nine miRNAs and eight miRNAs were differentially expressed under the ethylene and auxin signaling pathways, respectively. The network regulation diagram of miRNAs and target genes during fruit ripening was drawn. The expression of 11 miRNAs and 12 target genes was verified by RT-qPCR. The target gene verification showed that cpa-miR390a and cpa-miR396 target CpARF19-like and CpERF RAP2-12-like, respectively, affecting the ethylene and auxin signaling pathways and, therefore, papaya ripening.

Small RNAs and their targets are associated with the transgenerational effects of water-deficit stress in durum wheat

Water-deficit stress negatively affects wheat yield and quality. Abiotic stress on parental plants during reproduction may have transgenerational effects on progeny. Here we investigated the transgenerational influence of pre-anthesis water-deficit stress by detailed analysis of the yield components, grain quality traits, and physiological traits in durum wheat. Next-generation sequencing analysis profiled the small RNA-omics, mRNA transcriptomics, and mRNA degradomics in first generation progeny. Parental water-deficit stress had positive impacts on the progeny for traits including harvest index and protein content in the less stress-tolerant variety. Small RNA-seq identified 1739 conserved and 774 novel microRNAs (miRNAs). Transcriptome-seq characterised the expression of 66,559 genes while degradome-seq profiled the miRNA-guided mRNA cleavage dynamics. Differentially expressed miRNAs and genes were identified, with significant regulatory patterns subject to trans- and inter-generational stress. Integrated analysis using three omics platforms revealed significant biological interactions between stress-responsive miRNA and targets, with transgenerational stress tolerance potentially contributed via pathways such as hormone signalling and nutrient metabolism. Our study provides the first confirmation of the transgenerational effects of water-deficit stress in durum wheat. New insights gained at the molecular level indicate that key miRNA-mRNA modules are candidates for transgenerational stress improvement.

MiR396 regulatory network and its expression during grain development in wheat

Wheat contains the largest number of miR396 family with 17 miR396 in Poaceae. MiR396 regulatory network underlying wheat grain development has not comprehensively been explored. Our results showed that precursor miR396 family in Poaceae exhibited not only conservativeness but also diversification especially in wheat. Five haplotypes were detected in Poaceae species, while 4 haplotypes in wheat with Hap-4 (miR396a) and Hap-5 (miR396n) unique to wheat. GO enrichment analysis of target genes showed that the first 20 enrichment functions of miR396a and miR396n are completely different from each other, and also completely different from miR396(b-g), miR396(h-m), and miR396(o-q). Functional annotation on the 18 target genes shared by miR396(b-g), miR396(h-m), and miR396(o-q) found that 11 of the 18 target genes are growth-regulating factor (GRF) genes. Our results indicated that, during the grain filling stage of wheat, miR396 is involved in the development of grains by regulating the expression of GRF genes (GRF1, GRF6, and GRF9). Although the enrichment function of miR396(b-g), miR396(h-m), and miR396(o-q) is the same, the gene functional networks they formed differ greatly. Our results indicated that polyploidization enriches not only the diversity of miR396 family and its target genes but also gene functional networks in wheat. These results laid foundation for further elucidating function of miR396 gene family underlying wheat grain development.

Integrated Analysis of Small RNA, Transcriptome, and Degradome Sequencing Reveals the Water-Deficit and Heat Stress Response Network in Durum Wheat

Water-deficit and heat stress negatively impact crop production. Mechanisms underlying the response of durum wheat to such stresses are not well understood. With the new durum wheat genome assembly, we conducted the first multi-omics analysis with next-generation sequencing, providing a comprehensive description of the durum wheat small RNAome (sRNAome), mRNA transcriptome, and degradome. Single and combined water-deficit and heat stress were applied to stress-tolerant and -sensitive Australian genotypes to study their response at multiple time-points during reproduction. Analysis of 120 sRNA libraries identified 523 microRNAs (miRNAs), of which 55 were novel. Differentially expressed miRNAs (DEMs) were identified that had significantly altered expression subject to stress type, genotype, and time-point. Transcriptome sequencing identified 49,436 genes, with differentially expressed genes (DEGs) linked to processes associated with hormone homeostasis, photosynthesis, and signaling. With the first durum wheat degradome report, over 100,000 transcript target sites were characterized, and new miRNA-mRNA regulatory pairs were discovered. Integrated omics analysis identified key miRNA-mRNA modules (particularly, novel pairs of miRNAs and transcription factors) with antagonistic regulatory patterns subject to different stresses. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis revealed significant roles in plant growth and stress adaptation. Our research provides novel and fundamental knowledge, at the whole-genome level, for transcriptional and post-transcriptional stress regulation in durum wheat.

Different MicroRNA Families Involved in Regulating High Temperature Stress Response during Cotton (Gossypium hirsutum L.) Anther Development

MicroRNAs (miRNAs) are small molecule RNAs widely involved in responses to plant abiotic stresses. We performed small RNA sequencing of cotton anthers at four developmental stages under normal and high temperature (NT and HT, respectively) conditions to investigate the stress response characteristics of miRNA to HT. A total of 77 miRNAs, including 33 known miRNAs and 44 novel miRNAs, were identified, and 41 and 28 miRNAs were differentially expressed under NT and HT stress conditions, respectively. The sporogenous cell proliferation (SCP), meiotic phase (MP), microspore release period (MRP), and pollen maturity (PM) stages had 10 (including 12 miRNAs), four (including six miRNAs), four (including five miRNAs), and seven (including 11 miRNAs) HT stress-responsive miRNA families, respectively, which were identified after removing the changes in genotype-specific miRNAs under NT condition. Seven miRNA families (miR2949, miR167, and miR160 at the SCP stage; miR156 and miR172 at the MP stage; miR156 at the MRP stage; and miR393 and miR3476 at the PM stage), which had expression abundance of more than 10% of the total expression abundance, served as the main regulators responding to HT stress with positive or negative regulation patterns. These miRNAs orchestrated the expression of the corresponding target genes and led to different responses in the HT-tolerant and the HT-sensitive lines. The results revealed that the HT stress response of miRNAs in cotton anthers were stage-specific and differed with the development of anthers. Our study may enhance the understanding of the response of miRNAs to HT stress in cotton anthers and may clarify the mechanism of plant tolerance to HT stress.

Transgenerational Effects of Water-Deficit and Heat Stress on Germination and Seedling Vigour-New Insights from Durum Wheat microRNAs

Water deficiency and heat stress can severely limit crop production and quality. Stress imposed on the parents during reproduction could have transgenerational effects on their progeny. Seeds with different origins can vary significantly in their germination and early growth. Here, we investigated how water-deficit and heat stress on parental durum wheat plants affected seedling establishment of the subsequent generation. One stress-tolerant and one stress-sensitive Australian durum genotype were used. Seeds were collected from parents with or without exposure to stress during reproduction. Generally, stress on the previous generation negatively affected seed germination and seedling vigour, but to a lesser extent in the tolerant variety. Small RNA sequencing utilising the new durum genome assembly revealed significant differences in microRNA (miRNA) expression in the two genotypes. A bioinformatics approach was used to identify multiple miRNA targets which have critical molecular functions in stress adaptation and plant development and could therefore contribute to the phenotypic differences observed. Our data provide the first confirmation of the transgenerational effects of reproductive-stage stress on germination and seedling establishment in durum wheat. New insights gained on the epigenetic level indicate that durum miRNAs could be key factors in optimising seed vigour for breeding superior germplasm and/or varieties.

Identification and evolutionary analysis of polycistronic miRNA clusters in domesticated and wild wheat

MicroRNAs are ~22 nucleotide long non-coding RNAs that regulate gene expression at posttranscriptional level. Genome-wide analysis was performed to identify polycistronic miRNAs from wheat. Total 89 polycistronic miRNAs were identified in bread wheat which were distributed on three component sub-genomes (A = 26, B = 33 and D = 30). Except some, most of the identified polycistronic miRNAs were also present in other cultivated and wild wheat species. Expression of 11 identified polycistronic miRNAs could be validated using previously assembled transcriptomes, RNA-seq/s-RNA seq data of cultivated and wild wheats and RT-PCR. Polycistronic miRNAs orthologs were also localized on rice and Brachypodium genomes. As a case study, we also analyzed molecular evolution of miR395 family polycistrons in wheat. Both tandem and segmental duplications contributed to expansion of miR395 family polycistrons. Our findings provide a comprehensive view on wheat polycitronic miRNAs that will enable their in-depth functional analysis in the future.

MicroRNA528, a hub regulator modulating ROS homeostasis via targeting of a diverse set of genes encoding copper-containing proteins in monocots

Plant microRNAs (miRNAs) regulate vital cellular processes, including responses to extreme temperatures with which reactive oxygen species (ROS) are often closely associated. In the present study, it was found that aberrant temperatures caused extensive changes in abundance to numerous miRNAs in banana fruit, especially the copper (Cu)-associated miRNAs. Among them, miR528 was significantly downregulated under cold stress and it was found to target genes encoding polyphenol oxidase (PPO), different from those identified in rice and maize. Expression of PPO genes was upregulated by > 100-fold in cold conditions, leading to ROS surge and subsequent peel browning of banana fruit. Extensive comparative genomic analyses revealed that the monocot-specific miR528 can potentially target a large collection of genes encoding Cu-containing proteins. Most of them are actively involved in cellular ROS metabolism, including not only ROS generating oxidases, but also ROS scavenging enzymes. It also was demonstrated that miR528 has evolved a distinct preference of target genes in different monocots, with its target site varying in position among/within gene families, implying a highly dynamic process of target gene diversification. Its broad capacity to target genes encoding Cu-containing protein implicates miR528 as a key regulator for modulating the cellular ROS homeostasis in monocots.

Identification of Fusarium graminearum-responsive miRNAs and their targets in wheat by sRNA sequencing and degradome analysis

Fusarium head blight (FHB), a prevalent disease of bread wheat (Triticum aestivum L.) caused by Fusarium graminearum, leads to considerable losses of yield and quality in wheat production. MicroRNAs (miRNAs) are important regulators of plant defense responses. Here, to better understand the F. graminearum-responsive miRNAs, we constructed sRNA libraries for wheat cultivar Sumai 3 challenged with F. graminearum and sterile water, respectively. As a result, a total of 203 known miRNAs from 46 families and 68 novel miRNAs were identified. Among them, 18 known and six novel miRNAs were found to be differentially expressed between the F. graminearum-infected samples and the controls and thus were considered to be responsive to F. graminearum. The expression patterns of eight miRNAs were further validated by stem-loop qRT-PCR. Meanwhile, target genes were validated by degradome sequencing. Integrative analysis of the differentially expressed miRNAs and their targets revealed complex miRNA-mediated regulatory networks involved in the response of wheat to F. graminearum infection. Our findings are expected to facilitate a better understanding of the miRNA regulation in wheat-F. graminearum interaction.

Lineage-Specific Evolved MicroRNAs Regulating NB-LRR Defense Genes in Triticeae

Disease resistance genes encoding proteins with nucleotide binding sites and Leucine-Rich Repeat (NB-LRR) domains include many members involved in the effector-triggered immunity pathway in plants. The transcript levels of these defense genes are negatively regulated by diverse microRNAs (miRNAs) in angiosperms and gymnosperms. In wheat, using small RNA expression datasets and degradome datasets, we identified five miRNA families targeting NB-LRR defense genes in monocots, some of which arose in the Triticeae species era. These miRNAs regulate different types of NB-LRR genes, most of them with coil-coiled domains, and trigger the generation of secondary small interfering RNAs (siRNA) as a phased pattern in the target site regions. In addition to acting in response to biotic stresses, they are also responsive to abiotic stresses such as heat, drought, salt, and light stress. Their copy number and expression variation in Triticeae suggest a rapid birth and death frequency. Altogether, non-conserved miRNAs as conserved transcriptional regulators in gymnosperms and angiosperms regulating the disease resistance genes displayed quick plasticity including the variations of sequences, gene copy number, functions, and expression level, which accompanied with NB-LRR genes may be tune-regulated to plants in natural environments with various biotic and abiotic stresses.

MicroRNA-guided regulation of heat stress response in wheat

Background

With rising global temperature, understanding plants’ adaptation to heat stress has implications in plant breeding. MicroRNAs (miRNAs) are small, non-coding, regulatory RNAs guiding gene expression at the post-transcriptional level. In this study, small RNAs and the degradome (parallel analysis of RNA ends) of leaf tissues collected from control and heat-stressed wheat plants immediately at the end of the stress period and 1 and 4 days later were analysed.

Results

Sequencing of 24 small RNA libraries produced 55.2 M reads while 404 M reads were obtained from the corresponding 24 PARE libraries. From these, 202 miRNAs were ascertained, of which mature miRNA evidence was obtained for 104 and 36 were found to be differentially expressed after heat stress. The PARE analysis identified 589 transcripts targeted by 84 of the ascertained miRNAs. PARE sequencing validated the targets of the conserved members of miRNA156, miR166 and miR393 families as squamosa promoter-binding-like, homeobox leucine-zipper and transport inhibitor responsive proteins, respectively. Heat stress responsive miRNA targeted superoxide dismutases and an array of homeobox leucine-zipper proteins, F-box proteins and protein kinases. Query of miRNA targets to interactome databases revealed a predominant association of stress responses such as signalling, antioxidant activity and ubiquitination to superoxide dismutases, F-box proteins, pentatricopeptide repeat-containing proteins and mitochondrial transcription termination factor-like proteins.

Conclusion

The interlaced data set generated in this study identified and validated heat stress regulated miRNAs and their target genes associated with thermotolerance. Such accurate identification and validation of miRNAs and their target genes are essential to develop novel regulatory gene-based breeding strategies.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5799-6) contains supplementary material, which is available to authorized users.

OsmiR528 Enhances Cold Stress Tolerance by Repressing Expression of Stress Response-related Transcription Factor Genes in Plant Cells

Background:

MicroRNAs participate in many molecular mechanisms and signaling trans-duction pathways that are associated with plant stress tolerance by repressing expression of their target genes. However, how microRNAs enhance tolerance to low temperature stress in plant cells remains elusive.

Objective:

In this investigation, we demonstrated that overexpression of the rice microRNA528 (Os-miR528) increases cell viability, growth rate, antioxidants content, ascorbate peroxidase (APOX) activi-ty, and superoxide dismutase (SOD) activity and decreases ion leakage rate and thiobarbituric acid reac-tive substances (TBARS) under low temperature stress in Arabidopsis (Arabidopsis thaliana), pine (Pi-nus elliottii), and rice (Oryza sativa).

Methods:

To investigate the potential mechanism of OsmiR528 in increasing cold stress tolerance, we examined expression of stress-associated MYB transcription factors OsGAMYB-like1, OsMYBS3, OsMYB4, OsMYB3R-2, OsMYB5, OsMYB59, OsMYB30, OsMYB1R, and OsMYB20 in rice cells by qRT-PCR.

Results:

Our experiments demonstrated that OsmiR528 decreases expression of transcription factor OsMYB30 by targeting a F-box domain containing protein gene (Os06g06050), which is a positive regulator of OsMYB30. In OsmiR528 transgenic rice, reduced OsMYB30 expression results in in-creased expression of BMY genes OsBMY2, OsBMY6, and OsBMY10. The transcript levels of the OsBMY2, OsBMY6, and OsBMY10 were elevated by OsMYB30 knockdown, but decreased by Os-MYB30 overexpression in OsmiR528 transgenic cell lines, suggesting that OsmiR528 increases low temperature tolerance by modulating expression of stress response-related transcription factor.

Conclusion:

Our experiments provide novel information in increasing our understanding in molecular mechanisms of microRNAs-associated low temperature tolerance and are valuable in plant molecular breeding from monocotyledonous, dicotyledonous, and gymnosperm plants.

Facing Climate Change: Biotechnology of Iconic Mediterranean Woody Crops

The Mediterranean basin is especially sensitive to the adverse outcomes of climate change and especially to variations in rainfall patterns and the incidence of extremely high temperatures. These two concurring adverse environmental conditions will surely have a detrimental effect on crop performance and productivity that will be particularly severe on woody crops such as citrus, olive and grapevine that define the backbone of traditional Mediterranean agriculture. These woody species have been traditionally selected for traits such as improved fruit yield and quality or alteration in harvesting periods, leaving out traits related to plant field performance. This is currently a crucial aspect due to the progressive and imminent effects of global climate change. Although complete genome sequence exists for sweet orange (Citrus sinensis) and clementine (Citrus clementina), olive tree (Olea europaea) and grapevine (Vitis vinifera), the development of biotechnological tools to improve stress tolerance still relies on the study of the available genetic resources including interspecific hybrids, naturally occurring (or induced) polyploids and wild relatives under field conditions. To this respect, post-genomic era studies including transcriptomics, metabolomics and proteomics provide a wide and unbiased view of plant physiology and biochemistry under adverse environmental conditions that, along with high-throughput phenotyping, could contribute to the characterization of plant genotypes exhibiting physiological and/or genetic traits that are correlated to abiotic stress tolerance. The ultimate goal of precision agriculture is to improve crop productivity, in terms of yield and quality, making a sustainable use of land and water resources under adverse environmental conditions using all available biotechnological tools and high-throughput phenotyping. This review focuses on the current state-of-the-art of biotechnological tools such as high throughput -omics and phenotyping on grapevine, citrus and olive and their contribution to plant breeding programs.

Development of species specific putative miRNA and its target prediction tool in wheat (Triticum aestivum L.)

MicroRNA are 20-24 nt, non-coding, single stranded molecule regulating traits and stress response. Tissue and time specific expression limits its detection, thus is major challenge in their discovery. Wheat has limited 119 miRNAs in MiRBase due to limitation of conservation based methodology where old and new miRNA genes gets excluded. This is due to origin of hexaploid wheat by three successive hybridization, older AA, BB and younger DD subgenome. Species specific miRNA prediction (SMIRP concept) based on 152 thermodynamic features of training dataset using support vector machine learning approach has improved prediction accuracy to 97.7%. This has been implemented in TamiRPred ( http://webtom.cabgrid.res.in/tamirpred ). We also report highest number of putative miRNA genes (4464) of wheat from whole genome sequence populated in database developed in PHP and MySQL. TamiRPred has predicted 2092 (>45.10%) additional miRNA which was not predicted by miRLocator. Predicted miRNAs have been validated by miRBase, small RNA libraries, secondary structure, degradome dataset, star miRNA and binding sites in wheat coding region. This tool can accelerate miRNA polymorphism discovery to be used in wheat trait improvement. Since it predicts chromosome-wise miRNA genes with their respective physical location thus can be transferred using linked SSR markers. This prediction approach can be used as model even in other polyploid crops.

Exploring miRNAs for developing climate-resilient crops: A perspective review

Climate changes and environmental stresses have significant implications on global crop production and necessitate developing crops that can withstand an array of climate changes and environmental perturbations such as irregular water-supplies leading to drought or water-logging, hyper soil-salinity, extreme and variable temperatures, ultraviolet radiations and metal stress. Plants have intricate molecular mechanisms to cope with these dynamic environmental changes, one of the most common and effective being the reprogramming of expression of stress-responsive genes. Plant microRNAs (miRNAs) have emerged as key post-transcriptional and translational regulators of gene-expression for modulation of stress implications. Recent reports are establishing their key roles in epigenetic regulations of stress/adaptive responses as well as in providing plants genome-stability. Several stress responsive miRNAs are being identified from different crop plants and miRNA-driven RNA-interference (RNAi) is turning into a technology of choice for improving crop traits and providing phenotypic plasticity in challenging environments. Here we presents a perspective review on exploration of miRNAs as potent targets for engineering crops that can withstand multi-stress environments via loss-/gain-of-function approaches. This review also shed a light on potential roles plant miRNAs play in genome-stability and their emergence as potent target for genome-editing. Current knowledge on plant miRNAs, their biogenesis, function, their targets, and latest developments in bioinformatics approaches for plant miRNAs are discussed. Though there are recent reviews discussing primarily the individual miRNAs responsive to single stress factors, however, considering practical limitation of this approach, special emphasis is given in this review on miRNAs involved in responses and adaptation of plants to multi-stress environments including at epigenetic and/or epigenomic levels.

miRNomes involved in imparting thermotolerance to crop plants

Thermal stress is one of the challenges to crop plants that negatively impacts crop yield. To overcome this ever-growing problem, utilization of regulatory mechanisms, especially microRNAs (miRNAs), that provide efficient and precise regulation in a targeted manner have been found to play determining roles. Besides their roles in plant growth and development, many recent studies have shown differential regulation of several miRNAs during abiotic stresses including heat stress (HS). Thus, understanding the underlying mechanism of miRNA-mediated gene expression during HS will enable researchers to exploit this regulatory mechanism to address HS responses. This review focuses on the miRNAs and regulatory networks that were involved in physiological, metabolic and morphological adaptations during HS in plant, specifically in crops. Illustrated examples including, the miR156-SPL, miR169-NF-YA5, miR395-APS/AST, miR396-WRKY, etc., have been discussed in specific relation to the crop plants. Further, we have also discussed the available plant miRNA databases and bioinformatics tools useful for miRNA identification and study of their regulatory role in response to HS. Finally, we have briefly discussed the future prospects about the miRNA-related mechanisms of HS for improving thermotolerance in crop plants.

Characterization of miRNAs and their target genes in He-Ne laser pretreated wheat seedlings exposed to drought stress

Drought stress is considered a critical environmental factor that negatively affects wheat growth and development, which causes considerable losses in wheat yields worldwide. More recently, numerous microRNAs (miRNAs) have been found to be involved in wheat responses to drought stresses. However, there is little information regarding the effects of He-Ne laser irradiation on the expression traits of miRNAs and their targets in wheat seedlings exposed to drought stress. In the current study, therefore, a combination of physiological and molecular approaches was used to assess the effect of He-Ne laser irradiation on the expression of miRNAs and their targets in wheat seedlings exposed to drought stress. Our results demonstrated that drought stress significantly reduced plant height, root length, shoot and root fresh weight, relative water content, the expression level and activity of superoxide dismutase (SOD), enhanced malondialdehyde (MDA) concentration in the wheat seedlings. However, He-Ne laser irradiation significantly enhanced the activities of SOD, ascorbate peroxidase (APX), peroxidase (POD) and relative water content, and reduced MDA concentration of seedlings by regulating gene expression for SOD, POD, APX. In addition, in comparison with drought stress alone, miR160, miR164 and miR398 transcripts were down-regulated, and expression levels of its targets auxin response factor (ARF22), NAC domain transcription factor and Cu/Zn superoxide dismutases (CSD) were up-regulated in He-Ne laser irradiated seedlings exposed to drought stress. These results suggested that He-Ne laser irradiation could possible protection of drought stress, at least partially, by increasing the transcript levels and activities of SOD, POD and APX, and decreasing the transcript levels of miR160, miR164 and miR398. To the best of our knowledge, this is the first study to present biochemical and molecular evidence supporting the effect of He-Ne laser irradiation on the alleviation of drought stress in wheat seedlings mediated by miRNA expression.

De novo characterization of microRNAs in oriental fruit moth Grapholita molesta and selection of reference genes for normalization of microRNA expression

MicroRNAs (miRNAs) are a group of endogenous non-coding small RNAs that have critical regulatory functions in almost all known biological processes at the post-transcriptional level in a variety of organisms. The oriental fruit moth Grapholita molesta is one of the most serious pests in orchards worldwide and threatens the production of Rosacea fruits. In this study, a de novo small RNA library constructed from mixed stages of G. molesta was sequenced through Illumina sequencing platform and a total of 536 mature miRNAs consisting of 291 conserved and 245 novel miRNAs were identified. Most of the conserved and novel miRNAs were detected with moderate abundance. The miRNAs in the same cluster normally showed correlated expressional profiles. A comparative analysis of the 79 conserved miRNA families within 31 arthropod species indicated that these miRNA families were more conserved among insects and within orders of closer phylogenetic relationships. The KEGG pathway analysis and network prediction of target genes indicated that the complex composed of miRNAs, clock genes and developmental regulation genes may play vital roles to regulate the developmental circadian rhythm of G. molesta. Furthermore, based on the sRNA library of G. molesta, suitable reference genes were selected and validated for study of miRNA transcriptional profile in G. molesta under two biotic and six abiotic experimental conditions. This study systematically documented the miRNA profile in G. molesta, which could lay a foundation for further understanding of the regulatory roles of miRNAs in the development and metabolism in this pest and might also suggest clues to the development of genetic-based techniques for agricultural pest control.