Stress-induced alternative splicing provides a mechanism for the regulation of microRNA processing in Arabidopsis thaliana

Iso-Seq Allows Genome-Independent Transcriptome Profiling of Grape Berry Development

Transcriptomics has been widely applied to study grape berry development. With few exceptions, transcriptomic studies in grape are performed using the available genome sequence, PN40024, as reference. However, differences in gene content among grape accessions, which contribute to phenotypic differences among cultivars, suggest that a single reference genome does not represent the species' entire gene space. Though whole genome assembly and annotation can reveal the relatively unique or "private" gene space of any particular cultivar, transcriptome reconstruction is a more rapid, less costly, and less computationally intensive strategy to accomplish the same goal. In this study, we used single molecule-real time sequencing (SMRT) to sequence full-length cDNA (Iso-Seq) and reconstruct the transcriptome of Cabernet Sauvignon berries during berry ripening. In addition, short reads from ripening berries were used to error-correct low-expression isoforms and to profile isoform expression. By comparing the annotated gene space of Cabernet Sauvignon to other grape cultivars, we demonstrate that the transcriptome reference built with Iso-Seq data represents most of the expressed genes in the grape berries and includes 1,501 cultivar-specific genes. Iso-Seq produced transcriptome profiles similar to those obtained after mapping on a complete genome reference. Together, these results justify the application of Iso-Seq to identify cultivar-specific genes and build a comprehensive reference for transcriptional profiling that circumvents the necessity of a genome reference with its associated costs and computational weight.

Cold-Dependent Expression and Alternative Splicing of Arabidopsis Long Non-coding RNAs

Plants re-program their gene expression when responding to changing environmental conditions. Besides differential gene expression, extensive alternative splicing (AS) of pre-mRNAs and changes in expression of long non-coding RNAs (lncRNAs) are associated with stress responses. RNA-sequencing of a diel time-series of the initial response of Arabidopsis thaliana rosettes to low temperature showed massive and rapid waves of both transcriptional and AS activity in protein-coding genes. We exploited the high diversity of transcript isoforms in AtRTD2 to examine regulation and post-transcriptional regulation of lncRNA gene expression in response to cold stress. We identified 135 lncRNA genes with cold-dependent differential expression (DE) and/or differential alternative splicing (DAS) of lncRNAs including natural antisense RNAs, sORF lncRNAs, and precursors of microRNAs (miRNAs) and trans-acting small-interfering RNAs (tasiRNAs). The high resolution (HR) of the time-series allowed the dynamics of changes in transcription and AS to be determined and identified early and adaptive transcriptional and AS changes in the cold response. Some lncRNA genes were regulated only at the level of AS and using plants grown at different temperatures and a HR time-course of the first 3 h of temperature reduction, we demonstrated that the AS of some lncRNAs is highly sensitive to small temperature changes suggesting tight regulation of expression. In particular, a splicing event in TAS1a which removed an intron that contained the miR173 processing and phased siRNAs generation sites was differentially alternatively spliced in response to cold. The cold-induced reduction of the spliced form of TAS1a and of the tasiRNAs suggests that splicing may enhance production of the siRNAs. Our results identify candidate lncRNAs that may contribute to the regulation of expression that determines the physiological processes essential for acclimation and freezing tolerance.

The Involvement of Long Noncoding RNAs in Response to Plant Stress

Plant growth and productivity are greatly impacted by environmental stresses. Therefore, plants have evolved mechanisms which allow them to adapt to abiotic stresses through alterations in gene expression and metabolism. In recent years, studies have investigated the role of long noncoding RNA (lncRNA) in regulating gene expression in plants and characterized their involvement in various biological functions through their regulation of DNA methylation, DNA structural modifications, histone modifications, and RNA-RNA interactions. Genome-wide transcriptome analyses have identified various types of noncoding RNAs (ncRNAs) that respond to abiotic stress. These ncRNAs are in addition to the well-known housekeeping ncRNAs, such as rRNAs, tRNAs, snoRNAs, and snRNAs. In this review, recent research pertaining to the role of lncRNAs in the response of plants to abiotic stress is summarized and discussed.

The Plant microRNA Pathway: The Production and Action Stages

Plant microRNAs are an endogenous class of small regulatory RNA central to the posttranscriptional regulation of gene expression in plant development and environmental stress adaptation or in response to pathogen challenge. The plant microRNA pathway is readily separated into two distinct stages: (1) the production stage, which is localized to the plant cell nucleus and where the microRNA small RNA is processed from a double-stranded RNA precursor transcript, and (2) the action stage, which is localized to the plant cell cytoplasm and where the mature microRNA small RNA is loaded into an effector complex and is used by the complex as a sequence specificity guide to direct expression repression of target genes harboring highly complementary microRNA target sequences. Historical research indicated that the plant microRNA pathway was a highly structured, almost linear pathway requiring a small set of core machinery proteins. However, contemporary research has demonstrated that the plant microRNA pathway is highly dynamic, and to allow for this flexibility, a large and highly functionally diverse set of machinery proteins is now known to be required. For example, recent research has shown that plant microRNAs can regulate target gene expression via a translational repression mechanism of RNA silencing in addition to the standard messenger RNA cleavage-based mechanism of RNA silencing: a mode of RNA silencing originally assigned to all plant microRNAs. Using Arabidopsis thaliana as our model system, here we report on both the core and auxiliary sets of machinery proteins now known to be required for both microRNA production and microRNA action in plants.

miRNomes involved in imparting thermotolerance to crop plants

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

Alternatively spliced variants in Atlantic cod (Gadus morhua) support response to variable salinity environment

Analysis of gill transcriptome of the Atlantic cod from the Baltic Sea demonstrated that alternatively spliced (AS) variants may be actively involved in the process of adaptation to altered salinity. Some AS variants of different genes, like phospholipase A2 group IVC (PLA2G4C), appeared only in fish exposed to altered salinity, while other isoforms of the same genes were present in all experimental groups. Novel sequence arrangements represent 89% of all AS in the Baltic cod compared to the Atlantic population. Profiles of modified pathways suggest that regulation by AS can afford specific changes of genes expressed in response to the environment. The AS variants appear to be involved in the response to stress by modifications of signalling in apoptosis pathways, an innate immunological response and pro-inflammatory process. Present results support the hypothesis that developing new AS variants could support genome complexity and reinforce the ability to fast adapt to local environments.

A deeper view into the significance of simple sequence repeats in pre-miRNAs provides clues for its possible roles in determining the function of microRNAs

Background

The central tenet of ‘genome content’ has been that the ‘non-coding’ parts are highly enriched with ‘microsatellites’ or ‘Simple Sequence Repeats’ (SSRs). We presume that the presence and change in number of repeat unit (n) of SSRs in different genomic locations may or may not become beneficial, depending on the position of SSRs in a gene. Very few studies have looked into the existence of SSRs in the hair-pin precursors of miRNAs (pre-miRNAs). The interplay between SSRs and miRNAs is not yet clearly understood.

Results

Considering the potential significance of SSRs in pre-miRNAs, we analysed the miRNA hair-pin precursors of 171 organisms, which revealed a noticeable (29.8%) existence of SSRs in their pre-miRNAs. The maintenance of SSRs in pre-miRNAs even in the complex, highly evolved phyla like Chordata and Magnoliophyta shed light upon its diverse functions. Putative effects of SSRs in either regulating the biogenesis or function of miRNAs were more underlined based on computational and experimental analysis. A preliminary computational analysis to explore the relevance of such SSRs maintained in pre-miRNA sequences led to the detection of splicing regulatory elements (SREs) either in or near to the SSRs. The absence of SSRs correspondingly decreased the detection of SREs.

Conclusion

The present study is the first implication for the possible involvement of SSRs in shaping the SREs to undergo Alternative Splicing events to produce miRNA isoforms in accordance with different stress environments. This part of work well demonstrates the importance of studying such consistently maintained SSRs residing in pre-miRNAs and can enhance more and more research towards deciphering the exact function of SSRs in the near future.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0615-x) contains supplementary material, which is available to authorized users.

Thermopriming triggers splicing memory in Arabidopsis

Abiotic and biotic stresses limit crop productivity. Exposure to a non-lethal stress, referred to as priming, can allow plants to survive subsequent and otherwise lethal conditions; the priming effect persists even after a prolonged stress-free period. However, the molecular mechanisms underlying priming are not fully understood. Here, we investigated the molecular basis of heat-shock memory and the role of priming in Arabidopsis thaliana. Comprehensive analysis of transcriptome-wide changes in gene expression and alternative splicing in primed and non-primed plants revealed that alternative splicing functions as a novel component of heat-shock memory. We show that priming of plants with a non-lethal heat stress results in de-repression of splicing after a second exposure to heat stress. By contrast, non-primed plants showed significant repression of splicing. These observations link 'splicing memory' to the ability of plants to survive subsequent and otherwise lethal heat stress. This newly discovered priming-induced splicing memory may represent a general feature of heat-stress responses in plants and other organisms as many of the key components are conserved among eukaryotes. Furthermore, this finding could facilitate the development of novel approaches to improve plant survival under extreme heat stress.

Splicing regulation by long noncoding RNAs

Massive high-throughput sequencing techniques allowed the identification of thousands of noncoding RNAs (ncRNAs) and a plethora of different mRNA processing events occurring in higher organisms. Long ncRNAs can act directly as long transcripts or can be processed into active small si/miRNAs. They can modulate mRNA cleavage, translational repression or the epigenetic landscape of their target genes. Recently, certain long ncRNAs have been shown to play a crucial role in the regulation of alternative splicing in response to several stimuli or during disease. In this review, we focus on recent discoveries linking gene regulation by alternative splicing and its modulation by long and small ncRNAs.

UTR-Dependent Control of Gene Expression in Plants

Throughout their lives, plants sense many developmental and environmental stimuli, and activation of optimal responses against these stimuli requires extensive transcriptional reprogramming. To facilitate this activation, plant mRNA contains untranslated regions (UTRs) that significantly increase the coding capacity of the genome by producing multiple mRNA variants from the same gene. In this review we compare UTRs of arabidopsis (Arabidopsis thaliana) and rice (Oryza sativum) at the genome scale to highlight their complexity in crop plants. We discuss different modes of UTR-based regulation with emphasis on genes that regulate multiple plant processes, including flowering, stress responses, and nutrient homeostasis. We demonstrate functional specificity in genes with variable UTR length and propose future research directions.

Non-coding RNAs and plant male sterility: current knowledge and future prospects

Latest outcomes assign functional role to non-coding (nc) RNA molecules in regulatory networks that confer male sterility to plants. Male sterility in plants offers great opportunity for improving crop performance through application of hybrid technology. In this respect, cytoplasmic male sterility (CMS) and sterility induced by photoperiod (PGMS)/temperature (TGMS) have greatly facilitated development of high-yielding hybrids in crops. Participation of non-coding (nc) RNA molecules in plant reproductive development is increasingly becoming evident. Recent breakthroughs in rice definitively associate ncRNAs with PGMS and TGMS. In case of CMS, the exact mechanism through which the mitochondrial ORFs exert influence on the development of male gametophyte remains obscure in several crops. High-throughput sequencing has enabled genome-wide discovery and validation of these regulatory molecules and their target genes, describing their potential roles performed in relation to CMS. Discovery of ncRNA localized in plant mtDNA with its possible implication in CMS induction is intriguing in this respect. Still, conclusive evidences linking ncRNA with CMS phenotypes are currently unavailable, demanding complementing genetic approaches like transgenics to substantiate the preliminary findings. Here, we review the recent literature on the contribution of ncRNAs in conferring male sterility to plants, with an emphasis on microRNAs. Also, we present a perspective on improved understanding about ncRNA-mediated regulatory pathways that control male sterility in plants. A refined understanding of plant male sterility would strengthen crop hybrid industry to deliver hybrids with improved performance.

Heat stress alters genome-wide profiles of circular RNAs in Arabidopsis

1599 novel circRNAs and 1583 heat stress-specific circRNAs were identified in Arabidopsis. Heat stress enhanced accumulation of circRNAs remarkably. Heat stress altered the sizes of circRNAs, numbers of circularized exons and alterative circularization events. A putative circRNA-mediated ceRNA networks under heat stress was established. Heat stress retards plant growth and destabilizes crop yield. The noncoding RNAs were demonstrated to be involved in plant response to heat stress. As a newly-characterized class of noncoding RNAs, circular RNAs (circRNAs) play important roles in transcriptional and post-transcriptional regulation. A few recent investigations indicated that plant circRNAs were differentially expressed under abiotic stress. However, little is known about how heat stress mediates biogenesis of circRNAs in plants. Here, we uncovered 1599 previously-unknown circRNAs and 1583 heat-specific circRNAs, by RNA-sequencing and bioinformatic analysis. Our results indicated that much more circRNAs were expressed under heat stress than in control condition. Besides, heat stress also increased the length of circRNAs, the quantity of circularized exons, and alternative circularization events. Moreover, we observed a positive correlation between expression patterns of some circRNAs and their parental genes. The prediction of ceRNA (competing endogenous RNA) networks indicated that differentially-expressed circRNAs could influence expression of many important genes, that participate in response to heat stress, hydrogen peroxide, and phytohormone signaling pathways, by interacting with the corresponding microRNAs. Together, our observations indicated that heat stress had great impacts on the biogenesis of circRNAs. Heat-induced circRNAs might participate in plant response to heat stress through the circRNA-mediated ceRNA networks.

MicroRNA160 Modulates Plant Development and Heat Shock Protein Gene Expression to Mediate Heat Tolerance in Arabidopsis

Global warming is causing a negative impact on plant growth and adversely impacts on crop yield. MicroRNAs (miRNAs) are critical in regulating the expression of genes involved in plant development as well as defense responses. The effects of miRNAs on heat-stressed Arabidopsis warrants further investigation. Heat stress increased the expression of miR160 and its precursors but considerably reduced that of its targets, ARF10, ARF16, and ARF17. To study the roles of miR160 during heat stress, transgenic Arabidopsis plants overexpressing miR160 precursor a (160OE) and artificial miR160 (MIM160), which mimics an inhibitor of miR160, were created. T-DNA insertion mutants of miR160 targets were also used to examine their tolerances to heat stress. Results presented that overexpressing miR160 improved seed germination and seedling survival under heat stress. The lengths of hypocotyl elongation and rachis were also longer in 160OE than the wild-type (WT) plants under heat stress. Interestingly, MIM160 plants showed worse adaption to heat. In addition, arf10, arf16, and arf17 mutants presented similar phenotypes to 160OE under heat stress to advance abilities of thermotolerance. Moreover, transcriptome and qRT-PCR analyses revealed that HSP17.6A, HSP17.6II, HSP21, and HSP70B expression levels were regulated by heat in 160OE, MIM160, arf10, arf16, and arf17 plants. Hence, miR160 altered the expression of the heat shock proteins and plant development to allow plants to survive heat stress.

Functional intron-derived miRNAs and host-gene expression in plants

BACKGROUND

Recently, putative pre-miRNAs locations have been identified in the introns of plant genes, raising the question whether such genes can show a dual functionality by having both correct maturation of the host gene pre-mRNA and maturation of the miRNAs from the intron. Here, we demonstrated that such dual functionality is indeed possible, using as host gene the firefly luciferase gene with intron (ffgLUC), and different artificial intronic miRNAs (aimiRNA) placed within the intron of ffgLUC.

RESULTS

The miRNAs were based on the structure of the natural miR319a. Luciferase (LUC) activity in planta was used to evaluate a correct splicing of the ffgLUC mRNA. Different target sequences were inserted into the aimiRNA to monitor efficiency of silencing of different target mRNAs. After adjusting the insertion cloning strategy, the ffgLUC^aimiR-319a gene showed dual functionality with correct splicing of ffgLUC and efficient silencing of TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR1 transcription factor genes targeted in-trans by aimiR-319a or targeting the transgene ffLUC in-cis by an aimiR-LUC. Silencing of endogenous target genes by aimiRNA or amiRNA is efficient both in transient assays and stable transformants. A behave as strong phenotype the PHYTOCHROME B (PHYB) gene was also targeted by ffgLUC^aimiR-PHYB. The lack of silencing of the PHYB target was most likely due to an insensitive target site within the PHYB mRNA which can potentially form a double stranded stem structure.

CONCLUSION

The combination of an overexpression construct with an artificial intronic microRNA allows for a simultaneous dual function in plants. The concept therefore adds new options to engineering of plant traits that require multiple gene manipulations.

A global survey of alternative splicing in allopolyploid cotton: landscape, complexity and regulation

Alternative splicing (AS) is a crucial regulatory mechanism in eukaryotes, which acts by greatly increasing transcriptome diversity. The extent and complexity of AS has been revealed in model plants using high-throughput next-generation sequencing. However, this technique is less effective in accurately identifying transcript isoforms in polyploid species because of the high sequence similarity between coexisting subgenomes. Here we characterize AS in the polyploid species cotton. Using Pacific Biosciences single-molecule long-read isoform sequencing (Iso-Seq), we developed an integrated pipeline for Iso-Seq transcriptome data analysis (https://github.com/Nextomics/pipeline-for-isoseq). We identified 176 849 full-length transcript isoforms from 44 968 gene models and updated gene annotation. These data led us to identify 15 102 fibre-specific AS events and estimate that c. 51.4% of homoeologous genes produce divergent isoforms in each subgenome. We reveal that AS allows differential regulation of the same gene by miRNAs at the isoform level. We also show that nucleosome occupancy and DNA methylation play a role in defining exons at the chromatin level. This study provides new insights into the complexity and regulation of AS, and will enhance our understanding of AS in polyploid species. Our methodology for Iso-Seq data analysis will be a useful reference for the study of AS in other species.

Small RNAs in regulating temperature stress response in plants

Due to global climate change, temperature stress has become one of the primary causes of crop losses worldwide. Much progress has been made in unraveling the complex stress response mechanisms in plants, particularly in the identification of temperature stress responsive protein-coding genes. Recently discovered microRNAs (miRNAs) and endogenous small-interfering RNAs (siRN As) have also been demonstrated as important players in plant temperature stress response. Using high-throughput sequencing, many small RNAs, especially miRNAs, have been identified to be triggered by cold or heat. Subsequently, several studies have shown an important functional role for these small RNAs in cold or heat tolerance. These findings greatly broaden our understanding of endogenous small RNAs in plant stress response control. Here, we highlight new findings regarding the roles of miRNAs and siRNAs in plant temperature stress response and acclimation. We also review the current understanding of the regulatory mechanisms of small RNAs in temperature stress response, and explore the outlook for the use of these small RNAs in molecular breeding for improvement of temperature stress tolerance in plants.

Extensive transcriptomic and epigenomic remodelling occurs during Arabidopsis thaliana germination

BACKGROUND

Seed germination involves progression from complete metabolic dormancy to a highly active, growing seedling. Many factors regulate germination and these interact extensively, forming a complex network of inputs that control the seed-to-seedling transition. Our understanding of the direct regulation of gene expression and the dynamic changes in the epigenome and small RNAs during germination is limited. The interactions between genome, transcriptome and epigenome must be revealed in order to identify the regulatory mechanisms that control seed germination.

RESULTS

We present an integrated analysis of high-resolution RNA sequencing, small RNA sequencing and MethylC sequencing over ten developmental time points in Arabidopsis thaliana seeds, finding extensive transcriptomic and epigenomic transformations associated with seed germination. We identify previously unannotated loci from which messenger RNAs are expressed transiently during germination and find widespread alternative splicing and divergent isoform abundance of genes involved in RNA processing and splicing. We generate the first dynamic transcription factor network model of germination, identifying known and novel regulatory factors. Expression of both microRNA and short interfering RNA loci changes significantly during germination, particularly between the seed and the post-germinative seedling. These are associated with changes in gene expression and large-scale demethylation observed towards the end of germination, as the epigenome transitions from an embryo-like to a vegetative seedling state.

CONCLUSIONS

This study reveals the complex dynamics and interactions of the transcriptome and epigenome during seed germination, including the extensive remodelling of the seed DNA methylome from an embryo-like to vegetative-like state during the seed-to-seedling transition. Data are available for exploration in a user-friendly browser at https://jbrowse.latrobe.edu.au/germination_epigenome .

Molecular Basis Underlying Leaf Variegation of a Moth Orchid Mutant (Phalaenopsis aphrodite subsp. formosana)

Leaf variegation is often the focus of plant breeding. Here, we studied a variegated mutant of Phalaenopsis aphrodite subsp. formosana, which is usually used as a parent of horticultural breeding, to understand its anatomic and genetic regulatory mechanisms in variegation. Chloroplasts with well-organized thylakoids and starch grains were found only in the mesophyll cells of green sectors but not of yellow sectors, confirming that the variegation belongs to the chlorophyll type. The two-dimensional electrophoresis and LC/MS/MS also reveal differential expressions of PsbP and PsbO between the green and yellow leaf sectors. Full-length cDNA sequencing revealed that mutant transcripts were caused by intron retention. When conditioning on the total RNA expression, we found that the functional transcript of PsbO and mutant transcript of PsbP are higher expressed in the yellow sector than in the green sector, suggesting that the post-transcriptional regulation of PsbO and PsbP differentiates the performance between green and yellow sectors. Because PsbP plays an important role in the stability of thylakoid folding, we suggest that the negative regulation of PsbP may inhibit thylakoid development in the yellow sectors. This causes chlorophyll deficiency in the yellow sectors and results in leaf variegation. We also provide evidence of the link of virus CymMV and the formation of variegation according to the differential expression of CymMV between green and yellow sectors.

Genome-wide association mapping of latex yield and girth in Amazonian accessions of Hevea brasiliensis grown in a suboptimal climate zone

Latex yield and growth are the key complex traits in commercial rubber production. The present study is the first to report genome-wide association mapping of latex yield and girth, for 170 Amazonian accessions grown in a suboptimal area characterized by limited rainfall and a lengthy dry season. Targeted sequence enrichment to capture gene transcripts generated 14,155 high quality filtered single nucleotide polymorphisms (SNPs) of which 94.3% resided in coding regions. The rapid decay of linkage disequilibrium over physical and genetic distance found in the accessions was comparable to those previously reported for several outcrossing species. A mixed linear model detected three significant SNPs in three candidate genes involved in plant adaptation to drought stress, individually explaining 12.7-15.7% of the phenotypic variance. The SNPs identified in the study will help to extend understanding, and to support genetic improvement of rubber trees grown in drought-affected regions.

The SnRK2 kinases modulate miRNA accumulation in Arabidopsis

MicroRNAs (miRNAs) regulate gene expression and play critical roles in growth and development as well as stress responses in eukaryotes. miRNA biogenesis in plants requires a processing complex that consists of the core components DICER-LIKE 1 (DCL1), SERRATE (SE) and HYPONASTIC LEAVES (HYL1). Here we show that inactivation of functionally redundant members of the SnRK2 kinases, which are the core components of abscisic acid (ABA) and osmotic stress signaling pathways, leads to reduction in miRNA accumulation under stress conditions. Further analysis revealed that the steady state level of HYL1 protein in plants under osmotic stress is dependent on the SnRK2 kinases. Additionally, our results suggest that the SnRK2 kinases physically associate with the miRNA processing components SE and HYL1 and can phosphorylate these proteins in vitro. These findings reveal an important role for the SnRK2 kinases in the regulation of miRNA accumulation and establish a mechanism by which ABA and osmotic stress signaling is linked to miRNA biogenesis.

Identified of a novel cis-element regulating the alternative splicing of LcDREB2

Alternative splicing (AS) is an important gene regulation mechanism in plants. Despite the widespread use of AS in plant gene expression regulation, the identification of the cis-elements involved in the AS mechanism is rarely reported in plants. To explore the regulation mechanism of the AS of LcDREB2, a DREB2 ortholog from Sheepgrass (Leymus chinensis), the genomic sequences of LcDREB2 and its homologs in Poaceae were aligned, and six mutations were introduced in the conserved sequence of LcDREB2. By analyzing the distinct transcript patterns of the LcDREB2 mutants in transgenic Oryza sativa, a novel cis-element that affected the AS of LcDREB2 was identified as Exonic Splicing Enhancer 1 (ESE1). In addition, five serine-arginine rich (SR) proteins were confirmed to interact with ESE1 by electrophoretic mobility shift assay (EMSA). To further explore the expression regulation mechanism of the DREB subfamily, phylogenetic analysis of DREB2 paralogous genes was performed. The results strongly supported the hypothesis that AS is conserved in Poaceae plants and that it is an evolutionary strategy for the regulation of the functional expression of genes. The findings and methods of our study will promote a substantial step forward in understanding of the plant AS regulation mechanism.

Active 5' splice sites regulate the biogenesis efficiency of Arabidopsis microRNAs derived from intron-containing genes

Arabidopsis, miR402 that is encoded within the first intron of a protein-coding gene At1g77230, is induced by heat stress. Its upregulation correlates with splicing inhibition and intronic proximal polyA site selection. It suggests that miR402 is not processed from an intron, but rather from a shorter transcript after selection of the proximal polyA site within this intron. Recently, introns and active 5΄ splice sites (5΄ss’) have been shown to stimulate the accumulation of miRNAs encoded within the first exons of intron-containing MIR genes. In contrast, we have observed the opposite effect of splicing inhibition on intronic miR402 production. Transient expression experiments performed in tobacco leaves revealed a significant accumulation of the intronic mature miR402 when the 5΄ss of the miR402-hosting intron was inactivated. In contrast, when the miR402 stem-loop structure was moved into the first exon, mutation of the first-intron 5΄ss resulted in a decrease in the miRNA level. Thus, the 5΄ss controls the efficiency of miRNA biogenesis. We also show that the SERRATE protein (a key component of the plant microprocessor) colocalizes and interacts with several U1 snRNP auxiliary proteins. We postulate that SERRATE-spliceosome connections have a direct effect on miRNA maturation.

Alternative Splicing in Plant Genes: A Means of Regulating the Environmental Fitness of Plants

Gene expression can be regulated through transcriptional and post-transcriptional mechanisms. Transcription in eukaryotes produces pre-mRNA molecules, which are processed and spliced post-transcriptionally to create translatable mRNAs. More than one mRNA may be produced from a single pre-mRNA by alternative splicing (AS); thus, AS serves to diversify an organism’s transcriptome and proteome. Previous studies of gene expression in plants have focused on the role of transcriptional regulation in response to environmental changes. However, recent data suggest that post-transcriptional regulation, especially AS, is necessary for plants to adapt to a changing environment. In this review, we summarize recent advances in our understanding of AS during plant development in response to environmental changes. We suggest that alternative gene splicing is a novel means of regulating the environmental fitness of plants.

Regulation of FT splicing by an endogenous cue in temperate grasses

Appropriate flowering timing is crucial for plant reproductive success. The florigen, FLOWERING LOCUS T (FT), interacts with 14-3-3 proteins and the bZIP transcription factor FD, functioning at core nodes in multiple flowering pathways. There are two FT homologues, FT1 and FT2, in Brachypodium distachyon. Here we show that FT2 undergoes age-dependent alternative splicing (AS), resulting in two splice variants (FT2α and FT2β). The FT2β-encoded protein cannot interact with FD or 14-3-3s but is able to form heterodimers with FT2α and FT1, thereby interfering with the florigen-mediated assembly of the flowering initiation complex. Notably, transgenic plants overproducing FT2β exhibit delayed flowering, while transgenic plants in which FT2β is silenced by an artificial microRNA display accelerated flowering, demonstrating a dominant-negative role of FT2β in flowering induction. Furthermore, we show that the AS splicing of FT2 is conserved in important cereal crops, such as barley and wheat. Collectively, these findings reveal a novel posttranscriptional mode of FT regulation in temperate grasses.

Transcriptome Analysis of Flowering Time Genes under Drought Stress in Maize Leaves

Flowering time is an important factor determining yield and seed quality in maize. A change in flowering time is a strategy used to survive abiotic stresses. Among abiotic stresses, drought can increase anthesis-silking intervals (ASI), resulting in negative effects on maize yield. We have analyzed the correlation between flowering time and drought stress using RNA-seq and bioinformatics tools. Our results identified a total of 619 genes and 126 transcripts whose expression was altered by drought stress in the maize B73 leaves under short-day condition. Among drought responsive genes, we also identified 20 genes involved in flowering times. Gene Ontology (GO) enrichment analysis was used to predict the functions of the drought-responsive genes and transcripts. GO categories related to flowering time included reproduction, flower development, pollen-pistil interaction, and post-embryonic development. Transcript levels of several genes that have previously been shown to affect flowering time, such as PRR37, transcription factor HY5, and CONSTANS, were significantly altered by drought conditions. Furthermore, we also identified several drought-responsive transcripts containing C₂H₂ zinc finger, CCCH, and NAC domains, which are frequently involved in transcriptional regulation and may thus have potential to alter gene expression programs to change maize flowering time. Overall, our results provide a genome-wide analysis of differentially expressed genes (DEGs), novel transcripts, and isoform variants expressed during the reproductive stage of maize plants subjected to drought stress and short-day condition. Further characterization of the drought-responsive transcripts identified in this study has the potential to advance our understanding of the mechanisms that regulate flowering time under drought stress.

MicroRNAs Are Intensively Regulated during Induction of Somatic Embryogenesis in Arabidopsis

Several genes encoding transcription factors (TFs) were indicated to have a key role in the induction of somatic embryogenesis (SE), which is triggered in the somatic cells of plants. In order to further explore the genetic regulatory network that is involved in the embryogenic transition induced in plant somatic cells, micro-RNA (miRNAs) molecules, the products of MIRNA (MIR) genes and the common regulators of TF transcripts, were analyzed in an embryogenic culture of Arabidopsis thaliana. In total, the expression of 190 genes of the 114 MIRNA families was monitored during SE induction and the levels of the primary (pri-miRNAs) transcripts vs. the mature miRNAs were investigated. The results revealed that the majority (98%) of the MIR genes were active and that most of them (64%) were differentially expressed during SE. A distinct attribute of the MIR expression in SE was the strong repression of MIR transcripts at the early stage of SE followed by their significant up-regulation in the advanced stage of SE. Comparison of the mature miRNAs vs. pri-miRNAs suggested that the extensive post-transcriptional regulation of miRNA is associated with SE induction. Candidate miRNA molecules of the assumed function in the embryogenic response were identified among the mature miRNAs that had a differential expression in SE, including miR156, miR157, miR159, miR160, miR164, miR166, miR169, miR319, miR390, miR393, miR396, and miR398. Consistent with the central role of phytohormones and stress factors in SE induction, the functions of the candidate miRNAs were annotated to phytohormone and stress responses. To confirm the functions of the candidate miRNAs in SE, the expression patterns of the mature miRNAs and their presumed targets were compared and regulatory relation during SE was indicated for most of the analyzed miRNA-target pairs. The results of the study contribute to the refinement of the miRNA-controlled regulatory pathways that operate during embryogenic induction in plants and provide a valuable platform for the identification of the genes that are targeted by the candidate miRNAs in SE induction.

New technologies accelerate the exploration of non-coding RNAs in horticultural plants

Non-coding RNAs (ncRNAs), that is, RNAs not translated into proteins, are crucial regulators of a variety of biological processes in plants. While protein-encoding genes have been relatively well-annotated in sequenced genomes, accounting for a small portion of the genome space in plants, the universe of plant ncRNAs is rapidly expanding. Recent advances in experimental and computational technologies have generated a great momentum for discovery and functional characterization of ncRNAs. Here we summarize the classification and known biological functions of plant ncRNAs, review the application of next-generation sequencing (NGS) technology and ribosome profiling technology to ncRNA discovery in horticultural plants and discuss the application of new technologies, especially the new genome-editing tool clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9) systems, to functional characterization of plant ncRNAs.

miRNA Biogenesis: A Dynamic Pathway

MicroRNAs (miRNAs) modulate plant homeostasis through the inactivation of specific mRNAs, especially those encoding transcription factors. A delicate spatial/temporal balance between a miRNA and its targets is central to achieving the appropriate biological outcomes. In this review we discuss our growing understanding of the dynamic regulation of miRNA biogenesis. We put special emphasis on crosstalk between miRNA biogenesis and other cellular processes such as transcription and splicing. We also discuss how the pathway is regulated in specific tissues to achieve harmonious plant development through a subtle balance between gene expression and silencing.

Arabidopsis YL1/BPG2 Is Involved in Seedling Shoot Response to Salt Stress through ABI4

The chloroplast-localized proteins play roles in plant salt stress response, but their mechanisms remain largely unknown. In this study, we screened a yellow leaf mutant, yl1-1, whose shoots exhibited hypersensitivity to salt stress. We mapped YL1 to AT3G57180, which encodes a YqeH-type GTPase. YL1, as a chloroplast stroma-localized protein, could be markedly reduced by high salinity. Upon exposure to high salinity, seedling shoots of yl1-1 and yl1-2 accumulated significantly higher levels of Na⁺ than wild type. Expression analysis of factors involved in plant salt stress response showed that the expression of ABI4 was increased and HKT1 was evidently suppressed in mutant shoots compared with the wild type under normal growth conditions. Moreover, salinity effects on ABI4 and HKT1 were clearly weakened in the mutant shoots, suggesting that the loss of YL1 function impairs ABI4 and HKT1 expression. Notably, the shoots of yl1-2 abi4 double mutant exhibited stronger resistance to salt stress and accumulated less Na⁺ levels after salt treatment compared with the yl1-2 single mutant, suggesting the salt-sensitive phenotype of yl1-2 seedlings could be rescued via loss of ABI4 function. These results reveal that YL1 is involved in the salt stress response of seedling shoots through ABI4.

Small RNAs: essential regulators of gene expression and defenses against environmental stresses in plants

Eukaryotic genomes produce thousands of diverse small RNAs (smRNAs), which play vital roles in regulating gene expression in all conditions, including in survival of biotic and abiotic environmental stresses. SmRNA pathways intersect with most of the pathways regulating different steps in the life of a messenger RNA (mRNA), starting from transcription and ending at mRNA decay. SmRNAs function in both nuclear and cytoplasmic compartments; the regulation of mRNA stability and translation in the cytoplasm and the epigenetic regulation of gene expression in the nucleus are the main and best-known modes of smRNA action. However, recent evidence from animal systems indicates that smRNAs and RNA interference (RNAi) also participate in the regulation of alternative pre-mRNA splicing, one of the most crucial steps in the fast, efficient global reprogramming of gene expression required for survival under stress. Emerging evidence from bioinformatics studies indicates that a specific class of plant smRNAs, induced by various abiotic stresses, the sutr-siRNAs, has the potential to target regulatory regions within introns and thus may act in the regulation of splicing in response to stresses. This review summarizes the major types of plant smRNAs in the context of their mechanisms of action and also provides examples of their involvement in regulation of gene expression in response to environmental cues and developmental stresses. In addition, we describe current advances in our understanding of how smRNAs function in the regulation of pre-mRNA splicing. WIREs RNA 2016, 7:356-381. doi: 10.1002/wrna.1340 For further resources related to this article, please visit the WIREs website.

Identification and Characterization of Erysiphe necator-Responsive MicroRNAs in Chinese Wild Vitis pseudoreticulata by High-Throughput Sequencing

Grapevine powdery mildew is one of the most damaging fungal diseases. Therefore, a precise understanding of the grapevine disease resistance system becomes a subject of significant importance. Plant microRNAs(miRNAs) have been implicated to play regulatory roles in plant biotic stress responses. In this study, high-throughput sequencing and miRDeep-P were employed to identify miRNAs in Chinese wild Vitis pseudoreticulata leaves following inoculation with Erysiphe necator. Altogether, 126 previously identified microRNAs and 124 novel candidates of miRNA genes were detected. Among them, 43 conserved miRNAs belong to 20 families and 23 non-conserved but previously-known miRNAs belong to 15 families. Following E. necator inoculation, 119 miRNAs were down-regulated and 131 were up-regulated. Furthermore, the expression changes occurring in 32 miRNAs were significant. The expression patterns of some miRNAs were validated by semi-quantitative RT-PCR and qRT-PCR. A total of 485 target genes were predicted and categorized by Gene Ontology (GO). In addition, 14 vvi-miRNAs were screened with 36 targets which may be involved in powdery mildew resistance in grape. Highly accumulated vvi-NewmiR2118 was detected from accession "Baihe-35-1," whose targets were mostly NBS-LRR resistance genes. It was down-regulated rapidly and strongly in "Baihe-35-1" leaves after inoculated with E. necator, indicating its involvement in grape powdery mildew resistance. Finally, the study verified interaction between vvi-NewmiR2118 and RPP13 by histochemical staining and GUS fluorescence quantitative assay.

Regulation of Non-coding RNAs in Heat Stress Responses of Plants

Heat stress is an important factor limiting plant growth, development, and productivity; thus, plants have evolved special adaptive mechanisms to cope with high-temperature stress. Non-coding RNAs (ncRNAs) are a class of regulatory RNAs that play an important role in many biological processes. Recently developed advanced technologies, such as genome-wide transcriptomic analysis, have revealed that abundant ncRNAs are expressed under heat stress. Although this area of research is still in its infancy, an increasing number of several classes of regulatory ncRNA (i.e., miRNA, siRNA, and lncRNA) related to heat stress responses have been reported. In this mini-review, we discuss our current understanding of the role of ncRNAs in heat stress responses in plants, especially miRNAs, siRNAs, and their targets. For example, the miR398-CSD/CCS-HSF, miR396-WRKY6, miR159-GAMYB, and TAS1-HTT-HSF pathways regulate plant heat tolerance. We highlight the hormone/development-related miRNAs involved in heat stress, and discuss the regulatory networks of miRNA-targets. We also note that DNA methylation and alternative splicing could affect miRNA expression under heat stress, and some lncRNAs could respond to heat stress. Finally, we briefly discuss future prospects concerning the ncRNA-related mechanisms of heat stress responses in plants.

Salt-induced transcription factor MYB74 is regulated by the RNA-directed DNA methylation pathway in Arabidopsis

Salt stress is one of the major abiotic stresses in agriculture worldwide that causes crop failure by interfering with the profile of gene expression and cell metabolism. Transcription factors and RNA-directed DNA methylation (RdDM) play an important role in the regulation of gene activation under abiotic stress in plants. This work characterized AtMYB74, a member of the R2R3-MYB gene family, which is transcriptionally regulated mainly by RdDM as a response in salt stress in Arabidopsis. Bisulphite sequencing indicated that 24-nt siRNAs target a region approximately 500bp upstream of the transcription initiation site of AtMYB74, which is heavily methylated. Levels of DNA methylation in this region were significantly reduced in wild type plants under salt stress, whereas no changes were found in RdDM mutants. Northern blot and quantitative real-time reverse transcription PCR analysis showed that the accumulation of 24-nt siRNAs was decreased in WT plants under salt stress. Further promoter deletion analysis revealed that the siRNA target region is essential for maintaining AtMYB74 expression patterns. In addition, transgenic plants overexpressing AtMYB74 displayed hypersensitivity to NaCl during seed germination. These results suggest that changes in the levels of the five 24-nt siRNAs regulate the AtMYB74 transcription factor via RdDM in response to salt stress.

ΔNp63 intronic miR-944 is implicated in the ΔNp63-mediated induction of epidermal differentiation

ΔNp63 is required for both the proliferation and differentiation of keratinocytes, but its role in the differentiation of these cells is poorly understood. The corresponding gene, TP63, harbors the MIR944 sequence within its intron. However, the mechanism of biogenesis and the function of miR-944 are unknown. We found that miR-944 is highly expressed in keratinocytes, in a manner that is concordant with that of ΔNp63 mRNA, but the regulation of miR-944 expression under various conditions did not correspond with that of ΔNp63. Bioinformatics analysis and functional studies demonstrated that MIR944 has its own promoter. We demonstrate here that MIR944 is a target of ΔNp63. Promoter analysis revealed that the activity of the MIR944 promoter was markedly enhanced by the binding of ΔNp63, which was maintained by the supportive action of AP-2 during keratinocyte differentiation. Our results indicated that miR-944 biogenesis is dependent on ΔNp63 protein, even though it is generated from ΔNp63 mRNA-independent transcripts. We also demonstrated that miR-944 induces keratin 1 and keratin 10 expression by inhibiting ERK signaling and upregulating p53 expression. Our findings suggested that miR-944, as an intronic miRNA and a direct target of ΔNp63, contributes to the function of ΔNp63 in the induction of epidermal differentiation.

Nonadditive gene expression in polyploids

Allopolyploidy involves hybridization and duplication of divergent parental genomes and provides new avenues for gene expression. The expression levels of duplicated genes in polyploids can show deviation from parental additivity (the arithmetic average of the parental expression levels). Nonadditive expression has been widely observed in diverse polyploids and comprises at least three possible scenarios: (a) The total gene expression level in a polyploid is similar to that of one of its parents (expression-level dominance); (b) total gene expression is lower or higher than in both parents (transgressive expression); and (c) the relative contribution of the parental copies (homeologs) to the total gene expression is unequal (homeolog expression bias). Several factors may result in expression nonadditivity in polyploids, including maternal-paternal influence, gene dosage balance, cis- and/or trans-regulatory networks, and epigenetic regulation. As our understanding of nonadditive gene expression in polyploids remains limited, a new generation of investigators should explore additional phenomena (i.e., alternative splicing) and use other high-throughput "omics" technologies to measure the impact of nonadditive expression on phenotype, proteome, and metabolome.

New insights into pri-miRNA processing and accumulation in plants

MicroRNAs (miRNAs) regulate many biological processes such as development, metabolism, and others. They are processed from their primary transcripts called primary miRNA transcripts (pri-miRNAs) by the processor complex containing the RNAse III enzyme, DICER-LIKE1 (DCL1), in plants. Consequently, miRNA biogenesis is controlled through altering pri-miRNA accumulation and processing, which is crucial for plant development and adaptation to environmental changes. Plant pri-miRNAs are transcribed by DNA-dependent RNA polymerase II (Pol II) and their levels are determined through transcription and degradation, whereas pri-miRNA processing is affected by its structure, splicing, alternative splicing, loading to the processor and the processor activity, which involve in many accessory proteins. Here, we summarize recent progresses related to pri-miRNA transcription, stability, and processing in plants.

Alternative splicing in plants: directing traffic at the crossroads of adaptation and environmental stress

In recent years, high-throughput sequencing-based analysis of plant transcriptomes has suggested that up to ∼60% of plant gene loci encode alternatively spliced mature transcripts. These studies have also revealed that alternative splicing in plants can be regulated by cell type, developmental stage, the environment, and the circadian clock. Alternative splicing is coupled to RNA surveillance and processing mechanisms, including nonsense mediated decay. Recently, non-protein-coding transcripts have also been shown to undergo alternative splicing. These discoveries collectively describe a robust system of post-transcriptional regulatory feedback loops which influence RNA abundance. In this review, we summarize recent studies describing the specific roles alternative splicing and RNA surveillance play in plant adaptation to environmental stresses and the regulation of the circadian clock.

Bioinformatics analysis of small RNAs in pima (Gossypium barbadense L.)

Small RNAs (sRNAs) are ~20 to 24 nucleotide single-stranded RNAs that play crucial roles in regulation of gene expression. In plants, sRNAs are classified into microRNAs (miRNAs), repeat-associated siRNAs (ra-siRNAs), phased siRNAs (pha-siRNAs), cis and trans natural antisense transcript siRNAs (cis- and trans-nat siRNAs). Pima (Gossypium barbadense L.) is one of the most economically important fiber crops, producing the best and longest spinnable fiber. Although some miRNAs are profiled in Pima, little is known about siRNAs, the largest subclass of plant sRNAs. In order to profile these gene regulators in Pima, a comprehensive analysis of sRNAs was conducted by mining publicly available sRNA data, leading to identification of 678 miRNAs, 3,559,126 ra-siRNAs, 627 pha-siRNAs, 136,600 cis-nat siRNAs and 79,994 trans-nat siRNAs. The 678 miRNAs, belonging to 98 conserved and 402 lineage-specific families, were produced from 2,138 precursors, of which 297 arose from introns, exons, or intron/UTR-exon junctions of protein-coding genes. Ra-siRNAs were produced from various repeat loci, while most (97%) were yielded from retrotransposons, especially LTRs (long terminal repeats). The genes encoding auxin-signaling-related proteins, NBS-LRRs and transcription factors were major sources of pha-siRNAs, while two conserved TAS3 homologs were found as well. Most cis-NATs in Pima overlapped in enclosed and convergent orientations, while a few hybridized in divergent and coincided orientations. Most cis- and trans-nat siRNAs were produced from overlapping regions. Additionally, characteristics of length and the 5’-first nucleotide of each sRNA class were analyzed as well. Results in this study created a valuable molecular resource that would facilitate studies on mechanism of controlling gene expression.

Post-transcriptional and post-translational regulations of drought and heat response in plants: a spider's web of mechanisms

Drought and heat tolerance are complex quantitative traits. Moreover, the adaptive significance of some stress-related traits is more related to plant survival than to agronomic performance. A web of regulatory mechanisms fine-tunes the expression of stress-related traits and integrates both environmental and developmental signals. Both post-transcriptional and post-translational modifications contribute substantially to this network with a pivotal regulatory function of the transcriptional changes related to cellular and plant stress response. Alternative splicing and RNA-mediated silencing control the amount of specific transcripts, while ubiquitin and SUMO modify activity, sub-cellular localization and half-life of proteins. Interactions across these modification mechanisms ensure temporally and spatially appropriate patterns of downstream-gene expression. For key molecular components of these regulatory mechanisms, natural genetic diversity exists among genotypes with different behavior in terms of stress tolerance, with effects upon the expression of adaptive morphological and/or physiological target traits.

Stress-induced endogenous siRNAs targeting regulatory intron sequences in Brachypodium

Exposure to abiotic stresses triggers global changes in the expression of thousands of eukaryotic genes at the transcriptional and post-transcriptional levels. Small RNA (smRNA) pathways and splicing both function as crucial mechanisms regulating stress-responsive gene expression. However, examples of smRNAs regulating gene expression remain largely limited to effects on mRNA stability, translation, and epigenetic regulation. Also, our understanding of the networks controlling plant gene expression in response to environmental changes, and examples of these regulatory pathways intersecting, remains limited. Here, to investigate the role of smRNAs in stress responses we examined smRNA transcriptomes of Brachypodium distachyon plants subjected to various abiotic stresses. We found that exposure to different abiotic stresses specifically induced a group of novel, endogenous small interfering RNAs (stress-induced, UTR-derived siRNAs, or sutr-siRNAs) that originate from the 3' UTRs of a subset of coding genes. Our bioinformatics analyses predicted that sutr-siRNAs have potential regulatory functions and that over 90% of sutr-siRNAs target intronic regions of many mRNAs in trans. Importantly, a subgroup of these sutr-siRNAs target the important intron regulatory regions, such as branch point sequences, that could affect splicing. Our study indicates that in Brachypodium, sutr-siRNAs may affect splicing by masking or changing accessibility of specific cis-elements through base-pairing interactions to mediate gene expression in response to stresses. We hypothesize that this mode of regulation of gene expression may also serve as a general mechanism for regulation of gene expression in plants and potentially in other eukaryotes.

Dissection of miRNA pathways using arabidopsis mesophyll protoplasts

MicroRNAs (miRNAs) control gene expression mostly post-transcriptionally by guiding transcript cleavage and/or translational repression of complementary mRNA targets, thereby regulating developmental processes and stress responses. Despite the remarkable expansion of the field, the mechanisms underlying miRNA activity are not fully understood. In this article, we describe a transient expression system in Arabidopsis mesophyll protoplasts, which is highly amenable for the dissection of miRNA pathways. We show that by transiently overexpressing primary miRNAs and target mimics, we can manipulate miRNA levels and consequently impact on their targets. Furthermore, we developed a set of luciferase-based sensors for quantifying miRNA activity that respond specifically to both endogenous and overexpressed miRNAs and target mimics. We demonstrate that these miRNA sensors can be used to test the impact of putative components of the miRNA pathway on miRNA activity, as well as the impact of specific mutations, by either overexpression or the use of protoplasts from the corresponding mutants. We further show that our miRNA sensors can be used for investigating the effect of chemicals on miRNA activity. Our cell-based transient expression system is fast and easy to set up, and generates quantitative results, being a powerful tool for assaying miRNA activity in vivo.

Arabidopsis microRNA expression regulation in a wide range of abiotic stress responses

Arabidopsis microRNA expression regulation was studied in a wide array of abiotic stresses such as drought, heat, salinity, copper excess/deficiency, cadmium excess, and sulfur deficiency. A home-built RT-qPCR mirEX platform for the amplification of 289 Arabidopsis microRNA transcripts was used to study their response to abiotic stresses. Small RNA sequencing, Northern hybridization, and TaqMan® microRNA assays were performed to study the abundance of mature microRNAs. A broad response on the level of primary miRNAs (pri-miRNAs) was observed. However, stress response at the level of mature microRNAs was rather confined. The data presented show that in most instances, the level of a particular mature miRNA could not be predicted based on the level of its pri-miRNA. This points to an essential role of posttranscriptional regulation of microRNA expression. New Arabidopsis microRNAs responsive to abiotic stresses were discovered. Four microRNAs: miR319a/b, miR319b.2, and miR400 have been found to be responsive to several abiotic stresses and thus can be regarded as general stress-responsive microRNA species.

Genetic and epigenetic control of plant heat responses

Plants have evolved sophisticated genetic and epigenetic regulatory systems to respond quickly to unfavorable environmental conditions such as heat, cold, drought, and pathogen infections. In particular, heat greatly affects plant growth and development, immunity and circadian rhythm, and poses a serious threat to the global food supply. According to temperatures exposing, heat can be usually classified as warm ambient temperature (about 22-27°C), high temperature (27-30°C) and extremely high temperature (37-42°C, also known as heat stress) for the model plant Arabidopsis thaliana. The genetic mechanisms of plant responses to heat have been well studied, mainly focusing on elevated ambient temperature-mediated morphological acclimation and acceleration of flowering, modulation of circadian clock and plant immunity by high temperatures, and thermotolerance to heat stress. Recently, great progress has been achieved on epigenetic regulation of heat responses, including DNA methylation, histone modifications, histone variants, ATP-dependent chromatin remodeling, histone chaperones, small RNAs, long non-coding RNAs and other undefined epigenetic mechanisms. These epigenetic modifications regulate the expression of heat-responsive genes and function to prevent heat-related damages. This review focuses on recent progresses regarding the genetic and epigenetic control of heat responses in plants, and pays more attention to the role of the major epigenetic mechanisms in plant heat responses. Further research perspectives are also discussed.

microRNA biogenesis, degradation and activity in plants

microRNAs (miRNAs) are important regulators of gene expression. After excised from primary miRNA transcript by dicer-like1 (DCL1, an RNAse III enzyme), miRNAs bind and guide their effector protein named argonaute 1 (AGO1) to silence the expression of target RNAs containing their complementary sequences in plants. miRNA levels and activities are tightly controlled to ensure their functions in various biological processes such as development, metabolism and responses to abiotic and biotic stresses. Studies have identified many factors that involve in miRNA accumulation and activities. Characterization of these factors in turn greatly improves our understanding of the processes related to miRNAs. Here, we review recent progress of mechanisms underlying miRNA expression and functions in plants.

Comparisons of computational methods for differential alternative splicing detection using RNA-seq in plant systems

BACKGROUND

Alternative Splicing (AS) as a post-transcription regulation mechanism is an important application of RNA-seq studies in eukaryotes. A number of software and computational methods have been developed for detecting AS. Most of the methods, however, are designed and tested on animal data, such as human and mouse. Plants genes differ from those of animals in many ways, e.g., the average intron size and preferred AS types. These differences may require different computational approaches and raise questions about their effectiveness on plant data. The goal of this paper is to benchmark existing computational differential splicing (or transcription) detection methods so that biologists can choose the most suitable tools to accomplish their goals.

RESULTS

This study compares the eight popular public available software packages for differential splicing analysis using both simulated and real Arabidopsis thaliana RNA-seq data. All software are freely available. The study examines the effect of varying AS ratio, read depth, dispersion pattern, AS types, sample sizes and the influence of annotation. Using a real data, the study looks at the consistences between the packages and verifies a subset of the detected AS events using PCR studies.

CONCLUSIONS

No single method performs the best in all situations. The accuracy of annotation has a major impact on which method should be chosen for AS analysis. DEXSeq performs well in the simulated data when the AS signal is relative strong and annotation is accurate. Cufflinks achieve a better tradeoff between precision and recall and turns out to be the best one when incomplete annotation is provided. Some methods perform inconsistently for different AS types. Complex AS events that combine several simple AS events impose problems for most methods, especially for MATS. MATS stands out in the analysis of real RNA-seq data when all the AS events being evaluated are simple AS events.

Mutational bias is the driving force for shaping the synonymous codon usage pattern of alternatively spliced genes in rice (Oryza sativa L.)

Alternative splicing plays important roles in diverse aspects of plant development, metabolism, and stress responses. However, the regulatory mechanisms of alternative splicing of genes still remain incompletely elucidated, especially in plants. In this study, the synonymous codon usage pattern of alternatively spliced (AS) genes in rice was firstly explored using the combination of correspondence analysis (CA), internal CA, correlation and ANOVA analyses. The results show that alternatively and non-alternatively spliced (non-AS) genes have similar tendency for overall codon usage, but exhibit significant difference in 58 out of 64 codons. AS and non-AS genes are both under strong purifying selection, but the former ones have significant lower mutation rate and are prone to be enriched towards the chromosomal ends. In the group of AS genes, the variability in synonymous codon usage between genes is mainly due to the variations in GC content, CDS length, as well as gene functions. Mutational bias that accounts for 25.85 % of the total codon usage variability plays a major role in shaping the codon usage pattern of AS genes. In contrast, no obvious evidence is found for the contributions of translational selection, AS types, the conservation of AS events, and numbers of AS variants to the codon usage divergence between AS genes. These findings may be useful for further understanding the mechanisms of origination, differentiation and regulation of alternatively spliced genes in plants.

MicroRNA400-guided cleavage of Pentatricopeptide repeat protein mRNAs Renders Arabidopsis thaliana more susceptible to pathogenic bacteria and fungi

Although a large number of microRNAs (miRNAs) have been identified in different plant species, the functional roles and targets of the majority of miRNAs have not yet been determined. Here, Arabidopsis thaliana miRNA400 (miR400) was investigated for its functional role in the defense response to diverse pathogens. Transgenic Arabidopsis plants that overexpress MIR400 (35S::MIR400) displayed much more severe disease symptoms than the wild-type plants when infected with the bacterium Pseudomonas syringae pv. tomato DC3000 or the fungus Botrytis cinerea. MiR400 guided the cleavage of two genes (At1g06580 and At1g62720) encoding pentatricopeptide repeat (PPR) proteins. To confirm further that the miR400-mediated defense response was due to the cleavage of PPR mRNAs, loss-of-function mutant and artificial miRNA-mediated knockdown mutants of PPR were generated, and their disease responses were analyzed upon pathogen challenge. Similar to the 35S::MIR400 plants, the ppr mutants displayed much more severe disease symptoms than the wild-type plants when challenged with the pathogens, indicating that miR400 affects the defense response by cleaving PPR mRNAs. Expression of miR400 was down-regulated, whereas the PPR1 and PPR2 transcripts increased upon pathogen challenge. Collectively, the present study reveals that miR400-mediated dysfunction of PPR proteins renders Arabidopsis more susceptible to pathogenic bacteria and fungi, which emphasizes the importance of PPR proteins in plant defense against diverse pathogens.