SWAP1-SFPS-RRC1 splicing factor complex modulates pre-mRNA splicing to promote photomorphogenesis in Arabidopsis.

Light signals perceived by a group of photoreceptors have profound effects on the physiology, growth, and development of plants. The red/far-red light-absorbing phytochromes (phys) modulate these aspects by intricately regulating gene expression at multiple levels. Here, we report the identification and functional characterization of an RNA-binding splicing factor, SWAP1 (SUPPRESSOR-OF-WHITE-APRICOT/SURP RNA-BINDING DOMAIN-CONTAINING PROTEIN1). Loss-of-function swap1-1 mutant is hyposensitive to red light and exhibits a day length-independent early flowering phenotype. SWAP1 physically interacts with two other splicing factors, (SFPS) SPLICING FACTOR FOR PHYTOCHROME SIGNALING and (RRC1) REDUCED RED LIGHT RESPONSES IN CRY1CRY2 BACKGROUND 1 in a light-independent manner and forms a ternary complex. In addition, SWAP1 physically interacts with photoactivated phyB and colocalizes with nuclear phyB photobodies. Phenotypic analyses show that the swap1sfps, swap1rrc1, and sfpsrrc1 double mutants display hypocotyl lengths similar to that of the respective single mutants under red light, suggesting that they function in the same genetic pathway. The swap1sfps double and swap1sfpsrrc1 triple mutants display pleiotropic phenotypes, including sterility at the adult stage. Deep RNA sequencing (RNA-seq) analyses show that SWAP1 regulates the gene expression and pre-messenger RNA (mRNA) alternative splicing of a large number of genes, including those involved in plant responses to light signaling. A comparative analysis of alternative splicing among single, double, and triple mutants showed that all three splicing factors coordinately regulate the alternative splicing of a subset of genes. Our study uncovered the function of a splicing factor that modulates light-regulated alternative splicing by interacting with photoactivated phyB and other splicing factors.

An in-frame deletion mutation in the degron tail of auxin coreceptor IAA2 confers resistance to the herbicide 2,4-D in Sisymbrium orientale.

The natural auxin indole-3-acetic acid (IAA) is a key regulator of many aspects of plant growth and development. Synthetic auxin herbicides such as 2,4-D mimic the effects of IAA by inducing strong auxinic-signaling responses in plants. To determine the mechanism of 2,4-D resistance in a Sisymbrium orientale (Indian hedge mustard) weed population, we performed a transcriptome analysis of 2,4-D-resistant (R) and -susceptible (S) genotypes that revealed an in-frame 27-nucleotide deletion removing nine amino acids in the degron tail (DT) of the auxin coreceptor Aux/IAA2 (SoIAA2). The deletion allele cosegregated with 2,4-D resistance in recombinant inbred lines. Further, this deletion was also detected in several 2,4-D-resistant field populations of this species. Arabidopsis transgenic lines expressing the SoIAA2 mutant allele were resistant to 2,4-D and dicamba. The IAA2-DT deletion reduced binding to TIR1 in vitro with both natural and synthetic auxins, causing reduced association and increased dissociation rates. This mechanism of synthetic auxin herbicide resistance assigns an in planta function to the DT region of this Aux/IAA coreceptor for its role in synthetic auxin binding kinetics and reveals a potential biotechnological approach to produce synthetic auxin-resistant crops using gene-editing.

Drought induces epitranscriptome and proteome changes in stem-differentiating xylem of Populus trichocarpa.

Understanding gene expression and regulation requires insights into RNA transcription, processing, modification, and translation. However, the relationship between the epitranscriptome and the proteome under drought stress remains undetermined in poplar (Populus trichocarpa). In this study, we used Nanopore direct RNA sequencing and tandem mass tag-based proteomic analysis to examine epitranscriptomic and proteomic regulation induced by drought treatment in stem-differentiating xylem (SDX). Our results revealed a decreased full-length read ratio under drought treatment and, especially, a decreased association between transcriptome and proteome changes in response to drought. Epitranscriptome analysis of cellulose- and lignin-related genes revealed an increased N6-Methyladenosine (m6A) ratio, which was accompanied by decreased RNA abundance and translation, under drought stress. Interestingly, usage of the distal poly(A) site increased during drought stress. Finally, we found that transcripts of highly expressed genes tend to have shorter poly(A) tail length (PAL), and drought stress increased the percentage of transcripts with long PAL. These findings provide insights into the interplay among m6A, polyadenylation, PAL, and translation under drought stress in P. trichocarpa SDX.

Alternative Splicing Variation: Accessing and Exploiting in Crop Improvement Programs.

Alternative splicing (AS) is a gene regulatory mechanism modulating gene expression in multiple ways. AS is prevalent in all eukaryotes including plants. AS generates two or more mRNAs from the precursor mRNA (pre-mRNA) to regulate transcriptome complexity and proteome diversity. Advances in next-generation sequencing, omics technology, bioinformatics tools, and computational methods provide new opportunities to quantify and visualize AS-based quantitative trait variation associated with plant growth, development, reproduction, and stress tolerance. Domestication, polyploidization, and environmental perturbation may evolve novel splicing variants associated with agronomically beneficial traits. To date, pre-mRNAs from many genes are spliced into multiple transcripts that cause phenotypic variation for complex traits, both in model plant Arabidopsis and field crops. Cataloguing and exploiting such variation may provide new paths to enhance climate resilience, resource-use efficiency, productivity, and nutritional quality of staple food crops. This review provides insights into AS variation alongside a gene expression analysis to select for novel phenotypic diversity for use in breeding programs. AS contributes to heterosis, enhances plant symbiosis (mycorrhiza and rhizobium), and provides a mechanistic link between the core clock genes and diverse environmental clues.

The amino acid region from 448-517 of CAMTA3 transcription factor containing a part of the TIG domain represses the N-terminal repression module function.

CAMTA3, a Ca2+-regulated transcription factor, is a repressor of plant immune responses. A truncated version of CAMTA3; CAMTA3334 called N-terminal repression module (NRM), and its extended version (CAMTA447), which include the DNA binding domain, were previously reported to complement the camta3/2 mutant phenotype. Here, we generated a series of CAMTA3 truncated versions [the N-terminus (aa 1-517), C-terminus (aa 517-1032), R1 (aa 1-173), R2 (aa 174-345), R3 (aa 346-517), R4 (aa 517-689), R5 (aa 690-861) and R6 (aa 862-1032)], expressed in camta3 mutant and analyzed the phenotypes of the transgenic lines. Interestingly, unlike CAMTA447, extending the N-terminal region to 517 aa did not complement the camta3 phenotype, suggesting that the amino acid region from 448-517 (70 aa), which includes a part of the TIG domain suppresses the NRM activity. The C-terminus and other truncated versions (R1-R6) also failed to complement the camta3 mutant. Expressing the full length or NRM of CAMTA3 in camta3 plants suppressed the activation of immune-responsive genes and increased the expression of cold-induced genes. In contrast, the transgenic lines expressing the N- or C-terminus or R1-R6 of CAMTA3 showed expression patterns like those of the camta3 with enhanced expression of the defense genes and suppressed expression of the cold response genes. Furthermore, like camta3, the transgenic lines expressing the N- or C-terminus, or R1-R6 of CAMTA3 exhibited higher levels of H2O2 and increased resistance to a Pst DC3000 as compared to WT, NRM, or FL-CAMTA3 transgenic plants. Our studies identified a novel regulatory region in CAMTA3 that suppresses the NRM activity. Supplementary Information: The online version contains supplementary material available at 10.1007/s12298-023-01401-w.

RODAN: a fully convolutional architecture for basecalling nanopore RNA sequencing data.

BACKGROUND: Despite recent progress in basecalling of Oxford nanopore DNA sequencing data, its wide adoption is still being hampered by its relatively low accuracy compared to short read technologies. Furthermore, very little of the recent research was focused on basecalling of RNA data, which has different characteristics than its DNA counterpart. RESULTS: We fill this gap by benchmarking a fully convolutional deep learning basecalling architecture with improved performance compared to Oxford nanopore's RNA basecallers. AVAILABILITY: The source code for our basecaller is available at: https://github.com/biodlab/RODAN .

Wide-ranging transcriptome remodelling mediated by alternative polyadenylation in response to abiotic stresses in Sorghum.

Alternative polyadenylation (APA) regulates diverse developmental and physiological processes through its effects on gene expression, mRNA stability, translatability, and transport. Sorghum is a major cereal crop in the world and, despite its importance, not much is known about the role of post-transcriptional regulation in mediating responses to abiotic stresses in Sorghum. A genome-wide APA analysis unveiled widespread occurrence of APA in Sorghum in response to drought, heat, and salt stress. Abiotic stress treatments incited changes in poly(A) site choice in a large number of genes. Interestingly, abiotic stresses led to the re-directing of transcriptional output into non-productive pathways defined by the class of poly(A) site utilized. This result revealed APA to be part of a larger global response of Sorghum to abiotic stresses that involves the re-direction of transcriptional output into non-productive transcriptional and translational pathways. Large numbers of stress-inducible poly(A) sites could not be linked with known, annotated genes, suggestive of the existence of numerous unidentified genes whose expression is strongly regulated by abiotic stresses. Furthermore, we uncovered a novel stress-specific cis-element in intronic poly(A) sites used in drought- and heat-stressed plants that might play an important role in non-canonical poly(A) site choice in response to abiotic stresses.

Alternative splicing dynamics and evolutionary divergence during embryogenesis in wheat species.

Among polyploid species with complex genomic architecture, variations in the regulation of alternative splicing (AS) provide opportunities for transcriptional and proteomic plasticity and the potential for generating trait diversities. However, the evolution of AS and its influence on grain development in diploid grass and valuable polyploid wheat crops are poorly understood. To address this knowledge gap, we developed a pipeline for the analysis of alternatively spliced transcript isoforms, which takes the high sequence similarity among polyploid wheat subgenomes into account. Through analysis of synteny and detection of collinearity of homoeologous subgenomes, conserved and specific AS events across five wheat and grass species were identified. A global analysis of the regulation of AS in diploid grass and polyploid wheat grains revealed diversity in AS events not only between the endosperm, pericarp and embryo overdevelopment, but also between subgenomes. Analysis of AS in homoeologous triads of polyploid wheats revealed evolutionary divergence between gene-level and transcript-level regulation of embryogenesis. Evolutionary age analysis indicated that the generation of novel transcript isoforms has occurred in young genes at a more rapid rate than in ancient genes. These findings, together with the development of comprehensive AS resources for wheat and grass species, advance understanding of the evolution of regulatory features of AS during embryogenesis and grain development in wheat.

Differential nucleosome occupancy modulates alternative splicing in Arabidopsis thaliana.

Alternative splicing (AS) is a major gene regulatory mechanism in plants. Recent evidence supports co-transcriptional splicing in plants, hence the chromatin state can impact AS. However, how dynamic changes in the chromatin state such as nucleosome occupancy influence the cold-induced AS remains poorly understood. Here, we generated transcriptome (RNA-Seq) and nucleosome positioning (MNase-Seq) data for Arabidopsis thaliana to understand how nucleosome positioning modulates cold-induced AS. Our results show that characteristic nucleosome occupancy levels are strongly associated with the type and abundance of various AS events under normal and cold temperature conditions in Arabidopsis. Intriguingly, exitrons, alternatively spliced internal regions of protein-coding exons, exhibit distinctive nucleosome positioning pattern compared to other alternatively spliced regions. Likewise, nucleosome patterns differ between exitrons and retained introns, pointing to their distinct regulation. Collectively, our data show that characteristic changes in nucleosome positioning modulate AS in plants in response to cold.

Does co-transcriptional regulation of alternative splicing mediate plant stress responses?

Plants display exquisite control over gene expression to elicit appropriate responses under normal and stress conditions. Alternative splicing (AS) of pre-mRNAs, a process that generates two or more transcripts from multi-exon genes, adds another layer of regulation to fine-tune condition-specific gene expression in animals and plants. However, exactly how plants control splice isoform ratios and the timing of this regulation in response to environmental signals remains elusive. In mammals, recent evidence indicate that epigenetic and epitranscriptome changes, such as DNA methylation, chromatin modifications and RNA methylation, regulate RNA polymerase II processivity, co-transcriptional splicing, and stability and translation efficiency of splice isoforms. In plants, the role of epigenetic modifications in regulating transcription rate and mRNA abundance under stress is beginning to emerge. However, the mechanisms by which epigenetic and epitranscriptomic modifications regulate AS and translation efficiency require further research. Dynamic changes in the chromatin landscape in response to stress may provide a scaffold around which gene expression, AS and translation are orchestrated. Finally, we discuss CRISPR/Cas-based strategies for engineering chromatin architecture to manipulate AS patterns (or splice isoforms levels) to obtain insight into the epigenetic regulation of AS.

The Landscape of RNA-Protein Interactions in Plants: Approaches and Current Status.

RNAs transmit information from DNA to encode proteins that perform all cellular processes and regulate gene expression in multiple ways. From the time of synthesis to degradation, RNA molecules are associated with proteins called RNA-binding proteins (RBPs). The RBPs play diverse roles in many aspects of gene expression including pre-mRNA processing and post-transcriptional and translational regulation. In the last decade, the application of modern techniques to identify RNA-protein interactions with individual proteins, RNAs, and the whole transcriptome has led to the discovery of a hidden landscape of these interactions in plants. Global approaches such as RNA interactome capture (RIC) to identify proteins that bind protein-coding transcripts have led to the identification of close to 2000 putative RBPs in plants. Interestingly, many of these were found to be metabolic enzymes with no known canonical RNA-binding domains. Here, we review the methods used to analyze RNA-protein interactions in plants thus far and highlight the understanding of plant RNA-protein interactions these techniques have provided us. We also review some recent protein-centric, RNA-centric, and global approaches developed with non-plant systems and discuss their potential application to plants. We also provide an overview of results from classical studies of RNA-protein interaction in plants and discuss the significance of the increasingly evident ubiquity of RNA-protein interactions for the study of gene regulation and RNA biology in plants.

Genome-wide profiling of circular RNAs, alternative splicing, and R-loops in stem-differentiating xylem of Populus trichocarpa.

Circular RNAs (circRNAs) are a recently discovered type of non-coding RNA derived from pre-mRNAs. R-loops consist of a DNA:RNA hybrid and the associated single-stranded DNA. In Arabidopsis thaliana, circRNA:DNA R-loops regulate alternative splicing (AS) of SEPALLATA3 (SEP3). However, the occurrence and functions of circRNAs and R-loops in Populus trichocarpa are largely unexplored. Here, we performed circRNA-enriched sequencing in the stem-differentiating xylem (SDX) of P. trichocarpa and identified 2,742 distinct circRNAs, including circ-CESA4, circ-IRX7, and circ-GUX1, which are generated from genes involved in cellulose, and hemicellulose biosynthesis, respectively. To investigate the roles of circRNAs in modulating alternative splicing (AS), we detected 7,836 AS events using PacBio Iso-Seq and identified 634 circRNAs that overlapped with 699 AS events. Furthermore, using DNA:RNA hybrid immunoprecipitation followed by sequencing (DRIP-seq), we identified 8,932 R-loop peaks that overlapped with 181 circRNAs and 672 AS events. Notably, several SDX-related circRNAs overlapped with R-loop peaks, pointing to their possible roles in modulating AS in SDX. Indeed, overexpressing circ-IRX7 increased the levels of R-loop structures and decreased the frequency of intron retention in linear IRX7 transcripts. This study provides a valuable R-loop atlas resource and uncovers the interplay between circRNAs and AS in SDX of P. trichocarpa.

Salt-Induced Stability of SR1/CAMTA3 mRNA Is Mediated by Reactive Oxygen Species and Requires the 3' End of Its Open Reading Frame.

Soil salinity, a prevalent abiotic stress, causes enormous losses in global crop yields annually. Previous studies have shown that salt stress-induced reprogramming of gene expression contributes to the survival of plants under this stress. However, mechanisms regulating gene expression in response to salt stress at the posttranscriptional level are not well understood. In this study, we show that salt stress increases the level of Signal Responsive 1 (SR1) mRNA, a member of signal-responsive Ca2+/calmodulin-regulated transcription factors, by enhancing its stability. We present multiple lines of evidence indicating that reactive oxygen species generated by NADPH oxidase activity mediate salt-induced SR1 transcript stability. Using mutants impaired in either nonsense-mediated decay, XRN4 or mRNA decapping pathways, we show that neither the nonsense-mediated mRNA decay pathway, XRN4 nor the decapping of SR1 mRNA is required for its decay. We analyzed the salt-induced accumulation of eight truncated versions of the SR1 coding region (∼3 kb) in the sr1 mutant background. This analysis identified a 500-nt region at the 3' end of the SR1 coding region to be required for the salt-induced stability of SR1 mRNA. Potential mechanisms by which this region confers SR1 transcript stability in response to salt are discussed.

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.

Decoding co-/post-transcriptional complexities of plant transcriptomes and epitranscriptome using next-generation sequencing technologies.

Next-generation sequencing (NGS) technologies - Illumina RNA-seq, Pacific Biosciences isoform sequencing (PacBio Iso-seq), and Oxford Nanopore direct RNA sequencing (DRS) - have revealed the complexity of plant transcriptomes and their regulation at the co-/post-transcriptional level. Global analysis of mature mRNAs, transcripts from nuclear run-on assays, and nascent chromatin-bound mRNAs using short as well as full-length and single-molecule DRS reads have uncovered potential roles of different forms of RNA polymerase II during the transcription process, and the extent of co-transcriptional pre-mRNA splicing and polyadenylation. These tools have also allowed mapping of transcriptome-wide start sites in cap-containing RNAs, poly(A) site choice, poly(A) tail length, and RNA base modifications. The emerging theme from recent studies is that reprogramming of gene expression in response to developmental cues and stresses at the co-/post-transcriptional level likely plays a crucial role in eliciting appropriate responses for optimal growth and plant survival under adverse conditions. Although the mechanisms by which developmental cues and different stresses regulate co-/post-transcriptional splicing are largely unknown, a few recent studies indicate that the external cues target spliceosomal and splicing regulatory proteins to modulate alternative splicing. In this review, we provide an overview of recent discoveries on the dynamics and complexities of plant transcriptomes, mechanistic insights into splicing regulation, and discuss critical gaps in co-/post-transcriptional research that need to be addressed using diverse genomic and biochemical approaches.

Transcriptome Analysis of Drought-Resistant and Drought-Sensitive Sorghum (Sorghum bicolor) Genotypes in Response to PEG-Induced Drought Stress.

Drought is a major limiting factor of crop yields. In response to drought, plants reprogram their gene expression, which ultimately regulates a multitude of biochemical and physiological processes. The timing of this reprogramming and the nature of the drought-regulated genes in different genotypes are thought to confer differential tolerance to drought stress. Sorghum is a highly drought-tolerant crop and has been increasingly used as a model cereal to identify genes that confer tolerance. Also, there is considerable natural variation in resistance to drought in different sorghum genotypes. Here, we evaluated drought resistance in four genotypes to polyethylene glycol (PEG)-induced drought stress at the seedling stage and performed transcriptome analysis in seedlings of sorghum genotypes that are either drought-resistant or drought-sensitive to identify drought-regulated changes in gene expression that are unique to drought-resistant genotypes of sorghum. Our analysis revealed that about 180 genes are differentially regulated in response to drought stress only in drought-resistant genotypes and most of these (over 70%) are up-regulated in response to drought. Among these, about 70 genes are novel with no known function and the remaining are transcription factors, signaling and stress-related proteins implicated in drought tolerance in other crops. This study revealed a set of drought-regulated genes, including many genes encoding uncharacterized proteins that are associated with drought tolerance at the seedling stage.

Profiling of circular RNA N6 -methyladenosine in moso bamboo (Phyllostachys edulis) using nanopore-based direct RNA sequencing.

N6 -methyladenosine (m6 A) is a prevalent modification in messenger RNAs and circular RNAs that play important roles in regulating various aspects of RNA metabolism. However, the occurrence of the m6 A modification in plant circular RNAs has not been reported. A widely used method to identify m6 A modifications relies on m6 A-specific antibodies followed by next-generation sequencing of precipitated RNAs (MeRIP-Seq). However, one limitation of MeRIP-Seq is that it does not provide the precise location of m6 A at single-nucleotide resolution. Although more recent sequencing techniques such as Nanopore-based direct RNA sequencing (DRS) can overcome such limitations, the technology does not allow sequencing of circular RNAs, as these molecules lack a poly(A) tail. Here, we developed a novel method to detect the precise location of m6 A modifications in circular RNAs using Nanopore DRS. We first enriched our samples for circular RNAs, which we then fragmented and sequenced on the Nanopore platform with a customized protocol. Using this method, we identified 470 unique circular RNAs from DRS reads based on the back-spliced junction region. Among exonic circular RNAs, about 10% contained m6 A sites, which mainly occurred around acceptor and donor splice sites. This study demonstrates the utility of our antibody-independent method in identifying total and methylated circular RNAs using Nanopore DRS. This method has the additional advantage of providing the exact location of m6 A sites at single-base resolution in circular RNAs or linear transcripts from non-coding RNA without poly(A) tails.

The Arabidopsis splicing regulator SR45 confers salt tolerance in a splice isoform-dependent manner.

KEY MESSAGE: Functions of most splice isoforms that are generated by alternative splicing are unknown. We show that two splice variants that encode proteins differing in only eight amino acids have distinct functions in a stress response. Serine/arginine-rich (SR) and SR-like proteins, a conserved family of RNA binding proteins across eukaryotes, play important roles in pre-mRNA splicing and other post-transcriptional processes. Pre-mRNAs of SR and SR-like proteins undergo extensive alternative splicing in response to diverse stresses and produce multiple splice isoforms. However, the functions of most splice isoforms remain elusive. Alternative splicing of pre-mRNA of Arabidopsis SR45, which encodes an SR-like splicing regulator, generates two isoforms (long-SR45.1 and short-SR45.2). The proteins encoded by these two isoforms differ in eight amino acids. Here, we investigated the role of SR45 and its splice variants in salt stress tolerance. The loss of SR45 resulted in enhanced sensitivity to salt stress and changes in expression and splicing of genes involved in regulating salt stress response. Interestingly, only the long isoform (SR45.1) rescued the salt-sensitive phenotype as well as the altered gene expression and splicing patterns in the mutant. These results suggest that SR45 positively regulates salt tolerance. Furthermore, only the long isoform is required for SR45-mediated salt tolerance.

Exploring the relationship between intron retention and chromatin accessibility in plants.

BACKGROUND: Intron retention (IR) is the most prevalent form of alternative splicing in plants. IR, like other forms of alternative splicing, has an important role in increasing gene product diversity and regulating transcript functionality. Splicing is known to occur co-transcriptionally and is influenced by the speed of transcription which in turn, is affected by chromatin structure. It follows that chromatin structure may have an important role in the regulation of splicing, and there is preliminary evidence in metazoans to suggest that this is indeed the case; however, nothing is known about the role of chromatin structure in regulating IR in plants. DNase I-seq is a useful experimental tool for genome-wide interrogation of chromatin accessibility, providing information on regions of chromatin with very high likelihood of cleavage by the enzyme DNase I, known as DNase I Hypersensitive Sites (DHSs). While it is well-established that promoter regions are highly accessible and are over-represented with DHSs, not much is known about DHSs in the bodies of genes, and their relationship to splicing in general, and IR in particular. RESULTS: In this study we use publicly available DNase I-seq data in arabidopsis and rice to investigate the relationship between IR and chromatin structure. We find that IR events are highly enriched in DHSs in both species. This implies that chromatin is more open in retained introns, which is consistent with a kinetic model of the process whereby higher speeds of transcription in those regions give less time for the spliceosomal machinery to recognize and splice out those introns co-transcriptionally. The more open chromatin in IR can also be the result of regulation mediated by DNA-binding proteins. To test this, we performed an exhaustive search for footprints left by DNA-binding proteins that are associated with IR. We identified several hundred short sequence elements that exhibit footprints in their DNase I-seq coverage, the telltale sign for binding events of a regulatory protein, protecting its binding site from cleavage by DNase I. A highly significant fraction of those sequence elements are conserved between arabidopsis and rice, a strong indication of their functional importance. CONCLUSIONS: In this study we have established an association between IR and chromatin accessibility, and presented a mechanistic hypothesis that explains the observed association from the perspective of the co-transcriptional nature of splicing. Furthermore, we identified conserved sequence elements for DNA-binding proteins that affect splicing.