Transcriptome profiling of the murine testis during the first wave of spermatogenesis
Laiho A, Kotaja N, Gyenesei A, Sironen A., 17.04.2013
Abstract
Correct gene expression patterns form the basis for male germ cell differentiation and male fertility. Although previous studies have elucidated the importance of testis specific gene expression, the exact transcripts and comprehensive gene expression patterns remain unknown. Large scale sequencing techniques have enabled cost effective analysis of gene expression and isoform studies. Using the SOLiD 4 next-generation sequencing platform we have investigated the gene expression patterns at five different time points during the first wave on murine spermatogenesis. Our results highlight the upregulation of spermatogenesis related biological processes and associated cellular components. Elucidation of differential gene expression at important time points during the sperm development emphasizes the importance of correct timing of gene expression within biological processes. Differential gene level expression was analyzed with R/Bioconductor's Limma package and isoform analysis was conducted with the Cufflinks pipeline. At gene level total of 2494 differentially expressed genes were identified and Cufflinks characterized over 160,000 gene isoforms, of which 29% were novel transcripts assigned to known genes. Isoforms were detected for 57% of expressed genes and in a total over 26,000 genes were expressed in the testis. Differential promoter and transcription start site usage appears also to play a role in regulation of gene expression during spermatogenesis. Furthermore, we identified 947 upregulated long non-coding RNAs during the first wave of spermatogenesis. These RNAs appeared to be highly specific to different time points. Transcriptomic analysis of testis tissue samples is highly informative due to the large number of expressed genes and identified isoforms. Our study provides a very valuable basis for investigation of gene isoforms and regulation and factors contributing to male fertility.
Laiho A, Kotaja N, Gyenesei A, Sironen A. Transcriptome profiling of the murine testis during the first wave of spermatogenesis. PLoS One. 2013 Apr 17;8(4):e61558. doi: 10.1371/journal.pone.0061558. PMID: 23613874; PMCID: PMC3629203.
Publication: https://doi.org/10.1371/journal.pone.0061558 Repository: https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA326117
Disclaimer
The publication Transcriptome profiling of the murine testis during the first wave of spermatogenesis by Laiho A, Kotaja N, Gyenesei A, Sironen A. is published under an open access license: https://creativecommons.org/licenses/by/4.0/. Permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Curation by the MFGA team Relevant data sets presented in the publication have been identified. If possible, annotations (title, general information, conditions, processed tissue types and processed cell types) have been added based on information from the publication. Data tables and images that provide a good overview on the publication's findings on the data set have been extracted from the publication and/or supplement. If not stated otherwise, images are depicted with title and description exactly as in the publication. Tables have been adjusted to the MFGA table format. Conducted adjustments are explained in the detailed view of the tables. However, titles and descriptions have been adopted from the publication.
Data set 1: Genome‑wide DNA methylation levels in different cell stages
Methylome: Whole Genome Bisulfite Sequencing
Tissue Types
| BRENDA tissue ontology | Maturity | Description | Species | Replicates |
|---|---|---|---|---|
| BTO_0001363: testis | A typically paired male reproductive gland that produces sperm and that in most mammals is contained within the scrotum at sexual maturity. | Mouse | 2 |
Cell Types
| Cell ontology | Maturity | Description | Species | Replicates | Cells per replicate |
|---|---|---|---|---|---|
| CL_0000017: spermatocyte | A male germ cell that develops from spermatogonia. The euploid primary spermatocytes undergo meiosis and give rise to the haploid secondary spermatocytes which in turn give rise to spermatids. | Mouse | 2 | ||
| CL_0000020: spermatogonium | An euploid male germ cell of an early stage of spermatogenesis. | Mouse | 2 | ||
| CL_0000018: spermatid | A male germ cell that develops from the haploid secondary spermatocytes. Without further division, spermatids undergo structural changes and give rise to spermatozoa. | Mouse | 2 | ||
| -: epididymal spermatozoa | epididymal spermatozoa | Mouse | 2 |
Images
Figure 1: Clustering of sample groups and overlapping DE genes.
A. Hierarchical clustering for the normalized data. Biological replicate samples cluster closely together. Sample groups PND28 and PND7 are more distant to other groups and PND14 and PND17 are similar to each other based on the gene expression differences. B. Total number of overlapping DE genes between comparisons. C. Agreeing and overlapping upregulated genes between comparisons. D. Agreeing and overlapping downregulated genes between comparisons.
Licensed under: https://creativecommons.org/licenses/by/4.0/
Figure 2: Highest ranked GO terms (p<0.01) with differential expression between sample groups for GO biological processes (A) and cellular components (B) and for KEGG pathways (C).
The relative significance of term enrichment is indicated by levels of red shading (the lower the p-value the darker shade of red colour). Heat map was produced with GENE-E (http://www.broadinstitute.org/cancer/software/GENE-E/) using pearson’s correlation metrics and average linkage distance metrics.
Licensed under: https://creativecommons.org/licenses/by/4.0/
Figure 3: Enriched biological process and cellular component GO terms among upregulated genes.
Comparison of PND14/17 is not presented due to low number of DE genes. A. PND 7/14 biological processes, B. PND 7/14 cellular component, C. PND 17/21 biological processes, D. PND 17/21 cellular component, E. PND 21/28 biological processes and F. PND 21/28 cellular component. Only lower level GO terms are shown.
Licensed under: https://creativecommons.org/licenses/by/4.0/
Figure 4: Venn diagram and enriched cellular complex GO terms of upregulated genes in biological process GO term gamete generation.
Only lower level GO terms are shown.
Licensed under: https://creativecommons.org/licenses/by/4.0/
Figure 5: Gene and isoform counts and promoter and transcription start site usage in sample groups.
A. Expressed gene, known and novel isoform counts at different time points during the first wave of spermatogenesis. FPKM >1. B. Differential promoter usage in all sample group comparisons. C. Differential transcription start site usage between all sample groups. The comparisons in chronological order during the first wave of spermatogenesis are indicated in black (B and C).
Licensed under: https://creativecommons.org/licenses/by/4.0/
Figure 6: Enriched biological process (A) and cellular component (B) GO terms of genes with differential promoter usage between sample groups (p<10−6).
Only lower level terms are shown.
Licensed under: https://creativecommons.org/licenses/by/4.0/