The urothelial gene regulatory network: understanding biology to improve bladder cancer management.

The urothelium is a stratified epithelium composed of basal cells, one or more layers of intermediate cells, and an upper layer of differentiated umbrella cells. Most bladder cancers (BLCA) are urothelial carcinomas. Loss of urothelial lineage fidelity results in altered differentiation, highlighted by the taxonomic classification into basal and luminal tumors. There is a need to better understand the urothelial transcriptional networks. To systematically identify transcription factors (TFs) relevant for urothelial identity, we defined highly expressed TFs in normal human bladder using RNA-Seq data and inferred their genomic binding using ATAC-Seq data. To focus on epithelial TFs, we analyzed RNA-Seq data from patient-derived organoids recapitulating features of basal/luminal tumors. We classified TFs as "luminal-enriched", "basal-enriched" or "common" according to expression in organoids. We validated our classification by differential gene expression analysis in Luminal Papillary vs. Basal/Squamous tumors. Genomic analyses revealed well-known TFs associated with luminal (e.g., PPARG, GATA3, FOXA1) and basal (e.g., TP63, TFAP2) phenotypes and novel candidates to play a role in urothelial differentiation or BLCA (e.g., MECOM, TBX3). We also identified TF families (e.g., KLFs, AP1, circadian clock, sex hormone receptors) for which there is suggestive evidence of their involvement in urothelial differentiation and/or BLCA. Genomic alterations in these TFs are associated with BLCA. We uncover a TF network involved in urothelial cell identity and BLCA. We identify novel candidate TFs involved in differentiation and cancer that provide opportunities for a better understanding of the underlying biology and therapeutic intervention.

Genome-wide Association Study of Bladder Cancer Reveals New Biological and Translational Insights.

BACKGROUND: Genomic regions identified by genome-wide association studies (GWAS) for bladder cancer risk provide new insights into etiology. OBJECTIVE: To identify new susceptibility variants for bladder cancer in a meta-analysis of new and existing genome-wide genotype data. DESIGN, SETTING, AND PARTICIPANTS: Data from 32 studies that includes 13,790 bladder cancer cases and 343,502 controls of European ancestry were used for meta-analysis. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSES: Log-additive associations of genetic variants were assessed using logistic regression models. A fixed-effects model was used for meta-analysis of the results. Stratified analyses were conducted to evaluate effect modification by sex and smoking status. A polygenic risk score (PRS) was generated on the basis of known and novel susceptibility variants and tested for interaction with smoking. RESULTS AND LIMITATIONS: Multiple novel bladder cancer susceptibility loci (6p.22.3, 7q36.3, 8q21.13, 9p21.3, 10q22.1, 19q13.33) as well as improved signals in three known regions (4p16.3, 5p15.33, 11p15.5) were identified, bringing the number of independent markers at genome-wide significance (p < 5 × 10-8) to 24. The 4p16.3 (FGFR3/TACC3) locus was associated with a stronger risk for women than for men (p-interaction = 0.002). Bladder cancer risk was increased by interactions between smoking status and genetic variants at 8p22 (NAT2; multiplicative p value for interaction [pM-I] = 0.004), 8q21.13 (PAG1; pM-I = 0.01), and 9p21.3 (LOC107987026/MTAP/CDKN2A; pM-I = 0.02). The PRS based on the 24 independent GWAS markers (odds ratio per standard deviation increase 1.49, 95% confidence interval 1.44-1.53), which also showed comparable results in two prospective cohorts (UK Biobank, PLCO trial), revealed an approximately fourfold difference in the lifetime risk of bladder cancer according to the PRS (e.g., 1st vs 10th decile) for both smokers and nonsmokers. CONCLUSIONS: We report novel loci associated with risk of bladder cancer that provide clues to its biological underpinnings. Using 24 independent markers, we constructed a PRS to stratify lifetime risk. The PRS combined with smoking history, and other established risk factors, has the potential to inform future screening efforts for bladder cancer. PATIENT SUMMARY: We identified new genetic markers that provide biological insights into the genetic causes of bladder cancer. These genetic risk factors combined with lifestyle risk factors, such as smoking, may inform future preventive and screening strategies for bladder cancer.

Pseudoalignment tools as an efficient alternative to detect repeated transposable elements in scRNAseq data.

MOTIVATION: Transposable elements (TE) have played a major role in configuring the structures of mammalian genomes through evolution. In normal conditions, the expression of these elements is repressed by different epigenetic regulation mechanisms such as DNA methylation, histone modification and regulation by small RNAs. TE re-activation is associated with stemness potential acquisition, regulation of innate immunity and disease, such as cancer. However, the vast majority of current knowledge in the field is based on bulk expression studies, and very little is known on cell-type- or state-specific expression of TE-derived transcripts. Therefore, cost-efficient single-cell-resolution TE expression analytical approaches are needed. RESULTS: We have implemented an analytical approach based on pseudoalignment to consensus sequences to incorporate TE expression information to scRNAseq data. AVAILABILITY AND IMPLEMENTATION: All the data and code implemented are available as Supplementary data and in: https://github.com/jmzvillarreal/kallisto_TE_scRNAseq. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SCIRT lncRNA Restrains Tumorigenesis by Opposing Transcriptional Programs of Tumor-Initiating Cells.

In many tumors, cells transition reversibly between slow-proliferating tumor-initiating cells (TIC) and their differentiated, faster-growing progeny. Yet, how transcriptional regulation of cell-cycle and self-renewal genes is orchestrated during these conversions remains unclear. In this study, we show that as breast TIC form, a decrease in cell-cycle gene expression and increase in self-renewal gene expression are coregulated by SOX2 and EZH2, which colocalize at CpG islands. This pattern was negatively controlled by a novel long noncoding RNA (lncRNA) that we named Stem Cell Inhibitory RNA Transcript (SCIRT), which was markedly upregulated in tumorspheres but colocalized with and counteracted EZH2 and SOX2 during cell-cycle and self-renewal regulation to restrain tumorigenesis. SCIRT specifically interacted with EZH2 to increase EZH2 affinity to FOXM1 without binding the latter. In this manner, SCIRT induced transcription at cell-cycle gene promoters by recruiting FOXM1 through EZH2 to antagonize EZH2-mediated effects at target genes. Conversely, on stemness genes, FOXM1 was absent and SCIRT antagonized EZH2 and SOX2 activity, balancing toward repression. These data suggest that the interaction of an lncRNA with EZH2 can alter the affinity of EZH2 for its protein-binding partners to regulate cancer cell state transitions. SIGNIFICANCE: These findings show that a novel lncRNA SCIRT counteracts breast tumorigenesis by opposing transcriptional networks associated with cell cycle and self-renewal.See related commentary by Pardini and Dragomir, p. 535.

The GATA3 X308_Splice breast cancer mutation is a hormone context-dependent oncogenic driver.

As the catalog of oncogenic driver mutations is expanding, it becomes clear that alterations in a given gene might have different functions and should not be lumped into one class. The transcription factor GATA3 is a paradigm of this. We investigated the functions of the most common GATA3 mutation (X308_Splice) and five additional mutations, which converge into a neoprotein that we called "neoGATA3," associated with excellent prognosis in patients. Analysis of available molecular data from >3000 breast cancer patients revealed a dysregulation of the ER-dependent transcriptional response in tumors carrying neoGATA3-generating mutations. Mechanistic studies in vitro showed that neoGATA3 interferes with the transcriptional programs controlled by estrogen and progesterone receptors, without fully abrogating them. ChIP-Seq analysis indicated that ER binding is reduced in neoGATA3-expressing cells, especially at distal regions, suggesting that neoGATA3 interferes with the fine tuning of ER-dependent gene expression. This has opposite outputs in distinct hormonal context, having pro- or anti-proliferative effects, depending on the estrogen/progesterone ratio. Our data call for functional analyses of putative cancer drivers to guide clinical application.

HNF1A recruits KDM6A to activate differentiated acinar cell programs that suppress pancreatic cancer.

Defects in transcriptional regulators of pancreatic exocrine differentiation have been implicated in pancreatic tumorigenesis, but the molecular mechanisms are poorly understood. The locus encoding the transcription factor HNF1A harbors susceptibility variants for pancreatic ductal adenocarcinoma (PDAC), while KDM6A, encoding Lysine-specific demethylase 6A, carries somatic mutations in PDAC. Here, we show that pancreas-specific Hnf1a null mutant transcriptomes phenocopy those of Kdm6a mutations, and both defects synergize with KrasG12D to cause PDAC with sarcomatoid features. We combine genetic, epigenomic, and biochemical studies to show that HNF1A recruits KDM6A to genomic binding sites in pancreatic acinar cells. This remodels the acinar enhancer landscape, activates differentiated acinar cell programs, and indirectly suppresses oncogenic and epithelial-mesenchymal transition genes. We also identify a subset of non-classical PDAC samples that exhibit the HNF1A/KDM6A-deficient molecular phenotype. These findings provide direct genetic evidence that HNF1A deficiency promotes PDAC. They also connect the tumor-suppressive role of KDM6A deficiency with a cell-specific molecular mechanism that underlies PDAC subtype definition.

The EndoC-ßH1 cell line is a valid model of human beta cells and applicable for screenings to identify novel drug target candidates.

OBJECTIVE: To characterize the EndoC-ßH1 cell line as a model for human beta cells and evaluate its beta cell functionality, focusing on insulin secretion, proliferation, apoptosis and ER stress, with the objective to assess its potential as a screening platform for identification of novel anti-diabetic drug candidates. METHODS: EndoC-ßH1 was transplanted into mice for validation of in vivo functionality. Insulin secretion was evaluated in cells cultured as monolayer and as pseudoislets, as well as in diabetic mice. Cytokine induced apoptosis, glucolipotoxicity, and ER stress responses were assessed. Beta cell relevant mRNA and protein expression were investigated by qPCR and antibody staining. Hundreds of proteins or peptides were tested for their effect on insulin secretion and proliferation. RESULTS: Transplantation of EndoC-ßH1 cells restored normoglycemia in streptozotocin induced diabetic mice. Both in vitro and in vivo, we observed a clear insulin response to glucose, and, in vitro, we found a significant increase in insulin secretion from EndoC-ßH1 pseudoislets compared to monolayer cultures for both glucose and incretins. Apoptosis and ER stress were inducible in the cells and caspase 3/7 activity was elevated in response to cytokines, but not affected by the saturated fatty acid palmitate. By screening of various proteins and peptides, we found Bombesin (BB) receptor agonists and Pituitary Adenylate Cyclase-Activating Polypeptides (PACAP) to significantly induce insulin secretion and the proteins SerpinA6, STC1, and APOH to significantly stimulate proliferation. ER stress was readily induced by Tunicamycin and resulted in a reduction of insulin mRNA. Somatostatin (SST) was found to be expressed by 1% of the cells and manipulation of the SST receptors was found to significantly affect insulin secretion. CONCLUSIONS: Overall, the EndoC-ßH1 cells strongly resemble human islet beta cells in terms of glucose and incretin stimulated insulin secretion capabilities. The cell line has an active cytokine induced caspase 3/7 apoptotic pathway and is responsive to ER stress initiation factors. The cells' ability to proliferate can be further increased by already known compounds as well as by novel peptides and proteins. Based on its robust performance during the functionality assessment assays, the EndoC-ßH1 cell line was successfully used as a screening platform for identification of novel anti-diabetic drug candidates.

OCT4 and downstream factors are expressed in human somatic urogenital epithelia and in culture of epididymal spheres.

The transcription factor OCT4 plays a crucial role in the earliest differentiation of the mammalian embryo and in self-renewal of embryonic stem cells. However, it remains controversial whether this gene is also expressed in somatic tissues. Here, we use a combination of RT-PCR on whole and microdissected tissues, in situ hybridization, immunohistochemistry and western blotting to show that OCT4 and SOX2 together with downstream targets, UTF1 and REX1/ZFP42, are expressed in the human male urogenital tract. We further support these results by the analysis of DNA methylation of a region in the OCT4 promoter. In culture, human primary epididymal cells formed spheres that continued to express the investigated genes for at least 20 days. Transcriptomic analysis of cultured cells showed up-regulation of CD29, CD44 and CD133 that are normally associated with sphere-forming cancer stem cells. Furthermore, stimulation with retinoic acid resulted in down-regulation of OCT4 expression, however, without multilineage differentiation. Our results show that OCT4 and associated genes are expressed in somatic epithelial cells from the urogenital tract and that these cells can form spheres, a general marker of stem cell behaviour.

Cytokine signalling in embryonic stem cells.

Cytokines play a central role in maintaining self-renewal in mouse embryonic stem (ES) cells through a member of the interleukin-6 type cytokine family termed leukemia inhibitory factor (LIF). LIF activates the JAK-STAT3 pathway through the class I cytokine receptor gp130, which forms a trimeric complex with LIF and the class I cytokine receptor LIF receptor beta. STAT3 has been shown to play a crucial role in self-renewal in mouse ES cells probably by induction of c-myc expression. Thus, ablation of STAT3 activation leads to differentiation. However, important connections between STAT3 and other signalling pathways have been documented. In addition, gp130 activation leads to both PI3K and Src activation. The canonical Wnt pathway is sufficient to maintain self-renewal of both human ES cells and mouse ES cells. It seems quite possible that the main pathway maintaining self-renewal in ES cells is the Wnt pathway, while the LIF-JAK-STAT3 pathway is present in mouse cells as an adaptation for sustaining self-renewal during embryonic diapause, a condition of delayed implantation in mammals. In keeping with this scenario, the Wnt pathway has been shown to elevate the level of c-myc. Thus, the two pathways seem to converge on c-myc as a common target to promote self-renewal. Whereas LIF does not seem to stimulate self-renewal in human embryonic stem cells it cannot be excluded that other cytokines are involved. The pleiotropic actions of the increasing number of cytokines and receptors signalling via JAKs, STATs and SOCS exhibit considerable redundancy, compensation and plasticity in stem cells in accordance with the view that stem cells are governed by quantitative variations in strength and duration of signalling events known from other cell types rather than qualitatively different stem cell-specific factors.