Deleterious, protein-altering variants in the transcriptional coregulator ZMYM3 in 27 individuals with a neurodevelopmental delay phenotype.

Neurodevelopmental disorders (NDDs) result from highly penetrant variation in hundreds of different genes, some of which have not yet been identified. Using the MatchMaker Exchange, we assembled a cohort of 27 individuals with rare, protein-altering variation in the transcriptional coregulator ZMYM3, located on the X chromosome. Most (n = 24) individuals were males, 17 of which have a maternally inherited variant; six individuals (4 male, 2 female) harbor de novo variants. Overlapping features included developmental delay, intellectual disability, behavioral abnormalities, and a specific facial gestalt in a subset of males. Variants in almost all individuals (n = 26) are missense, including six that recurrently affect two residues. Four unrelated probands were identified with inherited variation affecting Arg441, a site at which variation has been previously seen in NDD-affected siblings, and two individuals have de novo variation resulting in p.Arg1294Cys (c.3880C>T). All variants affect evolutionarily conserved sites, and most are predicted to damage protein structure or function. ZMYM3 is relatively intolerant to variation in the general population, is widely expressed across human tissues, and encodes a component of the KDM1A-RCOR1 chromatin-modifying complex. ChIP-seq experiments on one variant, p.Arg1274Trp, indicate dramatically reduced genomic occupancy, supporting a hypomorphic effect. While we are unable to perform statistical evaluations to definitively support a causative role for variation in ZMYM3, the totality of the evidence, including 27 affected individuals, recurrent variation at two codons, overlapping phenotypic features, protein-modeling data, evolutionary constraint, and experimentally confirmed functional effects strongly support ZMYM3 as an NDD-associated gene.

Characterization of human transcription factor function and patterns of gene regulation in HepG2 cells.

Transcription factors (TFs) are trans-acting proteins that bind cis-regulatory elements (CREs) in DNA to control gene expression. Here, we analyzed the genomic localization profiles of 529 sequence-specific TFs and 151 cofactors and chromatin regulators in the human cancer cell line HepG2, for a total of 680 broadly termed DNA-associated proteins (DAPs). We used this deep collection to model each TF's impact on gene expression, and identified a cohort of 26 candidate transcriptional repressors. We examine high occupancy target (HOT) sites in the context of three-dimensional genome organization and show biased motif placement in distal-promoter connections involving HOT sites. We also found a substantial number of closed chromatin regions with multiple DAPs bound, and explored their properties, finding that a MAFF/MAFK TF pair correlates with transcriptional repression. Altogether, these analyses provide novel insights into the regulatory logic of the human cell line HepG2 genome and show the usefulness of large genomic analyses for elucidation of individual TF functions.

Clinical and cellular phenotypes resulting from a founder mutation in IL10RB.

Inborn errors of immunity are a group of rare genetically determined diseases that impair immune system development or function. Many of these diseases include immune dysregulation, autoimmunity or autoinflammation as prominent clinical features. In some children diagnosed with very early onset inflammatory bowel disease (VEOIBD), monogenic inborn errors of immune dysregulation underlie disease. We report a case of VEOIBD caused by a novel homozygous loss of function mutation in IL10RB. We use CyTOF with a broad panel of antibodies to interrogate the immunophenotype of this patient and detect reduced frequencies of CD4 and CD8 T cells with additional defects in some populations of T helper cells, innate-like T cells and memory B cells. Finally, we identify the patient's mutation as a founder allele in an isolated indigenous population and estimate the age of this variant by studying the shared ancestral haplotype.

Molecular and clinical characterization of a founder mutation causing G6PC3 deficiency.

Background: G6PC3 deficiency is a rare genetic disorder that causes syndromic congenital neutropenia. It is driven by the intracellular accumulation of a metabolite named 1,5-anhydroglucitol-6-phosphate (1,5-AG6P) that inhibits glycolysis. Patients display heterogeneous extra-hematological manifestations, contributing to delayed diagnosis. Objective: The G6PC3 c.210delC variant has been identified in patients of Mexican origin. We set out to study the origin and functional consequence of this mutation. Furthermore, we sought to characterize the clinical phenotypes caused by it. Methods: Using whole-genome sequencing data, we conducted haplotype analysis to estimate the age of this allele and traced its ancestral origin. We examined how this mutation affected G6PC3 protein expression and performed extracellular flux assays on patient-derived cells to characterize how this mutation impacts glycolysis. Finally, we compared the clinical presentations of patients with the c.210delC mutation relative to other G6PC3 deficient patients published to date. Results: Based on the length of haplotypes shared amongst ten carriers of the G6PC3 c.210delC mutation, we estimated that this variant originated in a common ancestor of indigenous American origin. The mutation causes a frameshift that introduces a premature stop codon, leading to a complete loss of G6PC3 protein expression. When treated with 1,5-anhydroglucitol (1,5-AG), the precursor to 1,5-AG6P, patient-derived cells exhibited markedly reduced engagement of glycolysis. Clinically, c.210delC carriers display all the clinical features of syndromic severe congenital neutropenia type 4 observed in prior reports of G6PC3 deficiency. Conclusion: The G6PC3 c.210delC is a loss-of-function mutation that arose from a founder effect in the indigenous Mexican population. These findings may facilitate the diagnosis of additional patients in this geographical area. Moreover, the in vitro 1,5-AG-dependent functional assay used in our study could be employed to assess the pathogenicity of additional G6PC3 variants.

Polygenic risk score for ulcerative colitis predicts immune checkpoint inhibitor-mediated colitis.

Immune checkpoint inhibitor-mediated colitis (IMC) is a common adverse event of treatment with immune checkpoint inhibitors (ICI). We hypothesize that genetic susceptibility to Crohn's disease (CD) and ulcerative colitis (UC) predisposes to IMC. In this study, we first develop a polygenic risk scores for CD (PRSCD) and UC (PRSUC) in cancer-free individuals and then test these PRSs on IMC in a cohort of 1316 patients with ICI-treated non-small cell lung cancer and perform a replication in 873 ICI-treated pan-cancer patients. In a meta-analysis, the PRSUC predicts all-grade IMC (ORmeta=1.35 per standard deviation [SD], 95% CI = 1.12-1.64, P = 2×10-03) and severe IMC (ORmeta=1.49 per SD, 95% CI = 1.18-1.88, P = 9×10-04). PRSCD is not associated with IMC. Furthermore, PRSUC predicts severe IMC among patients treated with combination ICIs (ORmeta=2.20 per SD, 95% CI = 1.07-4.53, P = 0.03). Overall, PRSUC can identify patients receiving ICI at risk of developing IMC and may be useful to monitor patients and improve patient outcomes.

Minimum entropy framework identifies a novel class of genomic functional elements and reveals regulatory mechanisms at human disease loci.

We introduce CoRE-BED, a framework trained using 19 epigenomic features in 33 major cell and tissue types to predict cell-type-specific regulatory function. CoRE-BED identifies nine functional classes de-novo, capturing both known and new regulatory categories. Notably, we describe a previously undercharacterized class that we term Development Associated Elements (DAEs), which are highly enriched in cell types with elevated regenerative potential and distinguished by the dual presence of either H3K4me2 and H3K9ac (an epigenetic signature associated with kinetochore assembly) or H3K79me3 and H4K20me1 (a signature associated with transcriptional pause release). Unlike bivalent promoters, which represent a transitory state between active and silenced promoters, DAEs transition directly to or from a non-functional state during stem cell differentiation and are proximal to highly expressed genes. CoRE-BED's interpretability facilitates causal inference and functional prioritization. Across 70 complex traits, distal insulators account for the largest mean proportion of SNP heritability (~49%) captured by the GWAS. Collectively, our results demonstrate the value of exploring non-conventional ways of regulatory classification that enrich for trait heritability, to complement existing approaches for cis-regulatory prediction.

Polygenic risk score for ulcerative colitis predicts immune checkpoint inhibitor-mediated colitis.

Immune checkpoint inhibitors (ICIs) are a remarkable advancement in cancer therapeutics; however, a substantial proportion of patients develop severe immune-related adverse events (irAEs). Understanding and predicting irAEs is a key to advancing precision immuno-oncology. Immune checkpoint inhibitor-mediated colitis (IMC) is a significant complication from ICI and can have life-threatening consequences. Based on clinical presentation, IMC mimics inflammatory bowel disease, however the link is poorly understood. We hypothesized that genetic susceptibility to Crohn's disease (CD) and ulcerative colitis (UC) may predispose to IMC. We developed and validated polygenic risk scores for CD (PRSCD) and UC (PRSUC) in cancer-free individuals and assessed the role of each of these PRSs on IMC in a cohort of 1,316 patients with non-small cell lung cancer who received ICIs. Prevalence of all-grade IMC in our cohort was 4% (55 cases), and for severe IMC, 2.5% (32 cases). The PRSUC predicted the development of all-grade IMC (HR=1.34 per standard deviation [SD], 95% CI=1.02-1.76, P=0.04) and severe IMC (HR=1.62 per SD, 95% CI=1.12-2.35, P=0.01). PRSCD was not associated with IMC or severe IMC. The association between PRSUC and IMC (all-grade and severe) was consistent in an independent pan-cancer cohort of patients treated with ICIs. Furthermore, PRSUC predicted severe IMC among patients treated with combination ICIs (OR = 2.20 per SD, 95% CI = 1.07-4.53, P=0.03). This is the first study to demonstrate the potential clinical utility of a PRS for ulcerative colitis in identifying patients receiving ICI at high risk of developing IMC, where risk reduction and close monitoring strategies could help improve overall patient outcomes.

Epitope Tagging ChIP-Seq of DNA Binding Proteins Using CETCh-Seq.

Chromatin immunoprecipitation followed by next-generation DNA sequencing (ChIP-seq) has been used to identify transcription factor (TF) binding proteins throughout the genome. Unfortunately, this approach traditionally requires commercially available, ChIP-seq grade antibodies that frequently fail to generate acceptable datasets. To obtain data for the many TFs for which there is no appropriate antibody, we recently developed a new method for performing ChIP-seq by epitope tagging endogenous TFs using CRISPR/Cas9 genome editing technology (CETCh-seq). Here, we describe our general protocol of CETCh-seq for both adherent and nonadherent cell lines using a commercially available FLAG antibody.