Search

Abstract Maple syrup urine disease (MSUD) is a disorder of branched‐chain amino acid metabolism caused by a defect in the branched‐chain α‐ketoacid dehydrogenase (BCKD) complex (OMIM #248600). The hallmark presentation is encephalopathic crisis in neonates, but can also present with metabolic decompensation, developmental delays, and feeding difficulties. Biochemical evidence for MSUD includes elevated branched‐chain amino acids (BCAA) and the pathognomonic presence of alloisoleucine. The BCKD complex contains several subunits associated with autosomal recessive MSUD, while its regulatory proteins have less well‐defined disease associations. We report on two families with the same BCKDK variant (c.1115C>G (p.Thr372Arg)). Probands were detected on newborn screening and demonstrated biochemical evidence of MSUD. The variant was identified in reportedly asymptomatic parents and additional family members who had elevated BCAA and alloisoleucine, following an autosomal dominant pattern of inheritance. To better define the functional effect of the variant on the kinase, we completed molecular modeling using sequence‐based (2D), structural‐based (3D), and dynamic‐based (4D) analyses. The BCKDK variant modeling indicated a gain‐of‐function which leads to impaired BCAA catabolism consistent with the biochemical evidence in this cohort. Combining the evidence gained from molecular modeling with the absence of metabolic decompensation in our patients and several adult family members, despite encountering stressors typically problematic in classic MSUD, we suggest that heterozygous gain‐of‐function variants in BCKDK may represent a novel biochemical phenotype of MSUD with a benign clinical course.

Abstract Multiple myeloma (MM) is a heterogenous plasma cell malignancy, for which the established prognostic models exhibit limitations in capturing the full spectrum of outcome variability. Leveraging single-cell RNA-sequencing data, we developed a novel plasma cell gene signature. We evaluated and validated the associations of the resulting plasma cell malignancy (PBM) score with disease state, progression and clinical outcomes using data from five independent myeloma studies consisting of 2115 samples (1978 MM, 65 monoclonal gammopathy of undetermined significance, 35 smoldering MM, and 37 healthy controls). Overall, a higher PBM score was significantly associated with a more advanced stage within the spectrum of plasma cell dyscrasias (all p

Objective:. The adverse effects of ischemia-reperfusion injury (IRI) remain a principal barrier to a successful outcome after lifesaving orthotopic liver transplantation (OLT). Gene expression during different phases of IRI is dynamic and modified by individual exposures, making it attractive for identifying potential therapeutic targets for improving the number of suitable organs for transplantation and patient outcomes. However, data remain limited on the functional landscape of gene expression during liver graft IRI, spanning procurement to reperfusion and recovery. Therefore, we sought to characterize transcriptomic profiles of IRI during multiple phases in human OLT. Methods:. We conducted clinical data analyses, histologic evaluation, and RNA sequencing of 17 consecutive human primary OLT. We performed liver allograft biopsies at 4 time points: baseline (B, before donor cross-clamp), at the end of cold ischemia (CI), during early reperfusion (ER, after revascularization), and during late reperfusion (LR). Data were generated and then recipients grouped by post-OLT outcomes categories: immediate allograft function (IAF; n = 11) versus early allograft dysfunction (EAD; n = 6) groups. Results:. We observed that CI (vs B) modified a transcriptomic landscape enriched for a metabolic and immune process. Expression levels of hallmark inflammatory response genes were higher transitioning from CI to ER and decreased from ER to LR. IAF group predominantly showed higher bile and fatty acid metabolism activity during LR compared with EAD group, while EAD group maintained more immunomodulatory activities. Throughout all time points, EAD specimens exhibited decreased metabolic activity in both bile and fatty acid pathways. Conclusions:. We report transcriptomic profiles of human liver allograft IRI from prepreservation in the donor to posttransplantation in the recipient. Immunomodulatory and metabolic landscapes across ER and LR phases were different between IAF and EAD allografts. Our study also highlights marker genes for these biological processes that we plan to explore as novel therapeutic targets or surrogate markers for severe allograft injury in clinical OLT.

Introduction: The Euchromatic Histone Methyl Transferase Protein 2 (EHMT2), also known as G9a, deposits transcriptionally repressive chromatin marks that play pivotal roles in the maturation and homeostasis of multiple organs. Recently, we have shown that Ehmt2 inactivation in the mouse pancreas alters growth and immune gene expression networks, antagonizing Kras-mediated pancreatic cancer initiation and promotion. Here, we elucidate the essential role of Ehmt2 in maintaining a transcriptional landscape that protects organs from inflammation.Methods: Comparative RNA-seq studies between normal postnatal and young adult pancreatic tissue from Ehmt2 conditional knockout animals (Ehmt2fl/fl) targeted to the exocrine pancreatic epithelial cells (Pdx1-Cre and P48Cre/+), reveal alterations in gene expression networks in the whole organ related to injury-inflammation-repair, suggesting an increased predisposition to damage. Thus, we induced an inflammation repair response in the Ehmt2fl/fl pancreas and used a data science-based approach to integrate RNA-seq-derived pathways and networks, deconvolution digital cytology, and spatial transcriptomics. We also analyzed the tissue response to damage at the morphological, biochemical, and molecular pathology levels.Results and discussion: The Ehmt2fl/fl pancreas displays an enhanced injury-inflammation-repair response, offering insights into fundamental molecular and cellular mechanisms involved in this process. More importantly, these data show that conditional Ehmt2 inactivation in exocrine cells reprograms the local environment to recruit mesenchymal and immunological cells needed to mount an increased inflammatory response. Mechanistically, this response is an enhanced injury-inflammation-repair reaction with a small contribution of specific Ehmt2-regulated transcripts. Thus, this new knowledge extends the mechanisms underlying the role of the Ehmt2-mediated pathway in suppressing pancreatic cancer initiation and modulating inflammatory pancreatic diseases.

Abstract We find that NUPR1, a stress-associated intrinsically disordered protein, induced droplet formation via liquid–liquid phase separation (LLPS). NUPR1-driven LLPS was crucial for the creation of NUPR1-dependent stress granules (SGs) in pancreatic cancer cells since genetic or pharmacological inhibition by ZZW-115 of NUPR1 activity impeded SGs formation. The KrasG12D mutation induced oncogenic stress, NUPR1 overexpression, and promoted SGs development. Notably, enforced NUPR1 expression induced SGs formation independently of mutated KrasG12D. Mechanistically, KrasG12D expression strengthened sensitivity to NUPR1 inactivation, inducing cell death, activating caspase 3 and releasing LDH. Remarkably, ZZW-115-mediated SG-formation inhibition hampered the development of pancreatic intraepithelial neoplasia (PanINs) in Pdx1-cre;LSL-Kras G12D (KC) mice. ZZW-115-treatment of KC mice triggered caspase 3 activation, DNA fragmentation, and formation of the apoptotic bodies, leading to cell death, specifically in KrasG12D-expressing cells. We further demonstrated that, in developed PanINs, short-term ZZW-115 treatment prevented NUPR1-associated SGs presence. Lastly, a four-week ZZW-115 treatment significantly reduced the number and size of PanINs in KC mice. This study proposes that targeting NUPR1-dependent SGs formation could be a therapeutic approach to induce cell death in KrasG12D-dependent tumors.

Abstract Despite mitochondrial DNA (mtDNA) mutations are common events in cancer, their global frequency and clinical impact have not been comprehensively characterized in patients with myelodysplastic neoplasia (also known as myelodysplastic syndromes, MDS). Here we performed whole-genome sequencing (WGS) on samples obtained before allogenic hematopoietic cell transplantation (allo-HCT) from 494 patients with MDS who were enrolled in the Center for International Blood and Marrow Transplant Research. We evaluated the impact of mtDNA mutations on transplantation outcomes, including overall survival (OS), relapse, relapse-free survival (RFS), and transplant-related mortality (TRM). A random survival forest algorithm was applied to evaluate the prognostic performance of models that include mtDNA mutations alone and combined with MDS- and HCT-related clinical factors. A total of 2666 mtDNA mutations were identified, including 411 potential pathogenic variants. We found that overall, an increased number of mtDNA mutations was associated with inferior transplantation outcomes. Mutations in several frequently mutated mtDNA genes (e.g., MT-CYB and MT-ND5) were identified as independent predictors of OS, RFS, relapse and/or TRM after allo-HCT. Integration of mtDNA mutations into the models based on the Revised International Prognostic Scores (IPSS-R) and clinical factors related to MDS and allo-HCT could capture more prognostic information and significantly improve the prognostic stratification efforts. Our study represents the first WGS effort in MDS receiving allo-HCT and shows that there may be clinical utility of mtDNA variants to predict allo-HCT outcomes in combination with more standard clinical parameters.

Introduction: Kleefstra Syndrome type 2 (KLEFS-2) is a genetic, neurodevelopmental disorder characterized by intellectual disability, infantile hypotonia, severe expressive language delay, and characteristic facial appearance, with a spectrum of other distinct clinical manifestations. Pathogenic mutations in the epigenetic modifier type 2 lysine methyltransferase KMT2C have been identified to be causative in KLEFS-2 individuals.Methods: This work reports a translational genomic study that applies a multidimensional computational approach for deep variant phenotyping, combining conventional genomic analyses, advanced protein bioinformatics, computational biophysics, biochemistry, and biostatistics-based modeling. We use standard variant annotation, paralog annotation analyses, molecular mechanics, and molecular dynamics simulations to evaluate damaging scores and provide potential mechanisms underlying KMT2C variant dysfunction.Results: We integrated data derived from the structure and dynamics of KMT2C to classify variants into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). When compared with controls, these variants show values reflecting alterations in molecular fitness in both structure and dynamics.Discussion: We demonstrate that our 3D models for KMT2C variants suggest distinct mechanisms that lead to their imbalance and are not predictable from sequence alone. Thus, the missense variants studied here cause destabilizing effects on KMT2C function by different biophysical and biochemical mechanisms which we adeptly describe. This new knowledge extends our understanding of how variations in the KMT2C gene cause the dysfunction of its methyltransferase enzyme product, thereby bearing significant biomedical relevance for carriers of KLEFS2-associated genomic mutations.

Summary: Interpreting genetic changes observed in individual patients is a critical challenge. The array of immune deficiency syndromes is typically caused by genetic variation unique to individuals. Therefore, new approaches are needed to interpret functional variation and accelerate genomics interpretation. We constructed the first full-length structural model of human RAG recombinase across four functional states of the recombination process. We functionally tested 182 clinically observed RAG missense mutations. These experiments revealed dysfunction due to recombinase dysfunction and altered chromatin interactions. Structural modeling identified mechanical and energetic roles for each mutation. We built regression models for RAG1 (R2 = 0.91) and RAG2 (R2 = 0.97) to predict RAG activity changes. We applied our model to 711 additional RAG variants observed in population studies and identified a subset that may impair RAG function. Thus, we demonstrated a fundamental advance in the mechanistic interpretation of human genetic variations spanning from rare and undiagnosed diseases to population health.

Current capabilities in genomic sequencing outpace functional interpretations. Our previous work showed that 3D protein structure calculations enhance mechanistic understanding of genetic variation in sequenced tumors and patients with rare diseases. The KRAS GTPase is among the critical genetic factors driving cancer and germline conditions. Because KRAS-altered tumors frequently harbor one of three classic hotspot mutations, nearly all studies have focused on these mutations, leaving significant functional ambiguity across the broader KRAS genomic landscape observed in cancer and non-cancer diseases. Herein, we extend structural bioinformatics with molecular simulations to study an expanded landscape of 86 KRAS mutations. We identify multiple coordinated changes strongly associated with experimentally established KRAS biophysical and biochemical properties. The patterns we observe span hotspot and non-hotspot alterations, which can all dysregulate Switch regions, producing mutation-restricted conformations with different effector binding propensities. We experimentally measured mutation thermostability and identified shared and distinct patterns with simulations. Our results indicate mutation-specific conformations, which show potential for future research into how these alterations reverberate into different molecular and cellular functions. The data we present is not predictable using current genomic tools, demonstrating the added functional information derived from molecular simulations for interpreting human genetic variation.

This study investigates the functional significance of assorted variants of uncertain significance (VUS) in euchromatic histone lysine methyltransferase 1 (EHMT1), which is critical for early development and normal physiology. EHMT1 mutations cause Kleefstra syndrome and are linked to various human cancers. However, accurate functional interpretations of these variants are yet to be made, limiting diagnoses and future research. To overcome this, we integrate conventional tools for variant calling with computational biophysics and biochemistry to conduct multi-layered mechanistic analyses of the SET catalytic domain of EHMT1, which is critical for this protein function. We use molecular mechanics and molecular dynamics (MD)-based metrics to analyze the SET domain structure and functional motions resulting from 97 Kleefstra syndrome missense variants within the domain. Our approach allows us to classify the variants in a mechanistic manner into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). Our findings reveal that the damaging variants are mostly mapped around the active site, substrate binding site, and pre-SET regions. Overall, we report an improvement for this method over conventional tools for variant interpretation and simultaneously provide a molecular mechanism for variant dysfunction.

Cancer-associated fibroblasts (CAFs) are the main component of the stroma in pancreatic cancer, but their tissue of origin remains to be defined. Here the authors perform lineage tracing and single cell RNA sequencing in mice and suggest the splanchnic mesenchyme as the tissue of origin for CAFs.

Summary: Background: Systems biology leveraging multi-OMICs technologies, is rapidly advancing development of precision therapies and matching patients to targeted therapies, leading to improved responses. A new pillar of precision oncology lies in the power of chemogenomics to discover drugs that sensitizes malignant cells to other therapies. Here, we test a chemogenomic approach using epigenomic inhibitors (epidrugs) to reset patterns of gene expression driving the malignant behavior of pancreatic tumors. Methods: We tested a targeted library of ten epidrugs targeting regulators of enhancers and super-enhancers on reprogramming gene expression networks in seventeen patient-derived primary pancreatic cancer cell cultures (PDPCCs), of both basal and classical subtypes. We subsequently evaluated the ability of these epidrugs to sensitize pancreatic cancer cells to five chemotherapeutic drugs that are clinically used for this malignancy. Findings: To comprehend the impact of epidrug priming at the molecular level, we evaluated the effect of each epidrugs at the transcriptomic level of PDPCCs. The activating epidrugs showed a higher number of upregulated genes than the repressive epidrugs (χ2 test p-value

The PI3K enzymes modify phospholipids to regulate cell growth and differentiation. Somatic variants in PI3K are recurrent in cancer and drive a proliferative phenotype. Somatic mosaicism of PIK3R1 and PIK3CA are associated with vascular anomalies and overgrowth syndromes. Germline PIK3R1 variants are associated with varying phenotypes, including immunodeficiency or facial dysmorphism with growth delay, lipoatrophy, and insulin resistance associated with SHORT syndrome. There has been limited study of the molecular mechanism to unify our understanding of how variants in PIK3R1 drive both undergrowth and overgrowth phenotypes. Thus, we compiled genomic variants from cancer and rare vascular anomalies and sought to interpret their effects using an unbiased physics-based simulation approach for the protein complex. We applied molecular dynamics simulations to mechanistically understand how genetic variants affect PIK3R1 and its interactions with PIK3CA. Notably, iSH2 genetic variants associated with undergrowth destabilize molecular interactions with the PIK3CA receptor binding domain in simulations, which is expected to decrease activity. On the other hand, overgrowth and cancer variants lead to loss of inhibitory interactions in simulations, which is expected to increase activity. We find that all disease variants display dysfunctions on either structural characteristics or intermolecular interaction energy. Thus, this comprehensive characterization of novel mosaic somatic variants associated with two opposing phenotypes has mechanistic importance and biomedical relevance and may aid in future therapeutic developments.

Abstract Pancreatic ductal adenocarcinoma (PDAC), has recently been found to be a heterogeneous disease, although the extension of its diversity remains to be fully understood. Here, we harmonize transcriptomic profiles derived from both PDAC epithelial and microenvironment cells to develop a Master Regulators (MR)-Gradient model that allows important inferences on transcriptional networks, epigenomic states, and metabolomics pathways that underlies this disease heterogeneity. This gradient model was generated by applying a blind source separation based on independent components analysis and robust principal component analyses (RPCA), following regulatory network inference. The result of these analyses reveals that PDAC prognosis strongly associates with the tumor epithelial cell phenotype and the immunological component. These studies were complemented by integration of methylome and metabolome datasets generated from patient-derived xenograft (PDX), together experimental measurements of metabolites, immunofluorescence microscopy, and western blot. At the metabolic level, PDAC favorable phenotype showed a positive correlation with enzymes implicated in complex lipid biosynthesis. In contrast, the unfavorable phenotype displayed an augmented OXPHOS independent metabolism centered on the Warburg effect and glutaminolysis. Epigenetically, we find that a global hypermethylation profile associates with the worst prognosis. Lastly, we report that, two antagonistic histone code writers, SUV39H1/SUV39H2 (H3K9Me3) and KAT2B (H3K9Ac) were identified key deregulated pathways in PDAC. Our analysis suggests that the PDAC phenotype, as it relates to prognosis, is determined by a complex interaction of transcriptomic, epigenomic, and metabolic features. Furthermore, we demonstrated that PDAC prognosis could be modulated through epigenetics.

NUPR1 binds to PARP1 and inhibits PARP1 activity in pancreatic cancer cells, with the NUPR1 inhibitor ZZW-115 abolishing the protective effect of NUPR1, inducing strong hyperPARylation and cell death by mitochondrial dysfunction.

Abstract Background Genomics enables individualized diagnosis and treatment, but large challenges remain to functionally interpret rare variants. To date, only one causative variant has been described for KCNK9 imprinting syndrome (KIS). The genotypic and phenotypic spectrum of KIS has yet to be described and the precise mechanism of disease fully understood. Methods This study discovers mechanisms underlying KCNK9 imprinting syndrome (KIS) by describing 15 novel KCNK9 alterations from 47 KIS-affected individuals. We use clinical genetics and computer-assisted facial phenotyping to describe the phenotypic spectrum of KIS. We then interrogate the functional effects of the variants in the encoded TASK3 channel using sequence-based analysis, 3D molecular mechanic and dynamic protein modeling, and in vitro electrophysiological and functional methodologies. Results We describe the broader genetic and phenotypic variability for KIS in a cohort of individuals identifying an additional mutational hotspot at p.Arg131 and demonstrating the common features of this neurodevelopmental disorder to include motor and speech delay, intellectual disability, early feeding difficulties, muscular hypotonia, behavioral abnormalities, and dysmorphic features. The computational protein modeling and in vitro electrophysiological studies discover variability of the impact of KCNK9 variants on TASK3 channel function identifying variants causing gain and others causing loss of conductance. The most consistent functional impact of KCNK9 genetic variants, however, was altered channel regulation. Conclusions This study extends our understanding of KIS mechanisms demonstrating its complex etiology including gain and loss of channel function and consistent loss of channel regulation. These data are rapidly applicable to diagnostic strategies, as KIS is not identifiable from clinical features alone and thus should be molecularly diagnosed. Furthermore, our data suggests unique therapeutic strategies may be needed to address the specific functional consequences of KCNK9 variation on channel function and regulation.

In the current study, we report computational scores for advancing genomic interpretation of disease-associated genomic variation in members of the RAS family of genes. For this purpose, we applied 31 sequence- and 3D structure-based computational scores, chosen by their breadth of biophysical properties. We parametrized our data by assembling a numerically homogenized experimentally-derived dataset, which when use in our calculations reveal that computational scores using 3D structure highly correlate with experimental measures (e.g., GAP-mediated hydrolysis RSpearman = 0.80 and RAF affinity Rspearman = 0.82), while sequence-based scores are discordant with this data. Performing all-against-all comparisons, we applied this parametrized modeling approach to the study of 935 RAS variants from 7 RAS genes, which led us to identify 4 groups of mutations according to distinct biochemical scores within each group. Each group was comprised of hotspot and non-hotspot KRAS variants, indicating that poorly characterized variants could functionally behave like pathogenic mutations. Combining computational scores using dimensionality reduction indicated that changes to local unfolding propensity associate with changes in enzyme activity by genomic variants. Hence, our systematic approach, combining methodologies from both clinical genomics and 3D structural bioinformatics, represents an expansion for interpreting genomic data, provides information of mechanistic value, and that is transferable to other proteins.

The histone demethylase KDM6A has recently elicited significant attention because its mutations are associated with a rare congenital disorder (Kabuki syndrome) and various types of human cancers. However, distinguishing KDM6A mutations that are deleterious to the enzyme and their underlying mechanisms of dysfunction remain to be fully understood. Here, we report the results from a multi-tiered approach evaluating the impact of 197 KDM6A somatic mutations using information derived from combining conventional genomics data with computational biophysics. This comprehensive approach incorporates multiple scores derived from alterations in protein sequence, structure, and molecular dynamics. Using this method, we classify the KDM6A mutations into 136 damaging variants (69.0%), 32 tolerated variants (16.2%), and 29 variants of uncertain significance (VUS, 14.7%), which is a significant improvement from the previous classification based on the conventional tools (over 40% VUS). We further classify the damaging variants into 15 structural variants (SV), 88 dynamic variants (DV), and 33 structural and dynamic variants (SDV). Comparison with variant scoring methods used in current clinical diagnosis guidelines demonstrates that our approach provides a more comprehensive evaluation of damaging potential and reveals mechanisms of dysfunction. Thus, these results should be taken into consideration for clinical assessment of the damaging potential of each mutation, as they provide hypotheses for experimental validation and critical information for the development of mutant-specific drugs to fight diseases caused by KDM6A dysfunctions.

Abstract Background Pancreatic ductal adenocarcinoma initiation is most frequently caused by Kras mutations. Results Here, we apply biological, biochemical, and network biology methods to validate GEMM-derived cell models using inducible KrasG12D expression. We describe the time-dependent, chromatin remodeling program that impacts function during early oncogenic signaling. We find that the KrasG12D-induced transcriptional response is dominated by downregulated expression concordant with layers of epigenetic events. More open chromatin characterizes the ATAC-seq profile associated with a smaller group of upregulated genes and epigenetic marks. RRBS demonstrates that promoter hypermethylation does not account for the silencing of the extensive gene promoter network. Moreover, ChIP-Seq reveals that heterochromatin reorganization plays little role in this early transcriptional program. Notably, both gene activation and silencing primarily depend on the marking of genes with a combination of H3K27ac, H3K4me3, and H3K36me3. Indeed, integrated modeling of all these datasets shows that KrasG12D regulates its transcriptional program primarily through unique super-enhancers and enhancers, and marking specific gene promoters and bodies. We also report chromatin remodeling across genomic areas that, although not contributing directly to cis-gene transcription, are likely important for KrasG12D functions. Conclusions In summary, we report a comprehensive, time-dependent, and coordinated early epigenomic program for KrasG12D in pancreatic cells, which is mechanistically relevant to understanding chromatin remodeling events underlying transcriptional outcomes needed for the function of this oncogene.

Abstract Background Kabuki syndrome is a genetic disorder that affects several body systems and presents with variations in symptoms and severity. The syndrome is named for a common phenotype of faces resembling stage makeup used in a Japanese traditional theatrical art named kabuki. The most frequent cause of this syndrome is mutations in the H3K4 family of histone methyltransferases while a smaller percentage results from genetic alterations affecting the histone demethylase, KDM6A. Because of the rare presentation of the latter form of the disease, little is known about how missense changes in the KDM6A protein sequence impact protein function. Results In this study, we use molecular mechanic and molecular dynamic simulations to enhance the annotation and mechanistic interpretation of the potential impact of eleven KDM6A missense variants found in Kabuki syndrome patients. These variants (N910S, D980V, S1025G, C1153R, C1153Y, P1195L, L1200F, Q1212R, Q1248R, R1255W, and R1351Q) are predicted to be pathogenic, likely pathogenic or of uncertain significance by sequence-based analysis. Here, we demonstrate, for the first time, that although Kabuki syndrome missense variants are found outside the functionally critical regions, they could affect overall function by significantly disrupting global and local conformation (C1153R, C1153Y, P1195L, L1200F, Q1212R, Q1248R, R1255W and R1351Q), chemical environment (C1153R, C1153Y, P1195L, L1200F, Q1212R, Q1248R, R1255W and R1351Q), and/or molecular dynamics of the catalytic domain (all variants). In addition, our approaches predict that many mutations, in particular C1153R, could allosterically disrupt the key enzymatic interactions of KDM6A. Conclusions Our study demonstrates that the KDM6A Kabuki syndrome variants may impair histone demethylase function through various mechanisms that include altered protein integrity, local environment, molecular interactions and protein dynamics. Molecular dynamics simulations of the wild type and the variants are critical to gain a better understanding of molecular dysfunction. This type of comprehensive structure- and MD-based analyses should help develop improved impact scoring systems to interpret the damaging effects of variants in this protein and other related proteins as well as provide detailed mechanistic insight that is not currently predictable from sequence alone.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive, painful disease with a 5-year survival rate of only 9%. Recent evidence indicates that distinct epigenomic landscapes underlie PDAC progression, identifying the H3K9me pathway as important to its pathobiology. Here, we delineate the role of Euchromatic Histone-lysine N-Methyltransferase 2 (EHMT2), the enzyme that generates H3K9me, as a downstream effector of oncogenic KRAS during PDAC initiation and pancreatitis-associated promotion. EHMT2 inactivation in pancreatic cells reduces H3K9me2 and antagonizes KrasG12D-mediated acinar-to-ductal metaplasia (ADM) and Pancreatic Intraepithelial Neoplasia (PanIN) formation in both the Pdx1-Cre and P48Cre/+KrasG12D mouse models. Ex vivo acinar explants also show impaired EGFR-KRAS-MAPK pathway-mediated ADM upon EHMT2 deletion. Notably, KrasG12D increases EHMT2 protein levels and EHMT2-EHMT1-WIZ complex formation. Transcriptome analysis reveals that EHMT2 inactivation upregulates a cell cycle inhibitory gene expression network that converges on the Cdkn1a/p21-Chek2 pathway. Congruently, pancreas tissue from KrasG12D animals with EHMT2 inactivation have increased P21 protein levels and enhanced senescence. Furthermore, loss of EHMT2 reduces inflammatory cell infiltration typically induced during KrasG12D-mediated initiation. The inhibitory effect on KrasG12D-induced growth is maintained in the pancreatitis-accelerated model, while simultaneously modifying immunoregulatory gene networks that also contribute to carcinogenesis. This study outlines the existence of a novel KRAS-EHMT2 pathway that is critical for mediating the growth-promoting and immunoregulatory effects of this oncogene in vivo, extending human observations to support a pathophysiological role for the H3K9me pathway in PDAC.

In pancreatic cancer, the epigenomic landscape can strongly impact the disease phenotype. Here, the authors discuss recent advances in our understanding of pancreatic cancer epigenomics, and how this knowledge can integrate with precision medicine approaches in this lethal disease.

Background & Aims: Transdifferentiation of hepatic stellate cells (HSCs) into myofibroblasts is a key event in the pathogenesis of liver fibrosis. Transforming growth factor β (TGF-β) and platelet-derived growth factor (PDGF) are canonical HSC activators after liver injury. The aim of this study was to analyze the epigenetic modulators that differentially control TGF-β and PDGF signaling pathways. Methods: We performed a transcriptomic comparison of HSCs treated with TGF-β or PDGF-BB using RNA sequencing. Among the targets that distinguish these 2 pathways, we focused on the histone methyltransferase class of epigenetic modulators. Results: Enhancer of zeste homolog 2 (EZH2) was expressed differentially, showing significant up-regulation in HSCs activated with TGF-β but not with PDGF-BB. Indeed, EZH2 inhibition using either a pharmacologic (GSK-503) or a genetic (small interfering RNA) approach caused a significant attenuation of TGF-β–induced fibronectin, collagen 1α1, and α-smooth muscle actin, both at messenger RNA and protein levels. Conversely, adenoviral overexpression of EZH2 in HSCs resulted in a significant stimulation of fibronectin protein and messenger RNA levels in TGF-β–treated cells. Finally, we conducted in vivo experiments with mice chronically treated with carbon tetrachloride or bile duct ligation. Administration of GSK-503 to mice receiving either carbon tetrachloride or bile duct ligation led to attenuated fibrosis as assessed by Trichrome and Sirius red stains, hydroxyproline, and α-smooth muscle actin/collagen protein assays. Conclusions: TGF-β and PDGF share redundant and distinct transcriptomic targets, with the former predominating in HSC activation. The EZH2 histone methyltransferase is preferentially involved in the TGF-β as opposed to the PDGF signaling pathway. Inhibition of EZH2 attenuates fibrogenic gene transcription in TGF-β–treated HSCs and reduces liver fibrosis in vivo. The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE119606 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE119606) Keywords: EZH2, Liver Fibrosis, Epigenetics, Histone Modifications

Background: Although significant advances have been made recently to characterize the biology of pancreatic ductal adenocarcinoma (PDAC), more efforts are needed to improve our understanding and to face challenges related to the aggressiveness, high mortality rate and chemoresistance of this disease. Methods: In this study, we perform the metabolomics profiling of 77 PDAC patient-derived tumor xenografts (PDTX) to investigate the relationship of metabolic profiles with overall survival (OS) in PDAC patients, tumor phenotypes and resistance to five anticancer drugs (gemcitabine, oxaliplatin, docetaxel, SN-38 and 5-Fluorouracil). Findings: We identified a metabolic signature that was able to predict the clinical outcome of PDAC patients (p < 0.001, HR=2.68 [95% CI: 1.5–4.9]). The correlation analysis showed that this metabolomic signature was significantly correlated with the PDAC molecular gradient (PAMG) (R = 0.44 and p < 0.001) indicating significant association to the transcriptomic phenotypes of tumors. Resistance score established, based on growth rate inhibition metrics using 35 PDTX-derived primary cells, allowed to identify several metabolites related to drug resistance which was globally accompanied by accumulation of several diacy-phospholipids and decrease in lysophospholipids. Interestingly, targeting glycerophospholipid synthesis improved sensitivity to the three tested cytotoxic drugs indicating that interfering with metabolism could be a promising therapeutic strategy to overcome the challenging resistance of PDAC. Interpretation: In conclusion, this study shows that the metabolomic profile of pancreatic PDTX models is strongly associated to clinical outcome, transcriptomic phenotypes and drug resistance. We also showed that targeting the lipidomic profile could be used in combinatory therapies against chemoresistance in PDAC.

We investigated germline variation in pancreatic ductal adenocarcinoma (PDAC) predisposition genes in 535 patients, using a custom-built panel and a new complementary bioinformatic approach. Our panel assessed genes belonging to DNA repair, cell cycle checkpoints, migration, and preneoplastic pancreatic conditions. Our bioinformatics approach integrated annotations of variants by using data derived from both germline and somatic references. This integrated approach with expanded evidence enabled us to consider patterns even among private mutations, supporting a functional role for certain alleles, which we believe enhances individualized medicine beyond classic gene-centric approaches. Concurrent evaluation of three levels of evidence, at the gene, sample, and cohort level, has not been previously done. Overall, we identified in PDAC patient germline samples, 12% with mutations previously observed in pancreatic cancers, 23% with mutations previously discovered by sequencing other human tumors, and 46% with mutations with germline associations to cancer. Non-polymorphic protein-coding pathogenic variants were found in 18.4% of patient samples. Moreover, among patients with metastatic PDAC, 16% carried at least one pathogenic variant, and this subgroup was found to have an improved overall survival (22.0 months versus 9.8; p=0.008) despite a higher pre-treatment CA19-9 level (p=0.02). Genetic alterations in DNA damage repair genes were associated with longer overall survival among patients who underwent resection surgery (92 months vs. 46; p=0.06). ATM alterations were associated with more frequent metastatic stage (p = 0.04) while patients with BRCA1 or BRCA2 alterations had improved overall survival (79 months vs. 39; p=0.05). We found that mutations in genes associated with chronic pancreatitis were more common in non-white patients (p

Gestational trophoblastic disease (GTD) is a heterogeneous group of lesions arising from placental tissue. Epithelioid trophoblastic tumor (ETT), derived from chorionic-type trophoblast, is the rarest form of GTD with only approximately 130 cases described in the literature. Due to its morphologic mimicry of epithelioid smooth muscle tumors and carcinoma, ETT can be misdiagnosed. To date, molecular characterization of ETTs is lacking. Furthermore, ETT is difficult to treat when disease spreads beyond the uterus. Here using RNA-Seq analysis in a cohort of ETTs and other gestational trophoblastic lesions we describe the discovery of LPCAT1-TERT fusion transcripts that occur in ETTs and coincide with underlying genomic deletions. Through cell-growth assays we demonstrate that LPCAT1-TERT fusion proteins can positively modulate cell proliferation and therefore may represent future treatment targets. Furthermore, we demonstrate that TERT upregulation appears to be a characteristic of ETTs, even in the absence of LPCAT1-TERT fusions, and that it appears linked to copy number gains of chromosome 5. No evidence of TERT upregulation was identified in other trophoblastic lesions tested, including placental site trophoblastic tumors and placental site nodules, which are thought to be the benign chorionic-type trophoblast counterpart to ETT. These findings indicate that LPCAT1-TERT fusions and copy-number driven TERT activation may represent novel markers for ETT, with the potential to improve the diagnosis, treatment, and outcome for women with this rare form of GTD.

Establishing the interactome of the cancer-associated stress protein Nuclear Protein 1 (NUPR1), we found that it binds to several hundreds of proteins, including proteins involved in nuclear translocation, DNA repair, and key factors of the SUMO pathway. We demonstrated that the NUPR1 inhibitor ZZW-115, an organic synthetic molecule, competes with importins for the binding to the NLS region of NUPR1, thereby inhibiting its nuclear translocation. We hypothesized, and then proved, that inhibition of NUPR1 by ZZW-115 sensitizes cancer cells to DNA damage induced by several genotoxic agents. Strikingly, we found that treatment with ZZW-115 reduced SUMOylation of several proteins involved in DNA damage response (DDR). We further report that the presence of recombinant NUPR1 improved the SUMOylation in a cell-free system, indicating that NUPR1 directly stimulates the SUMOylation machinery. We propose that ZZW-115 sensitizes cancer cells to genotoxic agents by inhibiting the nuclear translocation of NUPR1 and thereby decreasing the SUMOylation-dependent functions of key proteins involved in the DDR.

By reading the H3K9Me3 mark through their N-terminal chromodomain (CD), HP1 proteins play a significant role in cancer-associated processes, including cell proliferation, differentiation, chromosomal stability, and DNA repair. Here, we used a combination of bioinformatics-based methodologies, as well as experimentally-derived datasets, that reveal the existence of a novel short HP1γ (CBX3) isoform, named here sHP1γ, generated by alternative splicing of the CBX3 locus. The sHP1γ mRNA encodes a protein composed of 101 residues and lacks the C-terminal chromoshadow domain (CSD) that is required for dimerization and heterodimerization in the previously described 183 a. a HP1γ protein. Fold recognition, order-to-disorder calculations, threading, homology-based molecular modeling, docking, and molecular dynamic simulations show that the sHP1γ is comprised of a CD flanked by intrinsically disordered regions (IDRs) with an IDR-CD-IDR domain organization and likely retains the ability to bind to the H3K9Me3. Both qPCR analyses and mRNA-seq data derived from large-scale studies confirmed that sHP1γ mRNA is expressed in the majority of human tissues at approximately constant ratios with the chromoshadow domain containing isoform. However, sHP1γ mRNA levels appear to be dysregulated in different cancer types. Thus, our data supports the notion that, due to the existence of functionally different isoforms, the regulation of HP1γ-mediated functions is more complex than previously anticipated.

Pancreatic ductal adenocarcinoma (PDAC) is a complex disease and its underlying epigenomic heterogeneity is not fully understood. Here, the authors utilize patient-derived PDAC xenografts to define the epigenomic landscape of PDAC, highlighting chromatin states linked to differing disease aggressiveness and survival.

International experts gathered at the Mayo Clinic (Rochester MN, USA) on February 27th-28th, 2017 for a meeting entitled ‘Basic and Translational Facets of the Epigenetics of GI Diseases’. This workshop summarized recent advances on the role of epigenetics in the pathobiology of gastrointestinal (GI) diseases. Highlights of the meeting included recent advances on the involvement of different epigenetic mechanisms in malignant and nonmalignant GI disorders and the epigenetic heterogeneity exhibited in these diseases. The translational value of epigenetic drugs, as well as the current and future use of epigenetic changes (i.e., DNA methylation patterns) as biomarkers for early detection tools or disease stratification were also important topics of discussion.

Evidence-based medicine (EBM) relies on scientific data to guide diagnosis and treatment and is recognized as the current paradigm in medicine. Accordingly, every current and future physician should be knowledgeable about its principles and methodologies. Embracing research and EBM is essential to modern clinical practice; however, trainees and physicians still struggle with the value of research-related courses and knowledge on epidemiology and research methodology is often poor. In this article we provide a cogent discussion of the importance of research as an indispensable discipline in medical education through a detailed analysis of the literature. We review the evolution of medicine towards EBM and discuss the myriad of benefits that research has on medical careers, leadership roles, mentoring relationships, social networking, and personal growth and development. Participation in research contributes to medicine, public health, and society while simultaneously allowing the achievement of a high level of personal satisfaction

An amendment to this paper has been published and can be accessed via the original article.

Preclinical models based on patient-derived xenografts have remarkable specificity in distinguishing transformed human tumor cells from non-transformed murine stromal cells computationally. We obtained 29 pancreatic ductal adenocarcinoma (PDAC) xenografts from either resectable or non-resectable patients (surgery and endoscopic ultrasound-guided fine-needle aspirate, respectively). Extensive multiomic profiling revealed two subtypes with distinct clinical outcomes. These subtypes uncovered specific alterations in DNA methylation and transcription as well as in signaling pathways involved in tumor-stromal cross-talk. The analysis of these pathways indicates therapeutic opportunities for targeting both compartments and their interactions. In particular, we show that inhibiting NPC1L1 with Ezetimibe, a clinically available drug, might be an efficient approach for treating pancreatic cancers. These findings uncover the complex and diverse interplay between PDAC tumors and the stroma and demonstrate the pivotal role of xenografts for drug discovery and relevance to PDAC.

Abstract c‐MYC controls more than 15% of genes responsible for proliferation, differentiation, and cellular metabolism in pancreatic as well as other cancers making this transcription factor a prime target for treating patients. The transcriptome of 55 patient‐derived xenografts show that 30% of them share an exacerbated expression profile of MYC transcriptional targets (MYC‐high). This cohort is characterized by a high level of Ki67 staining, a lower differentiation state, and a shorter survival time compared to the MYC‐low subgroup. To define classifier expression signature, we selected a group of 10 MYC target transcripts which expression is increased in the MYC‐high group and six transcripts increased in the MYC‐low group. We validated the ability of these markers panel to identify MYC‐high patient‐derived xenografts from both: discovery and validation cohorts as well as primary cell cultures from the same patients. We then showed that cells from MYC‐high patients are more sensitive to JQ1 treatment compared to MYC‐low cells, in monolayer, 3D cultured spheroids and in vivo xenografted tumors, due to cell cycle arrest followed by apoptosis. Therefore, these results provide new markers and potentially novel therapeutic modalities for distinct subgroups of pancreatic tumors and may find application to the future management of these patients within the setting of individualized medicine clinics.

Graphical Abstract

Purpose: Demand is increasing for clinical genomic sequencing to provide diagnoses for patients presenting phenotypes indicative of genetic diseases, but for whom routine genetic testing failed to yield a diagnosis. DNA-based testing using high-throughput technologies often identifies variants with insufficient evidence to determine whether they are disease-causal or benign, leading to categorization as variants of uncertain significance (VUS).Methods: We used molecular modeling and simulation to generate specific hypotheses for the molecular effects of variants in the human glucose transporter, GLUT10 (SLC2A10). Similar to many disease-relevant membrane proteins, no experimentally derived 3D structure exists. An atomic model was generated and used to evaluate multiple variants, including pathogenic, benign, and VUS.Results: These analyses yielded detailed mechanistic data, not currently predictable from sequence, including altered protein stability, charge distribution of ligand binding surfaces, and shifts toward or away from transport-competent conformations. Consideration of the two major conformations of GLUT10 was important as variants have conformation-specific effects. We generated detailed molecular hypotheses for the functional impact of variants in GLUT10 and propose means to determine their pathogenicity.Conclusion: The type of workflow we present here is valuable for increasing the throughput and resolution with which VUS effects can be assessed and interpreted.

Variants in the TGFBR2 kinase domain cause several human diseases and can increase propensity for cancer. The widespread application of next generation sequencing within the setting of Individualized Medicine (IM) is increasing the rate at which TGFBR2 kinase domain variants are being identified. However, their clinical relevance is often uncertain. Consequently, we sought to evaluate the use of molecular modeling and molecular dynamics (MD) simulations for assessing the potential impact of variants within this domain. We documented the structural differences revealed by these models across 57 variants using independent MD simulations for each. Our simulations revealed various mechanisms by which variants may lead to functional alteration; some are revealed energetically, while others structurally or dynamically. We found that the ATP binding site and activation loop dynamics may be affected by variants at positions throughout the structure. This prediction cannot be made from the linear sequence alone. We present our structure-based analyses alongside those obtained using several commonly used genomics-based predictive algorithms. We believe the further mechanistic information revealed by molecular modeling will be useful in guiding the examination of clinically observed variants throughout the exome, as well as those likely to be discovered in the near future by clinical tests leveraging next-generation sequencing through IM efforts.

Cancer cells activate stress-response mechanisms to adapt themselves to a variety of stressful conditions. Among these protective mechanisms, those controlled by the stress-induced nuclear protein 1 (NUPR1) belong to the most conserved ones. NUPR1 is an 82-residue-long, monomeric, basic and intrinsically disordered protein (IDP), which was found to be invariably overexpressed in some, if not all, cancer tissues. Remarkably, we and others have previously showed that genetic inactivation of the Nupr1 gene antagonizes the growth of pancreatic cancer as well as several other tumors. With the use of a multidisciplinary strategy by combining biophysical, biochemical, bioinformatic, and biological approaches, a trifluoperazine-derived compound, named ZZW-115, has been identified as an inhibitor of the NUPR1 functions. The anticancer activity of the ZZW-115 was first validated on a large panel of cancer cells. Furthermore, ZZW-115 produced a dose-dependent tumor regression of the tumor size in xenografted mice. Mechanistically, we have demonstrated that NUPR1 binds to several importins. Because ZZW-115 binds NUPR1 through the region around the amino acid Thr68, which is located into the nuclear location signal (NLS) region of the protein, we demonstrated that treatment with ZZW-115 inhibits completely the translocation of NUPR1 from the cytoplasm to the nucleus by competing with importins.

During the last decades it has been established that breast cancer arises through the accumulation of genetic and epigenetic alterations in different cancer related genes. These alterations confer the tumor oncogenic abilities, which can be resumed as cancer hallmarks (CH). The purpose of this study was to establish the methylation profile of CpG sites located in cancer genes in breast tumors so as to infer their potential impact on 6 CH: i.e. sustained proliferative signaling, evasion of growth suppressors, resistance to cell death, induction of angiogenesis, genome instability and invasion and metastasis. For 51 breast carcinomas, MS-MLPA derived-methylation profiles of 81 CpG sites were converted into 6 CH profiles. CH profiles distribution was tested by different statistical methods and correlated with clinical-pathological data. Unsupervised Hierarchical Cluster Analysis revealed that CH profiles segregate in two main groups (bootstrapping 90-100%), which correlate with breast laterality (p = 0.05). For validating these observations, gene expression data was obtained by RealTime-PCR in a different cohort of 25 tumors and converted into CH profiles. This analyses confirmed the same clustering and a tendency of association with breast laterality (p = 0.15). In silico analyses on gene expression data from TCGA Breast dataset from left and right breast tumors showed that they differed significantly when data was previously converted into CH profiles (p = 0.033). We show here for the first time, that breast carcinomas arising on different sides of the body present differential cancer traits inferred from methylation and expression profiles. Our results indicate that by converting methylation or expression profiles in terms of Cancer Hallmarks, it would allow to uncover veiled associations with clinical features. These results contribute with a new finding to the better understanding of breast tumor behavior, and can moreover serve as proof of principle for other bilateral cancers like lung, testes or kidney.

Developmental morphogens play an important role in coordinating the ductular reaction and portal fibrosis occurring in the setting of cholangiopathies. However, little is known about how membrane signaling events in ductular reactive cells (DRCs) are transduced into nuclear transcriptional changes to drive cholangiocyte maturation and matrix deposition. Therefore, the aim of this study was to investigate potential mechanistic links between cell signaling events and epigenetic regulators in DRCs.Using directed differentiation of induced pluripotent stem cells (iPSC), isolated DRCs, and in vivo models, we examine the mechanisms whereby sonic hedgehog (Shh) overcomes an epigenetic barrier in biliary precursors and promotes both cholangiocyte maturation and deposition of fibronectin (FN).We demonstrate, for the first time, that Gli1 influences the differentiation state and fibrogenic capacity of iPSC-derived hepatic progenitors and isolated DRCs. We outline a novel pathway wherein Shh-mediated Gli1 binding in key cholangiocyte gene promoters overcomes an epigenetic barrier conferred by the polycomb protein, enhancer of zeste homolog 2 (EZH2) and initiates the transcriptional program of cholangiocyte maturation. We also define previously unknown functional Gli1 binding sites in the promoters of cytokeratin (CK)7, CK19, and FN. Our in vivo results show that EZH2 KO mice fed the choline-deficient, ethanolamine supplemented (CDE) diet have an exaggerated cholangiocyte expansion associated with more robust ductular reaction and increased peri-portal fibrosis.We conclude that Shh/Gli1 signaling plays an integral role in cholangiocyte maturation in vitro by overcoming an EZH2-dependent epigenetic barrier and this mechanism also promotes biliary expansion in vivo.

Abstract Pancreatic adenocarcinoma (PDAC) is an extremely deadly disease for which all treatments available have failed to improve life expectancy significantly. This may be explained by the high metastatic potential of PDAC cells, which results from their dedifferentiation towards a mesenchymal phenotype. Some PDAC present cell‐in‐cell structures whose origin and significance are currently unknown. We show here that cell‐in‐cells form after homotypic cell cannibalism (HoCC). We found PDAC patients whose tumours display HoCC develop less metastasis than those without. In vitro, HoCC was promoted by inactivation of the nuclear protein 1 (Nupr1), and was enhanced by treatment with transforming growth factor β. HoCC ends with death of PDAC cells, consistent with a metastasis suppressor role for this phenomenon. Hence, our data indicates a protective role for HoCC in PDAC and identifies Nupr1 as a molecular regulator of this process.

Krüppel-like factor 11 (KLF11) and the highly homologous KLF10 proteins are transcription factors originating from duplication of the Drosophila melanogaster ancestor cabut. The function of these proteins in epithelial cells has been previously characterized. In the current study, we report a functional role for KLF11 in mesenchymal cells and in mesenchymal cell dysfunction, namely, fibrosis, and subsequently perform a detailed cellular, molecular, and in vivo characterization of this phenomenon. We find that, in cultured mesenchymal cells, enhanced expression of KLF11 results in activated extracellular matrix pathways, including collagen gene silencing and matrix metalloproteinases activation without changes in tissue inhibitors of metalloproteinases. Combined, reporter and chromatin immunoprecipitation assays demonstrate that KLF11 interacts directly with the collagen 1a2 (COL1A2) promoter in mesenchymal cells to repress its activity. Mechanistically, KLF11 regulates collagen gene expression through the heterochromatin protein 1 gene-silencing pathway as mutants defective for coupling to this epigenetic modifier lose the ability to repress COL1A2. Expression studies reveal decreased levels of KLF11 during liver fibrogenesis after chemically induced injury in vivo. Congruently, KLF11(-/-) mice, which should be deficient in the hypothesized anti-fibrogenic brake imposed by this transcription factor, display an enhanced response to liver injury with increased collagen fibril deposition. Thus, KLFs expands the repertoire of transcription factors involved in the regulation of extracellular matrix proteins in mesenchymal cells and define a novel pathway that modulates the fibrogenic response during liver injury.

Endometriosis affects approximately 10% of young, reproductive-aged women. Disease associated pelvic pain; infertility and sexual dysfunction have a significant adverse clinical, social and financial impact. As precise disease etiology has remained elusive, current therapeutic strategies are empiric, unfocused and often unsatisfactory. Lack of a suitable genetic model has impaired further translational research in the field. In this study, we evaluated the role of the Sp/KLF transcription factor KLF11/Klf11 in the pathogenesis of endometriosis. KLF11, a human disease-associated gene is etiologically implicated in diabetes, uterine fibroids and cancer. We found that KLF11 expression was diminished in human endometriosis implants and further investigated its pathogenic role in Klf11-/- knockout mice with surgically induced endometriotic lesions. Lesions in Klf11-/- animals were large and associated with prolific fibrotic adhesions resembling advanced human disease in contrast to wildtype controls. To determine phenotype-specificity, endometriosis was also generated in Klf9-/- animals. Unlike in Klf11-/- mice, lesions in Klf9-/- animals were neither large, nor associated with a significant fibrotic response. KLF11 also bound to specific elements located in the promoter regions of key fibrosis-related genes from the Collagen, MMP and TGF-β families in endometrial stromal cells. KLF11 binding resulted in transcriptional repression of these genes. In summary, we identify a novel pathogenic role for KLF11 in preventing de novo disease-associated fibrosis in endometriosis. Our model validates in vivo the phenotypic consequences of dysregulated Klf11 signaling. Additionally, it provides a robust means not only for further detailed mechanistic investigation but also the ability to test any emergent translational ramifications thereof, so as to expand the scope and capability for treatment of endometriosis.

Class I MAGE proteins (MAGE I) are normally expressed only in developing germ cells but are aberrantly expressed in many cancers. They have been shown to promote tumor survival, aggressive growth, and chemoresistance but the underlying mechanisms and MAGE I functions have not been fully elucidated. KRAB domain zinc finger transcription factors (KZNFs) are the largest group of vertebrate transcription factors and regulate neoplastic transformation, tumor suppression, cellular proliferation, and apoptosis. KZNFs bind the KAP1 protein and direct KAP1 to specific DNA sequences where it suppresses gene expression by inducing localized heterochromatin characterized by histone 3 lysine 9 trimethylation (H3me3K9). Discovery that MAGE I proteins also bind to KAP1 prompted us to investigate whether MAGE I can affect KZNF and KAP1 mediated gene regulation. We found that expression of MAGE I proteins, MAGE-A3 or MAGE-C2, relieved repression of a reporter gene by ZNF382, a KZNF with tumor suppressor activity. ChIP of MAGE I (-) HEK293T cells showed KAP1 and H3me3K9 are normally bound to the ID1 gene, a target of ZNF382, but that binding is greatly reduced in the presence of MAGE I proteins. MAGE I expression relieved KAP1 mediated ID1 repression, causing increased expression of ID1 mRNA and ID1 chromatin relaxation characterized by loss of H3me3K9. MAGE I binding to KAP1 also induced ZNF382 poly-ubiquitination and degradation, consistent with loss of ZNF382 leading to decreased KAP1 binding to ID1. In contrast, MAGE I expression caused increased KAP1 binding to Ki67, another KAP1 target gene, with increased H3me3K9 and decreased Ki67 mRNA expression. Since KZNFs are required to direct KAP1 to specific genes, these results show that MAGE I proteins can differentially regulate members of the KZNF family and KAP1 mediated gene repression.