Histone methylation antagonism drives tumor immune evasion in squamous cell carcinomas.

How cancer-associated chromatin abnormalities shape tumor-immune interaction remains incompletely understood. Recent studies have linked DNA hypomethylation and de-repression of retrotransposons to anti-tumor immunity through the induction of interferon response. Here, we report that inactivation of the histone H3K36 methyltransferase NSD1, which is frequently found in squamous cell carcinomas (SCCs) and induces DNA hypomethylation, unexpectedly results in diminished tumor immune infiltration. In syngeneic and genetically engineered mouse models of head and neck SCCs, NSD1-deficient tumors exhibit immune exclusion and reduced interferon response despite high retrotransposon expression. Mechanistically, NSD1 loss results in silencing of innate immunity genes, including the type III interferon receptor IFNLR1, through depletion of H3K36 di-methylation (H3K36me2) and gain of H3K27 tri-methylation (H3K27me3). Inhibition of EZH2 restores immune infiltration and impairs the growth of Nsd1-mutant tumors. Thus, our work uncovers a druggable chromatin cross talk that regulates the viral mimicry response and enables immune evasion of DNA hypomethylated tumors.

H3K36 dimethylation shapes the epigenetic interaction landscape by directing repressive chromatin modifications in embryonic stem cells.

Epigenetic modifications on the chromatin do not occur in isolation. Chromatin-associated proteins and their modification products form a highly interconnected network, and disturbing one component may rearrange the entire system. We see this increasingly clearly in epigenetically dysregulated cancers. It is important to understand the rules governing epigenetic interactions. Here, we use the mouse embryonic stem cell (mESC) model to describe in detail the relationships within the H3K27-H3K36-DNA methylation subnetwork. In particular, we focus on the major epigenetic reorganization caused by deletion of the histone 3 lysine 36 methyltransferase NSD1, which in mESCs deposits nearly all of the intergenic H3K36me2. Although disturbing the H3K27 and DNA methylation (DNAme) components also affects this network to a certain extent, the removal of H3K36me2 has the most drastic effect on the epigenetic landscape, resulting in full intergenic spread of H3K27me3 and a substantial decrease in DNAme. By profiling DNMT3A and CHH methylation (mCHH), we show that H3K36me2 loss upon Nsd1-KO leads to a massive redistribution of DNMT3A and mCHH away from intergenic regions and toward active gene bodies, suggesting that DNAme reduction is at least in part caused by redistribution of de novo methylation. Additionally, we show that pervasive acetylation of H3K27 is regulated by the interplay of H3K36 and H3K27 methylation. Our analysis highlights the importance of H3K36me2 as a major determinant of the developmental epigenome and provides a framework for further consolidating our knowledge of epigenetic networks.

The histone mark H3K36me2 recruits DNMT3A and shapes the intergenic DNA methylation landscape.

Enzymes that catalyse CpG methylation in DNA, including the DNA methyltransferases 1 (DNMT1), 3A (DNMT3A) and 3B (DNMT3B), are indispensable for mammalian tissue development and homeostasis1-4. They are also implicated in human developmental disorders and cancers5-8, supporting the critical role of DNA methylation in the specification and maintenance of cell fate. Previous studies have suggested that post-translational modifications of histones are involved in specifying patterns of DNA methyltransferase localization and DNA methylation at promoters and actively transcribed gene bodies9-11. However, the mechanisms that control the establishment and maintenance of intergenic DNA methylation remain poorly understood. Tatton-Brown-Rahman syndrome (TBRS) is a childhood overgrowth disorder that is defined by germline mutations in DNMT3A. TBRS shares clinical features with Sotos syndrome (which is caused by haploinsufficiency of NSD1, a histone methyltransferase that catalyses the dimethylation of histone H3 at K36 (H3K36me2)8,12,13), which suggests that there is a mechanistic link between these two diseases. Here we report that NSD1-mediated H3K36me2 is required for the recruitment of DNMT3A and maintenance of DNA methylation at intergenic regions. Genome-wide analysis shows that the binding and activity of DNMT3A colocalize with H3K36me2 at non-coding regions of euchromatin. Genetic ablation of Nsd1 and its paralogue Nsd2 in mouse cells results in a redistribution of DNMT3A to H3K36me3-modified gene bodies and a reduction in the methylation of intergenic DNA. Blood samples from patients with Sotos syndrome and NSD1-mutant tumours also exhibit hypomethylation of intergenic DNA. The PWWP domain of DNMT3A shows dual recognition of H3K36me2 and H3K36me3 in vitro, with a higher binding affinity towards H3K36me2 that is abrogated by TBRS-derived missense mutations. Together, our study reveals a trans-chromatin regulatory pathway that connects aberrant intergenic CpG methylation to human neoplastic and developmental overgrowth.

Independent mutations at Arg181 and Arg274 of Vangl proteins that are associated with neural tube defects in humans decrease protein stability and impair membrane targeting.

In vertebrates, Vangl proteins play important roles during embryogenesis, including establishing planar polarity and coordinating convergent extension movements. In mice, homozygosity for mutations in the Vangl1 and Vangl2 genes or combined heterozygosity for Vangl1/Vangl2 mutations causes the very severe neural tube defect (NTD) craniorachischisis. Recently, a number of patient-specific VANGL1 and VANGL2 protein mutations have been identified in familial and sporadic cases of mild and severe forms of NTDs. The biochemical nature of pathological effects in these mutations remains unknown. Of interest are two arginine residues, R181 and R274, that are highly conserved in Vangl protein homologues and found to be independently mutated in VANGL1 (R181Q and R274Q) and VANGL2 (R177H and R270H) in human cases of NTDs. The cellular and biochemical properties of R181Q and R274Q were established in transfected MDCK kidney epithelial cells and compared to those of wild-type (WT) Vangl1. Compared to that of WT, these mutations displayed impaired targeting to the plasma membrane and were instead detected in an intracellular endomembrane compartment that was positive for the endoplasmic reticulum. R181Q and R274Q showed impaired stability with significant reductions in measured half-lives from >20 h for WT protein to 9 and 5 h, respectively. These mutations have a cellular and biochemical phenotype that is indistinguishable from that of Vangl mutations known to cause craniorachichisis in mice (Lp). These results strongly suggest that R181 and R274 play critical roles in Vangl protein function and that their mutations cause neural tube defects in humans.

H2A.Z histone variants facilitate HDACi-dependent removal of H3.3K27M mutant protein in pediatric high-grade glioma cells.

Diffuse intrinsic pontine gliomas (DIPGs) are deadly pediatric brain tumors, non-resectable due to brainstem localization and diffusive growth. Over 80% of DIPGs harbor a mutation in histone 3 (H3.3 or H3.1) resulting in a lysine-to-methionine substitution (H3K27M). Patients with DIPG have a dismal prognosis with no effective therapy. We show that histone deacetylase (HDAC) inhibitors lead to a significant reduction in the H3.3K27M protein (up to 80%) in multiple glioma cell lines. We discover that the SB939-mediated H3.3K27M loss is partially blocked by a lysosomal inhibitor, chloroquine. The H3.3K27M loss is facilitated by co-occurrence of H2A.Z, as evidenced by the knockdown of H2A.Z isoforms. Chromatin immunoprecipitation sequencing (ChIP-seq) analysis confirms the occupancy of H3.3K27M and H2A.Z at the same SB939-inducible genes. We discover a mechanism showing that HDAC inhibition in DIPG leads to pharmacological modulation of the oncogenic H3.3K27M protein levels. These findings show the possibility of directly targeting the H3.3K27M oncohistone.

Systematic perturbations of SETD2, NSD1, NSD2, NSD3 and ASH1L reveals their distinct contributions to H3K36 methylation.

Background: Methylation of histone 3 lysine 36 (H3K36me) has emerged as an essential epigenetic component for the faithful regulation of gene expression. Despite its importance in development, disease, and cancer, how the molecular agents collectively shape the H3K36me landscape is unclear. Results: We use a mouse mesenchymal stem cell model to perturb the H3K36me deposition machinery and infer the activities of the five most prominent players: SETD2, NSD1, NSD2, NSD3, and ASH1L. We find that H3K36me2 is the most abundant of the three methylation states and is predominantly deposited at intergenic regions by NSD1, and partly by NSD2. In contrast, H3K36me1/3 are most abundant within exons and are positively correlated with gene expression. We demonstrate that while SETD2 deposits most H3K36me3, it also deposits H3K36me2 within transcribed genes. Additionally, loss of SETD2 results in an increase of exonic H3K36me1, suggesting other H3K36 methyltransferases (K36MTs) prime gene bodies with lower methylation states ahead of transcription. Through a reductive approach, we uncover the distribution patterns of NSD3- and ASH1L-catalyzed H3K36me2. While NSD1/2 establish broad intergenic H3K36me2 domains, NSD3 deposits H3K36me2 peaks on active promoters and enhancers. Meanwhile, the activity of ASH1L is restricted to the regulatory elements of developmentally relevant genes, and our analyses implicate PBX2 as a potential recruitment factor. Conclusions: Within genes, SETD2 deposits both H3K36me2/3, while the other K36MTs are capable of depositing H3K36me1/2 independently of SETD2 activity. For the deposition of H3K36me1/2, we find a hierarchy of K36MT activities where NSD1>NSD2>NSD3>ASH1L. While NSD1 and NSD2 are responsible for most genome-wide propagation of H3K36me2, the activities of NSD3 and ASH1L are confined to active regulatory elements.

Transmembrane topology of mammalian planar cell polarity protein Vangl1.

Vangl1 and Vangl2 are membrane proteins that play an important role in neurogenesis, and Vangl1/Vangl2 mutations cause neural tube defects in mice and humans. At the cellular level, Vangl proteins regulate the establishment of planar cell polarity (PCP), a process requiring membrane assembly of asymmetrically distributed multiprotein complexes that transmit polarity information to neighboring cells. The membrane topology of Vangl proteins and the protein segments required for structural and functional aspects of multiprotein membrane PCP complexes is unknown. We have used epitope tagging and immunofluorescence to establish the secondary structure of Vangl proteins, including the number, position, and polarity of transmembrane domains. Antigenic hemagglutinin A (HA) peptides (YYDVPDYS) were inserted in predicted intra- or extracellular loops of Vangl1 at positions 18, 64, 139, 178, 213, and 314, and individual mutant variants were stably expressed at the membrane of MDCK polarized cells. The membrane topology of the exofacial HA tag was determined by a combination of immunofluorescence in intact (extracellular epitopes) and permeabilized (intracellular epitopes) cells and by surface labeling. Results indicate that Vangl proteins have a four-transmembrane domain structure with the N-terminal portion (HA18 and HA64) and the large C-terminal portion (HA314) of the protein being intracellular. Topology mapping and hydropathy profiling show that the loop delineated by TMD1-2 (HA139) and TMD3-4 (HA213) is extracellular while the segment separating predicted TMD2-3 (HA178) is intracellular. This secondary structure reveals a compact membrane-associated portion with very short predicted extra- and intracellular loops, while the protein harbors a large intracellular domain.

Loss of membrane targeting of Vangl proteins causes neural tube defects.

In the mouse, the loop-tail mutation (Lp) causes a very severe neural tube defect, which is caused by mutations in the Vangl2 gene. In mammals, Vangl1 and Vangl2 code for integral membrane proteins that assemble into asymmetrically distributed membrane complexes that establish planar cell polarity in epithelial cells and that regulate convergent extension movements during embryogenesis. To date, VANGL are the only genes in which mutations cause neural tube defects in humans. Three independently arising Lp alleles have been described for Vangl2: D255E, S464N, and R259L. Here we report a common mechanism for both the naturally occurring Lp (S464N) and a novel ENU-induced mutation Lp(m2Jus)(R259L). We show that the S464N and R259L variants stably expressed in polarized MDCK kidney cells fail to reach the plasma membrane, their site for biological function. The mutant variants are retained intracellularly in the endoplasmic reticulum, colocalizing with ER chaperone calreticulin. Furthermore, the mutants also show a dramatically reduced half-life of ∼3 h, compared to ∼22 h for the wild-type protein, and are rapidly degraded in a proteasome-dependent and MG132-sensitive fashion. Coexpressing individually the three known allelic Lp variants with the wild-type protein does not influence the localization of the WT at the plasma membrane, suggesting that the codominant nature of the Lp trait in vivo is due to haploid insufficiency caused by a partial loss of function in a gene dosage-dependent pathway, as opposed to a dominant negative phenotype. Our study provides a biochemical framework for the study of recently identified mutations in hVANGL1 and hVANGL2 in sporadic or familial cases of neural tube defects.

Chromatin dysregulation associated with NSD1 mutation in head and neck squamous cell carcinoma.

Chromatin dysregulation has emerged as an important mechanism of oncogenesis. To develop targeted treatments, it is important to understand the transcriptomic consequences of mutations in chromatin modifier genes. Recently, mutations in the histone methyltransferase gene nuclear receptor binding SET domain protein 1 (NSD1) have been identified in a subset of common and deadly head and neck squamous cell carcinomas (HNSCCs). Here, we use genome-wide approaches and genome editing to dissect the downstream effects of loss of NSD1 in HNSCC. We demonstrate that NSD1 mutations are responsible for loss of intergenic H3K36me2 domains, followed by loss of DNA methylation and gain of H3K27me3 in the affected genomic regions. In addition, those regions are enriched in cis-regulatory elements, and subsequent loss of H3K27ac correlates with reduced expression of their target genes. Our analysis identifies genes and pathways affected by the loss of NSD1 and paves the way to further understanding the interplay among chromatin modifications in cancer.

Two competing mechanisms of DNMT3A recruitment regulate the dynamics of de novo DNA methylation at PRC1-targeted CpG islands.

Precise deposition of CpG methylation is critical for mammalian development and tissue homeostasis and is often dysregulated in human diseases. The localization of de novo DNA methyltransferase DNMT3A is facilitated by its PWWP domain recognizing histone H3 lysine 36 (H3K36) methylation1,2 and is normally depleted at CpG islands (CGIs)3. However, methylation of CGIs regulated by Polycomb repressive complexes (PRCs) has also been observed4-8. Here, we report that DNMT3A PWWP domain mutations identified in paragangliomas9 and microcephalic dwarfism10 promote aberrant localization of DNMT3A to CGIs in a PRC1-dependent manner. DNMT3A PWWP mutants accumulate at regions containing PRC1-mediated formation of monoubiquitylated histone H2A lysine 119 (H2AK119ub), irrespective of the amounts of PRC2-catalyzed formation of trimethylated histone H3 lysine 27 (H3K27me3). DNMT3A interacts with H2AK119ub-modified nucleosomes through a putative amino-terminal ubiquitin-dependent recruitment region, providing an alternative form of DNMT3A genomic targeting that is augmented by the loss of PWWP reader function. Ablation of PRC1 abrogates localization of DNMT3A PWWP mutants to CGIs and prevents aberrant DNA hypermethylation. Our study implies that a balance between DNMT3A recruitment by distinct reader domains guides de novo CpG methylation and may underlie the abnormal DNA methylation landscapes observed in select human cancer subtypes and developmental disorders.

Molecular and cellular mechanisms underlying neural tube defects in the loop-tail mutant mouse.

Loop-tail (Lp) mice show a very severe neural tube defect (craniorachischisis) caused by mutations in the Vangl2 gene (D255E, S464N). Mammalian Vangl1 and Vangl2 are membrane proteins that play critical roles in development such as establishment of planar cell polarity (PCP) in epithelial layers and convergent extension movements during neurogenesis and cardiogenesis. Vangl proteins are thought to assemble with other PCP proteins (Dvl, Pk) to form membrane-bound PCP signaling complexes that provide polarity information to the cell. In the present study, we show that Vangl1 is expressed exclusively at the plasma membrane of transfected MDCK cells, where it is targeted to the basolateral membrane. Experiments with an inserted exofacial HA epitope indicate that the segment delimited by the predicted transmembrane domains 1 and 2 is exposed to the extracellular milieu. Comparative studies of the Lp-associated pathogenic mutation D255E indicate that the targeting of the mutant variant at the plasma membrane is greatly reduced; the mutant variant is predominantly retained intracellularly in endoplasmic reticulum (ER) vesicles colocalizing with the ER marker calreticulin. In addition, the D255E variant shows drastically reduced stability with a half-life of approximately 2 h, compared to >9 h for its wild type counterpart and is rapidly degraded in a proteasome-dependent and MG132 sensitive pathway. These findings highlight a critical role for D255 for normal folding and processing of Vangl proteins, with highly conservative substitutions not tolerated at that site. Our study provide an experimental framework for the analysis of human VANGL mutations recently identified in familial and sporadic cases of spina bifida.

H3K27M in Gliomas Causes a One-Step Decrease in H3K27 Methylation and Reduced Spreading within the Constraints of H3K36 Methylation.

The discovery of H3K27M mutations in pediatric gliomas marked a new chapter in cancer epigenomics. Numerous studies have investigated the effect of this mutation on H3K27 trimethylation, but only recently have we started to realize its additional effects on the epigenome. Here, we use isogenic glioma H3K27M+/- cell lines to investigate H3K27 methylation and its interaction with H3K36 and H3K9 modifications. We describe a "step down" effect of H3K27M on the distribution of H3K27 methylation: me3 is reduced to me2, me2 is reduced to me1, whereas H3K36me2/3 delineates the boundaries for the spread of H3K27me marks. We also observe a replacement of H3K27me2/3 silencing by H3K9me3. Using a computational simulation, we explain our observations by reduced effectiveness of PRC2 and constraints imposed on the deposition of H3K27me by antagonistic H3K36 modifications. Our work further elucidates the effects of H3K27M in gliomas as well as the general principles of deposition in H3K27 methylation.

Transcriptome and regulatory maps of decidua-derived stromal cells inform gene discovery in preterm birth.

While a genetic component of preterm birth (PTB) has long been recognized and recently mapped by genome-wide association studies (GWASs), the molecular determinants underlying PTB remain elusive. This stems in part from an incomplete availability of functional genomic annotations in human cell types relevant to pregnancy and PTB. We generated transcriptome (RNA-seq), epigenome (ChIP-seq of H3K27ac, H3K4me1, and H3K4me3 histone modifications), open chromatin (ATAC-seq), and chromatin interaction (promoter capture Hi-C) annotations of cultured primary decidua-derived mesenchymal stromal/stem cells and in vitro differentiated decidual stromal cells and developed a computational framework to integrate these functional annotations with results from a GWAS of gestational duration in 56,384 women. Using these resources, we uncovered additional loci associated with gestational duration and target genes of associated loci. Our strategy illustrates how functional annotations in pregnancy-relevant cell types aid in the experimental follow-up of GWAS for PTB and, likely, other pregnancy-related conditions.