Transcription of two long noncoding RNAs mediates mating-type control of gametogenesis in budding yeast.

The cell-fate decision leading to gametogenesis is essential for sexual reproduction. In S. cerevisiae, only diploid MATa/α but not haploid MATa or MATα cells undergo gametogenesis, known as sporulation. We find that transcription of two long noncoding RNAs (lncRNAs) mediates mating-type control of sporulation. In MATa or MATα haploids, expression of IME1, the central inducer of gametogenesis, is inhibited in cis by transcription of the lncRNA IRT1, located in the IME1 promoter. IRT1 transcription recruits the Set2 histone methyltransferase and the Set3 histone deacetylase complex to establish repressive chromatin at the IME1 promoter. Inhibiting expression of IRT1 and an antisense transcript that antagonizes the expression of the meiotic regulator IME4 allows cells expressing the haploid mating type to sporulate with kinetics that are indistinguishable from that of MATa/α diploids. Conversely, expression of the two lncRNAs abolishes sporulation in MATa/α diploids. Thus, transcription of two lncRNAs governs mating-type control of gametogenesis in yeast.

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.

The conserved histone deacetylase Rpd3 and its DNA binding subunit Ume6 control dynamic transcript architecture during mitotic growth and meiotic development.

It was recently reported that the sizes of many mRNAs change when budding yeast cells exit mitosis and enter the meiotic differentiation pathway. These differences were attributed to length variations of their untranslated regions. The function of UTRs in protein translation is well established. However, the mechanism controlling the expression of distinct transcript isoforms during mitotic growth and meiotic development is unknown. In this study, we order developmentally regulated transcript isoforms according to their expression at specific stages during meiosis and gametogenesis, as compared to vegetative growth and starvation. We employ regulatory motif prediction, in vivo protein-DNA binding assays, genetic analyses and monitoring of epigenetic amino acid modification patterns to identify a novel role for Rpd3 and Ume6, two components of a histone deacetylase complex already known to repress early meiosis-specific genes in dividing cells, in mitotic repression of meiosis-specific transcript isoforms. Our findings classify developmental stage-specific early, middle and late meiotic transcript isoforms, and they point to a novel HDAC-dependent control mechanism for flexible transcript architecture during cell growth and differentiation. Since Rpd3 is highly conserved and ubiquitously expressed in many tissues, our results are likely relevant for development and disease in higher eukaryotes.

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).

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.

Global alterations of the transcriptional landscape during yeast growth and development in the absence of Ume6-dependent chromatin modification.

Chromatin modification enzymes are important regulators of gene expression and some are evolutionarily conserved from yeast to human. Saccharomyces cerevisiae is a major model organism for genome-wide studies that aim at the identification of target genes under the control of conserved epigenetic regulators. Ume6 interacts with the upstream repressor site 1 (URS1) and represses transcription by recruiting both the conserved histone deacetylase Rpd3 (through the co-repressor Sin3) and the chromatin-remodeling factor Isw2. Cells lacking Ume6 are defective in growth, stress response, and meiotic development. RNA profiling studies and in vivo protein-DNA binding assays identified mRNAs or transcript isoforms that are directly repressed by Ume6 in mitosis. However, a comprehensive understanding of the transcriptional alterations, which underlie the complex ume6Δ mutant phenotype during fermentation, respiration, or sporulation, is lacking. We report the protein-coding transcriptome of a diploid MAT a/α wild-type and ume6/ume6 mutant strains cultured in rich media with glucose or acetate as a carbon source, or sporulation-inducing medium. We distinguished direct from indirect effects on mRNA levels by combining GeneChip data with URS1 motif predictions and published high-throughput in vivo Ume6-DNA binding data. To gain insight into the molecular interactions between successive waves of Ume6-dependent meiotic genes, we integrated expression data with information on protein networks. Our work identifies novel Ume6 repressed genes during growth and development and reveals a strong effect of the carbon source on the derepression pattern of transcripts in growing and developmentally arrested ume6/ume6 mutant cells. Since yeast is a useful model organism for chromatin-mediated effects on gene expression, our results provide a rich source for further genetic and molecular biological work on the regulation of cell growth and cell differentiation in eukaryotes.

Developmental stage dependent metabolic regulation during meiotic differentiation in budding yeast.

BACKGROUND: The meiotic developmental pathway in yeast enables both differentiation of vegetative cells into haploid spores that ensure long-term survival, and recombination of the parental DNA to create genetic diversity. Despite the importance of proper metabolic regulation for the supply of building blocks and energy, little is known about the reprogramming of central metabolic pathways in meiotically differentiating cells during passage through successive developmental stages. RESULTS: Metabolic regulation during meiotic differentiation in budding yeast was analysed by integrating information on genome-wide transcriptional activity, 26 enzymatic activities in the central metabolism, the dynamics of 67 metabolites, and a metabolic flux analysis at mid-stage meiosis. Analyses of mutants arresting sporulation at defined stages demonstrated that metabolic reprogramming is tightly controlled by the progression through the developmental pathway. The correlation between transcript levels and enzymatic activities in the central metabolism varies significantly in a developmental-stage dependent manner. The complete loss of phosphofructokinase activity at mid-stage meiosis enables a unique setup of the glycolytic pathway which facilitates carbon flux repartitioning into synthesis of spore-wall precursors during the co-assimilation of glycogen and acetate. The need for correct homeostasis of purine nucleotides during the meiotic differentiation was demonstrated by the sporulation defect of the AMP deaminase mutant, amd1, which exhibited hyper-accumulation of ATP accompanied by depletion of guanosine nucleotides. CONCLUSIONS: Our systems-level analysis shows that reprogramming of the central metabolism during the meiotic differentiation is controlled at different hierarchical levels to meet the metabolic and energetic needs at successive developmental stages.

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.

Direct iterative protein profiling (DIPP) - an innovative method for large-scale protein detection applied to budding yeast mitosis.

The budding yeast Saccharomyces cerevisiae is a major model organism for important biological processes such as mitotic growth and meiotic development, it can be a human pathogen, and it is widely used in the food-, and biotechnology industries. Consequently, the genomes of numerous strains have been sequenced and a very large amount of RNA profiling data is available. Moreover, it has recently become possible to quantitatively analyze the entire yeast proteome; however, efficient and cost-effective high-throughput protein profiling remains a challenge. We report here a new approach to direct and label-free large-scale yeast protein identification using a tandem buffer system for protein extraction, two-step protein prefractionation and enzymatic digestion, and detection of peptides by iterative mass spectrometry. Our profiling study of diploid cells undergoing rapid mitotic growth identified 86% of the known proteins and its output was found to be widely concordant with genome-wide mRNA concentrations and DNA variations between yeast strains. This paves the way for comprehensive and straightforward yeast proteome profiling across a wide variety of experimental conditions.

Execution of the meiotic noncoding RNA expression program and the onset of gametogenesis in yeast require the conserved exosome subunit Rrp6.

Budding yeast noncoding RNAs (ncRNAs) are pervasively transcribed during mitosis, and some regulate mitotic protein-coding genes. However, little is known about ncRNA expression during meiotic development. Using high-resolution profiling we identified an extensive meiotic ncRNA expression program interlaced with the protein-coding transcriptome via sense/antisense transcript pairs, bidirectional promoters, and ncRNAs that overlap the regulatory regions of genes. Meiotic unannotated transcripts (MUTs) are mitotic targets of the conserved exosome component Rrp6, which itself is degraded after the onset of meiosis when MUTs and other ncRNAs accumulate in successive waves. Diploid cells lacking Rrp6 fail to initiate premeiotic DNA replication normally and cannot undergo efficient meiotic development. The present study demonstrates a unique function for budding yeast Rrp6 in degrading distinct classes of meiotically induced ncRNAs during vegetative growth and the onset of meiosis and thus points to a critical role of differential ncRNA expression in the execution of a conserved developmental program.

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.

Enhanced identification of endocrine disruptors through integration of science-based regulatory practices and innovative methodologies: The MERLON Project.

The prevalence of hormone-related health issues caused by exposure to endocrine disrupting chemicals (EDCs) is a significant, and increasing, societal challenge. Declining fertility rates together with rising incidence rates of reproductive disorders and other endocrine-related diseases underscores the urgency in taking more action. Addressing the growing threat of EDCs in our environment demands robust and reliable test methods to assess a broad variety of endpoints relevant for endocrine disruption. EDCs also require effective regulatory frameworks, especially as the current move towards greater reliance on non-animal methods in chemical testing puts to test the current paradigm for EDC identification, which requires that an adverse effect is observed in an intact organism. Although great advances have been made in the field of predictive toxicology, disruption to the endocrine system and subsequent adverse health effects may prove particularly difficult to predict without traditional animal models. The MERLON project seeks to expedite progress by integrating multispecies molecular research, new approach methodologies (NAMs), human clinical epidemiology, and systems biology to furnish mechanistic insights and explore ways forward for NAM-based identification of EDCs. The focus is on sexual development and function, from foetal sex differentiation of the reproductive system through mini-puberty and puberty to sexual maturity. The project aims are geared towards closing existing knowledge gaps in understanding the effects of EDCs on human health to ultimately support effective regulation of EDCs in the European Union and beyond.

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.

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).

XY Sox9 embryonic loss-of-function mouse mutants show complete sex reversal and produce partially fertile XY oocytes.

Gonadal differentiation is the first step of mammalian sex determination. The expression of the Y chromosomal testis determining factor Sry leads to up-regulation of the transcription factor Sox9 which promotes testis differentiation. Previous studies showed that Sox9 deficiency induces expression of ovarian markers in XY mutant fetal gonads before they die. To better understand the genome-wide transcriptional profile underlying this process we compared samples from XY Sf1:Cre(Tg/+); Sox9(flox/flox) mutant gonads in which Sox9 is ablated in Sertoli-precursor cells during early stages of gonad development to XX Sox9(flox/flox) ovaries and XY Sox9(flox/flox) testes at E13.5. We found a complex mRNA signature that indicates wide-spread transcriptional de-regulation and revealed for XY mutants at E13.5 an intermediate transcript profile between male and female gonads. However, XY Sf1:Cre(Tg/+); Sox9(flox/flox) mutant gonads develop as ovaries containing XY developing follicles at P0 but less frequently so than in XX control ovaries. Furthermore, we studied the extent to which developing XY mutant ovaries are able to mediate adult fertility and observed that XY oocytes from XY mutant ovaries are competent for fertilization; however, two thirds of them fail to develop beyond two-cell stage embryos. Taken together, we found that XY Sf1:Cre(Tg/+); Sox9(flox/flox) females are capable of producing viable offspring albeit at a reduced level.

The homeobox gene CHX10/VSX2 regulates RdCVF promoter activity in the inner retina.

Rod-derived Cone Viability Factor (RdCVF) is a trophic factor with therapeutic potential for the treatment of retinitis pigmentosa, a retinal disease that commonly results in blindness. RdCVF is encoded by Nucleoredoxin-like 1 (Nxnl1), a gene homologous with the family of thioredoxins that participate in the defense against oxidative stress. RdCVF expression is lost after rod degeneration in the first phase of retinitis pigmentosa, and this loss has been implicated in the more clinically significant secondary cone degeneration that often occurs. Here, we describe a study of the Nxnl1 promoter using an approach that combines promoter and transcriptomic analysis. By transfection of selected candidate transcription factors, chosen based upon their expression pattern, we identified the homeodomain proteins CHX10/VSX2, VSX1 and PAX4, as well as the zinc finger protein SP3, as factors that can stimulate both the mouse and human Nxnl1 promoter. In addition, CHX10/VSX2 binds to the Nxnl1 promoter in vivo. Since CHX10/VSX2 is expressed predominantly in the inner retina, this finding motivated us to demonstrate that RdCVF is expressed in the inner as well as the outer retina. Interestingly, the loss of rods in the rd1 mouse, a model of retinitis pigmentosa, is associated with decreased expression of RdCVF by inner retinal cells as well as by rods. Based upon these results, we propose an alternative therapeutic strategy aimed at recapitulating RdCVF expression in the inner retina, where cell loss is not significant, to prevent secondary cone death and central vision loss in patients suffering from retinitis pigmentosa.

Global human tissue profiling and protein network analysis reveals distinct levels of transcriptional germline-specificity and identifies target genes for male infertility.

BACKGROUND: Mammalian spermatogenesis is a process that involves a complex expression program in both somatic and germ cells present in the male gonad. A number of studies have attempted to define the transcriptome of male meiosis and gametogenesis in rodents and primates. Few human transcripts, however, have been associated with testicular somatic cells and germ cells at different post-natal developmental stages and little is known about their level of germline-specificity compared with non-testicular tissues. METHODS: We quantified human transcripts using GeneChips and a total of 47 biopsies from prepubertal children diagnosed with undescended testis, infertile adult patients whose spermatogenesis is arrested at consecutive stages and fertile control individuals. These results were integrated with data from enriched normal germ cells, non-testicular expression data, phenotype information, predicted regulatory DNA-binding motifs and interactome data. RESULTS: Among 3580 genes for which we found differential transcript concentrations in somatic and germ cells present in human testis, 933 were undetectable in 45 embryonic and adult non-testicular tissues, including many that were corroborated at protein level by published gene annotation data and histological high-throughput protein immunodetection assays. Using motif enrichment analyses, we identified regulatory promoter elements likely involved in germline development. Finally, we constructed a regulatory disease network for human fertility by integrating expression signals, interactome information, phenotypes and functional annotation data. CONCLUSIONS: Our results provide broad insight into the post-natal human testicular transcriptome at the level of cell populations and in a global somatic tissular context. Furthermore, they yield clues for genetic causes of male infertility and will facilitate the identification of novel cancer/testis genes as targets for cancer immunotherapies.

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.

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.