Functions of microRNAs in plant stress responses

Identification of Micro Ribonucleic Acids and Their Targets in Response to Plasmodiophora brassicae Infection in Brassica napus

Clubroot disease, which is caused by the soil-borne pathogen Plasmodiophora brassicae War (P. brassicae), is one of the oldest and most destructive diseases of Brassica and cruciferous crops in the world. Plant microRNAs [micro ribonucleic acids (miRNAs)] play important regulatory roles in several developmental processes. Although the role of plant miRNAs in plant-microbe interaction has been extensively studied, there are only few reports on the specific functions of miRNAs in response to P. brassicae. This study investigated the roles of miRNAs and their targets during P. brassicae infection in a pair of Brassica napus near-isogenic lines (NILs), namely clubroot-resistant line 409R and clubroot-susceptible line 409S. Small RNA sequencing (sRNA-seq) and degradome-seq were performed on root samples of 409R and 409S with or without P. brassicae inoculation. sRNA-seq identified a total of 48 conserved and 72 novel miRNAs, among which 18 had a significant differential expression in the root of 409R, while only one miRNA was differentially expressed in the root of 409S after P. brassicae inoculation. The degradome-seq analysis identified 938 miRNA target transcripts, which are transcription factors, enzymes, and proteins involved in multiple biological processes and most significantly enriched in the plant hormone signal transduction pathway. Between 409R and 409S, we found eight different degradation pathways in response to P. brassicae infection, such as those related to fatty acids. By combining published transcriptome data, we identified a total of six antagonistic miRNA-target pairs in 409R that are responsive to P. brassicae infection and involved in pathways associated with root development, hypersensitive cell death, and chloroplast metabolic synthesis. Our results reveal that P. brassicae infection leads to great changes in miRNA pool and target transcripts. More interestingly, these changes are different between 409R and 409S. Clarification of the crosstalk between miRNAs and their targets may shed new light on the possible mechanisms underlying the pathogen resistance against P. brassicae.

Noncoding-RNA-Mediated Regulation in Response to Macronutrient Stress in Plants

Macronutrient elements including nitrogen (N), phosphorus (P), potassium (K), calcium (Ca), magnesium (Mg), and sulfur (S) are required in relatively large and steady amounts for plant growth and development. Deficient or excessive supply of macronutrients from external environments may trigger a series of plant responses at phenotypic and molecular levels during the entire life cycle. Among the intertwined molecular networks underlying plant responses to macronutrient stress, noncoding RNAs (ncRNAs), mainly microRNAs (miRNAs) and long ncRNAs (lncRNAs), may serve as pivotal regulators for the coordination between nutrient supply and plant demand, while the responsive ncRNA-target module and the interactive mechanism vary among elements and species. Towards a comprehensive identification and functional characterization of nutrient-responsive ncRNAs and their downstream molecules, high-throughput sequencing has produced massive omics data for comparative expression profiling as a first step. In this review, we highlight the recent findings of ncRNA-mediated regulation in response to macronutrient stress, with special emphasis on the large-scale sequencing efforts for screening out candidate nutrient-responsive ncRNAs in plants, and discuss potential improvements in theoretical study to provide better guidance for crop breeding practices.

A small cog in a large wheel: crucial role of miRNAs in root apical meristem patterning

In both animal and plants, establishment of body axes is fundamental for proper organ development. Plant roots show two main developmental axes: the proximo-distal axis, which spans from the hypocotyl-root junction to the root tip; and the radial axis, which traverses from the vascular tissue to the epidermis. Root axes are determined in the root meristem. The root meristem occupies the tip of the root and contains self-renewing stem cells, which continuously produce new root cells. An intricate network of signalling pathways regulates meristem function and patterning to ensure proper root development and growth. In the last decade, miRNAs, 20-21 nucleotide-long molecules with morphogenetic activity, emerged as central regulators of root cell patterning. Their activity intersects with master regulators of meristematic activity, including phytohormones. In this review, we discuss the latest findings about the activity of miRNAs and their interaction with other molecular networks in the formation of root meristem axes. Furthermore, we describe how these small molecules allow root growth to adapt to changes in the environment, while maintaining the correct patterning.

Cryptic promoter activation occurs by at least two different mechanisms in the Arabidopsis genome

In gene-trap screening of plant genomes, promoterless reporter constructs are often expressed without trapping of annotated gene promoters. The molecular basis of this phenomenon, which has been interpreted as the trapping of cryptic promoters, is poorly understood. Here, we found that cryptic promoter activation occurs by at least two different mechanisms using Arabidopsis gene-trap lines in which a firefly luciferase (LUC) open reading frame (ORF) without an apparent promoter sequence was expressed from intergenic regions: one mechanism is 'cryptic promoter capturing', in which the LUC ORF captured pre-existing promoter-like chromatin marked by H3K4me3 and H2A.Z, and the other is 'promoter de novo origination', in which the promoter chromatin was newly formed near the 5' end of the inserted LUC ORF. The latter finding raises a question as to how the inserted LUC ORF sequence is involved in this phenomenon. To examine this, we performed a model experiment with chimeric LUC genes in transgenic plants. Using Arabidopsis psaH1 promoter-LUC constructs, we found that the functional core promoter region, where transcription start sites (TSSs) occur, cannot simply be determined by the upstream nor core promoter sequences; rather, its positioning proximal to the inserted LUC ORF sequence was more critical. This result suggests that the insertion of the coding sequence alters the local distribution of TSSs in the plant genome. The possible impact of the two types of cryptic promoter activation mechanisms on plant genome evolution and endosymbiotic gene transfer is discussed.

Long noncoding RNA lncRNA354 functions as a competing endogenous RNA of miR160b to regulate ARF genes in response to salt stress in upland cotton

Long non-coding RNAs (lncRNAs) play important roles in response to biotic and abiotic stress through acting as competing endogenous RNAs (ceRNAs) to decoy mature miRNAs. However, whether this mechanism is involved in cotton salt stress response remains unknown. We report the characterization of an endogenous lncRNA, lncRNA354, whose expression was reduced in salt-treated cotton and was localized at the nucleus and cytoplasm. Using endogenous target mimic (eTM) analysis, we predicted that lncRNA354 had a potential binding site for miR160b. Transient expression in tobacco demonstrated that lncRNA354 was a miR160b eTM and attenuated miR160b suppression of its target genes, including auxin response factors (ARFs). Silencing or overexpressing lncRNA354 affected the expression of miR160b and target ARFs. Silencing lncRNA354 and targets GhARF17/18 resulted in taller cotton plants and enhanced the resistant to salt stress. Overexpression of lncRNA354 and targets GhARF17/18 in Arabidopsis led to dwarf plants, decreased root dry weight and reduced salt tolerance. Our results show that the lncRNA354-miR160b effect on GhARF17/18 expression may modulate auxin signalling and thus affect growth. These results also shed new light on a mechanism of lncRNA-associated responses to salt stress.

Comprehensive Mechanism of Gene Silencing and Its Role in Plant Growth and Development

Gene silencing is a negative feedback mechanism that regulates gene expression to define cell fate and also regulates metabolism and gene expression throughout the life of an organism. In plants, gene silencing occurs via transcriptional gene silencing (TGS) and post-transcriptional gene silencing (PTGS). TGS obscures transcription via the methylation of 5' untranslated region (5'UTR), whereas PTGS causes the methylation of a coding region to result in transcript degradation. In this review, we summarized the history and molecular mechanisms of gene silencing and underlined its specific role in plant growth and crop production.

A novel Transposable element-derived microRNA participates in plant immunity to rice blast disease

MicroRNAs (miRNAs) are small non-coding RNAs that direct post-transcriptional gene silencing in plant development and stress responses through cleavage or translational repression of target mRNAs. Here, we report the identification and functional characterization of a new member of the miR812 family in rice (named as miR812w) involved in disease resistance. miR812w is present in cultivated Oryza species, both japonica and indica subspecies, and wild rice species within the Oryza genus, but not in dicotyledonous species. miR812w is a 24nt-long that requires DCL3 for its biogenesis and is loaded into AGO4 proteins. Whereas overexpression of miR812w increased resistance to infection by the rice blast fungus Magnaporthe oryzae, CRISPR/Cas9-mediated MIR812w editing enhances disease susceptibility, supporting that miR812w plays a role in blast resistance. We show that miR812w derives from the Stowaway type of rice MITEs (Miniature Inverted-Repeat Transposable Elements). Moreover, miR812w directs DNA methylation in trans at target genes that have integrated a Stowaway MITE copy into their 3' or 5' untranslated region (ACO3, CIPK10, LRR genes), as well as in cis at the MIR812w locus. The target genes of miR812 were found to be hypo-methylated around the miR812 recognition site, their expression being up-regulated in transgene-free CRISPR/Cas9-edited miR812 plants. These findings further support that, in addition to post-transcriptional regulation of gene expression, miRNAs can exert their regulatory function at the transcriptional level. This relationship between miR812w and Stowaway MITEs integrated into multiple coding genes might eventually create a network for miR812w-mediated regulation of gene expression with implications in rice immunity.

Physiological characteristics and miRNA sequencing of two root zones with contrasting ammonium assimilation patterns in Populus

BACKGROUND

The net ammonium fluxes differ among the different root zones of Populus, but the physiological and microRNA regulatory mechanisms are unclear.

OBJECTIVE

To elucidate the physiological and miRNA regulatory mechanisms, we investigated the two root zones displaying significant differences in net NH₄⁺ effluxes of P. × canescens.

METHODS

Populus plantlets were cultivated with 500 μM NH₄Cl for 10 days. Six plants were randomly selected to determine the net NH₄⁺ fluxes using a noninvasive microtest technique. High-throughput sequencing were used to determine the dynamic expression profile of miRNA among the different root zones of Populus.

RESULTS

Net NH₄⁺ efflux in zone I (from 0 to 40 mm from the root apex) was - 19.64 pmol cm^-2 s^-1 and in zone II (from 40 to 80 mm) it was - 43.96 pmol cm^-2 s^-1. The expression of eleven miRNAs was significantly upregulated, whereas fifteen miRNAs were downregulated. Moreover, eighty-eight target genes of the significantly differentially expressed miRNAs were identified in root zone II compared with zone I. Particularly, ptc-miR171a/b/e and their target, SCL6, were found to be important for the difference in net NH₄⁺ effluxes in the two root zones. Moreover, the expression of the target of ptc-miR169d, NFYA3 was upregulated in root zone II compared with root zone I, contributing to increased NH₄⁺ efflux and decreased NH₄⁺ assimilation in root zone II.

CONCLUSION

These results indicate that miRNAs regulate the expression levels of their target genes and thus play key roles in net NH₄⁺ fluxes and NH₄⁺ assimilation in different poplar root zones.

Genome-wide analysis of the Universal stress protein A gene family in Vitis and expression in response to abiotic stress

Universal Stress Protein A (USPA) plays critical roles in the regulation of growth, development and response to abiotic stress in plants. To date, most research related to the role of USPA in plants has been carried out in herbaceous models such as Arabidopsis, rice and soybean. Here, we used bioinformatics approaches to identify 21 USPA genes in the genome of Vitis vinifera L. Phylogenetic analysis revealed that VvUSPAs could be divided into eight clades. Based on predicted chromosomal locations, we identified 16 pairs of syntenic, orthologous genes between A. thaliana and V. vinifera. Further promoter cis-elements analysis, together with identification of potential microRNA (miRNA) binding sites, suggested that at least some of the VvUSPAs participate in response to phytohormones and abiotic stress. To add support for this, we analyzed the developmental and stress-responsive expression patterns of the homologous USPA genes in the drought-resistant wild Vitis yeshanensis accession 'Yanshan-1' and the drought-sensitive Vitis riparia accession 'He'an'. Most of the USPA genes were upregulated in different degrees in the two genotypes after drought stress and exposure to ethephon (ETH), abscisic acid (ABA) and methyl jasmonate (MeJA). Individual USPA genes showed various tissue-specific expression patterns. Heterologous expression of five selected genes (VvUSPA2, VvUSPA3, VvUSPA11, VvUSPA13 and VvUSPA16) in Escherichia coli (E. coli) enhanced resistance to drought stress. Our study provides a model for mapping gene function in response to abiotic stress and identified three candidate genes, VvUSPA3, VvUSPA11 and VvUSPA16, as regulators of drought response in V. vinifera.

Hydrogen Sulfide Alleviates Alkaline Salt Stress by Regulating the Expression of MicroRNAs in Malus hupehensis Rehd. Roots

Malus hupehensis Rehd. var. pingyiensis Jiang (Pingyi Tiancha, PYTC) is an excellent apple rootstock and ornamental tree, but its tolerance to salt stress is weak. Our previous study showed that hydrogen sulfide (H₂S) could alleviate damage in M. hupehensis roots under alkaline salt stress. However, the molecular mechanism of H₂S mitigation alkaline salt remains to be elucidated. MicroRNAs (miRNAs) play important regulatory roles in plant response to salt stress. Whether miRNAs are involved in the mitigation of alkaline salt stress mediated by H₂S remains unclear. In the present study, through the expression analysis of miRNAs and target gene response to H₂S and alkaline salt stress in M. hupehensis roots, 115 known miRNAs (belonging to 37 miRNA families) and 15 predicted novel miRNAs were identified. In addition, we identified and analyzed 175 miRNA target genes. We certified the expression levels of 15 miRNAs and nine corresponding target genes by real-time quantitative PCR (qRT-PCR). Interestingly, H₂S pretreatment could specifically induce the downregulation of mhp-miR408a expression, and upregulated mhp-miR477a and mhp-miR827. Moreover, root architecture was improved by regulating the expression of mhp-miR159c and mhp-miR169 and their target genes. These results suggest that the miRNA-mediated regulatory network participates in the process of H₂S-mitigated alkaline salt stress in M. hupehensis roots. This study provides a further understanding of miRNA regulation in the H₂S mitigation of alkaline salt stress in M. hupehensis roots.

Heat stress in Marchantia polymorpha: Sensing and mechanisms underlying a dynamic response

Sensing and response to high temperatures are crucial to prevent heat-related damage and to preserve cellular and metabolic functions. The response to heat stress is a complex and coordinated process that involves several subcellular compartments and multi-level regulatory networks that are synchronized to avoid cell damage while maintaining cellular homeostasis. In this review, we provide an insight into the most recent advances in elucidating the molecular mechanisms involved in heat stress sensing and response in Marchantia polymorpha. Based on the signaling pathways and genes that were identified in Marchantia, our analyses indicate that although with specific particularities, the core components of the heat stress response seem conserved in bryophytes and angiosperms. Liverworts not only constitute a powerful tool to study heat stress response and signaling pathways during plant evolution, but also provide key and simple mechanisms to cope with extreme temperatures. Given the increasing prevalence of high temperatures around the world as a result of global warming, this knowledge provides a new set of molecular tools with potential agronomical applications.

Current status of microRNA-mediated regulation of drought stress responses in cereals

Drought is one of the most important abiotic stress factors impeding crop productivity. With the uncovering of their role as potential regulators of gene expression, microRNAs (miRNAs) have been recognized as new targets for developing stress resistance. MicroRNAs are small noncoding RNAs whose abundance is significantly altered under stress conditions. Interestingly, plant miRNAs predominantly targets transcription factors (TFs), and some of which are also the most critical drought-responsive genes that in turn could regulate the expression of numerous loci with drought-adaptive potential. The phytohormone ABA plays important roles in regulating stomatal conductance and in initiating an adaptive response to drought stress. miRNAs are implicated in regulating ABA-(abscisic acid) and non-ABA-mediated drought resistance pathways. For instance, miR159-MYB module and miR169-NFYA module participates in an ABA-dependent pathway, whereas several other ABA-independent miRNA-target modules (miR156-SPL; miR393-TIR1; miR160-ARF10, ARF16, ARF17; miR167-ARF6 and ARF8; miR390/TAS3siRNA-ARF2, ARF3, ARF4) collectively regulate drought responses in plants. Overall, miRNA-mediated drought response manifests diverse molecular, biochemical and physiological processes. Because of their immense role in controlling gene expression, miRNA manipulation has significant potential to augment plant tolerance to drought stress. This review compiles the current understanding of drought-responsive miRNAs in major cereals. Also, potential miRNA manipulation strategies currently in use along with the challenges and future perspectives are discussed.

MicroRNA Omics Analysis of Camellia sinesis Pollen Tubes in Response to Low-Temperature and Nitric Oxide

Nitric oxide (NO) as a momentous signal molecule participates in plant reproductive development and responds to various abiotic stresses. Here, the inhibitory effects of the NO-dominated signal network on the pollen tube growth of Camellia sinensis under low temperature (LT) were studied by microRNA (miRNA) omics analysis. The results showed that 77 and 71 differentially expressed miRNAs (DEMs) were induced by LT and NO treatment, respectively. Gene ontology (GO) analysis showed that DEM target genes related to microtubules and actin were enriched uniquely under LT treatment, while DEM target genes related to redox process were enriched uniquely under NO treatment. In addition, the target genes of miRNA co-regulated by LT and NO are only located on the cell membrane and cell wall, and most of them are enriched in metal ion binding and/or transport and cell wall organization. Furthermore, DEM and its target genes related to metal ion binding/transport, redox process, actin, cell wall organization and carbohydrate metabolism were identified and quantified by functional analysis and qRT-PCR. In conclusion, miRNA omics analysis provides a complex signal network regulated by NO-mediated miRNA, which changes cell structure and component distribution by adjusting Ca²⁺ gradient, thus affecting the polar growth of the C. sinensis pollen tube tip under LT.

Transcriptome and Degradome Profiling Reveals a Role of miR530 in the Circadian Regulation of Gene Expression in Kalanchoë marnieriana

Crassulacean acid metabolism (CAM) is an important photosynthetic pathway for plant adaptation to dry environments. CAM plants feature a coordinated interaction between mesophyll and epidermis functions that involves refined regulations of gene expression. Plant microRNAs (miRNAs) are crucial post-transcription regulators of gene expression, however, their roles underlying the CAM pathway remain poorly investigated. Here, we present a study characterizing the expression of miRNAs in an obligate CAM species Kalanchoë marnieriana. Through sequencing of transcriptome and degradome in mesophyll and epidermal tissues under the drought treatments, we identified differentially expressed miRNAs that were potentially involved in the regulation of CAM. In total, we obtained 84 miRNA genes, and eight of them were determined to be Kalanchoë-specific miRNAs. It is widely accepted that CAM pathway is regulated by circadian clock. We showed that miR530 was substantially downregulated in epidermal peels under drought conditions; miR530 targeted two tandem zinc knuckle/PLU3 domain encoding genes (TZPs) that were potentially involved in light signaling and circadian clock pathways. Our work suggests that the miR530-TZPs module might play a role of regulating CAM-related gene expression in Kalanchoë.

High-throughput miRNA deep sequencing in response to drought stress in sugarcane

Drought is a major factor which reduces cane growth and productivity. In the present study, we sequenced drought susceptible (V1) and drought tolerant (V2) sugarcane varieties using high-throughput miRNA deep sequencing method to study the regulation of gene expression by miRNAs during drought stress in sugarcane. A total of 1224 conserved miRNAs which belong to 89 miRNA families were identified and 38% of the differentially regulated miRNAs were common for both varieties. Additionally 435 novel miRNAs were also identified from four small RNA libraries. We identified 145 miRNAs that were differentially expressed in susceptible variety (V1-31) and 143 miRNAs differentially expressed in the tolerant variety (V2-31). Target prediction revealed that the genes mainly encoded transcription factors, proteins, phosphatase and kinases involved in signal transduction pathways, integral component of membrane and inorganic ion transport metabolism, enzymes involved in carbohydrate transport and metabolism and drought-stress-related proteins involved in defense mechanisms. Pathway analysis of targets revealed that "General function prediction only" was the most significant pathway observed in both tolerant and susceptible genotypes followed by "signal transduction mechanisms". Functional annotation of the transcripts revealed genes like calcium-dependent protein kinase, respiratory burst oxidase, caffeic acid 3-O-methyltransferase, peroxidase, calmodulin, glutathione S-transferase and transcription factors like MYB, WRKY that are involved in drought tolerant pathways. qRT-PCR was used to verify the expression levels of miRNAs and their potential targets obtained from RNA sequencing results.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-021-02857-x.

Molecular regulation of zinc deficiency responses in plants

Zinc (Zn) is an essential micronutrient for plants and animals. Because of its low availability in arable soils worldwide, Zn deficiency is becoming a serious agricultural problem resulting in decreases of crop yield and nutritional quality. Plants have evolved multiple responses to adapt to low levels of soil Zn supply, involving biochemical and physiological changes to improve Zn acquisition and utilization, and defend against Zn deficiency stress. In this review, we summarize the physiological and biochemical adaptations of plants to Zn deficiency, the roles of transporters and metal-binding compounds in Zn homeostasis regulation, and the recent progresses in understanding the sophisticated regulatory mechanisms of Zn deficiency responses that have been made by molecular and genetic analyses, as well as diverse 'omics' studies. Zn deficiency responses are tightly controlled by multiple layers of regulation, such as transcriptional regulation that is mediated by transcription factors like F-group bZIP proteins, epigenetic regulation at the level of chromatin, and post-transcriptional regulation mediated by small RNAs and alternative splicing. The insights into the regulatory network underlying Zn deficiency responses and the perspective for further understandings of molecular regulation of Zn deficiency responses have been discussed. The understandings of the regulatory mechanisms will be important for improving Zn deficiency tolerance, Zn use efficiency, and Zn biofortification in plants.

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

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

The miR172/IDS1 signaling module confers salt tolerance through maintaining ROS homeostasis in cereal crops

Salt stress triggers the overdose accumulation of reactive oxygen species (ROS) in crop plants, leading to severe oxidative damage to living tissues. MicroRNAs (miRNAs) act as master regulators orchestrating the stress responsive regulatory networks as well as salt tolerance. However, the fundamental roles of miRNAs in modulating salt tolerance in cereal crops, especially in salt-triggered ROS scavenging remain largely unknown. Through small RNA sequencing, a salt-responsive miRNA, miR172 was identified in rice. Further, by generating the miR172-overexpression or MIR172 gene loss-of-function mutant lines, the biological significance of miR172 and its downstream signaling pathways related to salt tolerance were defined. We demonstrated that miR172 is a positive regulator of salt tolerance in both rice and wheat. More interestingly, miR172a and miR172b, but not miR172c or miR172d are involved in salt stress response, emphasizing the functional differentiation within miR172 family members. Further evidence uncovers a novel miR172/IDS1 regulatory module that functions as a crucial molecular rheostat in maintaining ROS homeostasis during salt stress, mainly through balancing the expression of a group of ROS-scavenging genes. Our findings establish a direct molecular link between miRNAs and detoxification response in cereal crops for improving salt tolerance.

MIR172d Is Required for Floral Organ Identity and Number in Tomato

MicroRNA172 (miR172) functions as a central regulator of flowering time and flower development by post-transcriptional repression of APETALA2-LIKE transcription factors. In the model crop Solanum lycopersicum (tomato), the miR172 family is still poorly annotated and information about the functions of specific members is lacking. Here, de-novo prediction of tomato miR172 coding loci identified seven genes (SlMIR172a-g), that code for four unique species of miR172 (sly-miR172). During reproductive development, sly-miR172s are differentially expressed, with sly-miR172c and sly-miR172d being the most abundant members in developing flowers, and are predicted to guide the cleavage of eight APETALA2-LIKE transcription factors. By CRISPR-Cas9 co-targeting of SlMIR172c and SlMIR172d we have generated a battery of loss-of-function and hypomorphic mutants (slmir172c-d^CR). The slmir172c-d^CR plants developed normal shoot but their flowers displayed graded floral organ abnormalities. Whereas slmir172c^CR loss-of-function caused only a slight greening of petals and stamens, hypomorphic and loss-of-function slmir172d^CR alleles were associated with the conversion of petals and stamens to sepaloids, which were produced in excess. Interestingly, the degrees of floral organ identity alteration and proliferation were directly correlated with the reduction in sly-miR172d activity. These results suggest that sly-miR172d regulates in a dose-dependent manner floral organ identity and number, likely by negatively regulating its APETALA2-like targets.

Splice Variants of Superoxide Dismutases in Rice and Their Expression Profiles under Abiotic Stresses

The superoxide dismutases (SODs) play vital roles in controlling cellular reactive oxygen species (ROS) that are generated both under optimal as well as stress conditions in plants. The rice genome harbors seven SOD genes (CSD1, CSD2, CSD3, CSD4, FSD1, FSD2, and MSD) that encode seven constitutive transcripts. Of these, five (CSD2, CSD3, CSD4, FSD1, and MSD) utilizes an alternative splicing (AS) strategy and generate seven additional splice variants (SVs) or mRNA variants, i.e., three for CSD3, and one each for CSD2, CSD4, FSD1, and MSD. The exon-intron organization of these SVs revealed variations in the number and length of exons and/or untranslated regions (UTRs). We determined the expression patterns of SVs along with their constitutive forms of SODs in rice seedlings exposed to salt, osmotic, cold, heavy metal (Cu⁺²) stresses, as well as copper-deprivation. The results revealed that all seven SVs were transcriptionally active in both roots and shoots. When compared to their corresponding constitutive transcripts, the profiles of five SVs were almost similar, while two specific SVs (CSD3-SV4 and MSD-SV2) differed significantly, and the differences were also apparent between shoots and roots suggesting that the specific SVs are likely to play important roles in a tissue-specific and stress-specific manner. Overall, the present study has provided a comprehensive analysis of the SVs of SODs and their responses to stress conditions in shoots and roots of rice seedlings.

Genome-wide identification and expression analysis of detoxification efflux carriers (DTX) genes family under abiotic stresses in flax

The detoxification efflux carriers (DTX)/multidrug and toxic compound extrusion (MATE) transporters encompass an ancient gene family of secondary transporters involved in the process of plant detoxification. A genome-wide analysis of these transporters was carried out in order to better understand the transport of secondary metabolites in flaxseed genome (Linum usitassimum). A total of 73 genes coding for DTX/MATE transporters were identified. Gene structure, protein domain and motif organization were found to be notably conserved over the distinct phylogenetic groups, showing the evolutionary significant role of each class. Gene ontology (GO) annotation revealed a link to transporter activities, response to stimulus and localizations. The presence of various hormone and stress-responsive cis-regulatory elements in promoter regions could be directly correlated with the alteration of their transcripts. Tertiary structure showed conservation for pore size and constrains in the pore, which indicate their involvement in the exclusion of toxic substances from the cell. MicroRNA target analysis revealed that LuDTXs genes were targeted by different classes of miRNA families. Twelve LuDTX genes were chosen for further quantitative real-time polymerase chain reaction analysis in response to cold, salinity and cadmium stress at 0, 6, 12 and 24 hours after treatment. Altogether, the identified members of the DTX gene family, their expression profile, phylogenetic and miRNAs analysis might provide opportunities for future functional validation of this important gene family in flax.

MicroRNAs: Potential Targets for Developing Stress-Tolerant Crops

Crop yield is challenged every year worldwide by changing climatic conditions. The forecasted climatic scenario urgently demands stress-tolerant crop varieties to feed the ever-increasing global population. Molecular breeding and genetic engineering approaches have been frequently exploited for developing crops with desired agronomic traits. Recently, microRNAs (miRNAs) have emerged as powerful molecules, which potentially serve as expression markers during stress conditions. The miRNAs are small non-coding endogenous RNAs, usually 20-24 nucleotides long, which mediate post-transcriptional gene silencing and fine-tune the regulation of many abiotic- and biotic-stress responsive genes in plants. The miRNAs usually function by specifically pairing with the target mRNAs, inducing their cleavage or repressing their translation. This review focuses on the exploration of the functional role of miRNAs in regulating plant responses to abiotic and biotic stresses. Moreover, a methodology is also discussed to mine stress-responsive miRNAs from the enormous amount of transcriptome data available in the public domain generated using next-generation sequencing (NGS). Considering the functional role of miRNAs in mediating stress responses, these molecules may be explored as novel targets for engineering stress-tolerant crop varieties.

The role of miRNA in plant-virus interaction: a review

Plant viruses affect crop production both quantitatively and qualitatively. The viral genome consists of either DNA or RNA. However, most plant viruses are positive single-strand RNA viruses. MicroRNAs are involved in gene regulation and affect development as well as host-virus interaction. They are non-coding short with 20-24 nucleotides long capable of regulating gene expression. The miRNA gene is transcribed by RNA polymerase II to form pri-miRNA which will later cleaved by Dicer-like 1 to produce pre-miRNA with the help of HYPONASTIC LEAVES1 and SERRATE which finally methylated and exported via nucleopore with the help of HASTY. The outcome of plant virus interaction depends on the effectiveness of host defense and the ability of a virus counter-defense mechanism. In plants, miRNAs are involved in the repression of gene expression through transcript cleavage. On the other hand, viruses use viral suppressors of RNA silencing (VSRs) which affect RISC assembly and subsequent mRNA degradation. Passenger strands, miRNA*, have a significant biological function in plant defense response as well as plant development.

ASCORBATE PEROXIDASE6 delays the onset of age-dependent leaf senescence

Age-dependent changes in reactive oxygen species (ROS) levels are critical in leaf senescence. While H2O2-reducing enzymes such as catalases and cytosolic ASCORBATE PEROXIDASE1 (APX1) tightly control the oxidative load during senescence, their regulation and function are not specific to senescence. Previously, we identified the role of ASCORBATE PEROXIDASE6 (APX6) during seed maturation in Arabidopsis (Arabidopsis thaliana). Here, we show that APX6 is a bona fide senescence-associated gene. APX6 expression is specifically induced in aging leaves and in response to senescence-promoting stimuli such as abscisic acid (ABA), extended darkness, and osmotic stress. apx6 mutants showed early developmental senescence and increased sensitivity to dark stress. Reduced APX activity, increased H2O2 level, and altered redox state of the ascorbate pool in mature pre-senescing green leaves of the apx6 mutants correlated with the early onset of senescence. Using transient expression assays in Nicotiana benthamiana leaves, we unraveled the age-dependent post-transcriptional regulation of APX6. We then identified the coding sequence of APX6 as a potential target of miR398, which is a key regulator of copper redistribution. Furthermore, we showed that mutants of SQUAMOSA PROMOTER BINDING PROTEIN-LIKE7 (SPL7), the master regulator of copper homeostasis and miR398 expression, have a higher APX6 level compared with the wild type, which further increased under copper deficiency. Our study suggests that APX6 is a modulator of ROS/redox homeostasis and signaling in aging leaves that plays an important role in developmental- and stress-induced senescence programs.

Potential Mechanisms Linking Food-Derived MicroRNAs, Gut Microbiota and Intestinal Barrier Functions in the Context of Nutrition and Human Health

MicroRNAs (miRNAs) are non-coding single-stranded RNA molecules from 18 to 24 nucleotides that are produced by prokaryote and eukaryote organisms, which play a crucial role in regulating gene expression through binding to their mRNA targets. MiRNAs have acquired special attention for their potential in cross kingdom communication, notably food-derived microRNAs (xenomiRs), which could have an impact on microorganism and mammal physiology. In this review, we mainly aim to deal with new perspectives on: (1) The mechanism by which food-derived xenomiRs (mainly dietary plant xenomiRs) could be incorporated into humans through diet, in a free form, associated with proteins or encapsulated in exosome-like nanoparticles. (2) The impact of dietary plant-derived miRNAs in modulating gut microbiota composition, which in turn, could regulate intestinal barrier permeability and therefore, affect dietary metabolite, postbiotics or food-derived miRNAs uptake efficiency. Individual gut microbiota signature/composition could be also involved in xenomiR uptake efficiency through several mechanisms such us increasing the bioavailability of exosome-like nanoparticles miRNAs. (3) Gut microbiota dysbiosis has been proposed to contribute to disease development by affecting gut epithelial barrier permeability. For his reason, the availability and uptake of dietary plant xenomiRs might depend, among other factors, on this microbiota-related permeability of the intestine. We hypothesize and critically review that xenomiRs-microbiota interaction, which has been scarcely explored yet, could contribute to explain, at least in part, the current disparity of evidences found dealing with dietary miRNA uptake and function in humans. Furthermore, dietary plant xenomiRs could be involved in the establishment of the multiple gut microenvironments, in which microorganism would adapt in order to optimize the resources and thrive in them. Additionally, a particular xenomiR could preferentially accumulate in a specific region of the gastrointestinal tract and participate in the selection and functions of specific gut microbial communities.

Elucidating the morpho-physiological adaptations and molecular responses under long-term waterlogging stress in maize through gene expression analysis

Waterlogging stress in maize is one of the emerging abiotic stresses in the current climate change scenario. To gain insights in transcriptional reprogramming during late hours of waterlogging stress under field conditions, we aimed to elucidate the transcriptional and anatomical changes in two contrasting maize inbreds viz. I110 (susceptible) and I172 (tolerant). Waterlogging stress reduced dry matter translocations from leaves and stems to ears, resulting in a lack of sink capacity and inadequate grain filling in I110, thus decreased the grain yield drastically. The development of aerenchyma cells within 48 h in I172 enabled hypoxia tolerance. The upregulation of alanine aminotransferase, ubiquitin activating enzyme E1, putative mitogen activated protein kinase and pyruvate kinase in I172 suggested that genes involved in protein degradation, signal transduction and carbon metabolism provided adaptive mechanisms during waterlogging. Overexpression of alcohol dehydrogenase, sucrose synthase, aspartate aminotransferase, NADP dependent malic enzyme and many miRNA targets in I110 indicated that more oxygen and energy consumption might have shortened plant survival during long-term waterlogging exposure. To the best of our knowledge, this is the first report of transcript profiling at late stage (24-96 h) of waterlogging stress under field conditions and provides new visions to understand the molecular basis of waterlogging tolerance in maize.

GhARF16-1 modulates leaf development by transcriptionally regulating the GhKNOX2-1 gene in cotton

The leaf is a crucial organ evolved with remarkable morphological diversity to maximize plant photosynthesis. The leaf shape is a key trait that affects photosynthesis, flowering rates, disease resistance and yield. Although many genes regulating leaf development have been identified in the past years, the precise regulatory architecture underlying the generation of diverse leaf shapes remains to be elucidated. We used cotton as a reference model to probe the genetic framework underlying divergent leaf forms. Comparative transcriptome analysis revealed that the GhARF16-1 and GhKNOX2-1 genes might be potential regulators of leaf shape. We functionally characterized the auxin-responsive factor ARF16-1 acting upstream of GhKNOX2-1 to determine leaf morphology in cotton. The transcription of GhARF16-1 was significantly higher in lobed-leaved cotton than in smooth-leaved cotton. Furthermore, the overexpression of GhARF16-1 led to the up-regulation of GhKNOX2-1 and resulted in more and deeper serrations in cotton leaves, similar to the leaf shape of cotton plants overexpressing GhKNOX2-1. We found that GhARF16-1 specifically bound to the promoter of GhKNOX2-1 to induce its expression. The heterologous expression of GhARF16-1 and GhKNOX2-1 in Arabidopsis led to lobed and curly leaves, and a genetic analysis revealed that GhKNOX2-1 is epistatic to GhARF16-1 in Arabidopsis, suggesting that the GhARF16-1 and GhKNOX2-1 interaction paradigm also functions to regulate leaf shape in Arabidopsis. To our knowledge, our results uncover a novel mechanism by which auxin, through the key component ARF16-1 and its downstream-activated gene KNOX2-1, determines leaf morphology in eudicots.

Pu-miR172d regulates stomatal density and water-use efficiency via targeting PuGTL1 in poplar

miRNAs play essential regulatory roles in many aspects of plant development and in responses to abiotic and biotic stresses. Here, we characterize Pu-miR172d, which acts as a negative regulator of stomatal density by directly repressing the expression of PuGTL1 in Populus ussuriensis. Quantitative real-time PCR and GUS reporter analyses showed that Pu-miR172d was strongly expressed in the guard cells of young leaves. Overexpression of Pu-miR172d significantly decreased stomatal density, resulting in increases in water use efficiency (WUE) and drought tolerance by reducing net photosynthetic rate, stomatal conductance, and transpiration. Molecular analysis showed that PuGTL1 was a major target of Pu-miR172d cleavage. Moreover, PuGTL1-SRDX plants, in which PuGTL1 is suppressed, phenocopied Pu-miR172d-overexpression lines with reduced stomatal density and enhanced WUE. The expression of PuSDD1, a negative regulator of stomatal development, was significantly increased in young leaves of both Pu-miR172d-overexpression and PuGTL1-SRDX plants. RNA-seq analysis of mature leaves indicated that overexpression of Pu-miR172d decreased the expression of many genes related to photosynthesis. Our findings show that the Pu-miR172d/PuGTL1/PuSDD1 module plays an important role in stomatal differentiation, and hence it is a potential target for engineering improved drought tolerance in poplar.

Increased water use efficiency in miR396-downregulated tomato plants

MicroRNAs regulate plant development and responses to biotic and abiotic stresses but their impact on water use efficiency (WUE) is poorly known. Increasing WUE is a major task in crop improvement programs aimed to meet the challenges posed by the reduction in water availability associated with the ongoing climatic change. We have examined the physiological and molecular response to water stress of tomato (Solanum lycopersicum L.) plants downregulated for miR396 by target mimicry. In water stress conditions, miR396-downregulated plants displayed reduced transpiration and a less then proportional decrease in the photosynthetic rate that resulted in higher WUE. The increase in WUE was associated with faster foliar accumulation of abscisic acid (ABA), with the induction of several drought-protective genes and with the activation of the jasmonic acid (JA) and γ-aminobutyric acid (GABA) pathways. We propose a model in which the downregulation of miR396 leads to the activation of a complex molecular response to water stress. This response acts synergistically with a set of leaf morphological modifications to increase stomatal closure and preserve the efficiency of the photosynthetic activity, ultimately resulting in higher WUE.

Profiling of MicroRNAs and Their Targets in Roots and Shoots Reveals a Potential MiRNA-Mediated Interaction Network in Response to Phosphate Deficiency in the Forestry Tree Betula luminifera

Inorganic phosphate (Pi) is often lacking in natural and agro-climatic environments, which impedes the growth of economically important woody species. Plants have developed strategies to cope with low Pi (LP) availability. MicroRNAs (miRNAs) play important roles in responses to abiotic stresses, including nutrition stress, by regulating target gene expression. However, the miRNA-mediated regulation of these adaptive responses and their underlying coordinating signals are still poorly understood in forestry trees such as Betula luminifera. Transcriptomic libraries, small RNA (sRNA) libraries, and a mixed degradome cDNA library of B. luminifera roots and shoots treated under LP and normal conditions (CK) were constructed and sequenced using next-generation deep sequencing. A comprehensive B. luminifera transcriptome derived from its roots and shoots was constructed, and a total of 76,899 unigenes were generated. Analysis of the transcriptome identified 8,095 and 5,584 differentially expressed genes in roots and shoots, respectively, under LP conditions. sRNA sequencing analyses indicated that 66 and 60 miRNAs were differentially expressed in roots and shoots, respectively, under LP conditions. A total of 109 and 112 miRNA-target pairs were further validated in the roots and shoots, respectively, using degradome sequencing. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of differential miRNA targets indicated that the "ascorbate and aldarate metabolism" pathway responded to LP, suggesting miRNA-target pairs might participating in the removing of reactive oxidative species under LP stress. Moreover, a putative network of miRNA-target interactions involved in responses to LP stress in B. luminifera is proposed. Taken together, these findings provide useful information to decipher miRNA functions and establish a framework for exploring P signaling networks regulated by miRNAs in B. luminifera and other woody plants. It may provide new insights into the genetic engineering of high use efficiency of Pi in forestry trees.

Biotechnological strategies for enhancing heavy metal tolerance in neglected and underutilized legume crops: A comprehensive review

Contamination of agricultural land and water by heavy metals due to rapid industrialization and urbanization including various natural processes have become one of the major constraints to crop growth and productivity. Several studies have reported that to counteract heavy metal stress, plants should be able to maneuver various physiological, biochemical and molecular processes to improve their growth and development under heavy metal stress. With the advent of modern biotechnological tools and techniques it is now possible to tailor legume and other plants overexpressing stress-induced genes, transcription factors, proteins, and metabolites that are directly involved in heavy metal stress tolerance. This review provides an in-depth overview of various biotechnological approaches and/or strategies that can be used for enhancing detoxification of the heavy metals by stimulating phytoremediation processes. Synthetic biology tools involved in the engineering of legume and other crop plants against heavy metal stress tolerance are also discussed herewith some pioneering examples where synthetic biology tools that have been used to modify plants for specific traits. Also, CRISPR based genetic engineering of plants, including their role in modulating the expression of several genes/ transcription factors in the improvement of abiotic stress tolerance and phytoremediation ability using knockdown and knockout strategies has also been critically discussed.

Methods for Predicting CircRNA-miRNA-mRNA Regulatory Networks: GreenCircRNA and PlantCircNet Databases as Study Cases

Circular RNAs are molecules formed by 3'-5' ligation in a splicing reaction, the so-called backsplicing. Well described in other groups, especially in humans, circRNA studies that include prediction and validation in plants are recent. It has already been shown that circRNAs can interact with microRNAs, acting as sponges, and adding a new layer of complexity in regulating eukaryotic transcription. Here, we cover two up-to-date databases that allow the users to perform analyses of the circRNA-miRNA-mRNA interactions in plants. We choose two databases to demonstrate their functions and compare their approaches to obtain a more robust and reliable interaction network.

Rice-specific Argonaute 17 controls reproductive growth and yield-associated phenotypes

This manuscript describes the functions of an Argonaute protein named AGO17 in rice. AGO17 is required for the development of rice reproductive tissues. Argonaute (AGO) proteins are a well-conserved multigene family of regulators mediating gene silencing across eukaryotes. Monocot plants have additional members of AGO, the functions of which are poorly understood. Among the non-dicot AGO1 clade members in monocots, AGO17 expresses highly in reproductive tissues. Here we show that overexpression of Oryza sativa indica AGO17 in rice resulted in robust growth and increased yield, whereas its silencing resulted in reduced panicle length, less fertility, and poor growth. Small (s)RNA transcriptome analysis revealed misregulation of several miRNAs and other categories of sRNAs in silenced and overexpression lines, in agreement with its likely competition with other AGO1 clade members. Targets of differentially expressed miRNAs included previously unreported target RNAs coding for proteins involved in development, phase transition, and transport. Our results indicate a distinctive role for OsAGO17 in rice reproductive development that could be harnessed to improve yield.

Microarray analysis of Arabidopsis thaliana exposed to single and mixed infections with Cucumber mosaic virus and turnip viruses

Cucumber mosaic virus (CMV), Turnip mosaic virus (TuMV) and Turnip crinkle virus (TCV) are important plant infecting viruses. In the present study, whole transcriptome alteration of Arabidopsis thaliana in response to CMV, TuMV and TCV, individual as well as mixed infections of CMV and TuMV/CMV and TCV were investigated using microarray data. In response to CMV, TuMV and TCV infections, a total of 2517, 3985 and 277 specific differentially expressed genes (DEGs) were up-regulated, while 2615, 3620 and 243 specific DEGs were down-regulated, respectively. The number of 1222 and 30 common DEGs were up-regulated during CMV and TuMV as well as CMV and TCV infections, while 914 and 24 common DEGs were respectively down-regulated. Genes encoding immune response mediators, signal transducer activity, signaling and stress response functions were among the most significantly upregulated genes during CMV and TuMV or CMV and TCV mixed infections. The NAC, C3H, C2H2, WRKY and bZIP were the most commonly presented transcription factor (TF) families in CMV and TuMV infection, while AP2-EREBP and C3H were the TF families involved in CMV and TCV infections. Moreover, analysis of miRNAs during CMV and TuMV and CMV and TCV infections have demonstrated the role of miRNAs in the down regulation of host genes in response to viral infections. These results identified the commonly expressed virus-responsive genes and pathways during plant-virus interaction which might develop novel antiviral strategies for improving plant resistance to mixed viral infections.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-00925-3.

Response of Root Growth and Development to Nitrogen and Potassium Deficiency as well as microRNA-Mediated Mechanism in Peanut (Arachis hypogaea L.)

The mechanism of miRNA-mediated root growth and development in response to nutrient deficiency in peanut (Arachis hypogaea L.) is still unclear. In the present study, we found that both nitrogen (N) and potassium (K) deficiency resulted in a significant reduction in plant growth, as indicated by the significantly decreased dry weight of both shoot and root tissues under N or K deficiency. Both N and K deficiency significantly reduced the root length, root surface area, root volume, root vitality, and weakened root respiration, as indicated by the reduced O₂ consuming rate. N deficiency significantly decreased primary root length and lateral root number, which might be associated with the upregulation of miR160, miR167, miR393, and miR396, and the downregulation of AFB3 and GRF. The primary and lateral root responses to K deficiency were opposite to that of the N deficiency condition. The upregulated miR156, miR390, NAC4, ARF2, and AFB3, and the downregulated miR160, miR164, miR393, and SPL10 may have contributed to the growth of primary roots and lateral roots under K deficiency. Overall, roots responded differently to the N or K deficiency stresses in peanuts, potentially due to the miRNA-mediated pathway and mechanism.

Genome-wide identification, characterization and expression profile of glutaredoxin gene family in relation to fruit ripening and response to abiotic and biotic stresses in banana (Musa acuminata)

Glutaredoxins (GRXs) are disulfide oxidoreductases that are involved in various biological processes. However, little information on the role of GRXs in the regulation of fruit ripening and the response to stress is available. In this study, we isolated 64 GRX genes from banana genome. Their encoded GRX proteins could be classified into four classes: CC, CGFS, CPYC and GRL types. The distribution and synteny of these GRXs on chromosomes, the gene structures, the promoter sequences, and the possible protein subcellular localizations were characterized. Molecular interaction network analysis suggested that MaGRX might interact with glutathione reductase (GR), sulfiredoxin, peroxiredoxin (Prx), and NADPH-dependent thioredoxin reductase C (NTRC), contributing to the antioxidative defense of banana fruit. MicroRNA prediction showed that MaGRX genes might be targeted by different miRNAs. Transcriptome analysis characterized the expression profiles of different MaGRX genes during banana fruit ripening, and in response to different storage stresses. The results suggested that CC-type, CPYC-type and GRL-type MaGRXs might be more active than CGFS-type MaGRXs during banana fruit ripening and the response to stress. Moreover, MaGRX6/7/9/11/17/23/28 and MaGRL3/16/19 might play important roles in regulating fruit ripening or in response to low and high temperature, or Fusarium proliferatum infection.

RNAi Mediated Hypoxia Stress Tolerance in Plants

Small RNAs regulate various biological process involved in genome stability, development, and adaptive responses to biotic or abiotic stresses. Small RNAs include microRNAs (miRNAs) and small interfering RNAs (siRNAs). MicroRNAs (miRNAs) are regulators of gene expression that affect the transcriptional and post-transcriptional regulation in plants and animals through RNA interference (RNAi). miRNAs are endogenous small RNAs that originate from the processing of non-coding primary miRNA transcripts folding into hairpin-like structures. The mature miRNAs are incorporated into the RNA-induced silencing complex (RISC) and drive the Argonaute (AGO) proteins towards their mRNA targets. siRNAs are generated from a double-stranded RNA (dsRNA) of cellular or exogenous origin. siRNAs are also involved in the adaptive response to biotic or abiotic stresses. The response of plants to hypoxia includes a genome-wide transcription reprogramming. However, little is known about the involvement of RNA signaling in gene regulation under low oxygen availability. Interestingly, miRNAs have been shown to play a role in the responses to hypoxia in animals, and recent evidence suggests that hypoxia modulates the expression of various miRNAs in plant systems. In this review, we describe recent discoveries on the impact of RNAi on plant responses to hypoxic stress in plants.

Multiomics analyses unveil the involvement of microRNAs in pear fruit senescence under high- or low-temperature conditions

Senescence leads to declines in fruit quality and shortening of shelf life. It is known that low temperatures (LTs) efficiently delay fruit senescence and that high temperatures (HTs) accelerate senescence. However, the molecular mechanism by which temperature affects senescence is unclear. Herein, through multiomics analyses of fruits subjected to postharvest HT, LT, and room temperature treatments, a total of 56 metabolic compounds and 700 mRNAs were identified to be associated with fruit senescence under HT or LT conditions. These compounds could be divided into antisenescent (I→III) and prosenescent (IV→VI) types. HT affected the expression of 202 mRNAs to enhance the biosynthesis of prosenescent compounds of types V and VI and to inhibit the accumulation of antisenescent compounds of types II and III. LT affected the expression of 530 mRNAs to promote the accumulation of antisenescent compounds of types I and II and to impede the biosynthesis of prosenescent compounds of types IV and V. Moreover, 16 microRNAs were isolated in response to HT or LT conditions and interacted with the mRNAs associated with fruit senescence under HT or LT conditions. Transient transformation of pear fruit showed that one of these microRNAs, Novel_188, can mediate fruit senescence by interacting with its target Pbr027651.1. Thus, both HT and LT conditions can affect fruit senescence by affecting microRNA-mRNA interactions, but the molecular networks are different in pear fruit.

Characterizing microRNAs and their targets in different organs of Camellia sinensis var. assamica

To comprehensively annotate miRNAs and their targets in tea plant, Camellia sinensis, we sequenced small and messenger RNAs of 9 samples of Camellia sinensis var. assamica (YK-10), a diploid elite cultivar widely grown in southwest China. In order to identify targets of miRNAs, we sequenced two degradome sequencing profiles from leaves and roots of YK-10, respectively. By analyzing the small RNA-Seq profiles, we newly identified 137 conserved miRNAs and 23 species specific miRNAs in the genome of YK-10, which significantly improved the annotation of miRNAs in tea plant. Approximately 2000 differently expressed genes were identified when comparing RNA-Seq profiles of any two of the three organs selected in the study. Totally, more than 5000 targets of conserved miRNAs were identified in the two degradome profiles. Furthermore, our results suggest that a few miRNAs play roles in the biosynthesis pathways of theanine, caffeine and flavonoid. These results enhance our understanding of small RNA guided gene regulations in different organs of tea plant.

Aux/IAA14 Regulates microRNA-Mediated Cold Stress Response in Arabidopsis Roots

The phytohormone auxin and microRNA-mediated regulation of gene expressions are key regulators of plant growth and development at both optimal and under low-temperature stress conditions. However, the mechanistic link between microRNA and auxin in regulating plant cold stress response remains elusive. To better understand the role of microRNA (miR) in the crosstalk between auxin and cold stress responses, we took advantage of the mutants of Arabidopsis thaliana with altered response to auxin transport and signal. Screening of the mutants for root growth recovery after cold stress at 4 °C revealed that the auxin signaling mutant, solitary root 1 (slr1; mutation in Aux/IAA14), shows a hypersensitive response to cold stress. Genome-wide expression analysis of miRs in the wild-type and slr1 mutant roots using next-generation sequencing revealed 180 known and 71 novel cold-responsive microRNAs. Cold stress also increased the abundance of 26–31 nt small RNA population in slr1 compared with wild type. Comparative analysis of microRNA expression shows significant differential expression of 13 known and 7 novel miRs in slr1 at 4 °C compared with wild type. Target gene expression analysis of the members from one potential candidate miR, miR169, revealed the possible involvement of miR169/NF-YA module in the Aux/IAA14-mediated cold stress response. Taken together, these results indicate that SLR/IAA14, a transcriptional repressor of auxin signaling, plays a crucial role in integrating miRs in auxin and cold responses.

Plant Non-Coding RNAs: Origin, Biogenesis, Mode of Action and Their Roles in Abiotic Stress

As sessile species, plants have to deal with the rapidly changing environment. In response to these environmental conditions, plants employ a plethora of response mechanisms that provide broad phenotypic plasticity to allow the fine-tuning of the external cues related reactions. Molecular biology has been transformed by the major breakthroughs in high-throughput transcriptome sequencing and expression analysis using next-generation sequencing (NGS) technologies. These innovations have provided substantial progress in the identification of genomic regions as well as underlying basis influencing transcriptional and post-transcriptional regulation of abiotic stress response. Non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), short interfering RNAs (siRNAs), and long non-coding RNAs (lncRNAs), have emerged as essential regulators of plants abiotic stress response. However, shared traits in the biogenesis of ncRNAs and the coordinated cross-talk among ncRNAs mechanisms contribute to the complexity of these molecules and might play an essential part in regulating stress responses. Herein, we highlight the current knowledge of plant microRNAs, siRNAs, and lncRNAs, focusing on their origin, biogenesis, modes of action, and fundamental roles in plant response to abiotic stresses.

Genome-wide identification of low phosphorus responsive microRNAs in two soybean genotypes by high-throughput sequencing

MicroRNAs (miRNAs) have been reported to be correlated with various stress responses in soybean, but only a few miRNAs have been demonstrated to respond to low phosphorus (LP) stress. To unravel the response mechanisms of miRNAs to low-P stress, the roots of two representative soybean genotypes with different P efficiency, Nannong94-156 (a LP-tolerant genotype) and Bogao (a LP-sensitive genotype), were used for the construction of RNA sequencing (RNA-seq) libraries under low/normal-P treatment by high-throughput sequencing. In total, 603 existing miRNAs and 1699 novel miRNAs belonging to 248 and 1582 families in all samples were identified, respectively. Among these miRNAs, 777 miRNAs were differentially expressed (DE) across different P levels and genotypes. Furthermore, putative targets of DE miRNAs were predicted, and these miRNAs mainly targeted ERF (ethylene responsive factor), auxin response factors (ARF), zinc finger protein, MYB, and NAC domain transcription factors. Gene ontology (GO) analysis showed that targets of DE miRNAs were significantly enriched in binding, metabolic processes, biological regulation, response to stress, and phosphorus metabolic processes. In addition, the expression profiles of chosen P-responsive miRNAs and target genes were validated by quantitative real-time PCR (qRT-PCR). Our study focused on genome-wide miRNA identification in two representative soybean genotypes under low-P stress. Overall, the DE miRNAs across different P levels and genotypes and their putative target genes will provide useful information for further study of miRNAs mediating low-P response and facilitate improvements in soybean breeding.

Genome wide in-silico miRNA and target network prediction from stress responsive Horsegram (Macrotyloma uniflorum) accessions

Horsegram (Macrotyloma uniflorum (Lam.) Verdc.) is a drought hardy food and fodder legume of Indo-African continents with diverse germplasm sources demonstrating alternating mechanisms depicting contrasting adaptations to different climatic zones. Tissue specific expression of genes contributes substantially to location specific adaptations. Regulatory networks of such adaptive genes are elucidated for downstream translational research. MicroRNAs are small endogenous regulatory RNAs which alters the gene expression profiles at a particular time and type of tissue. Identification of such small regulatory RNAs in low moisture stress hardy crops can help in cross species transfer and validation confirming stress tolerance ability. This study outlined prediction of conserved miRNAs from transcriptome shotgun assembled sequences and EST sequences of horsegram. We could validate eight out of 15 of the identified miRNAs to demonstrate their role in deficit moisture stress tolerance mechanism of horsegram variety Paiyur1 with their target networks. The putative mumiRs were related to other food legumes indicating the presence of gene regulatory networks. Differential miRNA expression among drought specific tissues indicted the probable energy conservation mechanism. Targets were identified for functional characterization and regulatory network was constructed to find out the probable pathways of post-transcriptional regulation. The functional network revealed mechanism of biotic and abiotic stress tolerance, energy conservation and photoperiod responsiveness.

Abiotic stress signalling in extremophile land plants

Plant life relies on complex arrays of environmental stress sensing and signalling mechanisms. Extremophile plants develop and grow in harsh environments with extremes of cold, heat, drought, desiccation, or salinity, which have resulted in original adaptations. In accordance with their polyphyletic origins, extremophile plants likely possess core mechanisms of plant abiotic stress signalling. However, novel properties or regulations may have emerged in the context of extremophile adaptations. Comparative omics of extremophile genetic models, such as Arabidopsis lyrata, Craterostigma plantagineum, Eutrema salsugineum, and Physcomitrella patens, reveal diverse strategies of sensing and signalling that lead to a general improvement in abiotic stress responses. Current research points to putative differences of sensing and emphasizes significant modifications of regulatory mechanisms, at the level of secondary messengers (Ca2+, phospholipids, reactive oxygen species), signal transduction (intracellular sensors, protein kinases, transcription factors, ubiquitin-mediated proteolysis) or signalling crosstalk. Involvement of hormone signalling, especially ABA signalling, cell homeostasis surveillance, and epigenetic mechanisms, also shows that large-scale gene regulation, whole-plant integration, and probably stress memory are important features of adaptation to extreme conditions. This evolutionary and functional plasticity of signalling systems in extremophile plants may have important implications for plant biotechnology, crop improvement, and ecological risk assessment under conditions of climate change.

Physiological characteristics and RNA sequencing in two root zones with contrasting nitrate assimilation of Populus × canescens

Different root zones have distinct capacities for nitrate (NO3-) uptake in Populus species, but the underlying physiological and microRNA (miRNA) regulatory mechanisms remain largely unknown. To address this question, two root zones of Populus × canescens (Ait.) Smith. with contrasting capacities for NO3- uptake were investigated. The region of 0-40 mm (root zone I) to the root apex displayed net influxes, whereas the region of 40-80 mm (root zone II) exhibited net effluxes. Concentrations of NO3- and ammonium (NH4+) as well as nitrate reductase activity were lower in zone II than in zone I. Forty one upregulated and twenty three downregulated miRNAs, and 576 targets of these miRNAs were identified in zone II in comparison with zone I. Particularly, growth-regulating factor 4 (GRF4), a target of upregulated ptc-miR396g-5p and ptc-miR396f_L + 1R-1, was downregulated in zone II in comparison with zone I, probably contributing to lower NO3- uptake rates and assimilation in zone II. Furthermore, several miRNAs and their targets, members of C2H2 zinc finger family and APETALA2/ethylene-responsive element binding protein family, were found in root zones, which probably play important roles in regulating NO3- uptake. These results indicate that differentially expressed miRNA-target pairs play key roles in regulation of distinct NO3- uptake rates and assimilation in different root zones of poplars.