Temporal chromatin accessibility changes define transcriptional states essential for osteosarcoma metastasis.

The metastasis-invasion cascade describes the series of steps required for a cancer cell to successfully spread from its primary tumor and ultimately grow within a secondary organ. Despite metastasis being a dynamic, multistep process, most omics studies to date have focused on comparing primary tumors to the metastatic deposits that define end-stage disease. This static approach means we lack information about the genomic and epigenomic changes that occur during the majority of tumor progression. One particularly understudied phase of tumor progression is metastatic colonization, during which cells must adapt to the new microenvironment of the secondary organ. Through temporal profiling of chromatin accessibility and gene expression in vivo, we identify dynamic changes in the epigenome that occur as osteosarcoma tumors form and grow within the lung microenvironment. Furthermore, we show through paired in vivo and in vitro CRISPR drop-out screens and pharmacological validation that the upstream transcription factors represent a class of metastasis-specific dependency genes. While current models depict lung colonization as a discrete step within the metastatic cascade, our study shows it is a defined trajectory through multiple epigenetic states, revealing new therapeutic opportunities undetectable with standard approaches.

Epigenetic effects of high-fat diet on intestinal tumorigenesis in C57BL/6J-Apc Min/+ mice.

Aim: Obesity and obesogenic diets might partly accelerate cancer development through epigenetic mechanisms. To determine these early effects, we investigated the impact of three days of a high-fat diet on epigenomic and transcriptomic changes in Apc Min/+ murine intestinal epithelia. Method: ChIP-Seq and RNA-Seq were performed on small intestinal epithelia of WT and Apc Min/+ male mice fed high-fat diet (HFD) or low-fat diet (LFD) for three days to identify genomic regions associated with differential H3K27ac levels as a marker of variant enhancer loci (VELs) as well as differentially expressed genes (DEGs). Results: Regarding epigenetic and transcriptomic changes, diet type (LFD vs. HFD) showed a significant impact, and genotype (WT vs.Apc Min/+) showed a small impact. Compared to LFD, HFD resulted in 1306 gained VELs, 230 lost VELs, 133 upregulated genes, and 127 downregulated genes in WT mice, with 1056 gained VELs, 371 lost VELs, 222 upregulated genes, and 182 downregulated genes in Apc Min/+ mice. Compared to the WT genotype, the Apc Min/+ genotype resulted in zero changed VELs for either diet type group, 21 DEGs for LFD, and 48 DEGs for HFD. Most gained VELs, and upregulated genes were associated with lipid metabolic processes. Gained VELs were also associated with Wnt signaling. Downregulated genes were associated with antigen presentation and processing. Conclusion: Three days of HFD-induced epigenomic and transcriptomic changes involving metabolic and immunologic pathways that may promote tumor growth in the genetically predisposed murine intestine without affecting key cancer signaling pathways.

Distinct Cell Adhesion Signature Defines Glioblastoma Myeloid-Derived Suppressor Cell Subsets.

In multiple types of cancer, an increased frequency in myeloid-derived suppressor cells (MDSC) is associated with worse outcomes and poor therapeutic response. In the glioblastoma (GBM) microenvironment, monocytic (m) MDSCs represent the predominant subset. However, the molecular basis of mMDSC enrichment in the tumor microenvironment compared with granulocytic (g) MDSCs has yet to be determined. Here we performed the first broad epigenetic profiling of MDSC subsets to define underlying cell-intrinsic differences in behavior and found that enhanced gene accessibility of cell adhesion programs in mMDSCs is linked to their tumor-accelerating ability in GBM models upon adoptive transfer. Mouse and human mMDSCs expressed higher levels of integrin ß1 and dipeptidyl peptidase-4 (DPP-4) compared with gMDSCs as part of an enhanced cell adhesion signature. Integrin ß1 blockade abrogated the tumor-promoting phenotype of mMDSCs and altered the immune profile in the tumor microenvironment, whereas treatment with a DPP-4 inhibitor extended survival in preclinical GBM models. Targeting DPP-4 in mMDSCs reduced pERK signaling and their migration towards tumor cells. These findings uncover a fundamental difference in the molecular basis of MDSC subsets and suggest that integrin ß1 and DPP-4 represent putative immunotherapy targets to attenuate myeloid cell-driven immune suppression in GBM. SIGNIFICANCE: Epigenetic profiling uncovers cell adhesion programming as a regulator of the tumor-promoting functions of monocytic myeloid-derived suppressor cells in glioblastoma, identifying therapeutic targets that modulate the immune response and suppress tumor growth.

Convergence of case-specific epigenetic alterations identify a confluence of genetic vulnerabilities tied to opioid overdose.

Opioid use disorder is a highly heterogeneous disease driven by a variety of genetic and environmental risk factors which have yet to be fully elucidated. Opioid overdose, the most severe outcome of opioid use disorder, remains the leading cause of accidental death in the United States. We interrogated the effects of opioid overdose on the brain using ChIP-seq to quantify patterns of H3K27 acetylation in dorsolateral prefrontal cortical neurons isolated from 51 opioid-overdose cases and 51 accidental death controls. Among opioid cases, we observed global hypoacetylation and identified 388 putative enhancers consistently depleted for H3K27ac. Machine learning on H3K27ac patterns predicted case-control status with high accuracy. We focused on case-specific regulatory alterations, revealing 81,399 hypoacetylation events, uncovering vast inter-patient heterogeneity. We developed a strategy to decode this heterogeneity based on convergence analysis, which leveraged promoter-capture Hi-C to identify five genes over-burdened by alterations in their regulatory network or "plexus": ASTN2, KCNMA1, DUSP4, GABBR2, ENOX1. These convergent loci are enriched for opioid use disorder risk genes and heritability for generalized anxiety, number of sexual partners, and years of education. Overall, our multi-pronged approach uncovers neurobiological aspects of opioid use disorder and captures genetic and environmental factors perpetuating the opioid epidemic.

Non-canonical Targets of HIF1a Impair Oligodendrocyte Progenitor Cell Function.

Mammalian cells respond to insufficient oxygen through transcriptional regulators called hypoxia-inducible factors (HIFs). Although transiently protective, prolonged HIF activity drives distinct pathological responses in different tissues. Using a model of chronic HIF1a accumulation in pluripotent-stem-cell-derived oligodendrocyte progenitors (OPCs), we demonstrate that HIF1a activates non-canonical targets to impair generation of oligodendrocytes from OPCs. HIF1a activated a unique set of genes in OPCs through interaction with the OPC-specific transcription factor OLIG2. Non-canonical targets, including Ascl2 and Dlx3, were sufficient to block differentiation through suppression of the oligodendrocyte regulator Sox10. Chemical screening revealed that inhibition of MEK/ERK signaling overcame the HIF1a-mediated block in oligodendrocyte generation by restoring Sox10 expression without affecting canonical HIF1a activity. MEK/ERK inhibition also drove oligodendrocyte formation in hypoxic regions of human oligocortical spheroids. This work defines mechanisms by which HIF1a impairs oligodendrocyte formation and establishes that cell-type-specific HIF1a targets perturb cell function in response to low oxygen.

Robust Hi-C Maps of Enhancer-Promoter Interactions Reveal the Function of Non-coding Genome in Neural Development and Diseases.

Genome-wide mapping of chromatin interactions at high resolution remains experimentally and computationally challenging. Here we used a low-input "easy Hi-C" protocol to map the 3D genome architecture in human neurogenesis and brain tissues and also demonstrated that a rigorous Hi-C bias-correction pipeline (HiCorr) can significantly improve the sensitivity and robustness of Hi-C loop identification at sub-TAD level, especially the enhancer-promoter (E-P) interactions. We used HiCorr to compare the high-resolution maps of chromatin interactions from 10 tissue or cell types with a focus on neurogenesis and brain tissues. We found that dynamic chromatin loops are better hallmarks for cellular differentiation than compartment switching. HiCorr allowed direct observation of cell-type- and differentiation-specific E-P aggregates spanning large neighborhoods, suggesting a mechanism that stabilizes enhancer contacts during development. Interestingly, we concluded that Hi-C loop outperforms eQTL in explaining neurological GWAS results, revealing a unique value of high-resolution 3D genome maps in elucidating the disease etiology.

Cell Type-Specific Intralocus Interactions Reveal Oligodendrocyte Mechanisms in MS.

Multiple sclerosis (MS) is an autoimmune disease characterized by attack on oligodendrocytes within the central nervous system (CNS). Despite widespread use of immunomodulatory therapies, patients may still face progressive disability because of failure of myelin regeneration and loss of neurons, suggesting additional cellular pathologies. Here, we describe a general approach for identifying specific cell types in which a disease allele exerts a pathogenic effect. Applying this approach to MS risk loci, we pinpoint likely pathogenic cell types for 70%. In addition to T cell loci, we unexpectedly identified myeloid- and CNS-specific risk loci, including two sites that dysregulate transcriptional pause release in oligodendrocytes. Functional studies demonstrated inhibition of transcriptional elongation is a dominant pathway blocking oligodendrocyte maturation. Furthermore, pause release factors are frequently dysregulated in MS brain tissue. These data implicate cell-intrinsic aberrations outside of the immune system and suggest new avenues for therapeutic development. VIDEO ABSTRACT.

Functional Enhancers Shape Extrachromosomal Oncogene Amplifications.

Non-coding regions amplified beyond oncogene borders have largely been ignored. Using a computational approach, we find signatures of significant co-amplification of non-coding DNA beyond the boundaries of amplified oncogenes across five cancer types. In glioblastoma, EGFR is preferentially co-amplified with its two endogenous enhancer elements active in the cell type of origin. These regulatory elements, their contacts, and their contribution to cell fitness are preserved on high-level circular extrachromosomal DNA amplifications. Interrogating the locus with a CRISPR interference screening approach reveals a diversity of additional elements that impact cell fitness. The pattern of fitness dependencies mirrors the rearrangement of regulatory elements and accompanying rewiring of the chromatin topology on the extrachromosomal amplicon. Our studies indicate that oncogene amplifications are shaped by regulatory dependencies in the non-coding genome.

Mismatch repair-signature mutations activate gene enhancers across human colorectal cancer epigenomes.

Commonly-mutated genes have been found for many cancers, but less is known about mutations in cis-regulatory elements. We leverage gains in tumor-specific enhancer activity, coupled with allele-biased mutation detection from H3K27ac ChIP-seq data, to pinpoint potential enhancer-activating mutations in colorectal cancer (CRC). Analysis of a genetically-diverse cohort of CRC specimens revealed that microsatellite instable (MSI) samples have a high indel rate within active enhancers. Enhancers with indels show evidence of positive selection, increased target gene expression, and a subset is highly recurrent. The indels affect short homopolymer tracts of A/T and increase affinity for FOX transcription factors. We further demonstrate that signature mismatch-repair (MMR) mutations activate enhancers using a xenograft tumor metastasis model, where mutations are induced naturally via CRISPR/Cas9 inactivation of MLH1 prior to tumor cell injection. Our results suggest that MMR signature mutations activate enhancers in CRC tumor epigenomes to provide a selective advantage.

Corrigendum: Positively selected enhancer elements endow osteosarcoma cells with metastatic competence.

This corrects the article DOI: 10.1038/nm.4475.

CHD7 represses the retinoic acid synthesis enzyme ALDH1A3 during inner ear development.

CHD7, an ATP-dependent chromatin remodeler, is disrupted in CHARGE syndrome, an autosomal dominant disorder characterized by variably penetrant abnormalities in craniofacial, cardiac, and nervous system tissues. The inner ear is uniquely sensitive to CHD7 levels and is the most commonly affected organ in individuals with CHARGE. Interestingly, upregulation or downregulation of retinoic acid (RA) signaling during embryogenesis also leads to developmental defects similar to those in CHARGE syndrome, suggesting that CHD7 and RA may have common target genes or signaling pathways. Here, we tested three separate potential mechanisms for CHD7 and RA interaction: (a) direct binding of CHD7 with RA receptors, (b) regulation of CHD7 levels by RA, and (c) CHD7 binding and regulation of RA-related genes. We show that CHD7 directly regulates expression of Aldh1a3, the gene encoding the RA synthetic enzyme ALDH1A3 and that loss of Aldh1a3 partially rescues Chd7 mutant mouse inner ear defects. Together, these studies indicate that ALDH1A3 acts with CHD7 in a common genetic pathway to regulate inner ear development, providing insights into how CHD7 and RA regulate gene expression and morphogenesis in the developing embryo.

Positively selected enhancer elements endow osteosarcoma cells with metastatic competence.

Metastasis results from a complex set of traits acquired by tumor cells, distinct from those necessary for tumorigenesis. Here, we investigate the contribution of enhancer elements to the metastatic phenotype of osteosarcoma. Through epigenomic profiling, we identify substantial differences in enhancer activity between primary and metastatic human tumors and between near isogenic pairs of highly lung metastatic and nonmetastatic osteosarcoma cell lines. We term these regions metastatic variant enhancer loci (Met-VELs). Met-VELs drive coordinated waves of gene expression during metastatic colonization of the lung. Met-VELs cluster nonrandomly in the genome, indicating that activity of these enhancers and expression of their associated gene targets are positively selected. As evidence of this causal association, osteosarcoma lung metastasis is inhibited by global interruptions of Met-VEL-associated gene expression via pharmacologic BET inhibition, by knockdown of AP-1 transcription factors that occupy Met-VELs, and by knockdown or functional inhibition of individual genes activated by Met-VELs, such as that encoding coagulation factor III/tissue factor (F3). We further show that genetic deletion of a single Met-VEL at the F3 locus blocks metastatic cell outgrowth in the lung. These findings indicate that Met-VELs and the genes they regulate play a functional role in metastasis and may be suitable targets for antimetastatic therapies.

Hotspots of aberrant enhancer activity punctuate the colorectal cancer epigenome.

In addition to mutations in genes, aberrant enhancer element activity at non-coding regions of the genome is a key driver of tumorigenesis. Here, we perform epigenomic enhancer profiling of a cohort of more than forty genetically diverse human colorectal cancer (CRC) specimens. Using normal colonic crypt epithelium as a comparator, we identify enhancers with recurrently gained or lost activity across CRC specimens. Of the enhancers highly recurrently activated in CRC, most are constituents of super enhancers, are occupied by AP-1 and cohesin complex members, and originate from primed chromatin. Many activate known oncogenes, and CRC growth can be mitigated through pharmacologic inhibition or genome editing of these loci. Nearly half of all GWAS CRC risk loci co-localize to recurrently activated enhancers. These findings indicate that the CRC epigenome is defined by highly recurrent epigenetic alterations at enhancers which activate a common, aberrant transcriptional programme critical for CRC growth and survival.

Transcriptome-wide identification of mRNAs and lincRNAs associated with trastuzumab-resistance in HER2-positive breast cancer.

Approximately, 25-30% of early-stage breast tumors are classified at the molecular level as HER2-positive, which is an aggressive subtype of breast cancer. Amplification of the HER2 gene in these tumors results in a substantial increase in HER2 mRNA levels, and consequently, HER2 protein levels. HER2, a transmembrane receptor tyrosine kinase (RTK), is targeted therapeutically by a monoclonal antibody, trastuzumab (Tz), which has dramatically improved the prognosis of HER2-driven breast cancers. However, ~30% of patients develop resistance to trastuzumab and recur; and nearly all patients with advanced disease develop resistance over time and succumb to the disease. Mechanisms of trastuzumab resistance (TzR) are not well understood, although some studies suggest that growth factor signaling through other receptors may be responsible. However, these studies were based on cell culture models of the disease, and thus, it is not known which pathways are driving the resistance in vivo. Using an integrative transcriptomic approach of RNA isolated from trastuzumab-sensitive and trastuzumab-resistant HER2+ tumors, and isogenic cell culture models, we identified a small set of mRNAs and lincRNAs that are associated with trastuzumab-resistance (TzR). Functional analysis of a top candidate gene, S100P, demonstrated that inhibition of S100P results in reversing TzR. Mechanistically, S100P activates the RAS/MEK/MAPK pathway to compensate for HER2 inhibition by trastuzumab. Finally, we demonstrated that the upregulation of S100P appears to be driven by epigenomic changes at the enhancer level. Our current findings should pave the path toward new therapies for breast cancer patients.

Protective effects of C-type natriuretic peptide on linear growth and articular cartilage integrity in a mouse model of inflammatory arthritis.

OBJECTIVE: The C-type natriuretic peptide (CNP) signaling pathway is a major contributor to postnatal skeletal growth in humans. This study was undertaken to investigate whether CNP signaling could prevent growth delay and cartilage damage in an animal model of inflammatory arthritis. METHODS: We generated transgenic mice that overexpress CNP (B6.SJL-Col2a1-NPPC) in chondrocytes. We introduced the CNP transgene into mice with experimental systemic inflammatory arthritis (K/BxN T cell receptor [TCR]) and determined the effect of CNP overexpression in chondrocytes on the severity of arthritis, cartilage damage, and linear growth. We also examined primary chondrocyte cultures for changes in gene and protein expression resulting from CNP overexpression. RESULTS: K/BxN TCR mice exhibited linear growth delay (P < 0.01) compared to controls, and this growth delay was correlated with the severity of arthritis. Diminished chondrocyte proliferation and matrix production was also seen in K/BxN TCR mice. Compared to non-CNP-transgenic mice, K/BxN TCR mice with overexpressed CNP had milder arthritis, no growth delay, and less cartilage damage. Primary chondrocytes from mice overexpressing CNP were less sensitive to inflammatory cytokines than wild-type mouse chondrocytes. CONCLUSION: CNP overexpression in chondrocytes can prevent endochondral growth delay and protect against cartilage damage in a mouse model of inflammatory arthritis. Pharmacologic or biologic modulation of the CNP signaling pathway may prevent growth retardation and protect cartilage in patients with inflammatory joint diseases, such as juvenile idiopathic arthritis.

The diagnostic challenge of progressive pseudorheumatoid dysplasia (PPRD): a review of clinical features, radiographic features, and WISP3 mutations in 63 affected individuals.

Progressive pseudorheumatoid dysplasia (PPRD) is a genetic, non-inflammatory arthropathy caused by recessive loss of function mutations in WISP3 (Wnt1-inducible signaling pathway protein 3; MIM 603400), encoding for a signaling protein. The disease is clinically silent at birth and in infancy. It manifests between the age of 3 and 6 years with joint pain and progressive joint stiffness. Affected children are referred to pediatric rheumatologists and orthopedic surgeons; however, signs of inflammation are absent and anti-inflammatory treatment is of little help. Bony enlargement at the interphalangeal joints progresses leading to camptodactyly. Spine involvement develops in late childhood and adolescence leading to short trunk with thoracolumbar kyphosis. Adult height is usually below the 3rd percentile. Radiographic signs are relatively mild. Platyspondyly develops in late childhood and can be the first clue to the diagnosis. Enlargement of the phalangeal metaphyses develops subtly and is usually recognizable by 10 years. The femoral heads are large and the acetabulum forms a distinct "lip" overriding the femoral head. There is a progressive narrowing of all articular spaces as articular cartilage is lost. Medical management of PPRD remains symptomatic and relies on pain medication. Hip joint replacement surgery in early adulthood is effective in reducing pain and maintaining mobility and can be recommended. Subsequent knee joint replacement is a further option. Mutation analysis of WISP3 allowed the confirmation of the diagnosis in 63 out of 64 typical cases in our series. Intronic mutations in WISP3 leading to splicing aberrations can be detected only in cDNA from fibroblasts and therefore a skin biopsy is indicated when genomic analysis fails to reveal mutations in individuals with otherwise typical signs and symptoms. In spite of the first symptoms appearing in early childhood, the diagnosis of PPRD is most often made only in the second decade and affected children often receive unnecessary anti-inflammatory and immunosuppressive treatments. Increasing awareness of PPRD appears to be essential to allow for a timely diagnosis.

H3K4me3 inversely correlates with DNA methylation at a large class of non-CpG-island-containing start sites.

BACKGROUND: In addition to mutations, epigenetic silencing of genes has been recognized as a fundamental mechanism that promotes human carcinogenesis. To date, characterization of epigenetic gene silencing has largely focused on genes in which silencing is mediated by hypermethylation of promoter-associated CpG islands, associated with loss of the H3K4me3 chromatin mark. Far less is known about promoters lacking CpG-islands or genes that are repressed by alternative mechanisms. METHODS: We performed integrative ChIP-chip, DNase-seq, and global gene expression analyses in colon cancer cells and normal colon mucosa to characterize chromatin features of both CpG-rich and CpG-poor promoters of genes that undergo silencing in colon cancer. RESULTS: Epigenetically repressed genes in colon cancer separate into two classes based on retention or loss of H3K4me3 at transcription start sites. Quantitatively, of transcriptionally repressed genes that lose H3K4me3 in colon cancer (K4-dependent genes), a large fraction actually lacks CpG islands. Nonetheless, similar to CpG-island containing genes, cytosines located near the start sites of K4-dependent genes become DNA hypermethylated, and repressed K4-dependent genes can be reactivated with 5-azacytidine. Moreover, we also show that when the H3K4me3 mark is retained, silencing of CpG island-associated genes can proceed through an alternative mechanism in which repressive chromatin marks are recruited. CONCLUSIONS: H3K4me3 equally protects from DNA methylation at both CpG-island and non-CpG island start sites in colon cancer. Moreover, the results suggest that CpG-rich genes repressed by loss of H3K4me3 and DNA methylation represent special instances of a more general epigenetic mechanism of gene silencing, one in which gene silencing is mediated by loss of H3K4me3 and methylation of non-CpG island promoter-associated cytosines.

Epigenomic enhancer profiling defines a signature of colon cancer.

Cancer is characterized by gene expression aberrations. Studies have largely focused on coding sequences and promoters, even though distal regulatory elements play a central role in controlling transcription patterns. We used the histone mark H3K4me1 to analyze gain and loss of enhancer activity genome-wide in primary colon cancer lines relative to normal colon crypts. We identified thousands of variant enhancer loci (VELs) that comprise a signature that is robustly predictive of the in vivo colon cancer transcriptome. Furthermore, VELs are enriched in haplotype blocks containing colon cancer genetic risk variants, implicating these genomic regions in colon cancer pathogenesis. We propose that reproducible changes in the epigenome at enhancer elements drive a specific transcriptional program to promote colon carcinogenesis.

Mutations in LRP5 cause primary osteoporosis without features of OI by reducing Wnt signaling activity.

BACKGROUND: Primary osteoporosis is a rare childhood-onset skeletal condition whose pathogenesis has been largely unknown. We have previously shown that primary osteoporosis can be caused by heterozygous missense mutations in the Low-density lipoprotein receptor-related protein 5 (LRP5) gene, and the role of LRP5 is further investigated here. METHODS: LRP5 was analyzed in 18 otherwise healthy children and adolescents who had evidence of osteoporosis (manifested as reduced bone mineral density i.e. BMD, recurrent peripheral fractures and/or vertebral compression fractures) but who lacked the clinical features of osteogenesis imperfecta (OI) or other known syndromes linked to low BMD. Also 51 controls were analyzed. Methods used in the genetic analyses included direct sequencing and multiplex ligation-dependent probe amplification (MLPA). In vitro studies were performed using luciferase assay and quantitative real-time polymerase chain reaction (qPCR) to examine the effect of two novel and three previously identified mutations on the activity of canonical Wnt signaling and on expression of tryptophan hydroxylase 1 (Tph1) and 5-hydroxytryptamine (5-Htr1b). RESULTS: Two novel LRP5 mutations (c.3446 T > A; p.L1149Q and c.3553 G > A; p.G1185R) were identified in two patients and their affected family members. In vitro analyses showed that one of these novel mutations together with two previously reported mutations (p.C913fs, p.R1036Q) significantly reduced the activity of the canonical Wnt signaling pathway. Such reductions may lead to decreased bone formation, and could explain the bone phenotype. Gut-derived Lrp5 has been shown to regulate serotonin synthesis by controlling the production of serotonin rate-limiting enzyme, Tph1. LRP5 mutations did not affect Tph1 expression, and only one mutant (p.L1149Q) reduced expression of serotonin receptor 5-Htr1b (p < 0.002). CONCLUSIONS: Our results provide additional information on the role of LRP5 mutations and their effects on the development of juvenile-onset primary osteoporosis, and hence the pathogenesis of the disorder. The mutations causing primary osteoporosis reduce the signaling activity of the canonical Wnt signaling pathway and may therefore result in decreased bone formation. The specific mechanism affecting signaling activity remains to be resolved in future studies.