"Stripe" transcription factors provide accessibility to co-binding partners in mammalian genomes.

Regulatory elements activate promoters by recruiting transcription factors (TFs) to specific motifs. Notably, TF-DNA interactions often depend on cooperativity with colocalized partners, suggesting an underlying cis-regulatory syntax. To explore TF cooperativity in mammals, we analyze ∼500 mouse and human primary cells by combining an atlas of TF motifs, footprints, ChIP-seq, transcriptomes, and accessibility. We uncover two TF groups that colocalize with most expressed factors, forming stripes in hierarchical clustering maps. The first group includes lineage-determining factors that occupy DNA elements broadly, consistent with their key role in tissue-specific transcription. The second one, dubbed universal stripe factors (USFs), comprises ∼30 SP, KLF, EGR, and ZBTB family members that recognize overlapping GC-rich sequences in all tissues analyzed. Knockouts and single-molecule tracking reveal that USFs impart accessibility to colocalized partners and increase their residence time. Mammalian cells have thus evolved a TF superfamily with overlapping DNA binding that facilitate chromatin accessibility.

Cancer-associated polybromo-1 bromodomain 4 missense variants variably impact bromodomain ligand binding and cell growth suppression.

The polybromo, brahma-related gene 1-associated factors (PBAF) chromatin remodeling complex subunit polybromo-1 (PBRM1) contains six bromodomains that recognize and bind acetylated lysine residues on histone tails and other nuclear proteins. PBRM1 bromodomains thus provide a link between epigenetic posttranslational modifications and PBAF modulation of chromatin accessibility and transcription. As a putative tumor suppressor in several cancers, PBRM1 protein expression is often abrogated by truncations and deletions. However, ∼33% of PBRM1 mutations in cancer are missense and cluster within its bromodomains. Such mutations may generate full-length PBRM1 variant proteins with undetermined structural and functional characteristics. Here, we employed computational, biophysical, and cellular assays to interrogate the effects of PBRM1 bromodomain missense variants on bromodomain stability and function. Since mutations in the fourth bromodomain of PBRM1 (PBRM1-BD4) comprise nearly 20% of all cancer-associated PBRM1 missense mutations, we focused our analysis on PBRM1-BD4 missense protein variants. Selecting 16 potentially deleterious PBRM1-BD4 missense protein variants for further study based on high residue mutational frequency and/or conservation, we show that cancer-associated PBRM1-BD4 missense variants exhibit varied bromodomain stability and ability to bind acetylated histones. Our results demonstrate the effectiveness of identifying the unique impacts of individual PBRM1-BD4 missense variants on protein structure and function, based on affected residue location within the bromodomain. This knowledge provides a foundation for drawing correlations between specific cancer-associated PBRM1 missense variants and distinct alterations in PBRM1 function, informing future cancer personalized medicine approaches.

Advanced computational analysis of CD40LG variants in atypical X-linked hyper-IgM syndrome.

X-Linked Hyper-IgM Syndrome is caused by pathogenic variants in CD40LG. Three patients with atypical clinical and immunological features were identified with variants in CD40LG requiring further characterization. Flow cytometry was used to evaluate CD40L protein expression and binding capacity to a surrogate receptor, CD40-muIg. Though functional anomalies were observed, there was still a lack of clarity regarding the underlying mechanism. We developed structural models for wild-type and the three variants of CD40L protein observed in these patients (p. Lys143Asn, Leu225Ser and Met36Arg) to evaluate structural alterations by molecular mechanic calculations, and assess protein movement by molecular dynamic simulations. These studies demonstrate that functional analysis of variants of unknown significance in CD40LG can be supplemented by advanced computational analysis in atypical clinical contexts. These studies in combination identify the deleterious effects of these variants and potential mechanisms for protein dysfunction.

Effects of Subnormothermic Regulated Hepatic Reperfusion on Mitochondrial and Transcriptomic Profiles in a Porcine Model.

OBJECTIVE: We sought to investigate the biological effects of pre-reperfusion treatments of the liver after warm and cold ischemic injuries in a porcine donation after circulatory death model. SUMMARY OF BACKGROUND DATA: Donation after circulatory death represents a severe form of liver ischemia and reperfusion injury that has a profound impact on graft function after liver transplantation. METHODS: Twenty donor pig livers underwent 60 minutes of in situ warm ischemia after circulatory arrest and 120 minutes of cold static preservation prior to simulated transplantation using an ex vivo perfusion machine. Four reperfusion treatments were compared: Control-Normothermic (N), Control- Subnormothermic (S), regulated hepatic reperfusion (RHR)-N, and RHR-S (n = 5 each). The biochemical, metabolic, and transcriptomic profiles, as well as mitochondrial function were analyzed. RESULTS: Compared to the other groups, RHR-S treated group showed significantly lower post-reperfusion aspartate aminotransferase levels in the reperfusion effluent and histologic findings of hepatocyte viability and lesser degree of congestion and necrosis. RHR-S resulted in a significantly higher mitochondrial respiratory control index and calcium retention capacity. Transcriptomic profile analysis showed that treatment with RHR-S activated cell survival and viability, cellular homeostasis as well as other biological functions involved in tissue repair such as cytoskeleton or cytoplasm organization, cell migration, transcription, and microtubule dynamics. Furthermore, RHR-S inhibited organismal death, morbidity and mortality, necrosis, and apoptosis. CONCLUSION: Subnormothermic RHR mitigates IRI and preserves hepatic mitochondrial function after warm and cold hepatic ischemia. This organ resuscitative therapy may also trigger the activation of protective genes against IRI. Sub- normothermic RHR has potential applicability to clinical liver transplantation.

SASH3 variants cause a novel form of X-linked combined immunodeficiency with immune dysregulation.

Sterile alpha motif (SAM) and Src homology-3 (SH3) domain-containing 3 (SASH3), also called SH3-containing lymphocyte protein (SLY1), is a putative adaptor protein that is postulated to play an important role in the organization of signaling complexes and propagation of signal transduction cascades in lymphocytes. The SASH3 gene is located on the X-chromosome. Here, we identified 3 novel SASH3 deleterious variants in 4 unrelated male patients with a history of combined immunodeficiency and immune dysregulation that manifested as recurrent sinopulmonary, cutaneous, and mucosal infections and refractory autoimmune cytopenias. Patients exhibited CD4+ T-cell lymphopenia, decreased T-cell proliferation, cell cycle progression, and increased T-cell apoptosis in response to mitogens. In vitro T-cell differentiation of CD34+ cells and molecular signatures of rearrangements at the T-cell receptor α (TRA) locus were indicative of impaired thymocyte survival. These patients also manifested neutropenia and B-cell and natural killer (NK)-cell lymphopenia. Lentivirus-mediated transfer of the SASH3 complementary DNA-corrected protein expression, in vitro proliferation, and signaling in SASH3-deficient Jurkat and patient-derived T cells. These findings define a new type of X-linked combined immunodeficiency in humans that recapitulates many of the abnormalities reported in mice with Sly1-/- and Sly1Δ/Δ mutations, highlighting an important role of SASH3 in human lymphocyte function and survival.

Browning of human adipocytes requires KLF11 and reprogramming of PPARγ superenhancers.

Long-term exposure to peroxisome proliferator-activated receptor γ (PPARγ) agonists such as rosiglitazone induces browning of rodent and human adipocytes; however, the transcriptional mechanisms governing this phenotypic switch in adipocytes are largely unknown. Here we show that rosiglitazone-induced browning of human adipocytes activates a comprehensive gene program that leads to increased mitochondrial oxidative capacity. Once induced, this gene program and oxidative capacity are maintained independently of rosiglitazone, suggesting that additional browning factors are activated. Browning triggers reprogramming of PPARγ binding, leading to the formation of PPARγ "superenhancers" that are selective for brown-in-white (brite) adipocytes. These are highly associated with key brite-selective genes. Based on such an association, we identified an evolutionarily conserved metabolic regulator, Kruppel-like factor 11 (KLF11), as a novel browning transcription factor in human adipocytes that is required for rosiglitazone-induced browning, including the increase in mitochondrial oxidative capacity. KLF11 is directly induced by PPARγ and appears to cooperate with PPARγ in a feed-forward manner to activate and maintain the brite-selective gene program.

A GATA6-centred gene regulatory network involving HNFs and ΔNp63 controls plasticity and immune escape in pancreatic cancer.

OBJECTIVE: Molecular taxonomy of tumours is the foundation of personalised medicine and is becoming of paramount importance for therapeutic purposes. Four transcriptomics-based classification systems of pancreatic ductal adenocarcinoma (PDAC) exist, which consistently identified a subtype of highly aggressive PDACs with basal-like features, including ΔNp63 expression and loss of the epithelial master regulator GATA6. We investigated the precise molecular events driving PDAC progression and the emergence of the basal programme. DESIGN: We combined the analysis of patient-derived transcriptomics datasets and tissue samples with mechanistic experiments using a novel dual-recombinase mouse model for Gata6 deletion at late stages of KRasG12D-driven pancreatic tumorigenesis (Gata6LateKO). RESULTS: This comprehensive human-to-mouse approach showed that GATA6 loss is necessary, but not sufficient, for the expression of ΔNp63 and the basal programme in patients and in mice. The concomitant loss of HNF1A and HNF4A, likely through epigenetic silencing, is required for the full phenotype switch. Moreover, Gata6 deletion in mice dramatically increased the metastatic rate, with a propensity for lung metastases. Through RNA-Seq analysis of primary cells isolated from mouse tumours, we show that Gata6 inhibits tumour cell plasticity and immune evasion, consistent with patient-derived data, suggesting that GATA6 works as a barrier for acquiring the fully developed basal and metastatic phenotype. CONCLUSIONS: Our work provides both a mechanistic molecular link between the basal phenotype and metastasis and a valuable preclinical tool to investigate the most aggressive subtype of PDAC. These data, therefore, are important for understanding the pathobiological features underlying the heterogeneity of pancreatic cancer in both mice and human.

Computational modeling reveals key molecular properties and dynamic behavior of disruptor of telomeric silencing 1-like (DOT1L) and partnering complexes involved in leukemogenesis.

Disruptor of telomeric silencing 1-like (DOT1L) is the only non-SET domain histone lysine methyltransferase (KMT) and writer of H3K79 methylation on nucleosomes marked by H2B ubiquitination. DOT1L has elicited significant attention because of its interaction or fusion with members of the AF protein family in blood cell biology and leukemogenic transformation. Here, our goal was to extend previous structural information by performing a robust molecular dynamic study of DOT1L and its leukemogenic partners combined with mutational analysis. We show that statically and dynamically, D161, G163, E186, and F223 make frequent time-dependent interactions with SAM, while additional residues T139, K187, and N241 interact with SAM only under dynamics. Dynamics models reveal DOT1L, SAM, and H4 moving as one and show that more than twice the number of DOT1L residues interacts with these partners, relative to the static structure. Mutational analyses indicate that six of these residues are intolerant to substitution. We describe the dynamic behavior of DOT1L interacting with AF10 and AF9. Studies on the dynamics of a heterotrimeric complex of DOT1L1-AF10 illuminated describe coordinated motions that impact the relative position of the DOT1L HMT domain to the nucleosome. The molecular motions of the DOT1L-AF9 complex are less extensive and highly dynamic, resembling a swivel-like mechanics. Through molecular dynamics and mutational analysis, we extend the knowledge previous provided by static measurements. These results are important to consider when describing the biochemical properties of DOT1L, under normal and in disease conditions, as well as for the development of novel therapeutic agents.

Structural bioinformatics enhances mechanistic interpretation of genomic variation, demonstrated through the analyses of 935 distinct RAS family mutations.

MOTIVATION: Protein-coding genetic alterations are frequently observed in Clinical Genetics, but the high yield of variants of uncertain significance remains a limitation in decision making. RAS-family GTPases are cancer drivers, but only 54 variants, across all family members, fall within well-known hotspots. However, extensive sequencing has identified 881 non-hotspot variants for which significance remains to be investigated. RESULTS: Here, we evaluate 935 missense variants from seven RAS genes, observed in cancer, RASopathies and the healthy adult population. We characterized hotspot variants, previously studied experimentally, using 63 sequence- and 3D structure-based scores, chosen by their breadth of biophysical properties. Applying scores that display best correlation with experimental measures, we report new valuable mechanistic inferences for both hot-spot and non-hotspot variants. Moreover, we demonstrate that 3D scores have little-to-no correlation with those based on DNA sequence, which are commonly used in Clinical Genetics. Thus, combined, these new knowledge bear significant relevance. AVAILABILITY AND IMPLEMENTATION: All genomic and 3D scores, and markdown for generating figures, are provided in our supplemental data. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SOD2 acetylation on lysine 68 promotes stem cell reprogramming in breast cancer.

Mitochondrial superoxide dismutase (SOD2) suppresses tumor initiation but promotes invasion and dissemination of tumor cells at later stages of the disease. The mechanism of this functional switch remains poorly defined. Our results indicate that as SOD2 expression increases acetylation of lysine 68 ensues. Acetylated SOD2 promotes hypoxic signaling via increased mitochondrial reactive oxygen species (mtROS). mtROS, in turn, stabilize hypoxia-induced factor 2α (HIF2α), a transcription factor upstream of "stemness" genes such as Oct4, Sox2, and Nanog. In this sense, our findings indicate that SOD2K68Ac and mtROS are linked to stemness reprogramming in breast cancer cells via HIF2α signaling. Based on these findings we propose that, as tumors evolve, the accumulation of SOD2K68Ac turns on a mitochondrial pathway to stemness that depends on HIF2α and may be relevant for the progression of breast cancer toward poor outcomes.

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.

Inorganic arsenic promotes luminal to basal transition and metastasis of breast cancer.

Inorganic arsenic (iAs/As2 O32- ) is an environmental toxicant found in watersheds around the world including in densely populated areas. iAs is a class I carcinogen known to target the skin, lungs, bladder, and digestive organs, but its role as a primary breast carcinogen remains controversial. Here, we examined a different possibility: that exposure to iAs promotes the transition of well-differentiated epithelial breast cancer cells characterized by estrogen and progesterone receptor expression (ER+/PR+), to more basal phenotypes characterized by active proliferation, and propensity to metastasis in vivo. Our results indicate two clear phenotypic responses to low-level iAs that depend on the duration of the exposure. Short-term pulses of iAs activate ER signaling, consistent with its reported pseudo-estrogen activity, but longer-term, chronic treatments for over 6 months suppresses both ER and PR expression and signaling. In fact, washout of these chronically exposed cells for up to 1 month failed to fully reverse the transcriptional and phenotypic effects of prolonged treatments, indicating durable changes in cellular physiologic identity. RNA-seq studies found that chronic iAs drives the transition toward more basal phenotypes characterized by impaired hormone receptor signaling despite the conservation of estrogen receptor expression. Because treatments for breast cancer patients are largely designed based on the detection of hormone receptor expression, our results suggest greater scrutiny of ER+ cancers in patients exposed to iAs, because these tumors may spawn more aggressive phenotypes than unexposed ER+ tumors, in particular, basal subtypes that tend to develop therapy resistance and metastasis.

ErbB3-binding protein 1 (EBP1) represses HNF4α-mediated transcription and insulin secretion in pancreatic ß-cells.

HNF4α (hepatocyte nuclear factor 4α) is one of the master regulators of pancreatic ß-cell development and function, and mutations in the HNF4α gene are well-known monogenic causes of diabetes. As a member of the nuclear receptor family, HNF4α exerts its gene regulatory function through various molecular interactions; however, there is a paucity of knowledge of the different functional complexes in which HNF4α participates. Here, to find HNF4α-binding proteins in pancreatic ß-cells, we used yeast two-hybrid screening, a mammalian two-hybrid assay, and glutathione S-transferase pulldown approaches, which identified EBP1 (ErbB3-binding protein 1) as a factor that binds HNF4α in a LXXLL motif-mediated manner. In the ß-cells, EBP1 suppressed the expression of HNF4α target genes that are implicated in insulin secretion, which is impaired in HNF4α mutation-driven diabetes. The crystal structure of the HNF4α ligand-binding domain in complex with a peptide harboring the EBP1 LXXLL motif at 3.15Å resolution hinted at the molecular basis of the repression. The details of the structure suggested that EBP1's LXXLL motif competes with HNF4α coactivators for the same binding pocket and thereby prevents recruitment of additional transcriptional coactivators. These findings provide further evidence that EBP1 plays multiple cellular roles and is involved in nuclear receptor-mediated gene regulation. Selective disruption of the HNF4α-EBP1 interaction or tissue-specific EBP1 inactivation can enhance HNF4α activities and thereby improve insulin secretion in ß-cells, potentially representing a new strategy for managing diabetes and related metabolic disorders.

Intrinsically disordered chromatin protein NUPR1 binds to the C-terminal region of Polycomb RING1B.

Intrinsically disordered proteins (IDPs) are ubiquitous in eukaryotes, and they are often associated with diseases in humans. The protein NUPR1 is a multifunctional IDP involved in chromatin remodeling and in the development and progression of pancreatic cancer; however, the details of such functions are unknown. Polycomb proteins are involved in specific transcriptional cascades and gene silencing. One of the proteins of the Polycomb complex is the Ring finger protein 1 (RING1). RING1 is related to aggressive tumor features in multiple cancer types. In this work we characterized the interaction between NUPR1 and the paralogue RING1B in vitro, in silico, and in cellulo. The interaction occurred through the C-terminal region of RING1B (C-RING1B), with an affinity in the low micromolar range (∼10 µM). The binding region of NUPR1, mapped by NMR, was a hydrophobic polypeptide patch at the 30s region of its sequence, as pinpointed by computational results and site-directed mutagenesis at Ala33. The association between C-RING1B and wild-type NUPR1 also occurred in cellulo as tested by protein ligation assays; this interaction is inhibited by trifluoperazine, a drug known to hamper binding of wild-type NUPR1 with other proteins. Furthermore, the Thr68Gln and Ala33Gln/Thr68Gln mutants had a reduction in the binding toward C-RING1B as shown by in vitro, in silico, and in cellulo studies. This is an example of a well-folded partner of NUPR1, because its other interacting proteins are also unfolded. We hypothesize that NUPR1 plays an active role in chromatin remodeling and carcinogenesis, together with Polycomb proteins.

Variant anatomy of the biliary system as a cause of pancreatic and peri-ampullary cancers.

BACKGROUND: The cause of most pancreatic and periampullary cancers (PAC) is unknown. Recently, anatomic variations such as pancreatobiliary maljunction have been recognized as risk factors, similar to Barrett-related gastro-esophageal cancers. METHODS: Pre-operative MRI from 860 pancreatic/biliary resections, including 322 PACs, were evaluated for low-union (cystic duct joining the common hepatic duct inside of the pancreas or within 5 mm of the pancreatic border) RESULTS: Low-union, seen <10% of the population, was present in 44% of PACs (73% distal bile duct/cholangiocarcinoma, 42% pancreatic head, and 34% ampullary). It was significantly lower(11%) in conditions without connection to the ductal system (thus not exposed to the ductal/biliary tract contents), namely mucinous cystic neoplasms and intrahepatic cholangiocarcinomas(p < 0.0001). Intra-pancreatic type low-union was seen in 16% of PACs versus 2% of controls(p < 0.0001). DISCUSSION: This study establishes an association between low-union and PACs, and points to an anatomy-induced chemical/bilious carcinogenesis. This may explain why most pancreas cancers are in the head. It is possible that the same chemical milieu, caused by conditions other than low-union/insertion, may also play a role in the remaining half of PACs. This opens various treatment opportunities including milieu modifications (chemoprevention), focused screening of at-risk patients, and early detection with possible corrective actions.

Modeling post-translational modifications and cancer-associated mutations that impact the heterochromatin protein 1α-importin α heterodimers.

Heterochromatin protein 1α (HP1α) is a protein that mediates cancer-associated processes in the cell nucleus. Proteomic experiments, reported here, demonstrate that HP1α complexes with importin α (IMPα), a protein necessary for its nuclear transport. This data is congruent with Simple Linear Motif (SLiM) analyses that identify an IMPα-binding motif within the linker that joins the two globular domains of this protein. Using molecular modeling and dynamics simulations, we develop a model of the IMPα-HP1α complex and investigate the impact of phosphorylation and genomic variants on their interaction. We demonstrate that phosphorylation of the HP1α linker likely regulates its association with IMPα, which has implications for HP1α access to the nucleus, where it functions. Cancer-associated genomic variants do not abolish the interaction of HP1α but instead lead to rearrangements where the variant proteins maintain interaction with IMPα, but with less specificity. Combined, this new mechanistic insight bears biochemical, cell biological, and biomedical relevance.

P300 Acetyltransferase Mediates Stiffness-Induced Activation of Hepatic Stellate Cells Into Tumor-Promoting Myofibroblasts.

BACKGROUND & AIMS: Hepatic stellate cells (HSCs) contribute to desmoplasia and stiffness of liver metastases by differentiating into matrix-producing myofibroblasts. We investigated whether stiffness due to the presence of tumors increases activation of HSCs into myofibroblasts and their tumor-promoting effects, as well as the role of E1A binding protein p300, a histone acetyltransferase that regulates transcription, in these processes. METHODS: HSCs were isolated from liver tissues of patients, mice in which the p300 gene was flanked by 2 loxP sites (p300F/F mice), and p300+/+ mice (controls). The HSCs were placed on polyacrylamide gels with precisely defined stiffness, and their activation (differentiation into myofibroblasts) was assessed by immunofluorescence and immunoblot analyses for alpha-smooth muscle actin. In HSCs from mice, the p300 gene was disrupted by cre recombinase. In human HSCs, levels of p300 were knocked down with small hairpin RNAs or a mutant form of p300 that is not phosphorylated by AKT (p300S1834A) was overexpressed. Human HSCs were also cultured with inhibitors of p300 (C646), PI3K signaling to AKT (LY294002), or RHOA (C3 transferase) and effects on stiffness-induced activation were measured. RNA sequencing and chromatin immunoprecipitation-quantitative polymerase chain reaction were used to identify HSC genes that changed expression levels in response to stiffness. We measured effects of HSC-conditioned media on proliferation of HT29 colon cancer cells and growth of tumors following subcutaneous injection of these cells into mice. MC38 colon cancer cells were injected into portal veins of p300F/Fcre and control mice, and liver metastases were measured. p300F/Fcre and control mice were given intraperitoneal injections of CCl4 to induce liver fibrosis. Liver tissues were collected and analyzed by immunofluorescence, immunoblot, and histology. RESULTS: Substrate stiffness was sufficient to activate HSCs, leading to nuclear accumulation of p300. Disrupting p300 level or activity blocked stiffness-induced activation of HSCs. In HSCs, substrate stiffness activated AKT signaling via RHOA to induce phosphorylation of p300 at serine 1834; this caused p300 to translocate to the nucleus, where it up-regulated transcription of genes that increase activation of HSCs and metastasis, including CXCL12. MC38 cells, injected into portal veins, formed fewer metastases in livers of p300F/Fcre mice than control mice. Expression of p300 was increased in livers of mice following injection of CCl4; HSC activation and collagen deposition were reduced in livers of p300F/Fcre mice compared with control mice. CONCLUSIONS: In studies of mice, we found liver stiffness to activate HSC differentiation into myofibroblasts, which required nuclear accumulation of p300. p300 increases HSC expression of genes that promote metastasis.

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.

A Novel Kleefstra Syndrome-associated Variant That Affects the Conserved TPLX Motif within the Ankyrin Repeat of EHMT1 Leads to Abnormal Protein Folding.

Kleefstra syndrome (KS) (Mendelian Inheritance in Man (MIM) no. 610253), also known as 9q34 deletion syndrome, is an autosomal dominant disorder caused by haploinsufficiency of euchromatic histone methyltransferase-1 (EHMT1). The clinical phenotype of KS includes moderate to severe intellectual disability with absent speech, hypotonia, brachycephaly, congenital heart defects, and dysmorphic facial features with hypertelorism, synophrys, macroglossia, protruding tongue, and prognathism. Only a few cases of de novo missense mutations in EHMT1 giving rise to KS have been described. However, some EHMT1 variants have been described in individuals presenting with autism spectrum disorder or mild intellectual disability, suggesting that the phenotypic spectrum resulting from EHMT1 alterations may be quite broad. In this report, we describe two unrelated patients with complex medical histories consistent with KS in whom next generation sequencing identified the same novel c.2426C>T (p.P809L) missense variant in EHMT1 To examine the functional significance of this novel variant, we performed molecular dynamics simulations of the wild type and p.P809L variant, which predicted that the latter would have a propensity to misfold, leading to abnormal histone mark binding. Recombinant EHMT1 p.P809L was also studied using far UV circular dichroism spectroscopy and intrinsic protein fluorescence. These functional studies confirmed the model-based hypotheses and provided evidence for protein misfolding and aberrant target recognition as the underlying pathogenic mechanism for this novel KS-associated variant. This is the first report to suggest that missense variants in EHMT1 that lead to protein misfolding and disrupted histone mark binding can lead to KS.

The Role of the Histone Methyltransferase Enhancer of Zeste Homolog 2 (EZH2) in the Pathobiological Mechanisms Underlying Inflammatory Bowel Disease (IBD).

Regulatory T (Treg) cells expressing the transcription factor FOXP3 play a pivotal role in maintaining immunologic self-tolerance. We and others have shown previously that EZH2 is recruited to the FOXP3 promoter and its targets in Treg cells. To further address the role for EZH2 in Treg cellular function, we have now generated mice that lack EZH2 specifically in Treg cells (EZH2Δ/ΔFOXP3+). We find that EZH2 deficiency in FOXP3+ T cells results in lethal multiorgan autoimmunity. We further demonstrate that EZH2Δ/ΔFOXP3+ T cells lack a regulatory phenotype in vitro and secrete proinflammatory cytokines. Of special interest, EZH2Δ/ΔFOXP3+ mice develop spontaneous inflammatory bowel disease. Guided by these results, we assessed the FOXP3 and EZH2 gene networks by RNA sequencing in isolated intestinal CD4+ T cells from patients with Crohn's disease. Gene network analysis demonstrates that these CD4+ T cells display a Th1/Th17-like phenotype with an enrichment of gene targets shared by FOXP3 and EZH2. Combined, these results suggest that the inflammatory milieu found in Crohn's disease could lead to or result from deregulation of FOXP3/EZH2-enforced T cell gene networks contributing to the underlying intestinal inflammation.