Alternative splicing during fruit development among fleshy fruits

Evidence Supporting a Role of Alternative Splicing Participates in Melon (Cucumis melo L.) Fruit Ripening

One key post-transcriptional modification mechanism that dynamically controls a number of physiological processes in plants is alternative splicing (AS). However, the functional impacts of AS on fruit ripening remain unclear. In this research, we used RNA-seq data from climacteric (VED, Harukei 3) and non-climacteric (PI, PS) melon cultivars to explore alternative splicing (AS) in immature and mature fruit. The results revealed dramatic changes in differential AS genes (DAG) between the young and mature fruit stages, particularly in genes involved in fruit development/ripening, carotenoid and capsaicinoid biosynthesis, and starch and sucrose metabolism. Serine/arginine-rich (SR) family proteins are known as important splicing factors in AS events. From the melon genome, a total of 17 SR members were discovered in this study. These genes could be classified into eight distinct subfamilies based on gene structure and conserved motifs. Promoter analysis detected various cis-acting regulatory elements involved in hormone pathways and fruit development. Interestingly, these SR genes exhibited specific expression patterns in reproductive organs such as flowers and ovaries. Additionally, concurrent with the increase in AS levels in ripening fruit, the transcripts of these SR genes were activated during fruit maturation in both climacteric and non-climacteric melon varieties. We also found that most SR genes were under selection during domestication. These results represent a novel finding of increased AS levels and SR gene expression during fruit ripening, indicating that alternative splicing may play a role in fruit maturation.

Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

Transcriptomic Analysis of Alternative Splicing Events during Different Fruit Ripening Stages of Coffea arabica L

To date, genomic and transcriptomic data on Coffea arabica L. in public databases are very limited, and there has been no comprehensive integrated investigation conducted on alternative splicing (AS). Previously, we have constructed and sequenced eighteen RNA-seq libraries of C. arabica at different ripening stages of fruit development. From this dataset, a total of 3824, 2445, 2564, 2990, and 3162 DSGs were identified in a comparison of different fruit ripening stages. The largest proportion of DSGs, approximately 65%, were of the skipped exon (SE) type. Biologically, 9 and 29 differentially expressed DSGs in the spliceosome pathway and carbon metabolism pathway, respectively, were identified. These DSGs exhibited significant variations, primarily in S1 vs. S2 and S5 vs. S6, and they involve many aspects of organ development, hormone transduction, and the synthesis of flavor components. Through the examination of research findings regarding the biological functions and biochemical pathways associated with DSGs and DEGs, it was observed that six DSGs significantly enriched in ABC transporters, namely, LOC113712394, LOC113726618, LOC113739972, LOC113725240, LOC113730214, and LOC113707447, were continually down-regulated at the fruit ripening stage. In contrast, a total of four genes, which were LOC113732777, LOC113727880, LOC113690566, and LOC113711936, including those enriched in the cysteine and methionine metabolism, were continually up-regulated. Collectively, our findings may contribute to the exploration of alternative splicing mechanisms for focused investigations of potential genes associated with the ripening of fruits in C. arabica.

Comprehensive Profiling of Alternative Splicing and Alternative Polyadenylation during Fruit Ripening in Watermelon (Citrullus lanatus)

Fruit ripening is a highly complicated process that is accompanied by the formation of fruit quality. In recent years, a series of studies have demonstrated post-transcriptional control play important roles in fruit ripening and fruit quality formation. Till now, the post-transcriptional mechanisms for watermelon fruit ripening have not been comprehensively studied. In this study, we conducted PacBio single-molecule long-read sequencing to identify genome-wide alternative splicing (AS), alternative polyadenylation (APA) and long non-coding RNAs (lncRNAs) in watermelon fruit. In total, 6,921,295 error-corrected and mapped full-length non-chimeric (FLNC) reads were obtained. Notably, more than 42,285 distinct splicing isoforms were derived from 5,891,183 intron-containing full-length FLNC reads, including a large number of AS events associated with fruit ripening. In addition, we characterized 21,506 polyadenylation sites from 11,611 genes, 8703 of which have APA sites. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that fructose and mannose metabolism, starch and sucrose metabolism and carotenoid biosynthesis were both enriched in genes undergoing AS and APA. These results suggest that post-transcriptional regulation might potentially have a key role in regulation of fruit ripening in watermelon. Taken together, our comprehensive PacBio long-read sequencing results offer a valuable resource for watermelon research, and provide new insights into the molecular mechanisms underlying the complex regulatory networks of watermelon fruit ripening.

Development of Nanopore sequencing-based full-length transcriptome database toward functional genome annotation of the Pacific oyster, Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is a widely cultivated shellfish in the world, while its transcriptome diversity remains less unexplored due to the limitation of short reads. In this study, we used Oxford Nanopore sequencing to develop the full-length transcriptome database of C. gigas. We identified 77,920 full-length transcripts from 21,523 genes, and uncovered 9668 alternative splicing events and 87,468 alternative polyadenylation sites. Notably, a total of 16,721 novel transcripts were annotated in this work. Furthermore, integrative analysis of 25 publicly available RNA-seq datasets revealed the transcriptome diversity involved in post-transcriptional regulation in C. gigas. We further developed a Drupal based webserver, Cgtdb, which can be used for transcriptome visualization, sequence alignment, and functional genome annotation analyses. This work provides valuable resources and a useful tool for integrative analysis of various transcriptome datasets in C. gigas, which will serve as an essential reference for functional annotation of the oyster genome.

Alternative Splicing during Fiber Development in G. hirsutum

Cotton is a valuable cash crop in many countries. Cotton fiber is a trichome that develops from a single epidermal cell and serves as an excellent model for understanding cell differentiation and other life processes. Alternative splicing (AS) of genes is a common post-transcriptional regulatory process in plants that is essential for plant growth and development. The process of AS during cotton fiber formation, on the other hand, is mainly unknown. A substantial number of multi-exon genes were discovered to be alternatively spliced during cotton fiber formation in this study, accounting for 23.31% of the total number of genes in Gossypium hirsutum. Retention intron (RI) is not necessarily the most common AS type, indicating that AS genes and processes during fiber development are very temporal and tissue-specific. When compared to fiber samples, AS is more prevalent at the fiber initiation stages and in the ovule, indicating that development stages and tissues use different AS strategies. Genes involved in fiber development have gone through stage-specific AS, demonstrating that AS regulates cotton fiber development. Furthermore, AS can be regulated by trans-regulation elements such as splicing factor and cis-regulation elements such as gene length, exon numbers, and GC content, particularly at exon-intron junction sites. Our findings also suggest that increased DNA methylation may aid in the efficiency of AS, and that gene body methylation is key in AS control. Finally, our research will provide useful information about the roles of AS during the cotton fiber development process.

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.

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

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