Next Article in Journal
Importance and Antimicrobial Resistance of Mycoplasma bovis in Clinical Respiratory Disease in Feedlot Calves
Previous Article in Journal
Myostatin and Related Factors Are Involved in Skeletal Muscle Protein Breakdown in Growing Broilers Exposed to Constant Heat Stress
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 11
Issue 5
10.3390/ani11051469
animals-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

Characterization of Alternative Splicing (AS) Events during Chicken (Gallus gallus) Male Germ-Line Stem Cell Differentiation with Single-Cell RNA-seq

by
Changhua Sun
1,2,3,
Kai Jin
1,3,
Qisheng Zuo
1,3,
Hongyan Sun
1,3,
Jiuzhou Song
4,
Yani Zhang
1,3,
Guohong Chen
1,3 and
Bichun Li
1,3,*
1
Key Laboratory of Animal Breeding Reproduction and Molecular Design for Jiangsu Province, College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China
2
Department of Food Technology, College of Biochemical Engineering, Yangzhou Polytechnic College, Yangzhou 225009, China
3
Joint International Research Laboratory of Agriculture and Agri-Product Safety of the Ministry of Education of China, Yangzhou University, Yangzhou 225009, China
4
Animal & Avian Sciences, University of Maryland, College Park, MD 20741, USA
*
Author to whom correspondence should be addressed.
Animals 2021, 11(5), 1469; https://0-doi-org.brum.beds.ac.uk/10.3390/ani11051469
Submission received: 12 April 2021 / Revised: 17 May 2021 / Accepted: 18 May 2021 / Published: 20 May 2021
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Studies have shown that alternative splicing (AS) has been utilized in a wide variety of life processes. However, there are very few studies on AS during germ cell development. In this study, we preliminarily investigated the variation of variable shear events during the formation of chicken germ cells through the RNA-seq data analysis of embryonic stem cells (ESCs), gonad PGCs (gPGCs), and spermatogonia stem cells (SSCs), and the critical AS mode for several crucial stage-specific genes, which were identified during germ cell development. The results of this study lay a theoretical foundation for further analysis of the regulation mechanism of key genes involved in germ cell formation.

Abstract

Alternative splicing (AS) is a ubiquitous, co-transcriptional, and post-transcriptional regulation mechanism during certain developmental processes, such as germ cell differentiation. A thorough understanding of germ cell differentiation will help us to open new avenues for avian reproduction, stem cell biology, and advances in medicines for human consumption. Here, based on single-cell RNA-seq, we characterized genome-wide AS events in manifold chicken male germ cells: embryonic stem cells (ESCs), gonad primordial germ cells (gPGCs), and spermatogonia stem cells (SSCs). A total of 38,494 AS events from 15,338 genes were detected in ESCs, with a total of 48,955 events from 14,783 genes and 49,900 events from 15,089 genes observed in gPGCs and SSCs, respectively. Moreover, this distribution of AS events suggests the diverse splicing feature of ESCs, gPGCs, and SSCs. Finally, several crucial stage-specific genes, such as NANOG, POU5F3, LIN28B, BMP4, STRA8, and LHX9, were identified in AS events that were transmitted in ESCs, gPGCs, and SSCs. The gene expression results of the RNA-seq data were validated by qRT-PCR. In summary, we provided a comprehensive atlas of the genome-wide scale of the AS event landscape in male chicken germ-line cells and presented its distribution for the first time. This research may someday improve treatment options for men suffering from male infertility.

1. Introduction

Alternative splicing (AS) is an important, co-transcriptional, and post-transcriptional regulatory mechanism that results in multiple mRNA and protein isoforms from a single gene. Protein isoforms show different or mutually antagonistically functional and structural characteristics, which play critical roles in cellular processes, such as cell signaling and differentiation [1,2,3]. Recent high-throughput analyses have shown their ubiquitousness in nematode worm genes (Caenorhabditis elegans, 71%), human genes (Homo sapiens, 95%), mouse genes (Mus musculus, 57%), cow genes (Bos taurus, 21%), etc. [4,5,6]. The classical AS event types include the alternative 5′ splice site (A5′), alternative 3′ splice site (A3′), exon-skipping (ES), intron retention (IR), and alternative exon (AE), which were involved in the mRNA transcription from one gene via the spliced site of pre-mRNA [7]. Therefore, AS is the fundamental cause of significantly fewer genes than protein species in eukaryotes. Generally, AS events are involved in many physiological processes, as well as cell differentiation, apoptosis, and development. Abnormal AS events will cause deviant processes by changing gene expression products [8,9,10,11]. Hence, studying AS in biological processes is vitally important and essential.
Germ cells are exclusive cells which are capable of bringing genetic information to the next generation and sustaining life. Their regulation processes are subtle and intricate during cell differentiation [12,13,14]. It is noted that the formation, migration, and differentiation of germ cells are current hotspots in animal genetics, breeding, and reproduction research [15,16,17]. Moreover, germ cells can be isolated, cultured, and substituted for embryonic stem cells (ESCs) in the treatment of diseases (such as infertility and organ regeneration) without ethical and immune rejection complications [18,19]. Especially in chickens (Gallus gallus), the unique migration pattern of primordial germ cells (PGCs) and advantages of embryonic development in ovum make germ cell isolation easier, and abundant germ cells could be applied to advanced stem cell research and transgenic animal production [20,21]. In addition, chickens are a classic avian experimental model, and they have been used extensively in developmental animal biological research [22,23]. For these reasons, understanding of chicken germ cell biology has not only accelerated the practical application of avian reproduction biology and fundamental stem cell research, but also assisted human medicine and medical applications (including various birth defects, germ cell tumor, and drug target screening). The utilization of these methods allows researchers and scientists to avoid the ethical constraints of other methods [24]. Hence, the study of chicken germ cell differentiation is an extremely urgent need.
Recent evidence suggests that AS events are an essential and vital factor in developmental biology and have been shown to play important roles in testis proteome diversification and male germ cell differentiation via numerous mechanisms [25,26,27]. However, there has been a scarcity of reports regarding AS events studied in chicken germ cells, and the distribution and features of these reports are still unclear. Now, single RNA-seq and relevant bioinformatic analysis methods are powerful and accurate tools to characterize AS events in an entire genome [28,29]. Therefore, we screened relevant AS events related to the differentiation of chicken ESCs, gPGCs, and spermatogonia stem cells (SSCs) through single-cell RNA-seq. The primary aim of our study was to determine the distribution of AS events that participate in the development and differentiation of male germ cells. We went further by analyzing AS event genes to construct an interaction network and identify key pathways in bioinformatics to sharpen our insight into reinterpreting specific pathways. Finally, we investigated AS events and expression patterns of key genes in ESCs, gPGCs, and SSCs. This study not only provides a comprehensive view of AS events in chickens, but also highlights novel insights into potential roles of AS events in the study of male germ cell differentiation.

2. Materials and Methods

2.1. Materials

Chicken eggs were collected from Rugao Yellow chickens (from the Poultry Institute, Chinese Academy of Agricultural Sciences, Yangzhou, China). All eggs were incubated at 37 °C and 75% relative humidity for 4.5 days or 18.5 days.

2.2. Cells and RNA-seq

The method of chicken ESC, gPGC, and SSC isolation referred to the previous report [30,31]. Briefly, ESCs was separated from newly laid fertilized eggs (blastoderm, E0, HH stage X) and digested with trypsin-EDTA (0.25%) for 3 min. gPGCs were isolated from eggs which had been incubated for 4.5 days (gonads, E4.5, HH stage 28–30). Genital ridges from embryos were collected and cut up, followed by digestion with trypsin-EDTA (0.25%). SSCs were isolated from eggs which had been incubated for 18.5 days (testis, E18.5, HH stage 44)). Testes from male embryos were collected and cut up, followed by the digestion with Collagenase I and trypsin-EDTA (0.25%). Isolated cells were filtered and prepared for further RNA extraction.
Cells were collected and then treated with a SMART-Seq v4 ultra low input RNA kit for sequencing (Takara, Dalian, China, 634850). The RNA from each sample was used for cDNA synthesis and sequencing. The RNA library pools for three kinds of cells were established following the protocol of Illumina Hiseq2000 PE150, and the experiments were performed at the Novogene Company (Beijing, China).

2.3. RNA-Seq Reads Mapping

Filtering and quality control checks of raw reads from RNA-seq were performed utilizing the Fast QC program (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/, accessed on 1 August 2019). The clean reads were mapped to the chicken reference genome (Gallus gallus, galGal4.0) using Tophat (version 2.1.1) [32]. The mapped reads were assembled into contigs with Cufflinks (version 2.2.1). The expression value for each gene was calculated and normalized to FPKM (fragments per kilobase of exon per million fragments mapped) values using Cuffnorm (version 2.2.1) in the Cufflinks package [33].

2.4. Identification, Classification, and Quantification of AS Events

After utilizing Cufflinks, the identification and quantification of AS events was accomplished using ASprofile [34]. Initially, the extract AS program compares all pairs of transcripts within a gene to determine exon–intron structural differences that indicate an AS event. The following types of events have been implemented: exon skipping (SKIP), cassette exons (MSKIP), alternative transcript start and termination (TSS, TTS), retention of single or multiple introns (IR, MIR), and alternative exons (AEs). Additionally, variations on these classes have allowed for some wiggle room at the boundaries of the surrounding features, and they have been reported and prefixed with an X. Therefore, we classified these AS events into 5 types: alternative 5′ splice Site (A5′, including TSS and XTSS), alternative 3′ splice Site (A3′, including TTS and XTTS), exon skip (ES, including SKIP, MSKIP, XSKIP, and XMSKIP), retained intron (RI, including IR, XMIR, XIR, and XMIR), and alternative exon (AE, including AE and XAE). Finally, the extract-as-fpkm program calculates the fragments per kilobase million (FPKM) of each event from those of transcripts harboring an event in each sample. The results can be used to observe trends in dynamics of AS profiles.

2.5. Intersections and Interaction Network Analysis of the AS Event Gene

The intersections (Upset plot and Veen plot) between different types of AS events were also investigated and illustrated by employing the UpSetR package (version 1.4.0, https://www.r-project.org/web/packages/UpSetR, accessed on 22 May 2019) of R software (version 3.6.2, https://www.r-project.org/, accessed on 12 December 2019). After selection by Veen plot, the genes of AS events were submitted to String (version 11.0, https://string-db.org/, accessed on 19 January 2019) for the protein–protein interaction (PPI) analysis. Cystoscope (version 3.7.2, https://cytoscape.org/, accessed on 5 September 2019) was then applied to illustrate and identify neighborhoods where proteins were densely connected in the PPI network.

2.6. GO Enrichment and KEGG Pathway Analysis of AS Events Gene

Gene ontology (GO) enrichment analysis was performed using DAVID (v6.8, https://david.ncifcrf.gov/, accessed on 1 September 2020) [35,36], and GO terms were considered as statistically significant with p-values < 0.05. GO terms have been classified by the GO database (ftp://ftp.ncbi.nih.gov/gene/DATA/gene2-go.gz, accessed on 10 May 2017), which is an international standard gene functional classification system. In addition, BinGO, which illustrates relationships among GO terms, was explained using Cytoscape 3.7.2 version [37]. The Kyoto Encyclopedia of Genes and Genomes (KEGG) helps researchers to better understand the biological functions of genes based on large-scale molecular datasets (http://www.genome.jp/kegg/, accessed on 10 May 2017). In this study, the enrichment of KEGG was analyzed by KOBAS 3.0 (http://kobas.cbi.pku.edu.cn/kobas3/annotate/, accessed on 20 September 2018) and a p-value < 0.05 was set as the cutoff for significantly enriched KEGG terms. Illustrations were performed through the ggplot2 package in R [38].

2.7. RNA Preparation and RT-qPCR

Isolated cells were homogenized using TRIzol, and the total RNA was isolated according to the manufacturer’s instructions (QIAGEN, Beijing, China, DP424). qRT-PCR was performed using the FastKing One-Step RT-PCR Kit with SYBR green (QIAGEN, Beijing, China, KR123), and the mRNA levels were determined by a CFX-Connect Real-time PCR detection system (BIO-RAD, Hercules, CA, USA, 7500fast). Quantification results were normalized to the relative housing keeping gene (GAPDH) using the 2-ΔΔCt method. The sequences of qRT-PCR primers are listed in Table 1.

3. Results

3.1. The AS Events during Male Chicken Germ Cell Differentiation

To study AS events in male chicken germ cell differentiation, we performed genome-wide RNA-seq for ESCs, gPGCs, and SSCs isolated from chicken embryonic blastoderms (E0, HH stage X), gonads (E4.5, HH stage 28–30), and testes (E18.5, HH stage 44), respectively (Figure 1a). Our results showed that the sequencing of the three cells transcriptome provided 195.7 million reads of ESCs, 170.1 million reads of gPGCs, and 143.2 million reads of SSCs. After quality control and mapping to the entire annotated chicken genome (galGal 4.0), 28.7%, 31.47% and 28.13% mapped splice reads were observed in ESCs, gPGCs, and SSCs, respectively (see Figure 1b). To refine and explore AS event changes during germ cell differentiation, we divided AS events into alternative 5′ splice site (A5′), alternative 3′ splice site (A3′), exon skip (ES), retained intron (RI), and alternative exon (AE) (see Figure 1c). In ESCs, a total of 38,494 AS events from 15,338 genes were detected, indicating that each gene might have had 2–3 AS events on average. However, 48,955 AS events from 14,783 genes and 49,900 AS events from 15,089 genes were found in gPGCs and SSCs, respectively. This indicated an average of 3–4 AS events from each gene, which was a notably higher frequency than that in ESCs (see Table 2). In detail, most AS events detected during germ cell differentiation were A5′ and A3′ followed by ES, RI, and AE, which further implies that gPGCs and SSCs had a higher frequency of AS events than ESCs across all five kinds of AS events, especially ES, RI, and AE (see Figure 1d). Moreover, ES, RI, and AE were much less prevalent than A5′/A3′ in three germ cell types, which suggests that A5′/A3′ are major AS events during chicken male germ cell differentiation.

3.2. AS Event Distribution Analysis during Male Chicken Germ Cell Differentiation

To investigate the splicing feature of each germ cell, we generalized the common gene distribution of AS events during male germ cell differentiation. The results showed that the A5′/A3′ gene group contained 12,453, 8081, and 9031 genes in ESCs, gPGCs, and SSCs, respectively. This was the intersection group with the largest number of genes involved in two types of events. However, the A5′/A3′/ES and AT groups contained 749 and 390 genes which followed the A5′/A3′ group in ESCs. In gPGCs, the following groups were the A5′/A3′ groups, which contained 1593 and 1338 genes, respectively. Additionally, SSCs, accounting for 1413 and 1175 genes, were contained in the A5′/A3′/ES group, and the A3′ group followed the A5′/A3′ group. In each type of male germ cell, different splicing features were present; more details can be found in Figure 2 and Tables S1–S3. Overall, the male germ cell exhibited different distribution characteristics of AS events, suggesting that AS is involved in the male germ cell differentiation.

3.3. A Cell-Type-Specific Spliced Gene in ESCs, GPGCs, and SSCs

To find out the potential influence of AS in male germ cell gene expressions, we compared the total gene and spliced gene expression patterns in ESCs, gPGCs, and SSCs. They shared a similar expression pattern (see Figure 3a; Tables S4 and S5). In addition to the expression pattern comparison at the gene level, we studied differentially spliced genes during the male germ cell differentiation and compared AS genes of three germ cells. Their overlap had 14,168 genes, and 222, 83, and 79 of these genes were singlehandedly present in ESCs, gPGCs, and SSCs, respectively (see Figure 3b). To determine whether the genes were affected by male germ cell differentiation, we analyzed these cell-specific gene groups. The network results revealed that ESCs might contain four interaction networks, which include hub genes at MYC, HOXA, RARS, etc. However, the GPGCs and SSCs may exhibit fewer networks and hub genes (see Figure 3c). Taken together, these results imply that differentially spliced genes appear in male germ cell differentiations.

3.4. Function Enrichment Results of Spliced Genes during Male Germ Cell Differentiation

To verify the function of these spliced genes, we performed GO and KEGG pathway enrichment analyses universally for the process of spliced genes during male germ cell differentiation. In total, 22 first-class GO terms were significantly identified: 7 GO terms were enriched in biological processes, including cellular process (GO:0009987), developmental process (GO:0032502), and metabolic process (GO:0008152). Additionally, we found 10 and 5 GO terms in the cellular component and molecular function parts, respectively (see Figure 4a, Table S6). Moreover, GO terms such as cellular process (GO:0009987), developmental process (GO:0032502), and metabolic process (GO:0008152) stand out, which play important roles in the biological process (Figure 4b). A total of 166 KEGG pathways were significantly enriched by the overall spliced genes (Table S7). The top 30 pathways are shown in Figure 4c, containing metabolic pathways (gga01100), endocytosis (gga04144), MAPK signaling pathway (gga04010), calcium signaling pathway (gga04218), Wnt signaling pathway (gga04310), and spliceosome (gga03040) (Figure 4c). These results indicate that spliced genes widely exist in the cell differentiation process.

3.5. Validation of AS Events and Crucial Gene Expression

Even previous results showed that spliced genes have similar expression patterns to total genes, but details of the relationship between AS events and spliced genes are still unknown. We displayed AS events from six genes which play a crucial role in male chicken germ cell differentiation. NANOG and POU5F3, highly expressed in ESCs, exhibited A5′ and A3′ events [39,40]. The LIN28B had more types of AS events in gPGCs, and both of the BMP4 had specific expression in gPGCs [41,42] (see Figure 5a). STRA8 and LHX9 exhibited high expressions in SSCs and showed different patterns of AS events during male germ cell differentiation [43,44] (see Figure 5b). Transcript levels of these genes in ESCs, gPGCs, and SSCs were measured by qRT-PCR. The consistent results from both qRT-PCR and RNA-Seq suggested the reliability of RNA-Seq data and its accuracy in quantifying gene expressions in male germ cells (see Figure 5c).

4. Discussion

The phenomenon of AS was discovered in the twentieth century, but it has not been systemically analyzed [45]. Advances in high-throughput sequencing provide an efficient technology to acquire AS events at the genome-wide level [46]. In this study, we reported and characterized genome-wide AS events during male chicken germ cell differentiation through a systemic comparison of single-cell RNA-seq. Results showed that AS events, especially A5′/A3′ events, are comprehensively involved in the process of male germ cell differentiation, while showing different patterns in ESCs, gPGCs, and SSCs. Functional enrichments in GO and KEGG pathways are closely related to cell differentiation. Moreover, some crucial gene expression patterns are differential AS events. Collectively, all these results provide a glimpse of dynamic involvements and systematical profiling of AS events in chicken ESCs, gPGCs, and SSCs, which might arouse the research use of numerous, potential biomarkers and AS events in the study of male germ cell differentiation.
A diverse splicing pattern in one gene leads to a variety of isoforms and expression patterns, which makes AS and its regulation mechanism more complex in the germ cellular process [47]. In this study, a total of 38,494 AS events out of 15,338 genes in ESCs, 48,955 AS events out of 14,783 genes in gPGCs, and 49,900 out of 15,089 genes in SSCs were detected. This is significant, in that AS is a commonly used process during male chicken germ cell differentiation. In addition, all cells generated the largest number of A5′ and A3′, while the smallest number of AEs was generated. It is possible that the A5′/A3′ is major factor in determining the transcript splicing, thus each alternative transcript is somewhat defined by the A5′/A3′. After comparing the difference of AS preference through ESCs, gPGCs, and SSCs, five types of AS events were increased, and ES, RI, and AE were all significantly enhanced in gPGCs and SSCs. The difference in AS has also been described in mammalian spermatogenesis, giving support to the agreement that changes in AS may contribute to male germ cell differentiation [48,49]. Moreover, the distribution of AS events demonstrated differential AS event patterns in ESCs, gPGCs, and SSCs. It is obvious that the splicing features of AS events may lead to direct differentiation.
To further understand the gene function affected by AS through influencing male chicken germ cell differentiation, we performed PPI, GO, and KEGG analyses of AS genes, and found that they were enriched in some biological processes, such as the developmental process, cellular process, and metabolic process. These results further suggest that AS events play an important role during male chicken germ cell differentiation. An interesting finding is that the metabolic process was not only enriched in GO terms, but also highly enriched in KEGG pathways. Similarly, many reports conclude that metabolic changes accompany the SSC differentiation [50,51]. We suspect that AS events induce the pathway transmission during male germ cell formation via metabolic transform. In addition, the calcium signing pathway, a key player in cell differentiation and development, was found to be regulated by some spliced factors in our study [52,53]. Being aware of many pathways, such as Wnt, MAPK, and mTOR, involved in the chicken germ cell formation process, we observed that abundant genes within these crucial pathways exhibited AS events. The implication is that AS events may affect germ cell differentiation via these pathways [54,55,56]. Moreover, germ cell-specific genes such as NANOG, POU5F3, LIN28B, BMP4, STRA8, and LHX9, which act as stage markers of germ cells, were developmentally regulated and stage-specific. In our results, the expression of these genes was associated with AS events transmitted through male chicken germ cell differentiation [39,40,41,42,43,44]. AS events may regulate the expression of these genes and subsequently impact the differentiation and developmental cell process. Taken together, AS is a complex and important regulatory mechanism, which may play a crucial role during male germ cell differentiation.

5. Conclusions

In summary, we have provided a comprehensive atlas on the genome-wide scale of the AS event landscape in male chicken germ-line cells (ESCs, gPGCs, and SSCs) and presented the distribution of these AS events. AS events have reportedly been involved in the germ cell differentiation process since they were first discovered. This helps us understand how AS events are associated with the male germ cell differentiation process, and to possibly improve treatment options for male infertility in animals and humans.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ani11051469/s1, Table S1: The AS event analysis in chicken ESCs, Table S2: The AS event analysis in chicken gPGCs, Table S3: The AS event analysis in chicken SSCs, Table S4: The total gene mRNA expression (FPKM) of ESCs, gPGCs, and SSCs, Table S5: The spliced gene mRNA expression (FPKM) of ESCs, gPGCs, and SSCs, Table S6: The GO term analysis of AS event genes during male chicken germ cell differentiation, Table S7: The KEGG pathway analysis of AS event genes during male chicken germ cell differentiation.

Author Contributions

Conceptualization, B.L. and G.C.; methodology, C.S.; validation, C.S. and K.J.; formal analysis, C.S. and K.J.; data curation, C.S.; writing—original draft prepara-tion, C.S.; writing, review, and editing, C.S., Q.Z., Y.Z., H.S., J.S., and K.J.; funding acquisition, B.L. and G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Key Research and Development Program (2017YFE0108000), National Natural Science Foundation of China (31972547), Qing Lan project of Jiangsu Province (2020), Excellent Educational Talents Program of the Fourth Talent Cultivation Plan in Yangzhou City (2020), Natural Science Project of Yangzhou Polytechnic College (2018ZR19), Young Academic Leaders Project of Yangzhou Polytechnic College (yzd2020-08), Jiangsu Science and Technology Project (Youth Fund): BK20180918.

Institutional Review Board Statement

Animal experiments were approved by the Institutional Animal Care and Use Committee of the Yangzhou University Animal Experiments Ethics Committee (permit number: SYXK [Su] IACUC 2012-0029). All experimental procedures were performed in accordance with the Regulations for the Administration of Affairs Concerning Experimental Animals approved by the State Council of the People’s Republic of China.

Data Availability Statement

The Illumina reads are available in the Sequence Read Archive (SRA) database at NCBI under study accession number PRJNA608148. All data are presented in the main manuscript and are available to readers.

Acknowledgments

The authors would like to thank Bichun Li and Jiuzhou Song for their excellent assistance in editing this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ASAlternative Splicing
ESCsEmbryonic Stem Cells
gPGCsPrimordial Germ Cells
SSCsSpermatogonia Stem Cells
APAlternative Promoter
ATAlternative Terminator
ESExon Skipping
IRIntron Retention
AEAlternative Exon
FPKMFragments Per Kilobase Million
GOGene Ontology
KEGGKyoto Encyclopedia of Genes and Genomes

References

  1. Nilsen, T.W.; Graveley, B.R. Expansion of the eukaryotic proteome by alternative splicing. Nature 2010, 463, 457–463. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Lee, Y.; Rio, D.C. Mechanisms and regulation of alternative Pre-mRNA splicing. Annu. Rev. Biochem. 2015, 84, 291–323. [Google Scholar] [CrossRef] [Green Version]
  3. Chen, M.; Manley, J.L. Mechanisms of alternative splicing regulation: Insights from molecular and genomics approaches. Nat. Rev. Mol. Cell Biol. 2009, 10, 741–754. [Google Scholar] [CrossRef]
  4. Ramani, A.K.; Calarco, J.A.; Pan, Q. Genome-wide analysis of alternative splicing in caenorhabditis elegans. Genome Res. 2011, 21, 342–348. [Google Scholar] [CrossRef] [Green Version]
  5. Pan, Q.; Shai, O.; Lee, L.J.; Frey, B.J.; Blencowe, B.J. Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat. Genet. 2008, 40, 1413–1415. [Google Scholar] [CrossRef]
  6. Chacko, E.; Ranganathan, S. Genome-wide analysis of alternative splicing in cow: Implications in bovine as a model for human diseases. BMC Genom. 2009, 10, S11. [Google Scholar] [CrossRef] [Green Version]
  7. Keren, H.; Lev-Maor, G.; Ast, G. Alternative splicing and evolution: Diversification, exon definition and function. Nat. Rev. Genet. 2010, 11, 345–355. [Google Scholar] [CrossRef] [PubMed]
  8. Yabas, M.; Elliott, H.; Hoyne, G.F. The role of alternative splicing in the control of immune homeostasis and cellular differentiation. Int. J. Mol. Sci. 2016, 17, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Xu, Y.; Zhao, W.; Olson, S.D.; Prabhakara, K.S.; Zhou, X. Alternative splicing links histone modifications to stem cell fate decision. Genome Biol. 2018, 19, 1–21. [Google Scholar] [CrossRef] [Green Version]
  10. Ramanouskaya, T.V.; Grinev, V.V. The determinants of alternative RNA splicing in human cells. Mol. Genet. Genom. 2017, 292, 1175–1195. [Google Scholar] [CrossRef]
  11. Fiszbein, A.; Kornblihtt, A.R. Alternative splicing switches: Important players in cell differentiation. Bioessays 2017, 39. [Google Scholar] [CrossRef]
  12. Spiller, C.M.; Burnet, G.; Bowles, J. Mouse fetal germ cell isolation and culture techniques. Methods Mol. Biol. 2017, 1463, 173–183. [Google Scholar] [CrossRef]
  13. Smallwood, S.A.; Kelsey, G. De novo DNA methylation: A germ cell perspective. Trends Genet. 2012, 28, 33–42. [Google Scholar] [CrossRef] [PubMed]
  14. Saitou, M.; Yamaji, M. Primordial germ cells in mice. Cold Spring Harb. Perspect. Biol. 2012, 4, a008375. [Google Scholar] [CrossRef] [Green Version]
  15. Yoshizaki, G.; Lee, S. Production of live fish derived from frozen germ cells via germ cell transplantation. Stem Cell Res. 2018, 29, 103–110. [Google Scholar] [CrossRef]
  16. Tagami, T.; Miyahara, D.; Nakamura, Y. Avian primordial germ cells. Adv. Exp. Med. Biol. 2017, 1001, 1–18. [Google Scholar] [CrossRef] [PubMed]
  17. Costa, J.; Souza, G.; Soares, M.; Ribeiro, R.; van den Hurk, R.; Silva, J.V. In vitro differentiation of primordial germ cells and oocyte-like cells from stem cells. Histol. Histopathol. 2018, 33, 121–132. [Google Scholar] [CrossRef]
  18. Huch, M.; Koo, B.-K. Modeling mouse and human development using organoid cultures. Development 2015, 142, 3113–3125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Nuti, N.; Corallo, C.; Chan, B.M.F.; Ferrari, M.; Gerami-Naini, B. Multipotent differentiation of human dental pulp stem cells: A literature Review. Stem Cell Rev. Rep. 2016, 12, 511–523. [Google Scholar] [CrossRef]
  20. Farzaneh, M.; Attari, F.; Mozdziak, P.E.; Khoshnam, S.E. The evolution of chicken stem cell culture methods. Br. Poult. Sci. 2017, 58, 681–686. [Google Scholar] [CrossRef] [Green Version]
  21. Kim, Y.M.; Han, J.Y. The early development of germ cells in chicken. Int. J. Dev. Biol. 2018, 62, 145–152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Davey, M.G.; Tickle, C. The chicken as a model for embryonic development. Cytogenet. Genome Res. 2007, 117, 231–239. [Google Scholar] [CrossRef]
  23. Johnsson, M.; Henriksen, R.; Höglund, A.; Fogelholm, J.; Jensen, P.; Wright, D. Genetical genomics of growth in a chicken model. BMC Genom. 2018, 19, 1–12. [Google Scholar] [CrossRef] [PubMed]
  24. Nowak-Sliwinska, P.; Segura, T.; Iruela-Arispe, M.L. The chicken chorioallantoic membrane model in biology, medicine and bioengineering. Angiogenesis 2014, 17, 779–804. [Google Scholar] [CrossRef] [Green Version]
  25. Hannigan, M.M.; Zagore, L.L.; Licatalosi, D.D. Ptbp2 controls an alternative splicing network required for cell communication during spermatogenesis. Cell Rep. 2017, 19, 2598–2612. [Google Scholar] [CrossRef] [Green Version]
  26. Zagore, L.L.; Grabinski, S.E.; Sweet, T.J. RNA binding protein Ptbp2 is essential for male germ cell development. Mol. Cell Biol. 2015, 35, 4030–4042. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Chen, Y.; Zheng, Y.; Gao, Y. Single-Cell RNA-Seq uncovers dynamic processes and critical regulators in mouse spermatogenesis. Cell Res. 2018, 28, 879–896. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Marinov, G.K.; Williams, B.A.; McCue, K. From single-cell to cell-pool transcriptomes: Stochasticity in gene expression and RNA splicing. Genome Res. 2014, 24, 496–510. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Chen, G.; Ning, B.; Shi, T. Single-cell RNA-Seq technologies and related computational data analysis. Front. Genet. 2019, 10, 317. [Google Scholar] [CrossRef] [PubMed]
  30. Zhang, Z.; Elsayed, A.K.; Shi, Q. Crucial genes and pathways in chicken germ stem cell differentiation. J. Biol. Chem. 2015, 290, 13605–13621. [Google Scholar] [CrossRef] [Green Version]
  31. Gao, W.; Zhang, C.; Jin, K.; Zhang, Y.; Zuo, Q.; Li, B. Analysis of LncRNA expression profile during the formation of male germ cells in chickens. Animals 2020, 10, 1850. [Google Scholar] [CrossRef] [PubMed]
  32. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, R36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Trapnell, C.; Roberts, A.; Goff, L. Differential gene and transcript expression analysis of RNA-Seq experiments with TopHat and cufflinks. Nat. Protoc. 2012, 7, 562–578. [Google Scholar] [CrossRef] [Green Version]
  34. Florea, L.; Song, L.; Salzberg, S.L. Thousands of exon skipping events differentiate among splicing patterns in sixteen human tissues. F1000Research 2013, 2, 188. [Google Scholar] [CrossRef] [PubMed]
  35. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37, 1–13. [Google Scholar] [CrossRef] [Green Version]
  36. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  37. Maere, S.; Heymans, K.; Kuiper, M. BiNGO: A cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 2005, 21, 3448–3449. [Google Scholar] [CrossRef] [Green Version]
  38. Wu, J.; Mao, X.; Cai, T.; Luo, J.; Wei, L. KOBAS server: A web-based platform for automated annotation and pathway identification. Nucleic Acids Res. 2006, 34, W720–W724. [Google Scholar] [CrossRef]
  39. Onichtchouk, D. Evolution and functions of Oct4 homologs in non-mammalian vertebrates. Biochim. Biophys. Acta 2016, 1859, 770–779. [Google Scholar] [CrossRef]
  40. Blinka, S.; Rao, S. Nanog expression in embryonic stem cells—An ideal model system to dissect enhancer function. Bioessays 2017, 39, 1700086. [Google Scholar] [CrossRef] [Green Version]
  41. Lawson, K.A.; Dunn, N.R.; Roelen, B.A. Bmp4 is required for the generation of primordial germ cells in the mouse embryo. Genes Dev. 1999, 13, 424–436. [Google Scholar] [CrossRef]
  42. West, J.A.; Viswanathan, S.R.; Yabuuchi, A.; Cunniff, K.; Takeuchi, A.; Park, I.-H.; Sero, J.E.; Zhu, H.; Perez-Atayde, A.; Frazier, A.L.; et al. A role for Lin28 in primordial germ-cell development and germ-cell malignancy. Nature 2009, 460, 909–913. [Google Scholar] [CrossRef] [Green Version]
  43. Yamazaki, F.; Møller, M.; Fu, C. The Lhx9 homeobox gene controls pineal gland development and prevents postnatal hydrocephalus. Brain Struct. Funct. 2015, 220, 1497–1509. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Zhang, Y.; Wang, Y.; Zuo, Q. CRISPR/Cas9 mediated chicken Stra8 gene knockout and inhibition of male germ cell differentiation. PLoS ONE 2017, 12, e0172207. [Google Scholar] [CrossRef] [PubMed]
  45. Baralle, F.E.; Giudice, J. Alternative splicing as a regulator of development and tissue identity. Nat. Rev. Mol. Cell Biol. 2017, 18, 437–451. [Google Scholar] [CrossRef]
  46. Zhao, S. Alternative splicing, RNA-Seq and drug discovery. Drug Discov. Today 2019, 24, 1258–1267. [Google Scholar] [CrossRef] [PubMed]
  47. Dvinge, H. Regulation of alternative MRNA splicing: Old players and new perspectives. FEBS Lett. 2018, 592, 2987–3006. [Google Scholar] [CrossRef]
  48. Song, H.; Wang, L.; Chen, D.; Li, F. The function of Pre-MRNA alternative splicing in mammal spermatogenesis. Int. J. Biol. Sci. 2020, 16, 38–48. [Google Scholar] [CrossRef]
  49. Guan, Y.; Liang, G.; Martin, G.B.; Guan, L.L. Functional changes in MRNA expression and alternative Pre-MRNA splicing associated with the effects of nutrition on apoptosis and spermatogenesis in the adult testis. BMC Genom. 2017, 18, 64. [Google Scholar] [CrossRef] [Green Version]
  50. Sun, R.; Qi, H. Dynamic expression of combinatorial replication-dependent histone variant genes during mouse spermatogenesis. Gene Expr. Patterns 2014, 14, 30–41. [Google Scholar] [CrossRef]
  51. Lord, T.; Nixon, B. Metabolic changes accompanying spermatogonial stem cell differentiation. Dev. Cell 2020, 52, 399–411. [Google Scholar] [CrossRef] [PubMed]
  52. Jaffe, L.F. Calcium explosions as triggers of development. Ann. N. Y. Acad. Sci. 1980, 339, 86–101. [Google Scholar] [CrossRef] [PubMed]
  53. Brodskiy, P.A.; Zartman, J.J. Calcium as a signal integrator in developing epithelial tissues. Phys. Biol. 2018, 15, 051001. [Google Scholar] [CrossRef] [PubMed]
  54. Lee, H.C.; Lim, S.; Han, J.Y. Wnt/β-catenin signaling pathway activation is required for proliferation of chicken primordial germ cells in vitro. Sci. Rep. 2016, 6, 34510. [Google Scholar] [CrossRef] [Green Version]
  55. Zhang, T.; Du, X.; Zhao, L. SIRT1 facilitates primordial follicle recruitment independent of deacetylase activity through directly modulating Akt1 and MTOR transcription. FASEB J. 2019, 33, 14703–14716. [Google Scholar] [CrossRef] [Green Version]
  56. Jouni, F.J.; Abdolmaleki, P.; Behmanesh, M.; Movahedin, M. An in vitro study of the impact of 4mT static magnetic field to modify the differentiation rate of rat bone marrow stem cells into primordial germ cells. Differentiation 2014, 87, 230–237. [Google Scholar] [CrossRef]
Animals 11 01469 g001 550
Figure 1. AS events during chicken male germ cell differentiation. (a) ESCs were isolated from the blastoderm of fertile eggs at E0 (HH stage X). gPGCs were isolated from chicken gonads at E4.5 (HH stage 28–30), and SSCs were isolated from chicken testis at E18.5 (HH stage 44). (b) The overview of RNA-seq data analysis pipeline. (c) Representative model for five types of AS events, including alternative 5′ splice site (A5′), alternative 3′ splice site (A3′), exon skip (ES), retained intron (RI), and alternative exon (AE). (d) Numbers of AS events and involved genes during male chicken germ cell differentiation.
Figure 1. AS events during chicken male germ cell differentiation. (a) ESCs were isolated from the blastoderm of fertile eggs at E0 (HH stage X). gPGCs were isolated from chicken gonads at E4.5 (HH stage 28–30), and SSCs were isolated from chicken testis at E18.5 (HH stage 44). (b) The overview of RNA-seq data analysis pipeline. (c) Representative model for five types of AS events, including alternative 5′ splice site (A5′), alternative 3′ splice site (A3′), exon skip (ES), retained intron (RI), and alternative exon (AE). (d) Numbers of AS events and involved genes during male chicken germ cell differentiation.
Animals 11 01469 g001
Animals 11 01469 g002 550
Figure 2. Distribution analysis of AS events during male chicken germ cell differentiation. The Veen and Upset plots show the distribution of AS events in chicken ESCs, gPGCs, and SSCs. The bar chart above represents the number of genes contained in each type of group. The bar chart at the bottom-left represents the number of events included in each type of AS event. The dotted line at the bottom-right shows the types of events contained in the group.
Figure 2. Distribution analysis of AS events during male chicken germ cell differentiation. The Veen and Upset plots show the distribution of AS events in chicken ESCs, gPGCs, and SSCs. The bar chart above represents the number of genes contained in each type of group. The bar chart at the bottom-left represents the number of events included in each type of AS event. The dotted line at the bottom-right shows the types of events contained in the group.
Animals 11 01469 g002
Animals 11 01469 g003 550
Figure 3. A cell-type-specific spliced gene in ESCs, gPGCs, and SSCs. (a) The expression patterns of total genes and spliced genes during male chicken germ cell differentiation. (b) Differentially spliced genes among three cell types. (c) The gene network of cell-specific spliced gene groups. The red dotted line is the ESC-specific spliced gene–gene interaction network; the yellow dotted line is the gPGC-specific spliced gene–gene interaction network; and the green dotted line is the gPGC-specific spliced gene–gene interaction network.
Figure 3. A cell-type-specific spliced gene in ESCs, gPGCs, and SSCs. (a) The expression patterns of total genes and spliced genes during male chicken germ cell differentiation. (b) Differentially spliced genes among three cell types. (c) The gene network of cell-specific spliced gene groups. The red dotted line is the ESC-specific spliced gene–gene interaction network; the yellow dotted line is the gPGC-specific spliced gene–gene interaction network; and the green dotted line is the gPGC-specific spliced gene–gene interaction network.
Animals 11 01469 g003
Animals 11 01469 g004 550
Figure 4. Functional enrichment analysis results of spliced genes during male germ cell differentiation. (a) The significantly enriched (p < 0.05) first-class GO terms cluster of AS genes that were alternatively spliced in all cell types. The red line represents the biological process terms; the blue line represents cellular component terms; and the green line represents molecular function terms. (b) The enriched GO terms network. Each term is represented by a circle node, where its size is proportional to the number of AS genes, and its color represents its cluster identity. (c) The top 30 enriched KEGG pathway of AS genes that were alternatively spliced in all cell types.
Figure 4. Functional enrichment analysis results of spliced genes during male germ cell differentiation. (a) The significantly enriched (p < 0.05) first-class GO terms cluster of AS genes that were alternatively spliced in all cell types. The red line represents the biological process terms; the blue line represents cellular component terms; and the green line represents molecular function terms. (b) The enriched GO terms network. Each term is represented by a circle node, where its size is proportional to the number of AS genes, and its color represents its cluster identity. (c) The top 30 enriched KEGG pathway of AS genes that were alternatively spliced in all cell types.
Animals 11 01469 g004
Animals 11 01469 g005 550
Figure 5. Validation of AS events and crucial gene expression. (a) The gene transcript structure and AS events of crucial genes in male chicken germ cell differentiation. The red square is exon, the different color cycles are AS events, and the purple arrow represents the qRT-PCR primers at the constitutive exon loci. (b) The gene expression of RNA-seq data. (c) The gene expression of qRT-PCR results and correlation analysis.
Figure 5. Validation of AS events and crucial gene expression. (a) The gene transcript structure and AS events of crucial genes in male chicken germ cell differentiation. The red square is exon, the different color cycles are AS events, and the purple arrow represents the qRT-PCR primers at the constitutive exon loci. (b) The gene expression of RNA-seq data. (c) The gene expression of qRT-PCR results and correlation analysis.
Animals 11 01469 g005
Table
Table 1. Sequences of qRT-PCR primers.
Table 1. Sequences of qRT-PCR primers.
GeneTranscript IDPrimer (5′-3′)Title 3
NANOGENSGALT00000044516F: TACCCCAGACTCTGCCACTA163 bp
R: CTTCTGGCTCTGAAACCGC
POU5F3ENSGALT00000043197F: GATGCGCCGACCTCAGAG112 bp
R: CATAGAGCGTGCCCAGAG
LIN28BENSGALT00000044077F: CAATGTGAGGATGGGCTTCG220 bp
R: ACTTCCTAAACAGGGGCTCC
BMP4ENSGALT00000020316F: GATCTCTACCGGCTCCAGTC174 bp
R: GTTGAAGACGAAGCGGATCC
STRA8ENSGALT00000038151F: TCCACGGCTATTTCACACCT242 bp
R: TCAAGGAAACCAGCAGCAAC
LHX9ENSGALT00000003480F: GAACTCACCTGCTTTGCCAA181 bp
R: GAGTCTTGTTGCAGGTGGTG
GAPDHENSGALT00000023323F: TGGGAAGCTGTGGAGAGATG166 bp
R: GCAGGTCAGGTCAACAACAG
Table
Table 2. The identification of AS events in chicken ESCs, gPGCs, and SSCs.
Table 2. The identification of AS events in chicken ESCs, gPGCs, and SSCs.
AS ClassificationESCsgPGCsSSCs
AS EventsGeneAS EventsGeneAS EventsGene
Alternative 5′ Splice Site (A5′)16,74114,56018,69712,63919,39113,402
Alternative 3′ Splice Site (A3′)16,27314,60917,21712,87417,51513,601
Exon Skip (ES)296414127002248772122605
Retained Intron (RI)14766493564127131861169
Alternative Exon (AE)1040486247599525961113
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, C.; Jin, K.; Zuo, Q.; Sun, H.; Song, J.; Zhang, Y.; Chen, G.; Li, B. Characterization of Alternative Splicing (AS) Events during Chicken (Gallus gallus) Male Germ-Line Stem Cell Differentiation with Single-Cell RNA-seq. Animals 2021, 11, 1469. https://0-doi-org.brum.beds.ac.uk/10.3390/ani11051469

AMA Style

Sun C, Jin K, Zuo Q, Sun H, Song J, Zhang Y, Chen G, Li B. Characterization of Alternative Splicing (AS) Events during Chicken (Gallus gallus) Male Germ-Line Stem Cell Differentiation with Single-Cell RNA-seq. Animals. 2021; 11(5):1469. https://0-doi-org.brum.beds.ac.uk/10.3390/ani11051469

Chicago/Turabian Style

Sun, Changhua, Kai Jin, Qisheng Zuo, Hongyan Sun, Jiuzhou Song, Yani Zhang, Guohong Chen, and Bichun Li. 2021. "Characterization of Alternative Splicing (AS) Events during Chicken (Gallus gallus) Male Germ-Line Stem Cell Differentiation with Single-Cell RNA-seq" Animals 11, no. 5: 1469. https://0-doi-org.brum.beds.ac.uk/10.3390/ani11051469

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