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.

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.

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.

A tridimensional atlas of the developing human head

The evolution and development of the head have long captivated researchers due to the crucial role of the head as the gateway for sensory stimuli and the intricate structural complexity of the head. Although significant progress has been made in understanding head development in various vertebrate species, our knowledge of early human head ontogeny remains limited. Here, we used advanced whole-mount immunostaining and 3D imaging techniques to generate a comprehensive 3D cellular atlas of human head embryogenesis. We present detailed developmental series of diverse head tissues and cell types, including muscles, vasculature, cartilage, peripheral nerves, and exocrine glands. These datasets, accessible through a dedicated web interface, provide insights into human embryogenesis. We offer perspectives on the branching morphogenesis of human exocrine glands and unknown features of the development of neurovascular and skeletomuscular structures. These insights into human embryology have important implications for understanding craniofacial defects and neurological disorders and advancing diagnostic and therapeutic strategies.

Loss of NR5A1 in mouse Sertoli cells after sex determination changes cellular identity and induces cell death by anoikis

To investigate the role of the nuclear receptor NR5A1 in the testis after sex determination, we analyzed mice lacking NR5A1 in Sertoli cells (SCs) from embryonic day (E) 13.5 onwards. Ablation of Nr5a1 impaired the expression of genes characteristic of SC identity (e.g. Sox9 and Amh), caused SC death from E14.5 onwards through a Trp53-independent mechanism related to anoikis, and induced disorganization of the testis cords. Together, these effects caused germ cells to enter meiosis and die. Single-cell RNA-sequencing experiments revealed that NR5A1-deficient SCs changed their molecular identity: some acquired a 'pre-granulosa-like' cell identity, whereas other reverted to a 'supporting progenitor-like' cell identity, most of them being 'intersex' because they expressed both testicular and ovarian genes. Fetal Leydig cells (LCs) did not display significant changes, indicating that SCs are not required beyond E14.5 for their emergence or maintenance. In contrast, adult LCs were absent from postnatal testes. In addition, adult mutant males displayed persistence of Müllerian duct derivatives, decreased anogenital distance and reduced penis length, which could be explained by the loss of AMH and testosterone synthesis due to SC failure.

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.

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.

A new scaffold-free tumoroid model provides a robust preclinical tool to investigate invasion and drug response in Renal Cell Carcinoma

Clear cell Renal Cell Carcinoma (ccRCC) is one of the most prevalent kidney cancers, which is often asymptomatic and thus discovered at a metastatic state (mRCC). mRCC are highly heterogeneous tumors composed of subclonal populations that lead to poor treatment response rate. Several recent works explored the potential of ccRCC tumoroids culture derived from patients. However, these models were produced following a scaffold-based method using collagen I or Matrigel that exhibit lot variability and whose complexity could induce treatment response modifications and phenotypic alterations. Following the observation that ccRCC tumoroids can create their own niche by secreting extracellular matrix components, we developed the first scaffold-free tumoroid model of ccRCC tumors. Tumoroids from mice as well as from human tumors were generated with high success rate (≥90%) using a magnetic suspension method and standard culture media. Immunofluorescence analysis revealed their self-organization capacities to maintain multiple tumor-resident cell types, including endothelial progenitor cells. Transcriptomic analysis showed the reproducibility of the method highlighting that the majority of gene expression patterns was conserved in tumoroids compared to their matching tumor tissue. Moreover, this model enables to evaluate drug effects and invasiveness of renal cancer cells in a 3D context, providing a robust preclinical tool for drug screening and biomarker assessment in line with alternative ex vivo methods like tumor tissue slice culture or in vivo xenograft models.

Uterine defects and estradiol-dependent development of oviductal diverticula in mice lacking the SMAD4 C-terminal Mad homology 2 domain

In female mammals, the oviduct and uterus are essential sites for female and male gamete transport, fertilization, implantation, and maintenance of a successful pregnancy. To delineate the reproductive function of Mothers against decapentaplegic homolog 4 (Smad4), we specifically inactivated Smad4 in ovarian granulosa cells and, oviduct and uterine mesenchymal cells using the Amhr2-cre mouse line. Deletion of exon 8 of Smad4 results in the production of an MH2-truncated SMAD4 protein. These mutant mice are infertile due to the development of oviductal diverticula and defects during the implantation process. The ovaries are fully functional as demonstrated in an ovary transfer experiment. The development of oviductal diverticula occurs shortly after puberty and is dependent on estradiol. The diverticula interfere with sperm migration and embryo transit to the uterus, reducing the number of implantation sites. Analysis of the uterus shows that, even if implantation occurs, decidualization and vascularization are defective resulting in embryo resorption as early as the seventh day of pregnancy. Thus, Smad4 plays an important function in female reproduction by controlling the structural and functional integrity of the oviduct and uterus.

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.

New approach methods to improve human health risk assessment of thyroid hormone system disruption-a PARC project

Current test strategies to identify thyroid hormone (TH) system disruptors are inadequate for conducting robust chemical risk assessment required for regulation. The tests rely heavily on histopathological changes in rodent thyroid glands or measuring changes in systemic TH levels, but they lack specific new approach methodologies (NAMs) that can adequately detect TH-mediated effects. Such alternative test methods are needed to infer a causal relationship between molecular initiating events and adverse outcomes such as perturbed brain development. Although some NAMs that are relevant for TH system disruption are available-and are currently in the process of regulatory validation-there is still a need to develop more extensive alternative test batteries to cover the range of potential key events along the causal pathway between initial chemical disruption and adverse outcomes in humans. This project, funded under the Partnership for the Assessment of Risk from Chemicals (PARC) initiative, aims to facilitate the development of NAMs that are specific for TH system disruption by characterizing in vivo mechanisms of action that can be targeted by in embryo/in vitro/in silico/in chemico testing strategies. We will develop and improve human-relevant in vitro test systems to capture effects on important areas of the TH system. Furthermore, we will elaborate on important species differences in TH system disruption by incorporating non-mammalian vertebrate test species alongside classical laboratory rat species and human-derived in vitro assays.

Single-cell Deconvolution of a Specific Malignant Cell Population as a Poor Prognostic Biomarker in Low-risk Clear Cell Renal Cell Carcinoma Patients

BACKGROUND

Intratumor heterogeneity (ITH) is a key feature in clear cell renal cell carcinomas (ccRCCs) that impacts outcomes such as aggressiveness, response to treatments, or recurrence. In particular, it may explain tumor relapse after surgery in clinically low-risk patients who did not benefit from adjuvant therapy. Recently, single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool to unravel expression ITH (eITH) and might enable better assessment of clinical outcomes in ccRCC.

OBJECTIVE

To explore eITH in ccRCC with a focus on malignant cells (MCs) and assess its relevance to improve prognosis for low-risk patients.

DESIGN, SETTING, AND PARTICIPANTS

We performed scRNA-seq on tumor samples from five untreated ccRCC patients ranging from pT1a to pT3b. Data were complemented with a published dataset composed of pairs of matched normal and ccRCC samples.

INTERVENTION

Radical or partial nephrectomy on untreated ccRCC patients.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Viability and cell type proportions were determined by flow cytometry. Following scRNA-seq, a functional analysis was performed and tumor progression trajectories were inferred. A deconvolution approach was applied on an external cohort, and Kaplan-Meier survival curves were estimated with respect to the prevalence of malignant clusters.

RESULTS AND LIMITATIONS

We analyzed 54 812 cells and identified 35 cell subpopulations. The eITH analysis revealed that each tumor contained various degrees of clonal diversity. The transcriptomic signatures of MCs in one particularly heterogeneous sample were used to design a deconvolution-based strategy that allowed the risk stratification of 310 low-risk ccRCC patients.

CONCLUSIONS

We described eITH in ccRCCs, and used this information to establish significant cell population-based prognostic signatures and better discriminate ccRCC patients. This approach has the potential to improve the stratification of clinically low-risk patients and their therapeutic management.

PATIENT SUMMARY

We sequenced the RNA content of individual cell subpopulations composed of clear cell renal cell carcinomas and identified specific malignant cells the genetic information of which can be used to predict tumor progression.

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.

SMG6 localizes to the chromatoid body and shapes the male germ cell transcriptome to drive spermatogenesis

Nonsense-mediated RNA decay (NMD) is a highly conserved and selective RNA turnover pathway that depends on the endonuclease SMG6. Here, we show that SMG6 is essential for male germ cell differentiation in mice. Germ-cell conditional knockout (cKO) of Smg6 induces extensive transcriptome misregulation, including a failure to eliminate meiotically expressed transcripts in early haploid cells, and accumulation of NMD target mRNAs with long 3' untranslated regions (UTRs). Loss of SMG6 in the male germline results in complete arrest of spermatogenesis at the early haploid cell stage. We find that SMG6 is strikingly enriched in the chromatoid body (CB), a specialized cytoplasmic granule in male germ cells also harboring PIWI-interacting RNAs (piRNAs) and the piRNA-binding protein PIWIL1. This raises the possibility that SMG6 and the piRNA pathway function together, which is supported by several findings, including that Piwil1-KO mice phenocopy Smg6-cKO mice and that SMG6 and PIWIL1 co-regulate many genes in round spermatids. Together, our results demonstrate that SMG6 is an essential regulator of the male germline transcriptome, and highlight the CB as a molecular platform coordinating RNA regulatory pathways to control sperm production and fertility.

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.

Androglobin, a chimeric mammalian globin, is required for male fertility

Spermatogenesis is a highly specialized differentiation process driven by a dynamic gene expression program and ending with the production of mature spermatozoa. Whereas hundreds of genes are known to be essential for male germline proliferation and differentiation, the contribution of several genes remains uncharacterized. The predominant expression of the latest globin family member, androglobin (Adgb), in mammalian testis tissue prompted us to assess its physiological function in spermatogenesis. Adgb knockout mice display male infertility, reduced testis weight, impaired maturation of elongating spermatids, abnormal sperm shape, and ultrastructural defects in microtubule and mitochondrial organization. Epididymal sperm from Adgb knockout animals display multiple flagellar malformations including coiled, bifid or shortened flagella, and erratic acrosomal development. Following immunoprecipitation and mass spectrometry, we could identify septin 10 (Sept10) as interactor of Adgb. The Sept10-Adgb interaction was confirmed both in vivo using testis lysates and in vitro by reciprocal co-immunoprecipitation experiments. Furthermore, the absence of Adgb leads to mislocalization of Sept10 in sperm, indicating defective manchette and sperm annulus formation. Finally, in vitro data suggest that Adgb contributes to Sept10 proteolysis in a calmodulin-dependent manner. Collectively, our results provide evidence that Adgb is essential for murine spermatogenesis and further suggest that Adgb is required for sperm head shaping via the manchette and proper flagellum formation.

Transcriptomic analysis in zebrafish larvae identifies iron-dependent mitochondrial dysfunction as a possible key event of NAFLD progression induced by benzo[a]pyrene/ethanol co-exposure

Non-alcoholic fatty liver disease (NAFLD) is a worldwide epidemic for which environmental contaminants are increasingly recognized as important etiological factors. Among them, the combination of benzo[a]pyrene (B[a]P), a potent environmental carcinogen, with ethanol, was shown to induce the transition of steatosis toward steatohepatitis. However, the underlying mechanisms involved remain to be deciphered. In this context, we used high-fat diet fed zebrafish model, in which we previously observed progression of steatosis to a steatohepatitis-like state following a 7-day-co-exposure to 43 mM ethanol and 25 nM B[a]P. Transcriptomic analysis highlighted the potent role of mitochondrial dysfunction, alterations in heme and iron homeostasis, involvement of aryl hydrocarbon receptor (AhR) signaling, and oxidative stress. Most of these mRNA dysregulations were validated by RT-qPCR. Moreover, similar changes were observed using a human in vitro hepatocyte model, HepaRG cells. The mitochondria structural and functional alterations were confirmed by transmission electronic microscopy and Seahorse technology, respectively. Involvement of AhR signaling was evidenced by using in vivo an AhR antagonist, CH223191, and in vitro in AhR-knock-out HepaRG cells. Furthermore, as co-exposure was found to increase the levels of both heme and hemin, we investigated if mitochondrial iron could induce oxidative stress. We found that mitochondrial labile iron content was raised in toxicant-exposed larvae. This increase was prevented by the iron chelator, deferoxamine, which also inhibited liver co-exposure toxicity. Overall, these results suggest that the increase in mitochondrial iron content induced by B[a]P/ethanol co-exposure causes mitochondrial dysfunction that contributes to the pathological progression of NAFLD.

FlexDotPlot: a universal and modular dot plot visualization tool for complex multifaceted data

Motivation

Dot plots are heatmap-like charts that provide a compact way to simultaneously display two quantitative information by means of dots of different sizes and colors. Despite the popularity of this visualization method, particularly in single-cell RNA-sequencing (scRNA-seq) studies, existing tools used to make dot plots are limited in terms of functionality and usability.

Results

We developed FlexDotPlot, an R package for generating dot plots from multifaceted data, including scRNA-seq data. It provides a universal and easy-to-use solution with a high versatility. An interactive R Shiny application is also available allowing non-R users to easily generate dot plots with several tunable parameters.

Availability and implementation

Source code and detailed manual are available on CRAN (stable version) and at https://github.com/Simon-Leonard/FlexDotPlot (development version). Code to reproduce figures is available at https://github.com/Simon-Leonard/FlexDotPlot_paper. A Shiny app is available as a stand-alone application within the package.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

P–055 Methylation dynamics of the sperm epigenome after chemotherapy: a case study

Study question

What is the evolution of the sperm epigenome after chemotherapy in a patient with testicular cancer (TC)?

Summary answer

These new data on epigenetic recovery profil after TC are useful tools for counseling and reassuring these patients.

What is known already

An important issue for young men affected TC is how TC and its treatment will affect, transiently or permanently, their future reproductive health. The consequences of cancer treatment on the sperm epigenome during the recovery periods are topical issues of ascendant significance as epigenetic modifications to the paternal genome may have deleterious effects on the offspring.

Study design, size, duration

Here we report the epigenomic profiling of frozen sperm from a TC patient before and after the treatment at different time points (6, 9, 12 and 24 months) by using RRBS analysis (Reduced representation bisulfite sequencing method).

Participants/materials, setting, methods

A testicular tumor (testicular germ cell tumor) was diagnosed in a 30 years old patient. A cryopreservation of spermatozoa was proposed before treatment.Semen samples were obtained 2 times before treatment and 4 times after treatment (6, 9, 12 and 24 months following the initiation of treatment).

Main results and the role of chance

Upon collection, sampling after chemotherapy ranged from 0,6 to 4,2 million per sperm straw between 6 and 24 months after the treatment, always increasing.

In order to capture the direct effect of the treatment on the methylation changes, the DMR detection has been operated between pre-chemotherapy samples (pair-wise) and the time point of 6 months. Among the 179 hqDMRs, 74 are differentially methylated between the PreCT and PostCT6m samples (16 hyper- and 68 hypo-methylated) associated with 49 DMGs (15 hyper- and 34 hypo-methylated).

We further sub-clustered the 74 hqDMRs between PreCT and PostCT6m into 6 patterns, 3 hyper- and 3 hypo-methylated. Briefly, patterns P1 and P4 include hqDMRs that quickly get back to their pre-treatment methylation status just after 9th months onwards. Patterns P2 and P5 include hqDMRs that slowly get back to their pre-treatment methylation status between 12 and 24 months after treatment. Patterns P3 and P6 include hqDMRs that remain hyper- or hypo-methylated even after 24 months.

We have intersected the genes (DMGs) associated with the detected hqDMRs with those known to be important or expressed during embryogenesis. We thus detected that 7 hyper-methylated and 6 hypomethylated DMGs were involved (or expressed) during embryonic / fetal development.

Limitations, reasons for caution

This study involves a single patient. As the patient made no major changes in his personal way of life, we hypothesized that sperm parameter variations may be attributable to the BEP treatment.

Wider implications of the findings: The altered methylated status of those DMGs important for early development might modify their expression pattern and thus affect their function during key stages of embryogenesis leading to potential developmental disorders. It is important to notice that among the 110 DMGs none of them correspond to known imprinted genes.

Trial registration number

Not applicable

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.

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.

The ReproGenomics Viewer: a multi-omics and cross-species resource compatible with single-cell studies for the reproductive science community

MOTIVATION

Recent advances in transcriptomics have enabled unprecedented insight into gene expression analysis at a single-cell resolution. While it is anticipated that the number of publications based on such technologies will increase in the next decade, there is currently no public resource to centralize and enable scientists to explore single-cell datasets published in the field of reproductive biology.

RESULTS

Here, we present a major update of the ReproGenomics Viewer, a cross-species and cross-technology web-based resource of manually-curated sequencing datasets related to reproduction. The redesign of the ReproGenomics Viewer's architecture is accompanied by significant growth of the database content including several landmark single-cell RNA-sequencing datasets. The implementation of additional tools enables users to visualize and browse the complex, high-dimensional data now being generated in the reproductive field.

AVAILABILITY AND IMPLEMENTATION

The ReproGenomics Viewer resource is freely accessible at http://rgv.genouest.org. The website is implemented in Python, JavaScript and MongoDB, and is compatible with all major browsers. Source codes can be downloaded from https://github.com/fchalmel/RGV.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Interplay between transcription regulators RUNX1 and FUBP1 activates an enhancer of the oncogene c-KIT and amplifies cell proliferation

Runt-related transcription factor 1 (RUNX1) is a well-known master regulator of hematopoietic lineages but its mechanisms of action are still not fully understood. Here, we found that RUNX1 localizes on active chromatin together with Far Upstream Binding Protein 1 (FUBP1) in human B-cell precursor lymphoblasts, and that both factors interact in the same transcriptional regulatory complex. RUNX1 and FUBP1 chromatin localization identified c-KIT as a common target gene. We characterized two regulatory regions, at +700 bp and +30 kb within the first intron of c-KIT, bound by both RUNX1 and FUBP1, and that present active histone marks. Based on these regions, we proposed a novel FUBP1 FUSE-like DNA-binding sequence on the +30 kb enhancer. We demonstrated that FUBP1 and RUNX1 cooperate for the regulation of the expression of the oncogene c-KIT. Notably, upregulation of c-KIT expression by FUBP1 and RUNX1 promotes cell proliferation and renders cells more resistant to the c-KIT inhibitor imatinib mesylate, a common therapeutic drug. These results reveal a new mechanism of action of RUNX1 that implicates FUBP1, as a facilitator, to trigger transcriptional regulation of c-KIT and to regulate cell proliferation. Deregulation of this regulatory mechanism may explain some oncogenic function of RUNX1 and FUBP1.

New transcriptomic tools to understand testis development and functions

The testis plays a central role in the male reproductive system - secreting several hormones including male steroids and producing male gametes. A complex and coordinated molecular program is required for the proper differentiation of testicular cell types and maintenance of their functions in adulthood. The testicular transcriptome displays the highest levels of complexity and specificity across all tissues in a wide range of species. Many studies have used high-throughput sequencing technologies to define the molecular dynamics and regulatory networks in the testis as well as to identify novel genes or gene isoforms expressed in this organ. This review intends to highlight the complementarity of these transcriptomic studies and to show how the use of different sequencing protocols contribute to improve our global understanding of testicular biology.

Meiotic Genes Are Enriched in Regions of Reduced Archaic Ancestry

About 1-6% of the genetic ancestry of modern humans today originates from admixture with archaic humans. It has recently been shown that autosomal genomic regions with a reduced proportion of Neanderthal and Denisovan ancestries (NA and DA) are significantly enriched in genes that are more expressed in testis than in other tissues. To determine whether a cellular segregation pattern would exist, we combined maps of archaic introgression with a cross-analysis of three transcriptomic datasets deciphering the transcriptional landscape of human gonadal cell types. We reveal that the regions deficient in both NA and DA contain a significant enrichment of genes transcribed in meiotic germ cells. The interbreeding of anatomically modern humans with archaic humans may have introduced archaic-derived alleles that contributed to genetic incompatibilities affecting meiosis that were subsequently purged by natural selection.

Defining the human sperm microtubulome: an integrated genomics approach

Sperm motility notably depends on the structural integrity of the flagellum and the regulation of microtubule dynamics. Although researchers have started to use "omics" techniques to characterize the human sperm's molecular landscape, the constituents responsible for the assembly, organization, and dynamics of the flagellum microtubule have yet to be fully defined. In this study, we defined a core set of 116 gene products associated with the human sperm microtubulome (including products potentially involved in abnormal ciliary phenotypes and male infertility disorders). To this end, we designed and applied an integrated genomics workflow and combined relevant proteomics, transcriptomics, and interactomics datasets to reconstruct a dynamic interactome map. By further integrating phenotypic information, we defined a disease-interaction network; this enabled us to highlight a number of novel factors potentially associated with altered sperm motility and male fertility. Lastly, we experimentally validated the expression pattern of two candidate genes (CUL3 and DCDC2C) that had never previously been associated with male germline differentiation. Our analysis suggested that CUL3 and DCDC2C's products have important roles in the sperm flagellum. Taken as a whole, our results demonstrate that an integrated genomics strategy can highlight relevant molecular factors in specific sperm components. This approach could be easily extended by including other "omics" data (from asthenozoospermic men, for example) and identifying other critical proteins from the human sperm microtubulome.

Cytosine modifications modulate the chromatin architecture of transcriptional enhancers

Epigenetic mechanisms are believed to play key roles in the establishment of cell-specific transcription programs. Accordingly, the modified bases 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) have been observed in DNA of genomic regulatory regions such as enhancers, and oxidation of 5mC into 5hmC by Ten-eleven translocation (TET) proteins correlates with enhancer activation. However, the functional relationship between cytosine modifications and the chromatin architecture of enhancers remains elusive. To gain insights into their function, 5mC and 5hmC levels were perturbed by inhibiting DNA methyltransferases and TETs during differentiation of mouse embryonal carcinoma cells into neural progenitors, and chromatin characteristics of enhancers bound by the pioneer transcription factors FOXA1, MEIS1, and PBX1 were interrogated. In a large fraction of the tested enhancers, inhibition of DNA methylation was associated with a significant increase in monomethylation of H3K4, a characteristic mark of enhancer priming. In addition, at some specific enhancers, 5mC oxidation by TETs facilitated chromatin opening, a process that may stabilize MEIS1 binding to these genomic regions.

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.

Linking transcriptomics and proteomics in spermatogenesis

Spermatogenesis is a complex and tightly regulated process leading to the continuous production of male gametes, the spermatozoa. This developmental process requires the sequential and coordinated expression of thousands of genes, including many that are testis-specific. The molecular networks underlying normal and pathological spermatogenesis have been widely investigated in recent decades, and many high-throughput expression studies have studied genes and proteins involved in male fertility. In this review, we focus on studies that have attempted to correlate transcription and translation during spermatogenesis by comparing the testicular transcriptome and proteome. We also discuss the recent development and use of new transcriptomic approaches that provide a better proxy for the proteome, from both qualitative and quantitative perspectives. Finally, we provide illustrations of how testis-derived transcriptomic and proteomic data can be integrated to address new questions and how the 'proteomics informed by transcriptomics' technique, by combining RNA-seq and MS-based proteomics, can contribute significantly to the discovery of new protein-coding genes or new protein isoforms expressed during spermatogenesis.

PepPSy: a web server to prioritize gene products in experimental and biocuration workflows

Among the 20 000 human gene products predicted from genome annotation, about 3000 still lack validation at protein level. We developed PepPSy, a user-friendly gene expression-based prioritization system, to help investigators to determine in which human tissues they should look for an unseen protein. PepPSy can also be used by biocurators to revisit the annotation of specific categories of proteins based on the ‘omics’ data housed by the system. In this study, it was used to prioritize 21 dubious protein-coding genes among the 616 annotated in neXtProt for reannotation. PepPSy is freely available at http://peppsy.genouest.org.

Database URL:http://peppsy.genouest.org.

DICER Regulates the Formation and Maintenance of Cell-Cell Junctions in the Mouse Seminiferous Epithelium

The endonuclease DICER that processes micro-RNAs and small interfering RNAs is essential for normal spermatogenesis and male fertility. We previously showed that the deletion of Dicer1 gene in postnatal spermatogonia in mice using Ngn3 promoter-driven Cre expression caused severe defects in the morphogenesis of haploid spermatid to mature spermatozoon, including problems in cell polarization and nuclear elongation. In this study, we further analyzed the same mouse model and revealed that absence of functional DICER in differentiating male germ cells induces disorganization of the cell-cell junctions in the seminiferous epithelium. We detected discontinuous and irregular apical ectoplasmic specializations between elongating spermatids and Sertoli cells. The defective anchoring of spermatids to Sertoli cells caused a premature release of spermatids into the lumen. Our findings may help also explain the abnormal elongation process of remaining spermatids because these junctions and the correct positioning of germ cells in the epithelium are critically important for the progression of spermiogenesis. Interestingly, cell adhesion-related genes were generally upregulated in Dicer1 knockout germ cells. Claudin 5 ( Cldn5 ) was among the most upregulated genes and we show that the polarized localization of CLAUDIN5 in the apical ectoplasmic specializations was lost in Dicer1 knockout spermatids. Our results suggest that DICER-dependent pathways control the formation and organization of cell-cell junctions in the seminiferous epithelium via the regulation of cell adhesion-related genes.

The ReproGenomics Viewer: an integrative cross-species toolbox for the reproductive science community

We report the development of the ReproGenomics Viewer (RGV), a multi- and cross-species working environment for the visualization, mining and comparison of published omics data sets for the reproductive science community. The system currently embeds 15 published data sets related to gametogenesis from nine model organisms. Data sets have been curated and conveniently organized into broad categories including biological topics, technologies, species and publications. RGV's modular design for both organisms and genomic tools enables users to upload and compare their data with that from the data sets embedded in the system in a cross-species manner. The RGV is freely available at http://rgv.genouest.org.

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.