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

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

Catsnap: a user-friendly algorithm for determining the conservation of protein variants reveals extensive parallelisms in the evolution of alternative splicing

Understanding the evolutionary conservation of complex eukaryotic transcriptomes significantly illuminates the physiological relevance of alternative splicing (AS). Examining the evolutionary depth of a given AS event with ordinary homology searches is generally challenging and time-consuming. Here, we present Catsnap, an algorithmic pipeline for assessing the conservation of putative protein isoforms generated by AS. It employs a machine learning approach following a database search with the provided pair of protein sequences. We used the Catsnap algorithm for analyzing the conservation of emerging experimentally characterized alternative proteins from plants and animals. Indeed, most of them are conserved among other species. Catsnap can detect the conserved functional protein isoforms regardless of the AS type by which they are generated. Notably, we found that while the primary amino acid sequence is maintained, the type of AS determining the inclusion or exclusion of protein regions varies throughout plant phylogenetic lineages in these proteins. We also document that this phenomenon is less seen among animals. In sum, our algorithm highlights the presence of unexpectedly frequent hotspots where protein isoforms recurrently arise to carry physiologically relevant functions. The user web interface is available at https://catsnap.cesnet.cz/.

Alternative splicing: transcriptional regulatory network in agroforestry

Alternative splicing (AS) in plants plays a key role in regulating the expression of numerous transcripts from a single gene in a regulatory pathway. Variable concentrations of growth regulatory hormones and external stimuli trigger alternative splicing to switch among different growth stages and adapt to environmental stresses. In the AS phenomenon, a spliceosome causes differential transcriptional modifications in messenger RNA (mRNAs), resulting in partial or complete retention of one or more introns as compared to fully spliced mRNA. Differentially expressed proteins translated from intron-retaining messenger RNA (mRNA^ir) perform vital functions in the feedback mechanism. At the post-transcriptional level, AS causes the remodeling of transcription factors (TFs) by the addition or deletion of binding domains to activate and/or repress transcription. In this study, we have summarized the specific role of AS in the regulation of gene expression through repression and activation of the transcriptional regulatory network under external stimuli and switch among developmental stages.

Integrated omics analysis identified genes and their splice variants involved in fruit development and metabolites production in Capsicum species

To date, several transcriptomic studies during fruit development have been reported; however, no comprehensive integrated study on expression diversity, alternative splicing, and metabolomic profiling was reported in Capsicum. This study analyzed RNA-seq data and untargeted metabolomic profiling from early green (EG), mature green (MG), and breaker (Br) fruit stages from two Capsicum species, i.e., C. annuum (Cann) and C. frutescens (Cfrut) from Northeast India. A total of 117,416 and 96,802 alternatively spliced events (AltSpli-events) were identified from Cann and Cfrut, respectively. Among AltSpli-events, intron retention (IR; 32.2% Cann and 25.75% Cfrut) followed by alternative acceptor (AA; 15.4% Cann and 18.9% Cfrut) were the most abundant in Capsicum. Around 7600 genes expressed in at least one fruit stage of Cann and Cfrut were AltSpli. The study identified spliced variants of genes including transcription factors (TFs) potentially involved in fruit development/ripening (Aux/IAA 16-like, ETR, SGR1, ARF, CaGLK2, ETR, CaAGL1, MADS-RIN, FUL1, SEPALLATA1), carotenoid (PDS, CA1, CCD4, NCED3, xanthoxin dehydrogenase, CaERF82, CabHLH100, CaMYB3R-1, SGR1, CaWRKY28, CaWRKY48, CaWRKY54), and capsaicinoids or flavonoid biosynthesis (CaMYB48, CaWRKY51), which were significantly differentially spliced (DS) between consecutive Capsicum fruit stages. Also, this study observed that differentially expressed isoforms (DEiso) from 38 genes with differentially spliced events (DSE) were significantly enriched in various metabolic pathways such as starch and sucrose metabolism, amino acid metabolism, cysteine cutin suberin and wax biosynthesis, and carotenoid biosynthesis. Furthermore, the metabolomic profiling revealed that metabolites from aforementioned pathways such as carbohydrates (mainly sugars such as D-fructose, D-galactose, maltose, and sucrose), organic acids (carboxylic acids), and peptide groups significantly altered during fruit development. Taken together, our findings could help in alternative splicing-based targeted studies of candidate genes involved in fruit development and ripening in Capsicum crop.

Xylem Transcriptome Analysis in Contrasting Wood Phenotypes of Eucalyptus urophylla × tereticornis Hybrids

An investigation of the effects of two important post-transcriptional regulatory mechanisms, gene transcription and alternative splicing (AS), on the wood formation of Eucalyptusurophylla × tereticornis, an economic tree species widely planted in southern China, was carried out. We performed RNA-seq on E. urophylla × tereticornis hybrids with highly contrasting wood basic density (BD), cellulose content (CC), hemicellulose content (HC), and lignin content (LC). Signals of strong differentially expressed genes (DEGs) and differentially spliced genes (DSGs) were detected in all four groups of wood properties, suggesting that gene transcription and selective splicing may have important regulatory roles in wood properties. We found that there was little overlap between DEGs and DSGs in groups of the same trait. Furthermore, the key DEGs and DSGs that were detected simultaneously in the four groups tended to be enriched in different Gene Ontology terms, Kyoto Encyclopedia of Genes and Genomes pathways, and transcription factors. These results implied that regulation of gene transcription and AS is controlled by independent regulatory systems in wood formation. Lastly, we detected transcript levels of known wood biosynthetic genes and found that 79 genes encoding mainly enzymes or proteins such as UGT, LAC, CAD, and CESA may be involved in the positive or negative regulation of wood properties. This study reveals potential molecular mechanisms that may regulate wood formation and will contribute to the genetic improvement of Eucalyptus.

Alternative splicing during fruit development among fleshy fruits

BACKGROUND

Alternative splicing (AS) is an important mechanism of posttranscriptional modification and dynamically regulates multiple physiological processes in plants, including fruit ripening. However, little is known about alternative splicing during fruit development in fleshy fruits.

RESULTS

We studied the alternative splicing at the immature and ripe stages during fruit development in cucumber, melon, papaya and peach. We found that 14.96-17.48% of multiexon genes exhibited alternative splicing. Intron retention was not always the most frequent event, indicating that the alternative splicing pattern during different developmental process differs. Alternative splicing was significantly more prevalent at the ripe stage than at the immature stage in cucumber and melon, while the opposite trend was shown in papaya and peach, implying that developmental stages adopt different alternative splicing strategies for their specific functions. Some genes involved in fruit ripening underwent stage-specific alternative splicing, indicating that alternative splicing regulates fruits ripening. Conserved alternative splicing events did not appear to be stage-specific. Clustering fruit developmental stages across the four species based on alternative splicing profiles resulted in species-specific clustering, suggesting that diversification of alternative splicing contributes to lineage-specific evolution in fleshy fruits.

CONCLUSIONS

We obtained high quality transcriptomes and alternative splicing events during fruit development across the four species. Dynamics and nonconserved alternative splicing were discovered. The candidate stage-specific AS genes involved in fruit ripening will provide valuable insight into the roles of alternative splicing during the developmental processes of fleshy fruits.

Transcriptome-wide identification and characterization of microRNAs in diverse phases of wood formation in Populus trichocarpa

We applied miRNA expression profiling method to Populus trichocarpa stems of the three developmental stages, primary stem (PS), transitional stem (TS), and secondary stem (SS), to investigate miRNA species and their regulation on lignocellulosic synthesis and related processes. We obtained 892, 872, and 882 known miRNAs and 1727, 1723, and 1597 novel miRNAs, from PS, TS, and SS, respectively. Comparisons of these miRNA species among different developmental stages led to the identification of 114, 306, and 152 differentially expressed miRNAs (DE-miRNAs), which had 921, 2639, and 2042 candidate target genes (CTGs) in the three respective stages of the same order. Correlation analysis revealed 47, 439, and 71 DE-miRNA-CTG pairs of high negative correlation in PS, TS, and SS, respectively. Through biological process analysis, we finally identified 34, 6, and 76 miRNA-CTG pairs from PS, TS, and SS, respectively, and the miRNA target genes in these pairs regulate or participate lignocellulosic biosynthesis-related biological processes: cell division and differentiation, cell wall modification, secondary cell wall biosynthesis, lignification, and programmed cell death processes. This is the first report on an integrated analysis of genome-wide mRNA and miRNA profilings during multiple phases of poplar stem development. Our analysis results imply that individual miRNAs modulate secondary growth and lignocellulosic biosynthesis through regulating transcription factors and lignocellulosic biosynthetic pathway genes, resulting in more dynamic promotion, suppression, or regulatory circuits. This study advanced our understanding of many individual miRNAs and their essential, diversified roles in the dynamic regulation of secondary growth in woody tree species.

Investigation Into Different Wood Formation Mechanisms Between Angiosperm and Gymnosperm Tree Species at the Transcriptional and Post-transcriptional Level

Enormous distinctions of the stem structure and cell types between gymnosperms and angiosperms tree species are expected to cause quite different wood physical and mechanical attributes, however, the molecular mechanisms underlying the differing wood morphology are still unclear. In this study, we compared the transcriptomes obtained by RNA-Seq between Populus alba × P. glandulosa clone 84K, and Larix kaempferi (Lamb.) Carr trees. Available genome resource served as reference for P. alba × P. glandulosa and the Iso-Seq results of a three-tissues mixture (xylem, phloem, and leaf) were used as the reference for L. kaempferi to compare the xylem-specifically expressed genes and their alternative splicing model. Through screening, we obtained 13,907 xylem-specifically expressed genes (5,954 up-regulated, 7,953 down-regulated) in the xylem of P. alba × P. glandulosa, and 2,596 xylem-specifically expressed genes (1,648 up-regulated, 948 down-regulated) in the xylem of L. kaempferi. From the GO and KEGG analyses, some genes associated with two wood formation-related pathways, namely those for phenylpropanoid biosynthesis, and starch and sucrose metabolism, were successfully screened. Then the distributions and gene expression models between P. alba × P. glandulosa and L. kaempferi in those pathways were compared, which suggested differential wood formation processes between the angiosperm and gymnosperm trees. Furthermore, a Weight Gene Co-expression Network Analysis (WGCNA) for total xylem-specifically expressed genes in two species was conducted, from which wood formation-related modules were selected to build a co-expression network for the two tree species. The genes within this co-expression network showed different co-expression relationships between the angiosperm and gymnosperm woody species. Comparing the alternative splicing events for wood formation-related genes suggests a different post-transcriptional regulation process exists between the angiosperm and gymnosperm trees. Our research thus provides the foundation for the in-depth investigation of different wood formation mechanisms of angiosperm and gymnosperm species.

Breeding for Climate Change Resilience: A Case Study of Loblolly Pine (Pinus taeda L.) in North America

Earth's atmosphere is warming and the effects of climate change are becoming evident. A key observation is that both the average levels and the variability of temperature and precipitation are changing. Information and data from new technologies are developing in parallel to provide multidisciplinary opportunities to address and overcome the consequences of these changes in forest ecosystems. Changes in temperature and water availability impose multidimensional environmental constraints that trigger changes from the molecular to the forest stand level. These can represent a threat for the normal development of the tree from early seedling recruitment to adulthood both through direct mortality, and by increasing susceptibility to pathogens, insect attack, and fire damage. This review summarizes the strengths and shortcomings of previous work in the areas of genetic variation related to cold and drought stress in forest species with particular emphasis on loblolly pine (Pinus taeda L.), the most-planted tree species in North America. We describe and discuss the implementation of management and breeding strategies to increase resilience and adaptation, and discuss how new technologies in the areas of engineering and genomics are shaping the future of phenotype-genotype studies. Lessons learned from the study of species important in intensively-managed forest ecosystems may also prove to be of value in helping less-intensively managed forest ecosystems adapt to climate change, thereby increasing the sustainability and resilience of forestlands for the future.

Genome-wide profiling of alternative splicing genes in hybrid poplar (P.alba×P.glandulosa cv.84K) leaves

Alternative splicing (AS) is a post-transcriptional process common in plants and essential for regulation of environmental fitness of plants. In the present study, we focus on the AS events in poplar leaves to understand their effects on plant growth and development. The hybrid poplar (P.alba×P.glandulosa cv.84K) leaves were collected for RNA extraction. The extracted RNA was sequenced using on an Illumina HiSeq™ 2000 platform. Using the Populus trichocarpa genome as the reference, a total of 3810 AS genes were identified (9225 AS events), which accounted for 13.51% of all the expressed genes. Intron retention was the most common AS event, accounting for 43.86% of all the AS events, followed by alternative 3′ splice sites (23.75%), alternative 5′ splice sites (23.71%), and exon skipping (8.68%). Chromosomes 10 had the most condensed AS events (33.67 events/Mb) and chromosome 19 had the least (12.42 events/Mb). Association analysis showed that AS in the poplar leaves was positively correlated with intron length, exon number, exon length, and gene expression level, and was negatively correlated with GC content. AS genes in the poplar leaves were associated mainly with inositol phosphate metabolism and phosphatidylinositol signaling system pathways that would be significant on wooden plant production.

Alternative splicing and its regulatory role in woody plants

Alternative splicing (AS) is an important post-transcriptional process to enhance proteome diversity in eukaryotic organisms. In plants, numerous reports have primarily focused on AS analysis in model plant species or herbaceous plants, leading to a notable lack of research on AS in woody plants. More importantly, emerging evidence indicates that many important traits, including wood formation and stress resistance, in woody plants are controlled by AS. In this review article, we summarize the current progress of all kinds of AS studies in different tree species at various stages of development and in response to various stresses, revealing the significant role played by AS in woody plants, as well as the similar properties and differential regulation within their herbaceous counterparts. Furthermore, we propose several potential strategies to facilitate the functional characterization of splicing factors in woody plants and evaluate a general pipeline for the systematic characterization of splicing isoforms in a complex AS regulatory network. The utilization of genetic studies and high-throughput omics integration approaches to analyze AS genes and splicing factors is likely to further advance our understanding of AS modulation in woody plants.

Transcriptional and Post-transcriptional Regulation of Lignin Biosynthesis Pathway Genes in Populus

Lignin is a heterogeneous polymer of aromatic subunits derived from phenylalanine. It is polymerized in intimate proximity to the polysaccharide components in plant cell walls and provides additional rigidity and compressive strength for plants. Understanding the regulatory mechanisms of lignin biosynthesis is important for genetic modification of the plant cell wall for agricultural and industrial applications. Over the past 10 years the transcriptional regulatory model of lignin biosynthesis has been established in plants. However, the role of post-transcriptional regulation is still largely unknown. Increasing evidence suggests that lignin biosynthesis pathway genes are also regulated by alternative splicing, microRNA, and long non-coding RNA. In this review, we briefly summarize recent progress on the transcriptional regulation, then we focus on reviewing progress on the post-transcriptional regulation of lignin biosynthesis pathway genes in the woody model plant Populus.

Hybrid Sequencing of Full-Length cDNA Transcripts of the Medicinal Plant Scutellaria baicalensis

Scutellaria baicalensis is a well-known medicinal plant that produces biologically active flavonoids, such as baicalin, baicalein, and wogonin. Pharmacological studies have shown that these compounds have anti-inflammatory, anti-bacterial, and anti-cancer activities. Therefore, it is of great significance to investigate the genetic information of S. baicalensis, particularly the genes related to the biosynthetic pathways of these compounds. Here, we constructed the full-length transcriptome of S. baicalensis using a hybrid sequencing strategy and acquired 338,136 full-length sequences, accounting for 93.3% of the total reads. After the removal of redundancy and correction with Illumina short reads, 75,785 nonredundant transcripts were generated, among which approximately 98% were annotated with significant hits in the protein databases, and 11,135 sequences were classified as lncRNAs. Differentially expressed gene (DEG) analysis showed that most of the genes related to flavonoid biosynthesis were highly expressed in the roots, consistent with previous reports that the flavonoids were mainly synthesized and accumulated in the roots of S. baicalensis. By constructing unique transcription models, a total of 44,071 alternative splicing (AS) events were identified, with intron retention (IR) accounting for the highest proportion (44.5%). A total of 94 AS events were present in five key genes related to flavonoid biosynthesis, suggesting that AS may play important roles in the regulation of flavonoid biosynthesis in S. baicalensis. This study provided a large number of highly accurate full-length transcripts, which represents a valuable genetic resource for further research of the molecular biology of S. baicalensis, such as the development, breeding, and biosynthesis of active ingredients.

Computational analysis of alternative splicing in plant genomes

Computational analyses play crucial roles in characterizing splicing isoforms in plant genomes. In this review, we provide a survey of computational tools used in recently published, genome-scale splicing analyses in plants. We summarize the commonly used software and pipelines for read mapping, isoform reconstruction, isoform quantification, and differential expression analysis. We also discuss methods for analyzing long reads and the strategies to combine long and short reads in identifying splicing isoforms. We review several tools for characterizing local splicing events, splicing graphs, coding potential, and visualizing splicing isoforms. We further discuss the procedures for identifying conserved splicing isoforms across plant species. Finally, we discuss the outlook of integrating other genomic data with splicing analyses to identify regulatory mechanisms of AS on genome-wide scale.

Systems and Synthetic Biology of Forest Trees: A Bioengineering Paradigm for Woody Biomass Feedstocks

Fast-growing forest plantations are sustainable feedstocks of plant biomass that can serve as alternatives to fossil carbon resources for materials, chemicals, and energy. Their ability to efficiently harvest light energy and carbon from the atmosphere and sequester this into metabolic precursors for lignocellulosic biopolymers and a wide range of plant specialized metabolites make them excellent biochemical production platforms and living biorefineries. Their large sizes have facilitated multi-omics analyses and systems modeling of key biological processes such as lignin biosynthesis in trees. High-throughput 'omics' approaches have also been applied in segregating tree populations where genetic variation creates abundant genetic perturbations of system components allowing construction of systems genetics models linking genes and pathways to complex trait variation. With this information in hand, it is now possible to start using synthetic biology and genome editing techniques in a bioengineering approach based on a deeper understanding and rational design of biological parts, devices, and integrated systems. However, the complexity of the biology and interacting components will require investment in big data informatics, machine learning, and intuitive visualization to fully explore multi-dimensional patterns and identify emergent properties of biological systems. Predictive systems models could be tested rapidly through high-throughput synthetic biology approaches and multigene editing. Such a bioengineering paradigm, together with accelerated genomic breeding, will be crucial for the development of a new generation of woody biorefinery crops.

Expanding Alternative Splicing Identification by Integrating Multiple Sources of Transcription Data in Tomato

Tomato (Solanum lycopersicum) is an important vegetable and fruit crop. Its genome was completely sequenced and there are also a large amount of available expressed sequence tags (ESTs) and short reads generated by RNA sequencing (RNA-seq) technologies. Mapping transcripts including mRNA sequences, ESTs, and RNA-seq reads to the genome allows identifying pre-mRNA alternative splicing (AS), a post-transcriptional process generating two or more RNA isoforms from one pre-mRNA transcript. We comprehensively analyzed the AS landscape in tomato by integrating genome mapping information of all available mRNA and ESTs with mapping information of RNA-seq reads which were collected from 27 published projects. A total of 369,911 AS events were identified from 34,419 genomic loci involving 161,913 transcripts. Within the basic AS events, intron retention is the prevalent type (18.9%), followed by alternative acceptor site (12.9%) and alternative donor site (7.3%), with exon skipping as the least type (6.0%). Complex AS types having two or more basic event accounted for 54.9% of total AS events. Within 35,768 annotated protein-coding gene models, 23,233 gene models were found having pre-mRNAs generating AS isoform transcripts. Thus the estimated AS rate was 65.0% in tomato. The list of identified AS genes with their corresponding transcript isoforms serves as a catalog for further detailed examination of gene functions in tomato biology. The post-transcriptional information is also expected to be useful in improving the predicted gene models in tomato. The sequence and annotation information can be accessed at plant alternative splicing database (http://proteomics.ysu.edu/altsplice).

Abiotic Stresses Modulate Landscape of Poplar Transcriptome via Alternative Splicing, Differential Intron Retention, and Isoform Ratio Switching

Abiotic stresses affect plant physiology, development, growth, and alter pre-mRNA splicing. Western poplar is a model woody tree and a potential bioenergy feedstock. To investigate the extent of stress-regulated alternative splicing (AS), we conducted an in-depth survey of leaf, root, and stem xylem transcriptomes under drought, salt, or temperature stress. Analysis of approximately one billion of genome-aligned RNA-Seq reads from tissue- or stress-specific libraries revealed over fifteen millions of novel splice junctions. Transcript models supported by both RNA-Seq and single molecule isoform sequencing (Iso-Seq) data revealed a broad array of novel stress- and/or tissue-specific isoforms. Analysis of Iso-Seq data also resulted in the discovery of 15,087 novel transcribed regions of which 164 show AS. Our findings demonstrate that abiotic stresses profoundly perturb transcript isoform profiles and trigger widespread intron retention (IR) events. Stress treatments often increased or decreased retention of specific introns - a phenomenon described here as differential intron retention (DIR). Many differentially retained introns were regulated in a stress- and/or tissue-specific manner. A subset of transcripts harboring super stress-responsive DIR events showed persisting fluctuations in the degree of IR across all treatments and tissue types. To investigate coordinated dynamics of intron-containing transcripts in the study we quantified absolute copy number of isoforms of two conserved transcription factors (TFs) using Droplet Digital PCR. This case study suggests that stress treatments can be associated with coordinated switches in relative ratios between fully spliced and intron-retaining isoforms and may play a role in adjusting transcriptome to abiotic stresses.

Genome-wide gene expression analysis of amphioxus (Branchiostoma belcheri) following lipopolysaccharide challenge using strand-specific RNA-seq

Amphioxus is the closest living proxy for exploring the evolutionary origin of the immune system in vertebrates. To understand the immune responses of amphioxus to lipopolysaccharide (LPS), 5 ribosomal RNA (rRNA)-depleted libraries of amphioxus were constructed, including one control (0 h) library and 4 treatment libraries at 6, 12, 24, and 48 h post-injection (hpi) with LPS. The transcriptome of Branchiostoma belcheri was analyzed using strand-specific RNA sequencing technology (RNA-seq). A total of 6161, 6665, 7969, and 6447 differentially expressed genes (DEGs) were detected at 6, 12, 24, and 48 hpi, respectively, compared with expression levels at 0 h. We identified amphioxus genes active during the acute-phase response to LPS at different time points after stimulation. Moreover, to better visualize the resolution phase of the immune process during immune response, we identified 6057 and 5235 DEGs at 48 hpi by comparing with 6 and 24 hpi, respectively. Through real-time quantitative PCR (qRT-PCR) analysis of 12 selected DEGs, we demonstrated the accuracy of the RNA-seq data in this study. Functional enrichment analysis of DEGs demonstrated that most terms were related to defense and immune responses, disease and infection, cell apoptosis, and metabolism and catalysis. Subsequently, we identified 1330, 485, 670, 911, and 1624 time-specific genes (TSGs) at 0, 6, 12, 24, and 48 hpi. Time-specific terms at each of 5 time points were primarily involved in development, immune signaling, signal transduction, DNA repair and stability, and metabolism and catalysis, respectively. As this is the first study to report the transcriptome of an organism with primitive immunity following LPS challenge at multiple time points, it provides gene expression information for further research into the evolution of immunity in vertebrates.

Evolutionarily Conserved Alternative Splicing Across Monocots

One difficulty when identifying alternative splicing (AS) events in plants is distinguishing functional AS from splicing noise. One way to add confidence to the validity of a splice isoform is to observe that it is conserved across evolutionarily related species. We use a high throughput method to identify junction-based conserved AS events from RNA-Seq data across nine plant species, including five grass monocots (maize, sorghum, rice, Brachpodium, and foxtail millet), plus two nongrass monocots (banana and African oil palm), the eudicot Arabidopsis, and the basal angiosperm Amborella In total, 9804 AS events were found to be conserved between two or more species studied. In grasses containing large regions of conserved synteny, the frequency of conserved AS events is twice that observed for genes outside of conserved synteny blocks. In plant-specific RS and RS2Z subfamilies of the serine/arginine (SR) splice-factor proteins, we observe both conservation and divergence of AS events after the whole genome duplication in maize. In addition, plant-specific RS and RS2Z splice-factor subfamilies are highly connected with R2R3-MYB in STRING functional protein association networks built using genes exhibiting conserved AS. Furthermore, we discovered that functional protein association networks constructed around genes harboring conserved AS events are enriched for phosphatases, kinases, and ubiquitylation genes, which suggests that AS may participate in regulating signaling pathways. These data lay the foundation for identifying and studying conserved AS events in the monocots, particularly across grass species, and this conserved AS resource identifies an additional layer between genotype to phenotype that may impact future crop improvement efforts.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Transcript profiling of a novel plant meristem, the monocot cambium

While monocots lack the ability to produce a vascular cambium or woody growth, some monocot lineages evolved a novel lateral meristem, the monocot cambium, which supports secondary radial growth of stems. In contrast to the vascular cambium found in woody angiosperm and gymnosperm species, the monocot cambium produces secondary vascular bundles, which have an amphivasal organization of tracheids encircling a central strand of phloem. Currently there is no information concerning the molecular genetic basis of the development or evolution of the monocot cambium. Here we report high-quality transcriptomes for monocot cambium and early derivative tissues in two monocot genera, Yucca and Cordyline. Monocot cambium transcript profiles were compared to those of vascular cambia and secondary xylem tissues of two forest tree species, Populus trichocarpa and Eucalyptus grandis. Monocot cambium transcript levels showed that there are extensive overlaps between the regulation of monocot cambia and vascular cambia. Candidate regulatory genes that vary between the monocot and vascular cambia were also identified, and included members of the KANADI and CLE families involved in polarity and cell-cell signaling, respectively. We suggest that the monocot cambium may have evolved in part through reactivation of genetic mechanisms involved in vascular cambium regulation.

Alternative splicing enhances transcriptome complexity in desiccating seeds

Before being dispersed in the environment, mature seeds need to be dehydrated. The survival of seeds after dispersal depends on their low hydration in combination with high desiccation tolerance. These characteristics are established during seed maturation. Some key seed maturation genes have been reported to be regulated by alternative splicing (AS). However, so far AS was described only for single genes and a comprehensive analysis of AS during seed maturation has been lacking. We investigated gene expression and AS during Arabidopsis thaliana seed development at a global level, before and after desiccation. Bioinformatics tools were developed to identify differentially spliced regions within genes. Our data suggest the importance and shows the peculiar features of AS during seed desiccation. We identified AS in 34% of genes that are expressed at both timepoints before and after desiccation. Most of these AS transcript variants had not been found before in other tissues. Among the AS genes some seed master regulators could be found. Interestingly, 6% of all expressed transcripts were not transcriptionally regulated during desiccation, but only modified by AS. We propose that AS should be more routinely taken into account in the analysis of transcriptomic data to prevent overlooking potentially important regulators.

Evolutionary Character of Alternative Splicing in Plants

Alternative splicing (AS) is one of the most important ways to enhance the functional diversity of genes. Huge amounts of data have been produced by microarray, expressed sequence tag, and RNA-seq, and plenty of methods have been developed specifically for this task. The most frequently asked questions in previous research were as follows. What is the content rate of AS genes among the whole gene set? How many AS types are presented in the genome, and which type is dominant? How about the conservation ability of AS among different species? Which kinds of isoforms from some genes have the environmental response to help individual adaptation? Based on this background, we collected analysis results from 17 species to try to map out the landscape of AS studies in plants. We have noted the shortages of previous results, and we appeal to all scientists working in the AS field to make a standard protocol so that analyses between different projects are comparable.