Dynamics of the transcriptional landscape during human fetal testis and ovary development

STUDY QUESTION

Which transcriptional program triggers sex differentiation in bipotential gonads and downstream cellular events governing fetal testis and ovary development in humans?

SUMMARY ANSWER

The characterization of a dynamically regulated protein-coding and non-coding transcriptional landscape in developing human gonads of both sexes highlights a large number of potential key regulators that show an early sexually dimorphic expression pattern.

WHAT IS KNOWN ALREADY

Gonadal sex differentiation is orchestrated by a sexually dimorphic gene expression program in XX and XY developing fetal gonads. A comprehensive characterization of its non-coding counterpart offers promising perspectives for deciphering the molecular events underpinning gonad development and for a complete understanding of the etiology of disorders of sex development in humans.

STUDY DESIGN, SIZE, DURATION

To further investigate the protein-coding and non-coding transcriptional landscape during gonad differentiation, we used RNA-sequencing (RNA-seq) and characterized the RNA content of human fetal testis (N = 24) and ovaries (N = 24) from 6 to 17 postconceptional week (PCW), a key period in sex determination and gonad development.

PARTICIPANTS/MATERIALS, SETTING, METHODS

First trimester fetuses (6-12 PCW) and second trimester fetuses (13-14 and 17 PCW) were obtained from legally induced normally progressing terminations of pregnancy. Total RNA was extracted from whole human fetal gonads and sequenced as paired-end 2 × 50 base reads. Resulting sequences were mapped to the human genome, allowing for the assembly and quantification of corresponding transcripts.

MAIN RESULTS AND THE ROLE OF CHANCE

This RNA-seq analysis of human fetal testes and ovaries at seven key developmental stages led to the reconstruction of 22 080 transcripts differentially expressed during testicular and/or ovarian development. In addition to 8935 transcripts displaying sex-independent differential expression during gonad development, the comparison of testes and ovaries enabled the discrimination of 13 145 transcripts that show a sexually dimorphic expression profile. The latter include 1479 transcripts differentially expressed as early as 6 PCW, including 39 transcription factors, 40 long non-coding RNAs and 20 novel genes. Despite the use of stringent filtration criteria (expression cut-off of at least 1 fragment per kilobase of exon model per million reads mapped, fold change of at least 2 and false discovery rate adjusted P values of less than <1%), the possibility of assembly artifacts and of false-positive differentially expressed transcripts cannot be fully ruled out.

LARGE-SCALE DATA

Raw data files (fastq) and a searchable table (.xlss) containing information on genomic features and expression data for all refined transcripts have been submitted to the NCBI GEO under accession number GSE116278.

LIMITATIONS, REASONS FOR CAUTION

The intrinsic nature of this bulk analysis, i.e. the sequencing of transcripts from whole gonads, does not allow direct identification of the cellular origin(s) of the transcripts characterized. Potential cellular dilution effects (e.g. as a result of distinct proliferation rates in XX and XY gonads) may account for a few of the expression profiles identified as being sexually dimorphic. Finally, transcriptome alterations that would result from exposure to pre-abortive drugs cannot be completely excluded. Although we demonstrated the high quality of the sorted cell populations used for experimental validations using quantitative RT-PCR, it cannot be totally excluded that some germline expression may correspond to cell contamination by, for example, macrophages.

WIDER IMPLICATIONS OF THE FINDINGS

For the first time, this study has led to the identification of 1000 protein-coding and non-coding candidate genes showing an early, sexually dimorphic, expression pattern that have not previously been associated with sex differentiation. Collectively, these results increase our understanding of gonad development in humans, and contribute significantly to the identification of new candidate genes involved in fetal gonad differentiation. The results also provide a unique resource that may improve our understanding of the fetal origin of testicular and ovarian dysgenesis syndromes, including cryptorchidism and testicular cancers.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the French National Institute of Health and Medical Research (Inserm), the University of Rennes 1, the French School of Public Health (EHESP), the Swiss National Science Foundation [SNF n° CRS115_171007 to B.J.], the French National Research Agency [ANR n° 16-CE14-0017-02 and n° 18-CE14-0038-02 to F.C.], the Medical Research Council [MR/L010011/1 to P.A.F.] and the European Community's Seventh Framework Programme (FP7/2007-2013) [under grant agreement no 212885 to P.A.F.] and from the European Union's Horizon 2020 Research and Innovation Programme [under grant agreement no 825100 to P.A.F. and S.M.G.]. There are no competing interests related to this study.

RNA profiling of human testicular cells identifies syntenic lncRNAs associated with spermatogenesis

STUDY QUESTION

Is the noncoding transcriptional landscape during spermatogenesis conserved between human and rodents?

SUMMARY ANSWER

We identified a core group of 113 long noncoding RNAs (lncRNAs) and 20 novel genes dynamically and syntenically transcribed during spermatogenesis.

WHAT IS KNOWN ALREADY

Spermatogenesis is a complex differentiation process driven by a tightly regulated and highly specific gene expression program. Recently, several studies in various species have established that a large proportion of known lncRNAs are preferentially expressed during meiosis and spermiogenesis in a testis-specific manner.

STUDY DESIGN, SIZE, DURATION

To further investigate lncRNA expression in human spermatogenesis, we carried out a cross-species RNA profiling study using isolated testicular cells.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Human testes were obtained from post-mortem donors (N = 8, 51 years old on average) or from prostate cancer patients with no hormonal treatment (N = 9, 80 years old on average) and only patients with full spermatogenesis were used to prepare enriched populations of spermatocytes, spermatids, Leydig cells, peritubular cells and Sertoli cells. To minimize potential biases linked to inter-patient variations, RNAs from two or three donors were pooled prior to RNA-sequencing (paired-end, strand-specific). Resulting reads were mapped to the human genome, allowing for assembly and quantification of corresponding transcripts.

MAIN RESULTS AND THE ROLE OF CHANCE

Our RNA-sequencing analysis of pools of isolated human testicular cells enabled us to reconstruct over 25 000 transcripts. Among them we identified thousands of lncRNAs, as well as many previously unidentified genes (novel unannotated transcripts) that share many properties of lncRNAs. Of note is that although noncoding genes showed much lower synteny than protein-coding ones, a significant fraction of syntenic lncRNAs displayed conserved expression during spermatogenesis.

LARGE SCALE DATA

Raw data files (fastq) and a searchable table (.xlss) containing information on genomic features and expression data for all refined transcripts have been submitted to the NCBI Gene Expression Omnibus under accession number GSE74896.

LIMITATIONS, REASONS FOR CAUTION

Isolation procedures may alter the physiological state of testicular cells, especially for somatic cells, leading to substantial changes at the transcriptome level. We therefore cross-validated our findings with three previously published transcriptomic analyses of human spermatogenesis. Despite the use of stringent filtration criteria, i.e. expression cut-off of at least three fragments per kilobase of exon model per million reads mapped, fold-change of at least three and false discovery rate adjusted P-values of less than <1%, the possibility of assembly artifacts and false-positive transcripts cannot be fully ruled out.

WIDER IMPLICATIONS OF THE FINDINGS

For the first time, this study has led to the identification of a large number of conserved germline-associated lncRNAs that are potentially important for spermatogenesis and sexual reproduction. In addition to further substantiating the basis of the human testicular physiology, our study provides new candidate genes for male infertility of genetic origin. This is likely to be relevant for identifying interesting diagnostic and prognostic biomarkers and also potential novel therapeutic targets for male contraception.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by l'Institut national de la santé et de la recherche médicale (Inserm); l'Université de Rennes 1; l'Ecole des hautes études en santé publique (EHESP); INERIS-STORM to B.J. [N 10028NN]; Rennes Métropole 'Défis scientifiques émergents' to F.C (2011) and A.D.R (2013). The authors have no competing financial interests.