The Annotation, Mapping, Expression and Network (AMEN) suite of tools for molecular systems biology

Two thyroperoxidase-inhibiting chemicals induce shared transcriptional changes in hippocampus of developing rats

Thyroid hormone (TH) system disrupting compounds can impair brain development by perturbing TH action during critical life stages. Human exposure to TH system disrupting chemicals is therefore of great concern. To better protect humans against such chemicals, sensitive test methods that can detect effects on the developing brain are critical. Worryingly, however, current test methods are not sensitive and specific towards TH-mediated effects. To address this shortcoming, we performed RNA-sequencing of rat brains developmentally exposed to two different thyroperoxidase (TPO) inhibiting compounds, the medical drug methimazole (MMI) or the pesticide amitrole. Pregnant and lactating rats were exposed to 8 and 16 mg/kg/day(d) MMI or 25 and 50 mg/kg/d amitrole from gestational day 7 until postnatal day 16. Bulk-RNA-seq was performed on hippocampus from the 16-day old male pups. MMI and amitrole caused pronounced changes to the transcriptomes; 816 genes were differentially expressed, and 425 gene transcripts were similarly affected by both chemicals. Functional terms indicate effects from key cellular functions to changes in cell development, migration and differentiation of several cell populations. Of the total number of DEGs, 106 appeared to form a consistent transcriptional fingerprint of developmental hypothyroidism as they were similarly and dose-dependently expressed across all treatment groups. Using a filtering system, we identified 20 genes that appeared to represent the most sensitive, robust and dose-dependent markers of altered TH-mediated brain development. These markers provide inputs to the adverse outcome pathway (AOP) framework where they, in the context of linking TPO inhibiting compounds to adverse cognitive function, can be used to assess altered gene expression in the hippocampus in rat toxicity studies.

Transcriptional landscape of human keratinocyte models exposed to 60-GHz millimeter-waves

The use of millimeter waves (MMW) will exponentially grow in the coming years due to their future utilization in 5G/6G networks. The question of possible biological effects at these frequencies has been raised. In this present study, we aimed to investigate gene expression changes under exposure to MMW using the Bulk RNA Barcoding and sequencing (BRB-seq) technology. To address this issue, three exposure scenarios were performed aiming at: i) comparing the cellular response of two primary culture of keratinocytes (HEK and NHEK) and one keratinocyte derivate cell line (HaCaT) exposed to MMW; ii) exploring the incident power density dose-effect on gene expression in HaCaT cell line; and, iii) studying the exposure duration at the new ICNIRP exposure limit for the general population. With the exception of heat effect induced by high power MMW (over 10 mW/cm2), those exposure scenarios have not enabled us to demonstrate important gene expression changes in the different cell populations studied. Very few differentially genes were observed between MMW exposed samples and heat shock control, and most of them were significantly associated with heat shock response that may reflect small differences in the heat generation. Together these results show that acute exposure to MMW has no effects on the transcriptional landscape of human keratinocyte models under athermal conditions.

Sperm DNA methylation dynamics after chemotherapy: a longitudinal study of a patient with testicular germ cell tumor treatment

BACKGROUND

An important issue for young men affected by testicular germ cell tumor (TGCT) is how TGCT and its treatment will affect, transiently or permanently, their future reproductive health. Previous studies have reported that xenobiotics can induce changes on human sperm epigenome and have the potential to promote epigenetic alterations in the offspring.

OBJECTIVES

Here, we report the first longitudinal DNA methylation profiling of frozen sperm from a TGCT patient before and up to 2 years after a bleomycin, etoposide, and cisplatin (BEP) chemotherapy.

MATERIALS AND METHODS

A TGCT was diagnosed in a 30-year-old patient. A cryopreservation of spermatozoa was proposed before adjuvant BEP treatment. Semen samples were collected before and after chemotherapy at 6, 9, 12, and 24 months. The DNA methylation status was determined by RRBS to detect DNA differentially methylated regions (DMRs).

RESULTS

The analysis revealed that among the 74 DMRs showing modified methylation status 6 months after therapy, 17 remained altered 24 months after treatment. We next associated DMRs with differentially methylated genes (DMGs), which were subsequently intersected with loci known to be important or expressed during early development.

DISCUSSION AND CONCLUSION

The consequences of the cancer treatment on the sperm epigenome during the recovery periods are topical issues of increasing significance as epigenetic modifications to the paternal genome may have deleterious effects on the offspring. The altered methylated status of these DMGs important for early development might modify their expression pattern and thus affect their function during key stages of embryogenesis, potentially leading to developmental disorders or miscarriages.

The -KTS splice variant of WT1 is essential for ovarian determination in mice

Sex determination in mammals depends on the differentiation of the supporting lineage of the gonads into Sertoli or pregranulosa cells that govern testis and ovary development, respectively. Although the Y-linked testis-determining gene Sry has been identified, the ovarian-determining factor remains unknown. In this study, we identified -KTS, a major, alternatively spliced isoform of the Wilms tumor suppressor WT1, as a key determinant of female sex determination. Loss of -KTS variants blocked gonadal differentiation in mice, whereas increased expression, as found in Frasier syndrome, induced precocious differentiation of ovaries independently of their genetic sex. In XY embryos, this antagonized Sry expression, resulting in male-to-female sex reversal. Our results identify -KTS as an ovarian-determining factor and demonstrate that its time of activation is critical in gonadal sex differentiation.

Multimodal cartography of human lymphopoiesis reveals B and T/NK/ILC lineages are subjected to differential regulation

The developmental cartography of human lymphopoiesis remains incompletely understood. Here, we establish a multimodal map demonstrating that lymphoid specification follows independent direct or stepwise hierarchic routes converging toward the emergence of newly characterized CD117lo multi-lymphoid progenitors (MLPs) that undergo a proliferation arrest before entering the CD127- (NK/ILC/T) or CD127+ (B) lymphoid pathways. While the differentiation of CD127- early lymphoid progenitors is mainly driven by Flt3 signaling, emergence of their CD127+ counterparts is regulated cell-intrinsically and depends exclusively on the divisional history of their upstream precursors, including hematopoietic stem cells. Further, transcriptional mapping of differentiation trajectories reveals that whereas myeloid granulomonocytic lineages follow continuous differentiation pathways, lymphoid trajectories are intrinsically discontinuous and characterized by sequential waves of cell proliferation allowing pre-commitment amplification of lymphoid progenitor pools. Besides identifying new lymphoid specification pathways and regulatory checkpoints, our results demonstrate that NK/ILC/T and B lineages are under fundamentally distinct modes of regulation. (149 words).

Antiepileptic drugs are endocrine disruptors for the human fetal testis ex vivo

Valproic acid (VPA) has long been the most widely used antiepileptic drug (AED) for the treatment of epilepsy, bipolar psychiatric disorders, and migraine. However, long-term VPA treatment has several adverse effects on the male reproductive system notably on endocrine functions and/or spermatic parameters. In utero exposure of the fetus to VPA is well known to be associated with a higher risk of several congenital malformations including those of male reproductive organs. Subsequent generations of AEDs, such as carbamazepine (CARB) and lamotrigine (LAM), are considered safer and are currently recommended for women of child-bearing age with epilepsy. Because anomalies of the male genital tract mostly result from endocrine imbalance during fetal life, we hypothesized that AEDs could directly impair testis differentiation. We thus aimed at identifying and characterizing the effects of VPA, CARB, and LAM on the differentiation and function of the different testicular cell types, and at understanding the mechanisms underlying these effects. By using ex vivo culture of first-trimester human fetal testes, we show that VPA induces multiple endocrine disruptive effects, compared with the milder ones caused by CARB and LAM. AED also subtly altered the germ cell lineage in distinct manners. Transcriptomic analysis of VPA-induced alterations highlighted a very broad range of effects on the fetal testis. Overall, our results show that AEDs can behave as endocrine disruptors for the human fetal testis ex vivo. This is consistent with, and likely underlies, the VPA-induced male genital tract masculinization abnormalities observed in patients.

Expression of the endocannabinoid system and response to cannabinoid components by the human fetal testis

BACKGROUND

Cannabis consumption by pregnant women continues to increase worldwide, raising concerns about adverse effects on fetal growth and deleterious impacts on the newborn, in connection with evidence of placental transfer of cannabis compound. Cannabis action is mediated by the endocannabinoid system (ECS), which expression is well established in the brain but unknown in the developing testis. The fetal testis, whose endocrine function orchestrates the masculinization of many distant organs, is particularly sensitive to disruption by xenobiotics. In this context, we aimed to determine whether cannabis exposure has the potential to directly impact the human fetal testis.

METHODS

We determined the expression of components of the ECS in the human fetal testis from 6 to 17 developmental weeks and assessed the direct effects of phytocannabinoids Δ9-trans-tetrahydrocannabinol (THC) and cannabidiol (CBD) on the testis morphology and cell functions ex vivo.

RESULTS

We demonstrate the presence in the human fetal testis of two key endocannabinoids, 2-arachidonylglycerol (2-AG) and to a lower level anandamide (AEA), as well as a range of enzymes and receptors for the ECS. Ex vivo exposure of first trimester testes to CBD, THC, or CBD/THC [ratio 1:1] at 10^-7 to 10^-5 M altered testosterone secretion by Leydig cells, AMH secretion by Sertoli cells, and impacted testicular cell proliferation and viability as early as 72 h post-exposure. Transcriptomic analysis on 72 h-exposed fetal testis explants revealed 187 differentially expressed genes (DEGs), including genes involved in steroid synthesis and toxic substance response. Depending on the molecules and testis age, highly deleterious effects of phytocannabinoid exposure were observed on testis tissue after 14 days, including Sertoli and germ cell death.

CONCLUSIONS

Our study is the first to evidence the presence of the ECS in the human fetal testis and to highlight the potential adverse effect of cannabis consumption by pregnant women onto the development of the male gonad.

Distinct subsets of multi-lymphoid progenitors support ontogeny-related changes in human lymphopoiesis

Changes in lymphocyte production patterns occurring across human ontogeny remain poorly defined. In this study, we demonstrate that human lymphopoiesis is supported by three waves of embryonic, fetal, and postnatal multi-lymphoid progenitors (MLPs) differing in CD7 and CD10 expression and their output of CD127-/+ early lymphoid progenitors (ELPs). In addition, our results reveal that, like the fetal-to-adult switch in erythropoiesis, transition to postnatal life coincides with a shift from multilineage to B lineage-biased lymphopoiesis and an increase in production of CD127+ ELPs, which persists until puberty. A further developmental transition is observed in elderly individuals whereby B cell differentiation bypasses the CD127+ compartment and branches directly from CD10+ MLPs. Functional analyses indicate that these changes are determined at the level of hematopoietic stem cells. These findings provide insights for understanding identity and function of human MLPs and the establishment and maintenance of adaptative immunity.

Transcriptional profiling of the developing rat ovary following intrauterine exposure to the endocrine disruptors diethylstilbestrol and ketoconazole

Exposure to endocrine-disrupting chemicals (EDCs) during development may cause reproductive disorders in women. Although female reproductive endpoints are assessed in rodent toxicity studies, a concern is that typical endpoints are not sensitive enough to detect chemicals of concern to human health. If so, measured endpoints must be improved or new biomarkers of effects included. Herein, we have characterized the dynamic transcriptional landscape of developing rat ovaries exposed to two well-known EDCs, diethylstilbestrol (DES) and ketoconazole (KTZ), by 3' RNA sequencing. Rats were orally exposed from day 7 of gestation until birth, and from postnatal day 1 until days 6, 14 or 22. Three exposure doses for each chemical were used: 3, 6 and 12 µg/kg bw/day of DES; 3, 6, 12 mg/kg bw/day of KTZ. The transcriptome changed dynamically during perinatal development in control ovaries, with 1137 differentially expressed genes (DEGs) partitioned into 3 broad expression patterns. A cross-species deconvolution strategy based on a mouse ovary developmental cell atlas was used to map any changes to ovarian cellularity across the perinatal period to allow for characterization of actual changes to gene transcript levels. A total of 184 DEGs were observed across dose groups and developmental stages in DES-exposed ovaries, and 111 DEGs in KTZ-exposed ovaries across dose groups and developmental stages. Based on our analyses, we have identified new candidate biomarkers for female reproductive toxicity induced by EDC, including Kcne2, Calb2 and Insl3.

Throughout in vitro first spermatogenic wave: Next-generation sequencing gene expression patterns of fresh and cryopreserved prepubertal mice testicular tissue explants

INTRODUCTION

Suitable cryopreservation procedures of pre-pubertal testicular tissue associated with efficient culture conditions are crucial in the fields of fertility preservation and restoration. In vitro spermatogenesis remains a challenging technical procedure to undergo a complete spermatogenesis.The number of haploid cells and more specifically the spermatic yield produced in vitro in mice is still extremely low compared to age-matched in vivo controls and this procedure has never yet been successfully transferred to humans.

METHODS

To evaluate the impact of in vitro culture and freezing procedure, pre-pubertal testicular mice testes were directly cultured until day 4 (D4), D16 and D30 or cryopreserved by controlled slow freezing then cultured until D30. Testes composed of a panel of 6.5 dpp (days postpartum), 10.5 dpp, 22.5 dpp, and 36.5 dpp mice were used as in vivo controls. Testicular tissues were assessed by histological (HES) and immunofluorescence (stimulated by retinoic acid gene 8, STRA8) analyses. Moreover, a detailed transcriptome evaluation study has been carried out to study the gene expression patterns throughout the first in vitro spermatogenic wave.

RESULTS

Transcriptomic analyses reveal that cultured tissues expression profiles are almost comparable between D16 and D30; highlighting an abnormal kinetic throughout the second half of the first spermatogenesis during in vitro cultures. In addition, testicular explants have shown dysregulation of their transcriptomic profile compared to controls with genes related to inflammation response, insulin-like growth factor and genes involved in steroidogenesis.

DISCUSSION

The present work first shows that cryopreservation had very little impact on gene expression in testicular tissue, either directly after thawing or after 30 days in culture. Transcriptomic analysis of testis tissue samples is highly informative due to the large number of expressed genes and identified isoforms. This study provides a very valuable basis for future studies concerning in vitro spermatogenesis in mice.

The anti-androgenic fungicide triticonazole induces region-specific transcriptional changes in the developing rat perineum and phallus

Intrauterine exposure to endocrine disrupting chemicals can interfere with male reproductive development. This can lead to male reproductive disorders such as hypospadias, cryptorchidism and reduced fertility, as well as shorter anogenital distance (AGD) - a biomarker for incomplete androgen-dependent fetal masculinization. However, it remains challenging to predict adverse in vivo outcomes based on in vitro effect patterns for many chemicals. This is a challenge for modern toxicology, which aims to reduce animal testing for chemical safety assessments. To enable the transition towards higher reliance on alternative test methods, we need to better map underlying mechanisms leading to adverse effects. Herein, we have analyzed the transcriptome of the perineum and phallus of male fetal rats and defined the impacts of exposure to an anti-androgenic fungicide, triticonazole. Previously we have shown that developmental exposure to triticonazole can induce short male AGD, but without a marked effect on the transcriptome of the fetal testes. In contrast, we report here significant changes to the transcriptional landscape of the perineum and phallus, including regional differences between these adjacent tissues. This highlights the importance of analyzing the correct tissue when characterizing mechanisms of complex in vivo effect outcomes. Our results provide a rich resource for the spatiotemporal gene networks that are involved in the development of male external genitalia, and that can be disrupted upon exposure to chemicals that prevent normal masculinization of the perineum and phallus. Such data will be critical in the development of novel alternative test methods to determine the endocrine disrupting potential of existing and emerging chemicals.

Hypoxia and ERα Transcriptional Crosstalk Is Associated with Endocrine Resistance in Breast Cancer

Simple Summary

Hormone receptor positive breast cancer patients are treated with anti-hormone molecules as a standard of care. However, resistance frequently occurs, leading to hormone resistant metastatic relapses in foreign organs. Understanding the molecular mechanisms through which breast cancer cells evade therapeutic pressure is of paramount interest. Hypoxia, which refers to oxygen deprivation and is characterized by the activation of hypoxia inducible factors, is a common feature of the solid tumor microenvironment, yet its influence on estrogen receptor alpha activity remains elusive. Here, we investigate the consequence of hypoxia and the signaling of hypoxia inducible factors on hormone responsiveness in breast cancer cells and its clinical implications.

Abstract

Estrogen receptor-alpha (ERα) is the driving transcription factor in 70% of breast cancers and its activity is associated with hormone dependent tumor cell proliferation and survival. Given the recurrence of hormone resistant relapses, understanding the etiological factors fueling resistance is of major clinical interest. Hypoxia, a frequent feature of the solid tumor microenvironment, has been described to promote endocrine resistance by triggering ERα down-regulation in both in vitro and in vivo models. Yet, the consequences of hypoxia on ERα genomic activity remain largely elusive. In the present study, transcriptomic analysis shows that hypoxia regulates a fraction of ERα target genes, underlying an important regulatory overlap between hypoxic and estrogenic signaling. This gene expression reprogramming is associated with a massive reorganization of ERα cistrome, highlighted by a massive loss of ERα binding sites. Profiling of enhancer acetylation revealed a hormone independent enhancer activation at the vicinity of genes harboring hypoxia inducible factor (HIFα) binding sites, the major transcription factors governing hypoxic adaptation. This activation counterbalances the loss of ERα and sustains hormone-independent gene expression. We describe hypoxia in luminal ERα (+) breast cancer as a key factor interfering with endocrine therapies, associated with poor clinical prognosis in breast cancer patients.

PFOS disrupts key developmental pathways during hiPSC-derived cardiomyocyte differentiation in vitro

Exposure to perfluorooctanesulfonic acid (PFOS) has been associated with congenital heart disease (CHD) and decreased birth weight. PFOS exposure can disrupt signaling pathways relevant for cardiac development in stem cell-derived cardiomyocyte assays, such as the PluriBeat assay, where spheroids of human induced pluripotent stem cells (hiPSCs) differentiate into contracting cardiomyocytes. Notably, cell line origin can also affect how the assay responds to chemical exposure. Herein, we examined the effect of PFOS on cardiomyocyte differentiation by transcriptomics profiling of two different hiPSC lines to see if they exhibit a common pattern of disruption. Two stages of differentiation were investigated: the cardiac progenitor stage and the cardiomyocyte stage. Many differentially expressed genes (DEGs) were observed between cell lines independent of exposure. However, 135 DEGs were identified as common between the two cell lines. Of these, 10 DEGs were associated with GO-terms related to the heart. PFOS exposure disrupted multiple signaling pathways relevant to cardiac development, including WNT, TGF, HH, and EGF. Of these pathways, genes related to the non-canonical WNTCa²⁺ signaling was particularly affected. PFOS thus has the capacity to disrupt pathways important for cardiac development and function.

PFOS-induced thyroid hormone system disrupted rats display organ-specific changes in their transcriptomes

Perfluorooctanesulfonic acid (PFOS) is a persistent anthropogenic chemical that can affect the thyroid hormone system in humans and animals. In adults, thyroid hormones (THs) are regulated by the hypothalamic-pituitary-thyroid (HPT) axis, but also by organs such as the liver and potentially the gut microbiota. PFOS and other xenobiotics can therefore disrupt the TH system at various locations and through different mechanisms. To start addressing this, we exposed adult male rats to 3 mg PFOS/kg/day for 7 days and analysed effects on multiple organs and pathways simultaneously by transcriptomics. This included four primary organs involved in TH regulation, namely hypothalamus, pituitary, thyroid, and liver. To investigate a potential role of the gut microbiota in thyroid hormone regulation, two additional groups of animals were dosed with the antibiotic vancomycin (8 mg/kg/day), either with or without PFOS. PFOS exposure decreased thyroxine (T4) and triiodothyronine (T3) without affecting thyroid stimulating hormone (TSH), resembling a state of hypothyroxinemia. PFOS exposure resulted in 50 differentially expressed genes (DEGs) in the hypothalamus, 68 DEGs in the pituitary, 71 DEGs in the thyroid, and 181 DEGs in the liver. A concomitant compromised gut microbiota did not significantly change effects of PFOS exposure. Organ-specific DEGs did not align with TH regulating genes; however, genes associated with vesicle transport and neuronal signaling were affected in the hypothalamus, and phase I and phase II metabolism in the liver. This suggests that a decrease in systemic TH levels may activate the expression of factors altering trafficking, metabolism and excretion of TH. At the transcriptional level, little evidence suggests that the pituitary or thyroid gland is involved in PFOS-induced TH system disruption.

Understanding the Underlying Molecular Mechanisms of Meiotic Arrest during In Vitro Spermatogenesis in Rat Prepubertal Testicular Tissue

In vitro spermatogenesis appears to be a promising approach to restore the fertility of childhood cancer survivors. The rat model has proven to be challenging, since germ cell maturation is arrested in organotypic cultures. Here, we report that, despite a meiotic entry, abnormal synaptonemal complexes were found in spermatocytes, and in vitro matured rat prepubertal testicular tissues displayed an immature phenotype. RNA-sequencing analyses highlighted up to 600 differentially expressed genes between in vitro and in vivo conditions, including genes involved in blood-testis barrier (BTB) formation and steroidogenesis. BTB integrity, the expression of two steroidogenic enzymes, and androgen receptors were indeed altered in vitro. Moreover, most of the top 10 predicted upstream regulators of deregulated genes were involved in inflammatory processes or immune cell recruitment. However, none of the three anti-inflammatory molecules tested in this study promoted meiotic progression. By analysing for the first time in vitro matured rat prepubertal testicular tissues at the molecular level, we uncovered the deregulation of several genes and revealed that defective BTB function, altered steroidogenic pathway, and probably inflammation, could be at the origin of meiotic arrest.

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.

Multiple metastatic clones assessed by an integrative multiomics strategy in clear cell renal carcinoma: a case study

The dynamics of metastatic evolution in clear cell renal cell carcinoma (ccRCC) are complex. We report a case study where tumour heterogeneity resulting from clonal evolution is a frequent feature and could play a role in metastatic dissemination.We used an integrative multiomics strategy combining genomic and transcriptomic data to classify fourteen specimens from spatially different areas of a kidney tumour and three non-primary sites including a vein thrombus and two adrenal metastases.All sites were heterogeneous and polyclonal, each tumour site containing two different aggressive subclonal populations, with differentially expressed genes implicated in distinct biological functions. These are rare primary metastatic samples prior to any medical treatment, where we showed a multiple metastatic seeding of two subclonal populations.Multiple interdependent lineages could be the source of metastatic heterogeneity in ccRCC. By sampling metastases, patients with resistance to therapies could benefit a combination of targeted therapies based on more than one aggressive clone.

Transcriptome analysis of fetal rat testis following intrauterine exposure to the azole fungicides triticonazole and flusilazole reveals subtle changes despite adverse endocrine effects

Azoles are used in agriculture and medicine to combat fungal infections. We have previously examined the endocrine disrupting properties of the agricultural azole fungicides triticonazole and flusilazole. Triticonazole displayed strong androgen receptor (AR) antagonism in vitro, whereas in utero exposure resulted in anti-androgenic effects in vivo evidenced by shorter anogenital distance (AGD) in fetal male rats. Flusilazole displayed strong AR antagonism, but less potent than triticonazole, and disrupted steroidogenesis in vitro, whereas in utero exposure disrupted fetal male plasma hormone levels. To elaborate on how these azole fungicides can disrupt male reproductive development by different mechanisms, and to investigate whether feminization effects such as short AGD in males can also be detected at the transcript level in fetal testes, we profiled fetal testis transcriptomes after in utero exposure to triticonazole and flusilazole by 3'Digital Gene Expression (3'DGE). The analysis revealed few transcriptional changes after exposure to either compound at gestation day 17 and 21. This suggests that the observed influence of flusilazole on hormone production may be by directly targeting steroidogenic enzyme activity in the testis at the protein level, whereas observations of shorter AGD by triticonazole may primarily be due to disturbed androgen signaling in androgen-sensitive tissues. Expression of Calb2 and Gsta2 was altered by flusilazole but not triticonazole and may pinpoint novel pathways of disrupted testicular steroid synthesis. Our findings have wider implication for how we integrate omics data in chemical testing frameworks, including selection of non-animal test methods and building of Adverse Outcome Pathways for regulatory purposes.

Combined RNA/tissue profiling identifies novel Cancer/testis genes

Cancer/Testis (CT) genes are induced in germ cells, repressed in somatic cells, and derepressed in somatic tumors, where these genes can contribute to cancer progression. CT gene identification requires data obtained using standardized protocols and technologies. This is a challenge because data for germ cells, gonads, normal somatic tissues, and a wide range of cancer samples stem from multiple sources and were generated over substantial periods of time. We carried out a GeneChip‐based RNA profiling analysis using our own data for testis and enriched germ cells, data for somatic cancers from the Expression Project for Oncology, and data for normal somatic tissues from the Gene Omnibus Repository. We identified 478 candidate loci that include known CT genes, numerous genes associated with oncogenic processes, and novel candidates that are not referenced in the Cancer/Testis Database (www.cta.lncc.br). We complemented RNA expression data at the protein level for SPESP1, GALNTL5, PDCL2, and C11orf42 using cancer tissue microarrays covering malignant tumors of breast, uterus, thyroid, and kidney, as well as published RNA profiling and immunohistochemical data provided by the Human Protein Atlas (www.proteinatlas.org). We report that combined RNA/tissue profiling identifies novel CT genes that may be of clinical interest as therapeutical targets or biomarkers. Our findings also highlight the challenges of detecting truly germ cell‐specific mRNAs and the proteins they encode in highly heterogenous testicular, somatic, and tumor tissues.

Potential for Virus Endogenization in Humans through Testicular Germ Cell Infection: the Case of HIV

Viruses have colonized the germ line of our ancestors on several occasions during evolution, leading to the integration in the human genome of viral sequences from over 30 retroviral groups and a few nonretroviruses. Among the recently emerged viruses infecting humans, several target the testis (e.g., human immunodeficiency virus [HIV], Zika virus, and Ebola virus). Here, we aimed to investigate whether human testicular germ cells (TGCs) can support integration by HIV, a contemporary retrovirus that started to spread in the human population during the last century. We report that albeit alternative receptors enabled HIV-1 binding to TGCs, HIV virions failed to infect TGCs in vitro Nevertheless, exposure of TGCs to infected lymphocytes, naturally present in the testis from HIV+ men, led to HIV-1 entry, integration, and early protein expression. Similarly, cell-associated infection or bypassing viral entry led to HIV-1 integration in a spermatogonial cell line. Using DNAscope, HIV-1 and simian immunodeficiency virus (SIV) DNA were detected within a few TGCs in the testis from one infected patient, one rhesus macaque, and one African green monkey in vivo Molecular landscape analysis revealed that early TGCs were enriched in HIV early cofactors up to integration and had overall low antiviral defenses compared with testicular macrophages and Sertoli cells. In conclusion, our study reveals that TGCs can support the entry and integration of HIV upon cell-associated infection. This could represent a way for this contemporary virus to integrate into our germ line and become endogenous in the future, as happened during human evolution for a number of viruses.IMPORTANCE Viruses have colonized the host germ line on many occasions during evolution to eventually become endogenous. Here, we aimed at investigating whether human testicular germ cells (TGCs) can support such viral invasion by studying HIV interactions with TGCs in vitro Our results indicate that isolated primary TGCs express alternative HIV-1 receptors, allowing virion binding but not entry. However, HIV-1 entered and integrated into TGCs upon cell-associated infection and produced low levels of viral proteins. In vivo, HIV-1 and SIV DNA was detected in a few TGCs. Molecular landscape analysis showed that TGCs have overall weak antiviral defenses. Altogether, our results indicate that human TGCs can support HIV-1 early replication, including integration, suggesting potential for endogenization in future generations.

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.

Distinct Transcriptional Profiles of the Female, Male, and Finasteride-Induced Feminized Male Anogenital Region in Rat Fetuses

A short anogenital distance (AGD) in males is a marker for incomplete masculinization and a predictor of adverse effects on male reproductive health. For this reason, AGD is used to assess the endocrine disrupting potential of chemicals for risk assessment purposes. The molecular mechanisms underpinning this chemically induced shortening of the AGD, however, remains unclear. Although it is clear that androgen receptor-mediated signaling is essential, evidence also suggest the involvement of other signaling pathways. This study presents the first global transcriptional profile of the anogenital tissue in male rat fetuses with chemically induced short AGD, also including comparison to normal male and female control animals. The antiandrogenic drug finasteride (10 mg/kg bw/day) was used to induce short AGD by exposing time-mated Sprague Dawley rats at gestation days 7-21. The AGD was 37% shorter in exposed male fetuses compared with control males at gestation day 21. Transcriptomics analysis on anogenital tissues revealed a sexually dimorphic transcriptional profile. More than 350 genes were found to be differentially expressed between the 3 groups. The expression pattern of 4 genes of particular interest (Esr1, Padi2, Wnt2, and Sfrp4) was also tested by RT-qPCR analyses, indicating that estrogen and Wnt2 signaling play a role in the sexually dimorphic development of the anogenital region. Our transcriptomics profiles provide a stepping-stone for future studies aimed at characterizing the molecular events governing development of the anogenital tissues, as well as describing the detailed Adverse Outcome Pathways for short AGD; an accepted biomarker of endocrine effects for chemical risk assessment.

Characterization of Glyceollins as Novel Aryl Hydrocarbon Receptor Ligands and Their Role in Cell Migration

Recent studies strongly support the use of the aryl hydrocarbon receptor (AhR) as a therapeutic target in breast cancer. Glyceollins, a group of soybean phytoalexins, are known to exert therapeutic effects in chronic human diseases and also in cancer. To investigate the interaction between glyceollin I (GI), glyceollin II (GII) and AhR, a computational docking analysis, luciferase assays, immunofluorescence and transcriptome analyses were performed with different cancer cell lines. The docking experiments predicted that GI and GII can enter into the AhR binding pocket, but their interactions with the amino acids of the binding site differ, in part, from those interacting with 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Both GI and GII were able to weakly and partially activate AhR, with GII being more potent. The results from the transcriptome assays showed that approximately 10% of the genes regulated by TCDD were also modified by both GI and GII, which could have either antagonistic or synergistic effects upon TCDD activation. In addition, we report here, on the basis of phenotype, that GI and GII inhibit the migration of triple-negative (ER-, PgR-, HER2NEU-) MDA-MB-231 breast cancer cells, and that they inhibit the expression of genes which code for important regulators of cell migration and invasion in cancer tissues. In conclusion, GI and GII are AhR ligands that should be further investigated to determine their usefulness in cancer treatments.

Deciphering the Molecular Mechanisms Sustaining the Estrogenic Activity of the Two Major Dietary Compounds Zearalenone and Apigenin in ER-Positive Breast Cancer Cell Lines

The flavone apigenin and the mycotoxin zearalenone are two major compounds found in the human diet which bind estrogen receptors (ERs), and therefore influence ER activity. However, the underlying mechanisms are not well known. To unravel the molecular mechanisms that could explain the differential effect of zearalenone and apigenin on ER-positive breast cancer cell proliferation, gene-reporter assays, chromatin immunoprecipitation (ChIP) experiments, proliferation assays and transcriptomic analysis were performed. We found that zearalenone and apigenin transactivated ERs and promoted the expression of estradiol (E2)-responsive genes. However, zearalenone clearly enhanced cellular proliferation, while apigenin appeared to be antiestrogenic in the presence of E2 in both ER-positive breast cancer cell lines, MCF-7 and T47D. The transcriptomic analysis showed that both compounds regulate gene expression in the same way, but with differences in intensity. Two major sets of genes were identified; one set was linked to the cell cycle and the other set was linked to stress response and growth arrest. Our results show that the transcription dynamics in gene regulation induced by apigenin were somehow different with zearalenone and E2 and may explain the differential effect of these compounds on the phenotype of the breast cancer cell. Together, our results confirmed the potential health benefit effect of apigenin, while zearalenone appeared to be a true endocrine-disrupting compound.

Genome-Wide Transcriptional and Functional Analysis of Human T Lymphocytes Treated with Benzo[α]pyrene

Polycyclic aromatic hydrocarbons (PAHs) are widely distributed environmental contaminants, known to affect T lymphocytes. However, the molecular targets and pathways involved in their immunotoxic effects in human T lymphocytes remain unknown. Here, we analyzed the gene expression profile of primary human T lymphocytes treated with the prototypical PAH, benzo[α]pyrene (B[α]P), using a microarray-based transcriptome analysis. After a 48 h exposure to B[α]P, we identified 158 genes differentially expressed in T lymphocytes, including not only genes well-known to be affected by PAHs such as the cytochromes P450 (CYP) 1A1 and 1B1, but also others not previously shown to be targeted by B[α]P such as genes encoding the gap junction beta (GJB)-2 and 6 proteins. Functional enrichment analysis revealed that these candidates were significantly associated with the aryl hydrocarbon (AhR) and interferon (IFN) signaling pathways; a marked alteration in T lymphocyte recruitment was also observed. Using functional tests in transwell migration experiments, B[α]P was then shown to significantly decrease the chemokine (C-X-C motif) ligand 12-induced chemotaxis and transendothelial migration of T lymphocytes. In total, this study opens the way to unsuspected responsive pathway of interest, i.e., T lymphocyte migration, thus providing a more thorough understanding of the molecular basis of the immunotoxicity of PAHs.

funcExplorer: a tool for fast data-driven functional characterisation of high-throughput expression data

BACKGROUND

A widely applied approach to extract knowledge from high-throughput genomic data is clustering of gene expression profiles followed by functional enrichment analysis. This type of analysis, when done manually, is highly subjective and has limited reproducibility. Moreover, this pipeline can be very time-consuming and resource-demanding as enrichment analysis is done for tens to hundreds of clusters at a time. Thus, the task often needs programming skills to form a pipeline of different software tools or R packages to enable an automated approach. Furthermore, visualising the results can be challenging.

RESULTS

We developed a web tool, funcExplorer, which automatically combines hierarchical clustering and enrichment analysis to detect functionally related gene clusters. The functional characterisation is achieved using structured knowledge from data sources such as Gene Ontology, KEGG and Reactome pathways, Human Protein Atlas, and Human Phenotype Ontology. funcExplorer includes various measures for finding biologically meaningful clusters, provides a modern graphical user interface, and has wide-ranging data export and sharing options as well as software transparency by open-source code. The results are presented in a visually compact and interactive format, enabling users to explore the biological essence of the data. We compared our results with previously published gene clusters to demonstrate that funcExplorer can perform the data characterisation equally well, but without requiring labour-intensive manual interference.

CONCLUSIONS

The open-source web tool funcExplorer enables scientists with high-throughput genomic data to obtain a preliminary interactive overview of the expression patterns, gene names, and shared functionalities in their dataset in a visually pleasing format. funcExplorer is publicly available at https://biit.cs.ut.ee/funcexplorer.

Quantitative proteome profiling of dystrophic dog skeletal muscle reveals a stabilized muscular architecture and protection against oxidative stress after systemic delivery of MuStem cells

Proteomic profiling plays a decisive role in the elucidation of molecular signatures representative of a specific clinical context. MuStem cell based therapy represents a promising approach for clinical applications to cure Duchenne muscular dystrophy (DMD). To expand our previous studies collected in the clinically relevant DMD animal model, we decided to investigate the skeletal muscle proteome 4 months after systemic delivery of allogenic MuStem cells. Quantitative proteomics with isotope-coded protein labeling was used to compile quantitative changes in the protein expression profiles of muscle in transplanted Golden Retriever muscular dystrophy (GRMD) dogs as compared to Golden Retriever muscular dystrophy dogs. A total of 492 proteins were quantified, including 25 that were overrepresented and 46 that were underrepresented after MuStem cell transplantation. Interestingly, this study demonstrates that somatic stem cell therapy impacts on the structural integrity of the muscle fascicle by acting on fibers and its connections with the extracellular matrix. We also show that cell infusion promotes protective mechanisms against oxidative stress and favors the initial phase of muscle repair. This study allows us to identify putative candidates for tissue markers that might be of great value in objectively exploring the clinical benefits resulting from our cell-based therapy for DMD. All MS data have been deposited in the ProteomeXchange with identifier PXD001768 (http://proteomecentral.proteomexchange.org/dataset/PXD001768).

Glyceollins trigger anti-proliferative effects through estradiol-dependent and independent pathways in breast cancer cells

Background

Estrogen receptors (ER) α and β are found in both women and men in many tissues, where they have different functions, including having roles in cell proliferation and differentiation of the reproductive tract. In addition to estradiol (E2), a natural hormone, numerous compounds are able to bind ERs and modulate their activities. Among these compounds, phytoestrogens such as isoflavones, which are found in plants, are promising therapeutics for several pathologies. Glyceollins are second metabolites of isoflavones that are mainly produced in soybean in response to an elicitor. They have potentially therapeutic actions in breast cancer by reducing the proliferation of cancer cells. However, the molecular mechanisms driving these effects remain elusive.

Methods

First, to determine the proliferative or anti-proliferative effects of glyceollins, in vivo and in vitro approaches were used. The length of epithelial duct in mammary gland as well as uterotrophy after treatment by E2 and glyceollins and their effect on proliferation of different breast cell line were assessed. Secondly, the ability of glyceollin to activate ER was assessed by luciferase assay. Finally, to unravel molecular mechanisms involved by glyceollins, transcriptomic analysis was performed on MCF-7 breast cancer cells.

Results

In this study, we show that synthetic versions of glyceollin I and II exert anti-proliferative effects in vivo in mouse mammary glands and in vitro in different ER-positive and ER-negative breast cell lines. Using transcriptomic analysis, we produce for the first time an integrated view of gene regulation in response to glyceollins and reveal that these phytochemicals act through at least two major pathways. One pathway involving FOXM1 and ERα is directly linked to proliferation. The other involves the HIF family and reveals that stress is a potential factor in the anti-proliferative effects of glyceollins due to its role in increasing the expression of REDD1, an mTORC1 inhibitor.

Conclusion

Overall, our study clearly shows that glyceollins exert anti-proliferative effects by reducing the expression of genes encoding cell cycle and mitosis-associated factors and biomarkers overexpressed in cancers and by increasing the expression of growth arrest-related genes. These results reinforce the therapeutic potential of glyceollins for breast cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s12964-017-0182-1) contains supplementary material, which is available to authorized users.

Endogenous IL-33 Deficiency Exacerbates Liver Injury and Increases Hepatic Influx of Neutrophils in Acute Murine Viral Hepatitis

The alarmin IL-33 has been described to be upregulated in human and murine viral hepatitis. However, the role of endogenous IL-33 in viral hepatitis remains obscure. We aimed to decipher its function by infecting IL-33-deficient mice (IL-33 KO) and their wild-type (WT) littermates with pathogenic mouse hepatitis virus (L2-MHV3). The IL-33 KO mice were more sensitive to L2-MHV3 infection exhibiting higher levels of AST/ALT, higher tissue damage, significant weight loss, and earlier death. An increased depletion of B and T lymphocytes, NKT cells, dendritic cells, and macrophages was observed 48 h postinfection (PI) in IL-33 KO mice than that in WT mice. In contrast, a massive influx of neutrophils was observed in IL-33 KO mice at 48 h PI. A transcriptomic study of inflammatory and cell-signaling genes revealed the overexpression of IL-6, TNFα, and several chemokines involved in recruitment/activation of neutrophils (CXCL2, CXCL5, CCL2, and CCL6) at 72 h PI in IL-33 KO mice. However, the IFNγ was strongly induced in WT mice with less profound expression in IL-33 KO mice demonstrating that endogenous IL-33 regulated IFNγ expression during L2-MHV3 hepatitis. In conclusion, we demonstrated that endogenous IL-33 had multifaceted immunoregulatory effect during viral hepatitis via induction of IFNγ, survival effect on immune cells, and infiltration of neutrophils in the liver.

Cold/menthol TRPM8 receptors initiate the cold-shock response and protect germ cells from cold-shock-induced oxidation

Testes of most male mammals present the particularity of being externalized from the body and are consequently slightly cooler than core body temperature (4-8°C below). Although, hypothermia of the testis is known to increase germ cells apoptosis, little is known about the underlying molecular mechanisms, including cold sensors, transduction pathways, and apoptosis triggers. In this study, using a functional knockout mouse model of the cold and menthol receptors, dubbed transient receptor potential melastatine 8 (TRPM8) channels, we found that TRPM8 initiated the cold-shock response by differentially modulating cold- and heat-shock proteins. Besides, apoptosis of germ cells increased in proportion to the cooling level in control mice but was independent of temperature in knockout mice. We also observed that the rate of germ cell death correlated positively with the reactive oxygen species level and negatively with the expression of the detoxifying enzymes. This result suggests that the TRPM8 sensor is a key determinant of germ cell fate under hypothermic stimulation.-Borowiec, A.-S., Sion, B., Chalmel, F., Rolland, A. D., Lemonnier, L., De Clerck, T., Bokhobza, A., Derouiche, S., Dewailly, E., Slomianny, C., Mauduit, C., Benahmed, M., Roudbaraki, M., Jégou, B., Prevarskaya, N., Bidaux, G. Cold/menthol TRPM8 receptors initiate the cold-shock response and protect germ cells from cold-shock-induced oxidation.

The epigenetic processes of meiosis in male mice are broadly affected by the widely used herbicide atrazine

BACKGROUND

Environmental factors such as pesticides can cause phenotypic changes in various organisms, including mammals. We studied the effects of the widely used herbicide atrazine (ATZ) on meiosis, a key step of gametogenesis, in male mice.

METHODS

Gene expression pattern was analysed by Gene-Chip array. Genome-wide mapping of H3K4me3 marks distribution was done by ChIP-sequencing of testis tissue using Illumina technologies. RT-qPCR was used to validate differentially expressed genes or differential peaks.

RESULTS

We demonstrate that exposure to ATZ reduces testosterone levels and the number of spermatozoa in the epididymis and delays meiosis. Using Gene-Chip and ChIP-Seq analysis of H3K4me3 marks, we found that a broad range of cellular functions, including GTPase activity, mitochondrial function and steroid-hormone metabolism, are affected by ATZ. Furthermore, treated mice display enriched histone H3K4me3 marks in regions of strong recombination (double-strand break sites), within very large genes and reduced marks in the pseudoautosomal region of X chromosome.

CONCLUSIONS

Our data demonstrate that atrazine exposure interferes with normal meiosis, which affects spermatozoa production.

Differential Gene Expression Profiling of Dystrophic Dog Muscle after MuStem Cell Transplantation

BACKGROUND

Several adult stem cell populations exhibit myogenic regenerative potential, thus representing attractive candidates for therapeutic approaches of neuromuscular diseases such as Duchenne Muscular Dystrophy (DMD). We have recently shown that systemic delivery of MuStem cells, skeletal muscle-resident stem cells isolated in healthy dog, generates the remodelling of muscle tissue and gives rise to striking clinical benefits in Golden Retriever Muscular Dystrophy (GRMD) dog. This global effect, which is observed in the clinically relevant DMD animal model, leads us to question here the molecular pathways that are impacted by MuStem cell transplantation. To address this issue, we compare the global gene expression profile between healthy, GRMD and MuStem cell treated GRMD dog muscle, four months after allogenic MuStem cell transplantation.

RESULTS

In the dystrophic context of the GRMD dog, disease-related deregulation is observed in the case of 282 genes related to various processes such as inflammatory response, regeneration, calcium ion binding, extracellular matrix organization, metabolism and apoptosis regulation. Importantly, we reveal the impact of MuStem cell transplantation on several molecular and cellular pathways based on a selection of 31 genes displaying signals specifically modulated by the treatment. Concomitant with a diffuse dystrophin expression, a histological remodelling and a stabilization of GRMD dog clinical status, we show that cell delivery is associated with an up-regulation of genes reflecting a sustained enhancement of muscle regeneration. We also identify a decreased mRNA expression of a set of genes having metabolic functions associated with lipid homeostasis and energy. Interestingly, ubiquitin-mediated protein degradation is highly enhanced in GRMD dog muscle after systemic delivery of MuStem cells.

CONCLUSIONS

Overall, our results provide the first high-throughput characterization of GRMD dog muscle and throw new light on the complex molecular/cellular effects associated with muscle repair and the clinical efficacy of MuStem cell-based therapy.

Research resource: the dynamic transcriptional profile of sertoli cells during the progression of spermatogenesis

Sertoli cells (SCs), the only somatic cells within seminiferous tubules, associate intimately with developing germ cells. They not only provide physical and nutritional support but also secrete factors essential to the complex developmental processes of germ cell proliferation and differentiation. The SC transcriptome must therefore adapt rapidly during the different stages of spermatogenesis. We report comprehensive genome-wide expression profiles of pure populations of SCs isolated at 5 distinct stages of the first wave of mouse spermatogenesis, using RNA sequencing technology. We were able to reconstruct about 13 901 high-confidence, nonredundant coding and noncoding transcripts, characterized by complex alternative splicing patterns with more than 45% comprising novel isoforms of known genes. Interestingly, roughly one-fifth (2939) of these genes exhibited a dynamic expression profile reflecting the evolving role of SCs during the progression of spermatogenesis, with stage-specific expression of genes involved in biological processes such as cell cycle regulation, metabolism and energy production, retinoic acid synthesis, and blood-testis barrier biogenesis. Finally, regulatory network analysis identified the transcription factors endothelial PAS domain-containing protein 1 (EPAS1/Hif2α), aryl hydrocarbon receptor nuclear translocator (ARNT/Hif1β), and signal transducer and activator of transcription 1 (STAT1) as potential master regulators driving the SC transcriptional program. Our results highlight the plastic transcriptional landscape of SCs during the progression of spermatogenesis and provide valuable resources to better understand SC function and spermatogenesis and its related disorders, such as male infertility.

Combining RNA and protein profiling data with network interactions identifies genes associated with spermatogenesis in mouse and human

Genome-wide RNA profiling studies have identified hundreds of transcripts that are highly expressed in mammalian male germ cells, including many that are undetectable in somatic control tissues. Among them, genes important for spermatogenesis are significantly enriched. Information about mRNAs and their cognate proteins facilitates the identification of novel conserved target genes for functional studies in the mouse. By inspecting genome-wide RNA profiling data, we manually selected 81 genes for which RNA is detected almost exclusively in the human male germline and, in most cases, in rodent testicular germ cells. We observed corresponding mRNA/protein patterns in 43 cases using immunohistochemical data from the Human Protein Atlas and large-scale human protein profiling data obtained via mass spectroscopy. Protein network information enabled us to establish an interaction map of 38 proteins that points to potentially important testicular roles for some of them. We further characterized six candidate genes at the protein level in the mouse. We conclude that conserved genes induced in testis tend to show similar mRNA/protein expression patterns across species. Specifically, our results suggest roles during embryogenesis and adult spermatogenesis for Foxr1 and Sox30 and during spermiogenesis and fertility for Fam71b, 1700019N19Rik, Hmgb4, and Zfp597.

Forty-four novel protein-coding loci discovered using a proteomics informed by transcriptomics (PIT) approach in rat male germ cells

Spermatogenesis is a complex process, dependent upon the successive activation and/or repression of thousands of gene products, and ends with the production of haploid male gametes. RNA sequencing of male germ cells in the rat identified thousands of novel testicular unannotated transcripts (TUTs). Although such RNAs are usually annotated as long noncoding RNAs (lncRNAs), it is possible that some of these TUTs code for protein. To test this possibility, we used a "proteomics informed by transcriptomics" (PIT) strategy combining RNA sequencing data with shotgun proteomics analyses of spermatocytes and spermatids in the rat. Among 3559 TUTs and 506 lncRNAs found in meiotic and postmeiotic germ cells, 44 encoded at least one peptide. We showed that these novel high-confidence protein-coding loci exhibit several genomic features intermediate between those of lncRNAs and mRNAs. We experimentally validated the testicular expression pattern of two of these novel protein-coding gene candidates, both highly conserved in mammals: one for a vesicle-associated membrane protein we named VAMP-9, and the other for an enolase domain-containing protein. This study confirms the potential of PIT approaches for the discovery of protein-coding transcripts initially thought to be untranslated or unknown transcripts. Our results contribute to the understanding of spermatogenesis by characterizing two novel proteins, implicated by their strong expression in germ cells. The mass spectrometry proteomics data have been deposited with the ProteomeXchange Consortium under the data set identifier PXD000872.

An integrative omics strategy to assess the germ cell secretome and to decipher sertoli-germ cell crosstalk in the Mammalian testis

Mammalian spermatogenesis, which takes place in complex testicular structures called seminiferous tubules, is a highly specialized process controlled by the integration of juxtacrine, paracrine and endocrine information. Within the seminiferous tubules, the germ cells and Sertoli cells are surrounded by testicular fluid (TF), which probably contains most of the secreted proteins involved in crosstalk between these cells. It has already been established that germ cells can modulate somatic Sertoli cell function through the secretion of diffusible factors. We studied the germ cell secretome, which was previously considered inaccessible, by analyzing the TF collected by microsurgery in an “integrative omics” strategy combining proteomics, transcriptomics, genomics and interactomics data. This approach identified a set of proteins preferentially secreted by Sertoli cells or germ cells. An interaction network analysis revealed complex, interlaced cell-cell dialog between the secretome and membranome of seminiferous cells, mediated via the TF. We then focused on germ cell-secreted candidate proteins, and we identified several potential interacting partners located on the surface of Sertoli cells. Two interactions, APOH/CDC42 and APP/NGFR, were validated in situ, in a proximity ligation assay (PLA). Our results provide new insight into the crosstalk between germ cells and Sertoli cells occurring during spermatogenesis. Our findings also demonstrate that this “integrative omics” strategy is powerful enough for data mining and highlighting meaningful cell-cell communication events between different types of cells in a complex tissue, via a biological fluid. This integrative strategy could be applied more widely, to gain access to secretomes that have proved difficult to study whilst avoiding the limitations of in vitro culture.

Proteomics and integrative genomics for unraveling the mysteries of spermatogenesis: the strategies of a team

The strikingly complex structural organization of the mammalian testis in vivo creates particular difficulties for studies of its organization, function and regulation. These difficulties are particularly pronounced for investigations of the molecular communication networks within the seminiferous tubules that govern spermatogenesis. The use of classical molecular and cell biology approaches to unravel this complexity has proved problematic, due to difficulties in maintaining differentiated germ cells in vitro, in particular. The lack of a suitable testing ground has led to a greater reliance on high-quality proteomic and genomic analyses as a prelude to the in vitro antx1d in vivo testing of hypotheses. In this study, we highlight the options currently available for research, as used in our laboratory, in which proteomic and integrative genomic strategies are applied to the study of spermatogenesis in mammals. We will comment on results providing insight into the molecular mechanisms underlying normal and pathological spermatogenesis and new perspectives for the treatment of male infertility in humans. Finally, we will discuss the relevance of our strategies and the unexpected potential and perspectives they offer to teams involved in the study of male reproduction, within the framework of the Human Proteome Project.

SIGNIFICANCE

Integrative genomics is becoming a powerful strategy for discovering the biological significance hidden in proteomic datasets. This work introduces some of the integrative genomic concepts and works used by our team to gain new insight into mammalian spermatogenesis, a remarkably sophisticated process. We demonstrate the relevance of these integrative approaches to understand the cellular cross talks established between the somatic Sertoli cells and the germ cell lineage, within the seminiferous epithelium. Our work also contributes to new knowledge on the pathophysiology of testicular function, with promising clinical applications. This article is part of a Special Issue entitled: 20years of Proteomics in memory of Viatliano Pallini. Guest Editors: Luca Bini, Juan J. Calvete, Natacha Turck, Denis Hochstrasser and Jean-Charles Sanchez.

High-resolution profiling of novel transcribed regions during rat spermatogenesis

Mammalian spermatogenesis is a complex and highly orchestrated combination of processes in which male germline proliferation and differentiation result in the production of mature spermatozoa. If recent genome-wide studies have contributed to the in-depth analysis of the male germline protein-encoding transcriptome, little effort has yet been devoted to the systematic identification of novel unannotated transcribed regions expressed during mammalian spermatogenesis. We report high-resolution expression profiling of male germ cells in rat, using next-generation sequencing technology and highly enriched testicular cell populations. Among 20 424 high-confidence transcripts reconstructed, we defined a stringent set of 1419 long multi-exonic unannotated transcripts expressed in the testis (testis-expressed unannotated transcripts [TUTs]). TUTs were divided into 7 groups with different expression patterns. Most TUTs share many of the characteristics of vertebrate long noncoding RNAs (lncRNAs). We also markedly reinforced the finding that TUTs and known lncRNAs accumulate during the meiotic and postmeiotic stages of spermatogenesis in mammals and that X-linked meiotic TUTs do not escape the silencing effects of meiotic sex chromosome inactivation. Importantly, we discovered that TUTs and known lncRNAs with a peak expression during meiosis define a distinct class of noncoding transcripts that exhibit exons twice as long as those of other transcripts. Our study provides new insights in transcriptional profiling of the male germline and represents a high-quality resource for novel loci expressed during spermatogenesis that significantly contributes to rat genome annotation.

An atlas of chromatoid body components

The genome of male germ cells is actively transcribed during spermatogenesis to produce phase-specific protein-coding mRNAs and a considerable amount of different noncoding RNAs. Ribonucleoprotein (RNP) granule-mediated RNA regulation provides a powerful means to secure the quality and correct expression of the requisite transcripts. Haploid spermatids are characterized by a unique, unusually large cytoplasmic granule, the chromatoid body (CB), which emerges during the switch between the meiotic and post-meiotic phases of spermatogenesis. To better understand the role of the CB in male germ cell differentiation, we isolated CBs from mouse testes and revealed its full RNA and protein composition. We showed that the CB is mainly composed of RNA-binding proteins and other proteins involved RNA regulation. The CB was loaded with RNA, including pachytene piRNAs, a diverse set of mRNAs, and a number of uncharacterized long noncoding transcripts. The CB was demonstrated to accumulate nascent RNA during all the steps of round spermatid differentiation. Our results revealed the CB as a large germ cell-specific RNP platform that is involved in the control of the highly complex transcriptome of haploid male germ cells.

Expression of HLA-G by mast cells is associated with hepatitis C virus-induced liver fibrosis

BACKGROUND & AIMS

Infection with hepatitis C virus is a worldwide health problem. An inadequate Th2 cytokine response promotes the fibrosis-cirrhosis fate. Immune-modulating molecules favoring a Th2 profile, such as HLA-G molecules of the HLA class Ib family, may play a role in chronic hepatitis. HLA-G contributes to the escape of tumors, and their involvement in viral infections has been increasingly described. The aim of this work was to study the expression of HLA-G in the liver, its cellular source and its regulation in cases of chronic C hepatitis.

METHODS

HLA-G cells in blocks of liver derived from patients infected with HCV were labeled by immunohistochemistry and enumerated. Double immunofluorescence allowed the identification of the cellular source. HLA-G secretion by a human mast cell line was quantified by ELISA after various stimulations. After treatment with IFN-α, real-time PCR was performed to determine the kinetics of cytokine expression profiles, followed by heat map clustering analysis.

RESULTS

The number of HLA-G+ cells was significantly associated with the area of fibrosis. For the first time, we identify the HLA-G+ cells as being mast cells. HLA-G secretion was significantly induced in human mast cells stimulated by IL-10 or interferons of class I. The transcriptome of the secretome of this cell line stimulated by IFN-α revealed that (i) the HLA-G gene is upregulated late, and that (ii) T lymphocytes and NK cells are recruited.

CONCLUSIONS

These findings suggest an autocrine loop in the genesis of HCV liver fibrosis, based on mast cells expressing HLA-G.

Systemic compensatory response to neonatal estradiol exposure does not prevent depletion of the oocyte pool in the rat

The formation of ovarian follicles is a finely tuned process that takes place within a narrow time-window in rodents. Multiple factors and pathways have been proposed to contribute to the mechanisms triggering this process but the role of endocrine factors, especially estrogens, remains elusive. It is currently hypothesized that removal from the maternal hormonal environment permits follicle formation at birth. However, experimentally-induced maintenance of high 17β-estradiol (E2) levels leads to subtle, distinct, immediate effects on follicle formation and oocyte survival depending on the species and dose. In this study, we examined the immediate effects of neonatal E2 exposure from post-natal day (PND) 0 to PND2 on the whole organism and on ovarian follicle formation in rats. Measurements of plasma E2, estrone and their sulfate conjugates after E2 exposure showed that neonatal female rats rapidly acquire the capability to metabolize and clear excessive E2 levels. Concomitant modifications to the mRNA content of genes encoding selected E2 metabolism enzymes in the liver and the ovary in response to E2 exposure indicate that E2 may modify the neonatal maturation of these organs. In the liver, E2 treatment was associated with lower acquisition of the capability to metabolize E2. In the ovary, E2 depleted the oocyte pool in a dose dependent manner by PND3. In 10 µg/day E2-treated ovaries, apoptotic oocytes were observed in newly formed follicles in addition to areas of ovarian cord remodeling. At PND6, follicles without any visible oocyte were present and multi-oocyte follicles were not observed. Our study reveals a major species-difference. Indeed, neonatal exposure to E2 depletes the oocyte pool in the rat ovary, whereas in the mouse it is well known to increase oocyte survival.

Fsh controls gene expression in fish both independently of and through steroid mediation

The mechanisms and the mediators relaying Fsh action on testicular functions are poorly understood. Unlike in mammals, in fish both gonadotropins (Fsh and Lh) are able to efficiently stimulate steroidogenesis, likely through a direct interaction with their cognate receptors present on the Leydig cells. In this context, it is crucial to understand if Fsh effects are mediated through the production of steroids. To address this issue we performed transcriptome studies after in vitro incubations of rainbow trout testis explants in the presence of Fsh alone or in combination with trilostane, an inhibitor of Δ4- steroidogenesis. Trilostane significantly reduced or suppressed the response of many genes to Fsh (like wisp1, testis gapdhs, cldn11, inha, vt1 or dmrt1) showing that, in fish, important aspects of Fsh action follow indirect pathways and require the production of Δ4-steroids. What is more, most of the genes regulated by Fsh through steroid mediation were similarly regulated by Lh (and/or androgens). In contrast, the response to Fsh of other genes was not suppressed in the presence of trilostane. These latter included genes encoding for anti-mullerian hormone, midkine a (pleiotrophin related), angiopoietine-related protein, cyclins E1 and G1, hepatocyte growth factor activator, insulin-like growth factor 1b/3. A majority of those genes were preferentially regulated by Fsh, when compared to Lh, suggesting that specific regulatory effects of Fsh did not depend on steroid production. Finally, antagonistic effects between Fsh and steroids were found, in particular for genes encoding key factors of steroidogenesis (star, hsd3b1, cyp11b2-2) or for genes of the Igf system (igf1b/3). Our study provides the first clear evidence that, in fish, Fsh exerts Δ4-steroid-independent regulatory functions on many genes which are highly relevant for the onset of spermatogenesis.

Expression screening of cancer/testis genes in prostate cancer identifies NR6A1 as a novel marker of disease progression and aggressiveness

BACKGROUND

Cancer/Testis (CT) genes are expressed in male gonads, repressed in most healthy somatic tissues and de-repressed in various somatic malignancies including prostate cancers (PCa). Because of their specific expression signature and their associations with tumor aggressiveness and poor outcomes, CT genes are considered to be useful biomarkers and they are also targets for the development of new anti-cancer immunotherapies. The aim of this study was to identify novel CT genes associated with hormone-sensitive prostate cancer (HSPC), and castration-resistant prostate cancer (CRPC).

METHODS

To identify novel CT genes we screened genes for which transcripts were detected by RNA profiling specifically in normal testis and in either HSPC or CRPC as compared to normal prostate and 44 other healthy tissues using GeneChips. The expression and clinicopathological significance of a promising candidate--NR6A1--was examined in HSPC, CRPC, and metastatic site samples using tissue microarrays.

RESULTS

We report the identification of 98 genes detected in CRPC, HSPC and testicular samples but not in the normal controls. Among them, cellular levels of NR6A1 were found to be higher in HSPC compared to normal prostate and further increased in metastatic lesions and CRPC. Furthermore, increased NR6A1 immunoreactivity was significantly associated with a high Gleason score, advanced pT stage and cancer cell proliferation.

CONCLUSIONS

Our results show that cellular levels of NR6A1 are correlated with disease progression in PCa. We suggest that this essential orphan nuclear receptor is a potential therapeutic target as well as a biomarker of PCa aggressiveness.

FSH and LH have common and distinct effects on gene expression in rainbow trout testis

The general rules established from mammalian species for the regulation of spermatogenesis by gonadotropins may not be fully relevant in fish. Particularly, Fsh is as potent as Lh to stimulate steroidogenesis and the Fsh receptor is expressed in Leydig cells. In seasonal breeders, Fsh is likely the major gonadotropin involved in spermatogenesis onset and Lh is required to support spermatogenesis progression and gamete release. However, the genes that relay the action of Fsh and Lh have been poorly investigated in fish. The present study was aimed at identifying gonadotropin-dependent genes expressed in the testis during fish puberty. We cultured pubertal trout testicular explants for 96 h, with or without gonadotropin, and analyzed transcriptome variations using microarrays. Fsh and Lh had similar effects on a large group of genes while other genes were preferentially regulated by one or the other gonadotropin. We showed that most of the responsive genes were expressed in somatic cells and exhibited relevant patterns during the seasonal reproductive cycle. Some genes preferentially modulated by Lh could be involved in testicular cell fate (pvrl1 and bty) or sperm maturation (ehmt2 and racgap1) and will deserve further examination. Besides Fsh's effects on the steroidogenic pathway, our study demonstrates that Fsh coordinates relevant stimulatory and inhibitory paracrine factors known to regulate early germ cell proliferation and differentiation. Some of these genes belong to major regulatory pathways including the Igf pathway (igf1b/igf3 and igfbp6), the Tgfb pathway (amh, inha, inhba, and fstl3), the Wnt pathway (wisp1), and pleiotrophin (mdka).

Genome-wide identification of Sox8-, and Sox9-dependent genes during early post-natal testis development in the mouse

The SOX8 and SOX9 transcription factors are involved in, among others, sex differentiation, male gonad development and adult maintenance of spermatogenesis. Sox8(-/-) mice lacking Sox9 in Sertoli cells fail to form testis cords and cannot establish spermatogenesis. Although genetic and histological data show an important role for these transcription factors in regulating spermatogenesis, it is not clear which genes depend upon them at a genome-wide level. To identify transcripts that respond to the absence of Sox8 in all cells and Sox9 in Sertoli cells we measured mRNA concentrations in testicular samples from mice at 0, 6 and 18 days post-partum. In total, 621 and 629 transcripts were found at decreased or increased levels, respectively, at different time points in the mutant as compared to the control samples. These mRNAs were categorized as preferentially expressed in Sertoli cells or germ cells using data obtained with male and female gonad samples and enriched testicular cell populations. Five candidate genes were validated at the protein level. Furthermore, we identified putative direct SOX8 and SOX9 target genes by integrating predicted SOX-binding sites present in potential regulatory regions upstream of the transcription start site. Finally, we used protein network data to gain insight into the effects on regulatory interactions that occur when Sox8 and Sox9 are absent in developing Sertoli cells. The integration of testicular samples with enriched Sertoli cells, germ cells and female gonads enabled us to broadly distinguish transcripts directly affected in Sertoli cells from others that respond to secondary events in testicular cell types. Thus, combined RNA profiling signals, motif predictions and network data identified putative SOX8/SOX9 target genes in Sertoli cells and yielded insight into regulatory interactions that depend upon these transcription factors. In addition, our results will facilitate the interpretation of genome-wide in vivo SOX8 and SOX9 DNA binding data.

Profiling of androgen response in rainbow trout pubertal testis: relevance to male gonad development and spermatogenesis

The capacity of testicular somatic cells to promote and sustain germ cell differentiation is largely regulated by sexual steroids and notably androgens. In fish species the importance of androgens is emphasized by their ability to induce sex reversal of the developing fries and to trigger spermatogenesis. Here we studied the influence of androgens on testicular gene expression in trout testis using microarrays. Following treatment of immature males with physiological doses of testosterone or 11-ketotestosterone, 418 genes that exhibit changes in expression were identified. Interestingly, the activity of testosterone appeared stronger than that of 11-ketotestosterone. Expression profiles of responsive genes throughout testis development and in isolated germ cells confirmed androgens to mainly affect gene expression in somatic cells. Furthermore, specific clusters of genes that exhibit regulation coincidently with changes in the natural circulating levels of androgens during the reproductive cycle were highlighted, reinforcing the physiological significance of these data. Among somatic genes, a phylogenetic footprinting study identified putative androgen response elements within the proximal promoter regions of 42 potential direct androgen target genes. Finally, androgens were also found to alter the germ line towards meiotic expression profiles, supporting the hypothesis of a role for the somatic responsive genes in driving germ cell fate. This study significantly increases our understanding of molecular pathways regulated by androgens in vertebrates. The highly cyclic testicular development in trout together with functions associated with regulated genes reveal potential mechanisms for androgen actions in tubule formation, steroid production, germ cell development and sperm secretion.

Insulin and IGF1 receptors are essential for XX and XY gonadal differentiation and adrenal development in mice

Mouse sex determination provides an attractive model to study how regulatory genetic networks and signaling pathways control cell specification and cell fate decisions. This study characterizes in detail the essential role played by the insulin receptor (INSR) and the IGF type I receptor (IGF1R) in adrenogenital development and primary sex determination. Constitutive ablation of insulin/IGF signaling pathway led to reduced proliferation rate of somatic progenitor cells in both XX and XY gonads prior to sex determination together with the downregulation of hundreds of genes associated with the adrenal, testicular, and ovarian genetic programs. These findings indicate that prior to sex determination somatic progenitors in Insr;Igf1r mutant gonads are not lineage primed and thus incapable of upregulating/repressing the male and female genetic programs required for cell fate restriction. In consequence, embryos lacking functional insulin/IGF signaling exhibit (i) complete agenesis of the adrenal cortex, (ii) embryonic XY gonadal sex reversal, with a delay of Sry upregulation and the subsequent failure of the testicular genetic program, and (iii) a delay in ovarian differentiation so that Insr;Igf1r mutant gonads, irrespective of genetic sex, remained in an extended undifferentiated state, before the ovarian differentiation program ultimately is initiated at around E16.5.

Spermatogenesis in mammals: proteomic insights

Spermatogenesis is a highly sophisticated process involved in the transmission of genetic heritage. It includes halving ploidy, repackaging of the chromatin for transport, and the equipment of developing spermatids and eventually spermatozoa with the advanced apparatus (e.g., tightly packed mitochondrial sheat in the mid piece, elongating of the tail, reduction of cytoplasmic volume) to elicit motility once they reach the epididymis. Mammalian spermatogenesis is divided into three phases. In the first the primitive germ cells or spermatogonia undergo a series of mitotic divisions. In the second the spermatocytes undergo two consecutive divisions in meiosis to produce haploid spermatids. In the third the spermatids differentiate into spermatozoa in a process called spermiogenesis. Paracrine, autocrine, juxtacrine, and endocrine pathways all contribute to the regulation of the process. The array of structural elements and chemical factors modulating somatic and germ cell activity is such that the network linking the various cellular activities during spermatogenesis is unimaginably complex. Over the past two decades, advances in genomics have greatly improved our knowledge of spermatogenesis, by identifying numerous genes essential for the development of functional male gametes. Large-scale analyses of testicular function have deepened our insight into normal and pathological spermatogenesis. Progress in genome sequencing and microarray technology have been exploited for genome-wide expression studies, leading to the identification of hundreds of genes differentially expressed within the testis. However, although proteomics has now come of age, the proteomics-based investigation of spermatogenesis remains in its infancy. Here, we review the state-of-the-art of large-scale proteomic analyses of spermatogenesis, from germ cell development during sex determination to spermatogenesis in the adult. Indeed, a few laboratories have undertaken differential protein profiling expression studies and/or systematic analyses of testicular proteomes in entire organs or isolated cells from various species. We consider the pros and cons of proteomics for studying the testicular germ cell gene expression program. Finally, we address the use of protein datasets, through integrative genomics (i.e., combining genomics, transcriptomics, and proteomics), bioinformatics, and modelling.