Alternative splicing of transcription factors in plant responses to low temperature stress: mechanisms and functions

SMRT sequencing of a full-length transcriptome reveals cold induced alternative splicing in Vitis amurensis root

Alternative splicing enhances diversity at the transcriptional and protein levels that widely involved in plant response to biotic and abiotic stresses. V. amurensis is an extremely cold-tolerant wild grape variety, however, studies on alternative splicing (AS) in amur grape at low temperatures are currently poorly understood. In this study, we analyzed full-length transcriptome and RNA seq data at 0, 2, and 24 h after cold stress in V. amurensis roots. Following quality control and correction, 221,170 high-quality full-length non-concatemer (FLNC) reads were identified. A total of 16,181 loci and 30,733 isoforms were identified. These included 22,868 novel isoforms from annotated genes and 2815 isoforms from 2389 novel genes. Among the distinguished novel isoforms, 673 Long non-coding RNAs (LncRNAs) and 18,164 novel isoforms open reading frame (ORF) region were found. A total of 2958 genes produced 8797 AS events, of which 189 genes were involved in the low-temperature response. Twelve transcription factors show AS during cold treatment and VaMYB108 was selected for initial exploration. Two transcripts, Chr05.63.1 (VaMYB108short) and Chr05.63.2 (VaMYB108normal) of VaMYB108, display up-regulated expression after cold treatment in amur grape roots and are both localized in the nucleus. Only VaMYB108normal exhibits transcriptional activation activity. Overexpression of either VaMYB108short or VaMYB108normal in grape roots leads to increased expression of the other transcript and both increased chilling resistance of amur grape roots. The results improve and supplement the genome annotations and provide insights for further investigation into AS mechanisms during cold stress in V. amurensis.

Genome-wide characterization of post-transcriptional processes related to wood formation in Dalbergia odorifera

BACKGROUND

Alternative polyadenylation (APA), alternative splicing (AS), and long non-coding RNAs (lncRNAs) play regulatory roles in post-transcriptional processes in plants. However, little is known about their involvement in xylem development in Dalbergia odorifera, a valuable rosewood species with medicinal and commercial significance. We addressed this by conducting Isoform Sequencing (Iso-Seq) using PacBio's SMRT technology and combined it with RNA-seq analysis (RNA sequencing on Illumina platform) after collecting xylem samples from the transition zone and the sapwood of D. odorifera.

RESULTS

We identified 14,938 full-length transcripts, including 9,830 novel isoforms, which has updated the D. odorifera genome annotation. Our analysis has revealed that 4,164 genes undergo APA, whereas 3,084 genes encounter AS. We have also annotated 118 lncRNAs. Furthermore, RNA-seq analysis identified 170 differential alternative splicing (DAS) events, 344 genes with differential APA site usage (DE-APA), and 6 differentially expressed lncRNAs in the transition zone when compared to the sapwood. AS, APA, and lncRNAs are differentially regulated during xylem development. Differentially expressed APA genes were enriched for terpenoid and flavonoid metabolism, indicating their role in the heartwood formation. Additionally, DE-APA genes were associated with cell wall biosynthesis and terpenoid metabolism, implying an APA's role in wood formation. A DAS gene (involved in chalcone accumulation) with a significantly greater inclusion of the last exon in the transition zone than in the sapwood was identified. We also found that differentially expressed lncRNAs targeted the genes related to terpene synthesis.

CONCLUSIONS

This study enhances our understanding of the molecular regulatory mechanisms underlying wood formation in D. odorifera, and provides valuable genetic resources and insights for its molecular-assisted breeding.

Genomics for monitoring and understanding species responses to global climate change

All life forms across the globe are experiencing drastic changes in environmental conditions as a result of global climate change. These environmental changes are happening rapidly, incur substantial socioeconomic costs, pose threats to biodiversity and diminish a species' potential to adapt to future environments. Understanding and monitoring how organisms respond to human-driven climate change is therefore a major priority for the conservation of biodiversity in a rapidly changing environment. Recent developments in genomic, transcriptomic and epigenomic technologies are enabling unprecedented insights into the evolutionary processes and molecular bases of adaptation. This Review summarizes methods that apply and integrate omics tools to experimentally investigate, monitor and predict how species and communities in the wild cope with global climate change, which is by genetically adapting to new environmental conditions, through range shifts or through phenotypic plasticity. We identify advantages and limitations of each method and discuss future research avenues that would improve our understanding of species' evolutionary responses to global climate change, highlighting the need for holistic, multi-omics approaches to ecosystem monitoring during global climate change.

Regulatory mechanism of the miR172e-LbrAP2 module during the vegetative growth phase transition in Lilium

MAIN CONCLUSION

It was proved for the first time that the miR172e-LbrAP2 module regulated the vegetative growth phase transition in Lilium, which provided a new approach to shorten the juvenile stage of Lilium, improved the reproduction rate, and reduced the propagation cost of Lilium commercial bulbs. Lilium is an ornamental bulb plant that takes at least 3 years to cultivate into commercial seed bulbs under natural conditions. The aim of this study was to shorten the Lilium expansion cycle. In this study, the growth cycle of lily tubers induced by low temperature of 15 °C was significantly shorter than that of tubers grown at a conventional temperature. Quantitative real-time PCR analysis showed that the expression patterns of miR172e and LbrAP2 were negatively correlated. GUS histochemical staining confirmed that miR172e and LbrAP2 in tobacco leaves interacted with each other after co-transformation. The shear sites of miR172e and its target gene, LbrAP2, upon binding, were identified by RLM 5' RACE analysis. In addition, miR172e and LbrAP2 showed opposite expression patterns after the transformation of Arabidopsis. miR172e overexpression accelerated the transition from juvenile to adult plants, whereas LbrAP2 overexpression inhibited this process, thus indicating that miR172e negatively regulated the target gene LbrAP2. Upregulation of the transcription factor LbrAP2 delayed the phase transition of plants, whereas miR172 inhibited the transcriptional translation of LbrAP2, thereby accelerating the phase transition. Low-temperature treatment of Lilium bulbs can shorten Lilium development, which provides a new approach to accelerating Lilium commercial bulb breeding and reducing breeding costs.

Alternative Splicing under Cold Stress in Paper Mulberry

The paper mulberry is a commonly found tree species with a long history of cultivation. It also serves as a crucial case study for understanding how woody plants adapt to low temperatures. Under cold treatment, we observed a substantial number of alternative splicing (AS) genes, showcasing the intricate landscape of AS events. We have detected all seven types of AS events, with the alternative 3' splice site (A3) having the most. We observed that many genes that underwent differential AS were significantly enriched in starch and sucrose metabolism and circadian rhythm pathways. Moreover, a considerable proportion of differentially spliced genes (DSGs) also showed differential expression, with 20.38% and 25.65% under 12 h and 24 h cold treatments, respectively. This suggests a coordinated regulation between gene AS and expression, playing a pivotal role in the paper mulberry's adaptation to cold stress. We further investigated the regulatory mechanisms of AS, identifying 41 serine/arginine-rich (SR) splicing factors, among which 11 showed differential expression under cold treatment, while 29 underwent alternative splicing. Additionally, genes undergoing AS displayed significantly higher DNA methylation levels under cold stress, while normal splicing (non-AS) genes exhibited relatively lower methylation levels. These findings suggest that methylation may play an important role in governing gene AS. Finally, our research will provide useful information on the role of AS in the cold acclimation tolerance of the paper mulberry.

The Genome-Wide Characterization of Alternative Splicing and RNA Editing in the Development of Coprinopsis cinerea

Coprinopsis cinerea is one of the model species used in fungal developmental studies. This mushroom-forming Basidiomycetes fungus has several developmental destinies in response to changing environments, with dynamic developmental regulations of the organism. Although the gene expression in C. cinerea development has already been profiled broadly, previous studies have only focused on a specific stage or process of fungal development. A comprehensive perspective across different developmental paths is lacking, and a global view on the dynamic transcriptional regulations in the life cycle and the developmental paths is far from complete. In addition, knowledge on co- and post-transcriptional modifications in this fungus remains rare. In this study, we investigated the transcriptional changes and modifications in C. cinerea during the processes of spore germination, vegetative growth, oidiation, sclerotia formation, and fruiting body formation by inducing different developmental paths of the organism and profiling the transcriptomes using the high-throughput sequencing method. Transition in the identity and abundance of expressed genes drive the physiological and morphological alterations of the organism, including metabolism and multicellularity construction. Moreover, stage- and tissue-specific alternative splicing and RNA editing took place and functioned in C. cinerea. These modifications were negatively correlated to the conservation features of genes and could provide extra plasticity to the transcriptome during fungal development. We suggest that C. cinerea applies different molecular strategies in its developmental regulation, including shifts in expressed gene sets, diversifications of genetic information, and reversible diversifications of RNA molecules. Such features would increase the fungal adaptability in the rapidly changing environment, especially in the transition of developmental programs and the maintenance and balance of genetic and transcriptomic divergence. The multi-layer regulatory network of gene expression serves as the molecular basis of the functioning of developmental regulation.

Mapping intron retention events contributing to complex traits using splice quantitative trait locus

BACKGROUND

Alternative splicing (AS) of mRNA plays an important roles in transcriptome diversity, involving regulation of plant growth and stress response. Understanding the variation of AS events underlying GWAS loci in a crop population can provide insight into the molecular mechanisms of complex agronomic traits. To date, genome-wide association studies relating AS events to agronomic traits have rarely been conducted at the population level in crops.

RESULTS

Here, a pipeline was constructed to identify candidate AS events related to complex traits. Firstly, ovule transcriptome data were used to characterize intron retention (IR), the predominant type of AS in plants, on a genome-wide scale. This was done in a natural population consisting of 279 upland cotton lines. Secondly, splice quantitative trait locus (sQTL) analysis was carried out, which yielded a total of 2295 sQTLs involving 1607 genes. Of these, 14.25% (n = 427) were cis-regulatory loci. Integration with expression quantitative trait loci (eQTL) revealed that 53 (21.4%) cis-sGenes were regulated by both cis-sQTLs and cis-eQTLs. Finally, co-localization analysis integrated with GWAS loci in this population showed 32 cis-QTLs to be co-located with genetic regulatory loci related to fiber yield and quality traits, indicating that sQTLs are likely to participate in regulating cotton fiber yield and quality. An in-depth evaluation confirmed that differences in the IR rates of sQTL-regulated candidate genes such as GhLRRK1 and GhGC1 are associated with lint percentage (LP), which has potential in breeding applications.

CONCLUSION

This study provides a clue that AS of mRNA has an impact on crop yield, along with functional sQTLs are new genetic resources for cotton precision breeding.

Alternative splicing shapes the transcriptome complexity in blackgram [Vigna mungo (L.) Hepper]

Vigna mungo, a highly consumed crop in the pan-Asian countries, is vulnerable to several biotic and abiotic stresses. Understanding the post-transcriptional gene regulatory cascades, especially alternative splicing (AS), may underpin large-scale genetic improvements to develop stress-resilient varieties. Herein, a transcriptome based approach was undertaken to decipher the genome-wide AS landscape and splicing dynamics in order to establish the intricacies of their functional interactions in various tissues and stresses. RNA sequencing followed by high-throughput computational analyses identified 54,526 AS events involving 15,506 AS genes that generated 57,405 transcripts isoforms. Enrichment analysis revealed their involvement in diverse regulatory functions and demonstrated that transcription factors are splicing-intensive, splice variants of which are expressed differentially across tissues and environmental cues. Increased expression of a splicing regulator NHP2L1/SNU13 was found to co-occur with lower intron retention events. The host transcriptome is significantly impacted by differential isoform expression of 1172 and 765 AS genes that resulted in 1227 (46.8% up and 53.2% downregulated) and 831 (47.5% up and 52.5% downregulated) transcript isoforms under viral pathogenesis and Fe²⁺ stressed condition, respectively. However, genes experiencing AS operate differently from the differentially expressed genes, suggesting AS is a unique and independent mode of regulatory mechanism. Therefore, it can be inferred that AS mediates a crucial regulatory role across tissues and stressful situations and the results would provide an invaluable resource for future endeavours in V. mungo genomics.

SWATH-MS-based proteogenomic analysis reveals the involvement of alternative splicing in poplar upon lead stress

Alternative splicing (AS) regulates gene expression and increases proteomic diversity for the fine tuning of stress responses in plants, but the exact mechanism through which AS functions in plant stress responses is not thoroughly understood. Here, we investigated how AS functions in poplar (Populus trichocarpa), a popular plant for bioremediation, in response to lead (Pb) stress. Using a proteogenomic analysis, we determine that Pb stress induced alterations in AS patterns that are characterized by an increased use of nonconventional splice sites and a higher abundance of Pb-responsive splicing factors (SFs) associated with Pb-responsive transcription factors. A strong Pb(II)-inducible chaperone protein, PtHSP70, that undergoes AS was further characterized. Overexpression of its two spliced isoforms, PtHSP70-AS1 and PtHSP70-AS2, in poplar and Arabidopsis significantly enhances the tolerance to Pb. Further characterization shows that both isoforms can directly bind to Pb(II), and PtHSP70-AS2 exhibits 10-fold higher binding capacities and a greater increase in expression under Pb stress, thereby reducing cellular toxicity through Pb(II) extrusion and conferring Pb tolerance. AS of PtHSP70 is found to be regulated by PtU1-70K, a Pb(II)-inducible core SF involved in 5'-splice site recognition. Because the same splicing pattern is also found in HSP70 orthologs in other plant species, AS of HSP70 may be a common regulatory mechanism to cope with Pb(II) toxicity. Overall, we have revealed a novel post-transcriptional machinery that mediates heavy metal tolerance in diverse plant species. Our findings offer new molecular targets and bioengineering strategies for phytoremediation and provide new insight for future directions in AS research.

Alternative Splicing in the Regulatory Circuit of Plant Temperature Response

As sessile organisms, plants have evolved complex mechanisms to rapidly respond to ever-changing ambient temperatures. Temperature response in plants is modulated by a multilayer regulatory network, including transcriptional and post-transcriptional regulations. Alternative splicing (AS) is an essential post-transcriptional regulatory mechanism. Extensive studies have confirmed its key role in plant temperature response, from adjustment to diurnal and seasonal temperature changes to response to extreme temperatures, which has been well documented by previous reviews. As a key node in the temperature response regulatory network, AS can be modulated by various upstream regulations, such as chromatin modification, transcription rate, RNA binding proteins, RNA structure and RNA modifications. Meanwhile, a number of downstream mechanisms are affected by AS, such as nonsense-mediated mRNA decay (NMD) pathway, translation efficiency and production of different protein variants. In this review, we focus on the links between splicing regulation and other mechanisms in plant temperature response. Recent advances regarding how AS is regulated and the following consequences in gene functional modulation in plant temperature response will be discussed. Substantial evidence suggests that a multilayer regulatory network integrating AS in plant temperature response has been unveiled.

Unleashing the power within short-read RNA-seq for plant research: Beyond differential expression analysis and toward regulomics

RNA-seq has become a state-of-the-art technique for transcriptomic studies. Advances in both RNA-seq techniques and the corresponding analysis tools and pipelines have unprecedently shaped our understanding in almost every aspects of plant sciences. Notably, the integration of huge amount of RNA-seq with other omic data sets in the model plants and major crop species have facilitated plant regulomics, while the RNA-seq analysis has still been primarily used for differential expression analysis in many less-studied plant species. To unleash the analytical power of RNA-seq in plant species, especially less-studied species and biomass crops, we summarize recent achievements of RNA-seq analysis in the major plant species and representative tools in the four types of application: (1) transcriptome assembly, (2) construction of expression atlas, (3) network analysis, and (4) structural alteration. We emphasize the importance of expression atlas, coexpression networks and predictions of gene regulatory relationships in moving plant transcriptomes toward regulomics, an omic view of genome-wide transcription regulation. We highlight what can be achieved in plant research with RNA-seq by introducing a list of representative RNA-seq analysis tools and resources that are developed for certain minor species or suitable for the analysis without species limitation. In summary, we provide an updated digest on RNA-seq tools, resources and the diverse applications for plant research, and our perspective on the power and challenges of short-read RNA-seq analysis from a regulomic point view. A full utilization of these fruitful RNA-seq resources will promote plant omic research to a higher level, especially in those less studied species.

Genome-Wide Analysis of the Almond AP2/ERF Superfamily and Its Functional Prediction during Dormancy in Response to Freezing Stress

Simple Summary

The ethylene-responsive element (AP2/ERF) is one of the key and conserved transcription factors (TFs) in plants, and it plays a crucial role in regulating plant growth, development, and stress response. The cultivated almond in Xinjiang is often affected by short-term ultralow temperature freezing stress during the winter dormancy period, resulting in the death of large-scale almond plants. In this study, we conducted the first genome-wide analysis of the PdAP2/ERF family in almond, including protein physicochemical properties, phylogenetic relationships, motif types, gene structures, gene replication types, collinearity relationships, and cis-element types in promoter regions. We further analyzed the expression patterns of the PdAP2/ERF gene in different tissues of almond and under freezing stress at different temperatures in annual dormant branches using transcriptome data. In addition, we also analyzed the expression levels of 13 PdAP2/ERF genes in four tissues of almond and in annual dormant branches treated with freezing stress at different temperatures using fluorescence quantitative technology. This study laid the foundation for further exploring the function of the PdAP2/ERF gene in almond.

Abstract

The AP2/ERF transcription factor family is one of the largest transcription factor families in plants and plays an important role in regulating plant growth and development and the response to biotic and abiotic stresses. However, there is no report on the AP2/ERF gene family in almond (Prunus dulcis). In this study, a total of 136 PdAP2/ERF genes were identified from the almond genome, and their protein physicochemical properties were analyzed. The PdAP2/ERF members were divided into five subgroups: AP2, RAV, ERF, DREB, and Soloist. The PdAP2/ERF members in each subgroup had conserved motif types and exon/intron numbers. PdAP2/ERFS members are distributed on eight chromosomes, with 22 pairs of segmental duplications and 28 pairs of tandem duplications. We further explored the colinear relationship between almond and Arabidopsis thaliana, Oryza sativa, Malus domestica, and Prunus persicaAP2/ERF genes and their evolution. The results of cis-acting elements showed that PdAP2/ERF members are widely involved in various processes, such as growth and development, hormone regulation, and stress response. The results based on transcriptome expression patterns showed that PdAP2/ERF genes had significant tissue-specific expression characteristics and were involved in the response of annual dormant branches of almond to low-temperature freezing stress. In addition, the fluorescence quantitative relative expression results of 13 representative PdAP2/ERF genes in four tissues of ‘Wanfeng’ almond and under six low-temperature freezing treatments of annual dormant branches were consistent with the transcriptome results. It is worth noting that the fluorescence quantitative expression level showed that the PdERF24 gene was extremely significant at −30 °C, suggesting that this gene may play an important role in the response of almond dormancy to ultralow temperature freezing stress. Finally, we identified 7424 and 6971 target genes based on AP2 and ERF/DREB DNA-binding sites, respectively. The GO and KEGG enrichment results showed that these target genes play important roles in protein function and multiple pathways. In summary, we conducted bioinformatics and expression pattern studies on PdAP2/ERF genes, including 13 PdAP2/ERF genes, and performed fluorescence quantitative analysis of annual dormant shoots under different low-temperature freezing stress treatments to understand the tolerance of almond dormancy to freezing stress and suggest future improvements.

Alternative Splicing and Its Roles in Plant Metabolism

Plant metabolism, including primary metabolism such as tricarboxylic acid cycle, glycolysis, shikimate and amino acid pathways as well as specialized metabolism such as biosynthesis of phenolics, alkaloids and saponins, contributes to plant survival, growth, development and interactions with the environment. To this end, these metabolic processes are tightly and finely regulated transcriptionally, post-transcriptionally, translationally and post-translationally in response to different growth and developmental stages as well as the constantly changing environment. In this review, we summarize and describe the current knowledge of the regulation of plant metabolism by alternative splicing, a post-transcriptional regulatory mechanism that generates multiple protein isoforms from a single gene by using alternative splice sites during splicing. Numerous genes in plant metabolism have been shown to be alternatively spliced under different developmental stages and stress conditions. In particular, alternative splicing serves as a regulatory mechanism to fine-tune plant metabolism by altering biochemical activities, interaction and subcellular localization of proteins encoded by splice isoforms of various genes.

The Importance of a Genome-Wide Association Analysis in the Study of Alternative Splicing Mutations in Plants with a Special Focus on Maize

Alternative splicing is an important mechanism for regulating gene expressions at the post-transcriptional level. In eukaryotes, the genes are transcribed in the nucleus to produce pre-mRNAs and alternative splicing can splice a pre-mRNA to eventually form multiple different mature mRNAs, greatly increasing the number of genes and protein diversity. Alternative splicing is involved in the regulation of various plant life activities, especially the response of plants to abiotic stresses and is also an important process of plant growth and development. This review aims to clarify the usefulness of a genome-wide association analysis in the study of alternatively spliced variants by summarizing the application of alternative splicing, genome-wide association analyses and genome-wide association analyses in alternative splicing, as well as summarizing the related research progress.

Surviving and thriving: How plants perceive and respond to temperature stress

The dramatic temperature fluctuations spurred by climate change inhibit plant growth and threaten crop productivity. Unraveling how plants defend themselves against temperature-stress-induced cellular impairment is not only a crucial fundamental issue but is also of critical importance for agricultural sustainability and food security. Here, we review recent developments in elucidating the molecular mechanisms used by plants to sense and respond to cold and heat stress at multiple levels. We also describe the trade-off between plant growth and responses to high and low temperatures. Finally, we discuss possible strategies that could be used to engineer temperature-stress-tolerant, high-yielding crops.

ABA Mediates Plant Development and Abiotic Stress via Alternative Splicing

Alternative splicing (AS) exists in eukaryotes to increase the complexity and adaptability of systems under biophysiological conditions by increasing transcriptional and protein diversity. As a classic hormone, abscisic acid (ABA) can effectively control plant growth, improve stress resistance, and promote dormancy. At the transcriptional level, ABA helps plants respond to the outside world by regulating transcription factors through signal transduction pathways to regulate gene expression. However, at the post-transcriptional level, the mechanism by which ABA can regulate plant biological processes by mediating alternative splicing is not well understood. Therefore, this paper briefly introduces the mechanism of ABA-induced alternative splicing and the role of ABA mediating AS in plant response to the environment and its own growth.

CFI 25 Subunit of Cleavage Factor I is Important for Maintaining the Diversity of 3' UTR Lengths in Arabidopsis thaliana (L.) Heynh

Cleavage and polyadenylation at the 3' end of the pre-mRNA is essential for mRNA function, by regulating its translatability, stability and translocation to the cytoplasm. Cleavage factor I (CFI) is a multi-subunit component of the pre-mRNA 3' end processing machinery in eukaryotes. Here, we report that plant CFI 25 subunit of CFI plays an important role in maintaining the diversity of the 3' ends of mRNA. The genome of Arabidopsis thaliana (L.) Heynh. contained four genes encoding three putative CFI subunits (AtCFI 25, AtCFI 59 and AtCFI 68), orthologous to the mammalian CFI subunits. There were two CFI 25 paralogs (AtCFI 25a and AtCFI 25b) that shared homology with human CFI 25. Two null alleles of AtCFI 25a displayed smaller rosette leaves, longer stigmatic papilla, smaller anther, earlier flowering and lower fertility compared to wild-type plants. Null alleles of AtCFI 25b, as well as, plants ectopically expressing full-length cDNA of AtCFI 25a, displayed no obvious morphological defects. AtCFI 25a was shown to interact with AtCFI 25b, AtCFI 68 and itself, suggesting various forms of CFI in plants. Furthermore, we show that AtCFI 25a function was essential for maintaining proper diversity of the 3' end lengths of transcripts coding for CFI subunits, suggesting a self-regulation of the CFI machinery in plants. AtCFI 25a was also important to maintain 3' ends for other genes to different extent. Collectively, AtCFI 25a, but not AtCFI 25b, seemed to play important roles during Arabidopsis development by maintaining proper diversity of the 3' UTR lengths.

Transcriptome profiling reveals key genes in regulation of the tepal trichome development in Lilium pumilum D.C

A number of potential genes and pathways involved in tepal trichome development were identified in a natural lily mutant by transcriptome analysis and were confirmed with trichome and trichomeless species. Trichome is a specialized structure found on the surface of the plant with an important function in survival against abiotic and biotic stress. It is also an important economic trait in crop breeding. Extensive research has investigated the foliar trichome in model plants (Arabidopsis and tomato). However, the developmental mechanism of tepal trichome remains elusive. Lilium pumilum is an edible ornamental bulb and a good breeding parent possessing cold and salt-alkali resistance. Here, we found a natural mutant of Lilium pumilum grown on a highland whose tepals are covered by trichomes. Our data indicate that trichomes of the mutant are multicellular and branchless. Notably, stomata are also developed on the tepal of the mutant as well, suggesting there may be a correlation between trichome and stomata regulation. Furthermore, we isolated 27 differentially expressed genes (DEGs) by comparing the transcriptome profiling between the natural mutant and the wild type. These 27 genes belong to 4 groups: epidermal cell cycle and division, trichome morphogenesis, stress response, and transcription factors. Quantitative real-time PCR in Lilium pumilum (natural mutant and the wild type) and other lily species (Lilium leichtlinii var. maximowiczii/trichome; Lilium davidii var. willmottiae/, trichomeless) confirmed the validation of RNA-seq data and identified several trichome-related genes.

PacBio Iso-Seq Improves the Rainbow Trout Genome Annotation and Identifies Alternative Splicing Associated With Economically Important Phenotypes

Rainbow trout is an important model organism that has received concerted international efforts to study the transcriptome. For this purpose, short-read sequencing has been primarily used over the past decade. However, these sequences are too short of resolving the transcriptome complexity. This study reported a first full-length transcriptome assembly of the rainbow trout using single-molecule long-read isoform sequencing (Iso-Seq). Extensive computational approaches were used to refine and validate the reconstructed transcriptome. The study identified 10,640 high-confidence transcripts not previously annotated, in addition to 1,479 isoforms not mapped to the current Swanson reference genome. Most of the identified lncRNAs were non-coding variants of coding transcripts. The majority of genes had multiple transcript isoforms (average ∼3 isoforms/locus). Intron retention (IR) and exon skipping (ES) accounted for 56% of alternative splicing (AS) events. Iso-Seq improved the reference genome annotation, which allowed identification of characteristic AS associated with fish growth, muscle accretion, disease resistance, stress response, and fish migration. For instance, an ES in GVIN1 gene existed in fish susceptible to bacterial cold-water disease (BCWD). Besides, under five stress conditions, there was a commonly regulated exon in prolyl 4-hydroxylase subunit alpha-2 (P4HA2) gene. The reconstructed gene models and their posttranscriptional processing in rainbow trout provide invaluable resources that could be further used for future genetics and genomics studies. Additionally, the study identified characteristic transcription events associated with economically important phenotypes, which could be applied in selective breeding.

Beyond the Genetic Pathways, Flowering Regulation Complexity in Arabidopsis thaliana

Flowering is one of the most critical developmental transitions in plants’ life. The irreversible change from the vegetative to the reproductive stage is strictly controlled to ensure the progeny’s success. In Arabidopsis thaliana, seven flowering genetic pathways have been described under specific growth conditions. However, the evidence condensed here suggest that these pathways are tightly interconnected in a complex multilevel regulatory network. In this review, we pursue an integrative approach emphasizing the molecular interactions among the flowering regulatory network components. We also consider that the same regulatory network prevents or induces flowering phase change in response to internal cues modulated by environmental signals. In this sense, we describe how during the vegetative phase of development it is essential to prevent the expression of flowering promoting genes until they are required. Then, we mention flowering regulation under suboptimal growing temperatures, such as those in autumn and winter. We next expose the requirement of endogenous signals in flowering, and finally, the acceleration of this transition by long-day photoperiod and temperature rise signals allowing A. thaliana to bloom in spring and summer seasons. With this approach, we aim to provide an initial systemic view to help the reader integrate this complex developmental process.

PSDX: A Comprehensive Multi-Omics Association Database of Populus trichocarpa With a Focus on the Secondary Growth in Response to Stresses

Populus trichocarpa (P. trichocarpa) is a model tree for the investigation of wood formation. In recent years, researchers have generated a large number of high-throughput sequencing data in P. trichocarpa. However, no comprehensive database that provides multi-omics associations for the investigation of secondary growth in response to diverse stresses has been reported. Therefore, we developed a public repository that presents comprehensive measurements of gene expression and post-transcriptional regulation by integrating 144 RNA-Seq, 33 ChIP-seq, and six single-molecule real-time (SMRT) isoform sequencing (Iso-seq) libraries prepared from tissues subjected to different stresses. All the samples from different studies were analyzed to obtain gene expression, co-expression network, and differentially expressed genes (DEG) using unified parameters, which allowed comparison of results from different studies and treatments. In addition to gene expression, we also identified and deposited pre-processed data about alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI). The post-transcriptional regulation, differential expression, and co-expression network datasets were integrated into a new P. trichocarpa Stem Differentiating Xylem (PSDX) database (http://forestry.fafu.edu.cn/db/SDX), which further highlights gene families of RNA-binding proteins and stress-related genes. The PSDX also provides tools for data query, visualization, a genome browser, and the BLAST option for sequence-based query. Much of the data is also available for bulk download. The availability of PSDX contributes to the research related to the secondary growth in response to stresses in P. trichocarpa, which will provide new insights that can be useful for the improvement of stress tolerance in woody plants.

Temperature-Dependent Alternative Splicing of Precursor mRNAs and Its Biological Significance: A Review Focused on Post-Transcriptional Regulation of a Cold Shock Protein Gene in Hibernating Mammals

Multiple mRNA isoforms are often generated during processing such as alternative splicing of precursor mRNAs (pre-mRNA), resulting in a diversity of generated proteins. Alternative splicing is an essential mechanism for the functional complexity of eukaryotes. Temperature, which is involved in all life activities at various levels, is one of regulatory factors for controlling patterns of alternative splicing. Temperature-dependent alternative splicing is associated with various phenotypes such as flowering and circadian clock in plants and sex determination in poikilothermic animals. In some specific situations, temperature-dependent alternative splicing can be evoked even in homothermal animals. For example, the splicing pattern of mRNA for a cold shock protein, cold-inducible RNA-binding protein (CIRP or CIRBP), is changed in response to a marked drop in body temperature during hibernation of hamsters. In this review, we describe the current knowledge about mechanisms and functions of temperature-dependent alternative splicing in plants and animals. Then we discuss the physiological significance of hypothermia-induced alternative splicing of a cold shock protein gene in hibernating and non-hibernating animals.

Investigation for a multi-silique trait in Brassica napus by alternative splicing analysis

Background

Flower and fruit development are vital stages of the angiosperm lifecycle. We previously investigated the multi-silique trait in the rapeseed (Brassica napus) line zws-ms on a genomic and transcriptomic level, leading to the identification of two genomic regions and several candidate genes associated with this trait. However, some events on the transcriptome level, like alternative splicing, were poorly understood.

Methods

Plants from zws-ms and its near-isogenic line (NIL) zws-217 were both grown in Xindu with normal conditions and a colder area Ma'erkang. Buds from the two lines were sampled and RNA was isolated to perform the transcriptomic sequencing. The numbers and types of alternative splicing (AS) events from the two lines were counted and classified. Genes with AS events and expressed differentially between the two lines, as well as genes with AS events which occurred in only one line were emphasized. Their annotations were further studied.

Results

From the plants in Xindu District, an average of 205,496 AS events, which could be sorted into five AS types, were identified. zws-ms and zws-217 shared highly similar ratios of each AS type: The alternative 5' and 3' splice site types were the most common, while the exon skipping type was observed least often. Eleven differentially expressed AS genes were identified, of which four were upregulated and seven were downregulated in zws-ms. Their annotations implied that five of these genes were directly associated with the multi-silique trait. While samples from colder area Ma'erkang generated generally reduced number of each type of AS events except for Intron Retention; but the number of differentially expressed AS genes increased significantly. Further analysis found that among the 11 differentially expressed AS genes from Xindu, three of them maintained the same expression models, while the other eight genes did not show significant difference between the two lines in expression level. Additionally, the 205 line-specific expressed AS genes were analyzed, of which 187 could be annotated, and two were considered to be important.

Discussion

This study provides new insights into the molecular mechanism of the agronomically important multi-silique trait in rapeseed on the transcriptome level and screens out some environment-responding candidate genes.

Characterization of novel regulators for heat stress tolerance in tomato from Indian sub-continent

The footprint of tomato cultivation, a cool region crop that exhibits heat stress (HS) sensitivity, is increasing in the tropics/sub-tropics. Knowledge of novel regulatory hot spots from varieties growing in the Indian sub-continent climatic zones could be vital for developing HS-resilient crops. Comparative transcriptome-wide signatures of a tolerant (CLN1621L) and sensitive (CA4) cultivar pair shortlisted from a pool of varieties exhibiting variable thermo-sensitivity using physiological-, survival- and yield-related traits revealed redundant to cultivar-specific HS regulation. The antagonistically expressing genes encode enzymes and proteins that have roles in plant defence and abiotic stresses. Functional characterization of three antagonistic genes by overexpression and silencing established Solyc09g014280 (Acylsugar acyltransferase) and Solyc07g056570 (Notabilis) that are up-regulated in tolerant cultivar, as positive regulators of HS tolerance and Solyc03g020030 (Pin-II proteinase inhibitor), that are down-regulated in CLN1621L, as negative regulator of thermotolerance. Transcriptional assessment of promoters of these genes by SNPs in stress-responsive cis-elements and promoter swapping experiments in opposite cultivar background showed inherent cultivar-specific orchestration of transcription factors in regulating transcription. Moreover, overexpression of three ethylene response transcription factors (ERF.C1/F4/F5) also improved HS tolerance in tomato. This study identifies several novel HS tolerance genes and provides proof of their utility in tomato thermotolerance.

Properties of two spliceoforms of major ampullate spidroin 1 reveal unique functions of N-linker region

Spiders produce diverse silk fibers with distinct properties for daily survival. Among these silk fibers, dragline silk spun by major ampullate gland is used for bridgelines and web radii, exhibiting both outstanding tensile strength and extensibility. Although more and more full-length major ampullate spidroin gene sequences have been reported, the research regarding alternative splicing events of spidroins are rare. Here we describe two spliceoforms of major ampullate spidroin 1 (MaSp1) from Araneus ventricosus, and both of them are lack of central repetitive region. The minor isoform only has terminal regions. For the major isoform, however, the N-linker and terminal regions are all retained. Furthermore, we investigated the functions of N-linker structure of A. ventricosus MaSp1, based on the properties of the two spliceoforms. The dimer level of major isoform (MaSp1-2) is higher than that of the minor isoform (MaSp1-1). Moreover, the MaSp1-2 protein display higher melting temperature (T_m) than MaSp1-1, and the MaSp1-2 fibers exhibit higher tensile strength than MaSp1-1 fibers. These studies demonstrate that the N-linker region promotes the formation of intermolecular disulphide bond, suggesting a strategy to enhance the thermostabilization and mechanical properties of spidroins.

Gene expression and metabolic changes of Momordica charantia L. seedlings in response to low temperature stress

Low temperature is one of the abiotic factors limiting germination, growth and distribution of the plant in current plant-products industry, especially for the tropical vegetables in non-tropical area or other fields under cold temperature. Screening the plant with ability against cold temperature captured worldwide attention and exerted great importance. In our previous work, the anti-cold specie of Momordica Charantia L. seedlings was screened out. Yet, the molecular and physiological mechanisms underlying this adaptive process still remain unknown. This study was aimed to investigate adaption mechanism of anti-cold species of Momordica Charantia L. seedlings in genetical and metabolomics levels. Two species, cold-susceptible group (Y17) and cold-resistant group (Y54), were evaluated containing the indexes of malondialdehyde (MDA), hydrogen peroxide (H2O2), proline content, activities of antioxidant enzymes, metabolites changes and genes differentiation in plant tissues after cold treatment. It found that low temperature stress resulted in increased accumulation of MDA, H2O2 and proline content in two species, but less expressions in cold-resistant species Y54. As compared to Y17, cold-resistant species Y54 presented significantly enhanced antioxidant enzyme activities of POD (peroxidase), CAT (cataalase) and SOD (superoxide dismutase). Meanwhile, higher expressed genes encoded antioxidant enzymes and transcription factors when exposure to the low temperature were found in cold-resistant species Y54, and core genes were explored by Q-PCR validation, including McSOD1, McPDC1 and McCHS1. Moreover, plant metabolites containing amino acid, sugar, fatty acid and organic acid in Y54 were higher than Y17, indicating their important roles in cold acclimation. Meanwhile, initial metabolites, including amimo acids, polypeptides, sugars, organic acids and nucleobases, were apparently increased in cold resistant species Y54 than cold susceptible species Y17. Our results demonstrated that the Momordica Charantia L. seedlings achieved cold tolerance might be went through mobilization of antioxidant systems, adjustment of the transcription factors and accumulation of osmoregulation substance. This work presented meaning information for revealing the anti-cold mechanism of the Momordica Charantia L. seedlings and newsight for further screening of anti-cold species in other plant.

Comprehensive profiling of alternative splicing landscape during cold acclimation in tea plant

BACKGROUND

Alternative splicing (AS) may generate multiple mRNA splicing isoforms from a single mRNA precursor using different splicing sites, leading to enhanced diversity of transcripts and proteins. AS has been implicated in cold acclimation by affecting gene expression in various ways, yet little information is known about how AS influences cold responses in tea plant (Camellia sinensis).

RESULTS

In this study, the AS transcriptional landscape was characterized in the tea plant genome using high-throughput RNA-seq during cold acclimation. We found that more than 41% (14,103) of genes underwent AS events. We summarize the possible existence of 11 types of AS events, including the four common types of intron retention (IR), exon skipping (ES), alternative 5' splice site (A5SS), and alternative 3' splice site (A3SS); of these, IR was the major type in all samples. The number of AS events increased rapidly during cold treatment, but decreased significantly following de-acclimation (DA). It is notable that the number of differential AS genes gradually increased during cold acclimation, and these genes were enriched in pathways relating to oxidoreductase activity and sugar metabolism during acclimation and de-acclimation. Remarkably, the AS isoforms of bHLH transcription factors showed higher expression levels than their full-length ones during cold acclimation. Interestingly, the expression pattern of some AS transcripts of raffinose and sucrose synthase genes were significantly correlated with sugar contents.

CONCLUSION

Our findings demonstrated that changes in AS numbers and transcript expression may contribute to rapid changes in gene expression and metabolite profile during cold acclimation, suggesting that AS events play an important regulatory role in response to cold acclimation in tea plant.

The complexity of alternative splicing and landscape of tissue-specific expression in lotus (Nelumbo nucifera) unveiled by Illumina- and single-molecule real-time-based RNA-sequencing

Alternative splicing (AS) plays a critical role in regulating different physiological and developmental processes in eukaryotes, by dramatically increasing the diversity of the transcriptome and the proteome. However, the saturation and complexity of AS remain unclear in lotus due to its limitation of rare obtainment of full-length multiple-splice isoforms. In this study, we apply a hybrid assembly strategy by combining single-molecule real-time sequencing and Illumina RNA-seq to get a comprehensive insight into the lotus transcriptomic landscape. We identified 211,802 high-quality full-length non-chimeric reads, with 192,690 non-redundant isoforms, and updated the lotus reference gene model. Moreover, our analysis identified a total of 104,288 AS events from 16,543 genes, with alternative 3ʹ splice-site being the predominant model, following by intron retention. By exploring tissue datasets, 370 tissue-specific AS events were identified among 12 tissues. Both the tissue-specific genes and isoforms might play important roles in tissue or organ development, and are suitable for ‘ABCE’ model partly in floral tissues. A large number of AS events and isoform variants identified in our study enhance the understanding of transcriptional diversity in lotus, and provide valuable resource for further functional genomic studies.

Maize NAC-domain retained splice variants act as dominant negatives to interfere with the full-length NAC counterparts

The plant-specific NAC transcription factors play diverse roles in various stress signaling. Alternative splicing is particularly prevalent in plants under stress. However, the investigation of cadmium (Cd) on the differential expression of the splice variants of NACs is in its infancy. Here, we identified three Cd-induced intron retention splice NAC variants which only contained the canonical NAC domain, designated as nacDomains, derived from three Cd-upregulated maize NACs. Subcellular localization analysis indicated that both nacDomain and its full-length NAC counterpart co-localized in the nucleus as manifested in the BiFC assay, thus implied that nacDomains and their corresponding NACs form heterodimers through the identical NAC domain. Further chimeric reporter/effector transient expression assay and Cd-tolerance assay in tobacco leaves collectively indicated that nacDomain-NAC heterodimers were involved in the regulation of NAC function. The results obtained here were in accordance with the model of dominant negative, which suggested that nacDomain act as the dominant negative to antagonize the regulation of NAC on its target gene expression and the Cd-tolerance function performance of NAC transcription factor. These findings proposed a novel insight into understanding the molecular mechanisms of Cd response in plants.

A perspective on crassulacean acid metabolism photosynthesis evolution of orchids on different continents: Dendrobium as a case study

Members of the Orchidaceae, one of the largest families of flowering plants, evolved the crassulacean acid metabolism (CAM) photosynthesis strategy. It is thought that CAM triggers adaptive radiation into new niche spaces, yet very little is known about its origin and diversification on different continents. Here, we assess the prevalence of CAM in Dendrobium, which is one of the largest genera of flowering plants and found in a wide range of environments, from the high altitudes of the Himalayas to relatively arid habitats in Australia. Based on phylogenetic time trees, we estimated that CAM, as determined by δ 13C values less negative than -20.0‰, evolved independently at least eight times in Dendrobium. The oldest lineage appeared in the Asian clade during the middle Miocene, indicating the origin of CAM was associated with a pronounced climatic cooling that followed a period of aridity. Divergence of the four CAM lineages in the Asian clade appeared to be earlier than divergence of those in the Australasian clade. However, CAM species in the Asian clade are much less diverse (25.6%) than those in the Australasian clade (57.9%). These findings shed new light on CAM evolutionary history and the aridity levels of the paleoclimate on different continents.

Single-molecule real-time sequencing facilitates the analysis of transcripts and splice isoforms of anthers in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

BACKGROUND

Anther development has been extensively studied at the transcriptional level, but a systematic analysis of full-length transcripts on a genome-wide scale has not yet been published. Here, the Pacific Biosciences (PacBio) Sequel platform and next-generation sequencing (NGS) technology were combined to generate full-length sequences and completed structures of transcripts in anthers of Chinese cabbage.

RESULTS

Using single-molecule real-time sequencing (SMRT), a total of 1,098,119 circular consensus sequences (CCSs) were generated with a mean length of 2664 bp. More than 75% of the CCSs were considered full-length non-chimeric (FLNC) reads. After error correction, 725,731 high-quality FLNC reads were estimated to carry 51,501 isoforms from 19,503 loci, consisting of 38,992 novel isoforms from known genes and 3691 novel isoforms from novel genes. Of the novel isoforms, we identified 407 long non-coding RNAs (lncRNAs) and 37,549 open reading frames (ORFs). Furthermore, a total of 453,270 alternative splicing (AS) events were identified and the majority of AS models in anther were determined to be approximate exon skipping (XSKIP) events. Of the key genes regulated during anther development, AS events were mainly identified in the genes SERK1, CALS5, NEF1, and CESA1/3. Additionally, we identified 104 fusion transcripts and 5806 genes that had alternative polyadenylation (APA).

CONCLUSIONS

Our work demonstrated the transcriptome diversity and complexity of anther development in Chinese cabbage. The findings provide a basis for further genome annotation and transcriptome research in Chinese cabbage.

Low temperature upregulates cwp expression and modifies alternative splicing patterns, increasing the severity of cwp-induced tomato fruit cuticular microfissures

The cwp (cuticular water permeability) gene controls the development of cuticular microfissuring and subsequent fruit dehydration in tomato. The gene underwent silencing in the evolution of the fleshy cultivated tomato but is expressed in the primitive wild tomato relatives. The introgression of the expressed allele from the wild S. habrochaites (cwp ^h ) into the cultivated tomato (Solanum lycopersicum) leads to the phenotype of fruit water loss during and following ripening. In this report, we show that low temperature impacts on the severity of the cuticular microfissure phenotype via a combination of effects on both expression and alternative splicing of cwp ^h . The cwp gene, comprising four exons and three introns, undergoes post-transcriptional alternative splicing processes, leading to seven alternative transcripts that differ in reading-frame lengths. Transgenic plants expressing each of the alternative transcripts identified the longest reading frame (VAR1) as the functional splice variant. Low temperature led to a strong upregulation of cwp ^h expression, compounded by an increase in the relative proportion of the functional VAR1 transcript, leading to increased severity of microfissuring of the cuticle. In summary, we demonstrate the molecular mechanism behind the horticultural phenomenon of the low-temperature effect on cuticular microfissures in the dehydrating tomato.

Transcriptome-wide alternative splicing modulation during plant-pathogen interactions in wheat

Alternative splicing (AS) enhances the diversities of both transcripts and proteins in eukaryotes, which contribute to stress adaptation. To catalog wheat (Triticum aestivum L.) AS genes, we characterized 45 RNA-seq libraries from wheat seedlings infected by powdery mildew, Blumeria graminis f. sp. tritici (Bgt) or stripe rust fungus, Puccinia striiformis f. sp. tritici (Pst). We discovered that 11.2% and 10.4% of the multiexon genes had AS transcripts during Bgt and Pst infections, respectively. In response to fungal infection, wheat modulated AS not only in disease resistance proteins, but also in splicing related factors. Apart from the stress induced or activated splicing variants by pathogen, the differential expression profiles were fold increased through changing the ratio of full spliced transcripts versus intron retention (IR) transcripts. Comparing AS transcripts produced by the same gene in Bgt with Pst stress, the spliced terminal exons and the stranded introns are independent and different. This demonstrated that differential induction of specific splice variants were activated against two fungal pathogens. The specific induced AS genes in the Pst-resistant plants were enriched in improving the membrane permeability and protein modification ability, whereas gene expression involved in protein translation and transport were strengthened in Pst-susceptible plants.

Distinct types of short open reading frames are translated in plant cells

Genomes contain millions of short (<100 codons) open reading frames (sORFs), which are usually dismissed during gene annotation. Nevertheless, peptides encoded by such sORFs can play important biological roles, and their impact on cellular processes has long been underestimated. Here, we analyzed approximately 70,000 transcribed sORFs in the model plant Physcomitrella patens (moss). Several distinct classes of sORFs that differ in terms of their position on transcripts and the level of evolutionary conservation are present in the moss genome. Over 5000 sORFs were conserved in at least one of 10 plant species examined. Mass spectrometry analysis of proteomic and peptidomic data sets suggested that tens of sORFs located on distinct parts of mRNAs and long noncoding RNAs (lncRNAs) are translated, including conserved sORFs. Translational analysis of the sORFs and main ORFs at a single locus suggested the existence of genes that code for multiple proteins and peptides with tissue-specific expression. Functional analysis of four lncRNA-encoded peptides showed that sORFs-encoded peptides are involved in regulation of growth and differentiation in moss. Knocking out lncRNA-encoded peptides resulted in a decrease of moss growth. In contrast, the overexpression of these peptides resulted in a diverse range of phenotypic effects. Our results thus open new avenues for discovering novel, biologically active peptides in the plant kingdom.

Genome-wide analysis of alternative splicing divergences between Brassica hexaploid and its parents

Compared with its parents, Brassica hexaploid underwent significant AS changes, which may provide diversified gene expression regulation patterns and could enhance its adaptability during evolution Polyploidization is considered a significant evolution force that promotes species formation. Alternative splicing (AS) plays a crucial role in multiple biological processes during plant growth and development. To explore the effects of allopolyploidization on the AS patterns of genes, a genome-wide AS analysis was performed by RNA-seq in Brassica hexaploid and its parents. In total, we found 7913 (27540 AS events), 14447 (70179 AS events), and 13205 (60804 AS events) AS genes in Brassica rapa, Brassica carinata, and Brassica hexaploid, respectively. A total of 920 new AS genes were discovered in Brassica hexaploid. There were 56 differently spliced genes between Brassica hexaploid and its parents. In addition, most of the alternative 5' splice sites were located 4 bp upstream of the dominant 5' splice sites, and most of the alternative 3' splice sites were located 3 bp downstream of the dominant 3' splice sites in Brassica hexapliod, which was similar to B. carinata. Furthermore, we cloned and sequenced all amplicons from the RT-PCR products of GRP7/8, namely, Bol045859, Bol016025 and Bol02880. The three genes were found to produce AS transcripts in a new way. The AS patterns of genes were diverse between Brassica hexaploid and its parents, including the loss and gain of AS events. Allopolyploidization changed alternative splicing sites of pre-mRNAs in Brassica hexaploid, which brought about alterations in the sequences of transcripts. Our study provided novel insights into the AS patterns of genes in allopolyploid plants, which may provide a reference for the study of polyploidy adaptability.

Iso-Seq analysis of the Taxus cuspidata transcriptome reveals the complexity of Taxol biosynthesis

BACKGROUND

Taxus cuspidata is well known worldwide for its ability to produce Taxol, one of the top-selling natural anticancer drugs. However, current Taxol production cannot match the increasing needs of the market, and novel strategies should be considered to increase the supply of Taxol. Since the biosynthetic mechanism of Taxol remains largely unknown, elucidating this pathway in detail will be very helpful in exploring alternative methods for Taxol production.

RESULTS

Here, we sequenced Taxus cuspidata transcriptomes with next-generation sequencing (NGS) and third-generation sequencing (TGS) platforms. After correction with Illumina reads and removal of redundant reads, more than 180,000 nonredundant transcripts were generated from the raw Iso-Seq data. Using Cogent software and an alignment-based method, we identified a total of 139 cytochrome P450s (CYP450s), 31 BAHD acyltransferases (ACTs) and 1940 transcription factors (TFs). Based on phylogenetic and coexpression analysis, we identified 9 CYP450s and 7 BAHD ACTs as potential lead candidates for Taxol biosynthesis and 6 TFs that are possibly involved in the regulation of this process. Using coexpression analysis of genes known to be involved in Taxol biosynthesis, we elucidated the stem biosynthetic pathway. In addition, we analyzed the expression patterns of 12 characterized genes in the Taxol pathway and speculated that the isoprene precursors for Taxol biosynthesis were mainly synthesized via the MEP pathway. In addition, we found and confirmed that the alternative splicing patterns of some genes varied in different tissues, which may be an important tissue-specific method of posttranscriptional regulation.

CONCLUSIONS

A strategy was developed to generate corrected full-length or nearly full-length transcripts without assembly to ensure sequence accuracy, thus greatly improving the reliability of coexpression and phylogenetic analysis and greatly facilitating gene cloning and characterization. This strategy was successfully utilized to elucidate the Taxol biosynthetic pathway, which will greatly contribute to the goals of improving the Taxol content in Taxus spp. using molecular breeding or plant management strategies and synthesizing Taxol in microorganisms using synthetic biological technology.

Arabidopsis SME1 Regulates Plant Development and Response to Abiotic Stress by Determining Spliceosome Activity Specificity

The control of precursor-messenger RNA (pre-mRNA) splicing is emerging as an important layer of regulation in plant responses to endogenous and external cues. In eukaryotes, pre-mRNA splicing is governed by the activity of a large ribonucleoprotein machinery, the spliceosome, whose protein core is composed of the Sm ring and the related Sm-like 2-8 complex. Recently, the Arabidopsis (Arabidopsis thaliana) Sm-like 2-8 complex has been characterized. However, the role of plant Sm proteins in pre-mRNA splicing remains largely unknown. Here, we present the functional characterization of Sm protein E1 (SME1), an Arabidopsis homolog of the SME subunit of the eukaryotic Sm ring. Our results demonstrate that SME1 regulates the spliceosome activity and that this regulation is controlled by the environmental conditions. Indeed, depending on the conditions, SME1 ensures the efficiency of constitutive and alternative splicing of selected pre-mRNAs. Moreover, missplicing of most targeted pre-mRNAs leads to the generation of nonsense-mediated decay signatures, indicating that SME1 also guarantees adequate levels of the corresponding functional transcripts. In addition, we show that the selective function of SME1 in ensuring appropriate gene expression patterns through the regulation of specific pre-mRNA splicing is essential for adequate plant development and adaptation to freezing temperatures. These findings reveal that SME1 plays a critical role in plant development and interaction with the environment by providing spliceosome activity specificity.

Emerging Roles of LSM Complexes in Posttranscriptional Regulation of Plant Response to Abiotic Stress

It has long been assumed that the wide reprogramming of gene expression that modulates plant response to unfavorable environmental conditions is mainly controlled at the transcriptional level. A growing body of evidence, however, indicates that posttranscriptional regulatory mechanisms also play a relevant role in this control. Thus, the LSMs, a family of proteins involved in mRNA metabolism highly conserved in eukaryotes, have emerged as prominent regulators of plant tolerance to abiotic stress. Arabidopsis contains two main LSM ring-shaped heteroheptameric complexes, LSM1-7 and LSM2-8, with different subcellular localization and function. The LSM1-7 ring is part of the cytoplasmic decapping complex that regulates mRNA stability. On the other hand, the LSM2-8 complex accumulates in the nucleus to ensure appropriate levels of U6 snRNA and, therefore, correct pre-mRNA splicing. Recent studies reported unexpected results that led to a fundamental change in the assumed consideration that LSM complexes are mere components of the mRNA decapping and splicing cellular machineries. Indeed, these data have demonstrated that LSM1-7 and LSM2-8 rings operate in Arabidopsis by selecting specific RNA targets, depending on the environmental conditions. This specificity allows them to actively imposing particular gene expression patterns that fine-tune plant responses to abiotic stresses. In this review, we will summarize current and past knowledge on the role of LSM rings in modulating plant physiology, with special focus on their function in abiotic stress responses.

Protein interference for regulation of gene expression in plants

Transcription factors (TFs) play a central role in the gene regulation associated with a plant's development and its response to the environmental factors. The work of TFs is well regulated at each stage of their activities. TFs usually consist of three protein domains required for DNA binding, dimerization, and transcriptional regulation. Alternative splicing (AS) produces multiple proteins with varying composition of domains. Recent studies have shown that AS of some TF genes form small proteins (small interfering peptide/small interfering protein, siPEP/siPRoT), which lack one or more domains and negatively regulate target TFs by the mechanism of protein interference (peptide interference/protein interference, PEPi/PROTi). The presence of an alternative form for the transcription factor CCA1 of Arabidopsis thaliana, has been shown to be involved in the regulation of the response to cold stress. For the PtFLC protein, one of the isoforms was found, which is formed as a result of alternative splicing and acts as a negative repressor, binding to the full-length TF PtFLC and therefore regulating the development of the Poncirus trifoliata. For A. thaliana, a FLM gene was found forming the FLM-б isoform, which acts as a dominant negative regulator and stimulates the development of the flower formation process due to the formation of a heterodimer with SVP TF. Small interfering peptides and proteins can actively participate in the regulation of gene expression, for example, in situations of stress or at different stages of plant development. Moreover, small interfering peptides and proteins can be used as a tool for fundamental research on the function of genes as well as for applied research for permanent or temporary knockout of genes. In this review, we have demonstrated recent studies related to siPEP/siPROT and their involvement in the response to various stresses, as well as possible ways to obtain small proteins.

Molecular improvement of alfalfa for enhanced productivity and adaptability in a changing environment

Due to an expanding world population and increased buying power, the demand for ruminant products such as meat and milk is expected to grow substantially in coming years, and high levels of forage crop production will therefore be a necessity. Unfortunately, urbanization of agricultural land, intensive agricultural practices, and climate change are all predicted to limit crop production in the future, which means that the development of forage cultivars with improved productivity and adaptability will be essential. Because alfalfa (Medicago sativa L.) is one of the most widely cultivated perennial forage crops, it has been the target of much research in this field. In this review, we discuss progress that has been made towards the improvement of productivity, abiotic stress tolerance, and nutrient-use efficiency, as well as disease and pest resistance, in alfalfa using biotechnological techniques. Furthermore, we consider possible future priorities and avenues for attaining further enhancements in this crop as a means of contributing to the realization of food security in a changing environment.

Genome-wide identification and characterization of mRNAs and lncRNAs involved in cold stress in the wild banana (Musa itinerans)

Cold stress seriously affects banana growth, yield and fruit quality. Long noncoding RNAs (lncRNAs) have been demonstrated as key regulators of biotic and abiotic stress in plants, but the identification and prediction of cold responsive mRNAs and lncRNAs in wild banana remains unexplored. In present study, a cold resistant wild banana line from China was used to profile the cold-responsive mRNAs and lncRNAs by RNA-seq under cold stress conditions, i.e. 13°C (critical growth temperature), 4°C (chilling temperature), 0°C (freezing temperature) and normal growing condition, i.e. 28°C (control group). A total of 12,462 lncRNAs were identified in cold-stressed wild banana. In mRNA, much more alternative splicing events occurred in wild banana under the cold stress conditions compared with that in the normal growing condition. The GO analysis of differential expression genes (DEGs) showed the biochemical processes and membrane related genes responded positively to the cold stress. The KEGG pathway enrichment analysis of the DEGs showed that the pathways of photosynthesis, photosynthesis–antenna proteins, circadian rhythm–plant, glutathione metabolism, starch and sucrose metabolism, cutin/suberine/biosynthesis were altered or affected by the cold stress conditions. Our analyses of the generated transcriptome and lncRNAs provide new insights into regulating expression of genes and lncRNAs that respond to cold stress in the wild banana.

Genome-Wide Association Analyses Reveal the Importance of Alternative Splicing in Diversifying Gene Function and Regulating Phenotypic Variation in Maize

Alternative splicing (AS) enhances transcriptome diversity and plays important roles in regulating plant processes. Although widespread natural variation in AS has been observed in plants, how AS is regulated and contribute to phenotypic variation is poorly understood. Here, we report a population-level transcriptome assembly and genome-wide association study to identify splicing quantitative trait loci (sQTLs) in developing maize (Zea mays) kernels from 368 inbred lines. We detected 19,554 unique sQTLs for 6570 genes. Most sQTLs showed small isoform usage changes without involving major isoform switching between genotypes. The sQTL-affected isoforms tend to display distinct protein functions. We demonstrate that nonsense-mediated mRNA decay, microRNA-mediated regulation, and small interfering peptide-mediated peptide interference are frequently involved in sQTL regulation. The natural variation in AS and overall mRNA level appears to be independently regulated with different cis-sequences preferentially used. We identified 214 putative trans-acting splicing regulators, among which ZmGRP1, encoding an hnRNP-like glycine-rich RNA binding protein, regulates the largest trans-cluster. Knockout of ZmGRP1 by CRISPR/Cas9 altered splicing of numerous downstream genes. We found that 739 sQTLs colocalized with previous marker-trait associations, most of which occurred without changes in overall mRNA level. Our findings uncover the importance of AS in diversifying gene function and regulating phenotypic variation.

Splicing regulation by long noncoding RNAs

Massive high-throughput sequencing techniques allowed the identification of thousands of noncoding RNAs (ncRNAs) and a plethora of different mRNA processing events occurring in higher organisms. Long ncRNAs can act directly as long transcripts or can be processed into active small si/miRNAs. They can modulate mRNA cleavage, translational repression or the epigenetic landscape of their target genes. Recently, certain long ncRNAs have been shown to play a crucial role in the regulation of alternative splicing in response to several stimuli or during disease. In this review, we focus on recent discoveries linking gene regulation by alternative splicing and its modulation by long and small ncRNAs.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Reciprocal cross-regulation of VND and SND multigene TF families for wood formation in Populus trichocarpa

Secondary cell wall (SCW) biosynthesis is the biological process that generates wood, an important renewable feedstock for materials and energy. NAC domain transcription factors, particularly Vascular-Related NAC-Domain (VND) and Secondary Wall-Associated NAC Domain (SND) proteins, are known to regulate SCW differentiation. The regulation of VND and SND is important to maintain homeostasis for plants to avoid abnormal growth and development. We previously identified a splice variant, PtrSND1-A2^IR , derived from PtrSND1-A2 as a dominant-negative regulator, which suppresses the transactivation of all PtrSND1 family members. PtrSND1-A2^IR also suppresses the self-activation of the PtrSND1 family members except for its cognate transcription factor, PtrSND1-A2, suggesting the existence of an unknown factor needed to regulate PtrSND1-A2 Here, a splice variant, PtrVND6-C1^IR , derived from PtrVND6-C1 was discovered that suppresses the protein functions of all PtrVND6 family members. PtrVND6-C1^IR also suppresses the expression of all PtrSND1 members, including PtrSND1-A2, demonstrating that PtrVND6-C1^IR is the previously unidentified regulator of PtrSND1-A2 We also found that PtrVND6-C1^IR cannot suppress the expression of its cognate transcription factor, PtrVND6-C1PtrVND6-C1 is suppressed by PtrSND1-A2^IR Both PtrVND6-C1^IR and PtrSND1-A2^IR cannot suppress their cognate transcription factors but can suppress all members of the other family. The results indicate that the splice variants from the PtrVND6 and PtrSND1 family may exert reciprocal cross-regulation for complete transcriptional regulation of these two families in wood formation. This reciprocal cross-regulation between families suggests a general mechanism among NAC domain proteins and likely other transcription factors, where intron-retained splice variants provide an additional level of regulation.

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.

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.

Response of Gene Expression and Alternative Splicing to Distinct Growth Environments in Tomato

Phenotypic plasticity is the phenomenon that one particular genotype produces different phenotypes under different environmental conditions, but its underlying molecular and genetic mechanisms are poorly understood. Plastic traits may be under the control of genes whose expression is modulated by environmental cues. In this study, we investigated phenotypic plasticity in tomato (Solanum lycopersicum) and its ancestral species S. pimpinellifolium by comparing the global gene expression of young seedlings grown under two distinct growth conditions. Our results show that more than 7000 genes exhibited differential expression in response to environmental changes from phytotron to a plastic greenhouse, and 98 environmentally sensitive genes displayed the same patterns of expression response across the two tomato species. We also found that growth conditions had a remarkable impact on transcriptome complexity, attributable to alternative splicing (AS), in which 665 splice variants showed differential expression in response to the environmental changes. Moreover, more splice variants and AS events per gene were detected in plastic greenhouse-grown seedlings than their phytotron counterparts, and these seedlings also had higher percentages of intron retention events. The identification of the conserved environmentally-sensitive genes and the splice variants in this study will be useful for further analysis of gene regulation of environmental response in tomato and other crops.

A sugarcane R2R3-MYB transcription factor gene is alternatively spliced during drought stress

MYB transcription factors of the R2R3-MYB family have been shown to play important roles in many plant processes. A sugarcane R2R3-MYB gene (ScMYB2) and its two alternative forms of transcript (ScMYB2S1 and ScMYB2S2) were identified in this study. The deduced protein of ScMYB2S1 is a typical plant R2R3-MYB protein, while ScMYB2S2 encodes a truncated protein. Real-time qPCR analysis revealed that ScMYB2S1 is suppressed under PEG-simulated drought stress in sugarcane, while ScMYB2S2 is induced at later treatment stage. A senescence symptom was observed when ScMYB2S1 was injected into tobacco leaves mediated by Agrobacterium, but no symptom for ScMYB2S2. Further investigation showed that the expression levels of 4 senescence-associated genes, NtPR-1a, NtNYC1, NtCAT3 and NtABRE, were markedly induced in tobacco leaves after ScMYB2S1-injection, while they were not sensitive to ScMYB2S2-injection. Moreover, MDA and proline were also investigated after injection. Similarly, MDA and proline levels were induced by ABA and ScMYB2S1, while inhibited by ScMYB2S2. We propose that ScMYB2, by alternatively splicing two transcripts (ScMYB2S1 and ScMYB2S2), is involved in an ABA-mediated leaf senescence signaling pathway and play positive role in respond to drought-induced senescence in sugarcane. The results of this study provide information for further research in sugarcane stress processes.