Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

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.

Probing the diabetes and colorectal cancer relationship using gene - environment interaction analyses

BACKGROUND

Diabetes is an established risk factor for colorectal cancer. However, the mechanisms underlying this relationship still require investigation and it is not known if the association is modified by genetic variants. To address these questions, we undertook a genome-wide gene-environment interaction analysis.

METHODS

We used data from 3 genetic consortia (CCFR, CORECT, GECCO; 31,318 colorectal cancer cases/41,499 controls) and undertook genome-wide gene-environment interaction analyses with colorectal cancer risk, including interaction tests of genetics(G)xdiabetes (1-degree of freedom; d.f.) and joint testing of Gxdiabetes, G-colorectal cancer association (2-d.f. joint test) and G-diabetes correlation (3-d.f. joint test).

RESULTS

Based on the joint tests, we found that the association of diabetes with colorectal cancer risk is modified by loci on chromosomes 8q24.11 (rs3802177, SLC30A8 - ORAA: 1.62, 95% CI: 1.34-1.96; ORAG: 1.41, 95% CI: 1.30-1.54; ORGG: 1.22, 95% CI: 1.13-1.31; p-value3-d.f.: 5.46 × 10-11) and 13q14.13 (rs9526201, LRCH1 - ORGG: 2.11, 95% CI: 1.56-2.83; ORGA: 1.52, 95% CI: 1.38-1.68; ORAA: 1.13, 95% CI: 1.06-1.21; p-value2-d.f.: 7.84 × 10-09).

DISCUSSION

These results suggest that variation in genes related to insulin signaling (SLC30A8) and immune function (LRCH1) may modify the association of diabetes with colorectal cancer risk and provide novel insights into the biology underlying the diabetes and colorectal cancer relationship.

A Genetic Locus within the FMN1/GREM1 Gene Region Interacts with Body Mass Index in Colorectal Cancer Risk

Colorectal cancer risk can be impacted by genetic, environmental, and lifestyle factors, including diet and obesity. Gene-environment interactions (G × E) can provide biological insights into the effects of obesity on colorectal cancer risk. Here, we assessed potential genome-wide G × E interactions between body mass index (BMI) and common SNPs for colorectal cancer risk using data from 36,415 colorectal cancer cases and 48,451 controls from three international colorectal cancer consortia (CCFR, CORECT, and GECCO). The G × E tests included the conventional logistic regression using multiplicative terms (one degree of freedom, 1DF test), the two-step EDGE method, and the joint 3DF test, each of which is powerful for detecting G × E interactions under specific conditions. BMI was associated with higher colorectal cancer risk. The two-step approach revealed a statistically significant G×BMI interaction located within the Formin 1/Gremlin 1 (FMN1/GREM1) gene region (rs58349661). This SNP was also identified by the 3DF test, with a suggestive statistical significance in the 1DF test. Among participants with the CC genotype of rs58349661, overweight and obesity categories were associated with higher colorectal cancer risk, whereas null associations were observed across BMI categories in those with the TT genotype. Using data from three large international consortia, this study discovered a locus in the FMN1/GREM1 gene region that interacts with BMI on the association with colorectal cancer risk. Further studies should examine the potential mechanisms through which this locus modifies the etiologic link between obesity and colorectal cancer.

SIGNIFICANCE

This gene-environment interaction analysis revealed a genetic locus in FMN1/GREM1 that interacts with body mass index in colorectal cancer risk, suggesting potential implications for precision prevention strategies.

Genome-wide Interaction Study with Smoking for Colorectal Cancer Risk Identifies Novel Genetic Loci Related to Tumor Suppression, Inflammation, and Immune Response

BACKGROUND

Tobacco smoking is an established risk factor for colorectal cancer. However, genetically defined population subgroups may have increased susceptibility to smoking-related effects on colorectal cancer.

METHODS

A genome-wide interaction scan was performed including 33,756 colorectal cancer cases and 44,346 controls from three genetic consortia.

RESULTS

Evidence of an interaction was observed between smoking status (ever vs. never smokers) and a locus on 3p12.1 (rs9880919, P = 4.58 × 10-8), with higher associated risk in subjects carrying the GG genotype [OR, 1.25; 95% confidence interval (CI), 1.20-1.30] compared with the other genotypes (OR <1.17 for GA and AA). Among ever smokers, we observed interactions between smoking intensity (increase in 10 cigarettes smoked per day) and two loci on 6p21.33 (rs4151657, P = 1.72 × 10-8) and 8q24.23 (rs7005722, P = 2.88 × 10-8). Subjects carrying the rs4151657 TT genotype showed higher risk (OR, 1.12; 95% CI, 1.09-1.16) compared with the other genotypes (OR <1.06 for TC and CC). Similarly, higher risk was observed among subjects carrying the rs7005722 AA genotype (OR, 1.17; 95% CI, 1.07-1.28) compared with the other genotypes (OR <1.13 for AC and CC). Functional annotation revealed that SNPs in 3p12.1 and 6p21.33 loci were located in regulatory regions, and were associated with expression levels of nearby genes. Genetic models predicting gene expression revealed that smoking parameters were associated with lower colorectal cancer risk with higher expression levels of CADM2 (3p12.1) and ATF6B (6p21.33).

CONCLUSIONS

Our study identified novel genetic loci that may modulate the risk for colorectal cancer of smoking status and intensity, linked to tumor suppression and immune response.

IMPACT

These findings can guide potential prevention treatments.

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.

Abstract IA019: How to target lung metastasis: dynamic changes in chromatin accessibility during osteosarcoma lung metastasis expose vulnerabilities

The metastatic cascade describes the series of steps required for a cancer cell to successfully spread from its primary tumor to a secondary organ. Despite knowing that metastasis is a fluid process, almost all studies to date have relied on end-point assays. As a result, little is known about the continuum of molecular events required for successful growth within the secondary organ. By profiling accessible chromatin in osteosarcoma cells grown in vitro and isolated from an early and late timepoint during lung metastasis, we identify temporally distinct epigenomic reprogramming events that occur as tumors form and grow within the metastatic microenvironment. Furthermore, these data pinpoint key upstream transcription factors responsible for the dynamic chromatin changes. Through paired in vivo and in vitro CRISPR drop-out screens, we show a subset of these putative pro-metastatic transcription factors are essential for lung metastasis, but dispensable for in vitro growth. Targeting the same factors with chemical probes in vitro and ex vivo demonstrated similar context-dependent essentiality. Altogether, our study demonstrates the epigenomes of metastasizing cancer cells are tightly regulated in a dynamic and context-specific manner, and the transcription factors controlling this process may serve as novel vulnerabilities.

Citation Format: Peter C. Scacheri. How to target lung metastasis: dynamic changes in chromatin accessibility during osteosarcoma lung metastasis expose vulnerabilities [abstract]. In: Proceedings of the AACR Special Conference: Sarcomas; 2022 May 9-12; Montreal, QC, Canada. Philadelphia (PA): AACR; Clin Cancer Res 2022;28(18_Suppl):Abstract nr IA019.

Beyond GWAS of Colorectal Cancer: Evidence of Interaction with Alcohol Consumption and Putative Causal Variant for the 10q24.2 Region

BACKGROUND

Currently known associations between common genetic variants and colorectal cancer explain less than half of its heritability of 25%. As alcohol consumption has a J-shape association with colorectal cancer risk, nondrinking and heavy drinking are both risk factors for colorectal cancer.

METHODS

Individual-level data was pooled from the Colon Cancer Family Registry, Colorectal Transdisciplinary Study, and Genetics and Epidemiology of Colorectal Cancer Consortium to compare nondrinkers (≤1 g/day) and heavy drinkers (>28 g/day) with light-to-moderate drinkers (1-28 g/day) in GxE analyses. To improve power, we implemented joint 2df and 3df tests and a novel two-step method that modifies the weighted hypothesis testing framework. We prioritized putative causal variants by predicting allelic effects using support vector machine models.

RESULTS

For nondrinking as compared with light-to-moderate drinking, the hybrid two-step approach identified 13 significant SNPs with pairwise r2 > 0.9 in the 10q24.2/COX15 region. When stratified by alcohol intake, the A allele of lead SNP rs2300985 has a dose-response increase in risk of colorectal cancer as compared with the G allele in light-to-moderate drinkers [OR for GA genotype = 1.11; 95% confidence interval (CI), 1.06-1.17; OR for AA genotype = 1.22; 95% CI, 1.14-1.31], but not in nondrinkers or heavy drinkers. Among the correlated candidate SNPs in the 10q24.2/COX15 region, rs1318920 was predicted to disrupt an HNF4 transcription factor binding motif.

CONCLUSIONS

Our study suggests that the association with colorectal cancer in 10q24.2/COX15 observed in genome-wide association study is strongest in nondrinkers. We also identified rs1318920 as the putative causal regulatory variant for the region.

IMPACT

The study identifies multifaceted evidence of a possible functional effect for rs1318920.

Chromatin Mechanisms Driving Cancer

The change in cell state from normal to malignant is driven fundamentally by oncogenic mutations in cooperation with epigenetic alterations of chromatin. These alterations in chromatin can be a consequence of environmental stressors or germline and/or somatic mutations that directly alter the structure of chromatin machinery proteins, their levels, or their regulatory function. These changes can result in an inability of the cell to differentiate along a predefined lineage path, or drive a hyperactive, highly proliferative state with addiction to high levels of transcriptional output. We discuss how these genetic alterations hijack the chromatin machinery for the oncogenic process to reveal unique vulnerabilities and novel targets for cancer therapy.

Genetic architectures of proximal and distal colorectal cancer are partly distinct

OBJECTIVE

An understanding of the etiologic heterogeneity of colorectal cancer (CRC) is critical for improving precision prevention, including individualized screening recommendations and the discovery of novel drug targets and repurposable drug candidates for chemoprevention. Known differences in molecular characteristics and environmental risk factors among tumors arising in different locations of the colorectum suggest partly distinct mechanisms of carcinogenesis. The extent to which the contribution of inherited genetic risk factors for CRC differs by anatomical subsite of the primary tumor has not been examined.

DESIGN

To identify new anatomical subsite-specific risk loci, we performed genome-wide association study (GWAS) meta-analyses including data of 48 214 CRC cases and 64 159 controls of European ancestry. We characterised effect heterogeneity at CRC risk loci using multinomial modelling.

RESULTS

We identified 13 loci that reached genome-wide significance (p<5×10-8) and that were not reported by previous GWASs for overall CRC risk. Multiple lines of evidence support candidate genes at several of these loci. We detected substantial heterogeneity between anatomical subsites. Just over half (61) of 109 known and new risk variants showed no evidence for heterogeneity. In contrast, 22 variants showed association with distal CRC (including rectal cancer), but no evidence for association or an attenuated association with proximal CRC. For two loci, there was strong evidence for effects confined to proximal colon cancer.

CONCLUSION

Genetic architectures of proximal and distal CRC are partly distinct. Studies of risk factors and mechanisms of carcinogenesis, and precision prevention strategies should take into consideration the anatomical subsite of the tumour.

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.

Ex vivo screen identifies CDK12 as a metastatic vulnerability in osteosarcoma

Despite progress in intensification of therapy, outcomes for patients with metastatic osteosarcoma (OS) have not improved in thirty years. We developed a system that enabled preclinical screening of compounds against metastatic OS cells in the context of the native lung microenvironment. Using this strategy to screen a library of epigenetically targeted compounds, we identified inhibitors of CDK12 to be most effective, reducing OS cell outgrowth in the lung by more than 90% at submicromolar doses. We found that knockout of CDK12 in an in vivo model of lung metastasis significantly decreased the ability of OS to colonize the lung. CDK12 inhibition led to defects in transcription elongation in a gene length- and expression-dependent manner. These effects were accompanied by defects in RNA processing and altered the expression of genes involved in transcription regulation and the DNA damage response. We further identified OS models that differ in their sensitivity to CDK12 inhibition in the lung and provided evidence that upregulated MYC levels may mediate these differences. Our studies provided a framework for rapid preclinical testing of compounds with antimetastatic activity and highlighted CDK12 as a potential therapeutic target in OS.

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.

Chemotherapy-Induced Distal Enhancers Drive Transcriptional Programs to Maintain the Chemoresistant State in Ovarian Cancer

Chemoresistance is driven by unique regulatory networks in the genome that are distinct from those necessary for cancer development. Here, we investigate the contribution of enhancer elements to cisplatin resistance in ovarian cancers. Epigenome profiling of multiple cellular models of chemoresistance identified unique sets of distal enhancers, super-enhancers (SE), and their gene targets that coordinate and maintain the transcriptional program of the platinum-resistant state in ovarian cancer. Pharmacologic inhibition of distal enhancers through small-molecule epigenetic inhibitors suppressed the expression of their target genes and restored cisplatin sensitivity in vitro and in vivo. In addition to known drivers of chemoresistance, our findings identified SOX9 as a critical SE-regulated transcription factor that plays a critical role in acquiring and maintaining the chemoresistant state in ovarian cancer. The approach and findings presented here suggest that integrative analysis of epigenome and transcriptional programs could identify targetable key drivers of chemoresistance in cancers. SIGNIFICANCE: Integrative genome-wide epigenomic and transcriptomic analyses of platinum-sensitive and -resistant ovarian lines identify key distal regulatory regions and associated master regulator transcription factors that can be targeted by small-molecule epigenetic inhibitors.

Correction to: Genetic variant predictors of gene expression provide new insight into risk of colorectal cancer

Every author has erroneously been assigned to the affiliation "62". The affiliation 62 belongs to the author Graham Casey.

Enhancers: bridging the gap between gene control and human disease

Enhancers are a class of regulatory elements essential for precise spatio-temporal control of gene expression during development and in terminally differentiated cells. This review highlights signature features of enhancer elements as well as new advances that provide mechanistic insights into enhancer-mediated gene control in the context of three-dimensional chromatin. We detail the various ways in which non-coding mutations can instigate aberrant gene control and cause a variety of Mendelian disorders, common diseases and cancer.

Chromatin landscapes reveal developmentally encoded transcriptional states that define human glioblastoma

Mack et al. defined active chromatin landscapes of glioblastoma stem cells (GSCs) and primary tumor specimens, revealing novel transcriptional regulatory circuits and therapeutic targets. Super-enhancers identified essential transcription factors that underlie GSC identity and intertumoral diversity, potentially informing precision medicine.

Glioblastoma is an incurable brain cancer characterized by high genetic and pathological heterogeneity. Here, we mapped active chromatin landscapes with gene expression, whole exomes, copy number profiles, and DNA methylomes across 44 patient-derived glioblastoma stem cells (GSCs), 50 primary tumors, and 10 neural stem cells (NSCs) to identify essential super-enhancer (SE)–associated genes and the core transcription factors that establish SEs and maintain GSC identity. GSCs segregate into two groups dominated by distinct enhancer profiles and unique developmental core transcription factor regulatory programs. Group-specific transcription factors enforce GSC identity; they exhibit higher activity in glioblastomas versus NSCs, are associated with poor clinical outcomes, and are required for glioblastoma growth in vivo. Although transcription factors are commonly considered undruggable, group-specific enhancer regulation of the MAPK/ERK pathway predicts sensitivity to MEK inhibition. These data demonstrate that transcriptional identity can be leveraged to identify novel dependencies and therapeutic approaches.

Genetic variant predictors of gene expression provide new insight into risk of colorectal cancer

Genome-wide association studies have reported 56 independently associated colorectal cancer (CRC) risk variants, most of which are non-coding and believed to exert their effects by modulating gene expression. The computational method PrediXcan uses cis-regulatory variant predictors to impute expression and perform gene-level association tests in GWAS without directly measured transcriptomes. In this study, we used reference datasets from colon (n = 169) and whole blood (n = 922) transcriptomes to test CRC association with genetically determined expression levels in a genome-wide analysis of 12,186 cases and 14,718 controls. Three novel associations were discovered from colon transverse models at FDR ≤ 0.2 and further evaluated in an independent replication including 32,825 cases and 39,933 controls. After adjusting for multiple comparisons, we found statistically significant associations using colon transcriptome models with TRIM4 (discovery P = 2.2 × 10- 4, replication P = 0.01), and PYGL (discovery P = 2.3 × 10- 4, replication P = 6.7 × 10- 4). Interestingly, both genes encode proteins that influence redox homeostasis and are related to cellular metabolic reprogramming in tumors, implicating a novel CRC pathway linked to cell growth and proliferation. Defining CRC risk regions as one megabase up- and downstream of one of the 56 independent risk variants, we defined 44 non-overlapping CRC-risk regions. Among these risk regions, we identified genes associated with CRC (P < 0.05) in 34/44 CRC-risk regions. Importantly, CRC association was found for two genes in the previously reported 2q25 locus, CXCR1 and CXCR2, which are potential cancer therapeutic targets. These findings provide strong candidate genes to prioritize for subsequent laboratory follow-up of GWAS loci. This study is the first to implement PrediXcan in a large colorectal cancer study and findings highlight the utility of integrating transcriptome data in GWAS for discovery of, and biological insight into, risk loci.

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.

Discovery of common and rare genetic risk variants for colorectal cancer

To further dissect the genetic architecture of colorectal cancer (CRC), we performed whole-genome sequencing of 1,439 cases and 720 controls, imputed discovered sequence variants and Haplotype Reference Consortium panel variants into genome-wide association study data, and tested for association in 34,869 cases and 29,051 controls. Findings were followed up in an additional 23,262 cases and 38,296 controls. We discovered a strongly protective 0.3% frequency variant signal at CHD1. In a combined meta-analysis of 125,478 individuals, we identified 40 new independent signals at P < 5 × 10-8, bringing the number of known independent signals for CRC to ~100. New signals implicate lower-frequency variants, Krüppel-like factors, Hedgehog signaling, Hippo-YAP signaling, long noncoding RNAs and somatic drivers, and support a role for immune function. Heritability analyses suggest that CRC risk is highly polygenic, and larger, more comprehensive studies enabling rare variant analysis will improve understanding of biology underlying this risk and influence personalized screening strategies and drug development.

Targeting Epigenetics to Prevent Obesity Promoted Cancers

Epigenetic changes in DNA and associated chromatin proteins are increasingly being considered as important mediators of the linkage between obesity and cancer. Although multiple agents, targeted at epigenetic changes, are being tested for therapy of established cancers, this issue of Cancer Prevention Research carries two articles demonstrating that the bromodomain inhibitor I-BET-762 can attenuate adipose tissue-promoted cancers. Although I-BET-762 significantly delayed, rather than completely prevented, the onset of adiposity-promoted transformation and malignancy, these experiments provide important proof of principle for the strategies of targeting epigenetic changes to disrupt the obesity-cancer linkage. Because bromodomain proteins represent only one of multiple epigenetic mediators, it is probable that targeting other epigenetic processes, alone or in combination, may serve to even more effectively disrupt the obesity promotion of cancer. Given the magnitude of the current obesity pandemic and its impact on cancer, preventive measures to disrupt this linkage are critically important. Cancer Prev Res; 11(3); 125-8. ©2018 AACRSee related article by Chakraborty et al., p. 129.

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.

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.

Therapeutic targeting of ependymoma as informed by oncogenic enhancer profiling

Genomic sequencing has driven precision-based oncology therapy; however, genetic drivers remain unknown or non-targetable for many malignancies, demanding alternative approaches to identify therapeutic leads. Ependymomas are chemotherapy-resistant brain tumours, which, despite genomic sequencing, lack effective molecular targets. Intracranial ependymomas are segregated based on anatomical location – supratentorial region (ST) or posterior fossa (PF) – and further divided into distinct molecular subgroups that reflect differences in age of onset, gender predominance, and response to therapy^1–3. The most common and aggressive subgroup, Posterior Fossa Ependymoma Group A (PF-EPN-A), occurs in young children and appears to lack recurrent somatic mutations². Conversely, Posterior Fossa Ependymoma Group B (PF-EPN-B) tumours display frequent large-scale copy number gains and losses yet favourable clinical outcomes^1,3. Greater than 70% of supratentorial ependymomas are defined by highly recurrent gene fusions in the NFκB subunit RELA (ST-EPN-RELA), and less frequently involve fusion of the gene encoding the transcriptional activator YAP1 (ST-EPN-YAP1).^1,3,4 Subependymomas, a distinct histologic variant, can also be found within the ST and PF compartments accounting for the majority of tumours in the molecular subgroups ST-EPN-SE and PF-EPN-SE, respectively¹. Here, we mapped active chromatin landscapes in 42 primary ependymomas in two non-overlapping primary ependymoma cohorts with the goal of identifying essential super enhancer associated genes on which tumour cells were dependent. Enhancer regions revealed putative oncogenes, molecular targets, and pathways, which when subjected to small molecule inhibitor or shRNA treatment, diminished proliferation of patient-derived neurospheres and increased survival in mouse models of ependymomas. Through profiling of transcriptional enhancers, our study provides a framework for target and drug discovery in other cancers recalcitrant to therapeutic development because of their lack of known genetic drivers.

Mitotic Vulnerability in Triple-Negative Breast Cancer Associated with LIN9 Is Targetable with BET Inhibitors

Triple-negative breast cancers (TNBC) are highly aggressive, lack FDA-approved targeted therapies, and frequently recur, making the discovery of novel therapeutic targets for this disease imperative. Our previous analysis of the molecular mechanisms of action of bromodomain and extraterminal protein inhibitors (BETi) in TNBC revealed these drugs cause multinucleation, indicating BET proteins are essential for efficient mitosis and cytokinesis. Here, using live cell imaging, we show that BET inhibition prolonged mitotic progression and induced mitotic cell death, both of which are indicative of mitotic catastrophe. Mechanistically, the mitosis regulator LIN9 was a direct target of BET proteins that mediated the effects of BET proteins on mitosis in TNBC. Although BETi have been proposed to function by dismantling super-enhancers (SE), the LIN9 gene lacks an SE but was amplified or overexpressed in the majority of TNBCs. In addition, its mRNA expression predicted poor outcome across breast cancer subtypes. Together, these results provide a mechanism for cancer selectivity of BETi that extends beyond modulation of SE-associated genes and suggest that cancers dependent upon LIN9 overexpression may be particularly vulnerable to BETi. Cancer Res; 77(19); 5395-408. ©2017 AACR.

Transcription elongation factors represent in vivo cancer dependencies in glioblastoma

Glioblastoma is a universally lethal cancer with a median survival time of approximately 15 months. Despite substantial efforts to define druggable targets, there are no therapeutic options that notably extend the lifespan of patients with glioblastoma. While previous work has largely focused on in vitro cellular models, here we demonstrate a more physiologically relevant approach to target discovery in glioblastoma. We adapted pooled RNA interference (RNAi) screening technology for use in orthotopic patient-derived xenograft models, creating a high-throughput negative-selection screening platform in a functional in vivo tumour microenvironment. Using this approach, we performed parallel in vivo and in vitro screens and discovered that the chromatin and transcriptional regulators needed for cell survival in vivo are non-overlapping with those required in vitro. We identified transcription pause-release and elongation factors as one set of in vivo-specific cancer dependencies, and determined that these factors are necessary for enhancer-mediated transcriptional adaptations that enable cells to survive the tumour microenvironment. Our lead hit, JMJD6, mediates the upregulation of in vivo stress and stimulus response pathways through enhancer-mediated transcriptional pause-release, promoting cell survival specifically in vivo. Targeting JMJD6 or other identified elongation factors extends survival in orthotopic xenograft mouse models, suggesting that targeting transcription elongation machinery may be an effective therapeutic strategy for glioblastoma. More broadly, this study demonstrates the power of in vivo phenotypic screening to identify new classes of 'cancer dependencies' not identified by previous in vitro approaches, and could supply new opportunities for therapeutic intervention.

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.

Modeling disease risk through analysis of physical interactions between genetic variants within chromatin regulatory circuitry

SNPs associated with disease susceptibility often reside in enhancer clusters, or super-enhancers. Constituents of these enhancer clusters cooperate to regulate target genes and often extend beyond the linkage disequilibrium (LD) blocks containing risk SNPs identified in genome-wide association studies (GWAS). We identified 'outside variants', defined as SNPs in weak LD with GWAS risk SNPs that physically interact with risk SNPs as part of a target gene's regulatory circuitry. These outside variants further explain variation in target gene expression beyond that explained by GWAS-associated SNPs. Additionally, the clinical risk associated with GWAS SNPs is considerably modified by the genotype of outside variants. Collectively, these findings suggest a potential model in which outside variants and GWAS SNPs that physically interact in 3D chromatin collude to influence target transcript levels as well as clinical risk. This model offers an additional hypothesis for the source of missing heritability for complex traits.

Upregulation of Glucose-Regulated Protein 78 in Metastatic Cancer Cells Is Necessary for Lung Metastasis Progression

Metastasis is the cause of more than 90% of all cancer deaths. Despite this fact, most anticancer therapeutics currently in clinical use have limited efficacy in treating established metastases. Here, we identify the endoplasmic reticulum chaperone protein, glucose-regulated protein 78 (GRP78), as a metastatic dependency in several highly metastatic cancer cell models. We find that GRP78 is consistently upregulated when highly metastatic cancer cells colonize the lung microenvironment and that mitigation of GRP78 upregulation via short hairpin RNA or treatment with the small molecule IT-139, which is currently under clinical investigation for the treatment of primary tumors, inhibits metastatic growth in the lung microenvironment. Inhibition of GRP78 upregulation and an associated reduction in metastatic potential have been shown in four highly metastatic cell line models: three human osteosarcomas and one murine mammary adenocarcinoma. Lastly, we show that downmodulation of GRP78 in highly metastatic cancer cells significantly increases median survival times in our in vivo animal model of experimental metastasis. Collectively, our data indicate that GRP78 is an attractive target for the development of antimetastatic therapies.

"PEAR-ing" Genomic and Epigenomic Analyses for Cancer Gene Discovery

Somatic structural variants in tumor genomes can deregulate transcription through repositioning of enhancer elements. A new method, PEAR-ChIP, leverages paired-end H3K27ac chromatin immunoprecipitation combined with high-throughput sequencing and current computational methods to identify such events.

Abstract LB-151: Positively selected enhancer elements endow tumor 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 signature enhancer-histone marks between near-isogenic pairs of high and low lung-metastatic osteosarcoma cells. 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 non-randomly, indicating that activity of these enhancers and their associated gene targets is positively selected. 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 whose motifs are enriched in Met-VELs, and by knockdown of individual genes activated by Met-VELs. These observations have implications for the discovery and development of targeted anti-metastatic therapies.

Citation Format: James J. Morrow, Tyler E. Miller, Alina Saiakhova, Michael M. Lizardo, Cynthia F. Bartels, Ian Bayles, Stevephen Hung, Arnulfo Mendoza, Jay T. Myers, Frederick Allen, Analisa DiFeo, Brian P. Rubin, Alex Y. Huang, Paul S. Meltzer, Lee J. Helman, Chand Khanna, Peter C. Scacheri. Positively selected enhancer elements endow tumor cells with metastatic competence. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr LB-151.

mTOR Inhibition Mitigates Enhanced mRNA Translation Associated with the Metastatic Phenotype of Osteosarcoma Cells In Vivo

PURPOSE

To successfully metastasize, tumor cells must respond appropriately to biological stressors encountered during metastatic progression. We sought to test the hypothesis that enhanced efficiency of mRNA translation during periods of metastatic stress is required for metastatic competence of osteosarcoma and that this metastasis-specific adaptation is amenable to therapeutic intervention.

EXPERIMENTAL DESIGN

We employ novel reporter and proteomic systems that enable tracking of mRNA translation efficiency and output in metastatic osteosarcoma cells as they colonize the lungs. We test the potential to target mRNA translation as an antimetastatic therapeutic strategy through pharmacokinetic studies and preclinical assessment of the prototypic mTOR inhibitor, rapamycin, across multiple models of metastasis.

RESULTS

Metastatic osteosarcoma cells translate mRNA more efficiently than nonmetastatic cells during critical stressful periods of metastatic colonization of the lung. Rapamycin inhibits translational output during periods of metastatic stress, mitigates lung colonization, and prolongs survival. mTOR-inhibiting exposures of rapamycin are achievable in mice using treatment schedules that correspond to human doses well below the MTDs defined in human patients, and as such are very likely to be tolerated over long exposures alone and in combination with other agents.

CONCLUSIONS

Metastatic competence of osteosarcoma cells is dependent on efficient mRNA translation during stressful periods of metastatic progression, and the mTOR inhibitor, rapamycin, can mitigate this translation and inhibit metastasis in vivo Our data suggest that mTOR pathway inhibitors should be reconsidered in the clinic using rationally designed dosing schedules and clinical metrics related to metastatic progression. Clin Cancer Res; 22(24); 6129-41. ©2016 AACR.

12th International CHARGE syndrome conference proceedings

The CHARGE Syndrome Foundation holds an International conference for families and professionals every other summer. In July, 2015, the 12th meeting was held in Schaumburg, Illinois, at the Renaissance Schaumburg Hotel. Day one of the 4-day conference was dedicated to professionals caring for and researching various aspects of CHARGE, including education, medical management, animal models, and stem cell-based approaches to understanding and treating individuals with CHARGE. Here, we summarize presentations from the meeting, including a synopsis of each of the three different breakout sessions (Medical/Clinical, Basic Science/CHD7, and Education), followed by a list of abstracts and authors for both platform and poster presentations.

Long ncRNA expression associates with tissue-specific enhancers

Long non-coding RNAs (ncRNA) have recently been demonstrated to be expressed from a subset of enhancers and to be required for the distant regulation of gene expression. Several approaches to predict enhancers have been developed based on various chromatin marks and occupancy of enhancer-binding proteins. Despite the rapid advances in the field, no consensus how to define tissue specific enhancers yet exists. Here, we identify 2,695 long ncRNAs annotated by ENCODE (corresponding to 28% of all ENCODE annotated long ncRNAs) that overlap tissue-specific enhancers. We use a recently developed algorithm to predict tissue-specific enhancers, PreSTIGE, that is based on the H3K4me1 mark and tissue specific expression of mRNAs. The expression of the long ncRNAs overlapping enhancers is significantly higher when the enhancer is predicted as active in a specific cell line, suggesting a general interdependency of active enhancers and expression of long ncRNAs. This dependency is not identified using previous enhancer prediction algorithms that do not account for expression of their downstream targets. The predicted enhancers that overlap annotated long ncRNAs generally have a lower ratio of H3K4me1 to H3K4me3, suggesting that enhancers expressing long ncRNAs might be associated with specific epigenetic marks. In conclusion, we demonstrate the tissue-specific predictive power of PreSTIGE and provide evidence for thousands of long ncRNAs that are expressed from active tissue-specific enhancers, suggesting a particularly important functional relationship between long ncRNAs and enhancer activity in determining tissue-specific gene expression.

Epigenomic comparison reveals activation of "seed" enhancers during transition from naive to primed pluripotency

Naive mouse embryonic stem cells (mESCs) and primed epiblast stem cells (mEpiSCs) represent successive snapshots of pluripotency during embryogenesis. Using transcriptomic and epigenomic mapping we show that a small fraction of transcripts are differentially expressed between mESCs and mEpiSCs and that these genes show expected changes in chromatin at their promoters and enhancers. Unexpectedly, the cis-regulatory circuitry of genes that are expressed at identical levels between these cell states also differs dramatically. In mESCs, these genes are associated with dominant proximal enhancers and dormant distal enhancers, which we term seed enhancers. In mEpiSCs, the naive-dominant enhancers are lost, and the seed enhancers take up primary transcriptional control. Seed enhancers have increased sequence conservation and show preferential usage in downstream somatic tissues, often expanding into super enhancers. We propose that seed enhancers ensure proper enhancer utilization and transcriptional fidelity as mammalian cells transition from naive pluripotency to a somatic regulatory program.

ZNF143 provides sequence specificity to secure chromatin interactions at gene promoters

Chromatin interactions connect distal regulatory elements to target gene promoters guiding stimulus- and lineage-specific transcription. Few factors securing chromatin interactions have so far been identified. Here, by integrating chromatin interaction maps with the large collection of transcription factor-binding profiles provided by the ENCODE project, we demonstrate that the zinc-finger protein ZNF143 preferentially occupies anchors of chromatin interactions connecting promoters with distal regulatory elements. It binds directly to promoters and associates with lineage-specific chromatin interactions and gene expression. Silencing ZNF143 or modulating its DNA-binding affinity using single-nucleotide polymorphisms (SNPs) as a surrogate of site-directed mutagenesis reveals the sequence dependency of chromatin interactions at gene promoters. We also find that chromatin interactions alone do not regulate gene expression. Together, our results identify ZNF143 as a novel chromatin-looping factor that contributes to the architectural foundation of the genome by providing sequence specificity at promoters connected with distal regulatory elements.

Chromatin interactions can connect distal regulatory elements to promoters via protein factors, but few such factors have been identified. Here, the authors show that zinc-finger protein ZNF143 is a sequence-specific chromatin-looping factor that connects promoters with distal regulatory elements.

Enhancer variants: evaluating functions in common disease

Gene enhancer elements are noncoding segments of DNA that play a central role in regulating transcriptional programs that control development, cell identity, and evolutionary processes. Recent studies have shown that noncoding single nucleotide polymorphisms (SNPs) that have been associated with risk for numerous common diseases through genome-wide association studies frequently lie in cell-type-specific enhancer elements. These enhancer variants probably influence transcriptional output, thereby offering a mechanistic basis to explain their association with risk for many common diseases. This review focuses on the identification and interpretation of disease-susceptibility variants that influence enhancer function. We discuss strategies for prioritizing the study of functional enhancer SNPs over those likely to be benign, review experimental and computational approaches to identifying the gene targets of enhancer variants, and highlight efforts to quantify the impact of enhancer variants on target transcript levels and cellular phenotypes. These studies are beginning to provide insights into the mechanistic basis of many common diseases, as well as into how we might translate this knowledge for improved disease diagnosis, prevention and treatments. Finally, we highlight five major challenges often associated with interpreting enhancer variants, and discuss recent technical advances that may help to surmount these challenges.

Trithorax monomethylates histone H3K4 and interacts directly with CBP to promote H3K27 acetylation and antagonize Polycomb silencing

Trithorax (TRX) antagonizes epigenetic silencing by Polycomb group (PcG) proteins, stimulates enhancer-dependent transcription, and establishes a 'cellular memory' of active transcription of PcG-regulated genes. The mechanisms underlying these TRX functions remain largely unknown, but are presumed to involve its histone H3K4 methyltransferase activity. We report that the SET domains of TRX and TRX-related (TRR) have robust histone H3K4 monomethyltransferase activity in vitro and that Tyr3701 of TRX and Tyr2404 of TRR prevent them from being trimethyltransferases. The trx(Z11) missense mutation (G3601S), which abolishes H3K4 methyltransferase activity in vitro, reduces the H3K4me1 but not the H3K4me3 level in vivo. trx(Z11) also suppresses the impaired silencing phenotypes of the Pc(3) mutant, suggesting that H3K4me1 is involved in antagonizing Polycomb silencing. Polycomb silencing is also antagonized by TRX-dependent H3K27 acetylation by CREB-binding protein (CBP). We show that perturbation of Polycomb silencing by TRX overexpression requires CBP. We also show that TRX and TRR are each physically associated with CBP in vivo, that TRX binds directly to the CBP KIX domain, and that the chromatin binding patterns of TRX and TRR are highly correlated with CBP and H3K4me1 genome-wide. In vitro acetylation of H3K27 by CBP is enhanced on K4me1-containing H3 substrates, and independently altering the H3K4me1 level in vivo, via the H3K4 demethylase LSD1, produces concordant changes in H3K27ac. These data indicate that the catalytic activities of TRX and CBP are physically coupled and suggest that both activities play roles in antagonizing Polycomb silencing, stimulating enhancer activity and cellular memory.

Uracil DNA glycosylase (UNG) loss enhances DNA double strand break formation in human cancer cells exposed to pemetrexed

Misincorporation of genomic uracil and formation of DNA double strand breaks (DSBs) are known consequences of exposure to TS inhibitors such as pemetrexed. Uracil DNA glycosylase (UNG) catalyzes the excision of uracil from DNA and initiates DNA base excision repair (BER). To better define the relationship between UNG activity and pemetrexed anticancer activity, we have investigated DNA damage, DSB formation, DSB repair capacity, and replication fork stability in UNG^+/+ and UNG^−/− cells. We report that despite identical growth rates and DSB repair capacities, UNG^−/− cells accumulated significantly greater uracil and DSBs compared with UNG^+/+ cells when exposed to pemetrexed. ChIP-seq analysis of γ-H2AX enrichment confirmed fewer DSBs in UNG^+/+ cells. Furthermore, DSBs in UNG^+/+ and UNG^−/− cells occur at distinct genomic loci, supporting differential mechanisms of DSB formation in UNG-competent and UNG-deficient cells. UNG^−/− cells also showed increased evidence of replication fork instability (PCNA dispersal) when exposed to pemetrexed. Thymidine co-treatment rescues S-phase arrest in both UNG^+/+ and UNG^−/− cells treated with IC₅₀-level pemetrexed. However, following pemetrexed exposure, UNG^−/− but not UNG^+/+ cells are refractory to thymidine rescue, suggesting that deficient uracil excision rather than dTTP depletion is the barrier to cell cycle progression in UNG^−/− cells. Based on these findings we propose that pemetrexed-induced uracil misincorporation is genotoxic, contributing to replication fork instability, DSB formation and ultimately cell death.

Combinatorial effects of multiple enhancer variants in linkage disequilibrium dictate levels of gene expression to confer susceptibility to common traits

DNA variants (SNPs) that predispose to common traits often localize within noncoding regulatory elements such as enhancers. Moreover, loci identified by genome-wide association studies (GWAS) often contain multiple SNPs in linkage disequilibrium (LD), any of which may be causal. Thus, determining the effect of these multiple variant SNPs on target transcript levels has been a major challenge. Here, we provide evidence that for six common autoimmune disorders (rheumatoid arthritis, Crohn's disease, celiac disease, multiple sclerosis, lupus, and ulcerative colitis), the GWAS association arises from multiple polymorphisms in LD that map to clusters of enhancer elements active in the same cell type. This finding suggests a “multiple enhancer variant” hypothesis for common traits, where several variants in LD impact multiple enhancers and cooperatively affect gene expression. Using a novel method to delineate enhancer–gene interactions, we show that multiple enhancer variants within a given locus typically target the same gene. Using available data from HapMap and B lymphoblasts as a model system, we provide evidence at numerous loci that multiple enhancer variants cooperatively contribute to altered expression of their gene targets. The effects on target transcript levels tend to be modest and can be either gain- or loss-of-function. Additionally, the genes associated with multiple enhancer variants encode proteins that are often functionally related and enriched in common pathways. Overall, the multiple enhancer variant hypothesis offers a new paradigm by which noncoding variants can confer susceptibility to common traits.

CHD7 and retinoic acid signaling cooperate to regulate neural stem cell and inner ear development in mouse models of CHARGE syndrome

CHARGE syndrome is a multiple congenital anomaly disorder that leads to life-threatening birth defects, such as choanal atresia and cardiac malformations as well as multiple sensory impairments, that affect hearing, vision, olfaction and balance. CHARGE is caused by heterozygous mutations in CHD7, which encodes an ATP-dependent chromatin remodeling enzyme. Identification of the mechanisms underlying neurological and sensory defects in CHARGE is a first step toward developing treatments for CHARGE individuals. Here, we used mouse models of Chd7 deficiency to explore the function of CHD7 in the development of the subventricular zone (SVZ) neural stem cell niche and inner ear, structures that are important for olfactory bulb neurogenesis and hearing and balance, respectively. We found that loss of Chd7 results in cell-autonomous proliferative, neurogenic and self-renewal defects in the perinatal and mature mouse SVZ stem cell niche. Modulation of retinoic acid (RA) signaling prevented in vivo inner ear and in vitro neural stem cell defects caused by Chd7 deficiency. Our findings demonstrate critical, cooperative roles for RA and CHD7 in SVZ neural stem cell function and inner ear development, suggesting that altered RA signaling may be an effective method for treating Chd7 deficiency.

Knockdown of fbxl10/kdm2bb rescues chd7 morphant phenotype in a zebrafish model of CHARGE syndrome

CHARGE syndrome is a sporadic autosomal-dominant genetic disorder characterized by a complex array of birth defects so named for its cardinal features of ocular coloboma, heart defects, choanal atresia, growth retardation, genital abnormalities, and ear abnormalities. Approximately two-thirds of individuals clinically diagnosed with CHARGE syndrome have heterozygous loss-of-function mutations in the gene encoding chromodomain helicase DNA-binding protein 7 (CHD7), an ATP-dependent chromatin remodeler. To examine the role of Chd7 in development, a zebrafish model was generated through morpholino (MO)-mediated targeting of the zebrafish chd7 transcript. High doses of chd7 MO induce lethality early in embryonic development. However, low dose-injected embryos are viable, and by 4 days post-fertilization, morphant fish display multiple defects in organ systems analogous to those affected in humans with CHARGE syndrome. The chd7 morphants show elevated expression of several potent cell-cycle inhibitors including ink4ab (p16/p15), p21 and p27, accompanied by reduced cell proliferation. We also show that Chd7 is required for proper organization of neural crest-derived craniofacial cartilage structures. Strikingly, MO-mediated knockdown of the jumonji domain-containing histone demethylase fbxl10/kdm2bb, a repressor of ribosomal RNA (rRNA) genes, rescues cell proliferation and cartilage defects in chd7 morphant embryos and can lead to complete rescue of the CHARGE syndrome phenotype. These results indicate that CHARGE-like phenotypes in zebrafish can be mitigated through modulation of fbxl10 levels and implicate FBXL10 as a possible therapeutic target in CHARGE syndrome.

Twist1 mediates repression of chondrogenesis by β-catenin to promote cranial bone progenitor specification

The bones of the mammalian skull vault form through intramembranous ossification. Skull bones ossify directly, in a process regulated by β-catenin, instead of passing through a cartilage intermediate. We tested whether β-catenin is necessary for fate selection of intramembranous bone progenitors in the skull. Here, we show in mice that removal of β-catenin from skull bone progenitors results in the near complete transformation of the skull bones to cartilage, whereas constitutive β-catenin activation inhibits skull bone fate selection. β-catenin directly activated Twist1 expression in skull progenitors, conditional Twist1 deletion partially phenocopied the absence of β-catenin, and Twist1 deletion partially restored bone formation in the presence of constitutive β-catenin activation. Finally, Twist1 bound robustly to the 3'UTR of Sox9, the central initiator of chondrogenesis, suggesting that Twist1 might directly repress cartilage formation through Sox9. These findings provide insight into how β-catenin signaling via Twist1 actively suppresses the formation of cartilage and promotes intramembranous ossification in the skull.

The chromatin fingerprint of gene enhancer elements

Different cell types within a single organism are generally distinguished by strikingly different patterns of gene expression, which are dynamic throughout development and adult life. Distal enhancer elements are key drivers of spatiotemporal specificity in gene regulation. Often located tens of kilobases from their target promoters and functioning in an orientation-independent manner, the identification of bona fide enhancers has proved a formidable challenge. With the development of ChIP-seq, global cataloging of putative enhancers has become feasible. Here, we review the current understanding of the chromatin landscape at enhancers and how these chromatin features enable robust identification of tissue-specific enhancers.

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.

Otitis media in a new mouse model for CHARGE syndrome with a deletion in the Chd7 gene

Otitis media is a middle ear disease common in children under three years old. Otitis media can occur in normal individuals with no other symptoms or syndromes, but it is often seen in individuals clinically diagnosed with genetic diseases such as CHARGE syndrome, a complex genetic disease caused by mutation in the Chd7 gene and characterized by multiple birth defects. Although otitis media is common in human CHARGE syndrome patients, it has not been reported in mouse models of CHARGE syndrome. In this study, we report a mouse model with a spontaneous deletion mutation in the Chd7 gene and with chronic otitis media of early onset age accompanied by hearing loss. These mice also exhibit morphological alteration in the Eustachian tubes, dysregulation of epithelial proliferation, and decreased density of middle ear cilia. Gene expression profiling revealed up-regulation of Muc5ac, Muc5b and Tgf-β1 transcripts, the products of which are involved in mucin production and TGF pathway regulation. This is the first mouse model of CHARGE syndrome reported to show otitis media with effusion and it will be valuable for studying the etiology of otitis media and other symptoms in CHARGE syndrome.