Protocol to measure ribosome density along mRNA transcripts of Arabidopsis thaliana tissues using Ribo-seq.

Ribosome profiling (Ribo-seq) measures ribosome density along messenger RNA (mRNA) transcripts and is used to estimate the "translational fitness" of a given mRNA in response to environmental or developmental cues with high resolution. Here, we describe a protocol for Ribo-seq in plants adapted for the model plant Arabidopsis thaliana. We describe steps for lysis and nucleolytic digestion and ribosome footprinting. We then detail library construction, sequencing, and data analysis.

A GFP splicing reporter in a coilin mutant background reveals links between alternative splicing, siRNAs, and coilin function in Arabidopsis thaliana.

Coilin is a scaffold protein essential for the structure of Cajal bodies, which are nucleolar-associated, nonmembranous organelles that coordinate the assembly of nuclear ribonucleoproteins (RNPs) including spliceosomal snRNPs. To study coilin function in plants, we conducted a genetic suppressor screen using a coilin (coi1) mutant in Arabidopsis thaliana and performed an immunoprecipitation-mass spectrometry analysis on coilin protein. The coi1 mutations modify alternative splicing of a GFP reporter gene, resulting in a hyper-GFP phenotype in young coi1 seedlings relative to the intermediate wild-type level. As shown here, this hyper-GFP phenotype is extinguished in older coi1 seedlings by posttranscriptional gene silencing triggered by siRNAs derived from aberrant splice variants of GFP pre-mRNA. In the coi1 suppressor screen, we identified suppressor mutations in WRAP53, a putative coilin-interacting protein; SMU2, a predicted splicing factor; and ZCH1, an incompletely characterized zinc finger protein. These suppressor mutations return the hyper-GFP fluorescence of young coi1 seedlings to the intermediate wild-type level. Additionally, coi1 zch1 mutants display more extensive GFP silencing and elevated levels of GFP siRNAs, suggesting the involvement of wild-type ZCH1 in siRNA biogenesis or stability. The immunoprecipitation-mass spectrometry analysis reinforced the roles of coilin in pre-mRNA splicing, nucleolar chromatin structure, and rRNA processing. The participation of coilin in these processes, at least some of which incorporate small RNAs, supports the hypothesis that coilin provides a chaperone for small RNA trafficking. Our study demonstrates the usefulness of the GFP splicing reporter for investigating alternative splicing, ribosome biogenesis, and siRNA-mediated silencing in the context of coilin function.

The Arabidopsis Deubiquitylase OTU5 Suppresses Flowering by Histone Modification-Mediated Activation of the Major Flowering Repressors FLC, MAF4, and MAF5.

Distinct phylogeny and substrate specificities suggest that 12 Arabidopsis Ovarian Tumor domain-containing (OTU) deubiquitinases participate in conserved or plant-specific functions. The otu5-1 null mutant displayed a pleiotropic phenotype, including early flowering, mimicking that of mutants harboring defects in subunits (e.g., ARP6) of the SWR1 complex (SWR1c) involved in histone H2A.Z deposition. Transcriptome and RT-qPCR analyses suggest that downregulated FLC and MAF4-5 are responsible for the early flowering of otu5-1. qChIP analyses revealed a reduction and increase in activating and repressive histone marks, respectively, on FLC and MAF4-5 in otu5-1. Subcellular fractionation, GFP-fusion expression, and MNase treatment of chromatin showed that OTU5 is nucleus-enriched and chromatin-associated. Moreover, OTU5 was found to be associated with FLC and MAF4-5. The OTU5-associated protein complex(es) appears to be distinct from SWR1c, as the molecular weights of OTU5 complex(es) were unaltered in arp6-1 plants. Furthermore, the otu5-1 arp6-1 double mutant exhibited synergistic phenotypes, and H2A.Z levels on FLC/MAF4-5 were reduced in arp6-1 but not otu5-1. Our results support the proposition that Arabidopsis OTU5, acting independently of SWR1c, suppresses flowering by activating FLC and MAF4-5 through histone modification. Double-mutant analyses also indicate that OTU5 acts independently of the HUB1-mediated pathway, but it is partially required for FLC-mediated flowering suppression in autonomous pathway mutants and FRIGIDA-Col.

Transcriptomic Analysis Suggests Auxin Regulation in Dorsal-Ventral Petal Asymmetry of Wild Progenitor Sinningia speciosa.

The establishment of dorsal-ventral (DV) petal asymmetry is accompanied by differential growth of DV petal size, shape, and color differences, which enhance ornamental values. Genes involved in flower symmetry in Sinningia speciosa have been identified as CYCLOIDEA (SsCYC), but which gene regulatory network (GRN) is associated with SsCYC to establish DV petal asymmetry is still unknown. To uncover the GRN of DV petal asymmetry, we identified 630 DV differentially expressed genes (DV-DEGs) from the RNA-Seq of dorsal and ventral petals in the wild progenitor, S. speciosa 'ES'. Validated by qRT-PCR, genes in the auxin signaling transduction pathway, SsCYC, and a major regulator of anthocyanin biosynthesis were upregulated in dorsal petals. These genes correlated with a higher endogenous auxin level in dorsal petals, with longer tube length growth through cell expansion and a purple dorsal color. Over-expression of SsCYC in Nicotiana reduced petal size by regulating cell growth, suggesting that SsCYC also controls cell expansion. This suggests that auxin and SsCYC both regulate DV petal asymmetry. Transiently over-expressed SsCYC, however, could not activate most major auxin signaling genes, suggesting that SsCYC may not trigger auxin regulation. Whether auxin can activate SsCYC or whether they act independently to regulate DV petal asymmetry remains to be explored in the future.

An improved repertoire of splicing variants and their potential roles in Arabidopsis photomorphogenic development.

BACKGROUND: Light switches on the photomorphogenic development of young plant seedlings, allowing young seedlings to acquire photosynthetic capacities and gain survival fitness. Light regulates gene expression at all levels of the central dogma, including alternative splicing (AS) during the photomorphogenic development. However, accurate determination of full-length (FL) splicing variants has been greatly hampered by short-read RNA sequencing technologies. RESULT: In this study, we adopt PacBio isoform sequencing (Iso-seq) to overcome the limitation of the short-read RNA-seq technologies. Normalized cDNA libraries used for Iso-seq allows for comprehensive and effective identification of FL AS variants. Our analyses reveal more than 30,000 splicing variant models from approximately 16,500 gene loci and additionally identify approximately 700 previously unannotated genes. Among the variants, approximately 12,000 represent new gene models. Intron retention (IR) is the most frequently observed form of variants, and many IR-containing AS variants show evidence of engagement in translation. Our study reveals the formation of heterodimers of transcription factors composed of annotated and IR-containing AS variants. Moreover, transgenic plants overexpressing the IR forms of two B-BOX DOMAIN PROTEINs exhibits light-hypersensitive phenotypes, suggesting their regulatory roles in modulating optimal light responses. CONCLUSIONS: This study provides an accurate and comprehensive portrait of full-length transcript isoforms and experimentally confirms the presence of de novo synthesized AS variants that impose regulatory functions in photomorphogenic development in Arabidopsis.

Genomically Hardwired Regulation of Gene Activity Orchestrates Cellular Iron Homeostasis in Arabidopsis.

Iron (Fe) is an essential micronutrient which plays pivotal roles as electron donor and catalyst across organisms. In plants, variable, often insufficient Fe supply necessitates mechanisms that constantly attune Fe uptake rates and recalibrate cellular Fe homoeostasis. Here, we show that short-term (0.5, 6, and 12 h) exposure of Arabidopsis thaliana plants to Fe deficiency triggered massive changes in gene activity governed by transcription and alternative splicing (AS), regulatory layers that were to a large extent mutually exclusive. Such preclusion was not observed for genes that are directly involved in the acquisition of Fe, which appears to be concordantly regulated by both expression and AS. Generally, genes with lower splice site strengths and higher intron numbers were more likely to be regulated by AS, no dependence on gene architecture was observed for transcriptionally controlled genes. Conspicuously, specific processes were associated with particular genomic features and biased towards either regulatory mode, suggesting that genomic hardwiring is functionally biased. Early changes in splicing patterns were, in many cases, congruent with later changes in transcript or protein abundance, thus contributing to the pronounced transcriptome-proteome discordance observed in plants.

CFM6 is an Essential CRM Protein Required for the Splicing of nad5 Transcript in Arabidopsis Mitochondria.

Plant chloroplast RNA splicing and ribosome maturation (CRM)-domain-containing proteins are capable of binding RNA to facilitate the splicing of group I or II introns in chloroplasts, but their functions in mitochondria are less clear. In the present study, Arabidopsis thaliana CFM6, a protein with a single CRM domain, was expressed in most plant tissues, particularly in flower tissues, and restricted to mitochondria. Mutation of CFM6 causes severe growth defects, including stunted growth, curled leaves, delayed embryogenesis and pollen development. CFM6 functions specifically in the splicing of group II intron 4 of nad5, which encodes a subunit of mitochondrial complex I, as evidenced by the loss of nad5 intron 4 splicing and high accumulation of its pretranscripts in cfm6 mutants. The phenotypic and splicing defects of cfm6 were rescued in transgenic plants overexpressing 35S::CFM6-YFP. Splicing failure in cfm6 also led to the loss of complex I activity and to its improper assembly. Moreover, dysfunction of complex I induced the expression of proteins or genes involved in alternative respiratory pathways in cfm6. Collectively, CFM6, a previously uncharacterized CRM domain-containing protein, is specifically involved in the cis-splicing of nad5 intron 4 and plays a pivotal role in mitochondrial complex I biogenesis and normal plant growth.

From simple and specific zymographic detections to the annotation of a fungus Daldinia caldariorum D263 that encodes a wide range of highly bioactive cellulolytic enzymes.

BACKGROUND: Lignocellulolytic enzymes are essential for agricultural waste disposal and production of renewable bioenergy. Many commercialized cellulase mixtures have been developed, mostly from saprophytic or endophytic fungal species. The cost of complete cellulose digestion is considerable because a wide range of cellulolytic enzymes is needed. However, most fungi can only produce limited range of highly bioactive cellulolytic enzymes. We aimed to investigate a simple yet specific method for discovering unique enzymes so that fungal species producing a diverse group of cellulolytic enzymes can be identified. RESULTS: The culture medium of an endophytic fungus, Daldinia caldariorum D263, contained a complete set of cellulolytic enzymes capable of effectively digesting cellulose residues into glucose. By taking advantage of the unique product inhibition property of ß-glucosidases, we have established an improved zymography method that can easily distinguish ß-glucosidase and exoglucanase activity. Our zymography method revealed that D263 can secrete a wide range of highly bioactive cellulases. Analyzing the assembled genome of D263, we found over 100 potential genes for cellulolytic enzymes that are distinct from those of the commercially used fungal species Trichoderma reesei and Aspergillus niger. We further identified several of these cellulolytic enzymes by mass spectrometry. CONCLUSIONS: The genome of Daldinia caldariorum D263 has been sequenced and annotated taking advantage of a simple yet specific zymography method followed by mass spectrometry analysis, and it appears to encode and secrete a wide range of bioactive cellulolytic enzymes. The genome and cellulolytic enzyme secretion of this unique endophytic fungus should be of value for identifying active cellulolytic enzymes that can facilitate conversion of agricultural wastes to fermentable sugars for the industrial production of biofuels.

Evidence That Ion-Based Signaling Initiating at the Cell Surface Can Potentially Influence Chromatin Dynamics and Chromatin-Bound Proteins in the Nucleus.

We have developed tools and performed pilot experiments to test the hypothesis that an intracellular ion-based signaling pathway, provoked by an extracellular stimulus acting at the cell surface, can influence interphase chromosome dynamics and chromatin-bound proteins in the nucleus. The experimental system employs chromosome-specific fluorescent tags and the genome-encoded fluorescent pH sensor SEpHluorinA227D, which has been targeted to various intracellular membranes and soluble compartments in root cells of Arabidopsis thaliana. We are using this system and three-dimensional live cell imaging to visualize whether fluorescent-tagged interphase chromosome sites undergo changes in constrained motion concurrently with reductions in membrane-associated pH elicited by extracellular ATP, which is known to trigger a cascade of events in plant cells including changes in calcium ion concentrations, pH, and membrane potential. To examine possible effects of the proposed ion-based signaling pathway directly at the chromatin level, we generated a pH-sensitive fluorescent DNA-binding protein that allows pH changes to be monitored at specific genomic sites. Results obtained using these tools support the existence of a rapid, ion-based signaling pathway that initiates at the cell surface and reaches the nucleus to induce alterations in interphase chromatin mobility and the surrounding pH of chromatin-bound proteins. Such a pathway could conceivably act under natural circumstances to allow external stimuli to swiftly influence gene expression by affecting interphase chromosome movement and the structures and/or activities of chromatin-associated proteins.

Highly ABA-Induced 1 (HAI1)-Interacting protein HIN1 and drought acclimation-enhanced splicing efficiency at intron retention sites.

The Highly ABA-Induced 1 (HAI1) protein phosphatase is a central component of drought-related signaling. A screen for HAI1-interacting proteins identified HAI1-Interactor 1 (HIN1), a nuclear protein of unknown function which could be dephosphorylated by HAI1 in vitro. HIN1 colocalization and interaction with serine-arginine rich (SR) splicing factors and appearance of nuclear speckle-localized HIN1 during low water potential (ψw) stress suggested a pre-mRNA splicing-related function. RNA sequencing of Arabidopsis Col-0 wild type identified more than 500 introns where moderate severity low ψw altered intron retention (IR) frequency. Surprisingly, nearly 90% of these had increased splicing efficiency (decreased IR) during stress. For one-third of these introns, ectopic HIN1 expression (35S:HIN1) in unstressed plants mimicked the increased splicing efficiency seen in stress-treated wild type. HIN1 bound to a GAA-repeat, Exonic Splicing Enhancer-like RNA motif enriched in flanking sequence around HIN1-regulated introns. Genes with stress and HIN1-affected splicing efficiency were enriched for abiotic stress and signaling-related functions. The 35S:HIN1 plants had enhanced growth maintenance during low ψw, while hin1 mutants had reduced growth, further indicating the role of HIN1 in drought response. HIN1 is annotated as an MYB/SANT domain protein but has limited homology to other MYB/SANT proteins and is not related to known yeast or metazoan RNA-binding proteins or splicing regulators. Together these data identify HIN1 as a plant-specific RNA-binding protein, show a specific effect of drought acclimation to promote splicing efficiency of IR-prone introns, and also discover HAI1-HIN1 interaction and dephosphorylation that connects stress signaling to splicing regulation.

Phosphoproteomics of Arabidopsis Highly ABA-Induced1 identifies AT-Hook-Like10 phosphorylation required for stress growth regulation.

The clade A protein phosphatase 2C Highly ABA-Induced 1 (HAI1) plays an important role in stress signaling, yet little information is available on HAI1-regulated phosphoproteins. Quantitative phosphoproteomics identified phosphopeptides of increased abundance in hai1-2 in unstressed plants and in plants exposed to low-water potential (drought) stress. The identity and localization of the phosphoproteins as well as enrichment of specific phosphorylation motifs indicated that these phosphorylation sites may be regulated directly by HAI1 or by HAI1-regulated kinases including mitogen-activated protein kinases, sucrose non-fermenting-related kinase 2, or casein kinases. One of the phosphosites putatively regulated by HAI1 was S313/S314 of AT-Hook-Like10 (AHL10), a DNA-binding protein of unclear function. HAI1 could directly dephosphorylate AHL10 in vitro, and the level of HAI1 expression affected the abundance of phosphorylated AHL10 in vivo. AHL10 S314 phosphorylation was critical for restriction of plant growth under low-water potential stress and for regulation of jasmonic acid and auxin-related gene expression as well as expression of developmental regulators including Shootmeristemless These genes were also misregulated in hai1-2 AHL10 S314 phosphorylation was required for AHL10 complexes to form foci within the nucleoplasm, suggesting that S314 phosphorylation may control AHL10 association with the nuclear matrix or with other transcriptional regulators. These data identify a set of HAI1-affected phosphorylation sites, show that HAI1-regulated phosphorylation of AHL10 S314 controls AHL10 function and localization, and indicate that HAI1-AHL10 signaling coordinates growth with stress and defense responses.

PRP4KA, a Putative Spliceosomal Protein Kinase, Is Important for Alternative Splicing and Development in Arabidopsis thaliana.

Splicing of precursor messenger RNAs (pre-mRNAs) is an essential step in the expression of most eukaryotic genes. Both constitutive splicing and alternative splicing, which produces multiple messenger RNA (mRNA) isoforms from a single primary transcript, are modulated by reversible protein phosphorylation. Although the plant splicing machinery is known to be a target for phosphorylation, the protein kinases involved remain to be fully defined. We report here the identification of pre-mRNA processing 4 (PRP4) KINASE A (PRP4KA) in a forward genetic screen based on an alternatively spliced GFP reporter gene in Arabidopsis thaliana (Arabidopsis). Prp4 kinase is the first spliceosome-associated kinase shown to regulate splicing in fungi and mammals but it has not yet been studied in plants. In the same screen we identified mutants defective in SAC3A, a putative mRNA export factor that is highly coexpressed with PRP4KA in Arabidopsis Whereas the sac3a mutants appear normal, the prp4ka mutants display a pleiotropic phenotype featuring atypical rosettes, late flowering, tall final stature, reduced branching, and lowered seed set. Analysis of RNA-sequencing data from prp4ka and sac3a mutants identified widespread and partially overlapping perturbations in alternative splicing in the two mutants. Quantitative phosphoproteomic profiling of a prp4ka mutant detected phosphorylation changes in several serine/arginine-rich proteins, which regulate constitutive and alternative splicing, and other splicing-related factors. Tests of PRP4KB, the paralog of PRP4KA, indicated that the two genes are not functionally redundant. The results demonstrate the importance of PRP4KA for alternative splicing and plant phenotype, and suggest that PRP4KA may influence alternative splicing patterns by phosphorylating a subset of splicing regulators.

A Genetic Screen Identifies PRP18a, a Putative Second Step Splicing Factor Important for Alternative Splicing and a Normal Phenotype in Arabidopsis thaliana.

Splicing of pre-mRNA involves two consecutive trans-esterification steps that take place in the spliceosome, a large dynamic ribonucleoprotein complex situated in the nucleus. In addition to core spliceosomal proteins, each catalytic step requires step-specific factors. Although the Arabidopsis thaliana genome encodes around 430 predicted splicing factors, functional information about these proteins is limited. In a forward genetic screen based on an alternatively-spliced GFP reporter gene in Arabidopsis thaliana, we identified a mutant impaired in putative step II factor PRP18a, which has not yet been investigated for its role in pre-mRNA splicing in plants. Step II entails cleavage at the 3' splice site accompanied by ligation of the 5' and 3' exons and intron removal. In the prp18 mutant, splicing of a U2-type intron with non-canonical AT-AC splice sites in GFP pre-mRNA is reduced while splicing of a canonical GT-AG intron is enhanced, resulting in decreased levels of translatable GFP mRNA and GFP protein. These findings suggest that wild-type PRP18a may in some cases promote splicing at weak, non-canonical splice sites. Analysis of genome-wide changes in alternative splicing in the prp18a mutant identified numerous cases of intron retention and a preponderance of altered 3' splice sites, suggesting an influence of PRP18a on 3' splice site selection. The prp18a mutant featured short roots on synthetic medium and small siliques, illustrating that wild-type PRP18a function is needed for a normal phenotype. Our study expands knowledge of plant splicing factors and provides foundational information and resources for further functional studies of PRP18 proteins in plants.

A Genetic Screen for Pre-mRNA Splicing Mutants of Arabidopsis thaliana Identifies Putative U1 snRNP Components RBM25 and PRP39a.

In a genetic screen for mutants showing modified splicing of an alternatively spliced GFP reporter gene in Arabidopsis thaliana, we identified mutations in genes encoding the putative U1 small nuclear ribonucleoprotein (snRNP) factors RBM25 and PRP39a. The latter has not yet been studied for its role in pre-messenger RNA (pre-mRNA) splicing in plants. Both proteins contain predicted RNA-binding domains and have been implicated in 5' splice site selection in yeast and metazoan cells. In rbm25 mutants, splicing efficiency of GFP pre-mRNA was reduced and GFP protein levels lowered relative to wild-type plants. By contrast, prp39a mutants exhibited preferential splicing of a U2-type AT-AC intron in GFP pre-mRNA and elevated levels of GFP protein. These opposing findings indicate that impaired function of either RBM25 or PRP39a can differentially affect the same pre-mRNA substrate. Given a prior genome-wide analysis of alternative splicing in rbm25 mutants, we focused on examining the alternative splicing landscape in prp39a mutants. RNA-seq experiments performed using two independent prp39a alleles revealed hundreds of common genes undergoing changes in alternative splicing, including PRP39a itself, a second putative U1 snRNP component PRP40b, and genes encoding a number of general transcription-related proteins. The prp39a mutants displayed somewhat delayed flowering, shorter stature, and reduced seed set but no other obvious common defects under normal conditions. Mutations in PRP39b, the paralog of PRP39a, did not visibly alter GFP expression, indicating the paralogs are not functionally equivalent in this system. Our study provides new information on the contribution of PRP39a to alternative splicing and expands knowledge of plant splicing factors.

A genetic screen implicates a CWC16/Yju2/CCDC130 protein and SMU1 in alternative splicing in Arabidopsis thaliana.

To identify regulators of pre-mRNA splicing in plants, we developed a forward genetic screen based on an alternatively spliced GFP reporter gene in Arabidopsis thaliana In wild-type plants, three major splice variants issue from the GFP gene but only one represents a translatable GFP mRNA. Compared to wild-type seedlings, which exhibit an intermediate level of GFP expression, mutants identified in the screen feature either a "GFP-weak" or "Hyper-GFP" phenotype depending on the ratio of the three splice variants. GFP-weak mutants, including previously identified prp8 and rtf2, contain a higher proportion of unspliced transcript or canonically spliced transcript, neither of which is translatable into GFP protein. In contrast, the coilin-deficient hyper-gfp1 (hgf1) mutant displays a higher proportion of translatable GFP mRNA, which arises from enhanced splicing of a U2-type intron with noncanonical AT-AC splice sites. Here we report three new hgf mutants that are defective, respectively, in spliceosome-associated proteins SMU1, SmF, and CWC16, an Yju2/CCDC130-related protein that has not yet been described in plants. The smu1 and cwc16 mutants have substantially increased levels of translatable GFP transcript owing to preferential splicing of the U2-type AT-AC intron, suggesting that SMU1 and CWC16 influence splice site selection in GFP pre-mRNA. Genome-wide analyses of splicing in smu1 and cwc16 mutants revealed a number of introns that were variably spliced from endogenous pre-mRNAs. These results indicate that SMU1 and CWC16, which are predicted to act directly prior to and during the first catalytic step of splicing, respectively, function more generally to modulate splicing patterns in plants.

Identification of Coilin Mutants in a Screen for Enhanced Expression of an Alternatively Spliced GFP Reporter Gene in Arabidopsis thaliana.

Coilin is a marker protein for subnuclear organelles known as Cajal bodies, which are sites of various RNA metabolic processes including the biogenesis of spliceosomal small nuclear ribonucleoprotein particles. Through self-associations and interactions with other proteins and RNA, coilin provides a structural scaffold for Cajal body formation. However, despite a conspicuous presence in Cajal bodies, most coilin is dispersed in the nucleoplasm and expressed in cell types that lack these organelles. The molecular function of coilin, particularly of the substantial nucleoplasmic fraction, remains uncertain. We identified coilin loss-of-function mutations in a genetic screen for mutants showing either reduced or enhanced expression of an alternatively spliced GFP reporter gene in Arabidopsis thaliana The coilin mutants feature enhanced GFP fluorescence and diminished Cajal bodies compared with wild-type plants. The amount of GFP protein is several-fold higher in the coilin mutants owing to elevated GFP transcript levels and more efficient splicing to produce a translatable GFP mRNA. Genome-wide RNA-sequencing data from two distinct coilin mutants revealed a small, shared subset of differentially expressed genes, many encoding stress-related proteins, and, unexpectedly, a trend toward increased splicing efficiency. These results suggest that coilin attenuates splicing and modulates transcription of a select group of genes. The transcriptional and splicing changes observed in coilin mutants are not accompanied by gross phenotypic abnormalities or dramatically altered stress responses, supporting a role for coilin in fine tuning gene expression. Our GFP reporter gene provides a sensitive monitor of coilin activity that will facilitate further investigations into the functions of this enigmatic protein.

Discriminative gene co-expression network analysis uncovers novel modules involved in the formation of phosphate deficiency-induced root hairs in Arabidopsis.

Cell fate and differentiation in the Arabidopsis root epidermis are genetically defined but remain plastic to environmental signals such as limited availability of inorganic phosphate (Pi). Root hairs of Pi-deficient plants are more frequent and longer than those of plants grown under Pi-replete conditions. To dissect genes involved in Pi deficiency-induced root hair morphogenesis, we constructed a co-expression network of Pi-responsive genes against a customized database that was assembled from experiments in which differentially expressed genes that encode proteins with validated functions in root hair development were over-represented. To further filter out less relevant genes, we combined this procedure with a search for common cis-regulatory elements in the promoters of the selected genes. In addition to well-described players and processes such as auxin signalling and modifications of primary cell walls, we discovered several novel aspects in the biology of root hairs induced by Pi deficiency, including cell cycle control, putative plastid-to-nucleus signalling, pathogen defence, reprogramming of cell wall-related carbohydrate metabolism, and chromatin remodelling. This approach allows the discovery of novel of aspects of a biological process from transcriptional profiles with high sensitivity and accuracy.

Regulation of flowering time by the histone deacetylase HDA5 in Arabidopsis.

The acetylation level of histones on lysine residues regulated by histone acetyltransferases and histone deacetylases plays an important but under-studied role in the control of gene expression in plants. With the aim of characterizing the Arabidopsis RPD3/HDA1 family histone deacetylase HDA5, we present evidence showing that HDA5 displays deacetylase activity. Mutants defective in the expression of HDA5 displayed a late-flowering phenotype. Expression of the flowering repressor genes FLC and MAF1 was up-regulated in hda5 mutants. Furthermore, the gene activation markers, histone H3 acetylation and H3K4 trimethylation on FLC and MAF1 chromatin were increased in hda5-1 mutants. Chromatin immunoprecipitation analysis showed that HDA5 binds to the chromatin of FLC and MAF1. Bimolecular fluorescence complementation assays and co-immunoprecipitation assays showed that HDA5 interacts with FVE, FLD and HDA6, indicating that these proteins are present in a protein complex involved in the regulation of flowering time. Comparing gene expression profiles of hda5 and hda6 mutants by RNA-seq revealed that HDA5 and HDA6 co-regulate gene expression in multiple development processes and pathways.

An Rtf2 Domain-Containing Protein Influences Pre-mRNA Splicing and Is Essential for Embryonic Development in Arabidopsis thaliana.

Alternative splicing is prevalent in plants, but little is known about its regulation in the context of developmental and signaling pathways. We describe here a new factor that influences pre-messengerRNA (mRNA) splicing and is essential for embryonic development in Arabidopsis thaliana. This factor was retrieved in a genetic screen that identified mutants impaired in expression of an alternatively spliced GFP reporter gene. In addition to the known spliceosomal component PRP8, the screen recovered Arabidopsis RTF2 (AtRTF2), a previously uncharacterized, evolutionarily conserved protein containing a replication termination factor 2 (Rtf2) domain. A homozygous null mutation in AtRTF2 is embryo lethal, indicating that AtRTF2 is an essential protein. Quantitative RT-PCR demonstrated that impaired expression of GFP in atrtf2 and prp8 mutants is due to inefficient splicing of the GFP pre-mRNA. A genome-wide analysis using RNA sequencing indicated that 13-16% of total introns are retained to a significant degree in atrtf2 mutants. Considering these results and previous suggestions that Rtf2 represents an ubiquitin-related domain, we discuss the possible role of AtRTF2 in ubiquitin-based regulation of pre-mRNA splicing.

Genome-Wide Analysis of Heat-Sensitive Alternative Splicing in Physcomitrella patens.

Plant growth and development are constantly influenced by temperature fluctuations. To respond to temperature changes, different levels of gene regulation are modulated in the cell. Alternative splicing (AS) is a widespread mechanism increasing transcriptome complexity and proteome diversity. Although genome-wide studies have revealed complex AS patterns in plants, whether AS impacts the stress defense of plants is not known. We used heat shock (HS) treatments at nondamaging temperature and messenger RNA sequencing to obtain HS transcriptomes in the moss Physcomitrella patens. Data analysis identified a significant number of novel AS events in the moss protonema. Nearly 50% of genes are alternatively spliced. Intron retention (IR) is markedly repressed under elevated temperature but alternative donor/acceptor site and exon skipping are mainly induced, indicating differential regulation of AS in response to heat stress. Transcripts undergoing heat-sensitive IR are mostly involved in specific functions, which suggests that plants regulate AS with transcript specificity under elevated temperature. An exonic GAG-repeat motif in these IR regions may function as a regulatory cis-element in heat-mediated AS regulation. A conserved AS pattern for HS transcription factors in P. patens and Arabidopsis (Arabidopsis thaliana) reveals that heat regulation for AS evolved early during land colonization of green plants. Our results support that AS of specific genes, including key HS regulators, is fine-tuned under elevated temperature to modulate gene regulation and reorganize metabolic processes.