Since the discovery of 5-hydroxymethylcytosine (5hmC) as a prominent DNA modification found in mammalian genomes, an emergent question has been what role this mark plays in gene regulation. 5hmC is hypothesized to function as an intermediate in the demethylation of 5-methylcytosine (5mC) and in the reactivation of silenced promoters and enhancers. Further, weak positive correlations are observed between gene body 5hmC and gene expression. We previously demonstrated that ME-Class is an effective tool to understand relationships between whole-genome bisulfite sequencing data and expression. In this work, we present ME-Class2, a machine-learning based tool to perform integrative 5mCG, 5hmCG and expression analysis. Using ME-Class2 we analyze whole-genome single-base resolution 5mCG and 5hmCG datasets from 20 primary tissue and cell samples to reveal relationships between 5hmCG and expression. Our analysis indicates that conversion of 5mCG to 5hmCG within 2 kb of the transcription start site associates with distinct functions depending on the summed level of 5mCG + 5hmCG. Unchanged levels of 5mCG + 5hmCG (conversion from 5mCG to stable 5hmCG) associate with repression. Meanwhile, decreases in 5mCG + 5hmCG (5hmCG-mediated demethylation) associate with gene activation. Our results demonstrate that ME-Class2 will prove invaluable to interpret genome-wide 5mC and 5hmC datasets and guide mechanistic studies into the function of 5hmCG.

Epigenetic changes, such as aberrant DNA methylation, contribute to cancer clonal expansion and disease progression. However, identifying subpopulation-level changes in a heterogeneous sample remains challenging. Thus, we have developed a computational approach, DXM, to deconvolve the methylation profiles of major allelic subpopulations from the bisulfite sequencing data of a heterogeneous sample. DXM does not require prior knowledge of the number of subpopulations or types of cells to expect. We benchmark DXM's performance and demonstrate improvement over existing methods. We further experimentally validate DXM predicted allelic subpopulation-methylation profiles in four Diffuse Large B-Cell Lymphomas (DLBCLs). Lastly, as proof-of-concept, we apply DXM to a cohort of 31 DLBCLs and relate allelic subpopulation methylation profiles to relapse. We thus demonstrate that DXM can robustly find allelic subpopulation methylation profiles that may contribute to disease progression using bisulfite sequencing data of any heterogeneous sample.

Almost all CpG-rich promoters in the mammalian genome are bound by the multidomain FBXL10 protein (also known as KDM2B, JHDM1B, CXXC2, and NDY1). FBXL10 is expressed as two isoforms: FBXL10-1, a longer form that contains an N-terminal histone demethylase domain with C-terminal F-box, CXXC, PHD, RING, and leucine-rich repeat domains, and FBXL10-2, a shorter form that initiates at an alternative internal exon and which lacks the histone demethylase domain but retains all other annotated domains. Selective deletion of Fbxl10-1 had been reported to produce a low penetrance and variable phenotype; most of the mutant animals were essentially normal. We constructed mutant mouse strains that were either null for Fbxl10-2 but wild type for Fbxl10-1 or null for both Fbxl10-1 and Fbxl10-2.

Age is a significant risk factor for the development of cancer. However, the mechanisms that drive age-related increases in cancer remain poorly understood. To determine if senescent stromal cells influence tumorigenesis, we develop a mouse model that mimics the aged skin microenvironment. Using this model, here we find that senescent stromal cells are sufficient to drive localized increases in suppressive myeloid cells that contributed to tumour promotion. Further, we find that the stromal-derived senescence-associated secretory phenotype factor interleukin-6 orchestrates both increases in suppressive myeloid cells and their ability to inhibit anti-tumour T-cell responses. Significantly, in aged, cancer-free individuals, we find similar increases in immune cells that also localize near senescent stromal cells. This work provides evidence that the accumulation of senescent stromal cells is sufficient to establish a tumour-permissive, chronic inflammatory microenvironment that can shelter incipient tumour cells, thus allowing them to proliferate and progress unabated by the immune system.

Advances in sequencing technology allow researchers to map genome-wide changes in DNA methylation in development and disease. However, there is a lack of experimental tools to site-specifically manipulate DNA methylation to discern the functional consequences. We developed a CRISPR/Cas9 DNA methyltransferase 3A (DNMT3A) fusion to induce DNA methylation at specific loci in the genome. We induced DNA methylation at up to 50% of alleles for targeted CpG dinucleotides. DNA methylation levels peaked within 50 bp of the short guide RNA (sgRNA) binding site and between pairs of sgRNAs. We used our approach to target methylation across the entire CpG island at the CDKN2A promoter, three CpG dinucleotides at the ARF promoter, and the CpG island within the Cdkn1a promoter to decrease expression of the target gene. These tools permit mechanistic studies of DNA methylation and its role in guiding molecular processes that determine cellular fate.

Molecular analyses of Aplysia, a well-established model organism for cellular and systems neural science, have been seriously handicapped by a lack of adequate genomic information. By sequencing cDNA libraries from the central nervous system (CNS), we have identified over 175,000 expressed sequence tags (ESTs), of which 19,814 are unique neuronal gene products and represent 50%-70% of the total Aplysia neuronal transcriptome. We have characterized the transcriptome at three levels: (1) the central nervous system, (2) the elementary components of a simple behavior: the gill-withdrawal reflex-by analyzing sensory, motor, and serotonergic modulatory neurons, and (3) processes of individual neurons. In addition to increasing the amount of available gene sequences of Aplysia by two orders of magnitude, this collection represents the largest database available for any member of the Lophotrochozoa and therefore provides additional insights into evolutionary strategies used by this highly successful diversified lineage, one of the three proposed superclades of bilateral animals.

To understand the mechanisms that mediate germline genetic leukemia predisposition, we studied the inherited ribosomopathy Shwachman-Diamond syndrome (SDS), a bone marrow failure disorder with high risk of myeloid malignancies at an early age. To define the mechanistic basis of clonal hematopoiesis in SDS, we investigate somatic mutations acquired by patients with SDS followed longitudinally. Here we report that multiple independent somatic hematopoietic clones arise early in life, most commonly harboring heterozygous mutations in EIF6 or TP53. We show that germline SBDS deficiency establishes a fitness constraint that drives selection of somatic clones via two distinct mechanisms with different clinical consequences. EIF6 inactivation mediates a compensatory pathway with limited leukemic potential by ameliorating the underlying SDS ribosome defect and enhancing clone fitness. TP53 mutations define a maladaptive pathway with enhanced leukemic potential by inactivating tumor suppressor checkpoints without correcting the ribosome defect. Subsequent development of leukemia was associated with acquisition of biallelic TP53 alterations. These results mechanistically link leukemia predisposition to germline genetic constraints on cellular fitness, and provide a rational framework for clinical surveillance strategies.

Human T-cell leukemia virus type 1 (HTLV-1) is a retrovirus that causes an aggressive T-cell malignancy and a variety of inflammatory conditions. The integrated provirus includes a single binding site for the epigenomic insulator, CCCTC-binding protein (CTCF), but its function remains unclear. In the current study, a mutant virus was examined that eliminates the CTCF-binding site. The mutation did not disrupt the kinetics and levels of virus gene expression, or establishment of or reactivation from latency. However, the mutation disrupted the epigenetic barrier function, resulting in enhanced DNA CpG methylation downstream of the CTCF binding site on both strands of the integrated provirus and H3K4Me3, H3K36Me3, and H3K27Me3 chromatin modifications both up- and downstream of the site. A majority of clonal cell lines infected with wild type HTLV-1 exhibited increased plus strand gene expression with CTCF knockdown, while expression in mutant HTLV-1 clonal lines was unaffected. These findings indicate that CTCF binding regulates HTLV-1 gene expression, DNA and histone methylation in an integration site dependent fashion.

Numerous genomic studies are underway to determine which genes are abnormally regulated by DNA methylation in disease. However, we have a poor understanding of how disease-specific methylation changes affect expression. We thus developed an integrative analysis tool, Methylation-based Gene Expression Classification (ME-Class), to explain specific variation in methylation that associates with expression change. This model captures the complexity of methylation changes around a gene promoter. Using 17 whole-genome bisulfite sequencing and RNA-seq datasets from different tissues from the Roadmap Epigenomics Project, ME-Class significantly outperforms standard methods using methylation to predict differential gene expression change. To demonstrate its utility, we used ME-Class to analyze 32 datasets from different hematopoietic cell types from the Blueprint Epigenome project. Expression-associated methylation changes were predominantly found when comparing cells from distantly related lineages, implying that changes in the cell's transcriptional program precede associated methylation changes. Training ME-Class on normal-tumor pairs from The Cancer Genome Atlas indicated that cancer-specific expression-associated methylation changes differ from tissue-specific changes. We further show that ME-Class can detect functionally relevant cancer-specific, expression-associated methylation changes that are reversed upon the removal of methylation. ME-Class is thus a powerful tool to identify genes that are dysregulated by DNA methylation in disease.

Pericentromeric satellite repeats are enriched in 5-methylcytosine (5mC). Loss of 5mC at these sequences is common in cancer and is a hallmark of Immunodeficiency, Centromere and Facial abnormalities (ICF) syndrome. While the general importance of 5mC is well-established, the specific functions of 5mC at pericentromeres are less clear. To address this deficiency, we generated a viable animal model of pericentromeric hypomethylation through mutation of the ICF-gene ZBTB24. Deletion of zebrafish zbtb24 caused a progressive loss of 5mC at pericentromeres and ICF-like phenotypes. Hypomethylation of these repeats triggered derepression of pericentromeric transcripts and activation of an interferon-based innate immune response. Injection of pericentromeric RNA is sufficient to elicit this response in wild-type embryos, and mutation of the MDA5-MAVS dsRNA-sensing machinery blocks the response in mutants. These findings identify activation of the innate immune system as an early consequence of pericentromeric hypomethylation, implicating derepression of pericentromeric transcripts as a trigger of autoimmunity.

Promoter silencing by ectopic de novo methylation of tumor suppressor genes has been proposed as comparable or equivalent to inactivating mutations as a factor in carcinogenesis. However, this hypotheses had not previously been tested by high resolution, high-coverage whole-genome methylation profiling in primary carcinomas. We have determined the genomic methylation status of a series of primary mammary carcinomas and matched control tissues by examination of more than 2.7 billion CpG dinucleotides. Most of the tumors showed variable losses of DNA methylation from all sequence compartments, but increases in promoter methylation were infrequent, very small in extent, and were observed largely at CpG-poor promoters. De novo methylation at the promoters of proto-oncogenes and tumor suppressor genes occurred at approximately the same frequency. The findings indicate that tumor suppressor silencing by de novo methylation is much less common than currently believed. We put forward a hypothesis under which the demethylation commonly observed in carcinomas is a manifestation of a defensive system that kills incipient cancer cells.

Myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML) are neoplastic disorders of hematopoietic stem cells. DNA methyltransferase inhibitors, 5-azacytidine and 5-aza-2'-deoxycytidine (decitabine), benefit some MDS/AML patients. However, the role of DNA methyltransferase inhibitor-induced DNA hypomethylation in regulation of gene expression in AML is unclear.

To characterize patient-derived xenografts (PDXs) for functional studies, we made whole-genome comparisons with originating breast cancers representative of the major intrinsic subtypes. Structural and copy number aberrations were found to be retained with high fidelity. However, at the single-nucleotide level, variable numbers of PDX-specific somatic events were documented, although they were only rarely functionally significant. Variant allele frequencies were often preserved in the PDXs, demonstrating that clonal representation can be transplantable. Estrogen-receptor-positive PDXs were associated with ESR1 ligand-binding-domain mutations, gene amplification, or an ESR1/YAP1 translocation. These events produced different endocrine-therapy-response phenotypes in human, cell line, and PDX endocrine-response studies. Hence, deeply sequenced PDX models are an important resource for the search for genome-forward treatment options and capture endocrine-drug-resistance etiologies that are not observed in standard cell lines. The originating tumor genome provides a benchmark for assessing genetic drift and clonal representation after transplantation.

Despite advances in early diagnosis and treatment of cancer patients, metastasis remains the major cause of mortality. TP53 is one of the most frequently mutated genes in human cancer, and these alterations can occur during the early stages of oncogenesis or as later events as tumors progress to more aggressive forms. Previous studies have suggested that p53 plays a role in cellular pathways that govern metastasis. To investigate how p53 deficiency contributes to late-stage tumor growth and metastasis, we developed paired isogenic patient-derived xenograft (PDX) models of triple-negative breast cancer (TNBC) differing only in p53 status for longitudinal analysis.

Tumor cells disseminate early in tumor development making metastasis-prevention strategies difficult. Identifying proteins that promote the outgrowth of disseminated tumor cells may provide opportunities for novel therapeutic strategies. Despite multiple studies demonstrating that the mesenchymal-to-epithelial transition (MET) is critical for metastatic colonization, key regulators that initiate this transition remain unknown. We serially passaged lung metastases from a primary triple negative breast cancer xenograft to the mammary fat pads of recipient mice to enrich for gene expression changes that drive metastasis. An unbiased transcriptomic signature of potential metastatic drivers was generated, and a high throughput gain-of-function screen was performed in vivo to validate candidates. Carcinoembryonic antigen-related cell adhesion molecule 5 (CEACAM5) was identified as a metastatic driver. CEACAM5 overproduction enriched for an epithelial gene expression pattern and facilitated tumor outgrowth at metastatic sites. Tissues from patients with metastatic breast cancer confirmed elevated levels of CEACAM5 in lung metastases relative to breast tumors, and an inverse correlation between CEACAM5 and the mesenchymal marker vimentin was demonstrated. Thus, CEACAM5 facilitates tumor outgrowth at metastatic sites by promoting MET, warranting its investigation as a therapeutic target and biomarker of aggressiveness in breast cancer.

The genetic modules that contribute to human evolution are poorly understood. Here we investigate positive selection in the Epidermal Differentiation Complex locus for skin barrier adaptation in diverse HapMap human populations (CEU, JPT/CHB, and YRI). Using Composite of Multiple Signals and iSAFE, we identify selective sweeps for LCE1A-SMCP and involucrin (IVL) haplotypes associated with human migration out-of-Africa, reaching near fixation in European populations. CEU-IVL is associated with increased IVL expression and a known epidermis-specific enhancer. CRISPR/Cas9 deletion of the orthologous mouse enhancer in vivo reveals a functional requirement for the enhancer to regulate Ivl expression in cis. Reporter assays confirm increased regulatory and additive enhancer effects of CEU-specific polymorphisms identified at predicted IRF1 and NFIC binding sites in the IVL enhancer (rs4845327) and its promoter (rs1854779). Together, our results identify a selective sweep for a cis regulatory module for CEU-IVL, highlighting human skin barrier evolution for increased IVL expression out-of-Africa.

The androgen receptor (AR) is a ligand-inducible transcription factor that mediates androgen action in target tissues. Upon ligand binding, the AR binds to thousands of genomic loci and activates a cell-type specific gene program. Prostate cancer growth and progression depend on androgen-induced AR signaling. Treatment of advanced prostate cancer through medical or surgical castration leads to initial response and durable remission, but resistance inevitably develops. In castration-resistant prostate cancer (CRPC), AR activity remains critical for tumor growth despite androgen deprivation. Although previous studies have focused on ligand-dependent AR signaling, in this study we explore AR function under the androgen-deprived conditions characteristic of CRPC. Our data demonstrate that AR persistently occupies a distinct set of genomic loci after androgen deprivation in CRPC. These androgen-independent AR occupied regions have constitutively open chromatin structures that lack the canonical androgen response element and are independent of FoxA1, a transcription factor involved in ligand-dependent AR targeting. Many AR binding events occur at proximal promoters, which can act as enhancers to augment transcriptional activities of other promoters through DNA looping. We further show that androgen-independent AR binding directs a gene expression program in CRPC, which is necessary for the growth of CRPC after androgen withdrawal.

The mechanisms by which methylated mammalian promoters are transcriptionally silenced even in the presence of all of the factors required for their expression have long been a major unresolved issue in the field of epigenetics. Repression requires the assembly of a methylation-dependent silencing complex that contains the TRIM28 protein (also known as KAP1 and TIF1β), a scaffolding protein without intrinsic repressive or DNA-binding properties. The identity of the key effector within this complex that represses transcription is unknown. We developed a methylation-sensitized interaction screen which revealed that TRIM28 was complexed with O-linked β-N-acetylglucosamine transferase (OGT) only in cells that had normal genomic methylation patterns. OGT is the only glycosyltransferase that modifies cytoplasmic and nuclear protein by transfer of N-acetylglucosamine (O-GlcNAc) to serine and threonine hydroxyls. Whole-genome analysis showed that O-glycosylated proteins and TRIM28 were specifically bound to promoters of active retrotransposons and to imprinting control regions, the two major regulatory sequences controlled by DNA methylation. Furthermore, genome-wide loss of DNA methylation caused a loss of O-GlcNAc from multiple transcriptional repressor proteins associated with TRIM28. A newly developed Cas9-based editing method for targeted removal of O-GlcNAc was directed against retrotransposon promoters. Local chromatin de-GlcNAcylation specifically reactivated the expression of the targeted retrotransposon family without loss of DNA methylation. These data revealed that O-linked glycosylation of chromatin factors is essential for the transcriptional repression of methylated retrotransposons.

Primary immune regulatory disorders (PIRD) represent a group of disorders characterized by immune dysregulation, presenting with a wide range of clinical disease, including autoimmunity, autoinflammation, or lymphoproliferation. Autosomal dominant germline gain-of-function (GOF) variants in STAT3 result in a PIRD with a broad clinical spectrum. Studies in patients have documented a decreased frequency of FOXP3+ Tregs and an increased frequency of Th17 cells in some patients with active disease. However, the mechanisms of disease pathogenesis in STAT3 GOF syndrome remain largely unknown, and treatment is challenging. We developed a knock-in mouse model harboring a de novo pathogenic human STAT3 variant (p.G421R) and found these mice developed T cell dysregulation, lymphoproliferation, and CD4+ Th1 cell skewing. Surprisingly, Treg numbers, phenotype, and function remained largely intact; however, mice had a selective deficiency in the generation of iTregs. In parallel, we performed single-cell RNA-Seq on T cells from STAT3 GOF patients. We demonstrate only minor changes in the Treg transcriptional signature and an expanded, effector CD8+ T cell population. Together, these findings suggest that Tregs are not the primary driver of disease and highlight the importance of preclinical models in the study of disease mechanisms in rare PIRD.

Cardiac muscle has the highest mitochondrial density of any human tissue, but mitochondrial dysfunction is not a recognized cause of isolated cardiomyopathy. Here, we determined that the rare mitofusin (MFN) 2 R400Q mutation is 15-20× over-represented in clinical cardiomyopathy, whereas this specific mutation is not reported as a cause of MFN2 mutant-induced peripheral neuropathy, Charcot-Marie-Tooth disease type 2A (CMT2A). Accordingly, we interrogated the enzymatic, biophysical, and functional characteristics of MFN2 Q400 versus wild-type and CMT2A-causing MFN2 mutants. All MFN2 mutants had impaired mitochondrial fusion, the canonical MFN2 function. Compared to MFN2 T105M that lacked catalytic GTPase activity and exhibited normal activation-induced changes in conformation, MFN2 R400Q and M376A had normal GTPase activity with impaired conformational shifting. MFN2 R400Q did not suppress mitochondrial motility, provoke mitochondrial depolarization, or dominantly suppress mitochondrial respiration like MFN2 T105M. By contrast to MFN2 T105M and M376A, MFN2 R400Q was uniquely defective in recruiting Parkin to mitochondria. CRISPR editing of the R400Q mutation into the mouse Mfn2 gene induced perinatal cardiomyopathy with no other organ involvement; knock-in of Mfn2 T105M or M376V did not affect the heart. RNA sequencing and metabolomics of cardiomyopathic Mfn2 Q/Q400 hearts revealed signature abnormalities recapitulating experimental mitophagic cardiomyopathy. Indeed, cultured cardiomyoblasts and in vivo cardiomyocytes expressing MFN2 Q400 had mitophagy defects with increased sensitivity to doxorubicin. MFN2 R400Q is the first known natural mitophagy-defective MFN2 mutant. Its unique profile of dysfunction evokes mitophagic cardiomyopathy, suggesting a mechanism for enrichment in clinical cardiomyopathy.