RUNX transcription factors are essential in maintaining epididymal epithelial differentiation.

Apart from the androgen receptor, transcription factors (TFs) that are required for the development and formation of the different segments of the epididymis have remained unknown. We identified TF families expressed in the developing epididymides, of which many showed segment specificity. From these TFs, down-regulation of runt related transcription factors (RUNXs) 1 and 2 expression coincides with epithelial regression in Dicer1 cKO mice. Concomitant deletion of both Runx1 and Runx2 in a mouse epididymal epithelial cell line affected cell morphology, adhesion and mobility in vitro. Furthermore, lack of functional RUNXs severely disturbed the formation of 3D epididymal organoid-like structures. Transcriptomic analysis of the epididymal cell organoid-like structures indicated that RUNX1 and RUNX2 are involved in the regulation of MAPK signaling, NOTCH pathway activity, and EMT-related gene expression. This suggests that RUNXs are master regulators of several essential signaling pathways, and necessary for the maintenance of proper differentiation of the epididymal epithelium.

Epigenetic control of CD1D expression as a mechanism of resistance to immune checkpoint therapy in poorly immunogenic melanomas.

Immune Checkpoint Therapies (ICT) have revolutionized the treatment of metastatic melanoma. However, only a subset of patients reaches complete responses. Deficient ß2-microglobulin (ß2M) expression impacts antigen presentation to T cells, leading to ICT resistance. Here, we investigate alternative ß2M-correlated biomarkers that associate with ICT resistance. We shortlisted immune biomarkers interacting with human ß2M using the STRING database. Next, we profiled the transcriptomic expression of these biomarkers in association with clinical and survival outcomes in the melanoma GDC-TCGA-SKCM dataset and a collection of publicly available metastatic melanoma cohorts treated with ICT (anti-PD1). Epigenetic control of identified biomarkers was interrogated using the Illumina Human Methylation 450 dataset from the melanoma GDC-TCGA-SKCM study. We show that ß2M associates with CD1d, CD1b, and FCGRT at the protein level. Co-expression and correlation profile of B2M with CD1D, CD1B, and FCGRT dissociates in melanoma patients following B2M expression loss. Lower CD1D expression is typically found in patients with poor survival outcomes from the GDC-TCGA-SKCM dataset, in patients not responding to anti-PD1 immunotherapies, and in a resistant anti-PD1 pre-clinical model. Immune cell abundance study reveals that B2M and CD1D are both enriched in tumor cells and dendritic cells from patients responding to anti-PD1 immunotherapies. These patients also show increased levels of natural killer T (NKT) cell signatures in the tumor microenvironment (TME). Methylation reactions in the TME of melanoma impact the expression of B2M and SPI1, which controls CD1D expression. These findings suggest that epigenetic changes in the TME of melanoma may impact ß2M and CD1d-mediated functions, such as antigen presentation for T cells and NKT cells. Our hypothesis is grounded in comprehensive bioinformatic analyses of a large transcriptomic dataset from four clinical cohorts and mouse models. It will benefit from further development using well-established functional immune assays to support understanding the molecular processes leading to epigenetic control of ß2M and CD1d. This research line may lead to the rational development of new combinatorial treatments for metastatic melanoma patients that poorly respond to ICT.

Novel prognostication biomarker adipophilin reveals a metabolic shift in uveal melanoma and new therapeutic opportunities.

Metastatic uveal melanoma remains incurable at present. We previously demonstrated that loss of BAP1 gene expression in tumour cells triggers molecular mechanisms of immunosuppression in the tumour microenvironment (TME) of metastatic uveal melanoma. Adipophilin is a structural protein of lipid droplets involved in fat storage within mammalian cells, and its expression has been identified in uveal melanoma. We comprehensively evaluated adipophilin expression at the RNA (PLIN2) and protein levels of 80 patients of the GDC-TCGA-UM study and in a local cohort of 43 primary uveal melanoma samples respectively. PLIN2 expression is a survival prognosticator biomarker in uveal melanoma. Loss of adipophilin expression is significantly associated with monosomy 3 status and nuclear BAP1 losses in uveal melanoma tumours. Integrative transcriptomic and secretome studies show a relationship between transient loss of adipophilin expression and increased levels of tumour-associated macrophages and hypoxia genes, suggesting PLIN2-dependent changes in oxygen and lipid metabolism in the TME of low and high-metastatic risk uveal melanoma. We designed four adipophilin-based multigene signatures for uveal melanoma prognostication using a transcriptomic and secretome survival-functional network approach. Adipophilin-based multigene signatures were validated in BAP1-positive and BAP1-negative uveal melanoma cell lines using a next-generation RNA sequencing approach. We identified existing small molecules, mostly adrenergic, retinoid, and glucocorticoid receptor agonists, MEK, and RAF inhibitors, with the potential to reverse this multigene signature expression in uveal melanoma. Some of these molecules were able to impact tumour cell viability, and carvedilol, an adrenergic receptor antagonist, restored PLIN2 levels, mimicking the expression of normoxia/lipid storage signatures and reversing the expression of hypoxia/lipolysis signatures in co-cultures of uveal melanoma cells with human macrophages. These findings open up a new research line for understanding the lipid metabolic regulation of immune responses, with implications for therapeutic innovation in uveal melanoma. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

High intratumoral dihydrotestosterone is associated with antiandrogen resistance in VCaP prostate cancer xenografts in castrated mice.

Antiandrogen treatment resistance is a major clinical concern in castration-resistant prostate cancer (CRPC) treatment. Using xenografts of VCaP cells we showed that growth of antiandrogen resistant CRPC tumors were characterized by a higher intratumor dihydrotestosterone (DHT) concentration than that of treatment responsive tumors. Furthermore, the slow tumor growth after adrenalectomy was associated with a low intratumor DHT concentration. Reactivation of androgen signaling in enzalutamide-resistant tumors was further shown by the expression of several androgen-dependent genes. The data indicate that intratumor DHT concentration and expression of several androgen-dependent genes in CRPC lesions is an indication of enzalutamide treatment resistance and an indication of the need for further androgen blockade. The presence of an androgen synthesis, independent of CYP17A1 activity, has been shown to exist in prostate cancer cells, and thus, novel androgen synthesis inhibitors are needed for the treatment of enzalutamide-resistant CRPC tumors that do not respond to abiraterone.

Selective loss of kisspeptin signaling in oocytes causes progressive premature ovulatory failure.

STUDY QUESTION: Does direct kisspeptin signaling in the oocyte have a role in the control of follicular dynamics and ovulation? SUMMARY ANSWER: Kisspeptin signaling in the oocyte plays a relevant physiological role in the direct control of ovulation; oocyte-specific ablation of kisspeptin receptor, Gpr54, induces a state of premature ovulatory failure in mice that recapitulates some features of premature ovarian insufficiency (POI). WHAT IS KNOWN ALREADY: Kisspeptins, encoded by the Kiss1 gene, are essential for the control of ovulation and fertility, acting primarily on hypothalamic GnRH neurons to stimulate gonadotropin secretion. However, kisspeptins and their receptor, Gpr54, are also expressed in the ovary of different mammalian species, including humans, where their physiological roles remain contentious and poorly characterized. STUDY DESIGN, SIZE, DURATION: A novel mouse line with conditional ablation of Gpr54 in oocytes, named OoGpr54-/-, was generated and studied in terms of follicular and ovulatory dynamics at different age-points of postnatal maturation. A total of 59 OoGpr54-/- mice and 47 corresponding controls were analyzed. In addition, direct RNA sequencing was applied to ovarian samples from 8 OoGpr54-/- and 7 control mice at 6 months of age, and gonadotropin priming for ovulatory induction was conducted in mice (N = 7) from both genotypes. PARTICIPANTS/MATERIALS, SETTING, METHODS: Oocyte-selective ablation of Gpr54 in the oocyte was achieved in vivo by crossing a Gdf9-driven Cre-expressing transgenic mouse line with a Gpr54 LoxP mouse line. The resulting OoGpr54-/- mouse line was subjected to phenotypic, histological, hormonal and molecular analyses at different age-points of postnatal maturation (Day 45, and 2, 4, 6 and 10-11 months of age), in order to characterize the timing of puberty, ovarian follicular dynamics and ovulation, with particular attention to identification of features reminiscent of POI. The molecular signature of ovaries from OoGpr54-/- mice was defined by direct RNA sequencing. Ovulatory responses to gonadotropin priming were also assessed in OoGpr54-/- mice. MAIN RESULTS AND THE ROLE OF CHANCE: Oocyte-specific ablation of Gpr54 caused premature ovulatory failure, with some POI-like features. OoGpr54-/- mice had preserved puberty onset, without signs of hypogonadism. However, already at 2 months of age, 40% of OoGpr54-/- females showed histological features reminiscent of ovarian failure and anovulation. Penetrance of the phenotype progressed with age, with >80% and 100% of OoGpr54-/- females displaying complete ovulatory failure by 6- and 10 months, respectively. This occurred despite unaltered hypothalamic Gpr54 expression and gonadotropin levels. Yet, OoGpr54-/- mice had decreased sex steroid levels. While the RNA signature of OoGpr54-/- ovaries was dominated by the anovulatory state, oocyte-specific ablation of Gpr54 significantly up- or downregulated of a set of 21 genes, including those encoding pituitary adenylate cyclase-activating polypeptide, Wnt-10B, matrix-metalloprotease-12, vitamin A-related factors and calcium-activated chloride channel-2, which might contribute to the POI-like state. Notably, the anovulatory state of young OoGpr54-/- mice could be rescued by gonadotropin priming. LARGE SCALE DATA: N/A. . LIMITATIONS, REASONS FOR CAUTION: Conditional ablation of Gpr54 in oocytes unambiguously caused premature ovulatory failure in mice; yet, the ultimate molecular mechanisms for such state of POI can be only inferred on the basis of RNAseq data and need further elucidation, since some of the molecular changes observed in OoGpr54-/- ovaries were secondary to the anovulatory state. Direct translation of mouse findings to human disease should be made with caution since, despite the conserved expression of Kiss1/kisspeptin and Gpr54 in rodents and humans, our mouse model does not recapitulate all features of common forms of POI. WIDER IMPLICATIONS OF THE FINDINGS: Deregulation of kisspeptin signaling in the oocyte might be an underlying, and previously unnoticed, cause for some forms of POI in women. STUDY FUNDING/COMPETING INTEREST(S): This work was primarily supported by a grant to M.P. and M.T.-S. from the FiDiPro (Finnish Distinguished Professor) Program of the Academy of Finland. Additional financial support came from grant BFU2017-83934-P (M.T.-S.; Ministerio de Economía y Competitividad, Spain; co-funded with EU funds/FEDER Program), research funds from the IVIRMA International Award in Reproductive Medicine (M.T.-S.), and EFSD Albert Renold Fellowship Programme (S.T.R.). The authors have no conflicts of interest to declare in relation to the contents of this work. TRIAL REGISTRATION NUMBER: N/A.

Exon-level estimates improve the detection of differentially expressed genes in RNA-seq studies.

Detection of differentially expressed genes (DEGs) between different biological conditions is a key data analysis step of most RNA-sequencing studies. Conventionally, computational tools have used gene-level read counts as input to test for differential gene expression between sample condition groups. Recently, it has been suggested that statistical testing could be performed with increased power at a lower feature level prior to aggregating the results to the gene level. In this study, we systematically compared the performance of calling the DEGs when using read count data at different levels (gene, transcript, and exon) as input, in the context of two publicly available data sets. Additionally, we tested two different methods for aggregating the lower feature-level p-values to gene-level: Lancaster and empirical Brown's method. Our results show that detection of DEGs is improved compared to the conventional gene-level approach regardless of the lower feature-level used for statistical testing. The overall best balance between accuracy and false discovery rate was obtained using the exon-level approach with empirical Brown's aggregation method, which we provide as a freely available Bioconductor package EBSEA (https://bioconductor.org/packages/release/bioc/html/EBSEA.html).

Systematic evaluation of differential splicing tools for RNA-seq studies.

Differential splicing (DS) is a post-transcriptional biological process with critical, wide-ranging effects on a plethora of cellular activities and disease processes. To date, a number of computational approaches have been developed to identify and quantify differentially spliced genes from RNA-seq data, but a comprehensive intercomparison and appraisal of these approaches is currently lacking. In this study, we systematically evaluated 10 DS analysis tools for consistency and reproducibility, precision, recall and false discovery rate, agreement upon reported differentially spliced genes and functional enrichment. The tools were selected to represent the three different methodological categories: exon-based (DEXSeq, edgeR, JunctionSeq, limma), isoform-based (cuffdiff2, DiffSplice) and event-based methods (dSpliceType, MAJIQ, rMATS, SUPPA). Overall, all the exon-based methods and two event-based methods (MAJIQ and rMATS) scored well on the selected measures. Of the 10 tools tested, the exon-based methods performed generally better than the isoform-based and event-based methods. However, overall, the different data analysis tools performed strikingly differently across different data sets or numbers of samples.

Lack of androgen receptor SUMOylation results in male infertility due to epididymal dysfunction.

Androgen receptor (AR) is regulated by SUMOylation at its transactivation domain. In vitro, the SUMOylation is linked to transcriptional repression and/or target gene-selective regulation. Here, we generated a mouse model (ArKI) in which the conserved SUMO acceptor lysines of AR are permanently abolished (ArK381R, K500R). ArKI males develop normally, without apparent defects in their systemic androgen action in reproductive tissues. However, the ArKI males are infertile. Their spermatogenesis appears unaffected, but their epididymal sperm maturation is defective, shown by severely compromised motility and fertilization capacity of the sperm. Fittingly, their epididymal AR chromatin-binding and gene expression associated with sperm maturation and function are misregulated. AR is SUMOylated in the wild-type epididymis but not in the testis, which could explain the tissue-specific response to the lack of AR SUMOylation. Our studies thus indicate that epididymal AR SUMOylation is essential for the post-testicular sperm maturation and normal reproductive capability of male mice.

Adrenals Contribute to Growth of Castration-Resistant VCaP Prostate Cancer Xenografts.

The role of adrenal androgens as drivers for castration-resistant prostate cancer (CRPC) growth in humans is generally accepted; however, the value of preclinical mouse models of CRPC is debatable, because mouse adrenals do not produce steroids activating the androgen receptor. In this study, we confirmed the expression of enzymes essential for de novo synthesis of androgens in mouse adrenals, with high intratissue concentration of progesterone (P4) and moderate levels of androgens, such as androstenedione, testosterone, and dihydrotestosterone, in the adrenal glands of both intact and orchectomized (ORX) mice. ORX alone had no effect on serum P4 concentration, whereas orchectomized and adrenalectomized (ORX + ADX) resulted in a significant decrease in serum P4 and in a further reduction in the low levels of serum androgens (androstenedione, testosterone, and dihydrotestosterone), measured by mass spectrometry. In line with this, the serum prostate-specific antigen and growth of VCaP xenografts in mice after ORX + ADX were markedly reduced compared with ORX alone, and the growth difference was not abolished by a glucocorticoid treatment. Moreover, ORX + ADX altered the androgen-dependent gene expression in the tumors, similar to that recently shown for the enzalutamide treatment. These data indicate that in contrast to the current view, and similar to humans, mouse adrenals synthesize significant amounts of steroids that contribute to the androgen receptor-dependent growth of CRPC.

Intratumoral androgen levels are linked to TMPRSS2-ERG fusion in prostate cancer.

Intratumoral androgen biosynthesis is one of the mechanisms involved in the progression of prostate cancer, and an important target for novel prostate cancer therapies. Using gas chromatography-tandem mass spectrometry and genome-wide RNA sequencing, we have analyzed androgen concentrations and androgen-regulated gene expression in cancerous and morphologically benign prostate tissue specimens and serum samples obtained from 48 primary prostate cancer patients. Intratumoral dihydrotestosterone (DHT) concentrations were significantly higher in the cancerous tissues compared to benign prostate (P < 0.001). The tissue/serum ratios of androgens were highly variable between the patients, indicating individual patterns of androgen metabolism and/or uptake of androgens within the prostate tissue. An unsupervised hierarchical clustering analysis of intratissue androgen concentrations indicated that transmembrane protease, serine 2/ETS-related gene (TMPRSS2-ERG)-positive patients have different androgen profiles compared to TMPRSS2-ERG-negative patients. TMPRSS2-ERG gene fusion status was also associated with an enhanced androgen-regulated gene expression, along with altered intratumoral androgen metabolism, demonstrated by reduced testosterone concentrations and increased DHT/testosterone ratios in TMPRSS2-ERG-positive tumors. TMPRSS2-ERG-positive and -negative prostate cancer specimens have distinct intratumoral androgen profiles, possibly due to activation of testosterone-independent DHT biosynthesis via the alternative pathway in TMPRSS2-ERG-positive tumors. Thus, patients with TMPRSS2-ERG-positive prostate cancer may benefit from novel inhibitors targeting the alternative DHT biosynthesis.

Hydroxysteroid (17ß) dehydrogenase 13 deficiency triggers hepatic steatosis and inflammation in mice.

Hydroxysteroid (17ß) dehydrogenases (HSD17Bs) form an enzyme family characterized by their ability to catalyze reactions in steroid and lipid metabolism. In the present study, we characterized the phenotype of HSD17B13-knockout (HSD17B13KO) mice deficient in Hsd17b13. In these studies, hepatic steatosis was detected in HSD17B13KO male mice, indicated by histologic analysis and by the increased triglyceride concentration in the liver, whereas reproductive performance and serum steroid concentrations were normal in HSD17B13KO mice. In line with these changes, the expression of key proteins in fatty acid synthesis, such as FAS, acetyl-CoA carboxylase 1, and SCD1, was increased in the HSD17B13KO liver. Furthermore, the knockout liver showed an increase in 2 acylcarnitines, suggesting impaired mitochondrial ß-oxidation in the presence of unaltered malonyl CoA and AMPK expression. The glucose tolerance did not differ between wild-type and HSD17B13KO mice in the presence of lower levels of glucose 6-phosphatase in HSD17B13KO liver compared with wild-type liver. Furthermore, microgranulomas and increased portal inflammation together with up-regulation of immune response genes were observed in HSD17B13KO mice. Our data indicate that disruption of Hsd17b13 impairs hepatic-lipid metabolism in mice, resulting in liver steatosis and inflammation, but the enzyme does not play a major role in the regulation of reproductive functions.-Adam, M., Heikelä, H., Sobolewski, C., Portius, D., Mäki-Jouppila, J., Mehmood, A., Adhikari, P., Esposito, I., Elo, L. L., Zhang, F.-P., Ruohonen, S. T., Strauss, L., Foti, M., Poutanen, M. Hydroxysteroid (17ß) dehydrogenase 13 deficiency triggers hepatic steatosis and inflammation in mice.

Antiandrogens Reduce Intratumoral Androgen Concentrations and Induce Androgen Receptor Expression in Castration-Resistant Prostate Cancer Xenografts.

The development of castration-resistant prostate cancer (CRPC) is associated with the activation of intratumoral androgen biosynthesis and an increase in androgen receptor (AR) expression. We recently demonstrated that, similarly to the clinical CRPC, orthotopically grown castration-resistant VCaP (CR-VCaP) xenografts express high levels of AR and retain intratumoral androgen concentrations similar to tumors grown in intact mice. Herein, we show that antiandrogen treatment (enzalutamide or ARN-509) significantly reduced (10-fold, P < 0.01) intratumoral testosterone and dihydrotestosterone concentrations in the CR-VCaP tumors, indicating that the reduction in intratumoral androgens is a novel mechanism by which antiandrogens mediate their effects in CRPC. Antiandrogen treatment also altered the expression of multiple enzymes potentially involved in steroid metabolism. Identical to clinical CRPC, the expression levels of the full-length AR (twofold, P < 0.05) and the AR splice variants 1 (threefold, P < 0.05) and 7 (threefold, P < 0.01) were further increased in the antiandrogen-treated tumors. Nonsignificant effects were observed in the expression of certain classic androgen-regulated genes, such as TMPRSS2 and KLK3, despite the low levels of testosterone and dihydrotestosterone. However, other genes recently identified to be highly sensitive to androgen-regulated AR action, such as NOV and ST6GalNAc1, were markedly altered, which indicated reduced androgen action. Taken together, the data indicate that, besides blocking AR, antiandrogens modify androgen signaling in CR-VCaP xenografts at multiple levels.

Comparison of methods to detect differentially expressed genes between single-cell populations.

We compared five statistical methods to detect differentially expressed genes between two distinct single-cell populations. Currently, it remains unclear whether differential expression methods developed originally for conventional bulk RNA-seq data can also be applied to single-cell RNA-seq data analysis. Our results in three diverse comparison settings showed marked differences between the different methods in terms of the number of detections as well as their sensitivity and specificity. They, however, did not reveal systematic benefits of the currently available single-cell-specific methods. Instead, our previously introduced reproducibility-optimization method showed good performance in all comparison settings without any single-cell-specific modifications.