Genome-Wide Analysis of Alternative Splicing during Development and Drought Stress in Maize

Chemical Modulation of Alternative Splicing for Molecular-Target Identification by Potential Genetic Control in Agrochemical Research

Alternative splicing (AS), the process of removing introns from pre-mRNA and the rearrangement of exons to produce several types of mature transcripts, is a remarkable step preceding protein synthesis. In particular, it has now been conclusively shown that up to ∼95% of genes are alternatively spliced to generate a complex and diverse proteome in eukaryotic organisms. Consequently, AS is one of the determinants of the functional repertoire of cells. Many studies have revealed that AS in plants can be regulated by cell type, developmental stage, environmental stress, and the circadian clock. Moreover, increasing amounts of evidence reveal that chemical compounds can affect various steps during splicing to induce major effects on plant physiology. Hence, the chemical modulation of AS can serve as a good strategy for molecular-target identification in attempts to potentially control plant genetics. However, the kind of mechanisms involved in the chemical modulation of AS that can be used in agrochemical research remain largely unknown. This review introduces recent studies describing the specific roles AS plays in plant adaptation to environmental stressors and in the regulation of development. We also discuss recent advances in small molecules that induce alterations of AS and the possibility of using this strategy in agrochemical-target identification, giving a new direction for potential genetic control in agrochemical research.

Developmental and transcriptional responses of maize to drought stress under field conditions

Drought is a common abiotic stress which significantly limits global crop productivity. Maize is an important staple crop and its yield is determined by successful development of the female inflorescence, the ear. We investigated drought stress responses across several developmental stages of the maize B73 inbred line under field conditions. Drought suppressed plant growth, but had little impact on progression through developmental stages. While ear growth was suppressed by drought, the process of spikelet initiation was not significantly affected. Tassel growth was reduced to a lesser extent compared to the observed reduction in ear growth under stress. Parallel RNA-seq profiling of leaves, ears, and tassels at several developmental stages revealed tissue-specific differences in response to drought stress. High temperature fluctuation was an additional environmental factor that also likely influenced gene expression patterns in the field. Drought induced significant transcriptional changes in leaves and ears but only minor changes in the tassel. Additionally, more genes were drought responsive in ears compared to leaves over the course of drought treatment. Genes that control DNA replication, cell cycle, and cell division were significantly down-regulated in stressed ears, which was consistent with inhibition of ear growth under drought. Inflorescence meristem genes were affected by drought to a lesser degree which was consistent with the minimal impact of drought on spikelet initiation. In contrast, genes that are involved in floret and ovule development were sensitive to stress, which is consistent with the detrimental effect of drought on gynoecium development and kernel set.

Comparative transcriptomics analysis uncovers alternative splicing events and molecular markers in cabbage (Brassica oleracea L.)

Alternative splicing (AS) events were identified and verified in cabbage by comparative transcriptome analysis. The corresponding markers were developed and the germplasm resources were identified. Alternative splicing (AS) is a central regulatory mechanism that greatly contributes to plant gene expression and transcriptome diversity. A large body of evidence has shown that AS complexity is relevant for plant development, evolution, complexity, and adaptation. Both insertion/deletion (InDel) and single nucleotide polymorphism (SNP) are typically co-dominant inheritance markers and have abundant polymorphisms. These have been widely used for marker-assisted selection, genetic mapping, and germplasm identification in plants. However, little is known about the molecular mechanisms underlying AS events and the development of markers including SNP and InDel from the cabbage transcriptome. In this study, three cabbage transcriptome datasets were collected and aligned to the cabbage reference genome to analyze AS events and marker development. 31,524 AS events were identified from three cabbage genotypes, accounting for 20.8% of the total cabbage genes. Alternative 3' splice site donor (A3SS) was the most frequent type of the four main AS events in cabbage. 70,475 InDels and 706,269 SNPs were identified with average frequencies of 1 InDel/6.9 kb and 1 SNP/0.7 kb, respectively. 71,942 potential SSRs were identified in 53,129 assembled unigenes with a density of 1 SSR/6.8 kb. The ratio of SNPs with synonymous/non-synonymous mutations was 1:0.65. 142 InDels and 36 SNPs were randomly selected and validated via Sanger sequencing and polymorphism was found among 66.2% of the InDels and 78.6% of the SNPs. Furthermore, 35 informative InDel markers were successfully used for genetic diversity analysis on 36 cabbage accessions. These results facilitate understanding of the molecular regulation mechanism underlying AS events in cabbage. They also provide molecular marker resource data for genetic mapping construction and germplasm identification, and facilitate the genetic improvement of cabbage via breeding.

An event of alternative splicing affects the expression of two BnCYCD3-1-like genes in Brassica napus

Two full-length cDNAs of the cyclin-D3-1-like gene, named as BnCYCD3-1-like-1 and BnCYCD3-1-like-2 respectively were obtained from Brassica napus, both of which encoded a cell cycle protein CYCD3. Alternative splicing (AS) events of the two genes' transcripts were identified, assigned as BnCYCD3-1-like-1-1, BnCYCD3-1-like-1-2, BnCYCD3-1-like-2-1 and BnCYCD3-1-like-2-2 respectively. BnCYCD3-1-like-1-1 and BnCYCD3-1-like-2-1 were both fully-spliced transcripts which encoded a complete protein containing a LXCXE motif, two cyclin boxes and a PEST domain, while other two alternative splicing transcripts both resulted in the early termination of the protein translation. BnCYCD3-1-like-2-2 retained the third intron, lacking a PEST domain, while BnCYCD3-1-like-1-2 retained all the introns, lacking the C-terminal cyclin domain and a PEST domain. The expression pattern for tissue and development specification of the AS transcripts were investigated. The results showed that the standard splicing transcripts (BnCYCD3-1-like-1-1 and BnCYCD3-1-like-2-1) with complete structural domains were found with the most abundant expression in seeds, followed by leaves, and the least expression in stems. Both of BnCYCD3-1-like-2-1 and BnCYCD3-1-like-2-2 had the highest abundance in leaves, followed by roots. In addition, by applying various biotic and abiotic stresses on Brassica napus, the variations in the expression of each transcript under stress treatment were studied. Also, it was found that AS of the cyclin-D3-1-like gene may play an important role in helping Brassica napus respond to environmental stresses by coordinating the levels of transcripts of standard splicing and alternative splicing.

Spaceflight-induced alternative splicing during seedling development in Arabidopsis thaliana

Plants grown in spaceflight experience novel environmental signals, including those associated with microgravity and ionizing radiation. Spaceflight triggers a response involving transcriptional re-programming and altered cell morphology, though many aspects of this response remain uncharacterized. We analyzed the spaceflight-induced transcriptome with a focus on genes that undergo alternative splicing to examine differential splicing associated with spaceflight—an unstudied characteristic of the molecular response to spaceflight exposure. RNA sequence data obtained during the APEX03 spaceflight experiment that was collected from two Arabidopsis thaliana ecotypes at two seedling stages grown onboard the International Space Station, or as ground controls at Kennedy Space Center, were re-examined to detect alternative splicing differences induced by spaceflight. Presence/absence variation analysis was used to identify putative expression-level differences in alternatively spliced isoforms between spaceflight and ground controls and was followed by analysis of significant differential alternative splicing. This study provides the first evidence of a role for alternative splicing in the molecular processes of physiological adaptation to the spaceflight environment.

Genome-wide identification and comparative analysis of alternative splicing across four legume species

MAIN CONCLUSION

Alternative splicing EVENTS were genome-wide identified for four legume species, and nitrogen fixation-related gene families and evolutionary analysis was also performed. Alternative splicing (AS) is a key regulatory mechanism that contributes to transcriptome and proteome diversity. Investigation of the genome-wide conserved AS events across different species will help with the understanding of the evolution of the functional diversity in legumes, allowing for genetic improvement. Genome-wide identification and characterization of AS were performed using the publically available mRNA, EST, and RNA-Seq data for four important legume species. A total of 15,165 AS genes in Glycine max, 6077 in Cicer arietinum, 7240 in Medicago truncatula, and 7358 in Lotus japonicus were identified. Intron retention (IntronR) was the dominant AS type among the identified events, with IntronR occurring from 53.76% in M. truncatula to 43.91% in C. arietinum. We identified 1159 AS genes that were conserved among four species. Furthermore, nine nitrogen fixation-related gene families with 237 genes were identified, and 80 of them were AS, accounting for the 43.48% in G. max and 27.78% in C. arietinum. An evolutionary analysis showed that these AS genes tended to be located adjacent to each other in the evolutionary tree and are unbalanced in the distribution in the sub-family. This study provides a foundation for future studies on transcription complexity, evolution, and the role of AS on plant functional regulation.

Differential Alternative Splicing Genes in Response to Boron Deficiency in Brassica napus

Alternative splicing (AS) can increase transcriptome diversity, protein diversity and protein yield, and is an important mechanism to regulate plant responses to stress. Oilseed rape (Brassica napus L.), one of the main oil crops in China, shows higher sensitivity to boron (B) deficiency than other species. Here, we demonstrated AS changes that largely increased the diversity of the mRNA expressed in response to B deficiency in B. napus. Each gene had two or more transcripts on average. A total of 33.3% genes in both Qingyou10 (QY10, B-efficient cultivar) and Westar10 (W10, B-inefficient cultivar) showed AS in both B conditions. The types of AS events were mainly intron retention, 3′ alternative splice site, 5′ alternative splice site and exon skipping. The tolerance ability of QY10 was higher than that of W10, possibly because there were far more differential alternative splicing (DAS) genes identified in QY10 at low B conditions than in W10. The number of genes with both DAS and differentially expressed (DE) was far lower than that of the genes that were either with DAS or DE in QY10 and W10, suggesting that the DAS and DE genes were independent. Four Serine/Arginine-rich (SR) splicing factors, BnaC06g14780D, BnaA01g14750D, BnaA06g15930D and BnaC01g41640D, underwent differentially alternative splicing in both cultivars. There existed gene–gene interactions between BnaC06g14780D and the genes associated with the function of B in oilseed rape at low B supply. This suggests that oilseed rape could regulate the alterative pre-mRNA splicing of SR protein related genes to increase the plant tolerance to B deficiency.

Patterns of alternative splicing in response to cold acclimation in fish

Phenotypic plasticity is an important aspect of an organism's response to environmental change that often requires the modulation of gene expression. These changes in gene expression can be quantitative, as a result of increases or decreases in the amounts of specific transcripts, or qualitative, as a result of the expression of alternative transcripts from the same gene (e.g. via alternative splicing of pre-mRNAs). Although the role of quantitative changes in gene expression in phenotypic plasticity is well known, relatively few studies have examined the role of qualitative changes. Here, we use skeletal muscle RNA-seq data from Atlantic killifish (Fundulus heteroclitus), threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio) to investigate the extent of qualitative changes in gene expression in response to cold acclimation. Fewer genes demonstrated alternative splicing than differential expression as a result of cold acclimation; however, differences in splicing were detected for 426 to 866 genes depending on species, indicating that large numbers of qualitative changes in gene expression are associated with cold acclimation. Many of these alternatively spliced genes were also differentially expressed, and there was functional enrichment for involvement in muscle contraction among the genes demonstrating qualitative changes in response to cold acclimation. Additionally, there was a common group of 29 genes with cold-acclimation-mediated changes in splicing in all three species, suggesting that there may be a set of genes with expression patterns that respond qualitatively to prolonged exposure to cold temperatures across fishes.

Iso-Seq Allows Genome-Independent Transcriptome Profiling of Grape Berry Development

Transcriptomics has been widely applied to study grape berry development. With few exceptions, transcriptomic studies in grape are performed using the available genome sequence, PN40024, as reference. However, differences in gene content among grape accessions, which contribute to phenotypic differences among cultivars, suggest that a single reference genome does not represent the species' entire gene space. Though whole genome assembly and annotation can reveal the relatively unique or "private" gene space of any particular cultivar, transcriptome reconstruction is a more rapid, less costly, and less computationally intensive strategy to accomplish the same goal. In this study, we used single molecule-real time sequencing (SMRT) to sequence full-length cDNA (Iso-Seq) and reconstruct the transcriptome of Cabernet Sauvignon berries during berry ripening. In addition, short reads from ripening berries were used to error-correct low-expression isoforms and to profile isoform expression. By comparing the annotated gene space of Cabernet Sauvignon to other grape cultivars, we demonstrate that the transcriptome reference built with Iso-Seq data represents most of the expressed genes in the grape berries and includes 1,501 cultivar-specific genes. Iso-Seq produced transcriptome profiles similar to those obtained after mapping on a complete genome reference. Together, these results justify the application of Iso-Seq to identify cultivar-specific genes and build a comprehensive reference for transcriptional profiling that circumvents the necessity of a genome reference with its associated costs and computational weight.

Genome-Wide Identification, Expression Profile, and Alternative Splicing Analysis of the Brassinosteroid-Signaling Kinase (BSK) Family Genes in Arabidopsis

Brassinosteroids (BRs) are steroid hormones essential for different biological processes, ranging from growth to environmental adaptation in plants. The plant brassinosteroid-signaling kinase (BSK) proteins belong to a family of receptor-like cytoplasmic kinases, which have been reported to play an important role in BR signal transduction. However, the knowledge of BSK genes in plants is still quite limited. In the present study, a total of 143 BSK proteins were identified by a genome-wide search in 17 plant species. A phylogenetic analysis showed that the BSK gene originated in embryophytes, with no BSK found in green algae, and these BSK genes were divided into six groups by comparison with orthologs/paralogs. A further study using comparative analyses of gene structure, expression patterns and alternative splicing of BSK genes in Arabidopsis revealed that all BSK proteins shared similar protein structure with some exception and post-translation modifications including sumolyation and ubiquitination. An expression profile analysis showed that most Arabidopsis BSK genes were constitutively expressed in different tissues; of these, several BSK genes were significantly expressed in response to some hormones or abiotic stresses. Furthermore, reverse transcription-polymerase chain reaction (RT-PCR) assays showed that BSK5, BSK7, and BSK9 underwent alternative splicing in specific stress induced and tissue-dependent patterns. Collectively, these results lay the foundation for further functional analyses of these genes in plants.

Temporal Small RNA Expression Profiling under Drought Reveals a Potential Regulatory Role of Small Nucleolar RNAs in the Drought Responses of Maize

Small RNAs (sRNAs) are short noncoding RNAs that play roles in many biological processes, including drought responses in plants. However, how the expression of sRNAs dynamically changes with the gradual imposition of drought stress in plants is largely unknown. We generated time-series sRNA sequence data from maize ( L.) seedlings under drought stress (DS) and under well-watered (WW) conditions at the same time points. Analyses of length, functional annotation, and abundance of 736,372 nonredundant sRNAs from both DS and WW data, as well as genome copy numbers at the corresponding genomic regions, revealed distinct patterns of abundance and genome organization for different sRNA classes. The analysis identified 6646 sRNAs whose regulation was altered in response to drought stress. Among drought-responsive sRNAs, 1325 showed transient downregulation by the seventh day, coinciding with visible symptoms of drought stress. The profiles revealed drought-responsive microRNAs, as well as other sRNAs that originated from ribosomal RNAs (rRNAs), splicing small nuclear RNAs, and small nucleolar RNAs (snoRNA). Expression profiles of their sRNA derivers indicated that snoRNAs might play a regulatory role through regulating the stability of rRNAs and splicing small nuclear RNAs under drought condition.

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.

Spatio-Temporal Transcriptional Dynamics of Maize Long Non-Coding RNAs Responsive to Drought Stress

Long non-coding RNAs (lncRNAs) have emerged as important regulators in plant stress response. Here, we report a genome-wide lncRNA transcriptional analysis in response to drought stress using an expanded series of maize samples collected from three distinct tissues spanning four developmental stages. In total, 3488 high-confidence lncRNAs were identified, among which 1535 were characterized as drought responsive. By characterizing the genomic structure and expression pattern, we found that lncRNA structures were less complex than protein-coding genes, showing shorter transcripts and fewer exons. Moreover, drought-responsive lncRNAs exhibited higher tissue- and development-specificity than protein-coding genes. By exploring the temporal expression patterns of drought-responsive lncRNAs at different developmental stages, we discovered that the reproductive stage R1 was the most sensitive growth stage with more lncRNAs showing altered expression upon drought stress. Furthermore, lncRNA target prediction revealed 653 potential lncRNA-messenger RNA (mRNA) pairs, among which 124 pairs function in cis-acting mode and 529 in trans. Functional enrichment analysis showed that the targets were significantly enriched in molecular functions related to oxidoreductase activity, water binding, and electron carrier activity. Multiple promising targets of drought-responsive lncRNAs were discovered, including the V-ATPase encoding gene, vpp4. These findings extend our knowledge of lncRNAs as important regulators in maize drought response.

Alternative splicing and translation play important roles in hypoxic germination in rice

Post-transcriptional mechanisms (PTMs), including alternative splicing (AS) and alternative translation initiation (ATI), may explain the diversity of proteins involved in plant development and stress responses. Transcriptional regulation is important during the hypoxic germination of rice seeds, but the potential roles of PTMs in this process have not been characterized. We used a combination of proteomics and RNA sequencing to discover how AS and ATI contribute to plant responses to hypoxia. In total, 10 253 intron-containing genes were identified. Of these, ~1741 differentially expressed AS (DAS) events from 811 genes were identified in hypoxia-treated seeds compared with controls. Over 95% of these were not present in the list of differentially expressed genes. In particular, regulatory pathways such as the spliceosome, ribosome, endoplasmic reticulum protein processing and export, proteasome, phagosome, oxidative phosphorylation, and mRNA surveillance showed substantial AS changes under hypoxia, suggesting that AS responses are largely independent of transcriptional regulation. Considerable AS changes were identified, including the preferential usage of some non-conventional splice sites and enrichment of splicing factors in the DAS data sets. Taken together, these results not only demonstrate that AS and ATI function during hypoxic germination but they have also allowed the identification of numerous novel proteins/peptides produced via ATI.

The chromatin remodeler ZmCHB101 impacts alternative splicing contexts in response to osmotic stress

Maize SWI3-type chromatin remodeler impacts alternative splicing contexts in response to osmotic stress by altering nucleosome density and affecting transcriptional elongation rate. Alternative splicing (AS) is commonly found in higher eukaryotes and is an important posttranscriptional regulatory mechanism to generate transcript diversity. AS has been widely accepted as playing essential roles in different biological processes including growth, development, signal transduction and responses to biotic and abiotic stresses in plants. However, whether and how chromatin remodeling complex functions in AS in plant under osmotic stress remains unknown. Here, we show that a maize SWI3D protein, ZmCHB101, impacts AS contexts in response to osmotic stress. Genome-wide analysis of mRNA contexts in response to osmotic stress using ZmCHB101-RNAi lines reveals that ZmCHB101 impacts alternative splicing contexts of a subset of osmotic stress-responsive genes. Intriguingly, ZmCHB101-mediated regulation of gene expression and AS is largely uncoupled, pointing to diverse molecular functions of ZmCHB101 in transcriptional and posttranscriptional regulation. We further found ZmCHB101 impacts the alternative splicing contexts by influencing alteration of chromatin and histone modification status as well as transcriptional elongation rates mediated by RNA polymerase II. Taken together, our findings suggest a novel insight of how plant chromatin remodeling complex impacts AS under osmotic stress .

Identification of Novel lincRNA and Co-Expression Network Analysis Using RNA-Sequencing Data in Plants

Long intergenic noncoding RNA (lincRNA) plays important biological functions in plants. Identification and annotation of lincRNA in plants largely rely on RNA sequencing followed by computational analysis. In this protocol, we describe a multistep computational pipeline for lincRNA identification using RNA-sequencing data. This pipeline can also construct co-expression network that is made of both lincRNA and mRNA genes. The co-expression network generated by this pipeline can be used to provide putative annotation of lincRNAs that have no known biological functions.

Functional identification of lncRNAs in sweet cherry (Prunus avium) pollen tubes via transcriptome analysis using single-molecule long-read sequencing

Sweet cherry (Prunus avium) is a popular fruit with high nutritional value and excellent flavor. Although pollen plays an important role in the double fertilization and subsequent fruit production of this species, little is known about its pollen tube transcriptome. In this study, we identified 16,409 transcripts using single-molecule sequencing. After filtering 292 transposable elements, we conducted further analyses including mRNA classification, gene function prediction, alternative splicing (AS) analysis, and long noncoding RNA (lncRNA) identification to gain insight into the pollen transcriptome. The filtered transcripts could be matched with 3,438 coding region sequences from the sweet cherry genome. GO and KEGG analyses revealed complex biological processes during pollen tube elongation. A total of 2043 AS events were predicted, 7 of which were identified in different organs, such as the leaf, pistil and pollen tube. Using BLASTnt and the Coding-Potential Assessment Tool (CPAT), we distinguished a total of 284 lncRNAs, among which 154 qualified as natural antisense transcripts (NATs). As the NATs could be the reverse complements of coding mRNA sequences, they might bind to coding sequences. Antisense transfection assays showed that the NATs could regulate the expression levels of their complementary sequences and even affect the growth conditions of pollen tubes. In summary, this research characterizes the transcripts of P. avium pollen and lays the foundation for elucidating the physiological and biochemical mechanisms underlying sexual reproduction in the male gametes of this species.

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).

Insights into the transcriptional and post-transcriptional regulation of the rice SUMOylation machinery and into the role of two rice SUMO proteases

BACKGROUND

SUMOylation is an essential eukaryotic post-translation modification that, in plants, regulates numerous cellular processes, ranging from seed development to stress response. Using rice as a model crop plant, we searched for potential regulatory points that may influence the activity of the rice SUMOylation machinery genes.

RESULTS

We analyzed the presence of putative cis-acting regulatory elements (CREs) within the promoter regions of the rice SUMOylation machinery genes and found CREs related to different cellular processes, including hormone signaling. We confirmed that the transcript levels of genes involved in target-SUMOylation, containing ABA- and GA-related CREs, are responsive to treatments with these hormones. Transcriptional analysis in Nipponbare (spp. japonica) and LC-93-4 (spp. indica), showed that the transcript levels of all studied genes are maintained in the two subspecies, under normal growth. OsSUMO3 is an exceptional case since it is expressed at low levels or is not detectable at all in LC-93-4 roots and shoots, respectively. We revealed post-transcriptional regulation by alternative splicing (AS) for all genes studied, except for SUMO coding genes, OsSIZ2, OsOTS3, and OsELS2. Some AS forms have the potential to alter protein domains and catalytic centers. We also performed the molecular and phenotypic characterization of T-DNA insertion lines of some of the genes under study. Knockouts of OsFUG1 and OsELS1 showed increased SUMOylation levels and non-overlapping phenotypes. The fug1 line showed a dwarf phenotype, and significant defects in fertility, seed weight, and panicle architecture, while the els1 line showed early flowering and decreased plant height. We suggest that OsELS1 is an ortholog of AtEsd4, which was also supported by our phylogenetic analysis.

CONCLUSIONS

Overall, we provide a comprehensive analysis of the rice SUMOylation machinery and discuss possible effects of the regulation of these genes at the transcriptional and post-transcriptional level. We also contribute to the characterization of two rice SUMO proteases, OsELS1 and OsFUG1.

Alternative Splicing as a Regulator of Early Plant Development

Most plant genes are interrupted by introns and the corresponding transcripts need to undergo pre-mRNA splicing to remove these intervening sequences. Alternative splicing (AS) is an important posttranscriptional process that creates multiple mRNA variants from a single pre-mRNA molecule, thereby enhancing the coding and regulatory potential of genomes. In plants, this mechanism has been implicated in the response to environmental cues, including abiotic and biotic stresses, in the regulation of key developmental processes such as flowering, and in circadian timekeeping. The early plant development steps - from embryo formation and seed germination to skoto- and photomorphogenesis - are critical to both execute the correct body plan and initiate a new reproductive cycle. We review here the available evidence for the involvement of AS and various splicing factors in the initial stages of plant development, while highlighting recent findings as well as potential future challenges.

Transcriptomic analysis of drought stress responses of sea buckthorn (Hippophae rhamnoidessubsp. sinensis) by RNA-Seq

Sea buckthorn is one of the most important eco-economic tree species in China due to its ability to grow and produce acceptable yields under limited water and fertilizer availability. In this study, the differentially expressed genes under drought stress (DS) of sea buckthorn were identified and compared with control (CK) by RNA-Seq. A total of 122,803 unigenes were identified in sea buckthorn, and 70,025 unigenes significantly matched a sequence in at least one of the seven databases. A total of 24,060 (19.59%) unigenes can be assigned to 19 KEGG pathways, and 1,644 unigenes were differentially expressed between DS and CK, of which 519 unigenes were up-regulated and 1,125 unigenes down-regulated. Of the 47 significantly enriched GO terms, 14, 7 and 26 items were related to BP, CC and MF, respectively. KEGG enrichment analysis showed 398 DEGs involved in 97 different pathways, of which 119 DEGs were up-regulated and 279 DEGs were down-regulated under drought stress. In addition, we found 4438 transcriptor factors (TFs) in sea buckthorn, of which 100 were differentially expressed between DS and CK. These results lay a first foundation for further investigations of the very specific functions of these unigenes in sea buckthorn in response to drought stress.

Transcriptome analysis of alternative splicing in peanut (Arachis hypogaea L.)

BACKGROUND

Alternative splicing (AS) represents a mechanism widely used by eukaryotes for the post-transcriptional regulation of genes. The detailed exploration of AS in peanut has not been documented.

RESULTS

The strand-specific RNA-Seq technique was exploited to characterize the distribution of AS in the four samples of peanut (FH1-seed1, FH1-seed2, FH1-root and FH1-leaf). AS was detected as affecting around 37.2% of the full set of multi-exon genes. Some of these genes experienced AS throughout the plant, while in the case of others, the effect was organ-specific. Overall, AS was more frequent in the seed than in either the root or leaf. The predominant form of AS was intron retention, and AS in transcription start site and transcription terminal site were commonly identified in all the four samples. It is interesting that in genes affected by AS, the majority experienced only a single type of event. Not all of the in silico predicted transcripts appeared to be translated, implying that these are either degraded or sequestered away from the translation machinery. With respect to genes involved in fatty acid metabolism, about 61.6% were shown to experience AS.

CONCLUSION

Our report contributes significantly in AS analysis of peanut genes in general, and these results have not been mentioned before. The specific functions of different AS forms need further investigation.

Imprinted gene expression in maize starchy endosperm and aleurone tissues of reciprocal F1 hybrids at a defined developmental stage

Imprinted gene expression in flowering plants predominantly occurs in the triploid endosperm of developing seed. However, endosperm is composed of distinct tissue types. For example, the maize (Zea mays) endosperm is constituted by two major tissues, starchy endosperm and aleurone. Previous studies in imprinted gene expression have generally assumed that the different tissues constituting endosperm would behavior the same, and hence have not examined them separately. Here, to examine parental-specific expression of imprinted genes in different parts of the seed, eight previously reported maize protein-coding imprinted genes were selected, and analyzed by cleaved amplified polymorphic sequence (CAPS) coupled with Sanger sequencing for transcripts from the various seed tissues collected at 18 days after pollination (DAP). The studied tissues included seed coat, embryo, starchy endosperm and aleurone, which were collected from a pair of reciprocal F1 hybrids produced by crossing inbred lines B73 and Mo17. Six of these eight analyzed imprinted genes showed the same imprinted expression pattern between the starchy endosperm and aleurone, but two showed imprinted expression only in the starchy endosperm. Comparison of the expression pattern of 20 selected imprinted genes in multiple seed tissues and vegetative tissues indicated that the majority (~ 75%) of these imprinted genes exhibited seed-specific or endosperm-specific expression. Our results also uncovered that imprinted genes have a high propensity to be alternatively spliced via intron retention in the developing embryo compared with the other tissues.

Genetics of alternative splicing evolution during sunflower domestication

Alternative splicing enables organisms to produce the diversity of proteins necessary for multicellular life by using relatively few protein-coding genes. Although differences in splicing have been identified among divergent taxa, the shorter-term evolution of splicing is understudied. The origins of novel splice forms, and the contributions of alternative splicing to major evolutionary transitions, are largely unknown. This study used transcriptomes of wild and domesticated sunflowers to examine splice differentiation and regulation during domestication. We identified substantial splicing divergence between wild and domesticated sunflowers, mainly in the form of intron retention. Transcripts with divergent splicing were enriched for seed-development functions, suggesting that artificial selection impacted splicing patterns. Mapping of quantitative trait loci (QTLs) associated with 144 differential splicing cases revealed primarily trans-acting variation affecting splicing patterns. A large proportion of identified QTLs contain known spliceosome proteins and are associated with splicing variation in multiple genes. Examining a broader set of wild and domesticated sunflower genotypes revealed that most differential splicing patterns in domesticated sunflowers likely arose from standing variation in wild Helianthus annuus and gained frequency during the domestication process. However, several domesticate-associated splicing patterns appear to be introgressed from other Helianthus species. These results suggest that sunflower domestication involved selection on pleiotropic regulatory alleles. More generally, our findings indicate that substantial differences in isoform abundances arose rapidly during a recent evolutionary transition and appear to contribute to adaptation and population divergence.

Structural, Functional, and Evolutionary Characterization of Major Drought Transcription Factors Families in Maize

Drought is one of the major threats to the maize yield especially in subtropical production systems. Understanding the genes and regulatory mechanisms of drought tolerance is important to sustain the yield. Transcription factors (TFs) play a major role in gene regulation under drought stress. In the present study, a set of 15 major TF families comprising 1,436 genes was structurally and functionally characterized. The functional annotation indicated that the genes were involved in ABA signaling, ROS scavenging, photosynthesis, stomatal regulation, and sucrose metabolism. Duplication was identified as the primary force in divergence and expansion of TF families. Phylogenetic relationship was developed for individual TF and combined TF families. Phylogenetic analysis clustered the genes into specific and mixed groups. Gene structure analysis revealed that more number of genes were intron-rich as compared to intron-less. Drought-responsive cis-regulatory elements such as ABREA, ABREB, DRE1, and DRECRTCOREAT have been identified. Expression and interaction analyses identified leaf-specific bZIP TF, GRMZM2G140355, as a potential contributor toward drought tolerance in maize. Protein-protein interaction network of 269 drought-responsive genes belonging to different TFs has been provided. The information generated on structural and functional characteristics, expression, and interaction of the drought-related TF families will be useful to decipher the drought tolerance mechanisms and to breed drought-tolerant genotypes in maize.

Drought-tolerant and drought-sensitive genotypes of maize (Zea mays L.) differ in contents of endogenous brassinosteroids and their drought-induced changes

The contents of endogenous brassinosteroids (BRs) together with various aspects of plant morphology, water management, photosynthesis and protection against cell damage were assessed in two maize genotypes that differed in their drought sensitivity. The presence of 28-norbrassinolide in rather high quantities (1-2 pg mg-1 fresh mass) in the leaves of monocot plants is reported for the first time. The intraspecific variability in the presence/content of the individual BRs in drought-stressed plants is also described for the first time. The drought-resistant genotype was characterised by a significantly higher content of total endogenous BRs (particularly typhasterol and 28-norbrassinolide) compared with the drought-sensitive genotype. On the other hand, the drought-sensitive genotype showed higher levels of 28-norcastasterone. Both genotypes also differed in the drought-induced reduction/elevation of the levels of 28-norbrassinolide, 28-norcastasterone, 28-homocastasterone and 28-homodolichosterone. The differences observed between both genotypes in the endogenous BR content are probably correlated with their different degrees of drought sensitivity, which was demonstrated at various levels of plant morphology, physiology and biochemistry.

Capturing the Alternative Cleavage and Polyadenylation Sites of 14 NAC Genes in Populus Using a Combination of 3'-RACE and High-Throughput Sequencing

Detection of complex splice sites (SSs) and polyadenylation sites (PASs) of eukaryotic genes is essential for the elucidation of gene regulatory mechanisms. Transcriptome-wide studies using high-throughput sequencing (HTS) have revealed prevalent alternative splicing (AS) and alternative polyadenylation (APA) in plants. However, small-scale and high-depth HTS aimed at detecting genes or gene families are very few and limited. We explored a convenient and flexible method for profiling SSs and PASs, which combines rapid amplification of 3′-cDNA ends (3′-RACE) and HTS. Fourteen NAC (NAM, ATAF1/2, CUC2) transcription factor genes of Populus trichocarpa were analyzed by 3′-RACE-seq. Based on experimental reproducibility, boundary sequence analysis and reverse transcription PCR (RT-PCR) verification, only canonical SSs were considered to be authentic. Based on stringent criteria, candidate PASs without any internal priming features were chosen as authentic PASs and assumed to be PAS-rich markers. Thirty-four novel canonical SSs, six intronic/internal exons and thirty 3′-UTR PAS-rich markers were revealed by 3′-RACE-seq. Using 3′-RACE and real-time PCR, we confirmed that three APA transcripts ending in/around PAS-rich markers were differentially regulated in response to plant hormones. Our results indicate that 3′-RACE-seq is a robust and cost-effective method to discover SSs and label active regions subjected to APA for genes or gene families. The method is suitable for small-scale AS and APA research in the initial stage.

Evolution of intron-poor clades and expression patterns of the glycosyltransferase family 47

A large-scale bioinformatics analysis revealed the origin and evolution of GT47 gene family, and identified two clades of intron-poor genes with putative functions in drought stress responses and seed development in maize. Glycosyltransferase family 47 (GT47) genes encode β-galactosyltransferases and β-glucuronyltransferases that synthesize pectin, xyloglucans and xylan, which are important components of the plant cell wall. In this study, we performed a systematic and large-scale bioinformatics analysis of GT47 gene family using 352 GT47 proteins from 15 species ranging from cyanobacteria to seed plants. The analysis results showed that GT47 family may originate in cyanobacteria and expand along the evolutionary trajectory to moss. Further analysis of 47 GT47 genes in maize revealed that they can divide into five clades with diverse exon-intron structures. Among these five clades, two were mainly composed with intron-poor genes, which may originate in the moss. Gene duplication analysis revealed that the expansion of GT47 gene family in maize was significantly driven from tandem duplication events and segmental duplication events. Significantly, almost all duplicated genes are intron-poor genes. Expression analysis indicated that several intron-poor GT47 genes may be involved in the drought stress response and seed development in maize. This work provides insight into the origin and evolutionary process, expansion mechanisms and expression patterns of GT47 genes, thus facilitating their functional investigations in the future.

Global profiling of alternative splicing landscape responsive to drought, heat and their combination in wheat (Triticum aestivum L.)

Plant can acquire tolerance to environmental stresses via transcriptome reprogramming at transcriptional and alternative splicing (AS) levels. However, how AS coordinates with transcriptional regulation to contribute to abiotic stresses responses is still ambiguous. In this study, we performed genome-wide analyses of AS responses to drought stress (DS), heat stress (HS) and their combination (HD) in wheat seedlings, and further compared them with transcriptional responses. In total, we found 200, 3576 and 4056 genes exhibiting significant AS pattern changes in response to DS, HS and HD, respectively, and combined drought and heat stress can induce specific AS compared with individual one. In addition, wheat homeologous genes exhibited differential AS responses under stress conditions that more AS events occurred on B subgenome than on A and D genomes. Comparison of genes regulated at AS and transcriptional levels showed that only 12% of DS-induced AS genes were subjected to transcriptional regulation, whereas the proportion increased to ~40% under HS and HD. Functional enrichment analysis revealed that abiotic stress-responsive pathways tended to be highly overrepresented among these overlapped genes under HS and HD. Thus, we proposed that transcriptional regulation may play a major role in response to DS, which coordinates with AS regulation to contribute to HS and HD tolerance in wheat.

Comparative epigenomics reveals evolution of duplicated genes in potato and tomato

The evolution of duplicated genes after polyploidization has been the subject of many evolutionary biology studies. Potato (Solanum tuberosum) and tomato (Solanum lycopersicum) are the first two sequenced genomes of asterids, and share a common polyploidization event. However, the epigenetic role of DNA methylation on the evolution of duplicated genes derived from polyploidization is not fully understood. Here, we explore the role of the DNA methylation in the evolution of duplicated genes in potato and tomato. The overall levels of DNA methylation are different, although patterns of DNA methylation are similar in potato and tomato. Different types of duplicated genes can display different methylation patterns in potato and tomato. In addition, we found that differences in the methylation levels between duplicated genes were associated with gene expression divergence. In particular, for the majority of duplicated gene pairs, one copy is always hyper- or hypo-methylated compared with the other copy across different tomato fruit ripening stages, and these genes are enriched for specific function related to transcription factor (TF) activity. Furthermore, transcription of hundreds of duplicated TFs was shown to be regulated by DNA methylation during fruit ripening stages in tomato, some of which are well-known fruit ripening TFs. Taken together, our results support the notion that DNA methylation may facilitate divergent evolution of duplicated genes and play roles in important biological processes such as tomato fruit ripening.

Alternative Splicing Control of Abiotic Stress Responses

Alternative splicing, which generates multiple transcripts from the same gene, is an important modulator of gene expression that can increase proteome diversity and regulate mRNA levels. In plants, this post-transcriptional mechanism is markedly induced in response to environmental stress, and recent studies have identified alternative splicing events that allow rapid adjustment of the abundance and function of key stress-response components. In agreement, plant mutants defective in splicing factors are severely impaired in their response to abiotic stress. Notably, mounting evidence indicates that alternative splicing regulates stress responses largely by targeting the abscisic acid (ABA) pathway. We review here current understanding of post-transcriptional control of plant stress tolerance via alternative splicing and discuss research challenges for the near future.

RNA-seq analysis reveals alternative splicing under salt stress in cotton, Gossypium davidsonii

BACKGROUND

Numerous studies have focused on the regulation of gene expression in response to salt stress at the transcriptional level; however, little is known about this process at the post-transcriptional level.

RESULTS

Using a diploid D genome wild salinity-tolerant cotton species, Gossypium davidsonii, we analyzed alternative splicing (AS) of genes related to salt stress by comparing high-throughput transcriptomes from salt-treated and well-watered roots and leaves. A total of 14,172 AS events were identified involving 6798 genes, of which intron retention (35.73%) was the most frequent, being detected in 3492 genes. Under salt stress, 1287 and 1228 differential alternative splicing (DAS) events were identified in roots and leaves, respectively. These DAS genes were associated with specific functional pathways, such as "responses to stress", "metabolic process" and "RNA splicing", implying that AS represents an important pathway of gene regulation in response to salt stress. Several salt response genes, such as pyrroline-5-carboxylate synthase (P5CS), K⁺ channel outward (KCO1), plasma membrane intrinsic protein (PIP) and WRKY33 which were involved in osmotic balance, ion homeostasis, water transportation and transcriptional regulation, respectively, were identified with differential alternative splicing under salt stress. Moreover, we revealed that 13 genes encoding Ser/Arg-rich (SR) proteins related to AS regulation were differentially alternatively spliced under salt stress.

CONCLUSION

This study first provide a comprehensive view of AS in G. davidsonii, and highlight novel insights into the potential roles of AS in plant responses to salt stress.

The comparison of alternative splicing among the multiple tissues in cucumber

BACKGROUND

Alternative splicing (AS) is an important post-transcriptional process. It has been suggested that most AS events are subject to tissue-specific regulation. However, the global dynamics of AS in different tissues are poorly explored.

RESULTS

To analyse global changes in AS in multiple tissues, we identified the AS events and constructed a comprehensive catalogue of AS events within each tissue based on the genome-wide RNA-seq reads from ten tissues in cucumber. First, we found that 58% of the multi-exon genes underwent AS. We further obtained 565 genes with significantly more AS events compared with random genes. These genes were found significant enrichment in biological processes related to the regulation of actin filament length. Second, significantly different AS event profiles among ten tissues were found. The tissues with the same origin of development are more likely to have a relatively similar AS profile. Moreover, 7370 genes showed tissue-specific AS events and were highly enriched in biological processes related to the positive regulation of cellular component organization. Root-specificity AS genes were related to the cellular response to DNA damage stimulus. Third, the genes with different intron retention (IR) patterns among the ten tissues showed significant difference in GC percentages of the retained intron, and the number of exons and FPKM of the major transcripts.

CONCLUSIONS

Our study provided a comprehensive view of AS in multiple tissues. We revealed novel insights into the patterns of AS in multiple tissues and the tissue-specific AS in cucumber.

A global survey of alternative splicing in allopolyploid cotton: landscape, complexity and regulation

Alternative splicing (AS) is a crucial regulatory mechanism in eukaryotes, which acts by greatly increasing transcriptome diversity. The extent and complexity of AS has been revealed in model plants using high-throughput next-generation sequencing. However, this technique is less effective in accurately identifying transcript isoforms in polyploid species because of the high sequence similarity between coexisting subgenomes. Here we characterize AS in the polyploid species cotton. Using Pacific Biosciences single-molecule long-read isoform sequencing (Iso-Seq), we developed an integrated pipeline for Iso-Seq transcriptome data analysis (https://github.com/Nextomics/pipeline-for-isoseq). We identified 176 849 full-length transcript isoforms from 44 968 gene models and updated gene annotation. These data led us to identify 15 102 fibre-specific AS events and estimate that c. 51.4% of homoeologous genes produce divergent isoforms in each subgenome. We reveal that AS allows differential regulation of the same gene by miRNAs at the isoform level. We also show that nucleosome occupancy and DNA methylation play a role in defining exons at the chromatin level. This study provides new insights into the complexity and regulation of AS, and will enhance our understanding of AS in polyploid species. Our methodology for Iso-Seq data analysis will be a useful reference for the study of AS in other species.

Genome-wide characterization of differential transcript usage in Arabidopsis thaliana

Alternative splicing and the usage of alternate transcription start- or stop sites allows a single gene to produce multiple transcript isoforms. Most plant genes express certain isoforms at a significantly higher level than others, but under specific conditions this expression dominance can change, resulting in a different set of dominant isoforms. These events of differential transcript usage (DTU) have been observed for thousands of Arabidopsis thaliana, Zea mays and Vitis vinifera genes, and have been linked to development and stress response. However, neither the characteristics of these genes, nor the implications of DTU on their protein coding sequences or functions, are currently well understood. Here we present a dataset of isoform dominance and DTU for all genes in the AtRTD2 reference transcriptome based on a protocol that was benchmarked on simulated data and validated through comparison with a published reverse transciptase-polymerase chain reaction panel. We report DTU events for 8148 genes across 206 public RNA-Seq samples, and find that protein sequences are affected in 22% of the cases. The observed DTU events show high consistency across replicates, and reveal reproducible patterns in response to treatment and development. We also demonstrate that genes with different evolutionary ages, expression breadths and functions show large differences in the frequency at which they undergo DTU, and in the effect that these events have on their protein sequences. Finally, we showcase how the generated dataset can be used to explore DTU events for genes of interest or to find genes with specific DTU in samples of interest.

iTRAQ-based quantitative proteomic analysis in vernalization-treated faba bean (Vicia faba L.)

Vernalization is classically defined as the induction of flowering process by exposure of the plants to a prolonged cold condition. Normally, it is considered as a precondition of flowering. Vicia faba, commonly known as faba bean, belongs to family Fabaceae. It is one of the plant species that has been cultivated in the earliest human settlements. In this study, an iTRAQ-LC-MS/MS-based quantitative proteomic analysis has been conducted to compare the vernalized faba bean seedlings and its corresponding control. In total, 91 proteins from various functional categories were observed to be differentially accumulated in vernalized faba bean seedlings. Subsequent gene ontology analysis indicated that several biological processes or metabolic pathways including photosynthesis and phytic acid metabolism were differentially respond to vernalization in comparison to the control sample. Further investigation revealed that a family of proteins nominated as glycine-rich RNA-binding factor was accumulated in vernalized seedlings, indicating an extra layer of regulation by alternative splicing on transcript abundance in response to vernalization. These findings raise a possibility that these candidate proteins could be important to represent the responsive network under vernalization process. Therefore, we propose that the regulation of vernalization in faba bean not only occurs at the transcriptional level as previously reported, but also at the post-transcriptional level.

Multilevel Regulation of Abiotic Stress Responses in Plants

The sessile lifestyle of plants requires them to cope with stresses in situ. Plants overcome abiotic stresses by altering structure/morphology, and in some extreme conditions, by compressing the life cycle to survive the stresses in the form of seeds. Genetic and molecular studies have uncovered complex regulatory processes that coordinate stress adaptation and tolerance in plants, which are integrated at various levels. Investigating natural variation in stress responses has provided important insights into the evolutionary processes that shape the integrated regulation of adaptation and tolerance. This review primarily focuses on the current understanding of how transcriptional, post-transcriptional, post-translational, and epigenetic processes along with genetic variation orchestrate stress responses in plants. We also discuss the current and future development of computational tools to identify biologically meaningful factors from high dimensional, genome-scale data and construct the signaling networks consisting of these components.

Alternative splicing complexity contributes to genetic improvement of drought resistance in the rice maintainer HuHan2B

Water-saving and drought-resistantce rice (WDR) breeding practices have greatly increased grain yield and drought resistance. To study the genetic basis of adaptation to drought, transcriptome sequences from the WDR maintainer line HuHan2B and the recurrent parent HanFengB were analyzed for alternative splicing (AS) complexity. Intron retention, the dominant AS type, accounted for 42% of the observed AS events. Differential expression analysis revealed transcripts were preferentially expressed in different varieties and conditions. Based on gene ontology predictions, the biological functions of drought-induced transcripts were significantly enriched in genes involved in transcription regulation, chloroplast components and response to abiotic stimulus in HuHan2B, whereas developmental processes for reproduction were primarily enriched in HanFengB. The regulatory network of transcription factors was driven by cohorts of transcript splicing targets, resulting in more diversified regulatory relationships due to AS complexity than in our previous findings. Moreover, several genes were validated to accumulate novel splicing transcripts in a drought-induced manner. Together, these results suggest that HuHan2B and HanFengB share similar AS features but that a subset of genes with increased levels of AS involved in transcription regulatory networks may contribute an additional level of control for genetic improvement of drought resistance in rice maintainer HuHan2B through breeding.

Differential expression by chromatin modifications of alcohol dehydrogenase 1 of Chorispora bungeana in cold stress

Epigenetic modifications regulate plant genes to cope with a variety of environmental stresses. Chorispora bungeana is an alpine subnival plant with strong tolerance to multiple abiotic stresses, especially cold stress. In this study, we characterized the alcohol dehydrogenase 1 gene from Chorispora bungeana, CbADH1, that is up-regulated in cold conditions. Overexpression of CbADH1 in Arabidopsis thaliana improved cold tolerance, as indicated by a decreased lethal temperature (LT50). Chromatin immunoprecipitation assays showed that histone H3 is removed from the promoter region and the middle-coding region of the gene. H3K9 acetylation and H3K4 trimethylation increased throughout the gene and in the proximal promoter region, respectively. Moreover, increased Ser5P and Ser2P polymerase II accumulation further indicated changes in the transcription initiation and elongation of CbADH1 were due to the cold stress. Taken together, our results suggested that CbADH1 is highly expressed during cold stress, and is regulated by epigenetic modifications. This study expands our understanding of the regulation of gene expression by epigenetic modifications in response to environmental cues.

Domestication reduces alternative splicing expression variations in sorghum

Domestication is known to strongly reduce genomic diversity through population bottlenecks. The resulting loss of polymorphism has been thoroughly documented in numerous cultivated species. Here we investigate the impact of domestication on the diversity of alternative transcript expressions using RNAseq data obtained on cultivated and wild sorghum accessions (ten accessions for each pool). In that aim, we focus on genes expressing two isoforms in sorghum and estimate the ratio between expression levels of those isoforms in each accession. Noticeably, for a given gene, one isoform can either be overexpressed or underexpressed in some wild accessions, whereas in the cultivated accessions, the balance between the two isoforms of the same gene appears to be much more homogenous. Indeed, we observe in sorghum significantly more variation in isoform expression balance among wild accessions than among domesticated accessions. The possibility exists that the loss of nucleotide diversity due to domestication could affect regulatory elements, controlling transcription or degradation of these isoforms. Impact on the isoform expression balance is discussed. As far as we know, this is the first time that the impact of domestication on transcript isoform balance has been studied at the genomic scale. This could pave the way towards the identification of key domestication genes with finely tuned isoform expressions in domesticated accessions while being highly variable in their wild relatives.

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.

Proteogenomic analysis reveals alternative splicing and translation as part of the abscisic acid response in Arabidopsis seedlings

In eukaryotes, mechanisms such as alternative splicing (AS) and alternative translation initiation (ATI) contribute to organismal protein diversity. Specifically, splicing factors play crucial roles in responses to environment and development cues; however, the underlying mechanisms are not well investigated in plants. Here, we report the parallel employment of short-read RNA sequencing, single molecule long-read sequencing and proteomic identification to unravel AS isoforms and previously unannotated proteins in response to abscisic acid (ABA) treatment. Combining the data from the two sequencing methods, approximately 83.4% of intron-containing genes were alternatively spliced. Two AS types, which are referred to as alternative first exon (AFE) and alternative last exon (ALE), were more abundant than intron retention (IR); however, by contrast to AS events detected under normal conditions, differentially expressed AS isoforms were more likely to be translated. ABA extensively affects the AS pattern, indicated by the increasing number of non-conventional splicing sites. This work also identified thousands of unannotated peptides and proteins by ATI based on mass spectrometry and a virtual peptide library deduced from both strands of coding regions within the Arabidopsis genome. The results enhance our understanding of AS and alternative translation mechanisms under normal conditions, and in response to ABA treatment.

A systems approach to a spatio-temporal understanding of the drought stress response in maize

Crops are often subjected to periods of drought stress during their life cycle. However, how stress response mechanisms contribute to the crosstalk between stress signaling pathways and developmental signaling pathways is still unknown. We built a gene co-expression network from a spatio-temporal transcriptomic map of the drought stress response in maize (Zea mays), profiled from three tissues and four developmental stages and characterized hub genes associated with duplication events, selection, and regulatory networks. Co-expression analysis grouped drought-response genes into ten modules, covering 844 highly connected genes (hub genes). Of these, 15.4% hub genes had diverged by whole-genome duplication events and 2.5% might then have been selected during natural domestication and artificial improvement processes, successively. We identified key transcription factor hubs in a transcriptional regulatory network, which may function as a crosstalk mechanism between drought stress and developmental signalling pathways in maize. Understanding the evolutionary biases that have evolved to enhance drought adaptation lays the foundation for further dissection of crosstalk between stress signalling pathways and developmental signalling pathways in maize, towards molecular design of new cultivars with desirable yield and greater stress tolerance.

Genome-wide association mapping of latex yield and girth in Amazonian accessions of Hevea brasiliensis grown in a suboptimal climate zone

Latex yield and growth are the key complex traits in commercial rubber production. The present study is the first to report genome-wide association mapping of latex yield and girth, for 170 Amazonian accessions grown in a suboptimal area characterized by limited rainfall and a lengthy dry season. Targeted sequence enrichment to capture gene transcripts generated 14,155 high quality filtered single nucleotide polymorphisms (SNPs) of which 94.3% resided in coding regions. The rapid decay of linkage disequilibrium over physical and genetic distance found in the accessions was comparable to those previously reported for several outcrossing species. A mixed linear model detected three significant SNPs in three candidate genes involved in plant adaptation to drought stress, individually explaining 12.7-15.7% of the phenotypic variance. The SNPs identified in the study will help to extend understanding, and to support genetic improvement of rubber trees grown in drought-affected regions.

Alternative splicing shapes transcriptome but not proteome diversity in Physcomitrella patens

Alternative splicing (AS) can significantly impact the transcriptome and proteome of a eukaryotic cell. Here, using transcriptome and proteome profiling data, we analyzed AS in two life forms of the model moss Physcomitrella patens, namely protonemata and gametophores, as well as in protoplasts. We identified 12 043 genes subject to alternative splicing and analyzed the extent to which AS contributes to proteome diversity. We could distinguish a few examples that unambiguously indicated the presence of two or more splice isoforms from the same locus at the proteomic level. Our results indicate that alternative isoforms have a small effect on proteome diversity. We also revealed that mRNAs and pre-mRNAs have thousands of complementary binding sites for long non-coding RNAs (lncRNAs) that may lead to potential interactions in transcriptome. This finding points to an additional level of gene expression and AS regulation by non-coding transcripts in Physcomitrella patens. Among the differentially expressed and spliced genes we found serine/arginine-rich (SR) genes, which are known to regulate AS in cells. We found that treatment with abscisic (ABA) and methyl jasmonic acids (MeJA) led to an isoform-specific response and suggested that ABA in gametophores and MeJA in protoplasts regulate AS and the transcription of SR genes.

An event of alternative splicing affects the expression of the NTRC gene, encoding NADPH-thioredoxin reductase C, in seed plants

The NTRC gene encodes a NADPH-dependent thioredoxin reductase with a joint thioredoxin domain, exclusive of photosynthetic organisms. An updated search shows that although most species harbor a single copy of the NTRC gene, two copies were identified in different species of the genus Solanum, Glycine max and the moss Physcomitrella patens. The phylogenetic analysis of NTRCs from different sources produced a tree with the major groups of photosynthetic organisms: cyanobacteria, algae and land plants, indicating the evolutionary success of the NTRC gene among photosynthetic eukaryotes. An event of alternative splicing affecting the expression of the NTRC gene was identified, which is conserved in seed plants but not in algae, bryophytes and lycophytes. The alternative splicing event results in a transcript with premature stop codon, which would produce a truncated form of the enzyme. The standard splicing/alternative splicing (SS/AS) transcripts ratio was higher in photosynthetic tissues from Arabidopsis, Brachypodium and tomato, in line with the higher content of the NTRC polypeptide in these tissues. Moreover, environmental stresses such as cold or high salt affected the SS/AS ratio of the NTRC gene transcripts in Brachypodium seedlings. These results suggest that the alternative splicing of the NTRC gene might be an additional mechanism for modulating the content of NTRC in photosynthetic and non-photosynthetic tissues of seed plants.

Global identification of alternative splicing via comparative analysis of SMRT- and Illumina-based RNA-seq in strawberry

Alternative splicing (AS) is a key post-transcriptional regulatory mechanism, yet little information is known about its roles in fruit crops. Here, AS was globally analyzed in the wild strawberry Fragaria vesca genome with RNA-seq data derived from different stages of fruit development. The AS landscape was characterized and compared between the single-molecule, real-time (SMRT) and Illumina RNA-seq platform. While SMRT has a lower sequencing depth, it identifies more genes undergoing AS (57.67% of detected multiexon genes) when it is compared with Illumina (33.48%), illustrating the efficacy of SMRT in AS identification. We investigated different modes of AS in the context of fruit development; the percentage of intron retention (IR) is markedly reduced whereas that of alternative acceptor sites (AA) is significantly increased post-fertilization when compared with pre-fertilization. When all the identified transcripts were combined, a total of 66.43% detected multiexon genes in strawberry undergo AS, some of which lead to a gain or loss of conserved domains in the gene products. The work demonstrates that SMRT sequencing is highly powerful in AS discovery and provides a rich data resource for later functional studies of different isoforms. Further, shifting AS modes may contribute to rapid changes of gene expression during fruit set.

RNA-seq Analysis of Cold and Drought Responsive Transcriptomes of Zea mays ssp. mexicana L

The annual Zea mays ssp. mexicana L. is a member of teosinte, a wild relative of the Zea mays spp. mays L. This subspecies has strong growth and regeneration ability, high tiller numbers, high protein and lysine content as well as resistance to many fungal diseases, and it can be effectively used in maize improvement. In this study, we reported a Zea mays ssp. mexicana L. transcriptome by merging data from untreated control (CK), cold (4°C) and drought (PEG2000, 20%) treated plant samples. A total of 251,145 transcripts (N50 = 1,269 bp) and 184,280 unigenes (N50 = 923 bp) were predicted, which code for homologs of near 47% of the published maize proteome. Under cold conditions, 2,232 and 817 genes were up-regulated and down-regulated, respectively, while fewer genes were up-regulated (532) and down-regulated (82) under drought stress, indicating that Zea mays ssp. mexicana L. is more sensitive to the applied cold rather than to the applied drought stresses. Functional enrichment analyses identified many common or specific biological processes and gene sets in response to drought and cold stresses. The ABA dependent pathway, trehalose synthetic pathway and the ICE1-CBF pathway were up-regulated by both stresses. GA associated genes have been shown to differentially regulate the responses to cold in close subspecies in Zea mays. These findings and the identified functional genes can provide useful clues for improving abiotic stress tolerance of maize.