Alternative mRNA processing increases the complexity of microRNA-based gene regulation in Arabidopsis

microRNA858 represses the transcription factor gene SbMYB47 and regulates flavonoid biosynthesis in Scutellaria baicalensis

MicroRNAs (miRNAs) are noncoding endogenous single-stranded RNAs that regulate target gene expression by reducing their transcription and translation. Several miRNAs in plants function in secondary metabolism. The dried root of Scutellaria baicalensis Georgi is a traditional Chinese medicine that contains flavonoids (baicalin, wogonoside, and baicalein) as its main active ingredients. Although the S. baicalensis genome sequence has been published, information regarding its miRNAs is lacking. In this study, 12 small RNA libraries of different S. baicalensis tissues were compiled, including roots, stems, leaves, and flowers. A total of 129 miRNAs were identified, including 99 miRNAs from 27 miRNA families and 30 predicted miRNAs. Furthermore, 46 reliable target genes of 15 miRNA families were revealed using psRNATarget and confirmed by degradome sequencing. It was speculated that the microRNA858 (miR858)-SbMYB47 module might be involved in flavonoid biosynthesis. Transient assays in Nicotiana benthamiana leaves indicated that miR858 targets SbMYB47 and suppresses its expression. Artificial miRNA-mediated knockdown of miR858 and overexpression of SbMYB47 significantly increased the flavonoid content in S. baicalensis hairy roots, while SbMYB47 knockdown inhibited flavonoid accumulation. Yeast one-hybrid and dual-luciferase assays indicated that SbMYB47 directly binds to and activates the S. baicalensis phenylalanine ammonia-lyase 3 (SbPAL-3) and flavone synthase II (SbFNSⅡ-2) promoters. Our findings reveal the link between the miR858-SbMYB47 module and flavonoid biosynthesis, providing a potential strategy for the production of flavonoids with important pharmacological activities through metabolic engineering.

Phytopathogens Reprogram Host Alternative mRNA Splicing

Alternative splicing (AS) is an evolutionarily conserved cellular process in eukaryotes in which multiple messenger RNA (mRNA) transcripts are produced from a single gene. The concept that AS adds to transcriptome complexity and proteome diversity introduces a new perspective for understanding how phytopathogen-induced alterations in host AS cause diseases. Recently, it has been recognized that AS represents an integral component of the plant immune system during parasitic, commensalistic, and symbiotic interactions. Here, I provide an overview of recent progress detailing the reprogramming of plant AS by phytopathogens and the functional implications on disease phenotypes. Additionally, I discuss the vital function of AS of immune receptors in regulating plant immunity and how phytopathogens use effector proteins to target key components of the splicing machinery and exploit alternatively spliced variants of immune regulators to negate defense responses. Finally, the functional association between AS and nonsense-mediated mRNA decay in the context of plant-pathogen interface is recapitulated.

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

Introduction

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

Methods

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

Results and discussion

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

Conclusion

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

The Splicing Factor SR45 Negatively Regulates Anthocyanin Accumulation under High-Light Stress in Arabidopsis thaliana

High-intensity light (HL) greatly induces the accumulation of anthocyanin, a fundamental compound in photoprotection and antioxidation. Many mechanisms regulating anthocyanin biosynthesis are well-characterized across developmental and environmental conditions; however, post-transcriptional regulation of its biosynthesis remains unclear. RNA splicing is one mechanism of post-transcriptional control and reprogramming in response to different developmental cues and stress conditions. The Arabidopsis splicing modulator SR45 regulates a number of developmental and environmental stress responses. Here, we investigated the role of SR45 and its isoforms in HL-induced anthocyanin accumulation. We found that the SR45 promoter contains light-responsive cis-elements, and that light stress significantly increases SR45 expression. Furthermore, we found that mutant plants lacking SR45 function (sr45) accumulate significantly more anthocyanin under HL. SR45 is alternatively spliced to produce two proteins, SR45.1 and SR45.2, which differ by seven amino acids. Intriguingly, these isoforms exhibited distinct functions, with only SR45.1 reversing anthocyanin accumulation in the sr45 plants. We also identified possible SR45 target genes that are involved in anthocyanin synthesis. Consistent with the antioxidant role of anthocyanin, we found that sr45 mutants and SR45.2 overexpression lines accumulate anthocyanin and better tolerate paraquat which induces oxidative stress. Collectively, our results reveal that the Arabidopsis splicing regulator SR45 inhibits anthocyanin accumulation under HL, which may negatively affect oxidative stress tolerance. This study illuminates splicing-level regulation of anthocyanin production in response to light stress and offers a possible target for genetic modification to increase plant stress tolerance.

Lipid-based nucleic acid therapeutics with in vivo efficacy

Synthetic vectors for therapeutic nucleic acid delivery are currently competing significantly with their viral counter parts due to their reduced immunogenicity, large payload capacity, and ease of manufacture under GMP-compliant norms. The approval of Onpattro, a lipid-based siRNA therapeutic, and the proven clinical success of two lipid-based COVID-19 vaccines from Pfizer-BioNTech, and Moderna heralded the specific advantages of lipid-based systems among all other synthetic nucleic acid carriers. Lipid-based systems with diverse payloads-plasmid DNA (pDNA), antisense oligonucleotide (ASO), small interfering RNA (siRNA), microRNA (miRNA), small activating RNA (saRNA), and messenger RNA (mRNA)-are now becoming a mature technology, with growing impact in the clinic. Research over four decades identified the key factors determining the therapeutic success of these multi-component systems. Here, we discuss the main nucleic acid-based technologies, presenting their mechanism of action, delivery barriers facing them, the structural properties of the payload as well as the component lipids that regulate physicochemical properties, pharmacokinetics and biodistribution, efficacy, and toxicity of the resultant nanoparticles. We further detail on the formulation parameters, evolution of the manufacturing techniques that generate reproducible and scalable outputs, and key manufacturing aspects that enable control over physicochemical properties of the resultant particles. Preclinical applications of some of these formulations that were successfully translated from in vitro studies to animal models are subsequently discussed. Finally, clinical success and failure of these systems starting from 1993 to present are highlighted, in a holistic literature review focused on lipid-based nucleic acid delivery systems. This article is categorized under: Therapeutic Approaches and Drug Discovery > Emerging Technologies Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic Disease Toxicology and Regulatory Issues in Nanomedicine > Toxicology of Nanomaterials.

AGO1 and HSP90 buffer different genetic variants in Arabidopsis thaliana

Argonaute 1 (AGO1), the principal protein component of microRNA-mediated regulation, plays a key role in plant growth and development. AGO1 physically interacts with the chaperone HSP90, which buffers cryptic genetic variation in plants and animals. We sought to determine whether genetic perturbation of AGO1 in Arabidopsis thaliana would also reveal cryptic genetic variation, and if so, whether AGO1-dependent loci overlap with those dependent on HSP90. To address these questions, we introgressed a hypomorphic mutant allele of AGO1 into a set of mapping lines derived from the commonly used Arabidopsis strains Col-0 and Ler. Although we identified several cases in which AGO1 buffered genetic variation, none of the AGO1-dependent loci overlapped with those buffered by HSP90 for the traits assayed. We focused on 1 buffered locus where AGO1 perturbation uncoupled the traits days to flowering and rosette leaf number, which are otherwise closely correlated. Using a bulk segregant approach, we identified a nonfunctional Ler hua2 mutant allele as the causal AGO1-buffered polymorphism. Introduction of a nonfunctional hua2 allele into a Col-0 ago1 mutant background recapitulated the Ler-dependent ago1 phenotype, implying that coupling of these traits involves different molecular players in these closely related strains. Taken together, our findings demonstrate that even though AGO1 and HSP90 buffer genetic variation in the same traits, these robustness regulators interact epistatically with different genetic loci, suggesting that higher-order epistasis is uncommon. Plain Language Summary Argonaute 1 (AGO1), a key player in plant development, interacts with the chaperone HSP90, which buffers environmental and genetic variation. We found that AGO1 buffers environmental and genetic variation in the same traits; however, AGO1-dependent and HSP90-dependent loci do not overlap. Detailed analysis of a buffered locus found that a nonfunctional HUA2 allele decouples days to flowering and rosette leaf number in an AGO1-dependent manner, suggesting that the AGO1-dependent buffering acts at the network level.

Structural and Functional Analyses of Hub MicroRNAs in An Integrated Gene Regulatory Network of Arabidopsis

MicroRNAs (miRNAs) are trans-acting small regulatory RNAs that work coordinately with transcription factors (TFs) to shape the repertoire of cellular mRNAs available for translation. Despite our growing knowledge of individual plant miRNAs, their global roles in gene regulatory networks remain mostly unassessed. Based on interactions obtained from public databases and curated from the literature, we reconstructed an integrated miRNA network in Arabidopsis that includes 66 core TFs, 318 miRNAs, and 1712 downstream genes. We found that miRNAs occupy distinct niches and enrich miRNA-containing feed-forward loops (FFLs), particularly those with miRNAs as intermediate nodes. Further analyses revealed that miRNA-containing FFLs coordinate TFs located in different hierarchical layers and that intertwined miRNA-containing FFLs are associated with party and date miRNA hubs. Using the date hub MIR858A as an example, we performed detailed molecular and genetic analyses of three interconnected miRNA-containing FFLs. These analyses revealed individual functions of the selected miRNA-containing FFLs and elucidated how the date hub miRNA fulfills multiple regulatory roles. Collectively, our findings highlight the prevalence and importance of miRNA-containing FFLs, and provide new insights into the design principles and control logics of miRNA regulatory networks governing gene expression programs in plants.

Moderate Soil Drying-Induced Alternative Splicing Provides a Potential Novel Approach for the Regulation of Grain Filling in Rice Inferior Spikelets

Poor grain filling of inferior spikelets, especially in some large-panicle rice varieties, is becoming a major limitation in breaking the ceiling of rice production. In our previous studies, we proved that post-anthesis moderate soil drying (MD) was an effective way to promote starch synthesis and inferior grain filling. As one of the most important regulatory processes in response to environmental cues and at different developmental stages, the function of alternative splicing (AS) has not yet been revealed in regulating grain filling under MD conditions. In this study, AS events at the most active grain-filling stage were identified in inferior spikelets under well-watered control (CK) and MD treatments. Of 16,089 AS events, 1840 AS events involving 1392 genes occurred differentially between the CK and MD treatments, many of which function on spliceosome, ncRNA metabolic process, starch, and sucrose metabolism, and other functions. Some of the splicing factors and starch synthesis-related genes, such as SR protein, hnRNP protein, OsAGPL2, OsAPS2, OsSSIVa, OsSSIVb, OsGBSSII, and OsISA1 showed differential AS changes under MD treatment. The expression of miR439f and miR444b was reduced due to an AS event which occurred in the intron where miRNAs were located in the MD-treated inferior spikelets. On the contrary, OsAGPL2, an AGPase encoding gene, was alternatively spliced, resulting in different transcripts with or without the miR393b binding site, suggesting a potential mechanism for miRNA-mediated gene regulation on grain filling of inferior spikelets in response to MD treatment. This study provides some new insights into the function of AS on the MD-promoted grain filling of inferior spikelets, and potential application in agriculture to increase rice yields by genetic approaches.

44 Current Challenges in miRNomics

Mature microRNAs (miRNAs) are short RNA sequences about 18-24 nucleotide long, which provide the recognition key within RISC for the posttranscriptional regulation of target RNAs. Considering the canonical pathway, mature miRNAs are produced via a multistep process. Their transcription (pri-miRNAs) and first processing step via the microprocessor complex (pre-miRNAs) occur in the nucleus. Then they are exported into the cytosol, processed again by Dicer (dsRNA) and finally a single strand (mature miRNA) is incorporated into RISC (miRISC). The sequence of the incorporated miRNA provides the function of RNA target recognition via hybridization. Following binding of the target, the mRNA is either degraded or translation is inhibited, which ultimately leads to less protein production. Conversely, it has been shown that binding within the 5' UTR of the mRNA can lead to an increase in protein product. Regulation of homeostasis is very important for a cell; therefore, all steps in the miRNA-based regulation pathway, from transcription to the incorporation of the mature miRNA into RISC, are under tight control. While much research effort has been exerted in this area, the knowledgebase is not sufficient for accurately modelling miRNA regulation computationally. The computational prediction of miRNAs is, however, necessary because it is not feasible to investigate all possible pairs of a miRNA and its target, let alone miRNAs and their targets. We here point out open challenges important for computational modelling or for our general understanding of miRNA-based regulation and show how their investigation is beneficial. It is our hope that this collection of challenges will lead to their resolution in the near future.

Genome-Wide Analysis of Alternative Splicing and Non-Coding RNAs Reveal Complicated Transcriptional Regulation in Cannabis sativa L

It is of significance to mine the structural genes related to the biosynthetic pathway of fatty acid (FA) and cellulose as well as explore the regulatory mechanism of alternative splicing (AS), microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) in the biosynthesis of cannabinoids, FA and cellulose, which would enhance the knowledge of gene expression and regulation at post-transcriptional level in Cannabis sativa L. In this study, transcriptome, small RNA and degradome libraries of hemp 'Yunma No.1' were established, and comprehensive analysis was performed. As a result, a total of 154, 32 and 331 transcripts encoding key enzymes involved in the biosynthesis of cannabinoids, FA and cellulose were predicted, respectively, among which AS occurred in 368 transcripts. Moreover, 183 conserved miRNAs, 380 C. sativa-specific miRNAs and 7783 lncRNAs were predicted. Among them, 70 miRNAs and 17 lncRNAs potentially targeted 13 and 17 transcripts, respectively, encoding key enzymes or transporters involved in the biosynthesis of cannabinoids, cellulose or FA. Finally, the crosstalk between AS and miRNAs or lncRNAs involved in cannabinoids and cellulose was also predicted. In summary, all these results provided insights into the complicated network of gene expression and regulation in C. sativa.

Alternative splicing regulation appears to play a crucial role in grape berry development and is also potentially involved in adaptation responses to the environment

BACKGROUND

Alternative splicing (AS) produces transcript variants playing potential roles in proteome diversification and gene expression regulation. AS modulation is thus essential to respond to developmental and environmental stimuli. In grapevine, a better understanding of berry development is crucial for implementing breeding and viticultural strategies allowing adaptation to climate changes. Although profound changes in gene transcription have been shown to occur in the course of berry ripening, no detailed study on splicing modifications during this period has been published so far. We report here on the regulation of gene AS in developing berries of two grapevine (Vitis vinifera L.) varieties, Gewurztraminer (Gw) and Riesling (Ri), showing distinctive phenotypic characteristics. Using the software rMATS, the transcriptomes of berries at four developmental steps, from the green stage to mid-ripening, were analysed in pairwise comparisons between stages and varieties.

RESULTS

A total of 305 differential AS (DAS) events, affecting 258 genes, were identified. Interestingly, 22% of these AS events had not been reported before. Among the 80 genes that underwent the most significant variations during ripening, 22 showed a similar splicing profile in Gw and Ri, which suggests their involvement in berry development. Conversely, 23 genes were subjected to splicing regulation in only one variety. In addition, the ratios of alternative isoforms were different in Gw and Ri for 35 other genes, without any change during ripening. This last result indicates substantial AS differences between the two varieties. Remarkably, 8 AS events were specific to one variety, due to the lack of a splice site in the other variety. Furthermore, the transcription rates of the genes affected by stage-dependent splicing regulation were mostly unchanged, identifying AS modulation as an independent way of shaping the transcriptome.

CONCLUSIONS

The analysis of AS profiles in grapevine varieties with contrasting phenotypes revealed some similarity in the regulation of several genes with developmental functions, suggesting their involvement in berry ripening. Additionally, many splicing differences were discovered between the two varieties, that could be linked to phenotypic specificities and distinct adaptive capacities. Together, these findings open perspectives for a better understanding of berry development and for the selection of grapevine genotypes adapted to climate change.

MicroRNA775 regulates intrinsic leaf size and reduces cell wall pectin levels by targeting a galactosyltransferase gene in Arabidopsis

Plants possess unique primary cell walls made of complex polysaccharides that play critical roles in determining intrinsic cell and organ size. How genes responsible for synthesizing and modifying the polysaccharides in the cell wall are regulated by microRNAs (miRNAs) to control plant size remains largely unexplored. Here we identified 23 putative cell wall-related miRNAs, termed as CW-miRNAs, in Arabidopsis thaliana and characterized miR775 as an example. We showed that miR775 post-transcriptionally silences GALT9, which encodes an endomembrane-located galactosyltransferase belonging to the glycosyltransferase 31 family. Over-expression of miR775 and deletion of GALT9 led to significantly enlarged leaf-related organs, primarily due to increased cell size. Monosaccharide quantification, confocal Raman imaging, and immunolabeling combined with atomic force microscopy revealed that the MIR775A-GALT9 circuit modulates pectin levels and the elastic modulus of the cell wall. We also showed that MIR775A is directly repressed by the transcription factor ELONGATED HYPOCOTYL5 (HY5). Genetic analysis confirmed that HY5 is a negative regulator of leaf size that acts through the HY5-MIR775A-GALT9 repression cascade to control pectin levels. These findings demonstrate that miR775-regulated cell wall remodeling is an integral determinant of intrinsic leaf size in A. thaliana. Studying other CW-miRNAs would provide more insights into cell wall biology.

Molecular Tools for Adapting Viticulture to Climate Change

Adaptation of viticulture to climate change includes exploration of new geographical areas, new training systems, new management practices, or new varieties, both for rootstocks and scions. Molecular tools can be defined as molecular approaches used to study DNAs, RNAs, and proteins in all living organisms. We present here the current knowledge about molecular tools and their potential usefulness in three aspects of grapevine adaptation to the ongoing climate change. (i) Molecular tools for understanding grapevine response to environmental stresses. A fine description of the regulation of gene expression is a powerful tool to understand the physiological mechanisms set up by the grapevine to respond to abiotic stress such as high temperatures or drought. The current knowledge on gene expression is continuously evolving with increasing evidence of the role of alternative splicing, small RNAs, long non-coding RNAs, DNA methylation, or chromatin activity. (ii) Genetics and genomics of grapevine stress tolerance. The description of the grapevine genome is more and more precise. The genetic variations among genotypes are now revealed with new technologies with the sequencing of very long DNA molecules. High throughput technologies for DNA sequencing also allow now the genetic characterization at the same time of hundreds of genotypes for thousands of points in the genome, which provides unprecedented datasets for genotype-phenotype associations studies. We review the current knowledge on the genetic determinism of traits for the adaptation to climate change. We focus on quantitative trait loci and molecular markers available for developmental stages, tolerance to water stress/water use efficiency, sugar content, acidity, and secondary metabolism of the berries. (iii) Controlling the genome and its expression to allow breeding of better-adapted genotypes. High-density DNA genotyping can be used to select genotypes with specific interesting alleles but genomic selection is also a powerful method able to take into account the genetic information along the whole genome to predict a phenotype. Modern technologies are also able to generate mutations that are possibly interesting for generating new phenotypes but the most promising one is the direct editing of the genome at a precise location.

Targeting alternative splicing by RNAi: from the differential impact on splice variants to triggering artificial pre-mRNA splicing

Alternative splicing generates multiple transcript and protein isoforms from a single gene and controls transcript intracellular localization and stability by coupling to mRNA export and nonsense-mediated mRNA decay (NMD). RNA interference (RNAi) is a potent mechanism to modulate gene expression. However, its interactions with alternative splicing are poorly understood. We used artificial microRNAs (amiRNAs, also termed shRNAmiR) to knockdown all splice variants of selected target genes in Arabidopsis thaliana. We found that splice variants, which vary by their protein-coding capacity, subcellular localization and sensitivity to NMD, are affected differentially by an amiRNA, although all of them contain the target site. Particular transcript isoforms escape amiRNA-mediated degradation due to their nuclear localization. The nuclear and NMD-sensitive isoforms mask RNAi action in alternatively spliced genes. Interestingly, Arabidopsis SPL genes, which undergo alternative splicing and are targets of miR156, are regulated in the same manner. Moreover, similar results were obtained in mammalian cells using siRNAs, indicating cross-kingdom conservation of these interactions among RNAi and splicing isoforms. Furthermore, we report that amiRNA can trigger artificial alternative splicing, thus expanding the RNAi functional repertoire. Our findings unveil novel interactions between different post-transcriptional processes in defining transcript fates and regulating gene expression.

miRNA Targets: From Prediction Tools to Experimental Validation

MicroRNAs (miRNAs) are post-transcriptional regulators of gene expression in both animals and plants. By pairing to microRNA responsive elements (mREs) on target mRNAs, miRNAs play gene-regulatory roles, producing remarkable changes in several physiological and pathological processes. Thus, the identification of miRNA-mRNA target interactions is fundamental for discovering the regulatory network governed by miRNAs. The best way to achieve this goal is usually by computational prediction followed by experimental validation of these miRNA-mRNA interactions. This review summarizes the key strategies for miRNA target identification. Several tools for computational analysis exist, each with different approaches to predict miRNA targets, and their number is constantly increasing. The major algorithms available for this aim, including Machine Learning methods, are discussed, to provide practical tips for familiarizing with their assumptions and understanding how to interpret the results. Then, all the experimental procedures for verifying the authenticity of the identified miRNA-mRNA target pairs are described, including High-Throughput technologies, in order to find the best approach for miRNA validation. For each strategy, strengths and weaknesses are discussed, to enable users to evaluate and select the right approach for their interests.

Strontium stress disrupts miRNA biogenesis by reducing HYL1 protein levels in Arabidopsis

Strontium (Sr) is an emerging environmental pollutant that has become a major global concern after the nuclear accident at the Fukushima Daiichi Nuclear Power Plant in 2011. Although many studies have demonstrated the harmful effects of Sr on plant growth and development at the physiological level, knowledge regarding how plants sense and respond to Sr stress at the molecular level is limited. Recent studies have suggested that microRNAs (miRNAs) function as key regulators of plant growth and development as well as in the responses of plants to environmental stresses, including salinity, drought, cold, nutrient starvation, and heavy metals. In this study, we examined the global expression profile of miRNAs under Sr stress using small RNA sequencing analysis in Arabidopsis to better understand the molecular basis of plant responses to Sr stress. To identify specific Sr-responsive miRNAs, we performed comparative miRNA expression profiling analysis using control, CaCl₂-, and SrCl₂-treated seedlings. Compared to the control treatment, the expressions of most miRNAs were considerably decreased in the Sr-treated seedlings. However, under Sr stress, the expressions of primary miRNAs (pri-miRNAs) and their target genes were significantly increased; the protein levels of HYPONASTIC LEAVES 1 (HYL1), one of the core components of the microprocessor complex, were strongly reduced despite the increased HYL1 mRNA expression. In addition, hyl1-2 mutant plants were shown to be more sensitive to Sr stress than wild-type plants. Collectively, our results strongly suggested that Sr stress may be associated with the disruption of miRNA biogenesis by reducing the protein level of HYL1, which is required to maintain proper growth and development for plants. Our findings further indicated that some miRNAs may play important roles in plant responses to Sr stress.

Interactions of gene expression, alternative splicing, and DNA methylation in determining nodule identity

Soybean nodulation is a highly controlled process that involves complex gene regulation at both transcriptional and post-transcriptional levels. In the present study, we profiled gene expression changes, alternative splicing events, and DNA methylation patterns during nodule formation, development, and senescence. The transcriptome data uncovered key transcription patterns of nodule development that included 9669 core genes and 7302 stage-specific genes. Alternative splicing analysis uncovered a total of 2323 genes that undergo alternative splicing events in at least one nodule developmental stage, with activation of exon skipping and repression of intron retention being the most common splicing events in nodules compared to roots. Approximately 40% of the differentially spliced genes were also differentially expressed at the same nodule developmental stage, implying a substantial association between gene expression and alternative splicing. Genome-wide-DNA methylation analysis revealed dynamic changes in nodule methylomes that were specific to each nodule stage, occurred in a sequence-specific manner, and impacted the expression of 1864 genes. An attractive hypothesis raised by our data is that increased DNA methylation may contribute to the efficiency of alternative splicing. Together, our results provide intriguing insights into the associations between gene expression, alternative splicing, and DNA methylation that may shape transcriptome complexity and proteome specificity in developing soybean nodules.

Jasmonate induced alternative splicing responses in Arabidopsis

Jasmonate is an essential phytohormone regulating plant growth, development, and defense. Alternative splicing (AS) in jasmonate ZIM-domain (JAZ) repressors is well-characterized and plays an important role in jasmonate signaling regulation. However, it is unknown whether other genes in the jasmonate signaling pathway are regulated by AS. We explore the potential for AS regulation in three Arabidopsis genotypes (WT, jaz2, jaz7) in response to methyl jasmonate (MeJA) treatment with respect to: (a) differential AS, (b) differential miRNA targeted AS, and (c) AS isoforms with novel functions. AS events identified from transcriptomic data were validated with proteomic data. Protein interaction networks identified two genes, SKIP and ALY4 whose products have both DNA- and RNA-binding affinities, as potential key regulators mediating jasmonate signaling and AS regulation. We observed cases where AS alone, or AS and transcriptional regulation together, can influence gene expression in response to MeJA. Twenty-one genes contain predicted miRNA target sites subjected to AS, which implies that AS is coupled to miRNA regulation. We identified 30 cases where alternatively spliced isoforms may have novel functions. For example, AS of bHLH160 generates an isoform without a basic domain, which may convert it from an activator to a repressor. Our study identified potential key regulators in AS regulation of jasmonate signaling pathway. These findings highlight the importance of AS regulation in the jasmonate signaling pathway, both alone and in collaboration with other regulators.

SIGNIFICANCE STATEMENT

By exploring alternative splicing, we demonstrate its regulation in the jasmonate signaling pathway alone or in collaboration with other posttranscriptional regulations such as nonsense and microRNA-mediated decay. A signal transduction network model for alternative splicing in jasmonate signaling pathway was generated, contributing to our understanding for this important, prevalent, but relatively unexplored regulatory mechanism in plants.

Pituitary miRNAs target GHRHR splice variants to regulate GH synthesis by mediating different intracellular signalling pathways

Growth hormone (GH), whose synthesis and release are mainly regulated by intracellular signals mediated by growth hormone-releasing hormone receptor (GHRHR), is one of the major pituitary hormones and critical regulators of organism growth, metabolism, and immunoregulation. Pig GHRHR splice variants (SVs) may activate different signalling pathways via the variable C-terminal by alternative splicing, and SVs have the potential to change microRNA (miRNA) binding sites. In this study, we first confirmed the existence of pig GHRHR SVs (i.e., GHRHR, GHRHR SV1 and SV2) and demonstrated the inhibitory effects of critical pituitary miRNAs (i.e., let-7e and miR-328-5p) on GH synthesis and cell proliferation of primary pituitary cells. The SVs of GHRHR targeted by let-7e and miR-328-5p were predicted via bioinformatics analysis and verified by performing dual-luciferase reporter assays and detecting the expression of target transcripts. The differential responses of let-7e, and miR-328-5p to GH-releasing hormone and the changes in signalling pathways mediated by GHRHR suggested that let-7e and miR-328-5p were involved in GH synthesis mediated by GHRHR SVs, indicating that the two miRNAs played different roles by different ways. Finally, results showed that the protein coded by the GHRHR transcript regulated GH through the NO/NOS signalling pathway, whereas that coded by SV1 and SV2 regulated GH through the PKA/CREB signalling pathway, which was confirmed by the changes in signalling pathways after transfecting the expression vectors of GHRHR SVs to GH3 cells. To the best of our knowledge, this paper is the first to report pituitary miRNAs regulate GH synthesis by targeting the different SVs of GHRHR.

The interplay between microRNA and alternative splicing of linear and circular RNAs in eleven plant species

MOTIVATION

MicroRNA (miRNA) and alternative splicing (AS)-mediated post-transcriptional regulation has been extensively studied in most eukaryotes. However, the interplay between AS and miRNAs has not been explored in plants. To our knowledge, the overall profile of miRNA target sites in circular RNAs (circRNA) generated by alternative back splicing has never been reported previously. To address the challenge, we identified miRNA target sites located in alternatively spliced regions of the linear and circular splice isoforms using the up-to-date single-molecule real-time (SMRT) isoform sequencing (Iso-Seq) and Illumina sequencing data in eleven plant species.

RESULTS

In total, we identified 399 401 and 114 574 AS events from linear and circular RNAs, respectively. Among them, there were 64 781 and 41 146 miRNA target sites located in linear and circular AS region, respectively. In addition, we found 38 913 circRNAs to be overlapping with 45 648 AS events of its own parent isoforms, suggesting circRNA regulation of AS of linear RNAs by forming R-loop with the genomic locus. Here, we present a comprehensive database of miRNA targets in alternatively spliced linear and circRNAs (ASmiR) and a web server for deposition and identification of miRNA target sites located in the alternatively spliced region of linear and circular RNAs. This database is accompanied by an easy-to-use web query interface for meaningful downstream analysis. Plant research community can submit user-defined datasets to the web service to search AS regions harboring small RNA target sites. In conclusion, this study provides an unprecedented resource to understand regulatory relationships between miRNAs and AS in both gymnosperms and angiosperms.

AVAILABILITY AND IMPLEMENTATION

The readily accessible database and web-based tools are available at http://forestry.fafu.edu.cn/bioinfor/db/ASmiR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The expanding roles of Argonautes: RNA interference, splicing and beyond

Argonaute (AGO) protein family is highly conserved in eukaryotes and prokaryotes, reflecting its evolutionarily indispensible role in maintaining normal life cycle of the organisms. Small RNA-guided, AGO-dependent RNA interference (RNAi) is a well-studied pathway for gene expression regulation, which can be performed at transcriptional, posttranscriptional or translational level. In addition to RNAi, growing pieces of evidence point to a novel role of AGOs in pre-mRNA (messenger RNA precursor) splicing in animals. Many noncoding RNAs (ncRNAs) share common structural features with protein-coding genes, indicating that these ncRNAs might be subject to AGO-mediated splicing. Finally, we provide a comprehensive view that RNAi, transcription and RNA splicing are highly interactive processes, all of which involve several key factors such as AGOs. In this regard, the AGO proteins contribute to orchestrate an exquisite gene regulatory network in vivo. However, more research efforts are needed to reach a thorough understanding of the AGO activities.

Casting Away the Shadows: Elucidating the Role of Light-mediated Posttranscriptional Control in Plants

Light signals trigger precise changes in gene expression networks that activate distinctive developmental programs in plants. The transcriptome is shaped at different stages, both by the regulation of gene expression and also by posttranscriptional mechanisms that alter the sequence or abundance of the transcripts generated. Posttranscriptional mechanisms have attracted much interest in recent years with the advent of high-throughput technologies and bioinformatics tools. One such posttranscriptional process, alternative splicing, increases proteome diversity without increasing gene number by changing the function of individual proteins, while another, miRNA-mediated gene silencing, fine-tunes the amount of mRNA produced. The manner in which plants make use of these two crucial posttranscriptional mechanisms to respond to light and adapt to their environment is the focus of active research. In this review, we summarize the current knowledge of light-mediated posttranscriptional control in Arabidopsis thaliana and focus on the biological impact of the various posttranscriptional processes. We also discuss a potential cross talk between the alternative splicing and miRNA pathways, highlighting the complexity of light responsiveness.

Coupling of microRNA-directed phased small interfering RNA generation from long noncoding genes with alternative splicing and alternative polyadenylation in small RNA-mediated gene silencing

MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play vital regulatory roles in plant growth and development. Little is known about these small RNAs in litchi (Litchi chinensis), an economically important fruit crop widely cultivated in Southeast Asia. We profiled the litchi small RNA population with various deep-sequencing techniques and in-depth bioinformatic analyses. The genome-wide identification of miRNAs, their target genes, and phasiRNA-generating (PHAS) genes/loci showed that the function of miR482/2118 has expanded, relative to its canonical function. We also discovered that, for 29 PHAS loci, miRNA-mediated phasiRNA production was coupled with alternative splicing (AS) and alternative polyadenylation (APA). Most of these loci encoded long noncoding RNAs. An miR482/2118 targeted locus gave rise to four main transcript isoforms through AS/APA, and diverse phasiRNAs generated from these isoforms appeared to target long terminal repeat (LTR) retrotransposons and other unrelated genes. This coupling enables phasiRNA production from different exons of noncoding PHAS genes and yields diverse phasiRNA populations, both broadening and altering the range of downstream phasiRNA-regulated genes. Our results reveal the diversity of miRNA and phasiRNA in litchi, and demonstrate AS/APA as a new layer of regulation in small RNA-mediated gene silencing.

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.

MicroRNA1917 targets CTR4 splice variants to regulate ethylene responses in tomato

Ethylene perception is regulated by receptors, and the downstream protein CONSTITUTIVE TRIPLE RESPONSE1 is a key suppressor of ethylene signalling. The non-conserved tomato (Solanum lycopersicum) microRNA1917 (Sly-miR1917) mediates degradation of SlCTR4 splice variants (SlCTR4sv) but the molecular details of this pathway remain unknown. Sly-miR1917 and the targeted SlCTR4sv are ubiquitously expressed in all tomato organs. Overexpression of Sly-miR1917 enhances ethylene responses, including the triple response in etiolated seedlings, in the absence of ethylene, as well as epinastic petiole growth, accelerated pedicel abscission, and fruit ripening. Enhanced ethylene signalling in Sly-miR1917-overexpressing plants (1917-OE) is accompanied by up-regulation of ethylene biosynthesis and signalling genes, and increased ethylene emission. These phenotypes were recovered by repressing the positive ethylene regulator EIN2. Moreover, the Sly-miR1917-targeted SlCTR4 splice variant SlCTR4sv3, expressed specifically in the abscission zone, exhibited the opposite expression pattern to Sly-miR1917. Complementation of the Arabidopsis thaliana ctr-1 mutant and yeast two-hybrid and bimolecular fluorescence complementation assays suggested that SlCTR4sv3 functions in ethylene signalling. Co-expression of Sly-miR1917 and SlCTR4sv3 in Nicotiana benthamiana further suggested that Sly-miR1917 cleaves SlCTR4sv3 in vivo. Database homology searching revealed a Solanum tuberosum CTR-like splice variant containing a Sly-miR1917 binding sequence, and a homologue of mature Sly-miR1917 in potato, indicating a conserved function for miR1917 and the regulatory miRNA-mediated ethylene network in solanaceous species.

Bypassing miRNA-mediated gene regulation under drought stress: alternative splicing affects CSD1 gene expression

The binding site for miR398 in an isoform of Cu/Zn superoxide dismutase (CSD1) is eliminated by alternative splicing to bypass miR398-mediated gene down-regulation under drought stress. MicroRNA (miRNA) binding sites (MBSs) are frequently interrupted by introns and therefore require proper splicing to generate functional MBSs in target transcripts. MBSs can also be excluded during splicing of pre-messenger RNA, leading to different regulation among isoforms. Previous studies have shown that levels of Cu/Zn superoxide dismutase (CSD) are down-regulated by miR398. In this study, sequences and transcript levels of peanut CSD1 isoforms (AhCSD1-1, AhCSD1-2.1, and AhCSD1-2.2) were analyzed under the drought stress. Results demonstrated that a miR398 binding site is eliminated in AhCSD1-2.2 as a consequence of alternative splicing, which bypasses miRNA-mediated down-regulation under drought stress. This alternative isoform was not only identified in peanut but also in soybean and Arabidopsis. In addition, transgenic Arabidopsis plants expressing AhCSD1 were more tolerant to osmotic stress. We hypothesize that the level of AhCSD1 is increased to allow diverse plant responses to overcome environmental challenges even in the presence of increased miR398 levels. These findings suggest that studies on the role of alternatively spliced MBSs affecting transcript levels are important for understanding plant stress responses.

Identification and analysis of alternative splicing events in Phaseolus vulgaris and Glycine max

BACKGROUND

The vast diversification of proteins in eukaryotic cells has been related with multiple transcript isoforms from a single gene that result in alternative splicing (AS) of primary transcripts. Analysis of RNA sequencing data from expressed sequence tags and next generation RNA sequencing has been crucial for AS identification and genome-wide AS studies. For the identification of AS events from the related legume species Phaseolus vulgaris and Glycine max, 157 and 88 publicly available RNA-seq libraries, respectively, were analyzed.

RESULTS

We identified 85,570 AS events from P. vulgaris in 72% of expressed genes and 134,316 AS events in 70% of expressed genes from G. max. These were categorized in seven AS event types with intron retention being the most abundant followed by alternative acceptor and alternative donor, representing ~75% of all AS events in both plants. Conservation of AS events in homologous genes between the two species was analyzed where an overrepresentation of AS affecting 5'UTR regions was observed for certain types of AS events. The conservation of AS events was experimentally validated for 8 selected genes, through RT-PCR analysis. The different types of AS events also varied by relative position in the genes. The results were consistent in both species.

CONCLUSIONS

The identification and analysis of AS events are first steps to understand their biological relevance. The results presented here from two related legume species reveal high conservation, over ~15-20 MY of divergence, and may point to the biological relevance of AS.

Light Regulation of Alternative Pre-mRNA Splicing in Plants

Alternative splicing (AS) is a major post-transcriptional mechanism to enhance the diversity of proteome in response to environmental signals. Among the numerous external signals perceived by plants, light is the most crucial one. Plants utilize complex photoreceptor signaling networks to sense different light conditions and adjust their growth and development accordingly. Although light-mediated gene expression has been widely investigated, little is known regarding the mechanism of light affecting AS to modulate mRNA at the post-transcriptional level. In this minireview, we summarize current progresses on how light affects AS, and how sensory photoreceptors and retrograde signaling pathways may coordinately regulate AS of pre-mRNAs. In addition, we also discuss the possibility that AS of the mRNAs encoding photoreceptors may be involved in feedback control of AS. We hypothesize that light regulation of the expression and activity of splicing factors would be a major mechanism of light-mediated AS. The combination of genetic study and high-throughput analyses of AS and splicing complexes in response to light is likely to further advance our understanding of the molecular mechanisms underlying light control of AS and plant development.

Computational Prediction of MicroRNA Target Genes, Target Prediction Databases, and Web Resources

MicroRNA (miRNA) mediated silencing and repression of mRNA molecules requires complementary base pairing between the "seed" region of the miRNA and the "seed match" region of target mRNAs. While this mechanism is fairly well understood, accurate prediction of valid miRNA targets remains challenging due to factors such as imperfect sequence specificity, target site availability, and the thermodynamic stability of the mRNA structure itself. As knowledge of what genes are being targeted by each miRNA is arguably the most important facet of miRNA biology, many approaches have been developed to address the need for reliable prediction and ranking of putative targets, with most using a combination of various strategies such as evolutionary conservation, statistical inference, and distinct features of the target sequences themselves. This chapter reviews the pros and cons of a number of different prediction algorithms, showcases some databases that store experimentally validated miRNA targets, and also provides a case study that profiles some of the potential microRNA-mRNA interactions predicted by each methodology for various human genes.

Cross-talk between freezing response and signaling for regulatory transcriptions of MIR475b and its targets by miR475b promoter in Populus suaveolens

MicroRNAs (miRNAs) are small, non-coding RNAs that play important roles in post-transcriptional regulation of their target genes, yet the transcriptional regulation of plant miRNAs by promoter is poorly understood. Here, we firstly clone pri-miR475b cDNA and its native promoter from P. suaveolens, and characterize Psu-MIR475b as class-II gene transcribed by RNA polymerase II. By 5′ deletion analysis of Psu-miR475b promoter in a series of promoter-GUS chimeric vectors, we functionally identify three positive regulatory regions and multiple cis-acting elements responsible for Psu-miR475b promoter activity in response to freezing stress and exogenous hormone treatment. Moreover, the Psu-miR475b promoter activity displays a tissue-specific manner, negatively regulated by freezing stress and positively by MeJA, SA or GA treatment. Importantly, we comparatively analyze the time-course transcriptional profiles of Psu-miR475b and its targets in Psu-miR475b over-expression transgenic plants controlled by Psu-miR475b-specific promoter or CaMV 35S constitutive promoter, and explore the regulatory mechanism of Psu-miR475b promoter controlling transcriptional expressions of Psu-MIR475b and its targets in response to freezing stress and exogenous hormone treatment. Our results reveal that Psu-miR475b promoter-mediated transcriptions of Psu-MIR475b and its targets in response to freezing stress may be involved in a cross-talk between freezing response and stress signaling process.

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.

Genome-wide cataloging and analysis of alternatively spliced genes in cereal crops

Background

Protein functional diversity at the post-transcriptional level is regulated through spliceosome mediated pre-mRNA alternative splicing (AS) events and that has been widely demonstrated to be a key player in regulating the functional diversity in plants. Identification and analysis of AS genes in cereal crop plants are critical for crop improvement and understanding regulatory mechanisms.

Results

We carried out the comparative analyses of the functional landscapes of the AS using the consensus assembly of expressed sequence tags and available mRNA sequences in four cereal plants. We identified a total of 8,734 in Oryza sativa subspecies (ssp) japonica, 2,657 in O. sativa ssp indica, 3,971 in Sorghum bicolor, and 10,687 in Zea mays AS genes. Among the identified AS events, intron retention remains to be the dominant type accounting for 23.5 % in S. bicolor, and up to 55.8 % in O. sativa ssp indica. We identified a total of 887 AS genes that were conserved among Z. mays, S. bicolor, and O. sativa ssp japonica; and 248 AS genes were found to be conserved among all four studied species or ssp. Furthermore, we identified 53 AS genes conserved with Brachypodium distachyon. Gene Ontology classification of AS genes revealed functional assignment of these genes in many biological processes with diverse molecular functions.

Conclusions

AS is common in cereal plants. The AS genes identified in four cereal crops in this work provide the foundation for further studying the roles of AS in regulation of cereal plant growth and development. The data can be accessed at Plant Alternative Splicing Database (http://proteomics.ysu.edu/altsplice/).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1914-5) contains supplementary material, which is available to authorized users.

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.

Differential gene expression and alternative splicing between diploid and tetraploid watermelon

The exploitation of synthetic polyploids for producing seedless fruits is well known in watermelon. Tetraploid progenitors of triploid watermelon plants, compared with their diploid counterparts, exhibit wide phenotypic differences. Although many factors modulate alternative splicing (AS) in plants, the effects of autopolyploidization on AS are still unknown. In this study, we used tissues of leaf, stem, and fruit of diploid and tetraploid sweet watermelon to understand changes in gene expression and the occurrence of AS. RNA-sequencing analysis was performed along with reverse transcription quantitative PCR and rapid amplification of cDNA ends (RACE)-PCR to demonstrate changes in expression and splicing. All vegetative tissues except fruit showed an increased level of AS in the tetraploid watermelon throughout the growth period. The ploidy levels of diploids and the tetraploid were confirmed using a ploidy analyser. We identified 5362 and 1288 genes that were up- and downregulated, respectively, in tetraploid as compared with diploid plants. We further confirmed that 22 genes underwent AS events across tissues, indicating possibilities of generating different protein isoforms with altered functions of important transcription factors and transporters. Arginine biosynthesis, chlorophyllide synthesis, GDP mannose biosynthesis, trehalose biosynthesis, and starch and sucrose degradation pathways were upregulated in autotetraploids. Phloem protein 2, chloroplastic PGR5-like protein, zinc-finger protein, fructokinase-like 2, MYB transcription factor, and nodulin MtN21 showed AS in fruit tissues. These results should help in developing high-quality seedless watermelon and provide additional transcriptomic information related to other cucurbits.

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.

Alternative splicing at GYNNGY 5' splice sites: more noise, less regulation

Numerous eukaryotic genes are alternatively spliced. Recently, deep transcriptome sequencing has skyrocketed proportion of alternatively spliced genes; over 95% human multi-exon genes are alternatively spliced. One fundamental question is: are all these alternative splicing (AS) events functional? To look into this issue, we studied the most common form of alternative 5' splice sites-GYNNGYs (Y = C/T), where both GYs can function as splice sites. Global analyses suggest that splicing noise (due to stochasticity of splicing process) can cause AS at GYNNGYs, evidenced by higher AS frequency in non-coding than in coding regions, in non-conserved than in conserved genes and in lowly expressed than in highly expressed genes. However, ∼20% AS GYNNGYs in humans and ∼3% in mice exhibit tissue-dependent regulation. Consistent with being functional, regulated GYNNGYs are more conserved than unregulated ones. And regulated GYNNGYs have distinctive sequence features which may confer regulation. Particularly, each regulated GYNNGY comprises two splice sites more resembling each other than unregulated GYNNGYs, and has more conserved downstream flanking intron. Intriguingly, most regulated GYNNGYs may tune gene expression through coupling with nonsense-mediated mRNA decay, rather than encode different proteins. In summary, AS at GYNNGY 5' splice sites is primarily splicing noise, and secondarily a way of regulation.

Conservation and functional influence of alternative splicing in wood formation of Populus and Eucalyptus

BACKGROUND

Wood formation in tree species is regulated by multiple factors at various layers. Alternative splicing (AS) occurs within a large number of genes in wood formation. However, the functional implications and conservation of the AS occurrence are not well understood.

RESULTS

In this study, we profiled AS events in wood-forming tissues of Populus and Eucalyptus, and analyzed their functional implications as well as inter-species conservation. 28.3% and 20.7% of highly expressed transcripts in the developing xylem of Populus and Eucalyptus respectively were affected by AS events. Around 42% of the AS events resulted in changes to the original reading frame. 25.0% (in Populus) and 26.8% (in Eucalyptus) of the AS events may cause protein domain modification. In the process of wood formation, about 28% of AS-occurring genes were putative orthologs and 71 conserved AS events were identified in the two species.

CONCLUSION

Through analysis of AS events in developing xylem of two tree species, this study reveals an array of new information regarding AS occurrence and function in tree development.

Genome-wide analysis of alternative splicing in Zea mays: landscape and genetic regulation

Alternative splicing enhances transcriptome diversity in all eukaryotes and plays a role in plant tissue identity and stress adaptation. To catalog new maize (Zea mays) transcripts and identify genomic loci that regulate alternative splicing, we analyzed over 90 RNA-seq libraries from maize inbred lines B73 and Mo17, as well as Syn10 doubled haploid lines (progenies from B73 × Mo17). Transcript discovery was augmented with publicly available data from 14 maize tissues, expanding the maize transcriptome by more than 30,000 and increasing the percentage of intron-containing genes that undergo alternative splicing to 40%. These newly identified transcripts greatly increase the diversity of the maize proteome, sometimes coding for entirely different proteins compared with their most similar annotated isoform. In addition to increasing proteome diversity, many genes encoding novel transcripts gained an additional layer of regulation by microRNAs, often in a tissue-specific manner. We also demonstrate that the majority of genotype-specific alternative splicing can be genetically mapped, with cis-acting quantitative trait loci (QTLs) predominating. A large number of trans-acting QTLs were also apparent, with nearly half located in regions not shown to contain genes associated with splicing. Taken together, these results highlight the currently underappreciated role that alternative splicing plays in tissue identity and genotypic variation in maize.

Identification of differentially expressed genes that potentially confer pest resistance in transgenic ChIFN-γ tobacco

Chicken interferon-γ (ChIFN-γ) is both an inhibitor of viral replication and a regulator of numerous immunological functions. However, since little is known about the mechanisms underlying the insect-resistance of transgenic ChIFN-γ, a transgenic ChIFN-γ tobacco line was employed in the present study to explore this mechanism. A cDNA microarray (with 43,760 unigenes) was used to analyze the gene expression profiles of transgenic and wild-type (WT) tobacco leaves at two different growth stages. Compared with the WT, 1529 and 405 expressed sequence tags were significantly up- or downregulated on days 119 and 147, respectively. The differentially expressed genes (DEGs) are involved in metabolic regulation, cell division and differentiation, material synthesis and transport, signal transduction, and protein synthesis and degradation. Candidate genes that may increase cell density, thicken cell walls, promote secondary metabolite synthesis, and mediate plant hormone-induced resistance responses were used to identify the ChIFN-γ-mediated insect-resistance mechanisms. The insect-resistance of transgenic ChIFN-γ tobacco possibly involves unknown signaling pathways, which may directly or indirectly affect DEG expression-mediating genes. The degree of pest resistance increased as the plants grew. Three genes likely to be related to jasmonic acid- or salicylic acid-dependent plant defense responses, including CAF 1, Cop 8/CSN, and HD, are implicated in the insect-resistance of the transgenic plants. The mechanism of transgenic ChIFN-γ tobacco resistance also involves RPS20 and other genes that induce microRNA-based gene regulation. The ChIFN-γ-mediated DGEs contribute to insect-resistance in transgenic ChIFN-γ tobacco, which provides new insight into the role of ChIFN-γ.

A deep survey of alternative splicing in grape reveals changes in the splicing machinery related to tissue, stress condition and genotype

BACKGROUND

Alternative splicing (AS) significantly enhances transcriptome complexity. It is differentially regulated in a wide variety of cell types and plays a role in several cellular processes. Here we describe a detailed survey of alternative splicing in grape based on 124 SOLiD RNAseq analyses from different tissues, stress conditions and genotypes.

RESULTS

We used the RNAseq data to update the existing grape gene prediction with 2,258 new coding genes and 3,336 putative long non-coding RNAs. Several gene structures have been improved and alternative splicing was described for about 30% of the genes. A link between AS and miRNAs was shown in 139 genes where we found that AS affects the miRNA target site. A quantitative analysis of the isoforms indicated that most of the spliced genes have one major isoform and tend to simultaneously co-express a low number of isoforms, typically two, with intron retention being the most frequent alternative splicing event.

CONCLUSIONS

As described in Arabidopsis, also grape displays a marked AS tissue-specificity, while stress conditions produce splicing changes to a minor extent. Surprisingly, some distinctive splicing features were also observed between genotypes. This was further supported by the observation that the panel of Serine/Arginine-rich splicing factors show a few, but very marked differences between genotypes. The finding that a part the splicing machinery can change in closely related organisms can lead to some interesting hypotheses for evolutionary adaptation, that could be particularly relevant in the response to sudden and strong selective pressures.

Complexity of the alternative splicing landscape in plants

Alternative splicing (AS) of precursor mRNAs (pre-mRNAs) from multiexon genes allows organisms to increase their coding potential and regulate gene expression through multiple mechanisms. Recent transcriptome-wide analysis of AS using RNA sequencing has revealed that AS is highly pervasive in plants. Pre-mRNAs from over 60% of intron-containing genes undergo AS to produce a vast repertoire of mRNA isoforms. The functions of most splice variants are unknown. However, emerging evidence indicates that splice variants increase the functional diversity of proteins. Furthermore, AS is coupled to transcript stability and translation through nonsense-mediated decay and microRNA-mediated gene regulation. Widespread changes in AS in response to developmental cues and stresses suggest a role for regulated splicing in plant development and stress responses. Here, we review recent progress in uncovering the extent and complexity of the AS landscape in plants, its regulation, and the roles of AS in gene regulation. The prevalence of AS in plants has raised many new questions that require additional studies. New tools based on recent technological advances are allowing genome-wide analysis of RNA elements in transcripts and of chromatin modifications that regulate AS. Application of these tools in plants will provide significant new insights into AS regulation and crosstalk between AS and other layers of gene regulation.

Fine-tuning of microRNA-mediated repression of mRNA by splicing-regulated and highly repressive microRNA recognition element

BACKGROUND

MicroRNAs are very small non-coding RNAs that interact with microRNA recognition elements (MREs) on their target messenger RNAs. Varying the concentration of a given microRNA may influence the expression of many target proteins. Yet, the expression of a specific target protein can be fine-tuned by alternative cleavage and polyadenylation to the corresponding mRNA.

RESULTS

This study showed that alternative splicing of mRNA is a fine-tuning mechanism in the cellular regulatory network. The splicing-regulated MREs are often highly repressive MREs. This phenomenon was observed not only in the hsa-miR-148a-regulated DNMT3B gene, but also in many target genes regulated by hsa-miR-124, hsa-miR-1, and hsa-miR-181a. When a gene contains multiple MREs in transcripts, such as the VEGF gene, the splicing-regulated MREs are again the highly repressive MREs. Approximately one-third of the analysable human MREs in MiRTarBase and TarBase can potentially perform the splicing-regulated fine-tuning. Interestingly, the high (+30%) repression ratios observed in most of these splicing-regulated MREs indicate associations with functions. For example, the MRE-free transcripts of many oncogenes, such as N-RAS and others may escape microRNA-mediated suppression in cancer tissues.

CONCLUSIONS

This fine-tuning mechanism revealed associations with highly repressive MRE. Since high-repression MREs are involved in many important biological phenomena, the described association implies that splicing-regulated MREs are functional. A possible application of this observed association is in distinguishing functionally relevant MREs from predicted MREs.

Small changes in ambient temperature affect alternative splicing in Arabidopsis thaliana

Alternative splicing (AS) gives rise to multiple mRNA isoforms from the same gene, providing possibilities to regulate gene expression beyond the level of transcription. In a recent paper in Nucleic Acids Research we used a high resolution RT-PCR based panel to study changes in AS patterns in plants with altered levels of an hnRNP-like RNA-binding protein in Arabidopsis thaliana. Furthermore, we detected significant changes in AS patterns between different Arabidopsis ecotypes. Here we investigated how small changes in ambient temperature affect AS. We found significant changes in AS for 12 of 28 investigated events (43%) upon transfer of Arabidopsis plants from 20°C to 16°C and for 6 of the 28 investigated events (21%) upon transfer from 20°C to 24°C.

Identification of a novel microRNA (miRNA) from rice that targets an alternatively spliced transcript of the Nramp6 (Natural resistance-associated macrophage protein 6) gene involved in pathogen resistance

Plants have evolved efficient defence mechanisms to defend themselves from pathogen attack. Although many studies have focused on the transcriptional regulation of defence responses, less is known about the involvement of microRNAs (miRNAs) as post-transcriptional regulators of gene expression in plant immunity. This work investigates miRNAs that are regulated by elicitors from the blast fungus Magnaporthe oryzae in rice (Oryza sativa). Small RNA libraries were constructed from rice tissues and subjected to high-throughput sequencing for the identification of elicitor-responsive miRNAs. Target gene expression was examined by microarray analysis. Transgenic lines were used for the analysis of miRNA functioning in disease resistance. Elicitor treatment is accompanied by dynamic alterations in the expression of a significant number of miRNAs, including new members of annotated miRNAs. Novel miRNAs from rice are proposed. We report a new rice miRNA, osa-miR7695, which negatively regulates an alternatively spliced transcript of OsNramp6 (Natural resistance-associated macrophage protein 6). This novel miRNA experienced natural and domestication selection events during evolution, and its overexpression in rice confers pathogen resistance. This study highlights an miRNA-mediated regulation of OsNramp6 in disease resistance, whilst illustrating the existence of a novel regulatory network that integrates miRNA function and mRNA processing in plant immunity.

Comparative analysis of microRNA promoters in Arabidopsis and rice

Endogenously-encoded microRNAs (miRNAs) are a class of small regulatory RNAs that modulate gene expression at the post-transcriptional level. In plants, miRNAs have increasingly been identified by experiments based on next-generation sequencing (NGS). However, promoter organization is currently unknown for most plant miRNAs, which are transcribed by RNA polymerase II. This deficiency prevents a comprehensive understanding of miRNA-mediated gene networks. In this study, by analyzing full-length cDNA sequences related to miRNAs, we mapped transcription start sites (TSSs) for 62 and 55 miRNAs in Arabidopsis and rice, respectively. The average free energy (AFE) profiles in the vicinity of TSSs were studied for both species. By employing position weight matrices (PWM) for 99 plant cis-elements, we discovered that three cis-elements were over-represented in the miRNA promoters of both species, while four and ten cis-elements were over-represented in Arabidopsis only and in rice only. Thus, comparison of miRNA promoters between Arabidopsis and rice provides a new perspective for studying miRNA regulation in plants.

Identification and analysis of the proximal promoters of microRNA genes in Arabidopsis

Endogenous microRNAs (miRNAs) modulate gene expression at the post-transcription level. In plants, a vast majority of MIR genes are thought to be transcribed by RNA Polymerase II (Pol II). However, promoter organization is currently unknown for most plant MIR genes. This deficiency prevents a comprehensive understanding of miRNA-mediated gene networks. In this study, we performed Pol II chromatin immunoprecipitation (ChIP) analysis in Arabidopsis using a genome tiling microarray. Distinct Pol II binding was found at most MIR loci, which allowed prediction of the transcription start sites (TSSs) for 167 MIR genes in Arabidopsis that was validated by average free energy profiling. By employing 99 position weight matrices (PWM), we systematically scanned the regulatory regions upstream of the TSSs. We discovered eleven and ten cis-elements that are statistically over- and under-represented in MIR promoters, respectively. Thus, analysis of Pol II binding provides a new perspective for studying miRNAs in plants.

Construction of gene regulatory networks mediated by vegetative and reproductive stage-specific small RNAs in rice (Oryza sativa)

Although huge amounts of high-throughput sequencing (HTS) data are available, limited systematic analyses have been performed by integrating these valuable resources. Based on small RNA (sRNA), RNA and degradome HTS data, the sRNAs specifically expressed at vegetative and reproductive stages were identified separately in rice. Two distinct groups of sRNA HTS data, which were prepared during the vegetative and the reproductive stages, were utilized to extract stage-specific sRNAs. Degradome sequencing data were employed for sRNA target validation. RNA sequencing data were used to construct expression-based, sRNA-mediated networks. As a result, 26 microRNAs and 413 sRNAs were specifically expressed at the vegetative stage, and 79 microRNAs and 539 sRNAs were specifically expressed at the reproductive stage. In addition to the microRNAs, numerous stage-specific sRNAs enriched in ARGONAUTE1 showed great potential to perform cleavage-based repression on the targets. Several stage-specific sRNAs were indicated to result from the wobble effect of Dicer-like 1-mediated processing of microRNA precursors. The expression patterns of the sRNA targets, and the stage-specific cleavage signals strongly indicated the reliability of the constructed networks. A set of rice stage-specific sRNAs along with the regulatory cascades, which have great potential in regulating specific developmental stages, were provided for further investigation.

"And yet it moves": cell-to-cell and long-distance signaling by plant microRNAs

MicroRNAs (miRNAs) are key regulators of numerous genes in many eukaryotes. Some plant miRNAs are involved in developmental and physiological processes that require intercellular or inter-organ signaling. Movement of other small RNAs within plants has been established. Recent findings also demonstrate intercellular signaling by miRNAs and strongly support that a subset of these regulatory molecules move from one cell to another or over long distances. Phloem exudates contain diverse miRNAs and at least two of them, involved in responses to nutrient availability, are transmitted through grafts, indicating long-distance movement. Two miRNAs that regulate developmental processes are present in cells outside their domains of expression. Several results strongly support that one of them moves from cell to cell. Research on a mutant affected in plasmodesmata trafficking indicates that these intercellular channels are required for transmission of miRNA activity to adjacent cells. Moreover, ARGONAUTE proteins might be involved in the regulation of miRNA trafficking. Hypothesis on the features and mechanisms that may determine miRNA mobility are presented. Future challenges include identifying other mobile miRNAs; demonstrating that miRNA movement is required for non-cell autonomous action; and characterizing the mechanisms of translocation and genetic pathways that regulate miRNA movement.

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

MicroRNAs (miRNAs) have emerged as a class of regulators of gene expression through posttranscriptional degradation or translational repression in living cells. Increasing evidence points to the important relationship between miRNAs and environmental stress responses, but the regulatory mechanisms in plants are poorly understood. Here, we found that Arabidopsis thaliana intronic miR400 was cotranscribed with its host gene (At1g32583) and downregulated by heat treatment. Intriguingly, an alternative splicing (AS) event that occurred in the intron (306 bp) where MIR400 was located was specifically induced by heat stress. A 100 bp fragment was excised, and the remaining 206 bp intron containing MIR400 transcripts was retained in the host gene. The stress-induced AS event thus resulted in greater accumulation of miR400 primary transcripts and a low level of mature miR400. Together, these results provide the direct evidence that AS acts as a regulatory mechanism linking miRNAs and environmental stress in plants.

A reversed framework for the identification of microRNA-target pairs in plants

Most plant microRNAs (miRNAs) perform their repressive regulation through target cleavages. The resulting slicing sites on the target transcripts could be mapped by sequencing of the 3'-cleavage remnants, called degradome sequencing. The high sequence complementarity between miRNAs and their targets has greatly facilitated the development of the target prediction tools for plant miRNAs. The prediction results were then subjected to degradome sequencing data-based validation, through which numerous miRNA-target interactions have been extracted. However, some drawbacks are unavoidable when using this forward approach. Essentially, a known list of plant miRNAs should be obtained in advance of target prediction and validation. This becomes an obstacle to discover novel miRNAs and their targets. Here, after reviewing the current available algorithms for reverse identification of miRNA-target pairs in plants, a case study was performed by using a newly established framework with adjustable parameters. In this workflow, integration of degradome and ARGONAUTE 1-enriched small RNA sequencing data was recommended to do a relatively comprehensive and reliable search. Besides, several computational algorithms such as BLAST, target plots and RNA secondary structure prediction were used. The results demonstrated the prevalent utility of the reversed approach for uncovering miRNA-target interactions in plants.