Light Regulates Plant Alternative Splicing through the Control of Transcriptional Elongation.

Light makes carbon fixation possible, allowing plant and animal life on Earth. We have previously shown that light regulates alternative splicing in plants. Light initiates a chloroplast retrograde signaling that regulates nuclear alternative splicing of a subset of Arabidopsis thaliana transcripts. Here, we show that light promotes RNA polymerase II (Pol II) elongation in the affected genes, whereas in darkness, elongation is lower. These changes in transcription are consistent with elongation causing the observed changes in alternative splicing, as revealed by different drug treatments and genetic evidence. The light control of splicing and elongation is abolished in an Arabidopsis mutant defective in the transcription factor IIS (TFIIS). We report that the chloroplast control of nuclear alternative splicing in plants responds to the kinetic coupling mechanism found in mammalian cells, providing unique evidence that coupling is important for a whole organism to respond to environmental cues.

Genome-Wide Identification of microRNAs Associated with Starch Biosynthesis and Endosperm Development in Foxtail Millet.

Foxtail millet is one of the oldest crops, and its endosperm contains up to 70% of starch. Grain filling is an important starch accumulation process associated with foxtail millet yield and quality. However, the molecular mechanisms of grain filling in foxtail millet are relatively unclear. Here, we investigate the genes and regulated miRNAs associated with starch synthesis and metabolism in foxtail millet using high-throughput small RNA, mRNA and degradome sequencing. The regulation of starch synthesis and quality is carried out mainly at the 15 DAA to 35 DAA stage during grain filling. The DEGs between waxy and non-waxy foxtail millet were significant, especially for GBSS. Additionally, ptc-miR169i_R+2_1ss21GA, fve-miR396e_L-1R+1, mtr-miR162 and PC-5p-221_23413 regulate the expression of genes associated with the starch synthesis pathway in foxtail millet. This study provides new insights into the molecular mechanisms of starch synthesis and quality formation in foxtail millet.

Computational methods for annotation of plant regulatory non-coding RNAs using RNA-seq.

Plant transcriptome encompasses numerous endogenous, regulatory non-coding RNAs (ncRNAs) that play a major biological role in regulating key physiological mechanisms. While studies have shown that ncRNAs are extremely diverse and ubiquitous, the functions of the vast majority of ncRNAs are still unknown. With ever-increasing ncRNAs under study, it is essential to identify, categorize and annotate these ncRNAs on a genome-wide scale. The use of high-throughput RNA sequencing (RNA-seq) technologies provides a broader picture of the non-coding component of transcriptome, enabling the comprehensive identification and annotation of all major ncRNAs across samples. However, the detection of known and emerging class of ncRNAs from RNA-seq data demands complex computational methods owing to their unique as well as similar characteristics. Here, we discuss major plant endogenous, regulatory ncRNAs in an RNA sample followed by computational strategies applied to discover each class of ncRNAs using RNA-seq. We also provide a collection of relevant software packages and databases to present a comprehensive bioinformatics toolbox for plant ncRNA researchers. We assume that the discussions in this review will provide a rationale for the discovery of all major categories of plant ncRNAs.

Oryza sativa RNA-Dependent RNA Polymerase 6 Contributes to Double-Strand Break Formation in Meiosis.

RNA-dependent RNA polymerase 6 (RDR6) is a core component of the small RNA biogenesis pathway, but its function in meiosis is unclear. Here, we report a new allele of OsRDR6 (Osrdr6-meiosis [Osrdr6-mei]), which causes meiosis-specific phenotypes in rice (Oryza sativa). In Osrdr6-mei, meiotic double-strand break (DSB) formation is partially blocked. We created a biallelic mutant with more severe phenotypes, Osrdr6-bi, by crossing Osrdr6-mei with a knockout mutant, Osrdr6-edit In Osrdr6-bi meiocytes, 24 univalents were observed, and no histone H2AX phosphorylation foci were detected. Compared with the wild type, the number of 21-nucleotide small RNAs in Osrdr6-mei was dramatically lower, while the number of 24-nucleotide small RNAs was significantly higher. Thousands of differentially methylated regions (DMRs) were discovered in Osrdr6-mei, implying that OsRDR6 plays an important role in DNA methylation. There were 457 genes downregulated in Osrdr6-mei, including three genes, CENTRAL REGION COMPONENT1, P31 comet , and O. sativa SOLO DANCERS, related to DSB formation. Interestingly, the downregulated genes were associated with a high level of 24-nucleotide small RNAs but less strongly associated with DMRs. Therefore, we speculate that the alteration in expression of small RNAs in Osrdr6 mutants leads to the defects in DSB formation during meiosis, which might not be directly dependent on RNA-directed DNA methylation.

Evolution of plant telomerase RNAs: farther to the past, deeper to the roots.

The enormous sequence heterogeneity of telomerase RNA (TR) subunits has thus far complicated their characterization in a wider phylogenetic range. Our recent finding that land plant TRs are, similarly to known ciliate TRs, transcribed by RNA polymerase III and under the control of the type-3 promoter, allowed us to design a novel strategy to characterize TRs in early diverging Viridiplantae taxa, as well as in ciliates and other Diaphoretickes lineages. Starting with the characterization of the upstream sequence element of the type 3 promoter that is conserved in a number of small nuclear RNAs, and the expected minimum TR template region as search features, we identified candidate TRs in selected Diaphoretickes genomes. Homologous TRs were then used to build covariance models to identify TRs in more distant species. Transcripts of the identified TRs were confirmed by transcriptomic data, RT-PCR and Northern hybridization. A templating role for one of our candidates was validated in Physcomitrium patens. Analysis of secondary structure demonstrated a deep conservation of motifs (pseudoknot and template boundary element) observed in all published TRs. These results elucidate the evolution of the earliest eukaryotic TRs, linking the common origin of TRs across Diaphoretickes, and underlying evolutionary transitions in telomere repeats.

Histone methyltransferase ATX1 dynamically regulates fiber secondary cell wall biosynthesis in Arabidopsis inflorescence stem.

Secondary wall thickening in the sclerenchyma cells is strictly controlled by a complex network of transcription factors in vascular plants. However, little is known about the epigenetic mechanism regulating secondary wall biosynthesis. In this study, we identified that ARABIDOPSIS HOMOLOG of TRITHORAX1 (ATX1), a H3K4-histone methyltransferase, mediates the regulation of fiber cell wall development in inflorescence stems of Arabidopsis thaliana. Genome-wide analysis revealed that the up-regulation of genes involved in secondary wall formation during stem development is largely coordinated by increasing level of H3K4 tri-methylation. Among all histone methyltransferases for H3K4me3 in Arabidopsis, ATX1 is markedly increased during the inflorescence stem development and loss-of-function mutant atx1 was impaired in secondary wall thickening in interfascicular fibers. Genetic analysis showed that ATX1 positively regulates secondary wall deposition through activating the expression of secondary wall NAC master switch genes, SECONDARY WALL-ASSOCIATED NAC DOMAIN PROTEIN1 (SND1) and NAC SECONDARY WALL THICKENING PROMOTING FACTOR1 (NST1). We further identified that ATX1 directly binds the loci of SND1 and NST1, and activates their expression by increasing H3K4me3 levels at these loci. Taken together, our results reveal that ATX1 plays a key role in the regulation of secondary wall biosynthesis in interfascicular fibers during inflorescence stem development of Arabidopsis.

DNA nanostructures coordinate gene silencing in mature plants.

Delivery of biomolecules to plants relies on Agrobacterium infection or biolistic particle delivery, the former of which is amenable only to DNA delivery. The difficulty in delivering functional biomolecules such as RNA to plant cells is due to the plant cell wall, which is absent in mammalian cells and poses the dominant physical barrier to biomolecule delivery in plants. DNA nanostructure-mediated biomolecule delivery is an effective strategy to deliver cargoes across the lipid bilayer of mammalian cells; however, nanoparticle-mediated delivery without external mechanical aid remains unexplored for biomolecule delivery across the cell wall in plants. Herein, we report a systematic assessment of different DNA nanostructures for their ability to internalize into cells of mature plants, deliver siRNAs, and effectively silence a constitutively expressed gene in Nicotiana benthamiana leaves. We show that nanostructure internalization into plant cells and corresponding gene silencing efficiency depends on the DNA nanostructure size, shape, compactness, stiffness, and location of the siRNA attachment locus on the nanostructure. We further confirm that the internalization efficiency of DNA nanostructures correlates with their respective gene silencing efficiencies but that the endogenous gene silencing pathway depends on the siRNA attachment locus. Our work establishes the feasibility of biomolecule delivery to plants with DNA nanostructures and both details the design parameters of importance for plant cell internalization and also assesses the impact of DNA nanostructure geometry for gene silencing mechanisms.

RNA polymerase II-independent recruitment of SPT6L at transcription start sites in Arabidopsis.

SPT6 is a conserved elongation factor that is associated with phosphorylated RNA polymerase II (RNAPII) during transcription. Recent transcriptome analysis in yeast mutants revealed its potential role in the control of transcription initiation at genic promoters. However, the mechanism by which this is achieved and how this is linked to elongation remains to be elucidated. Here, we present the genome-wide occupancy of Arabidopsis SPT6-like (SPT6L) and demonstrate its conserved role in facilitating RNAPII occupancy across transcribed genes. We also further demonstrate that SPT6L enrichment is unexpectedly shifted, from gene body to transcription start site (TSS), when its association with RNAPII is disrupted. Protein domains, required for proper function and enrichment of SPT6L on chromatin, are subsequently identified. Finally, our results suggest that recruitment of SPT6L at TSS is indispensable for its spreading along the gene body during transcription. These findings provide new insights into the mechanisms underlying SPT6L recruitment in transcription and shed light on the coordination between transcription initiation and elongation.

Genome-Wide Identification of Copper Stress-Regulated and Novel MicroRNAs in Mulberry Leaf.

Copper (Cu) is an essential trace element for plant growth and development. It is widely involved in respiration, photosynthesis, pollen formation, and other biological processes. Therefore, low or excessive copper causes damage to plants. Mulberry is an essential perennial economic tree. At present, research on the abiotic stress responses in mulberry is mainly focused on the identification of resistant germplasm resources and cloning of resistant genes. In contrast, studies on the resistance function of microRNAs and the regulatory gene responses to stress are rare. In this study, small RNA libraries (control and copper stressed) were constructed from mulberry leaf RNA. High-throughput sequencing and screening were employed, a total of 65 known miRNAs and 78 predicted novel mature miRNAs were identified, among which 40 miRNAs were differentially expressed under copper stress. Subsequently, expression patterns were verified for 14 miRNAs by real-time fluorescence quantitative PCR (qPCR). The target genes of miRNAs were validated by 5' RLM-RACE. Our results provide the bases for further study on the molecular mechanism of copper stress regulation in mulberry.

Grape Small Auxin Upregulated RNA (SAUR) 041 Is a Candidate Regulator of Berry Size in Grape.

Grape (Vitis vinifera) is an important horticultural crop that can be used to make juice and wine. However, the small size of the berry limits its yield. Cultivating larger berry varieties can be an effective way to solve this problem. As the largest family of auxin early response genes, SAUR (small auxin upregulated RNA) plays an important role in the growth and development of plants. Berry size is one of the important factors that determine grape quality. However, the SAUR gene family's function in berry size of grape has not been studied systematically. We identified 60 SAUR members in the grape genome and divided them into 12 subfamilies based on phylogenetic analysis. Subsequently, we conducted a comprehensive and systematic analysis on the SAUR gene family by analyzing distribution of key amino acid residues in the domain, structural features, conserved motifs, and protein interaction network, and combined with the heterologous expression in Arabidopsis and tomato. Finally, the member related to grape berry size in SAUR gene family were screened. This genome-wide study provides a systematic analysis of grape SAUR gene family, further understanding the potential functions of candidate genes, and provides a new idea for grape breeding.

Plant transfer RNA-derived fragments: Biogenesis and functions.

Processing of mature transfer RNAs (tRNAs) produces complex populations of tRNA-derived fragments (tRFs). Emerging evidence shows that tRFs have important functions in bacteria, animals, and plants. Here, we review recent advances in understanding plant tRFs, focusing on their biological and cellular functions, such as regulating stress responses, mediating plant-pathogen interactions, and modulating post-transcriptional gene silencing and translation. We also review sequencing strategies and bioinformatics resources for studying tRFs in plants. Finally, we discuss future directions for plant tRF research, which will expand our knowledge of plant non-coding RNAs.

Keep calm and carry on: miRNA biogenesis under stress.

MicroRNAs (miRNAs) are major post-transcriptional regulators of gene expression. Their biogenesis relies on the cleavage of longer precursors by a nuclear localized processing machinery. The evolutionary preference of plant miRNAs to silence transcription factors turned these small molecules into key actors during growth and adaptive responses. Furthermore, during their life cycle plants are subject to changes in the environmental conditions surrounding them. In order to face these changes, plants display unique adaptive capacities based on an enormous developmental plasticity, where miRNAs play central roles. Many individual miRNAs have been shown to modulate the plant response to different environmental cues and stresses. In the last few years, increasing evidence has shown that not only individual genes encoding miRNAs but also the miRNA pathway as a whole is subject to regulation in response to external stimulus. In this review, we discuss the current knowledge about the miRNA pathway. We dissect the pathway to analyze the events leading to the generation of these small RNAs and emphasize the regulation of core components of the miRNA biogenesis machinery.

Genomic identification and expression analysis of the BBX transcription factor gene family in Petunia hybrida.

The B-box proteins (BBXs) are a class of zinc finger transcription factors containing one or two B-BOX domains that play important roles in plant growth, development and stress response. The petunia (Petunia hybrida) is a model ornamental plant, and its draft genome has been published. However, no systematic study of the BBX gene family in Petunia has been reported. In this study, a total of 28 BBX members from the Petunia genome were identified. We performed analyses of their phylogenetic relationships, structures, conserved motifs, promoter regions, and expression patterns. Based on the phylogenetic relationship, the PhBBXs were divided into six groups. Analysis of the gene structures and conserved motifs further confirmed the closer relationships in each group. Based on the RNA-seq data, the transcript abundance of PhBBXs in different tissues were divided into two major groups. The analysis of cis-elements showed that many stress responsive elements appeared in the promoter region of most PhBBX genes. The stress response patterns of PhBBXs were detected under drought, salinity, cold and heat treatments. Based on the RNA-seq data, we found that 3 genes responded to drought, 8 genes responded to salt, 18 genes responded to cold, and 15 genes responded to heat. In conclusion, this study may facilitate further functional studies of BBXs in Petunia.

SRA- and SET-domain-containing proteins link RNA polymerase V occupancy to DNA methylation.

RNA-directed DNA methylation in Arabidopsis thaliana depends on the upstream synthesis of 24-nucleotide small interfering RNAs (siRNAs) by RNA POLYMERASE IV (Pol IV) and downstream synthesis of non-coding transcripts by Pol V. Pol V transcripts are thought to interact with siRNAs which then recruit DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2) to methylate DNA. The SU(VAR)3-9 homologues SUVH2 and SUVH9 act in this downstream step but the mechanism of their action is unknown. Here we show that genome-wide Pol V association with chromatin redundantly requires SUVH2 and SUVH9. Although SUVH2 and SUVH9 resemble histone methyltransferases, a crystal structure reveals that SUVH9 lacks a peptide-substrate binding cleft and lacks a properly formed S-adenosyl methionine (SAM)-binding pocket necessary for normal catalysis, consistent with a lack of methyltransferase activity for these proteins. SUVH2 and SUVH9 both contain SRA (SET- and RING-ASSOCIATED) domains capable of binding methylated DNA, suggesting that they function to recruit Pol V through DNA methylation. Consistent with this model, mutation of DNA METHYLTRANSFERASE 1 (MET1) causes loss of DNA methylation, a nearly complete loss of Pol V at its normal locations, and redistribution of Pol V to sites that become hypermethylated. Furthermore, tethering SUVH9 [corrected] with a zinc finger to an unmethylated site is sufficient to recruit Pol V and establish DNA methylation and gene silencing. These results indicate that Pol V is recruited to DNA methylation through the methyl-DNA binding SUVH2 and SUVH9 proteins, and our mechanistic findings suggest a means for selectively targeting regions of plant genomes for epigenetic silencing.

In vitro transcription activities of Pol IV, Pol V, and RDR2 reveal coupling of Pol IV and RDR2 for dsRNA synthesis in plant RNA silencing.

In Arabidopsis, RNA-dependent DNA methylation and transcriptional silencing involves three nuclear RNA polymerases that are biochemically undefined: the presumptive DNA-dependent RNA polymerases Pol IV and Pol V and the putative RNA-dependent RNA polymerase RDR2. Here we demonstrate their RNA polymerase activities in vitro. Unlike Pol II, Pols IV and V require an RNA primer, are insensitive to α-amanitin, and differ in their ability to displace the nontemplate DNA strand during transcription. Biogenesis of 24 nt small interfering RNAs (siRNAs), which guide cytosine methylation to corresponding sequences, requires both Pol IV and RDR2, which physically associate in vivo. Whereas Pol IV does not require RDR2 for activity, RDR2 is nonfunctional in the absence of associated Pol IV. These results suggest that the physical and mechanistic coupling of Pol IV and RDR2 results in the channeled synthesis of double-stranded precursors for 24 nt siRNA biogenesis.

Identification of novel phasiRNAs loci on long non-coding RNAs in Arabidopsis thaliana.

Long non-coding RNAs (lncRNAs) are the "dark matters"involved in gene regulation with complex mechanisms. However, the functions of most lncRNAs remain to be determined. Our previous work revealed a massive number of degradome-supported cleavage signatures on Arabidopsis lncRNAs. Some of them have been confirmed associated with miRNAs-like sRNAs production, while others without long stem structure remain unexplored. A systematical search for phasiRNAs generating ability of these lncRNAs was conducted. Eight novel small RNA triggered lncRNA-phasiRNA pathways were discovered and three of them were found to be conserved in Arabidopsis, Oryza sativa, Glycine max and Gossypium hirsutum. Besides, Five novel ta-siRNAs derived from these lncRNAs were further identified to be involved in the regulation of plant development, stress responses and aromatic amino acids synthesis. These results substantially expanded the gene regulation mechanisms of lncRNAs.

Plant Elongator-Protein Complex of Diverse Activities Regulates Growth, Development, and Immune Responses.

Contrary to the conserved Elongator composition in yeast, animals, and plants, molecular functions and catalytic activities of the complex remain controversial. Elongator was identified as a component of elongating RNA polymerase II holoenzyme in yeast, animals, and plants. Furthermore, it was suggested that Elonagtor facilitates elongation of transcription via histone acetyl transferase activity. Accordingly, phenotypes of Arabidopsis elo mutants, which show development, growth, or immune response defects, correlate with transcriptional downregulation and the decreased histone acetylation in the coding regions of crucial genes. Plant Elongator was also implicated in other processes: transcription and processing of miRNA, regulation of DNA replication by histone acetylation, and acetylation of alpha-tubulin. Moreover, tRNA modification, discovered first in yeast and confirmed in plants, was claimed as the main activity of Elongator, leading to specificity in translation that might also result indirectly in a deficiency in transcription. Heterologous overexpression of individual Arabidopsis Elongator subunits and their respective phenotypes suggest that single Elongator subunits might also have another function next to being a part of the complex. In this review, we shall present the experimental evidence of all molecular mechanisms and catalytic activities performed by Elongator in nucleus and cytoplasm of plant cells, which might explain how Elongator regulates growth, development, and immune responses.

Natural antisense transcripts are significantly involved in regulation of drought stress in maize.

Natural antisense transcripts (NATs) are a prominent and complex class of regulatory RNAs. Using strand-specific RNA sequencing, we identified 1769 sense and antisense transcript pairs (NAT pairs) in two maize inbreds with different sensitivity to drought, as well as in two derivative recombination inbred lines (RILs). A significantly higher proportion of NATs relative to non-NATs are specifically expressed under water stress (WS). Surprisingly, expression of sense and antisense transcripts produced by NAT pairs is significantly correlated, particularly under WS. We found an unexpected large proportion of NATs with protein coding potential, as estimated by ribosome release scores. Small RNAs significantly accumulate within NAT pairs, with 21 nt smRNA particularly enriched in overlapping regions of these pairs of genes. The abundance of these smRNAs is significantly altered in the leafbladeless1 mutant, suggesting that these genes may be regulated by the tasiRNA pathway. Further, NATs are significantly hypomethylated and include fewer transposable element sequences relative to non-NAT genes. NAT gene regions also exhibit higher levels of H3K36me3, H3K9ac, and H3K4me3, but lower levels of H3K27me3, indicating that NAT gene pairs generally exhibit an open chromatin configuration. Finally, NAT pairs in 368 diverse maize inbreds and 19 segregating populations were specifically enriched for polymorphisms associated with drought tolerance. Taken together, the data highlight the potential impact of that small RNAs and histone modifications have in regulation of NAT expression, and the significance of NATs in response to WS.

Unravelling the MicroRNA-Mediated Gene Regulation in Developing Pongamia Seeds by High-Throughput Small RNA Profiling.

Pongamia (Millettia pinnata syn. Pongamia pinnata) is a multipurpose biofuel tree which can withstand a variety of abiotic stresses. Commercial applications of Pongamia trees may substantially benefit from improvements in their oil-seed productivity, which is governed by complex regulatory mechanisms underlying seed development. MicroRNAs (miRNAs) are important molecular regulators of plant development, while relatively little is known about their roles in seed development, especially for woody plants. In this study, we identified 236 conserved miRNAs within 49 families and 143 novel miRNAs via deep sequencing of Pongamia seeds sampled at three developmental phases. For these miRNAs, 1327 target genes were computationally predicted. Furthermore, 115 differentially expressed miRNAs (DEmiRs) between successive developmental phases were sorted out. The DEmiR-targeted genes were preferentially enriched in the functional categories associated with DNA damage repair and photosynthesis. The combined analyses of expression profiles for DEmiRs and functional annotations for their target genes revealed the involvements of both conserved and novel miRNA-target modules in Pongamia seed development. Quantitative Real-Time PCR validated the expression changes of 15 DEmiRs as well as the opposite expression changes of six targets. These results provide valuable miRNA candidates for further functional characterization and breeding practice in Pongamia and other oilseed plants.

Enhanced Senescence Process is the Major Factor Stopping Spike Differentiation of Wheat Mutant ptsd1.

Complete differentiation of the spikes guarantees the final wheat (Triticum aestivum L.) grain yield. A unique wheat mutant that prematurely terminated spike differentiation (ptsd1) was obtained from cultivar Guomai 301 treated with ethyl methane sulfonate (EMS). The molecular mechanism study on ptsd1 showed that the senescence-associated genes (SAGs) were highly expressed, and spike differentiation related homeotic genes were depressed. Cytokinin signal transduction was weakened and ethylene signal transduction was enhanced. The enhanced expression of Ca2+ signal transduction related genes and the accumulation of reactive oxygen species (ROS) caused the upper spikelet cell death. Many genes in the WRKY, NAC and ethylene response factor (ERF) transcription factor (TF) families were highly expressed. Senescence related metabolisms, including macromolecule degradation, nutrient recycling, as well as anthocyanin and lignin biosynthesis, were activated. A conserved tae-miR164 and a novel-miR49 and their target genes were extensively involved in the senescence related biological processes in ptsd1. Overall, the abnormal phytohormone homeostasis, enhanced Ca2+ signaling and activated senescence related metabolisms led to the spikelet primordia absent their typical meristem characteristics, and ultimately resulted in the phenotype of ptsd1.