Glioma intratumoral heterogeneity enables adaptation to challenging microenvironments and contributes to therapeutic resistance. We integrated 914 single-cell DNA methylomes, 55,284 single-cell transcriptomes, and bulk multi-omic profiles across 11 adult IDH-mutant or IDH-wild-type gliomas to delineate sources of intratumoral heterogeneity. We show that local DNA methylation disorder associates with cell-to-cell DNA methylation differences, is elevated in more aggressive tumors, links with transcriptional disruption, and is altered during environmental stress response. Glioma cells under in vitro hypoxic and irradiation stress increased local DNA methylation disorder and shifted cell states. We identified a positive association between genetic and epigenetic instability that was supported in bulk longitudinally collected DNA methylation data. Increased DNA methylation disorder associated with accelerated disease progression, and recurrently selected DNA methylation changes were enriched for environmental stress response pathways. Our work identifies an epigenetically facilitated adaptive stress response process and highlights the importance of epigenetic heterogeneity in shaping therapeutic outcomes.

Lamin B1 is a component of the nuclear lamina and plays a critical role in maintaining nuclear architecture, regulating gene expression and modulating chromatin positioning. We have previously shown that LMNB1 gene duplications cause autosomal dominant leukodystrophy (ADLD), a fatal adult onset demyelinating disease. The mechanisms by which increased LMNB1 levels cause ADLD are unclear. To address this, we used a transgenic mouse model where Lamin B1 overexpression is targeted to oligodendrocytes. These mice showed severe vacuolar degeneration of the spinal cord white matter together with marked astrogliosis, microglial infiltration, and secondary axonal damage. Oligodendrocytes in the transgenic mice revealed alterations in histone modifications favoring a transcriptionally repressed state. Chromatin changes were accompanied by reduced expression of genes involved in lipid synthesis pathways, many of which are known to play important roles in myelin regulation and are preferentially expressed in oligodendrocytes. Decreased lipogenic gene expression resulted in a significant reduction in multiple classes of lipids involved in myelin formation. Many of these gene expression changes and lipid alterations were observed even before the onset of the phenotype, suggesting a causal role. Our findings establish, for the first time, a link between LMNB1 and lipid synthesis in oligodendrocytes, and provide a mechanistic framework to explain the age dependence and white matter involvement of the disease phenotype. These results have implications for disease pathogenesis and may also shed light on the regulation of lipid synthesis pathways in myelin maintenance and turnover. Key words: chromatin; demyelination; gene expression; inflammation; Lamin B1; lipid Significance StatementAutosomal dominant leukodystrophy (ADLD) is fatal neurological disorder caused by increased levels of the nuclear protein, Lamin B1. The disease is characterized by an age-dependent loss of myelin, the fatty sheath that covers nerve fibers. We have studied a mouse model where Lamin B1 level are increased in oligodendrocytes, the cell type that produces myelin in the CNS. We demonstrate that destruction of myelin in the spinal cord is responsible for the degenerative phenotype in our mouse model. We show that this degeneration is mediated by reduced expression of lipid synthesis genes and the subsequent reduction in myelin enriched lipids. These findings provide a mechanistic framework to explain the age dependence and tissue specificity of the ADLD disease phenotype.

BACKGROUND & AIMS Subgroups of colorectal carcinomas (CRCs) characterized by DNA methylation anomalies are termed CpG island methylator phenotype (CIMP)1, CIMP2, or CIMP-negative. The pathogenesis of CIMP1 colorectal carcinomas, and their effects on patients’ prognoses and responses to treatment, differ from those of other CRCs. We sought to identify genetic somatic alterations associated with CIMP1 CRCs. METHODS We examined genomic DNA samples from 100 primary CRCs, 10 adenomas, and adjacent normal-appearing mucosae from patients undergoing surgery or colonoscopy at 3 tertiary medical centers. We performed exome sequencing of 16 colorectal tumors and their adjacent normal tissues. Extensive comparison with known somatic alterations in CRCs allowed segregation of CIMP1-exclusive alterations. The prevalence of mutations in selected genes was determined from an independent cohort. RESULTS We found that genes that regulate chromatin were mutated in CIMP1 CRCs; the highest rates of mutation were observed in CHD7 and CHD8, which encode members of the chromodomain helicase/adenosine triphosphate—dependent chromatin remodeling family. Somaticmutations in these 2 genes were detected in 5 of 9 CIMP1 CRCs. A prevalence screen showed that nonsilencing mutations in CHD7 and CHD8 occurred significantly more frequently in CIMP1 tumors (18 of 42 [43%]) than in CIMP2 (3 of 34 [9%]; P < .01) or CIMP-negative tumors (2 of 34 [6%]; P < .001). CIMP1 markers had increased binding by CHD7, compared with all genes. Genes altered in patients with CHARGE syndrome (congenital malformations involving the central nervous system, eye, ear, nose, and mediastinal organs) who had CHD7 mutations were also altered in CRCs with mutations in CHD7. CONCLUSIONS Aberrations in chromatin remodeling could contribute to the development of CIMP1 CRCs. A better understanding of the biological determinants of CRCs can be achieved when these tumors are categorized according to their epigenetic status.

Immunotherapies revolutionized cancer treatment by harnessing the immune system to target cancer cells. However, most patients are resistant to immunotherapies and the mechanisms underlying this resistant is still poorly understood. Here, we report that overexpression of BMP7, a member of the TGFB superfamily, represents a mechanism for resistance to anti-PD1 therapy in preclinical models and in patients with disease progression while on immunotherapies. BMP7 secreted by tumor cells acts on macrophages and CD4+ T cells in the tumor microenvironment, inhibiting MAPK14 expression and impairing pro-inflammatory responses. Knockdown of BMP7 or its neutralization via follistatin in combination with anti-PD1 re-sensitizes resistant tumors to immunotherapies. Thus, we identify the BMP7 signaling pathway as a potential immunotherapeutic target in cancer.

In myotonic dystrophy type 1 (DM1), somatic mosaicism of the (CTG)n repeat expansion is age-dependent, tissue-specific, and expansion-biased. These features contribute toward variation in disease severity and confound genotype to phenotype analyses. To investigate how the (CTG)n repeat expansion changes over time, we collected three longitudinal blood DNA samples separated by 8 to 15 years and used small pool and single molecule PCR in 43 DM1 patients. We used the lower boundary of the allele length distribution as the best estimate for the inherited progenitor allele length (ePAL), which is itself the best predictor of disease severity. Although in most patients the lower boundary of the allele length distribution was conserved over time, in many this estimate also increased with age, suggesting samples for research studies and clinical trials should be obtained as early as possible. As expected, the modal allele length increased over time, driven primarily by ePAL, age-at-sampling and the time interval. As expected, small expansions &lt; 100 repeats, did not expand as rapidly as larger alleles. However, the rate of expansion of very large alleles was not obviously proportionally higher. This may, at least in part, be a result of the allele length-dependent increase in large contractions that we also observed. We also determined that individual-specific variation in the increase of modal allele length over time not accounted for by ePAL, age-at-sampling and time was inversely associated with individual-specific variation in age-at-onset not accounted for by ePAL, further highlighting somatic expansion as a therapeutic target in DM1.