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.
Display omitted
•HiCorr allows robust mapping of sub-TAD chromatin interactions with Hi-C•Low-input “easy Hi-C” protocol compatible with 50–100k cells•Enhancer loops and aggregates are better marks of cell identity than compartments•Chromatin loops outperform eQTLs in defining neurological GWAS target genes
Lu et al. developed a rigorous Hi-C bias-correction pipeline to significantly improve the robustness of high-resolution chromatin interaction maps. With a new low-input “easy Hi-C” protocol, they mapped chromatin interactions in neural samples, defined cell-type-specific enhancer loops and aggregates, and concluded that Hi-C outperforms eQTL in explaining GWAS results.
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.
Display omitted
•Enhancers active in the cell of origin are co-amplified with oncogenes•Circular extrachromosomal amplicons are associated with enhancer rewiring•Endogenous and new enhancers on amplicons contribute to cell proliferation•Skewed co-amplification that selects enhancers is found across several tumor types
Extrachromosomal oncogene amplification is a common occurrence across a broad range of cancers, yet the chromatin landscape of these high-level amplifications is poorly understood. Using a combination of approaches to explore their chromatin topology and enhancer landscape, Morton et al. observed that these oncogene amplifications are shaped by regulatory dependencies in the non-coding genome.
CHD7 is one of nine members of the chromodomain helicase DNA-binding domain family of ATP-dependent chromatin remodeling enzymes found in mammalian cells. De novo mutation of CHD7 is a major cause of ...CHARGE syndrome, a genetic condition characterized by multiple congenital anomalies. To gain insights to the function of CHD7, we used the technique of chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP-Seq) to map CHD7 sites in mouse ES cells. We identified 10,483 sites on chromatin bound by CHD7 at high confidence. Most of the CHD7 sites show features of gene enhancer elements. Specifically, CHD7 sites are predominantly located distal to transcription start sites, contain high levels of H3K4 mono-methylation, found within open chromatin that is hypersensitive to DNase I digestion, and correlate with ES cell-specific gene expression. Moreover, CHD7 co-localizes with P300, a known enhancer-binding protein and strong predictor of enhancer activity. Correlations with 18 other factors mapped by ChIP-seq in mouse ES cells indicate that CHD7 also co-localizes with ES cell master regulators OCT4, SOX2, and NANOG. Correlations between CHD7 sites and global gene expression profiles obtained from Chd7(+/+), Chd7(+/-), and Chd7(-/-) ES cells indicate that CHD7 functions at enhancers as a transcriptional rheostat to modulate, or fine-tune the expression levels of ES-specific genes. CHD7 can modulate genes in either the positive or negative direction, although negative regulation appears to be the more direct effect of CHD7 binding. These data indicate that enhancer-binding proteins can limit gene expression and are not necessarily co-activators. Although ES cells are not likely to be affected in CHARGE syndrome, we propose that enhancer-mediated gene dysregulation contributes to disease pathogenesis and that the critical CHD7 target genes may be subject to positive or negative regulation.
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.
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.
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.
CHD7 is a member of the chromodomain helicase DNA binding domain family of ATP-dependent chromatin remodeling enzymes. De novo mutation of the CHD7 gene is a major cause of CHARGE syndrome, a genetic ...disease characterized by a complex constellation of birth defects (Coloboma of the eye, Heart defects, Atresia of the choanae, severe Retardation of growth and development, Genital abnormalities, and Ear abnormalities). To gain insight into the function of CHD7, we mapped the distribution of the CHD7 protein on chromatin using the approach of chromatin immunoprecipitation on tiled microarrays (ChIP-chip). These studies were performed in human colorectal carcinoma cells, human neuroblastoma cells, and mouse embryonic stem (ES) cells before and after differentiation into neural precursor cells. The results indicate that CHD7 localizes to discrete locations along chromatin that are specific to each cell type, and that the cell-specific binding of CHD7 correlates with a subset of histone H3 methylated at lysine 4 (H3K4me). The CHD7 sites change concomitantly with H3K4me patterns during ES cell differentiation, suggesting that H3K4me is part of the epigenetic signature that defines lineage-specific association of CHD7 with specific sites on chromatin. Furthermore, the CHD7 sites are predominantly located distal to transcription start sites, most often contained within DNase hypersensitive sites, frequently conserved, and near genes expressed at relatively high levels. These features are similar to those of gene enhancer elements, raising the possibility that CHD7 functions in enhancer mediated transcription, and that the congenital anomalies in CHARGE syndrome are due to alterations in transcription of tissue-specific genes normally regulated by CHD7 during development.
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.
Mutations in the low-density lipoprotein receptor–related protein 5 gene (LRP5) cause autosomal recessive osteoporosis-pseudoglioma syndrome (OPPG). We sequenced the coding exons of LRP5 in 37 ...probands suspected of having OPPG on the basis of the co-occurrence of severe congenital or childhood-onset visual impairment with bone fragility or osteoporosis recognized by young adulthood. We found two putative mutant alleles in 26 probands, only one mutant allele in 4 probands, and no mutant alleles in 7 probands. Looking for digenic inheritance, we sequenced the genes encoding the functionally related receptor LRP6, an LRP5 coreceptor FZD4, and an LRP5 ligand, NDP, in the four probands with one mutant allele, and, looking for locus heterogeneity, we sequenced FZD4 and NDP in the seven probands with no mutations, but we found no additional mutations. When we compared clinical features between probands with and without LRP5 mutations, we found no difference in the severity of skeletal disease, prevalence of cognitive impairment, or family history of consanguinity. However, four of the seven probands without detectable mutations had eye pathology that differed from pathology previously described for OPPG. Since many LRP5 mutations are missense changes, to differentiate between a disease-causing mutation and a benign variant, we measured the ability of wild-type and mutant LRP5 to transduce Wnt and Norrin signal ex vivo. Each of the seven OPPG mutations tested, had reduced signal transduction compared with wild-type mutations. These results indicate that early bilateral vitreoretinal eye pathology coupled with skeletal fragility is a strong predictor of LRP5 mutation and that mutations in LRP5 cause OPPG by impairing Wnt and Norrin signal transduction.
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.
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).
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).
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.