Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Urinary Long Non-Coding RNA Levels as Biomarkers of Lupus Nephritis
Previous Article in Journal
Effect of Substituents on Molecular Reactivity during Lignin Oxidation by Chlorine Dioxide: A Density Functional Theory Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 14
10.3390/ijms241411812
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Alternative Splicing during Fiber Development in G. hirsutum

by
Jing Zheng
†,
Shuhan Wen
†,
Zhipeng Yu
,
Keyan Luo
,
Junkang Rong
* and
Mingquan Ding
*
The Key Laboratory for Quality Improvement of Agricultural Products of Zhejiang Province, College of Advanced Agricultural Sciences, Zhejiang A&F University, Linan, Hangzhou 311300, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(14), 11812; https://doi.org/10.3390/ijms241411812
Submission received: 29 May 2023 / Revised: 12 July 2023 / Accepted: 20 July 2023 / Published: 22 July 2023
(This article belongs to the Section Molecular Informatics)
Browse Figures
Versions Notes

Abstract

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

1. Introduction

Cotton is an important fiber crop and plays an important role in the textile industry [1]. Cotton fiber is a linear single cell, originated and differentiated from epidermal cells of cotton ovule [2]. The fiber development process is complicated and overlapping, with five stages: cell initiation, rapid elongation, primary cell wall remodeling, secondary cell wall biosynthesis, and maturation [3]. A number of important genes responsible for cotton fiber development were identified with the whole-genome sequencing of the Gossypium hirsutum and RNA-seq technology [4,5]. These genes include transcription factors, genes related to the cell wall and phytohormones, and so on [6]. Although various dynamic transcriptome studies have revealed the molecular basis for fiber development, the post-transcriptional modification regulation mechanism of genes related to cotton fiber development is still unknown.
Post-transcriptional modifications include alternative splicing (AS), RNA editing, and mRNA decay [7]. AS is a pre-mRNA process mechanism that enables the production of many transcript isoforms from a single gene, and each transcript can generate a different protein [8]. More than 60%–85% of genes in plants are AS genes [9,10]. Many studies have indicated that AS is highly related to plant defense response, such as high temperature, cold, salt, and drought stress [11,12,13]. AS has also been recently proven to be related to plant growth and development, such as seed germination, circadian clock control, photomorphogenesis, flowing time control, and ovule development in Arabidopsis and rice [14,15,16,17,18,19]. Whether AS event regulation in polyploid crops, such as wheat, cotton, and oilseed rape, shares the same mechanism as the diploid remains unclear compared with the AS analysis in Arabidopsis and rice [20]. A study has shown that AS variation primarily occurs after polyploidization and domestication in wheat [21]. During embryogenesis and grain development in wheat, a global analysis of the regulation of AS in diploid grass and polyploid wheat grains revealed diversity in AS events not only between the endosperm, pericarp, and embryo overdevelopment but also between subgenomes [22]. Researchers have found that AS in cotton occurs in many gene families, such as peroxiredoxin [23] and the MADS-box family of genes [24]. However, a limited number of whole-genome-wide analyses of the AS genes and AS events exist in the cotton species, where in most of them are from diploid cotton [25,26,27] and one is from tetraploid cotton Gossypium barbadense [28]. To date, few systematic studies have been conducted on whether AS plays an important role in G. hirsutum, especially in cotton fiber development.
AS is mainly regulated by the activity of splicing regulators. Changes in the abundance of these proteins through transcription and AS, post-translational modifications, and interactions with exonic and intronic cis-elements and core elements of the spliceosomes modulate the outcome of pre-mRNA splicing [13]. The AS genes could also be regulated by miRNA, nucleosome occupancy, and DNA methylation [28,29,30]. Little research has been undertaken on the cloning and function of splicing factors in cotton because research on AS in cotton is still in the early stages, and the regulation mechanism of AS is still unclear.
Because little research has been conducted on AS in plant development, particularly in the polar development phase of cotton fiber cell initiation and elongation, the variations and potential roles of AS genes in cotton fiber development are unknown. We used publicly accessible, high-density time point transcriptome data during cotton fiber development to conduct a complete time-course AS analysis in G. hirsutum. According to our findings, AS variation occurs often during fiber production in G. hirsutum. Furthermore, we discovered large sources of other regulatory mechanisms related to cotton fiber AS in our investigation.

2. Results

2.1. Statistical Analysis of AS Genes and AS Events

The AS events during the ovule and fiber development in G. hirsutum were identified by using the transcriptome data with rMATs software (version:4.0.2) [31]. Given that the fiber development is a dynamic and sophisticated development process with overlapping stages, we roughly divided our samples into three stages: namely, the fiber initiation and primary elongation period (stage 1: −3 DPA–5 DPA), fiber elongation period (stage 2: 10 DPA–15 DPA), and secondary wall thickening period (stage 3: 20 DPA–25 DPA). By customizing a user-friendly program, five major AS events were identified, including skipped exon (SE), alternative 5′ splicing site (A5SS), alternative 3′ splicing site(A3SS), mutually exclusive exons (MXE), and retained intron(RI). Figure 1A shows the corresponding types of AS events in G. hirsutum. A total of 28,768 pairs of At and Dt homologous genes were extracted via a bidirectional blast and used as the research objects to compare the differences between AS events in G. hirsutum.
When the screening thresholds of the AS events were set to 0.1 < PSI (percentage spliced in index) < 0.9 and SJC (skipping junction counts) ≥ 516,283 AS genes (47,751 AS events) were detected during the cotton fiber development in G. hirsutum, accounting for 23.31% of the total number of genes in the G. hirsutum (Figure 1B and Table 1). Each AS gene produced 2.93 AS events on average. Approximately 41.01% of the total AS genes produced only one AS event. Moreover, the number of AS hotspot genes is 2872, accounting for 17.63% of the total AS genes. These genes could produce more than four AS events during fiber development. Gh_A11G2975 was the most variable splicing gene in G. hirsutum, which produced 39 different AS events. Furthermore, 41,253 genes did not have any AS events in all samples, accounting for 71.69% of the total number of G. hirsutum genes.
The AS genes and AS events were not evenly distributed across all of the G. hirsutum chromosomes. The chrA05 chromosome has the largest number of AS genes and AS events, with a total of 1114 AS genes (3372 AS events). Meanwhile, the chrA04 chromosome has the smallest number of AS genes and AS events, with a total of 308 AS genes (930 AS events). The AS genes and AS events are mainly distributed in the euchromatin region at both ends of each chromosome but less distributed near the centromere in the heterochromatin region (Figure 2).
AS events can be divided into five types according to the splicing sites: variable 3′-end site (A3SS), variable 5′-end site (A5SS), exon mutually exclusion (MXE), intron retention (RI), and exon jumping (SE). Among the AS genes, 206 AS genes (1.27%) produced five types of AS events (A3SS, A5SS, MXE, RI, and SE), while 8453 AS genes (51.91%) produced only one type of AS event. The number of AS events revealed that the A3SS type was the most abundant AS event type (30.67–38.66%) among all the time course samples of different periods and tissues, followed by SE (20.34–29.67%), RI (6.79–27.64%), A5SS (6.79–27.64%), and MXE (3.29–8.17%). The number of RI types of AS events fluctuated in samples of different periods and tissues, indicating that it has great temporal and tissue specificity (Figure 3).
G. hirsutum is an allotetraploid that contains At and Dt subgenomes. In this study, the differences between AS genes and AS events were compared between the two subgenomes. Totals of 8, 176 and 8107 AS genes were detected in the At and Dt subgenomes of G. hirsutum involving 24,301 and 23,450 AS events, accounting for 28.42% and 28.18% of the total At and Dt subgenomes, respectively. Each At- and Dt-subgenome-splicing gene produced 2.97 and 2.89 splicing events on average, respectively. The statistical analysis showed significant differences in the number of AS genes and AS events among G. hirsutum subgenomes. The At subgenome had more AS genes (p-value = 8.24 × 10−11) and AS events (p-value = 2.03 × 10−13) compared with the Dt subgenome (Figure 4).
According to the cotton fiber’s morphological changes, a large number of AS genes and AS events occurred during the three stages of the fiber development of G. hirsutum. In ovule samples, 14,154, 9441, and 5705 AS genes (Figure 5A and Supplemental Table S2) involving 37,956, 19,639, and 9174 AS events (Figure 5B and Supplemental Table S3) were identified for stages 1, 2, and 3, respectively. In fiber samples, the number of AS genes and AS events in the ovule samples of the fiber initiation stage (stage 1) was significantly higher than that of the ovule and fiber samples of the fiber elongation (stage 2) and secondary wall thickening stage (stage 3). The 0 DPA ovule samples had the highest number of AS genes (9360) and AS events (19,920), while the 25 DPA fiber samples have the lowest number of AS genes (2363) and AS events (3187).

2.2. Identification of the Fiber-Development-Related AS Genes and AS Events

In our study, the AS genes and AS events detected in the ovule and fiber samples at each time point were defined as conservative AS genes and AS events. A total of 2910 conservative AS genes were found, accounting for 13.71% of the total number of AS genes, and 3018 conservative AS events were involved, accounting for 6.32% of the total number of AS events in G. hirsutum (Figure 5). The gene ontology (GO) enrichment analysis indicated that these genes are mostly enriched in biological processes (GO terms, such as GO:0006397: mRNA processing, GO:0015031: protein transport, GO:0042752: regulation of circadian rhythm, and GO:0045292: mRNA cis-splicing) via spliceosome (Supplemental Table S4). The large number of universally enriched GO terms indicates that AS is a common phenomenon in fiber development due to the substantial number of genes enriched in various biological process categories.
We discovered that 2384 unique AS genes (3130 AS events) existed in the fiber samples compared with that in the ovule (Supplemental Table S5). The GO enrichment analysis results revealed that these genes can mainly be categorized into GO terms such as “GO:0043412 macromolecule modification, GO:0016192 vesicle-mediated transport, GO:0006468 protein phosphorylation, GO:0016579 protein deubiquitination, GO:0006885 regulation of pH, GO:0006886 intracellular protein transport, and GO:0006511 ubiquitin-dependent protein catabolic process”. The enrichment analysis revealed that the rapid elongation of fibers may be related to the transportation of proteins in the cell, signal transductions, and changes in intracellular pH. Specific AS genes in the above-mentioned biological processes may play important roles in fiber elongation.
We also identified the specific AS genes in different stages of fiber development. A total of 37,956 AS events were detected in the samples at fiber initiation and primary elongation (stage 1), of which 18,424 were specific, involving 9506 AS genes. Moreover, 25,095 AS events were found in the ovule and fiber samples during fiber elongation (stage 2), of which 5602 AS events were specific, involving 4011 AS genes. Approximately 15,664 AS events were found in the samples at the secondary wall thickening period (stage 3), among which 2, 356 AS events were unique to this stage, involving 1868 AS genes (Figure 5).
The specific AS genes in stage 1 were mostly enriched in GO terms, such as “GO:0006396 RNA processing, GO:0006259 DNA metabolic process, GO:0016071 mRNA metabolic process, GO:0044265 cellular macromolecule catabolic process, and GO:0006281 DNA repair”, indicating that genes responsible for the genetic material’s metabolism are prone to be alternatively spliced during stage 1, which may play an important role in fiber initiation (Supplemental Table S6). The specific AS genes in stage 2 were mainly enriched in GO terms, such as “GO:0006644 phospholipid metabolic process, GO:0006654 phosphatidic acid biosynthetic process, GO:0046474 glycerophospholipid biosynthetic process, and GO:0044255 cellular lipid metabolic process”, indicating that AS of the dynamic lipid signal is a key cellular signaling mechanism for fiber elongation (Supplemental Table S7). The GO terms enriched in specific AS genes of stage 3 contained “GO:0043412 macromolecule modification, GO:0006464 cellular protein modification process, GO:0008064 regulation of actin polymerization or depolymerization, GO:0032956 regulation of actin cytoskeleton organization, and GO:0006996 organelle organization”, indicating that protein modification and cytoskeleton organization genes would play important roles in the secondary wall thickening stage of cotton fiber development (Supplemental Table S8).

2.3. Regulation of AS in Cotton

The AS is regulated by cis-acting elements within pre-mRNAs and trans-acting factors. The cis-acting elements are the splicing site and the branchpoint sequence that contains partially conserved motifs that are recognized by trans-acting factors. To identify the influencing factors of the AS, we first researched the expression of the pre-mRNA splicing factors (transacting factors) in our samples. Then, the gene length, exon/intron numbers, and GC content of the AS genes and the splicing site were deeply mined.
We first queried 59 splicing factors in Arabidopsis, which are responsible for mRNA processing, then 108 homologous genes in G. hirsutum were identified through the blast. The majority of the gene’s transcription levels remained stable during −3 DPA to 20 DPA. However, many of the splicing factors’ expression was significantly down-regulated in 25 DPA ovules and 10 DPA–25 DPA fibers. This finding is consistent with the result that the number of specific splicing genes in cotton fiber is less than that in the ovule (Figure 6).
After analyzing the splicing factor, we investigated whether AS genes have a special gene sequence characteristic. To simplify the study, only 3351 “RI group”- and 4896 “SE group”-type AS genes were included for further research. Genes that do not have any AS events were used as the “control group”. These data were statistically analyzed and compared between RI, SE, and normal groups to explore whether AS events have a bias on AS genes’ lengths, exon numbers, GC content of the entire gene, and the splicing site. The results revealed that the average length of normal-group genes was 1868.18 bp, while SE- and RI-group genes were 4517.96 and 4407.29 bp, respectively (Figure 7A). Moreover, each group in the normal group contained an average of 3.30 exons, while every gene in the SE and RI groups contained 9.44 and 10.03 exons on average, respectively (Figure 7B). The average GC content in the normal gene group was 40.15%. The GC contents of the SE and RI groups were 36.27% and 36.74%, respectively (Figure 7C). The statistical analysis showed that the SE and RI splicing events tended to occur in genes with longer sequence lengths, more exon numbers, and lower gene GC content.
The exon–intron splicing site is important for the correctness of the gene splicing. To explore the relationship between the occurrence of SE, RI, and the sequences near the splicing site, the 25 bp upstream and downstream sequences with a total length of 51 bp were extracted from the splicing site as the center to identify the sequence characteristics of the splicing site. The GC content of the sequence with the length of 20 bp on the intron side of the SE type of genes was not significantly different between the normal group’s introns (I1 = I1′, I2 = I2′) in the two splicing sites involved in the SE type (A’ and B’ splicing points) (A: p-value = 0.8553 and B: p-value = 0.4432). However, the GC content of the sequence of the SE-type exons was significantly lower than that of the normal type (E1′ < E1, E2′ < E2) (A: p-value = 2.51 × 10−11 and B: p-value = 6.31 × 10−11). Of the two splicing sites involved in the RI type (C and D splicing sites), the GC content on the exon side was different from the normal group’s exon (C: p-value = 1.24 × 10−4 and D: p-value = 0.0137). The GC content of the 20 bp sequence on the intron side was significantly higher than the normal group’s intron (I1′ > I1 and I2′ > I2) (C: p-value = 1.02 × 10−8 and D: p-value = 1.44 × 10−9) (Figure 8).
Due to significant differences detected in the GC content of the gene body regions and splicing sites between AS genes and normal splicing genes, this study also investigated the DNA methylation levels between the gene regions and splicing sites of AS genes. By analyzing the DNA methylome data of cotton 0 DPA ovule samples, it was found that in the gene body region, AS genes (SE, RI) showed a significant increase in m5C levels compared to non-splicing genes. The level of CpG methylation showed the same trend as m5C, while there was no significant change in the methylation levels of CHG and CHH. At the region of gene splicing sites, the m5C level of the SE type was slightly lower than that of RI type and normal splicing genes. However, the level of CpG methylation exhibited the same trend as the methylation level of the gene body region of AS genes (SE, RI), with higher levels of CpG methylation. Interestingly, the levels of CHG and CHH methylation were significantly reduced compared to normal splicing genes.

2.4. Different Types of AS Genes Methylation in Cotton

Due to significant differences detected in the GC content of the gene body regions and splicing sites between AS genes and normal splicing genes, this study also investigated the DNA methylation levels between the gene regions and splicing sites of AS genes. By analyzing the DNA methylome data of cotton 0 DPA ovule samples, it was found that in the gene body region, AS genes (SE, RI) showed a significant increase in m5C levels compared to non-splicing genes. The level of CpG methylation showed the same trend as m5C, while there was no significant change in the methylation levels of CHG and CHH. (Figure 9A) At the region of gene splicing sites, the m5C level of SE type is slightly lower than that of RI type and normal splicing genes. However, the level of CpG methylation is the same trend as the methylation level of the gene body region of AS genes (SE, RI), with higher levels of CpG methylation. Interestingly, the levels of CHG and CHH methylation were significantly reduced compared to normal splicing genes (Figure 9B).

3. Discussion

Because tolerance to abiotic stress is often characterized as polygenic adaptation, many adaptive loci contain genes with subtle effects, and AS has been extensively researched in plant response to abiotic stress. This could be connected to genes that are alternatively spliced at numerous loci and have distinct activities that contribute to abiotic stress [32,33]. In cotton, AS has been investigated in Gossypium davidsonii, Gossypium raimondii, and Gossypium barbadense. Similar in other plants, previous studies have mainly focused on the response of cotton to abiotic and biotic stresses [26,27,28]. The research material of the plant tissue is mainly on the leaf, with fewer studies on the AS genes and events of the ovule and fiber tissues in G. hirsutum. In this study, we explored the complicated post-transcriptional regulation during cotton fiber development in G. hirsutum via RNA-seq. The considerable AS events suggested that post-transcriptional regulation plays a crucial role in the fiber development process. Our data also indicate that many AS events have taken place in G. hirsutum ovules and fibers, which may contribute to the increase of the transcriptome diversity or proteomic complexity, consistent with the previous study [34]. Although AS is largely reported to relate to stress responses in the plant, a number of studies have demonstrated that proper AS of pre-mRNAs and functional genes controlled by the spliceosome is critical for the development of various plant tissues and organs [35,36]. Recently, AS has proven to be involved in JA and ABA signal transductions, indicating that the gene’s AS could participate in the hormone regulation of plant development [9,37,38]. Our study has revealed that the number of AS genes and AS events are highly spatiotemporal and tissue-specific during cotton fiber development. The trend is consistent with the previously defined fiber development stages, such as fiber initiation, elongation, secondary wall thickening stage, and cotton fiber maturation and dehydration stage [3]. The number of AS genes significantly increased at 0 DPA, 5 DPA, and 15 DPA. Given that 0 DPA is the key period in the process of fiber initiation, 5 DPA is the stage of fiber transition from initiation to elongation, and 15 DPA is the timepoint of transition from elongation to secondary wall thickening. Our results suggest that AS may be an important mechanism for G. hirsutum to regulate the growth and development of the fiber cells.
In contrast with the results that the number of RI-type AS genes accounts for the vast majority of all the types of AS genes in the relevant studies in cotton, this study found that the A3SS type has the largest number of AS genes, followed by the SE, RI, A5SS, and MXE types. In previous studies, the leaf has been used as the research material. Accordingly, we hypothesize that the distribution of the AS is related to tissues in cotton. In different samples of cotton fiber development, the proportion of RI-type events varies greatly (by up to 4.07 times). In the ovule samples of 10 DPA and 20 DPA, the RI type of AS genes accounted for less than 10%. In the ovules of 15 DPA and fiber samples of 10 DPA and 20 DPA, the RI type of AS genes accounted for more than 20%. The content of RI-type AS events during the fiber development has significant time and tissue specificities. Our results are highly consistent with recent research on AS during fruit development among fresh fruit, and they also have the same finding that the RI type of AS is tissue- and genus-specific [39].
AS has proven to be a critical regulatory mechanism for plant development throughout a plant’s entire lifespan [20]. Many genetic and gene function analyses suggest that the AS has been involved in seed dormancy, seed germination regulation, circadian clock control, and phytohormone signaling via DELAY OF GERMINATION 1 (DOG1), PHYTOCHROME-INTERACTING FACTOR 6 (PIF6), CIRCADIAN-CLOCK-ASSOCIATED 1 (CCA1), and HYPERSENSITIVE TO ABA1 (HAB1) gene [14,15,40,41]. Although a large number of AS genes and events have been found in cotton by our and the other groups, the key role of AS genes in cotton fiber development is still unknown. Here, the specific AS genes identified in the ovule and fiber tissues are enriched in GO terms highly related to fiber development. Hence, the AS genes’ function and their roles during cotton fiber development must be intensively explored further.
AS is regulated at several levels; to date, multiple sequence-conserved splicing factors and regulators have been functionally characterized in Arabidopsis, which regulate many aspects of development and stress responses [42]. We were able to identify 108 splicing factor genes in G. hirsutum. Some of these genes’ lower expression may lead to a smaller number of AS genes being present in cotton fiber samples. Furthermore, the gene structure and GC content may affect the preference for AS. The RI and SE types of AS genes have relatively longer gene lengths, lower GC content, and more exons at the gene level, indicating that AS genes have their characteristics. Another interesting result is that the GC content of the AS genes’ splicing site is different from the constitutive splicing genes, indicating that the junction site’s GC content and DNA methylation level would be important factors.
We discovered that increasing DNA methylation, namely CpG methylation in gene body areas, could improve the efficiency of AS in cotton. Our findings are congruent with those of other crops and plants, including maize and linseed [43,44]. Furthermore, a novel discovery was that the splicing junction’s methylation level contributes to the AS, with high CpG methylation associated with low CHH and CHG methylation, indicating that cis-element regulation via DNA methylation is crucial for gene post-transcriptional regulation such as AS.
In summary, our research has identified an intriguing phenomenon in which many critical functional genes involved in cotton fiber development undergo AS at different phAS events of fiber development. However, whether distinct spliced isoform proteins have diverse activities must be confirmed by cloning each isoform gene and validating overexpression experiments in cotton. Splicing is a complex process involving multiple RNA and protein components. Our findings indicate that during cotton fiber development, the expression of splicing factors significantly regulates the degree of AS. However, it is unclear which splicing factors are the key AS regulatory genes and which types of AS are regulated. In the future, precisely knocking out a splicing factor member with CRISPR-Cas9 could be used to study the biological functions of splicing factors and how to govern cotton fiber growth. Furthermore, our findings suggest that there are links between DNA methylation and AS. The mechanism through which DNA methylation influences fiber development by modulating the AS gene is currently unknown. Using whole genome bisulfite sequencing analysis (WGBS) to study the profile of DNA methylation of AS genes at various stages of cotton fiber development can be a good way to systematically study how DNA methylation maintains correct splicing of cotton fiber genes and thus regulates natural fiber development.

4. Materials and Methods

4.1. Data Availability

All of the transcriptome data used for this study came from NCBI’s SRA public database (accession numbers in Supplemental Table S1) [5]. In this study, a total of 34 relevant fastq samples were retrieved. The ovule samples ranged from −3 DPA to 25 DPA (−3 DPA, 0 DPA, 1 DPA, 3 DPA, 5 DPA, 10 DPA, 15 DPA, 20 DPA, 25 DPA). Samples of fiber were from the 10 DPA, 15 DPA, 20 DPA, and 25 DPA time periods. The downloads of the reference genome and associated annotation data are available on the CottonGen website [45].

4.2. Processing of the Raw Data and Alignment

SRAToolkit [46] was used to download the raw fastq data and filter it using fastp software(version: 0.23.4). The Hisat2 program [47] was then used to align all of the paired-end clean reads to the G. hirsutum reference genome. To construct the GTF file to be imported for AS analysis, the resultant SAM files were sorted and assembled using Cufflinks software (version: 2.2.1) [48]. Using CuffMerge tool (in Cufflink), the final GTF data were compared and merged with the original GTF file of the TM1 reference genome.

4.3. Performing AS and Statistical Analysis

The AS genes and events were identified with the software, rMATS (version:4.0.2) [31] according to the author’s manual. The IJC (inclusion junction counts) and SJC (skipping junction counts) were used to calculate the PSI values of the AS genes and AS events. AS events with support PSI > 0.1, and PSI < 0.9 and SJC ≥ 5 were retained and counted. Gene ontology enrichment analysis was performed using the OmicShare tools, a free online platform for data analysis (https://www.omicshare.com/).

4.4. Gene Expression Analysis

Taking G. hirsutum acc. TM-1 (version 2.1) as the reference genome, Htseq-counts [49] were used to calculate the number of reads of a gene in different samples, and the counts of each gene were obtained and normalized to be TPM. The R package, DESeq2, was used for differential expression between samples with the threshold of 2-fold expression changes [|log2(fold-change)| ≥ 1] and false discovery rate (FDR) < 0.05. The heatmap of gene’s expression was generated via the tbtools [50].

4.5. AS Genes’ Sequence Analysis

In order to analyze the structural and sequence differences of normal splicing genes (genes never undergoing AS) and AS genes in TM-1, the gene length, exon, and intron number, and GC content were extracted and calculated from the GTF file with the in-house Perl scripts. In order to identify the GC content of the splicing junction site, in this study, we defined the splice site as the center, and the upstream and downstream of each 25 bp sequence and the total length of 51 bp were extracted for the in-depth analysis of splice site’s GC characteristics.

4.6. Methylation’s Sequence Analysis

The cotton ovules’ (0 DPA) whole-genome bisulfite sequencing (WGBS) data were obtained from the NCBI SRA database (Accession number: GSE61774). The fastp software(version: 0.23.2) was used to filter the raw fastq file [51]. Using Bismark software (version: 0.23.1), the filtered reads were mapped on the Gossypium hirsutum TM1 genomic sequences [52]. The methylation level at each site was determined using the calmeth and methGff tool in the Batmeth2 software [53]. A sliding window approach was used to calculate the methylation level of each gene’s gene body region and 2000 bp upstream and downstream regulatory regions (Window size: 0.02 times the gene length or regulation regions, step size: 0.01 times the gene length or regulation regions). The sequences of the 50 bp upstream and downstream splicing sites of SE and RI genes, as well as each exon splicing site of normal splicing genes, were retrieved. These sites’ methylation levels were determined using sliding windows (Window size: 10 bp, Step size: 5 bp).

5. Conclusions

The AS genes and events during various phAS events of the development of upland cotton fibers were thoroughly investigated in this study. Although AS occurs regularly throughout the growth of fibers, the genes and events involved in this process were highly spatiotemporally and tissue-specific. With minimal impact on chromosomal length, the distribution of AS genes and AS events on each particular chromosome was strongly correlated with the gene density. A significant number of At and Dt subgenome AS genes and AS events was identified. In addition, we found that AS may be controlled by the way of the splicing factors’ expression and changes in the length, exon count, and GC content of the gene, particularly at the splicing junction location. This research contributes to a better understanding of the gene expression and AS regulatory mechanisms of genes involved in fiber development, as well as the AS profiles of G. hirsutum. The discovery of a potential link between AS and the fiber development in G. hirsutum may open up a new direction of research into the identification of the molecular mechanisms underlying cotton fiber development in the future.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms241411812/s1.

Author Contributions

Conceptualization, J.Z. and M.D.; methodology, data curation, S.W.; writing—original draft preparation, S.W.; writing—review and editing, Z.Y. and K.L.; project administration, J.R. and M.D.; funding acquisition, J.R. and M.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key R&D Program of China (No. 2021YFD1000100) and National Natural Science Foundation of China (No. 32171923).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All sequence data used in the analysis are available from NCBI under PRJNA490626 and PRJNA262137, and all scripts used to analyze the data are available from Multivariate Analysis of Transcript Splicing (MATS) under (Multivariate Analysis of Transcript Splicing (MATS) (sourceforge.net)).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, G.; Huang, J.Q.; Chen, X.Y.; Zhu, Y.X. Recent Advances and Future Perspectives in Cotton Research. Annu. Rev. Plant Biol. 2021, 72, 437–462. [Google Scholar] [CrossRef]
  2. Wang, S.; Wang, J.W.; Yu, N.; Li, C.H.; Luo, B.; Gou, J.Y.; Wang, L.J.; Chen, X.Y. Control of plant trichome development by a cotton fiber MYB gene. Plant Cell 2004, 16, 2323–2334. [Google Scholar] [CrossRef] [Green Version]
  3. Haigler, C.; Betancur, L.; Stiff, M.; Tuttle, J. Cotton fiber: A powerful single-cell model for cell wall and cellulose research. Front. Plant Sci. 2012, 3, 104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Li, F.; Fan, G.; Lu, C.; Xiao, G.; Zou, C.; Kohel, R.J.; Ma, Z.; Shang, H.; Ma, X.; Wu, J.; et al. Genome sequence of cultivated Upland cotton (Gossypium hirsutum TM-1) provides insights into genome evolution. Nat. Biotechnol. 2015, 33, 524–530. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Zhang, T.; Hu, Y.; Jiang, W.; Fang, L.; Guan, X.; Chen, J.; Zhang, J.; Saski, C.A.; Scheffler, B.E.; Stelly, D.M.; et al. Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement. Nat. Biotechnol. 2015, 33, 531–537. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Zhang, J.; Mei, H.; Lu, H.; Chen, R.; Hu, Y.; Zhang, T. Transcriptome Time-Course Analysis in the Whole Period of Cotton Fiber Development. Front. Plant Sci. 2022, 13, 864529. [Google Scholar] [CrossRef]
  7. Lam, P.Y.; Wang, L.; Lo, C.; Zhu, F.Y. Alternative Splicing and Its Roles in Plant Metabolism. Int. J. Mol. Sci. 2022, 23, 7355. [Google Scholar] [CrossRef] [PubMed]
  8. Syed, N.H.; Kalyna, M.; Marquez, Y.; Barta, A.; Brown, J.W. Alternative splicing in plants--coming of age. Trends Plant Sci. 2012, 17, 616–623. [Google Scholar] [CrossRef] [Green Version]
  9. Yang, X.; Jia, Z.; Pu, Q.; Tian, Y.; Zhu, F.; Liu, Y. ABA Mediates Plant Development and Abiotic Stress via Alternative Splicing. Int. J. Mol. Sci. 2022, 23, 3796. [Google Scholar] [CrossRef]
  10. Liu, J.; Chen, S.; Liu, M.; Chen, Y.; Fan, W.; Lee, S.; Xiao, H.; Kudrna, D.; Li, Z.; Chen, X.; et al. Full-Length Transcriptome Sequencing Reveals Alternative Splicing and lncRNA Regulation during Nodule Development in Glycine max. Int. J. Mol. Sci. 2022, 23, 7371. [Google Scholar] [CrossRef]
  11. Liu, X.X.; Guo, Q.H.; Xu, W.B.; Liu, P.; Yan, K. Rapid Regulation of Alternative Splicing in Response to Environmental Stresses. Front. Plant Sci. 2022, 13, 832177. [Google Scholar] [CrossRef] [PubMed]
  12. Zhou, Y.; Li, X.H.; Guo, Q.H.; Liu, P.; Li, Y.; Wu, C.A.; Yang, G.D.; Huang, J.G.; Zhang, S.Z.; Zheng, C.C.; et al. Salt responsive alternative splicing of a RING finger E3 ligase modulates the salt stress tolerance by fine-tuning the balance of COP9 signalosome subunit 5A. PLoS Genet. 2021, 17, e1009898. [Google Scholar] [CrossRef]
  13. Rosenkranz, R.R.E.; Ullrich, S.; Löchli, K.; Simm, S.; Fragkostefanakis, S. Relevance and Regulation of Alternative Splicing in Plant Heat Stress Response: Current Understanding and Future Directions. Front. Plant Sci. 2022, 13, 911277. [Google Scholar] [CrossRef] [PubMed]
  14. Nakabayashi, K.; Bartsch, M.; Ding, J.; Soppe, W.J. Seed Dormancy in Arabidopsis Requires Self-Binding Ability of DOG1 Protein and the Presence of Multiple Isoforms Generated by Alternative Splicing. PLoS Genet. 2015, 11, e1005737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. James, A.B.; Syed, N.H.; Bordage, S.; Marshall, J.; Nimmo, G.A.; Jenkins, G.I.; Herzyk, P.; Brown, J.W.; Nimmo, H.G. Alternative splicing mediates responses of the Arabidopsis circadian clock to temperature changes. Plant Cell 2012, 24, 961–981. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Dong, J.; Chen, H.; Deng, X.W.; Irish, V.F.; Wei, N. Phytochrome B Induces Intron Retention and Translational Inhibition of PHYTOCHROME-INTERACTING FACTOR3. Plant Physiol. 2020, 182, 159–166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Lutz, U.; Posé, D.; Pfeifer, M.; Gundlach, H.; Hagmann, J.; Wang, C.; Weigel, D.; Mayer, K.F.; Schmid, M.; Schwechheimer, C. Modulation of Ambient Temperature-Dependent Flowering in Arabidopsis thaliana by Natural Variation of FLOWERING LOCUS M. PLoS Genet. 2015, 11, e1005588. [Google Scholar] [CrossRef]
  18. Zhan, X.; Qian, B.; Cao, F.; Wu, W.; Yang, L.; Guan, Q.; Gu, X.; Wang, P.; Okusolubo, T.A.; Dunn, S.L.; et al. An Arabidopsis PWI and RRM motif-containing protein is critical for pre-mRNA splicing and ABA responses. Nat. Commun. 2015, 6, 8139. [Google Scholar] [CrossRef] [Green Version]
  19. Cucinotta, M.; Cavalleri, A.; Guazzotti, A.; Astori, C.; Manrique, S.; Bombarely, A.; Oliveto, S.; Biffo, S.; Weijers, D.; Kater, M.M.; et al. Alternative Splicing Generates a MONOPTEROS Isoform Required for Ovule Development. Curr. Biol. 2021, 31, 892–899.e893. [Google Scholar] [CrossRef]
  20. Liu, L.; Tang, Z.; Liu, F.; Mao, F.; Yujuan, G.; Wang, Z.; Zhao, X. Normal, novel or none: Versatile regulation from alternative splicing. Plant Signal Behav. 2021, 16, 1917170. [Google Scholar] [CrossRef] [PubMed]
  21. Yu, K.; Feng, M.; Yang, G.; Sun, L.; Qin, Z.; Cao, J.; Wen, J.; Li, H.; Zhou, Y.; Chen, X.; et al. Changes in Alternative Splicing in Response to Domestication and Polyploidization in Wheat. Plant Physiol. 2020, 184, 1955–1968. [Google Scholar] [CrossRef]
  22. Gao, P.; Quilichini, T.D.; Zhai, C.; Qin, L.; Nilsen, K.T.; Li, Q.; Sharpe, A.G.; Kochian, L.V.; Zou, J.; Reddy, A.S.N.; et al. Alternative splicing dynamics and evolutionary divergence during embryogenesis in wheat species. Plant Biotechnol. J. 2021, 19, 1624–1643. [Google Scholar] [CrossRef]
  23. Feng, Y.; Wei, R.; Liu, A.; Fan, S.; Che, J.; Zhang, Z.; Tian, B.; Yuan, Y.; Shi, G.; Shang, H. Genome-wide identification, evolution, expression, and alternative splicing profiles of peroxiredoxin genes in cotton. PeerJ 2021, 9, e10685. [Google Scholar] [CrossRef]
  24. Lightfoot, D.J.; Malone, K.M.; Timmis, J.N.; Orford, S.J. Evidence for alternative splicing of MADS-box transcripts in developing cotton fibre cells. Mol. Genet. Genom. 2008, 279, 75–85. [Google Scholar] [CrossRef]
  25. Feng, S.; Xu, M.; Liu, F.; Cui, C.; Zhou, B. Reconstruction of the full-length transcriptome atlas using PacBio Iso-Seq provides insight into the alternative splicing in Gossypium australe. BMC Plant Biol. 2019, 19, 365. [Google Scholar] [CrossRef] [Green Version]
  26. Li, Q.; Xiao, G.; Zhu, Y.X. Single-nucleotide resolution mapping of the Gossypium raimondii transcriptome reveals a new mechanism for alternative splicing of introns. Mol. Plant 2014, 7, 829–840. [Google Scholar] [CrossRef] [Green Version]
  27. Zhu, G.; Li, W.; Zhang, F.; Guo, W. RNA-seq analysis reveals alternative splicing under salt stress in cotton, Gossypium davidsonii. BMC Genom. 2018, 19, 73. [Google Scholar] [CrossRef]
  28. Wang, M.; Wang, P.; Liang, F.; Ye, Z.; Li, J.; Shen, C.; Pei, L.; Wang, F.; Hu, J.; Tu, L.; et al. A global survey of alternative splicing in allopolyploid cotton: Landscape, complexity and regulation. New Phytol. 2018, 217, 163–178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Thatcher, S.R.; Zhou, W.; Leonard, A.; Wang, B.B.; Beatty, M.; Zastrow-Hayes, G.; Zhao, X.; Baumgarten, A.; Li, B. Genome-wide analysis of alternative splicing in Zea mays: Landscape and genetic regulation. Plant Cell 2014, 26, 3472–3487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Chodavarapu, R.K.; Feng, S.; Bernatavichute, Y.V.; Chen, P.Y.; Stroud, H.; Yu, Y.; Hetzel, J.A.; Kuo, F.; Kim, J.; Cokus, S.J.; et al. Relationship between nucleosome positioning and DNA methylation. Nature 2010, 466, 388–392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Shen, S.; Park, J.W.; Lu Z-x Lin, L.; Henry, M.D.; Wu, Y.N.; Zhou, Q.; Xing, Y. rMATS: Robust and flexible detection of differ-ential alternative splicing from replicate RNA-Seq data. Proc. Natl. Acad. Sci. USA 2014, 111, E5593–E5601. [Google Scholar] [CrossRef]
  32. Cortés, A.J.; Skeen, P.; Blair, M.W.; Chacón-Sánchez, M.I. Does the Genomic Landscape of Species Divergence in Phaseolus Beans Coerce Parallel Signatures of Adaptation and Domestication? Front. Plant Sci. 2018, 9, 1816. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Cortés, A.J.; Blair, M.W. Genotyping by Sequencing and Genome–Environment Associations in Wild Common Bean Predict Widespread Divergent Adaptation to Drought. Front. Plant Sci. 2018, 9, 128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Singh, P.; Ahi, E.P. The importance of alternative splicing in adaptive evolution. Mol. Ecol. 2022, 31, 1928–1938. [Google Scholar] [CrossRef]
  35. Wang, Z.; Hong, Y.; Yao, J.; Huang, H.; Qian, B.; Liu, X.; Chen, Y.; Pang, J.; Zhan, X.; Zhu, J.K.; et al. Modulation of plant development and chilling stress responses by alternative splicing events under control of the spliceosome protein SmEb in Arabidopsis. Plant Cell Env. 2022, 45, 2762–2779. [Google Scholar] [CrossRef] [PubMed]
  36. Liu, X.; Tian, Y.; Chi, W.; Zhang, H.; Yu, J.; Chen, G.; Wu, W.; Jiang, X.; Wang, S.; Lin, Z.; et al. Alternative splicing of OsGS1;1 affects nitrogen-use efficiency, grain development, and amylose content in rice. Plant J. 2022, 110, 1751–1762. [Google Scholar] [CrossRef]
  37. Zhu, J.; Yan, X.; Liu, S.; Xia, X.; An, Y.; Xu, Q.; Zhao, S.; Liu, L.; Guo, R.; Zhang, Z.; et al. Alternative splicing of CsJAZ1 negatively regulates flavan-3-ol biosynthesis in tea plants. Plant J. 2022, 110, 243–261. [Google Scholar] [CrossRef]
  38. Zhang, Q.; Zhang, W.; Wei, J.; Gao, Z.; Guan, J.; Cui, Z.; Wang, X. SKIP Regulates ABA Signaling through Alternative Splicing in Arabidopsis. Plant Cell Physiol. 2022, 63, 494–507. [Google Scholar] [CrossRef]
  39. Yan, X.; Bai, D.; Song, H.; Lin, K.; Pang, E. Alternative splicing during fruit development among fleshy fruits. BMC Genom. 2021, 22, 762. [Google Scholar] [CrossRef]
  40. Wang, Z.; Ji, H.; Yuan, B.; Wang, S.; Su, C.; Yao, B.; Zhao, H.; Li, X. ABA signalling is fine-tuned by antagonistic HAB1 variants. Nat. Commun. 2015, 6, 8138. [Google Scholar] [CrossRef] [Green Version]
  41. Penfield, S.; Josse, E.M.; Halliday, K.J. A role for an alternative splice variant of PIF6 in the control of Arabidopsis primary seed dormancy. Plant Mol. Biol. 2010, 73, 89–95. [Google Scholar] [CrossRef]
  42. Meyer, K.; Koester, T.; Staiger, D. Pre-mRNA Splicing in Plants: In Vivo Functions of RNA-Binding Proteins Implicated in the Splicing Process. Biomolecules 2015, 5, 1717–1740. [Google Scholar] [CrossRef] [Green Version]
  43. Wang, L.; Wang, L.; Tan, M.; Wang, L.; Zhao, W.; You, J.; Wang, L.; Yan, X.; Wang, W. The pattern of alternative splicing and DNA methylation alteration and their interaction in linseed (Linum usitatissimum L.) response to repeated drought stresses. Biol. Res. 2023, 56, 12. [Google Scholar] [CrossRef] [PubMed]
  44. Wang, Q.; Xu, J.; Pu, X.; Lv, H.; Liu, Y.; Ma, H.; Wu, F.; Wang, Q.; Feng, X.; Liu, T.; et al. Maize DNA Methylation in Response to Drought Stress Is Involved in Target Gene Expression and Alternative Splicing. Int. J. Mol. Sci. 2021, 22, 8285. [Google Scholar] [CrossRef]
  45. Yu, J.; Jung, S.; Cheng, C.-H.; Ficklin, S.P.; Lee, T.; Zheng, P.; Jones, D.; Percy, R.G.; Main, D. CottonGen: A genomics, genetics and breeding database for cotton research. Nucleic Acids Res. 2013, 42, D1229–D1236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef] [PubMed]
  47. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differ-entiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Anders, S.; Pyl, P.T.; Huber, W. HTSeq–A Python framework to work with high-throughput sequencing data. bioRxiv 2014, 2824. [Google Scholar] [CrossRef] [Green Version]
  50. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  51. Song, Q.; Guan, X.; Chen, Z.J. Dynamic Roles for Small RNAs and DNA Methylation during Ovule and Fiber Development in Allotetraploid Cotton. PLoS Genet. 2015, 11, e1005724. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Krueger, F.; Andrews, S.R. Bismark: A flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics 2011, 27, 1571–1572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Zhou, Q.; Lim, J.-Q.; Sung, W.-K.; Li, G. An integrated package for bisulfite DNA methylation data analysis with Indel-sensitive mapping. BMC Bioinform. 2019, 20, 47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijms 24 11812 g001 550
Figure 1. (A) Different types of AS events in G. hirsutum. (B) The bar chart displays the number of the AS events (A3SS, A5SS, MXE, RI, and SE) during cotton fiber development; AS events were statistically identified in ovule samples from nine periods (−3 DPA–25 DPA) and in fiber samples from four periods (10 DPA–25 DPA). The black arc represents the fiber samples, and the gray arc represents the ovule samples. At and Dt subgenomes of G. hirsutum were separately counted; Stage 1 Initiation: −3 DPA–5 DPA, Stage 2 fiber elongation period: 10 DPA–15 DPA, Stage 3 secondary wall thickening period (SCW): 20 DPA–25 DPA; skipped exon (SE), Alternative 5′ splicing site (A5SS), Alternative 3′ splicing site (A3SS), mutually exclusive exons (MXE), retained intron (RI)).
Figure 1. (A) Different types of AS events in G. hirsutum. (B) The bar chart displays the number of the AS events (A3SS, A5SS, MXE, RI, and SE) during cotton fiber development; AS events were statistically identified in ovule samples from nine periods (−3 DPA–25 DPA) and in fiber samples from four periods (10 DPA–25 DPA). The black arc represents the fiber samples, and the gray arc represents the ovule samples. At and Dt subgenomes of G. hirsutum were separately counted; Stage 1 Initiation: −3 DPA–5 DPA, Stage 2 fiber elongation period: 10 DPA–15 DPA, Stage 3 secondary wall thickening period (SCW): 20 DPA–25 DPA; skipped exon (SE), Alternative 5′ splicing site (A5SS), Alternative 3′ splicing site (A3SS), mutually exclusive exons (MXE), retained intron (RI)).
Ijms 24 11812 g001
Ijms 24 11812 g002 550
Figure 2. Distribution of the AS genes and AS events across the whole chromosome (A: AS events, B: AS genes, and C: gene densities).
Figure 2. Distribution of the AS genes and AS events across the whole chromosome (A: AS events, B: AS genes, and C: gene densities).
Ijms 24 11812 g002
Ijms 24 11812 g003 550
Figure 3. Statistical analysis of different AS events of the samples (skipped exon (SE), Alternative 5′ splicing site (A5SS), Alternative 3′ splicing site(A3SS), mutually exclusive exons (MXE), retained intron(RI)).
Figure 3. Statistical analysis of different AS events of the samples (skipped exon (SE), Alternative 5′ splicing site (A5SS), Alternative 3′ splicing site(A3SS), mutually exclusive exons (MXE), retained intron(RI)).
Ijms 24 11812 g003
Ijms 24 11812 g004 550
Figure 4. Number of AS genes and AS events between the At and the Dt subgenomes in G. hirsutum. This box plot represents the significant difference in the number of AS genes and AS events between At and Dt subgenomes. (*** p  <  0.001, based on Student’s t-test. Skipped exon (SE), Alternative 5′ splicing site (A5SS), Alternative 3′ splicing site(A3SS), mutually exclusive exons (MXE) and retained intron (RI)).
Figure 4. Number of AS genes and AS events between the At and the Dt subgenomes in G. hirsutum. This box plot represents the significant difference in the number of AS genes and AS events between At and Dt subgenomes. (*** p  <  0.001, based on Student’s t-test. Skipped exon (SE), Alternative 5′ splicing site (A5SS), Alternative 3′ splicing site(A3SS), mutually exclusive exons (MXE) and retained intron (RI)).
Ijms 24 11812 g004
Ijms 24 11812 g005 550
Figure 5. Conservative and specific AS genes (A) and AS events (B) of different cotton fiber development stages. This Upset plot shows the conservative and specific AS events detected in cotton fiber development stages: ovule initiation, ovule elongation, ovule secondary cell wall biosynthesis, fiber elongation, and fiber secondary cell wall biosynthesis, with the different points representing the AS genes and AS events occurring at each period, red dots representing the number of AS genes and AS events occurring at all five periods, blue dots representing the number of AS genes and AS events occurring at only a single period, and black dots representing the number of AS genes and AS events occurring at 2, 3 or 4 periods.
Figure 5. Conservative and specific AS genes (A) and AS events (B) of different cotton fiber development stages. This Upset plot shows the conservative and specific AS events detected in cotton fiber development stages: ovule initiation, ovule elongation, ovule secondary cell wall biosynthesis, fiber elongation, and fiber secondary cell wall biosynthesis, with the different points representing the AS genes and AS events occurring at each period, red dots representing the number of AS genes and AS events occurring at all five periods, blue dots representing the number of AS genes and AS events occurring at only a single period, and black dots representing the number of AS genes and AS events occurring at 2, 3 or 4 periods.
Ijms 24 11812 g005
Ijms 24 11812 g006 550
Figure 6. The expression of 108 splicing factors during fiber development in G. hirsutum.
Figure 6. The expression of 108 splicing factors during fiber development in G. hirsutum.
Ijms 24 11812 g006
Ijms 24 11812 g007 550
Figure 7. Statistical comparison analysis of the sequence characteristics of AS genes. (A): gene length, (B) number of exons, and (C) GC content of the genes. (a represents Normal genes, b represents Skipped exon (SE) genes and c represent retained intron (RI)genes).
Figure 7. Statistical comparison analysis of the sequence characteristics of AS genes. (A): gene length, (B) number of exons, and (C) GC content of the genes. (a represents Normal genes, b represents Skipped exon (SE) genes and c represent retained intron (RI)genes).
Ijms 24 11812 g007
Ijms 24 11812 g008 550
Figure 8. Comparison analysis of the GC content of the junction sites of AS.(*** p  <  0.001, based on Student’s t-test).
Figure 8. Comparison analysis of the GC content of the junction sites of AS.(*** p  <  0.001, based on Student’s t-test).
Ijms 24 11812 g008
Ijms 24 11812 g009 550
Figure 9. (A) Methylation level of the gene body region of AS genes (SE, RI). (B) Methylation levels of 50 bp upstream and downstream from the gene splicing sites of AS genes (SE, RI) exon. Skipped exon (SE), retained intron (RI).
Figure 9. (A) Methylation level of the gene body region of AS genes (SE, RI). (B) Methylation levels of 50 bp upstream and downstream from the gene splicing sites of AS genes (SE, RI) exon. Skipped exon (SE), retained intron (RI).
Ijms 24 11812 g009
Table
Table 1. Statistics on the number of AS genes in different periods (unit: pieces; skipped exon (SE), Alternative 5′ splicing site (A5SS), Alternative 3′ splicing site (A3SS), mutually exclusive exons (MXE) and retained intron (RI)).
Table 1. Statistics on the number of AS genes in different periods (unit: pieces; skipped exon (SE), Alternative 5′ splicing site (A5SS), Alternative 3′ splicing site (A3SS), mutually exclusive exons (MXE) and retained intron (RI)).
A3SSA5SSMXERISESUM
−3 DPA_Ovule43692343438196030728437
0 DPA_Ovule48322656563233936429360
1 DPA_Ovule43622386527213326148309
3 DPA_Ovule32531664394143117336354
5 DPA_Ovule45992514799217633749237
10 DPA_Ovule2268111447545313164605
15 DPA_Ovule41772206489318826458669
20 DPA_Ovule197794941737112404091
25 DPA_Ovule11026021418048052894
10 DPA_Fiber30401636494261216766849
15 DPA_Fiber30531528434121516756067
20 DPA_Fiber24021213392183313215420
25 DPA_Fiber10535081644115992363
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zheng, J.; Wen, S.; Yu, Z.; Luo, K.; Rong, J.; Ding, M. Alternative Splicing during Fiber Development in G. hirsutum. Int. J. Mol. Sci. 2023, 24, 11812. https://doi.org/10.3390/ijms241411812

AMA Style

Zheng J, Wen S, Yu Z, Luo K, Rong J, Ding M. Alternative Splicing during Fiber Development in G. hirsutum. International Journal of Molecular Sciences. 2023; 24(14):11812. https://doi.org/10.3390/ijms241411812

Chicago/Turabian Style

Zheng, Jing, Shuhan Wen, Zhipeng Yu, Keyan Luo, Junkang Rong, and Mingquan Ding. 2023. "Alternative Splicing during Fiber Development in G. hirsutum" International Journal of Molecular Sciences 24, no. 14: 11812. https://doi.org/10.3390/ijms241411812

APA Style

Zheng, J., Wen, S., Yu, Z., Luo, K., Rong, J., & Ding, M. (2023). Alternative Splicing during Fiber Development in G. hirsutum. International Journal of Molecular Sciences, 24(14), 11812. https://doi.org/10.3390/ijms241411812

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop