KRAS Promotes GLI2-Dependent Transcription during Pancreatic Carcinogenesis.

Aberrant activation of GLI transcription factors has been implicated in the pathogenesis of different tumor types including pancreatic ductal adenocarcinoma. However, the mechanistic link with established drivers of this disease remains in part elusive. In this study, using a new genetically engineered mouse model overexpressing constitutively active mouse form of GLI2 and a combination of genome-wide assays, we provide evidence of a novel mechanism underlying the interplay between KRAS, a major driver of pancreatic ductal adenocarcinoma development, and GLI2 to control oncogenic gene expression. These mice, also expressing KrasG12D, show significantly reduced median survival rate and accelerated tumorigenesis compared with the KrasG12D only expressing mice. Analysis of the mechanism using RNA sequencing demonstrate higher levels of GLI2 targets, particularly tumor growth-promoting genes, including Ccnd1, N-Myc, and Bcl2, in KrasG12D mutant cells. Furthermore, chromatin immunoprecipitation sequencing studies showed that in these cells KrasG12D increases the levels of trimethylation of lysine 4 of the histone 3 (H3K4me3) at the promoter of GLI2 targets without affecting significantly the levels of other major active chromatin marks. Importantly, Gli2 knockdown reduces H3K4me3 enrichment and gene expression induced by mutant Kras. In summary, we demonstrate that Gli2 plays a significant role in pancreatic carcinogenesis by acting as a downstream effector of KrasG12D to control gene expression.

SUFU promotes GLI activity in a Hedgehog-independent manner in pancreatic cancer.

Aberrant activation of the Hedgehog (Hh) signaling pathway, through which the GLI family of transcription factors (TF) is stimulated, is commonly observed in cancer cells. One well-established mechanism of this increased activity is through the inactivation of Suppressor of Fused (SUFU), a negative regulator of the Hh pathway. Relief from negative regulation by SUFU facilitates GLI activity and induction of target gene expression. Here, we demonstrate a novel role for SUFU as a promoter of GLI activity in pancreatic ductal adenocarcinoma (PDAC). In non-ciliated PDAC cells unresponsive to Smoothened agonism, SUFU overexpression increases GLI transcriptional activity. Conversely, knockdown (KD) of SUFU reduces the activity of GLI in PDAC cells. Through array PCR analysis of GLI target genes, we identified B-cell lymphoma 2 (BCL2) among the top candidates down-regulated by SUFU KD. We demonstrate that SUFU KD results in reduced PDAC cell viability, and overexpression of BCL2 partially rescues the effect of reduced cell viability by SUFU KD. Further analysis using as a model GLI1, a major TF activator of the GLI family in PDAC cells, shows the interaction of SUFU and GLI1 in the nucleus through previously characterized domains. Chromatin immunoprecipitation (ChIP) assay shows the binding of both SUFU and GLI1 at the promoter of BCL2 in PDAC cells. Finally, we demonstrate that SUFU promotes GLI1 activity without affecting its protein stability. Through our findings, we propose a novel role of SUFU as a positive regulator of GLI1 in PDAC, adding a new mechanism of Hh/GLI signaling pathway regulation in cancer cells.

GLI1 interaction with p300 modulates SDF1 expression in cancer-associated fibroblasts to promote pancreatic cancer cells migration.

Carcinoma-associated fibroblasts (CAFs) play an important role in the progression of multiple malignancies. Secretion of cytokines and growth factors underlies the pro-tumoral effect of CAFs. Although this paracrine function has been extensively documented, the molecular mechanisms controlling the expression of these factors remain elusive. In this study, we provide evidence of a novel CAF transcriptional axis regulating the expression of SDF1, a major driver of cancer cell migration, involving the transcription factor GLI1 and histone acetyltransferase p300. We demonstrate that conditioned media from CAFs overexpressing GLI1 induce the migration of pancreatic cancer cells, and this effect is impaired by an SDF1-neutralizing antibody. Using a combination of co-immunoprecipitation, proximity ligation assay and chromatin immunoprecipitation assay, we further demonstrate that GLI1 and p300 physically interact in CAFs to co-occupy and drive SDF1 promoter activity. Mapping experiments highlight the requirement of GLI1 N-terminal for the interaction with p300. Importantly, knockdowns of both GLI1 and p300 reduce SDF1 expression. Further analysis shows that knockdown of GLI1 decreases SDF1 promoter activity, p300 recruitment, and levels of its associated histone marks (H4ac, H3K27ac, and H3K14ac). Finally, we show that the integrity of two GLI binding sites in the SDF1 promoter is required for p300 recruitment. Our findings define a new role for the p300-GLI1 complex in the regulation of SDF1, providing new mechanistic insight into the molecular events controlling pancreatic cancer cells migration.

GLI1, a novel target of the ER stress regulator p97/VCP, promotes ATF6f-mediated activation of XBP1.

Upon accumulation of improperly folded proteins in the Endoplasmic Reticulum (ER), the Unfolded Protein Response (UPR) is triggered to restore ER homeostasis. The induction of stress genes is a sine qua non condition for effective adaptive UPR. Although this requirement has been extensively described, the mechanisms underlying this process remain in part uncharacterized. Here, we show that p97/VCP, an AAA+ ATPase known to contribute to ER stress-induced gene expression, regulates the transcription factor GLI1, a primary effector of Hedgehog (Hh) signaling. Under basal (non-ER stress) conditions, GLI1 is repressed by a p97/VCP-HDAC1 complex while upon ER stress GLI1 is induced through a mechanism requiring both USF2 binding and increase histone acetylation at its promoter. Interestingly, the induction of GLI1 was independent of ligand-regulated Hh signaling. Further analysis showed that GLI1 cooperates with ATF6f to induce promoter activity and expression of XBP1, a key transcription factor driving UPR. Overall, our work demonstrates a novel role for GLI1 in the regulation of ER stress gene expression and defines the interplay between p97/VCP, HDAC1 and USF2 as essential players in this process.

Culture media composition influences patient-derived organoid ability to predict therapeutic responses in gastrointestinal cancers.

BACKGROUNDA patient-derived organoid (PDO) platform may serve as a promising tool for translational cancer research. In this study, we evaluated PDO's ability to predict clinical response to gastrointestinal (GI) cancers.METHODSWe generated PDOs from primary and metastatic lesions of patients with GI cancers, including pancreatic ductal adenocarcinoma, colorectal adenocarcinoma, and cholangiocarcinoma. We compared PDO response with the observed clinical response for donor patients to the same treatments.RESULTSWe report an approximately 80% concordance rate between PDO and donor tumor response. Importantly, we found a profound influence of culture media on PDO phenotype, where we showed a significant difference in response to standard-of-care chemotherapies, distinct morphologies, and transcriptomes between media within the same PDO cultures.CONCLUSIONWhile we demonstrate a high concordance rate between donor tumor and PDO, these studies also showed the important role of culture media when using PDOs to inform treatment selection and predict response across a spectrum of GI cancers.TRIAL REGISTRATIONNot applicable.FUNDINGThe Joan F. & Richard A. Abdoo Family Fund in Colorectal Cancer Research, GI Cancer program of the Mayo Clinic Cancer Center, Mayo Clinic SPORE in Pancreatic Cancer, Center of Individualized Medicine (Mayo Clinic), Department of Laboratory Medicine and Pathology (Mayo Clinic), Incyte Pharmaceuticals and Mayo Clinic Hepatobiliary SPORE, University of Minnesota-Mayo Clinic Partnership, and the Early Therapeutic program (Department of Oncology, Mayo Clinic).

Nuclear Dynamics and Chromatin Structure: Implications for Pancreatic Cancer.

Changes in nuclear shape have been extensively associated with the dynamics and functionality of cancer cells. In most normal cells, nuclei have a regular ellipsoid shape and minimal variation in nuclear size; however, an irregular nuclear contour and abnormal nuclear size is often observed in cancer, including pancreatic cancer. Furthermore, alterations in nuclear morphology have become the 'gold standard' for tumor staging and grading. Beyond the utility of altered nuclear morphology as a diagnostic tool in cancer, the implications of altered nuclear structure for the biology and behavior of cancer cells are profound as changes in nuclear morphology could impact cellular responses to physical strain, adaptation during migration, chromatin organization, and gene expression. Here, we aim to highlight and discuss the factors that regulate nuclear dynamics and their implications for pancreatic cancer biology.

Hedgehog/GLI Signaling Pathway: Transduction, Regulation, and Implications for Disease.

The Hh/GLI signaling pathway was originally discovered in Drosophila as a major regulator of segment patterning in development. This pathway consists of a series of ligands (Shh, Ihh, and Dhh), transmembrane receptors (Ptch1 and Ptch2), transcription factors (GLI1-3), and signaling regulators (SMO, HHIP, SUFU, PKA, CK1, GSK3ß, etc.) that work in concert to repress (Ptch1, Ptch2, SUFU, PKA, CK1, GSK3ß) or activate (Shh, Ihh, Dhh, SMO, GLI1-3) the signaling cascade. Not long after the initial discovery, dysregulation of the Hh/GLI signaling pathway was implicated in human disease. Activation of this signaling pathway is observed in many types of cancer, including basal cell carcinoma, medulloblastoma, colorectal, prostate, pancreatic, and many more. Most often, the activation of the Hh/GLI pathway in cancer occurs through a ligand-independent mechanism. However, in benign disease, this activation is mostly ligand-dependent. The upstream signaling component of the receptor complex, SMO, is bypassed, and the GLI family of transcription factors can be activated regardless of ligand binding. Additional mechanisms of pathway activation exist whereby the entirety of the downstream signaling pathway is bypassed, and PTCH1 promotes cell cycle progression and prevents caspase-mediated apoptosis. Throughout this review, we summarize each component of the signaling cascade, non-canonical modes of pathway activation, and the implications in human disease, including cancer.

Impact of integrated translational research on clinical exome sequencing.

PURPOSE: Exome sequencing often identifies pathogenic genetic variants in patients with undiagnosed diseases. Nevertheless, frequent findings of variants of uncertain significance necessitate additional efforts to establish causality before reaching a conclusive diagnosis. To provide comprehensive genomic testing to patients with undiagnosed disease, we established an Individualized Medicine Clinic, which offered clinical exome testing and included a Translational Omics Program (TOP) that provided variant curation, research activities, or research exome sequencing. METHODS: From 2012 to 2018, 1101 unselected patients with undiagnosed diseases received exome testing. Outcomes were reviewed to assess impact of the TOP and patient characteristics on diagnostic rates through descriptive and multivariate analyses. RESULTS: The overall diagnostic yield was 24.9% (274 of 1101 patients), with 174 (15.8% of 1101) diagnosed on the basis of clinical exome sequencing alone. Four hundred twenty-three patients with nondiagnostic or without access to clinical exome sequencing were evaluated by the TOP, with 100 (9% of 1101) patients receiving a diagnosis, accounting for 36.5% of the diagnostic yield. The identification of a genetic diagnosis was influenced by the age at time of testing and the disease phenotype of the patient. CONCLUSION: Integration of translational research activities into clinical practice of a tertiary medical center can significantly increase the diagnostic yield of patients with undiagnosed disease.

GLI1/GLI2 functional interplay is required to control Hedgehog/GLI targets gene expression.

The Hedgehog-regulated transcription factors GLI1 and GLI2 play overlapping roles in development and disease; however, the mechanisms underlying their interplay remain elusive. We report for the first time that GLI1 and GLI2 physically and functionally interact in cancer cells. GLI1 and GLI2 were shown to co-immunoprecipitate in PANC1 pancreatic cancer cells and RMS13 rhabdomyosarcoma cells. Mapping analysis demonstrated that the zinc finger domains of both proteins are required for their heteromerization. RNAi knockdown of either GLI1 or GLI2 inhibited expression of many well-characterized GLI target genes (BCL2, MYCN, PTCH2, IL7 and CCND1) in PANC1 cells, whereas PTCH1 expression was only inhibited by GLI1 depletion. qPCR screening of a large set of putative canonical and non-canonical Hedgehog/GLI targets identified further genes (e.g. E2F1, BMP1, CDK2) strongly down-regulated by GLI1 and/or GLI2 depletion in PANC1 cells, and demonstrated that ANO1, AQP1 and SOCS1 are up-regulated by knockdown of either GLI1 or GLI2. Chromatin immunoprecipitation showed that GLI1 and GLI2 occupied the same regions at the BCL2, MYCN and CCND1 promoters. Furthermore, depletion of GLI1 inhibited GLI2 occupancy at these promoters, suggesting that GLI1/GLI2 interaction is required for the recruitment of GLI2 to these sites. Together, these findings indicate that GLI1 and GLI2 co-ordinately regulate the transcription of some genes, and provide mechanistic insight into the roles of GLI proteins in carcinogenesis.

TFII-I-mediated polymerase pausing antagonizes GLI2 induction by TGFß.

The modulation of GLI2, an oncogenic transcription factor commonly upregulated in cancer, is in many cases not due to genetic defects, suggesting dysregulation through alternative mechanisms. The identity of these molecular events remains for the most part unknown. Here, we identified TFII-I as a novel repressor of GLI2 expression. Mapping experiments suggest that the INR region of the GLI2 promoter is necessary for GLI2 repression. ChIP studies showed that TFII-I binds to this INR. TFII-I knockdown decreased the binding of NELF-A, a component of the promoter-proximal pausing complex at this site, and enriched phosphorylated RNAPII serine 2 in the GLI2 gene body. Immunoprecipitation studies demonstrate TFII-I interaction with SPT5, another pausing complex component. TFII-I overexpression antagonized GLI2 induction by TGFß, a known activator of GLI2 in cancer cells. TGFß reduced endogenous TFII-I binding to the INR and increased RNAPII SerP2 in the gene body. We demonstrate that this regulatory mechanism is not exclusive of GLI2. TGFß-induced genes CCR7, TGFß1 and EGR3 showed similar decreased TFII-I and NELF-A INR binding and increased RNAPII SerP2 in the gene body post-TGFß treatment. Together these results identify TFII-I as a novel repressor of a subset of TGFß-responsive genes through the regulation of RNAPII pausing.

De novo variants of NR4A2 are associated with neurodevelopmental disorder and epilepsy.

PURPOSE: This study characterizes the clinical and genetic features of nine unrelated patients with de novo variants in the NR4A2 gene. METHODS: Variants were identified and de novo origins were confirmed through trio exome sequencing in all but one patient. Targeted RNA sequencing was performed for one variant to confirm its splicing effect. Independent discoveries were shared through GeneMatcher. RESULTS: Missense and loss-of-function variants in NR4A2 were identified in patients from eight unrelated families. One patient carried a larger deletion including adjacent genes. The cases presented with developmental delay, hypotonia (six cases), and epilepsy (six cases). De novo status was confirmed for eight patients. One variant was demonstrated to affect splicing and result in expression of abnormal transcripts likely subject to nonsense-mediated decay. CONCLUSION: Our study underscores the importance of NR4A2 as a disease gene for neurodevelopmental disorders and epilepsy. The identified variants are likely causative of the seizures and additional developmental phenotypes in these patients.

Molecular characterization of known and novel ACVR1 variants in phenotypes of aberrant ossification.

Diffuse idiopathic skeletal hyperostosis (DISH) is a disorder principally characterized by calcification and ossification of spinal ligaments and entheses. Fibrodysplasia ossificans progressiva (FOP) is a rare autosomal dominant disabling disorder characterized by progressive ossification of skeletal muscle, fascia, tendons, and ligaments. These conditions manifest phenotypic overlap in the ossification of tendons and ligaments. We describe herein a patient with DISH, exhibiting heterotopic ossification of the posterior longitudinal ligament where clinical whole exome sequencing identified a variant within ACVR1, a gene implicated in FOP. This variant, p.K400E, is a novel variant, not identified previously, and occurs in a highly conserved region across orthologs. We used sequence-based predicative algorithms, molecular modeling, and molecular dynamics simulations, to test the potential for p.K400E to alter the structure and dynamics of ACVR1. We applied the same modeling and simulation methods to established FOP variants, to identify the detailed effects that they have on the ACVR1 protein, as well as to act as positive controls against which the effects of p.K400E could be evaluated. Our in silico molecular analyses support p.K400E as altering the behavior of ACVR1. In addition, functional testing to measure the effect of this variant on BMP-pSMAD 1/5/8 target genes was carried out which revealed this variant to cause increased ID1 and Msx2 expression compared with the wild-type receptor. This analysis supports the potential for the variant of uncertain significance to contribute to the patient's phenotype.

Novel zebrafish behavioral assay to identify modifiers of the rapid, nongenomic stress response.

When vertebrates face acute stressors, their bodies rapidly undergo a repertoire of physiological and behavioral adaptations, which is termed the stress response. Rapid changes in heart rate and blood glucose levels occur via the interaction of glucocorticoids and their cognate receptors following hypothalamic-pituitary-adrenal axis activation. These physiological changes are observed within minutes of encountering a stressor and the rapid time domain rules out genomic responses that require gene expression changes. Although behavioral changes corresponding to physiological changes are commonly observed, it is not clearly understood to what extent hypothalamic-pituitary-adrenal axis activation dictates adaptive behavior. We hypothesized that rapid locomotor response to acute stressors in zebrafish requires hypothalamic-pituitary-interrenal (HPI) axis activation. In teleost fish, interrenal cells are functionally homologous to the adrenocortical layer. We derived eight frameshift mutants in genes involved in HPI axis function: two mutants in exon 2 of mc2r (adrenocorticotropic hormone receptor), five in exon 2 or 5 of nr3c1 (glucocorticoid receptor [GR]) and two in exon 2 of nr3c2 (mineralocorticoid receptor [MR]). Exposing larval zebrafish to mild environmental stressors, acute changes in salinity or light illumination, results in a rapid locomotor response. We show that this locomotor response requires a functioning HPI axis via the action of mc2r and the canonical GR encoded by nr3c1 gene, but not MR (nr3c2). Our rapid behavioral assay paradigm based on HPI axis biology can be used to screen for genetic and environmental modifiers of the hypothalamic-pituitary-adrenal axis and to investigate the effects of corticosteroids and their cognate receptor interactions on behavior.

PCNT point mutations and familial intracranial aneurysms.

OBJECTIVE: To identify novel genes involved in the etiology of intracranial aneurysms (IAs) or subarachnoid hemorrhages (SAHs) using whole-exome sequencing. METHODS: We performed whole-exome sequencing in 13 individuals from 3 families with an autosomal dominant IA/SAH inheritance pattern to look for candidate genes for disease. In addition, we sequenced PCNT exon 38 in a further 161 idiopathic patients with IA/SAH to find additional carriers of potential pathogenic variants. RESULTS: We identified 2 different variants in exon 38 from the PCNT gene shared between affected members from 2 different families with either IA or SAH (p.R2728C and p.V2811L). One hundred sixty-four samples with either SAH or IA were Sanger sequenced for the PCNT exon 38. Five additional missense mutations were identified. We also found a second p.V2811L carrier in a family with a history of neurovascular diseases. CONCLUSION: The PCNT gene encodes a protein that is involved in the process of microtubule nucleation and organization in interphase and mitosis. Biallelic loss-of-function mutations in PCNT cause a form of primordial dwarfism (microcephalic osteodysplastic primordial dwarfism type II), and ≈50% of these patients will develop neurovascular abnormalities, including IAs and SAHs. In addition, a complete Pcnt knockout mouse model (Pcnt -/-) published previously showed general vascular abnormalities, including intracranial hemorrhage. The variants in our families lie in the highly conserved PCNT protein-protein interaction domain, making PCNT a highly plausible candidate gene in cerebrovascular disease.

Utility of DNA, RNA, Protein, and Functional Approaches to Solve Cryptic Immunodeficiencies.

PURPOSE: We report a female infant identified by newborn screening for severe combined immunodeficiencies (NBS SCID) with T cell lymphopenia (TCL). The patient had persistently elevated alpha-fetoprotein (AFP) with IgA deficiency, and elevated IgM. Gene sequencing for a SCID panel was uninformative. We sought to determine the cause of the immunodeficiency in this infant. METHODS: We performed whole-exome sequencing (WES) on the patient and parents to identify a genetic diagnosis. Based on the WES result, we developed a novel flow cytometric panel for rapid assessment of DNA repair defects using blood samples. We also performed whole transcriptome sequencing (WTS) on fibroblast RNA from the patient and father for abnormal transcript analysis. RESULTS: WES revealed a pathogenic paternally inherited indel in ATM. We used the flow panel to assess several proteins in the DNA repair pathway in lymphocyte subsets. The patient had absent phosphorylation of ATM, resulting in absent or aberrant phosphorylation of downstream proteins, including γH2AX. However, ataxia-telangiectasia (AT) is an autosomal recessive condition, and the abnormal functional data did not correspond with a single ATM variant. WTS revealed in-frame reciprocal fusion transcripts involving ATM and SLC35F2 indicating a chromosome 11 inversion within 11q22.3, of maternal origin. Inversion breakpoints were identified within ATM intron 16 and SLC35F2 intron 7. CONCLUSIONS: We identified a novel ATM-breaking chromosome 11 inversion in trans with a pathogenic indel (compound heterozygote) resulting in non-functional ATM protein, consistent with a diagnosis of AT. Utilization of several molecular and functional assays allowed successful resolution of this case.

The endocannabinoid gene faah2a modulates stress-associated behavior in zebrafish.

The ability to orchestrate appropriate physiological and behavioral responses to stress is important for survival, and is often dysfunctional in neuropsychiatric disorders that account for leading causes of global disability burden. Numerous studies have shown that the endocannabinoid neurotransmitter system is able to regulate stress responses and could serve as a therapeutic target for the management of these disorders. We used quantitative reverse transcriptase-polymerase chain reactions to show that genes encoding enzymes that synthesize (abhd4, gde1, napepld), enzymes that degrade (faah, faah2a, faah2b), and receptors that bind (cnr1, cnr2, gpr55-like) endocannabinoids are expressed in zebrafish (Danio rerio). These genes are conserved in many other vertebrates, including humans, but fatty acid amide hydrolase 2 has been lost in mice and rats. We engineered transcription activator-like effector nucleases to create zebrafish with mutations in cnr1 and faah2a to test the role of these genes in modulating stress-associated behavior. We showed that disruption of cnr1 potentiated locomotor responses to hyperosmotic stress. The increased response to stress was consistent with rodent literature and served to validate the use of zebrafish in this field. Moreover, we showed for the first time that disruption of faah2a attenuated the locomotor responses to hyperosmotic stress. This later finding suggests that FAAH2 may be an important mediator of stress responses in non-rodent vertebrates. Accordingly, FAAH and FAAH2 modulators could provide distinct therapeutic options for stress-aggravated disorders.

Molecular Modeling and Functional Analysis of Exome Sequencing-Derived Variants of Unknown Significance Identify a Novel, Constitutively Active FGFR2 Mutant in Cholangiocarcinoma.

PURPOSE: Genomic testing has increased the quantity of information available to oncologists. Unfortunately, many identified sequence alterations are variants of unknown significance (VUSs), which thus limit the clinician's ability to use these findings to inform treatment. We applied a combination of in silico prediction and molecular modeling tools and laboratory techniques to rapidly define actionable VUSs. MATERIALS AND METHODS: Exome sequencing was conducted on 308 tumors from various origins. Most single nucleotide alterations within gene coding regions were VUSs. These VUSs were filtered to identify a subset of therapeutically targetable genes that were predicted with in silico tools to be altered in function by their variant sequence. A subset of receptor tyrosine kinase VUSs was characterized by laboratory comparison of each VUS versus its wild-type counterpart in terms of expression and signaling activity. RESULTS: The study identified 4,327 point mutations of which 3,833 were VUSs. Filtering for mutations in genes that were therapeutically targetable and predicted to affect protein function reduced these to 522VUSs of interest, including a large number of kinases. Ten receptortyrosine kinase VUSs were selected to explore in the laboratory. Of these, seven were found to be functionally altered. Three VUSs (FGFR2 F276C, FGFR4 R78H, and KDR G539R) showed increased basal or ligand-stimulated ERK phosphorylation compared with their wild-type counterparts, which suggests that they support transformation. Treatment of a patient who carried FGFR2 F276C with an FGFR inhibitor resulted in significant and sustained tumor response with clinical benefit. CONCLUSION: The findings demonstrate the feasibility of rapid identification of the biologic relevance of somatic mutations, which thus advances clinicians' ability to make informed treatment decisions.