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.

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.

Mind bomb 1 is required for pancreatic ß-cell formation.

During early pancreatic development, Notch signaling represses differentiation of endocrine cells and promotes proliferation of Nkx6-1(+)Ptf1a(+) multipotent progenitor cells (MPCs). Later, antagonistic interactions between Nkx6 transcription factors and Ptf1a function to segregate MPCs into distal Nkx6-1(-)Ptf1a(+) acinar progenitors and proximal Nkx6-1(+)Ptf1a(-) duct and ß-cell progenitors. Distal cells are initially multipotent, but evolve into unipotent, acinar cell progenitors. Conversely, proximal cells are bipotent and give rise to duct cells and late-born endocrine cells, including the insulin producing ß-cells. However, signals that regulate proximodistal (P-D) patterning and thus formation of ß-cell progenitors are unknown. Here we show that Mind bomb 1 (Mib1) is required for correct P-D patterning of the developing pancreas and ß-cell formation. We found that endoderm-specific inactivation of Mib1 caused a loss of Nkx6-1(+)Ptf1a(-) and Hnf1ß(+) cells and a corresponding loss of Neurog3(+) endocrine progenitors and ß-cells. An accompanying increase in Nkx6-1(-)Ptf1a(+) and amylase(+) cells, occupying the proximal domain, suggests that proximal cells adopt a distal fate in the absence of Mib1 activity. Impeding Notch-mediated transcriptional activation by conditional expression of dominant negative Mastermind-like 1 (Maml1) resulted in a similarly distorted P-D patterning and suppressed ß-cell formation, as did conditional inactivation of the Notch target gene Hes1. Our results reveal iterative use of Notch in pancreatic development to ensure correct P-D patterning and adequate ß-cell formation.

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.

EVEN-SKIPPED HOMEOBOX 1 controls human ES cell differentiation by directly repressing GOOSECOID expression.

TGFß signaling patterns the primitive streak, yet little is known about transcriptional effectors that mediate the cell fate choices during streak-like development in mammalian embryos and in embryonic stem (ES) cells. Here we demonstrate that cross-antagonistic actions of EVEN-SKIPPED HOMEOBOX 1 (EVX1) and GOOSECOID (GSC) regulate cell fate decisions in streak-like progenitors derived from human ES cells exposed to BMP4 and/or activin. We found that EVX1 repressed GSC expression and promoted formation of posterior streak-like progeny in response to BMP4, and conversely that GSC repressed EVX1 expression and was required for development of anterior streak-like progeny in response to activin. Chromatin immunoprecipitation assays showed that EVX1 bound to the GSC 5'-flanking region in BMP4 treated human ES cells, and band shift assays identified two EVX1 binding sites in the GSC 5'-region. Significantly, we found that intact EVX1 binding sites were required for BMP4-mediated repression of GSC reporter constructs. We conclude that BMP4-induced EVX1 repress GSC directly and the two genes form the core of a gene regulatory network (GRN) controlling cell fates in streak-like human ES cell progeny.

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.

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.

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.

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.

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.

Bioenergetic Changes during Differentiation of Human Embryonic Stem Cells along the Hepatic Lineage.

Mitochondrial dysfunction has been demonstrated to result in premature aging due to its effects on stem cells. Nevertheless, a full understanding of the role of mitochondrial bioenergetics through differentiation is still lacking. Here we show the bioenergetics profile of human stem cells of embryonic origin differentiating along the hepatic lineage. Our study reveals especially the transition between hepatic specification and hepatic maturation as dependent on mitochondrial respiration and demonstrates that even though differentiating cells are primarily dependent on glycolysis until induction of hepatocyte maturation, oxidative phosphorylation is essential at all stages of differentiation.

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.

Evaluation in mammalian oocytes of gene transcripts linked to epigenetic reprogramming.

The mature mammalian metaphase II (MII) oocyte has a unique ability to reprogram sperm chromatin and support early embryonic development. This feature even extends to the epigenetic reprogramming of a terminally differentiated cell nucleus as observed in connection with somatic cell nuclear transfer. Epigenetic nuclear reprogramming is highly linked to chromatin structure and includes covalent modifications of DNA and core histone proteins as well as reorganization of higher-order chromatin structure. A group of conserved enzymes mediating DNA methylation, methyl-CpG-binding protein (MeCP), histone acetylation and methylation, and chromatin remodeling are extensively involved in epigenetic reprogramming in mammalian cells. Using the oligonucleotide microarray technique, the present study compared the expression levels of 86 genes associated with epigenetic reprogramming in murine in vivo matured MII oocytes with that of germinal vesicle oocytes. Correlation between biological replicates was high. A total of 57 genes with potential reprogramming effect were detected. In MII oocytes, four genes were significant up-regulated, whereas 18 were down-regulated and 35 unchanged. The significantly regulated genes were validated by real-time quantitative RT-PCR. For example, MII oocytes showed a significant down-regulation of oocyte-specific maintenance DNA methyltransferase, Dnmt1o, and up-regulation of MeCP transcript, methyl-CpG binding domain protein 2. Furthermore, histone acetyltransferases were proportionally overrepresented when compared with histone deacetylases. These data elucidate for the first time some of the mechanisms that the oocyte may employ to reprogram a foreign genome either in form of a spermatozoa or a somatic nucleus and may thus be of importance for advancing the fields of stem cell research and regenerative medicine.