MicroRNA and Transcription Factor: Key Players in Plant Regulatory Network

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.

Proteomic and Transcriptomic Analyses Provide Novel Insights into the Crucial Roles of Host-Induced Carbohydrate Metabolism Enzymes in Xanthomonas oryzae pv. oryzae Virulence and Rice-Xoo Interaction

BACKGROUND

Xanthomonas oryzae pv. oryzae (Xoo) causes bacterial leaf blight, a devastating rice disease. The Xoo-rice interaction, wherein wide ranging host- and pathogen-derived proteins and genes wage molecular arms race, is a research hotspot. Hence, the identification of novel rice-induced Xoo virulence factors and characterization of their roles affecting rice global gene expression profiles will provide an integrated and better understanding of Xoo-rice interactions from the molecular perspective.

RESULTS

Using comparative proteomics and an in vitro interaction system, we revealed that 5 protein spots from Xoo exhibited significantly different expression patterns (|fold change| > 1.5) at 3, 6, 12 h after susceptible rice leaf extract (RLX) treatment. MALDI-TOF MS analysis and pathogenicity tests showed that 4 host-induced proteins, including phosphohexose mutase, inositol monophosphatase, arginase and septum site-determining protein, affected Xoo virulence. Among them, mutants of two host-induced carbohydrate metabolism enzyme-encoding genes, ΔxanA and Δimp, elicited enhanced defense responses and nearly abolished Xoo virulence in rice. To decipher rice differentially expressed genes (DEGs) associated with xanA and imp, transcriptomic responses of ΔxanA-treated and Δimp-treated susceptible rice were compared to those in rice treated with PXO99^A at 1 and 3 dpi. A total of 1521 and 227 DEGs were identified for PXO99^A vs Δimp at 1 and 3 dpi, while for PXO99^A vs ΔxanA, there were 131 and 106 DEGs, respectively. GO, KEGG and MapMan analyses revealed that the DEGs for PXO99^A vs Δimp were mainly involved in photosynthesis, signal transduction, transcription, oxidation-reduction, hydrogen peroxide catabolism, ion transport, phenylpropanoid biosynthesis and metabolism of carbohydrates, lipids, amino acids, secondary metabolites, hormones, and nucleotides, while the DEGs from PXO99^A vs ΔxanA were predominantly associated with photosynthesis, signal transduction, oxidation-reduction, phenylpropanoid biosynthesis, cytochrome P450 and metabolism of carbohydrates, lipids, amino acids, secondary metabolites and hormones. Although most pathways were associated with both the Δimp and ΔxanA treatments, the underlying genes were not the same.

CONCLUSION

Our study identified two novel host-induced virulence factors XanA and Imp in Xoo, and revealed their roles in global gene expression in susceptible rice. These results provide valuable insights into the molecular mechanisms of pathogen infection strategies and plant immunity.

Long Non-Coding RNA and Its Regulatory Network Response to Cold Stress in Eucalyptus urophylla S.T.Blake

Long non-coding RNA (lncRNA) plays an important regulatory role in plant growth and development, but its systematic identification and analysis in Eucalyptus has not yet been reported. Cold stress has a huge impact on the survival and yield of Eucalyptus seedlings, but the regulatory mechanism of lncRNA in Eucalyptus in response to cold stress is still unclear. In this study, the transcriptomes of young leaves of Eucalyptus urophylla S.T.Blake under low-temperature treatment and restoration were analyzed by RNA-seq. A total of 11,394 lncRNAs and 46,276 mRNAs were identified, of which 300 were differentially expressed lncRNAs (DE_lncRNAs) and 5606 were differentially expressed target genes of lncRNAs under cold stress, with the total number of target genes of DE_lncRNAs being 1681. A total of 677 differentially expressed transcription factors (TFs) were also identified, mainly including ERF, MYB and the NAC transcription factor family. Gene ontology (GO) analysis of the differentially expressed genes (DEGs) and target genes of DE_lncRNAs was mostly related to the response to cold stress and external stimuli. Furthermore, lncRNA–miRNA–mRNA regulatory networks were constructed, and 22 DE_lncRNAs were predicted to be targets or targeting mimics of 20 miRNAs. A qRT-PCR was used to verify the relative expression of genes in the regulatory EuGBF3-EUC_00002677-MSTRG.7690 network, and it matched the transcriptome data, indicating that it may play an important role in the response to cold stress in E. urophylla. This study provides a new insight into lncRNA and its regulatory network under abiotic stress, especially cold stress in E. urophylla.

Genome-Wide Identification and Characterization of Small Peptides in Maize

Small peptides (sPeptides), <100 amino acids (aa) long, are encoded by small open reading frames (sORFs) often found in the 5' and 3' untranslated regions (or other parts) of mRNAs, in long non-coding RNAs, or transcripts from introns and intergenic regions; various sPeptides play important roles in multiple biological processes. In this study, we conducted a comprehensive study of maize (Zea mays) sPeptides using mRNA sequencing, ribosome profiling (Ribo-seq), and mass spectrometry (MS) on six tissues (each with at least two replicates). To identify maize sORFs and sPeptides from these data, we set up a robust bioinformatics pipeline and performed a genome-wide scan. This scan uncovered 9,388 sORFs encoding peptides of 2-100 aa. These sORFs showed distinct genomic features, such as different Kozak region sequences, higher specificity of translation, and high translational efficiency, compared with the canonical protein-coding genes. Furthermore, the MS data verified 2,695 sPeptides. These sPeptides perfectly discriminated all the tissues and were highly associated with their parental genes. Interestingly, the parental genes of sPeptides were significantly enriched in multiple functional gene ontology terms related to abiotic stress and development, suggesting the potential roles of sPeptides in the regulation of their parental genes. Overall, this study lays out the guidelines for genome-wide scans of sORFs and sPeptides in plants by integrating Ribo-seq and MS data and provides a more comprehensive resource of functional sPeptides in maize and gives a new perspective on the complex biological systems of plants.

Investigation of long non-coding RNAs as regulatory players of grapevine response to powdery and downy mildew infection

BACKGROUND

Long non-coding RNAs (lncRNAs) are regulatory transcripts of length > 200 nt. Owing to the rapidly progressing RNA-sequencing technologies, lncRNAs are emerging as considerable nodes in the plant antifungal defense networks. Therefore, we investigated their role in Vitis vinifera (grapevine) in response to obligate biotrophic fungal phytopathogens, Erysiphe necator (powdery mildew, PM) and Plasmopara viticola (downy mildew, DM), which impose huge agro-economic burden on grape-growers worldwide.

RESULTS

Using computational approach based on RNA-seq data, 71 PM- and 83 DM-responsive V. vinifera lncRNAs were identified and comprehensively examined for their putative functional roles in plant defense response. V. vinifera protein coding sequences (CDS) were also profiled based on expression levels, and 1037 PM-responsive and 670 DM-responsive CDS were identified. Next, co-expression analysis-based functional annotation revealed their association with gene ontology (GO) terms for 'response to stress', 'response to biotic stimulus', 'immune system process', etc. Further investigation based on analysis of domains, enzyme classification, pathways enrichment, transcription factors (TFs), interactions with microRNAs (miRNAs), and real-time quantitative PCR of lncRNAs and co-expressing CDS pairs suggested their involvement in modulation of basal and specific defense responses such as: Ca²⁺-dependent signaling, cell wall reinforcement, reactive oxygen species metabolism, pathogenesis related proteins accumulation, phytohormonal signal transduction, and secondary metabolism.

CONCLUSIONS

Overall, the identified lncRNAs provide insights into the underlying intricacy of grapevine transcriptional reprogramming/post-transcriptional regulation to delay or seize the living cell-dependent pathogen growth. Therefore, in addition to defense-responsive genes such as TFs, the identified lncRNAs can be further examined and leveraged to candidates for biotechnological improvement/breeding to enhance fungal stress resistance in this susceptible fruit crop of economic and nutritional importance.

Genome-Wide Identification and Analysis of the APETALA2 (AP2) Transcription Factor in Dendrobium officinale

The APETALA2 (AP2) transcription factors (TFs) play crucial roles in regulating development in plants. However, a comprehensive analysis of the AP2 family members in a valuable Chinese herbal orchid, Dendrobium officinale, or in other orchids, is limited. In this study, the 14 DoAP2 TFs that were identified from the D. officinale genome and named DoAP2-1 to DoAP2-14 were divided into three clades: euAP2, euANT, and basalANT. The promoters of all DoAP2 genes contained cis-regulatory elements related to plant development and also responsive to plant hormones and stress. qRT-PCR analysis showed the abundant expression of DoAP2-2, DoAP2-5, DoAP2-7, DoAP2-8 and DoAP2-12 genes in protocorm-like bodies (PLBs), while DoAP2-3, DoAP2-4, DoAP2-6, DoAP2-9, DoAP2-10 and DoAP2-11 expression was strong in plantlets. In addition, the expression of some DoAP2 genes was down-regulated during flower development. These results suggest that DoAP2 genes may play roles in plant regeneration and flower development in D. officinale. Four DoAP2 genes (DoAP2-1 from euAP2, DoAP2-2 from euANT, and DoAP2-6 and DoAP2-11 from basal ANT) were selected for further analyses. The transcriptional activation of DoAP2-1, DoAP2-2, DoAP2-6 and DoAP2-11 proteins, which were localized in the nucleus of Arabidopsis thaliana mesophyll protoplasts, was further analyzed by a dual-luciferase reporter gene system in Nicotiana benthamiana leaves. Our data showed that pBD-DoAP2-1, pBD-DoAP2-2, pBD-DoAP2-6 and pBD-DoAP2-11 significantly repressed the expression of the LUC reporter compared with the negative control (pBD), suggesting that these DoAP2 proteins may act as transcriptional repressors in the nucleus of plant cells. Our findings on AP2 genes in D. officinale shed light on the function of AP2 genes in this orchid and other plant species.

Identification of miRNAs and Their Targets Involved in Flower and Fruit Development across Domesticated and Wild Capsicum Species

MicroRNAs (miRNAs) are regulators of the post-transcription stage of gene activity documented to play central roles in flower and fruit development in model plant species. However, little is known about their roles and differences in domesticated and wild Capsicum species. In this study, we used high-throughput sequencing to analyze the miRNA content at three developmental stages (flower, small fruit, and middle fruit) from two cultivated (C. baccatum and C. annuum) and two wild (C. chacoense and C. eximium) pepper species. This analysis revealed 22 known and 27 novel miRNAs differentially expressed across species and tissues. A number of stage- and species-specific miRNAs were identified, and Gene Ontology terms were assigned to 138 genes targeted by the miRNAs. Most Gene Ontology terms were for the categories "genetic information processing", "signaling and cellular processes", "amino acid metabolism", and "carbohydrate metabolism". Enriched KEGG analysis revealed the pathways amino acids, sugar and nucleotide metabolism, starch and sucrose metabolism, and fructose-mannose metabolism among the principal ones regulated by miRNAs during pepper fruit ripening. We predicted miRNA-target gene interactions regulating flowering time and fruit development, including miR156/157 with SPL genes, miR159 with GaMYB proteins, miR160 with ARF genes, miR172 with AP2-like transcription factors, and miR408 with CLAVATA1 gene across the different Capsicum species. In addition, novel miRNAs play an important role in regulating interactions potentially controlling plant pathogen defense and fruit quality via fructokinase, alpha-L-arabinofuranosidase, and aromatic and neutral amino acid transporter. Overall, the small RNA-sequencing results from this study represent valuable information that provides a solid foundation for uncovering the miRNA-mediated mechanisms of flower and fruit development between domesticated and wild Capsicum species.

The role of RNA-binding protein, microRNA and alternative splicing in seed germination: a field need to be discovered

Seed germination is the process through which a quiescent organ reactivates its metabolism culminating with the resumption cell divisions. It is usually the growth of a plant contained within a seed and results in the formation of a seedling. Post-transcriptional regulation plays an important role in gene expression. In cells, post-transcriptional regulation is mediated by many factors, such as RNA-binding proteins, microRNAs, and the spliceosome. This review provides an overview of the relationship between seed germination and post-transcriptional regulation. It addresses the relationship between seed germination and RNA-binding proteins, microRNAs and alternative splicing. This presentation of the current state of the knowledge will promote new investigations into the relevance of the interactions between seed germination and post-transcriptional regulation in plants.

Differential miRNA profiles in South African cassava mosaic virus-infected cassava landraces reveal clues to susceptibility and tolerance to cassava mosaic disease

Specific miRNA families are involved in susceptibility or antiviral immunity in plants. Manihot esculenta Crantz (cassava) is a perennial plant that is an important food security crop in sub-Saharan Africa. Cassava is susceptible to several begomoviruses that cause cassava mosaic disease (CMD). In this study, we investigated the leaf miRNAome response in a tolerant (TME3) and susceptible (T200) cassava landrace challenged with South African cassava mosaic virus. RNAseq was performed on leaf samples at 12, 32 and 67 days post infection (dpi), representing early, symptomatic and late persistent stages of CMD infection. Significantly, distinct profiles of conserved miRNA family expression between the T200 and TME3 landraces at the three infection stages were observed. Notably at 12 days post SACMV infection, TME3 exhibited significant downregulation (log2fold<2.0) of 42 %, compared to 9% in T200, of the conserved miRNA families. This demonstrates an overall early response to SACMV in TME3 prior to symptom appearance not observed in T200, and expression of a large cohort of miRNA-regulated genes. Notably, at early infection, downregulation of mes-miR162 and 168 that target antiviral posttransriptional gene silencing (PTGS) regulators DCL1 and AGO1, respectively, was observed in TME3, and AGO1 and DCL1 expression was higher compared to T200 post infection. Early rapid responses prior to symptom development, including RNA silencing, may be key to establishing the tolerance/recovery phenotype exhibited by TME3 landrace later on at 67 dpi. At recovery, TME3 was hallmarked by a highly significant down-regulation of mes-miR167. MiR167 targets an auxin responsive factor which plays a role in auxin signaling and adaptive responses to stress, suggesting the importance of the auxin signaling in recovery of SACMV-induced symptoms. The gene targets of these miRNAs and their associated networks may provide clues to the molecular basis of CMD tolerance in perennial hosts such as cassava.

Abiotic Stress-Responsive miRNA and Transcription Factor-Mediated Gene Regulatory Network in Oryza sativa: Construction and Structural Measure Study

Climate changes and environmental stresses have a consequential association with crop plant growth and yield, meaning it is necessary to cultivate crops that have tolerance toward the changing climate and environmental disturbances such as water stress, temperature fluctuation, and salt toxicity. Recent studies have shown that trans-acting regulatory elements, including microRNAs (miRNAs) and transcription factors (TFs), are emerging as promising tools for engineering naive improved crop varieties with tolerance for multiple environmental stresses and enhanced quality as well as yield. However, the interwoven complex regulatory function of TFs and miRNAs at transcriptional and post-transcriptional levels is unexplored in Oryza sativa. To this end, we have constructed a multiple abiotic stress responsive TF-miRNA-gene regulatory network for O. sativa using a transcriptome and degradome sequencing data meta-analysis approach. The theoretical network approach has shown the networks to be dense, scale-free, and small-world, which makes the network stable. They are also invariant to scale change where an efficient, quick transmission of biological signals occurs within the network on extrinsic hindrance. The analysis also deciphered the existence of communities (cluster of TF, miRNA, and genes) working together to help plants in acclimatizing to multiple stresses. It highlighted that genes, TFs, and miRNAs shared by multiple stress conditions that work as hubs or bottlenecks for signal propagation, for example, during the interaction between stress-responsive genes (TFs/miRNAs/other genes) and genes involved in floral development pathways under multiple environmental stresses. This study further highlights how the fine-tuning feedback mechanism works for balancing stress tolerance and how timely flowering enable crops to survive in adverse conditions. This study developed the abiotic stress-responsive regulatory network, APRegNet database (http://lms.snu.edu.in/APRegNet), which may help researchers studying the roles of miRNAs and TFs. Furthermore, it advances current understanding of multiple abiotic stress tolerance mechanisms.

Insight of transcriptional regulators reveals the tolerance mechanism of carpet-grass (Axonopus compressus) against drought

BACKGROUND

Carpet grass [Axonopus compressus (L.)] is an important warm-season perennial grass around the world and is known for its adaptability to varied environmental conditions. However, Carpet grass lacks enough data in public data banks, which confined our comprehension of the mechanisms of environmental adaptations, gene discovery, and development of molecular markers. In current study, the DEGs (differentially expressed genes) in Axonopus compressus under drought stress (DS) were identified and compared with CK (control) by RNA-Seq.

RESULTS

A total of 263,835 unigenes were identified in Axonopus compressus, and 201,303 (also added to the numbers of the remaining 2 databases) a sequence of unigenes significantly matched in at least one of the seven databases. A total of 153,697 (58.25%) unigenes classified to 144 KEGG pathways, and 7444 unigenes were expressed differentially between DS and CK, of which 4249 were up-regulated and 3195 were down-regulated unigenes. Of the 50 significantly enriched GO terms, 18, 6, and 14 items were related to BP, CC, and MF respectively. Analysis of KEGG enrichment revealed 2569 DEGs involved in 143 different pathways, under drought stress. 2747 DEGs were up-regulated and 2502 DEGs were down-regulated. Moreover, we identified 352 transcription factors (TFs) in Axonopus compressus, of which 270 were differentially expressed between CK and DS. The qRT-PCR validation experiment also supports the transcriptional response of Axonopus compressus against drought. Accuracy of transcriptome unigenes of Axonopus compressus was assessed with BLAST, which showed 3300 sequences of Axonopus compressus in the NCBI.

CONCLUSION

The 7444 unigenes were found to be between DS and CK treatments, which indicate the existence of a strong mechanism of drought tolerance in Axonopus compressus. The current findings provide the first framework for further investigations for the particular roles of these unigenes in Axonopus compressus in response to drought.

Systematic analyses of the MIR172 family members of Arabidopsis define their distinct roles in regulation of APETALA2 during floral transition

MicroRNAs (miRNAs) play important roles in regulating flowering and reproduction of angiosperms. Mature miRNAs are encoded by multiple MIRNA genes that can differ in their spatiotemporal activities and their contributions to gene regulatory networks, but the functions of individual MIRNA genes are poorly defined. We functionally analyzed the activity of all 5 Arabidopsis thaliana MIR172 genes, which encode miR172 and promote the floral transition by inhibiting the accumulation of APETALA2 (AP2) and APETALA2-LIKE (AP2-LIKE) transcription factors (TFs). Through genome editing and detailed confocal microscopy, we show that the activity of miR172 at the shoot apex is encoded by 3 MIR172 genes, is critical for floral transition of the shoot meristem under noninductive photoperiods, and reduces accumulation of AP2 and TARGET OF EAT2 (TOE2), an AP2-LIKE TF, at the shoot meristem. Utilizing the genetic resources generated here, we show that the promotion of flowering by miR172 is enhanced by the MADS-domain TF FRUITFULL, which may facilitate long-term silencing of AP2-LIKE transcription, and that their activities are partially coordinated by the TF SQUAMOSA PROMOTER-BINDING-LIKE PROTEIN 15. Thus, we present a genetic framework for the depletion of AP2 and AP2-LIKE TFs at the shoot apex during floral transition and demonstrate that this plays a central role in floral induction.

Comprehensive analysis of NAC transcription factor family uncovers drought and salinity stress response in pearl millet (Pennisetum glaucum)

BACKGROUND

Pearl millet (Pennisetum glaucum) is a cereal crop that possesses the ability to withstand drought, salinity and high temperature stresses. The NAC [NAM (No Apical Meristem), ATAF1 (Arabidopsis thaliana Activation Factor 1), and CUC2 (Cup-shaped Cotyledon)] transcription factor family is one of the largest transcription factor families in plants. NAC family members are known to regulate plant growth and abiotic stress response. Currently, no reports are available on the functions of the NAC family in pearl millet.

RESULTS

Our genome-wide analysis found 151 NAC transcription factor genes (PgNACs) in the pearl millet genome. Thirty-eight and 76 PgNACs were found to be segmental and dispersed duplicated respectively. Phylogenetic analysis divided these NAC transcription factors into 11 groups (A-K). Three PgNACs (- 073, - 29, and - 151) were found to be membrane-associated transcription factors. Seventeen other conserved motifs were found in PgNACs. Based on the similarity of PgNACs to NAC proteins in other species, the functions of PgNACs were predicted. In total, 88 microRNA target sites were predicted in 59 PgNACs. A previously performed transcriptome analysis suggests that the expression of 30 and 42 PgNACs are affected by salinity stress and drought stress, respectively. The expression of 36 randomly selected PgNACs were examined by quantitative reverse transcription-PCR. Many of these genes showed diverse salt- and drought-responsive expression patterns in roots and leaves. These results confirm that PgNACs are potentially involved in regulating abiotic stress tolerance in pearl millet.

CONCLUSION

The pearl millet genome contains 151 NAC transcription factor genes that can be classified into 11 groups. Many of these genes are either upregulated or downregulated by either salinity or drought stress and may therefore contribute to establishing stress tolerance in pearl millet.

Comparative profiling of the resistance of different genotypes of mannose-binding lectin to Mycoplasma pneumoniae infection in Chinese Merino sheep based on high-throughput sequencing technology

Mannose-binding lectin (MBL) glycoproteins in blood can selectively recognise lectins on the surface of bacteria, and play an important role in natural immunity. Micro RNAs (miRNAs) are key molecules that regulate gene expression at the post-transcriptional level in vivo, and their pathways are specific and effective. Previous studies indicate that small RNAs such as miRNAs perform regulatory roles in immunology. Herein, we investigated differential expression of miRNAs during MBL protein immunotherapy in sheep following treatment with different MBL genotypes (resistant and susceptible), and identified miRNAs linked to different target genes and pathways. RNA was extracted from liver tissue of resistant and susceptible sheep, miRNAs were identified by high-throughput sequencing, and differentially expressed miRNAs were analysed by SOAP to predict target genes and biological pathways. Results: Some miRNAs (oar-mir-143, oar-mir-10b, oar-mir-382, oar-mir-432 and oar-mir-379) were up-regulated, while others were down-regulated. GPATCH3 and DNAJC5 were predicted target genes of oar-mir-379, DMRT1 and GATA4 were linked to oar-mir-382, and oar-mir-432 was associated with STAT2, DMRT1 and ATG16L1. Identification of miRNAs differentially expressed in resistant and susceptible sheep may expand our understanding of miRNAs in immune regulation, and the role of MBL in innate immunity.

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

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

The BTB-TAZ protein MdBT2 negatively regulates the drought stress response by interacting with the transcription factor MdNAC143 in apple

Drought stress is a severe source of abiotic stress that can affect apple yield and quality, yet the underlying molecular mechanism of the drought stress response and the role of MdBT2 in the process remain unclear. Here, we find that MdBT2 negatively regulates the drought stress response. Both in vivo and in vitro assays indicated that MdBT2 interacted physically with and ubiquitinated MdNAC143, a member of the NAC TF family that is a positive regulator under drought stress. In addition, MdBT2 promotes the degradation of MdNAC143 via the 26S proteasome system. A series of transgenic assays in apple calli and Arabidopsis verify that MdBT2 confers susceptibility to drought stress at least in part by the regulation of MdNAC143. Overall, our findings provide new insight into the mechanism of MdBT2, which functions antagonistically to MdNAC143 in regulating drought stress by regulating the potential downstream target protein MdNAC143 for proteasomal degradation in apple.

Identification of genuine and novel miRNAs in Amaranthus hypochondriacus from high-throughput sequencing data

Amaranth has been proposed as an exceptional alternative for food security and climate change mitigation. Information about the distribution, abundance, or specificity of miRNAs in amaranth species is scare. Here, small RNAs from seedlings under control, drought, heat, and cold stress conditions of the Amaranthus hypocondriacus variety "Gabriela" were sequenced and miRNA loci identified in the amaranth genome using the ShortStack software. Fifty-three genuine miRNA clustersthirty-nine belonging to conserved families, and fourteen novel, were identified. Identification of their target genes suggests that conserved amaranth miRNAs are involved in growth and developmental processes, as well as stress responses. MiR0005, an amaranth-specific miRNA, exhibited an unusual high level of expression, akin to that of conserved miRNAs. Overall, our results broaden our knowledge regarding the distribution, abundance and expression of miRNAs in amaranth, providing the basis for future research on miRNAs and their functions in this important species.

Plant Volatile Organic Compounds Evolution: Transcriptional Regulation, Epigenetics and Polyploidy

Volatile organic compounds (VOCs) are emitted by plants as a consequence of their interaction with biotic and abiotic factors, and have a very important role in plant evolution. Floral VOCs are often involved in defense and pollinator attraction. These interactions often change rapidly over time, so a quick response to those changes is required. Epigenetic factors, such as DNA methylation and histone modification, which regulate both genes and transcription factors, might trigger adaptive responses to these evolutionary pressures as well as regulating the rhythmic emission of VOCs through circadian clock regulation. In addition, transgenerational epigenetic effects and whole genome polyploidy could modify the generation of VOCs’ profiles of offspring, contributing to long-term evolutionary shifts. In this article, we review the available knowledge about the mechanisms that may act as epigenetic regulators of the main VOC biosynthetic pathways, and their importance in plant evolution.

Underpinning miRNA-miRNA co-functional interaction patterns in the metabolism of Oryza sativa by genome-scale network analysis

MicroRNA (miRNA) is a class of non-coding small RNAs, which post-transcriptionally regulate a large number of genes and are now known to be important regulators in a wide variety of biological processes including metabolism. Thus, for better understanding these complex biological networks, and to derive their significance and inter-dependency, a systems biology approach enables us to explore and draw vital insights into these molecular network architectures. In this study, we aimed to understand the significance of synergistic miRNA-miRNA interactions in rice by constructing and analysing metabolic networks. The construction of the network involves target gene prediction of experimentally verified miRNAs of rice and then appending associated metabolic pathways to the network. A genome-scale miRNA-miRNA co-functional network (MFSN) is constructed based on co-regulatory interactions among the miRNAs and common target genes by applying transformational procedures. The analysis of the extracted MFSN modules identifies co-regulated target genes that are associated with corresponding interconnected metabolic pathways such as VALDEG-PWY (L-valine degradation I pathway was found to be targeted by multiple miRNA families, such as osa-miR812, osa-miR818, osa-miR821, and osa-miR5799 families while another pathway that was found to be associated with multiple miRNA families was PWY-6952 (glycerophosphodiester degradation pathway), PWY-6952 was found to be targeted by osa-miR812, osa-miR11344 and osa-miR5801 families. Such extensive study will help in systematically elucidating the regulatory networks in metabolism of rice, which in turn can be utilised to devise strategies for crop improvement and novel cultivar development.

OsmiR535, a Potential Genetic Editing Target for Drought and Salinity Stress Tolerance in Oryza sativa

OsmiR535 belongs to the miR156/miR529/miR535 superfamily, a highly conserved miRNA family in plants. OsmiR535 is involved in regulating the cold-stress response, modulating plant development, and determining panicle architecture and grain length. However, the role that OsmiR535 plays in plant responses to drought and salinity are elusive. In the current study, molecular and genetic engineering techniques were used to elucidate the possible role of OsmiR535 in response to NaCl, PEG(Poly ethylene glycol), ABA(Abscisic acid), and dehydration stresses. Our results showed that OsmiR535 is induced under stressed conditions as compared to control. With transgenic and CRISPR/Cas9 knockout system techniques, our results verified that either inhibition or knockout of OsmiR535 in rice could enhance the tolerance of plants to NaCl, ABA, dehydration and PEG stresses. In addition, the overexpression of OsmiR535 significantly reduced the survival rate of rice seedlings during PEG and dehydration post-stress recovery. Our results demonstrated that OsmiR535 negatively regulates the stress response in rice. Moreover, our practical application of CRISPR/Cas9 mediated genome editing created a homozygous 5 bp deletion in the coding sequence of OsmiR535, demonstrating that OsmiR535 could be a useful genetic editing target for drought and salinity tolerance and a new marker for molecular breeding of Oryza sativa.

Computational Identification and Comparative Analysis of Conserved miRNAs and Their Putative Target Genes in the Juglans regia and J. microcarpa Genomes

MicroRNAs (miRNAs) are important factors for the post-transcriptional regulation of protein-coding genes in plants and animals. They are discovered either by sequencing small RNAs or computationally. We employed a sequence-homology-based computational approach to identify conserved miRNAs and their target genes in Persian (English) walnut, Juglans regia, and its North American wild relative, J. microcarpa. A total of 119 miRNA precursors (pre-miRNAs) were detected in the J. regia genome and 121 in the J. microcarpa genome and miRNA target genes were predicted and their functional annotations were performed in both genomes. In the J. regia genome, 325 different genes were targets; 87.08% were regulated by transcript cleavage and 12.92% by translation repression. In the J. microcarpa genome, 316 different genes were targets; 88.92% were regulated by transcript cleavage and 11.08% were regulated by translation repression. Totals of 1.3% and 2.0% of all resistance gene analogues (RGA) and 2.7% and 2.6% of all transcription factors (TFs) were regulated by miRNAs in the J. regia and J. microcarpa genomes, respectively. Juglans genomes evolved by a whole genome duplication (WGD) and consist of eight pairs of fractionated homoeologous chromosomes. Within each pair, the chromosome that has more genes with greater average transcription also harbors more pre-miRNAs and more target genes than its homoeologue. While only minor differences were detected in pre-miRNAs between the J. regia and J. microcarpa genomes, about one-third of the pre-miRNA loci were not conserved between homoeologous chromosome within each genome. Pre-miRNA and their corresponding target genes showed a tendency to be collocated within a subgenome.

smi-miR396b targeted SmGRFs, SmHDT1, and SmMYB37/4 synergistically regulates cell growth and active ingredient accumulation in Salvia miltiorrhiza hairy roots

MIR396b had been cloned and overexpressed in Salvia miltiorrhiza hairy roots. MiR396b targets SmGRFs, SmHDT1, and SmMYB37/4 to regulate cell growth and secondary metabolism in S. miltiorrhiza hairy roots. Danshen (Salvia miltiorrhiza Bunge) is a valuable medicinal herb with two kinds of clinically used natural products, salvianolic acids and tanshinones. miR396 is a conserved microRNA and plays extensive roles in plants. However, it is still unclear how miR396 works in S. miltiorrhiza. In this study, an smi-MIR396b has been cloned from S. miltiorrhiza. Overexpression of miR396b in danshen hairy roots inhibited hairy root growth, reduced salvianolic acid concentration, but enhanced tanshinone accumulation, resulting in the biomass and total salvianolic acids respectively reduced to 55.5 and 72.1% of the control and total tanshinones increased up to 1.91-fold of the control. Applied degradome sequencing, 5'RLM-RACE, and qRT-PCR, 13 targets for miR396b were identified including seven conserved SmGRF1-7 and six novel ones. Comparative transcriptomics and microRNomics analysis together with qRT-PCR results confirmed that miR396b targets SmGRFs, SmHDT1, and SmMYB37/4 to mediate the phytohormone, especially gibberellin signaling pathways and consequentially resulted in the phenotype variation of miR396b-OE hairy roots. Furthermore, miR396b could be activated by methyl jasmonate, abscisic acid, gibberellin, salt, and drought stresses. The findings in this study indicated that smi-miR396b acts as an upstream and central regulator in cell growth and the biosynthesis of tanshinones and salvianolic acids, shedding light on the coordinated regulation of plant growth and biosynthesis of active ingredients in S. miltiorrhiza.

Characterizing the Role of the miR156-SPL Network in Plant Development and Stress Response

MicroRNA (miRNA) is a short, single-stranded, non-coding RNA found in eukaryotic cells that can regulate the expression of many genes at the post-transcriptional level. Among various plant miRNAs with diverse functions, miR156 plays a key role in biological processes, including developmental regulation, immune response, metabolic regulation, and abiotic stress. MiRNAs have become the regulatory center for plant growth and development. MicroRNA156 (miR156) is a highly conserved and emerging tool for the improvement of plant traits, including crop productivity and stress tolerance. Fine-tuning of squamosa promoter biding-like (SPL) gene expression might be a useful strategy for crop improvement. Here, we studied the regulation of the miR156 module and its interaction with SPL factors to understand the developmental transition of various plant species. Furthermore, this review provides a strong background for plant biotechnology and is an important source of information for further molecular breeding to optimize farming productivity.

Hybrid transcriptome sequencing approach improved assembly and gene annotation in Cynara cardunculus (L.)

Background

The investigation of transcriptome profiles using short reads in non-model organisms, which lack of well-annotated genomes, is limited by partial gene reconstruction and isoform detection. In contrast, long-reads sequencing techniques revealed their potential to generate complete transcript assemblies even when a reference genome is lacking. Cynara cardunculus var. altilis (DC) (cultivated cardoon) is a perennial hardy crop adapted to dry environments with many industrial and nutraceutical applications due to the richness of secondary metabolites mostly produced in flower heads. The investigation of this species benefited from the recent release of a draft genome, but the transcriptome profile during the capitula formation still remains unexplored. In the present study we show a transcriptome analysis of vegetative and inflorescence organs of cultivated cardoon through a novel hybrid RNA-seq assembly approach utilizing both long and short RNA-seq reads.

Results

The inclusion of a single Nanopore flow-cell output in a hybrid sequencing approach determined an increase of 15% complete assembled genes and 18% transcript isoforms respect to short reads alone. Among 25,463 assembled unigenes, we identified 578 new genes and updated 13,039 gene models, 11,169 of which were alternatively spliced isoforms. During capitulum development, 3424 genes were differentially expressed and approximately two-thirds were identified as transcription factors including bHLH, MYB, NAC, C2H2 and MADS-box which were highly expressed especially after capitulum opening. We also show the expression dynamics of key genes involved in the production of valuable secondary metabolites of which capitulum is rich such as phenylpropanoids, flavonoids and sesquiterpene lactones. Most of their biosynthetic genes were strongly transcribed in the flower heads with alternative isoforms exhibiting differentially expression levels across the tissues.

Conclusions

This novel hybrid sequencing approach allowed to improve the transcriptome assembly, to update more than half of annotated genes and to identify many novel genes and different alternatively spliced isoforms. This study provides new insights on the flowering cycle in an Asteraceae plant, a valuable resource for plant biology and breeding in Cynara and an effective method for improving gene annotation.

Homology modeling and molecular dynamics based insights into Chalcone synthase and Chalcone isomerase in Phyllanthus emblica L

Chalcone synthase (CHS) and chalcone isomerase (CHI) plays a major role in the biosynthesis of flavonoid in plants. In this study, we made extensive bioinformatics analysis to gain functional and structural insight into PeCHS and PeCHI proteins. The phylogenetic distribution of PeCHS and PeCHI genes encoding proteins demonstrated the close evolutionary relationship with different CHS and CHI proteins of other dicot plants. MicroRNA target analysis showed miR169n and 3p miR5053 targeting PeCHS gene while miR169c-3p and miR4248 are targeting PeCHI gene, respectively. Three-dimensional structural models of PeCHS and PeCHI proteins were elucidated by homology modeling with Ramachandran plots showing the excellent geometry of the proteins structure. Molecular docking revealed that cinnamoyl-coa and naringenin chalcone substrates are strongly bound to PeCHS and PeCHI proteins, respectively. Finally, molecular dynamics (MD) simulation for 30 ns, further yielded stability checks of ligands in the binding pocket and behavior of protein complexes. Thus MD simulation and interaction fraction analysis showed the stable conformation of PeCHS and PeCHI proteins with their respective substrates during theee simulation. Our study provides first-hand structural prospective of PeCHS and PeCHI proteins towards understanding the mechanism of flavonoid biosynthetic pathway in P. emblica.

miRNAs involved in transcriptome remodeling during pollen development and heat stress response in Solanum lycopersicum

Cellular transitions during development and stress response depend on coordinated transcriptomic and proteomic alterations. Pollen is particular because its development is a complex process that includes meiotic and mitotic divisions which causes a high heat sensitivity of these cells. Development and stress response are accompanied by a reprogramming of the transcriptome, e.g. by post-transcriptional regulation via miRNAs. We identified known and potentially novel miRNAs in the transcriptome of developing and heat-stressed pollen of Solanum lycopersicum (tomato). The prediction of target mRNAs yielded an equal number of predicted target-sites in CDS and 3'UTR regions of target mRNAs. The result enabled the postulation of a possible link between miRNAs and a fine-tuning of transcription factor abundance during pollen development. miRNAs seem to play a role in the pollen heat stress response as well. We identified several heat stress transcription factors and heat shock proteins as putative targets of miRNAs in response to heat stress, thereby placing these miRNAs as important elements of thermotolerance. Moreover, for members of the AP2, SBP and ARF family members we could predict a miRNA-mediated regulation during development via the miR172, mir156 and mir160-family strengthening the current concept of a cross-connection between development and stress response in plants.

Transcriptional regulation and expression network responding to cadmium stress in a Cd-tolerant perennial grass Poa Pratensis

Kentucky bluegrass has good capability to absorb and accumulate cadmium (Cd) through developed root system, thus having potential phytoremediation function in Cd contaminated soils. Understanding the molecular mechanisms of Cd tolerance and accumulation in this species will be crucial to generating novel Cd-tolerance cultivars through genetic improvement, while it has not well documented yet. In the present study, comparative transcriptome analysis was performed for the seedlings of high Cd-tolerant genotype (M) and low Cd-tolerant genotype (R) under Cd stress. A total of 7022 up-regulated and 1033 down-regulated transcripts were identified in M genotype, whereas, only 850 up-regulated and 846 down-regulated transcripts were detected in R. Further transcriptional regulation analysis in M genotype showed that Dof, MADS25, BBR-BPC, B3, bZIP23 and MYB30 might be the hub transcription factors in response to Cd stress due to the orchestrated multiple functional genes associated with carbohydrate, lipid and secondary metabolism, as well as signal transduction. Differential expressed genes involved in auxin, ethylene, brassinosteroid and ABA signalling formed signal transduction cascades, which interacted with hub transcription factors, thereby finally orchestrated the expression of multiple genes associated with cell wall and membrane stability, cell elongation and Cd tolerance, including IAAs, ARFs, SnRK2, PP2C, PIFs, BES1/BZR1, CCR, CAD, FATB, fabF and HACD. Additionally, post-transcriptional modification of CIPKs, MAPKs, WAXs, UBCs, and E3 ubiquitin ligases were identified and also involved in plant signalling pathways and abiotic resistance. The study could contribute to our understanding the transcriptional regulation and complex internal network associated with Cd tolerance in Kentucky bluegrass.

Identification of small RNAs during cold acclimation in Arabidopsis thaliana

BACKGROUND

Cold stress causes dynamic changes in gene expression that are partially caused by small non-coding RNAs since they regulate protein coding transcripts and act in epigenetic gene silencing pathways. Thus, a detailed analysis of transcriptional changes of small RNAs (sRNAs) belonging to all known sRNA classes such as microRNAs (miRNA) and small interfering RNA (siRNAs) in response to cold contributes to an understanding of cold-related transcriptome changes.

RESULT

We subjected A. thaliana plants to cold acclimation conditions (4 °C) and analyzed the sRNA transcriptomes after 3 h, 6 h and 2 d. We found 93 cold responsive differentially expressed miRNAs and only 14 of these were previously shown to be cold responsive. We performed miRNA target prediction for all differentially expressed miRNAs and a GO analysis revealed the overrepresentation of miRNA-targeted transcripts that code for proteins acting in transcriptional regulation. We also identified a large number of differentially expressed cis- and trans-nat-siRNAs, as well as sRNAs that are derived from long non-coding RNAs. By combining the results of sRNA and mRNA profiling with miRNA target predictions and publicly available information on transcription factors, we reconstructed a cold-specific, miRNA and transcription factor dependent gene regulatory network. We verified the validity of links in the network by testing its ability to predict target gene expression under cold acclimation.

CONCLUSION

In A. thaliana, miRNAs and sRNAs derived from cis- and trans-NAT gene pairs and sRNAs derived from lncRNAs play an important role in regulating gene expression in cold acclimation conditions. This study provides a fundamental database to deepen our knowledge and understanding of regulatory networks in cold acclimation.

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

Background

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

Results

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

Conclusions

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

The Biosynthesis of Main Taste Compounds Is Coordinately Regulated by miRNAs and Phytohormones in Tea Plant (Camellia sinensis)

Based on the abundance of taste compounds in leaves at different leaf positions on the same shoot, green tea made from one bud and one leaf, or even just one bud, has the best quality. To elucidate the mechanism underlying the regulation of the biosynthesis of these compounds, we profiled the metabolome, transcriptome, sRNA, degradome, and WGCNA using leaves from five leaf positions of shoots. Through this analysis, we found 139 miRNA-target pairs related to taste compound biosynthesis and 96 miRNA-target pairs involved in phytohormone synthesis or signal transduction. Moreover, miR166-HD-ZIP, miR169-NF-YA, IAA, ZA, ABA, and JA were positively related to the accumulation of gallated catechin, caffeine, and theanine. However, miR396-GRF, miR393-bHLH, miR156-SBP, and SA were negatively correlated with these compounds. Among these important pairs, the miR396-GRF and miR156-SBP pairs were further validated by using qRT-PCR, Northern blots, and cotransformation. This is the first report describing that miRNA-TF pairs and phytohormones might synergistically regulate the biosynthesis of taste compounds in the leaves of tea plants.

PmiREN: a comprehensive encyclopedia of plant miRNAs

MicroRNAs (miRNAs) are small non-coding RNA molecules that function as diverse endogenous gene regulators at the post-transcriptional level. In the past two decades, as research effort on miRNA identification, function and evolution has soared, so has the demand for miRNA databases. However, the current plant miRNA databases suffer from several typical drawbacks, including a lack of entries for many important species, uneven annotation standards across different species, abundant questionable entries, and limited annotation. To address these issues, we developed a knowledge-based database called Plant miRNA Encyclopedia (PmiREN, http://www.pmiren.com/), which was based on uniform processing of sequenced small RNA libraries using miRDeep-P2, followed by manual curation using newly updated plant miRNA identification criteria, and comprehensive annotation. PmiREN currently contains 16,422 high confidence novel miRNA loci in 88 plant species and 3,966 retrieved from miRBase. For every miRNA entry, information on precursor sequence, precursor secondary structure, expression pattern, clusters and synteny in the genome, potential targets supported by Parallel Analysis of RNA Ends (PARE) sequencing, and references is attached whenever possible. PmiREN is hierarchically accessible and has eight built-in search engines. We believe PmiREN is useful for plant miRNA cataloguing and data mining, therefore a resource for data-driven miRNA research in plants.

Early Drought-Responsive Genes Are Variable and Relevant to Drought Tolerance

Drought stress is an important crop yield limiting factor worldwide. Plant physiological responses to drought stress are driven by changes in gene expression. While drought-responsive genes (DRGs) have been identified in maize, regulation patterns of gene expression during progressive water deficits remain to be elucidated. In this study, we generated time-series transcriptomic data from the maize inbred line B73 under well-watered and drought conditions. Comparisons between the two conditions identified 8,626 DRGs and the stages (early, middle, and late drought) at which DRGs occurred. Different functional groups of genes were regulated at the three stages. Specifically, early and middle DRGs display higher copy number variation among diverse Zea mays lines, and they exhibited stronger associations with drought tolerance as compared to late DRGs. In addition, correlation of expression between small RNAs (sRNAs) and DRGs from the same samples identified 201 negatively sRNA/DRG correlated pairs, including genes showing high levels of association with drought tolerance, such as two glutamine synthetase genes, gln2 and gln6 The characterization of dynamic gene responses to progressive drought stresses indicates important adaptive roles of early and middle DRGs, as well as roles played by sRNAs in gene expression regulation upon drought stress.

Plant stress biology in epigenomic era

Recent progress in "omics" methodologies allow us to gain insight into the complex molecular regulatory networks underlying plant responses to environmental stresses. Among the different genome-wide analysis, epigenomics is the most under-investigated "omic" approach requiring more critical and speculative discussion about approaches, methods and experimental designs. Epigenomics allows us to gain insight into the molecular adaptation of plants in response to environmental stresses. The identification of epigenetic marks transmitted during filial generations enables new theories to be developed on the evolution of living organisms in relation to environmental changes. The molecular mechanisms driving the capacity of plants to memorize a stress and to generate stress-resistant progenies are still unclear and scarcely investigated. The elucidation of these cryptic molecular switches will assist breeders in designing crops characterized by minimally compromised productivity in relation to stresses caused by climate change. The aim of this review is to briefly describe the most uptodate epigenomic approaches, update recent progresses in crop epigenomics in plant stress biology, and to stimulate the discussion of new epigenomic methods and approaches in the new era of "omic" sciences.

Identification of MicroRNAs and Their Targets That Respond to Powdery Mildew Infection in Cucumber by Small RNA and Degradome Sequencing

Powdery mildew (PM) is a prevalent disease known to limit cucumber production worldwide. MicroRNAs (miRNAs) are single-stranded molecules that regulate host defense responses through posttranscriptional gene regulation. However, which specific miRNAs are involved and how they regulate cucumber PM resistance remain elusive. A PM-resistant single-segment substitution line, SSSL508-28, was developed previously using marker-assisted backcrossing of the PM-susceptible cucumber inbred D8 line. In this study, we applied small RNA and degradome sequencing to identify PM-responsive miRNAs and their target genes in the D8 and SSSL508-28 lines. The deep sequencing resulted in the identification of 156 known and 147 novel miRNAs. Among them, 32 and six differentially expressed miRNAs (DEMs) were detected in D8 and SSSL508-28, respectively. The positive correlation between DEMs measured by small RNA sequencing and stem-loop quantitative real-time reverse transcription-polymerase chain reaction confirmed the accuracy of the observed miRNA abundances. The 32 DEMs identified in the PM-susceptible D8 were all upregulated, whereas four of the six DEMs identified in the PM-resistant SSSL508-28 were downregulated. Using in silico and degradome sequencing approaches, 517 and 20 target genes were predicted for the D8 and SSSL508-28 DEMs, respectively. Comparison of the DEM expression profiles with the corresponding mRNA expression profiles obtained in a previous study with the same experimental design identified 60 and three target genes in D8 and SSSL508-28, respectively, which exhibited inverse expression patterns with their respective miRNAs. In particular, five DEMs were located in the substituted segment that contained two upregulated DEMs, Csa-miR172c-3p and Csa-miR395a-3p, in D8 and two downregulated DEMs, Csa-miR395d-3p and Csa-miR398b-3p, in SSSL508-28. One gene encoding L-aspartate oxidase, which was targeted by Csa-miR162a, was also located on the same segment and was specifically downregulated in PM-inoculated D8 leaves. Our results will facilitate the future use of miRNAs in breeding cucumber varieties with enhanced resistance to PM.

Dynamic changes in the expression pattern of miRNAs and associated target genes during coconut somatic embryogenesis

Genome-wide analysis of small RNAs identifies somatic embryogenesis- specific miRNAs and their targets and provides novel insights into the mechanisms governing somatic embryogenesis in coconut, a highly in vitro recalcitrant species. Coconut, a major plantation crop of the tropics is recalcitrant to in vitro culture with a very low rate of somatic embryo turnover. Clonal propagation to enhance the production of high yielding, disease-free planting material in coconut has remained a distant reality. To better understand the molecular basis of this recalcitrance and to throw light on the complex regulatory network involved in the transition of coconut somatic cells to embryogenic calli, genome-wide profiling of small RNAs from embryogenic (EC) and non-embryogenic calli (NEC) was undertaken using Illumina Hiseq 2000 platform. We have identified a total of 110 conserved miRNAs (representing 46 known miRNA families) in both types of calli. In addition, 97 novel miRNAs (48 specific to EC, 21 specific to NEC and 28 common to both the libraries) were also identified. Among the conserved miRNAs, 10 were found to be differentially expressed between NEC and EC libraries with a log2 fold change > 2 following RPM-based normalization. miR156f, miR167c, miR169a, miR319a, miR535a, and miR5179 are upregulated and miR160a, miR166a, miR171a, and miR319b are down-regulated in NEC. To confirm the differential expression pattern and their regulatory role in SE, the expression patterns of miRNAs and their putative targets were analyzed using qRT- PCR and most of the analyzed miRNA-target pairs showed inverse correlation during somatic embryogenesis. Selected targets were further validated by RNA ligase mediated rapid amplification of 5' cDNA ends (5'RLM-RACE). Our data suggest that a few conserved miRNAs and species-specific miRNAs act in concert to regulate the process of somatic embryogenesis in coconut. The results of this study provide the first overview into the regulatory landscape of somatic embryogenesis in coconut and possible strategies for fine-tuning or reprogramming to enhance somatic embryo turn over in coconut.

Diverse and dynamic roles of F-box proteins in plant biology

The SCF complex is a widely studied multi-subunit ring E3 ubiquitin ligase that tags targeted proteins with ubiquitin for protein degradation by the ubiquitin 26S-proteasome system (UPS). The UPS is an important system that generally keeps cellular events tightly regulated by purging misfolded or damaged proteins and selectively degrading important regulatory proteins. The specificity of this post-translational regulation is controlled by F-box proteins (FBPs) via selective recognition of a protein-protein interaction motif at the C-terminal domain. Hence, FBPs are pivotal proteins in determining the plant response in multiple scenarios. It is not surprising that the FBP family is one of the largest protein families in the plant kingdom. In this review, the roles of FBPs, specifically in plants, are compiled to provide insights into their involvement in secondary metabolites, plant stresses, phytohormone signalling, plant developmental processes and miRNA biogenesis.

Computational screening of miRNAs and their targets in leaves of Hypericum spp. by transcriptome-mining: a pilot study

Our work provides a survey of mature miRNAs, their target genes and primary precursors identified by in-silico approach in leaf transcriptomes of five selected Hypericum species. MiRNAs are small non-coding RNA molecules found in animals, terrestrial plants, several algae and molds. As their role lies in the post-transcriptional gene silencing, these tiny molecules regulate many biological processes. Phyto-miRNAs are considered the important regulators of secondary metabolism in medicinal plants. The genus Hypericum comprises many producers of bioactive compounds, mainly unique naphtodianthrones with a great therapeutic potential. The main goal of our work was to identify genetically conserved miRNAs, characterize their primary precursors and target sequences in the leaf transcriptomes of five Hypericum species using in-silico approach. We found 20 sequences of potential Hypericum pri-miRNAs, and predicted and computationally validated their secondary structures. The mature miRNAs were identified by target genes screening analysis. Whereas predicted miRNA profiles differed in less genetically conserved families, the highly conserved miRNAs were found in almost all studied species. Moreover, we detected several novel highly likely miRNA-mRNA interactions, such as mir1171 with predicted regulatory role in the biosynthesis of melatonin in plants. Our work contributes to the knowledge of Hypericum miRNAome and miRNA-mRNA interactions.

Small RNA Bidirectional Crosstalk During the Interaction Between Wheat and Zymoseptoria tritici

Cross-kingdom RNA interference (RNAi) has been shown to play important roles during plant-pathogen interactions, and both plants and pathogens can use small RNAs (sRNAs) to silence genes in each other. This bidirectional cross-kingdom RNAi was still unexplored in the wheat-Zymoseptoria tritici pathosystem. Here, we performed a detailed analysis of the sRNA bidirectional crosstalk between wheat and Z. tritici. Using a combination of small RNA sequencing (sRNA-seq) and microRNA sequencing (mRNA-seq), we were able to identify known and novel sRNAs and study their expression and their action on putative targets in both wheat and Z. tritici. We predicted the target genes of all the sRNAs in either wheat or Z. tritici transcriptome and used degradome analysis to validate the cleavage of these gene transcripts. We could not find any clear evidence of a cross-kingdom RNAi acting by mRNA cleavage in this pathosystem. We also found that the fungal sRNA enrichment was lower in planta than during in vitro growth, probably due to the lower expression of the only Dicer gene of the fungus during plant infection. Our results support the recent finding that Z. tritici sRNAs cannot play important roles during wheat infection. However, we also found that the fungal infection induced wheat sRNAs regulating the expression of specific wheat genes, including auxin-related genes, as an immune response. These results indicate a role of sRNAs in the regulation of wheat defenses during Z. tritici infection. Our findings contribute to improve our understanding of the interactions between wheat and Z. tritici.

Genome-Wide Identification of Cassava Serine/Arginine-Rich Proteins: Insights into Alternative Splicing of Pre-mRNAs and Response to Abiotic Stress

Serine/arginine-rich (SR) proteins have an essential role in the splicing of pre-messenger RNA (pre-mRNA) in eukaryote. Pre-mRNA with introns can be alternatively spliced to generate multiple transcripts, thereby increasing adaptation to the external stress conditions in planta. However, pre-mRNA of SR proteins can also be alternatively spliced in different plant tissues and in response to diverse stress treatments, indicating that SR proteins might be involved in regulating plant development and adaptation to environmental changes. We identified and named 18 SR proteins in cassava and systematically studied their splicing and transcriptional changes under tissue-specific and abiotic stress conditions. Fifteen out of 18 SR genes showed alternative splicing in the tissues. 45 transcripts were found from 18 SR genes under normal conditions, whereas 55 transcripts were identified, and 21 transcripts were alternate spliced in some SR genes under salt stress, suggesting that SR proteins might participate in the plant adaptation to salt stress. We then found that overexpression of MeSR34 in Arabidopsis enhanced the tolerance to salt stress through maintaining reactive oxygen species homeostasis and increasing the expression of calcineurin B-like proteins (CBL)-CBL-interacting protein kinases and osmotic stress-related genes. Therefore, our findings highlight the critical role of cassava SR proteins as regulators of RNA splicing and salt tolerance in planta.

Sugar Signaling and Post-transcriptional Regulation in Plants: An Overlooked or an Emerging Topic?

Plants are autotrophic organisms that self-produce sugars through photosynthesis. These sugars serve as an energy source, carbon skeletons, and signaling entities throughout plants' life. Post-transcriptional regulation of gene expression plays an important role in various sugar-related processes. In cells, it is regulated by many factors, such as RNA-binding proteins (RBPs), microRNAs, the spliceosome, etc. To date, most of the investigations into sugar-related gene expression have been focused on the transcriptional level in plants, while only a few studies have been conducted on post-transcriptional mechanisms. The present review provides an overview of the relationships between sugar and post-transcriptional regulation in plants. It addresses the relationships between sugar signaling and RBPs, microRNAs, and mRNA stability. These new items insights will help to reach a comprehensive understanding of the diversity of sugar signaling regulatory networks, and open onto new investigations into the relevance of these regulations for plant growth and development.

Dynamic architecture and regulatory implications of the miRNA network underlying the response to stress in melon

miRNAs are small RNAs that regulate mRNAs at both transcriptional and posttranscriptional level. In plants, miRNAs are involved in the regulation of different processes including development and stress-response. Elucidating how stress-responsive miRNAs are regulated is key to understand the global response to stress but also to develop efficient biotechnological tools that could help to cope with stress. Here, we describe a computational approach based on sRNA sequencing, transcript quantification and degradome data to analyse the accumulation, function and structural organization of melon miRNAs reactivated under seven biotic and abiotic stress conditions at two and four days post-treatment. Our pipeline allowed us to identify fourteen stress-responsive miRNAs (including evolutionary conserved such as miR156, miR166, miR172, miR319, miR398, miR399, miR894 and miR408) at both analysed times. According to our analysis miRNAs were categorized in three groups showing a broad-, intermediate- or narrow- response range. miRNAs reactive to a broad range of environmental cues appear as central components in the stress-response network. The strictly coordinated response of miR398 and miR408 (broad response-range) to the seven stress treatments during the period analysed here reinforces this notion. Although both, the amplitude and diversity of the miRNA-related response to stress changes during the exposition time, the architecture of the miRNA-network is conserved. This organization of miRNA response to stress is also conserved in rice and soybean supporting the conservation of miRNA-network organization in other crops. Overall, our work sheds light into how miRNA networks in plants organize and function during stress.

Diurnal Regulation of Plant Epidermal Wax Synthesis through Antagonistic Roles of the Transcription Factors SPL9 and DEWAX

Plant surface waxes form an outer barrier that protects the plant from many forms of environmental stress. The deposition of cuticular waxes on the plant surface is regulated by external environmental changes, including light and dark cycles. However, the underlying molecular mechanisms controlling light regulation of wax production are still poorly understood, especially at the posttranscriptional level. In this paper, we report the regulation of cuticular wax production by the miR156-SPL9 (SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 9) module in Arabidopsis (Arabidopsis thaliana). When compared with wild-type plants, miR156 and SPL9 mutants showed significantly altered cuticular wax amounts in both stems and leaves. Furthermore, it was found that SPL9 positively regulates gene expression of the alkane-forming enzyme ECERIFERUM1 (CER1), as well as the primary (1-) alcohol-forming enzyme ECERIFERUM4 (CER4), to enhance alkane and 1-alcohol synthesis, respectively. Our results indicate that complex formation of SPL9 with a negative regulator of wax synthesis, DEWAX, will hamper SPL9 DNA binding ability, possibly by interfering with SPL9 homodimerization. Combined with their diurnal gene and protein expressions, this dynamic repression-activation transcriptional module defines a dynamic mechanism that may allow plants to optimize wax synthesis during daily cycles. These findings provide a regulatory framework for environmental signal integration in the regulation of wax synthesis.

Plant microRNAs: biogenesis, gene silencing, web-based analysis tools and their use as molecular markers

MicroRNAs (miRNAs) are tiny (20-24 nt bp) regulatory non-protein-coding RNA molecules that have been extensively characterized and found important for many physiological and developmental processes. The miss-expression of miRNAs leads to various defects in plants. MicroRNAs repress gene expression by directing mRNA degradation or translational arrest. Several proteins such as PP43A, HYL1, DCL, HST are indispensable role players in promoting miRNA biogenesis in plants. During miRNA biogenesis, lariat RNAs are produced as by-products of pre-mRNA splicing which have a negative role in regulation of miRNA homeostasis. By acting as a decoy and by sequestering to the dicing complex, lariat RNA can prevent the processing of miRNAs. A number of bioinformatic tools with different methodologies are available to identify and validate miRNAs and their targets. Many miRNAs have been reported in different crops for different traits; however, no reports are available on their use in plant breeding. Recently, researchers have developed trait specific miRNA-based molecular markers (miRNA-SSRs/SNP) for many quantitative traits in different plant species. In the future, these molecular markers can be used for plant breeding programs. In this review, a comprehensive up-to-date information is provided on the bioinformatic tools used for analysis of plant miRNAs and their targets, the number of miRNAs, their biogenesis, gene silencing mechanism and miRNA-based molecular markers.

Profiling miRNA expression in photo-thermo-sensitive male genic sterility line (PTGMS) PA64S under high and low temperature

Peiai64S (PA64S) is a photo-thermo-sensitive genic male sterile line (PTGMS), with wide application in hybrid seed production in rice (Oryza sativa L.). Micro-RNAs are 21-24 nt, endogenously expressed small RNAs that have been characterized in various developmental stages of rice, but none have been studied with respect to the regulation of TGMS in rice. Here, we employed high-throughput sequencing to identify expression profiles of miRNAs in the anthers of PA64S at high (PA64S-H) and low temperature (PA64S-L). Two small RNA libraries from PA64S-H and PA64-L anthers were sequenced, and 263 known and 321 novel candidate miRNAs were identified. Based on the number of sequencing reads, a total of 133 known miRNAs were found to be differentially expressed between PA64S-H and PA64S-L. Target prediction showed that the target genes encode MYB and TCP transcription factors, and bHLH proteins. These target genes are related to pollen development and male sterility, suggesting that miRNA/targets may play roles in regulating TGMS in rice. Further, starch and sucrose metabolism pathways, sphingolipid metabolism, arginine and proline metabolism, and plant hormone signal transduction pathways were enriched by KEGG pathway annotation. These findings contribute to our understanding of the role of miRNAs during anther development and TGMS occurrence in rice.

Transcription Factors Associated with Abiotic and Biotic Stress Tolerance and Their Potential for Crops Improvement

In field conditions, crops are adversely affected by a wide range of abiotic stresses including drought, cold, salt, and heat, as well as biotic stresses including pests and pathogens. These stresses can have a marked effect on crop yield. The present and future effects of climate change necessitate the improvement of crop stress tolerance. Plants have evolved sophisticated stress response strategies, and genes that encode transcription factors (TFs) that are master regulators of stress-responsive genes are excellent candidates for crop improvement. Related examples in recent studies include TF gene modulation and overexpression approaches in crop species to enhance stress tolerance. However, much remains to be discovered about the diverse plant TFs. Of the >80 TF families, only a few, such as NAC, MYB, WRKY, bZIP, and ERF/DREB, with vital roles in abiotic and biotic stress responses have been intensively studied. Moreover, although significant progress has been made in deciphering the roles of TFs in important cereal crops, fewer TF genes have been elucidated in sorghum. As a model drought-tolerant crop, sorghum research warrants further focus. This review summarizes recent progress on major TF families associated with abiotic and biotic stress tolerance and their potential for crop improvement, particularly in sorghum. Other TF families and non-coding RNAs that regulate gene expression are discussed briefly. Despite the emphasis on sorghum, numerous examples from wheat, rice, maize, and barley are included. Collectively, the aim of this review is to illustrate the potential application of TF genes for stress tolerance improvement and the engineering of resistant crops, with an emphasis on sorghum.

Roles of the INO80 and SWR1 Chromatin Remodeling Complexes in Plants

Eukaryotic genes are packed into a dynamic but stable nucleoprotein structure called chromatin. Chromatin-remodeling and modifying complexes generate a dynamic chromatin environment that ensures appropriate DNA processing and metabolism in various processes such as gene expression, as well as DNA replication, repair, and recombination. The INO80 and SWR1 chromatin remodeling complexes (INO80-c and SWR1-c) are ATP-dependent complexes that modulate the incorporation of the histone variant H2A.Z into nucleosomes, which is a critical step in eukaryotic gene regulation. Although SWR1-c has been identified in plants, plant INO80-c has not been successfully isolated and characterized. In this review, we will focus on the functions of the SWR1-c and putative INO80-c (SWR1/INO80-c) multi-subunits and multifunctional complexes in Arabidopsis thaliana. We will describe the subunit compositions of the SWR1/INO80-c and the recent findings from the standpoint of each subunit and discuss their involvement in regulating development and environmental responses in Arabidopsis.

Drought Response in Rice: The miRNA Story

As a semi-aquatic plant, rice requires water for proper growth, development, and orientation of physiological processes. Stress is induced at the cellular and molecular level when rice is exposed to drought or periods of low water availability. Plants have existing defense mechanisms in planta that respond to stress. In this review we examine the role played by miRNAs in the regulation and control of drought stress in rice through a summary of molecular studies conducted on miRNAs with emphasis on their contribution to drought regulatory networks in comparison to other plant systems. The interaction between miRNAs, target genes, transcription factors and their respective roles in drought-induced stresses is elaborated. The cross talk involved in controlling drought stress responses through the up and down regulation of targets encoding regulatory and functional proteins is highlighted. The information contained herein can further be explored to identify targets for crop improvement in the future.