Mice lacking a Myc enhancer that includes human SNP rs6983267 are resistant to intestinal tumors

The interplay of long non-coding RNAs and MYC in cancer

Long non-coding RNAs (lncRNAs) are a class of RNA molecules that are changing how researchers view eukaryotic gene regulation. Once considered to be non-functional products of low-level aberrant transcription from non-coding regions of the genome, lncRNAs are now viewed as important epigenetic regulators and several lncRNAs have now been demonstrated to be critical players in the development and/or maintenance of cancer. Similarly, the emerging variety of interactions between lncRNAs and MYC, a well-known oncogenic transcription factor linked to most types of cancer, have caught the attention of many biomedical researchers. Investigations exploring the dynamic interactions between lncRNAs and MYC, referred to as the lncRNA-MYC network, have proven to be especially complex. Genome-wide studies have shown that MYC transcriptionally regulates many lncRNA genes. Conversely, recent reports identified lncRNAs that regulate MYC expression both at the transcriptional and post-transcriptional levels. These findings are of particular interest because they suggest roles of lncRNAs as regulators of MYC oncogenic functions and the possibility that targeting lncRNAs could represent a novel avenue to cancer treatment. Here, we briefly review the current understanding of how lncRNAs regulate chromatin structure and gene transcription, and then focus on the new developments in the emerging field exploring the lncRNA-MYC network in cancer.

Regulation of MYC gene expression by aberrant Wnt/β-catenin signaling in colorectal cancer

The Wnt/β-catenin signaling pathway controls intestinal homeostasis and mutations in components of this pathway are prevalent in human colorectal cancers (CRCs). These mutations lead to inappropriate expression of genes controlled by Wnt responsive DNA elements (WREs). T-cell factor/Lymphoid enhancer factor transcription factors bind WREs and recruit the β-catenin transcriptional co-activator to activate target gene expression. Deregulated expression of the c-MYC proto-oncogene (MYC) by aberrant Wnt/β-catenin signaling drives colorectal carcinogenesis. In this review, we discuss the current literature pertaining to the identification and characterization of WREs that control oncogenic MYC expression in CRCs. A common theme has emerged whereby these WREs often map distally to the MYC genomic locus and control MYC gene expression through long-range chromatin loops with the MYC proximal promoter. We propose that by determining which of these WREs is critical for CRC pathogenesis, novel strategies can be developed to treat individuals suffering from this disease.

The complexity of epigenetic diseases

Over the past 30 years, a plethora of pathogenic mutations affecting enhancer regions and epigenetic regulators have been identified. Coupled with more recent genome‐wide association studies (GWAS) and epigenome‐wide association studies (EWAS) implicating major roles for regulatory mutations in disease, it is clear that epigenetic mechanisms represent important biomarkers for disease development and perhaps even therapeutic targets. Here, we discuss the diversity of disease‐causing mutations in enhancers and epigenetic regulators, with a particular focus on cancer. © 2015 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Demystifying the secret mission of enhancers: linking distal regulatory elements to target genes

Enhancers are short regulatory sequences bound by sequence-specific transcription factors and play a major role in the spatiotemporal specificity of gene expression patterns in development and disease. While it is now possible to identify enhancer regions genomewide in both cultured cells and primary tissues using epigenomic approaches, it has been more challenging to develop methods to understand the function of individual enhancers because enhancers are located far from the gene(s) that they regulate. However, it is essential to identify target genes of enhancers not only so that we can understand the role of enhancers in disease but also because this information will assist in the development of future therapeutic options. After reviewing models of enhancer function, we discuss recent methods for identifying target genes of enhancers. First, we describe chromatin structure-based approaches for directly mapping interactions between enhancers and promoters. Second, we describe the use of correlation-based approaches to link enhancer state with the activity of nearby promoters and/or gene expression. Third, we describe how to test the function of specific enhancers experimentally by perturbing enhancer–target relationships using high-throughput reporter assays and genome editing. Finally, we conclude by discussing as yet unanswered questions concerning how enhancers function, how target genes can be identified, and how to distinguish direct from indirect changes in gene expression mediated by individual enhancers.

Long non-coding RNA CARLo-5 expression is associated with disease progression and predicts outcome in hepatocellular carcinoma patients

Recently, many studies show that long non-coding RNAs (lncRNAs) play important roles in cancer biology. Although its expression was reported dysregulated during tumorigenesis, the contributions of lncRNAs to hepatocellular carcinoma (HCC) are still largely unknown. In particular, the lncRNA CARLo-5 has a functional role in cell-cycle regulation in colon cancer, while the clinical significance and biological function of CARLo-5 in HCC remain unelucidated. In order to fill those study blanks, the expression level of CARLo-5 in human HCC specimens was tested, and its correlation with clinicopathologic features as well as the prognosis for patients with HCC was analyzed. Additionally, MTT, wound healing and transwell assays were employed to investigate the biological function of CARLo-5. The results showed that CARLo-5 levels were significantly overexpressed in HCC tissues compared to ANLT. Besides, high expression of CARLo-5 was associated with liver cirrhosis (P = 0.001), tumor number (P < 0.001), vascular invasion (P = 0.001), capsular formation (P = 0.014) and Edmondson-Steiner grade (P < 0.001), which proved that CARLo-5 was an independent risk factor for overall survival and disease-free survival. In addition, in highly metastatic HCC cell lines (HCCLM3 and MHCC97-L), CARLo-5 was up-regulated, but in lowly metastatic HCC cell lines (HepG2, SNU387), it showed down-regulated. Besides, by using gain and loss of function experiments in HCC cell lines (HCCLM3 and HepG2), the results showed that CARLo-5 overexpression significantly enhanced cell proliferation, migration and invasion in vitro. Our study also revealed that CARLo-5 was prominently up-regulated in HCC specimens and its high expression was associated with poor prognosis of HCC patients. Totally, those findings together indicate that CARLo-5 promotes proliferation and metastasis of HCC and potentially emerged as a novel therapeutic target.

Genetic architecture of colorectal cancer

Colorectal cancer (CRC) is a complex disease that develops as a consequence of both genetic and environmental risk factors. A small proportion (3-5%) of cases arise from hereditary syndromes predisposing to early onset CRC as a result of mutations in over a dozen well defined genes. In contrast, CRC is predominantly a late onset 'sporadic' disease, developing in individuals with no obvious hereditary syndrome. In recent years, genome wide association studies have discovered that over 40 genetic regions are associated with weak effects on sporadic CRC, and it has been estimated that increasingly large genome wide scans will identify many additional novel genetic regions. Subsequent experimental validations have identified the causally related variant(s) in a limited number of these genetic regions. Further biological insight could be obtained through ethnically diverse study populations, larger genetic sequencing studies and development of higher throughput functional experiments. Along with inherited variation, integration of the tumour genome may shed light on the carcinogenic processes in CRC. In addition to summarising the genetic architecture of CRC, this review discusses genetic factors that modify environmental predictors of CRC, as well as examples of how genetic insight has improved clinical surveillance, prevention and treatment strategies. In summary, substantial progress has been made in uncovering the genetic architecture of CRC, and continued research efforts are expected to identify additional genetic risk factors that further our biological understanding of this disease. Subsequently these new insights will lead to improved treatment and prevention of colorectal cancer.

Network-Based Integration of GWAS and Gene Expression Identifies a HOX-Centric Network Associated with Serous Ovarian Cancer Risk

BACKGROUND

Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations.

METHODS

We selected TF genes within 1 Mb of the top signal at the 12 genome-wide significant risk loci. Mutual information, a form of correlation, was used to build networks of genes strongly coexpressed with each selected TF gene in the unified microarray dataset of 489 serous EOC tumors from The Cancer Genome Atlas. Genes represented in this dataset were subsequently ranked using a gene-level test based on results for germline SNPs from a serous EOC GWAS meta-analysis (2,196 cases/4,396 controls).

RESULTS

Gene set enrichment analysis identified six networks centered on TF genes (HOXB2, HOXB5, HOXB6, HOXB7 at 17q21.32 and HOXD1, HOXD3 at 2q31) that were significantly enriched for genes from the risk-associated end of the ranked list (P < 0.05 and FDR < 0.05). These results were replicated (P < 0.05) using an independent association study (7,035 cases/21,693 controls). Genes underlying enrichment in the six networks were pooled into a combined network.

CONCLUSION

We identified a HOX-centric network associated with serous EOC risk containing several genes with known or emerging roles in serous EOC development.

IMPACT

Network analysis integrating large, context-specific datasets has the potential to offer mechanistic insights into cancer susceptibility and prioritize genes for experimental characterization.

The long noncoding RNA regulation at the MYC locus

Aberrant expression of long noncoding RNAs (lncRNAs) has been linked to cancers. The MYC oncoprotein is a key contributor to the development of many human tumors. Recent studies have revealed that a number of lncRNAs originating from the human 8q24 locus previously known to corresponding to a 'gene desert' are transcribed and play important roles in MYC regulation. In this review, we highlight recent progress in how these lncRNAs participate in control of MYC levels in normal and tumor cells.

Evolutionary determinants of cancer

Our understanding of cancer is being transformed by exploring clonal diversity, drug resistance, and causation within an evolutionary framework. The therapeutic resilience of advanced cancer is a consequence of its character as a complex, dynamic, and adaptive ecosystem engendering robustness, underpinned by genetic diversity and epigenetic plasticity. The risk of mutation-driven escape by self-renewing cells is intrinsic to multicellularity but is countered by multiple restraints, facilitating increasing complexity and longevity of species. But our own species has disrupted this historical narrative by rapidly escalating intrinsic risk. Evolutionary principles illuminate these challenges and provide new avenues to explore for more effective control.

SIGNIFICANCE

Lifetime risk of cancer now approximates to 50% in Western societies. And, despite many advances, the outcome for patients with disseminated disease remains poor, with drug resistance the norm. An evolutionary perspective may provide a clearer understanding of how cancer clones develop robustness and why, for us as a species, risk is now off the scale. And, perhaps, of what we might best do to achieve more effective control.

Mouse models of colorectal cancer as preclinical models

In this review, we discuss the application of mouse models to the identification and pre-clinical validation of novel therapeutic targets in colorectal cancer, and to the search for early disease biomarkers. Large-scale genomic, transcriptomic and epigenomic profiling of colorectal carcinomas has led to the identification of many candidate genes whose direct contribution to tumourigenesis is yet to be defined; we discuss the utility of cross-species comparative 'omics-based approaches to this problem. We highlight recent progress in modelling late-stage disease using mice, and discuss ways in which mouse models could better recapitulate the complexity of human cancers to tackle the problem of therapeutic resistance and recurrence after surgical resection.

Can genome engineering be used to target cancer-associated enhancers?

Transcriptional misregulation is involved in the development of many diseases, especially neoplastic transformation. Distal regulatory elements, such as enhancers, play a major role in specifying cell-specific transcription patterns in both normal and diseased tissues, suggesting that enhancers may be prime targets for therapeutic intervention. By focusing on modulating gene regulation mediated by cell type-specific enhancers, there is hope that normal epigenetic patterning in an affected tissue could be restored with fewer side effects than observed with treatments employing relatively nonspecific inhibitors such as epigenetic drugs. New methods employing genomic nucleases and site-specific epigenetic regulators targeted to specific genomic regions, using either artificial DNA-binding proteins or RNA-DNA interactions, may allow precise genome engineering at enhancers. However, this field is still in its infancy and further refinements that increase specificity and efficiency are clearly required.

Evolutionary determinants of cancer

Our understanding of cancer is being transformed by exploring clonal diversity, drug resistance, and causation within an evolutionary framework. The therapeutic resilience of advanced cancer is a consequence of its character as a complex, dynamic, and adaptive ecosystem engendering robustness, underpinned by genetic diversity and epigenetic plasticity. The risk of mutation-driven escape by self-renewing cells is intrinsic to multicellularity but is countered by multiple restraints, facilitating increasing complexity and longevity of species. But our own species has disrupted this historical narrative by rapidly escalating intrinsic risk. Evolutionary principles illuminate these challenges and provide new avenues to explore for more effective control.

SIGNIFICANCE

Lifetime risk of cancer now approximates to 50% in Western societies. And, despite many advances, the outcome for patients with disseminated disease remains poor, with drug resistance the norm. An evolutionary perspective may provide a clearer understanding of how cancer clones develop robustness and why, for us as a species, risk is now off the scale. And, perhaps, of what we might best do to achieve more effective control.

TRIP through the chromatin: a high throughput exploration of enhancer regulatory landscapes

Enhancers are regulatory elements that promote gene expression in a spatio-temporal way and are involved in a wide range of developmental and disease processes. Both the identification and subsequent functional dissection of enhancers are key steps in understanding these processes. Several high-throughput approaches were recently developed for these purposes; however, in almost all cases enhancers are being tested outside their native chromatin context. Until recently, the analysis of enhancer activities at their native genomic locations was low throughput, laborious and time-consuming. Here, we discuss the potential of a powerful approach, TRIP, to study the functioning of enhancers in their native chromatin environments by introducing sensor constructs directly in the genome. TRIP allows for simultaneously analyzing the quantitative readout of numerous sensor constructs integrated at random locations in the genome. The high-throughput and flexible nature of TRIP opens up potential to study different aspects of enhancer biology at an unprecedented level.

Learning about mammalian gene regulation from functional enhancer assays in the mouse

Enhancer biology is emerging as a critical area of research that informs studies of evolution, development, and disease. From early experiments that defined and mapped the first enhancers to recent enhancer models of human disease, functional experiments in the mouse have played a central role in revealing enhancer biology. Three decades of in vivo enhancer studies in mouse have laid the groundwork for the current understanding of mammalian enhancers, demonstrating the developmental and tissue-specific activity of enhancers and illuminating general features of enhancer evolution and function. Recent studies offer an emerging perspective on the importance of chromosomal context, combinatorial enhancer activity, and the functional consequences of enhancer sequence variation. This review describes the basic principles of functional testing in mouse, summarizes the contributions these studies have made to our understanding of enhancer biology, and describes limitations and future outlook of in vivo mouse enhancer studies.

Mapping long-range promoter contacts in human cells with high-resolution capture Hi-C

Transcriptional control in large genomes often requires looping interactions between distal DNA elements, such as enhancers and target promoters. Current chromosome conformation capture techniques do not offer sufficiently high resolution to interrogate these regulatory interactions on a genomic scale. Here we use Capture Hi-C (CHi-C), an adapted genome conformation assay, to examine the long-range interactions of almost 22,000 promoters in 2 human blood cell types. We identify over 1.6 million shared and cell type-restricted interactions spanning hundreds of kilobases between promoters and distal loci. Transcriptionally active genes contact enhancer-like elements, whereas transcriptionally inactive genes interact with previously uncharacterized elements marked by repressive features that may act as long-range silencers. Finally, we show that interacting loci are enriched for disease-associated SNPs, suggesting how distal mutations may disrupt the regulation of relevant genes. This study provides new insights and accessible tools to dissect the regulatory interactions that underlie normal and aberrant gene regulation.

Inferring regulatory element landscapes and transcription factor networks from cancer methylomes

Recent studies indicate that DNA methylation can be used to identify transcriptional enhancers, but no systematic approach has been developed for genome-wide identification and analysis of enhancers based on DNA methylation. We describe ELMER (Enhancer Linking by Methylation/Expression Relationships), an R-based tool that uses DNA methylation to identify enhancers and correlates enhancer state with expression of nearby genes to identify transcriptional targets. Transcription factor motif analysis of enhancers is coupled with expression analysis of transcription factors to infer upstream regulators. Using ELMER, we investigated more than 2,000 tumor samples from The Cancer Genome Atlas. We identified networks regulated by known cancer drivers such as GATA3 and FOXA1 (breast cancer), SOX17 and FOXA2 (endometrial cancer), and NFE2L2, SOX2, and TP63 (squamous cell lung cancer). We also identified novel networks with prognostic associations, including RUNX1 in kidney cancer. We propose ELMER as a powerful new paradigm for understanding the cis-regulatory interface between cancer-associated transcription factors and their functional target genes.

Heterospecific SNP diversity in humans and rhesus macaque (Macaca mulatta)

BACKGROUND

Conservation of single nucleotide polymorphisms (SNPs) between human and other primates (i.e., heterospecific SNPs) in candidate genes can be used to assess the utility of those organisms as models for human biomedical research.

METHODS

A total of 59,691 heterospecific SNPs in 22 rhesus macaques and 20 humans were analyzed for human trait associations and 4207 heterospecific SNPs biallelic in both taxa were compared for genetic variation.

RESULTS

Variation comparisons at the 4207 SNPs showed that humans were more genetically diverse than rhesus macaques with observed and expected heterozygosities of 0.337 and 0.323 vs. 0.119 and 0.102, and minor allele frequencies of 0.239 and 0.063, respectively. In total, 431 of the 59,691 heterospecific SNPs are reportedly associated with human-specific traits.

CONCLUSION

While comparisons between human and rhesus macaque genomes are plausible, functional studies of heterospecific SNPs are necessary to determine whether rhesus macaque alleles are associated with the same phenotypes as their corresponding human alleles.

Gene analysis techniques and susceptibility gene discovery in non-BRCA1/BRCA2 familial breast cancer

Breast cancer is the leading cause of cancer deaths in females worldwide occurring in both hereditary and sporadic forms. Women with inherited pathogenic mutations in the BRCA1 or BRCA2 genes have up to an 85% risk of developing breast cancer in their lifetimes. These patients are candidates for risk-reduction measures such as intensive radiological screening, prophylactic surgery or chemoprevention. However, only about 20% of familial breast cancer cases are attributed to mutations in BRCA1 and BRCA2, while a further 5-10% are attributed to mutations in other rare susceptibility genes such as TP53, STK11, PTEN, ATM and CHEK2. A multitude of genome wide association studies (GWAS) have been conducted confirming low-risk common variants associated with breast cancer in excess of 90 loci, which may contribute to a further 23% of the heritability. We currently find ourselves in "the next generation", with technologies offering deep sequencing at a fraction of the cost. Starting off primarily in a research setting, multi-gene panel testing is now utilized in the clinic to sequence multiple predisposing genes simultaneously (otherwise known as multi-gene panel testing). In this review, we focus on the hereditary breast cancer discoveries, techniques and the challenges we face in this complex disease, especially in the light of the vast amount of data we now have at hand. It has been 20 years since the first breast cancer susceptibility gene has been discovered and there has been substantial progress in unraveling the genetic component of the disease. However, hereditary breast cancer remains a challenging topic subject to common debate.

The type 2 diabetes associated rs7903146 T allele within TCF7L2 is significantly under-represented in Hereditary Multiple Exostoses: insights into pathogenesis

Hereditary Multiple Exostoses (HME) is an autosomal-dominant disorder characterized by benign cartilage tumors (exostoses) forming near the growth plates, leading to severe health problems. EXT1 and EXT2 are the two genes known to harbor heterozygous loss-of-function mutations that account for the vast majority of the primary genetic component of HME. However, patients present with wide clinical heterogeneity, suggesting that modifier genes play a role in determining severity. Our previous work has pointed to an imbalance of β-catenin signaling being involved in the pathogenesis of osteochondroma formation. TCF7L2 is one of the key 'gate-keeper' TCF family members for Wnt/β-catenin signaling pathway, and TCF7L2 and EXT2 are among the earliest associated loci reported in genome wide appraisals of type 2 diabetes (T2D). Thus we investigated if the key T allele of single nucleotide polymorphism (SNP) rs7903146 within the TCF7L2 locus, which is strongly over-represented among T2D cases, was also associated with HME. We leveraged genotype data available from ongoing GWAS efforts from genomics and orthopedic centers in the US, Canada and Italy. Collectively 213 cases and 1890 controls were analyzed and, surprisingly, the T allele was in fact significantly under-represented in the HME patient group [P = 0.009; odds ratio = 0.737 (95% C.I. 0.587-0.926)]; in addition, the direction of effect was consistent within each individual cohort. Immunohistochemical analyses revealed that TCF7L2 is differentially expressed and distributed in normal human growth plate zones, and exhibits substantial variability in human exostoses in terms of staining intensity and distribution. In summary, the data indicate that there is a putative genetic connection between TCF7L2 and EXT in the context of HME. Given this observation, we suggest that these loci could possibly modulate shared pathways, in particular with respect to β-catenin, and their respective variants interplay to influence HME pathogenesis as well as T2D.

The selection and function of cell type-specific enhancers

The human body contains several hundred cell types, all of which share the same genome. In metazoans, much of the regulatory code that drives cell type-specific gene expression is located in distal elements called enhancers. Although mammalian genomes contain millions of potential enhancers, only a small subset of them is active in a given cell type. Cell type-specific enhancer selection involves the binding of lineage-determining transcription factors that prime enhancers. Signal-dependent transcription factors bind to primed enhancers, which enables these broadly expressed factors to regulate gene expression in a cell type-specific manner. The expression of genes that specify cell type identity and function is associated with densely spaced clusters of active enhancers known as super-enhancers. The functions of enhancers and super-enhancers are influenced by, and affect, higher-order genomic organization.

CRISPR-Cas9 genome editing of a single regulatory element nearly abolishes target gene expression in mice--brief report

OBJECTIVE

To ascertain the importance of a single regulatory element in the control of Cnn1 expression using CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9) genome editing.

APPROACH AND RESULTS

The CRISPR/Cas9 system was used to produce 3 of 18 founder mice carrying point mutations in an intronic CArG box of the smooth muscle cell-restricted Cnn1 gene. Each founder was bred for germline transmission of the mutant CArG box and littermate interbreeding to generate homozygous mutant (Cnn1(ΔCArG/ΔCArG)) mice. Quantitative reverse transcription polymerase chain reaction, Western blotting, and confocal immunofluorescence microscopy showed dramatic reductions in Cnn1 mRNA and CNN1 protein expression in Cnn1(ΔCArG/ΔCArG) mice with no change in other smooth muscle cell-restricted genes and little evidence of off-target edits elsewhere in the genome. In vivo chromatin immunoprecipitation assay revealed a sharp decrease in binding of serum response factor to the mutant CArG box. Loss of CNN1 expression was coincident with an increase in Ki-67 positive cells in the normal vessel wall.

CONCLUSIONS

CRISPR/Cas9 genome editing of a single CArG box nearly abolishes Cnn1 expression in vivo and evokes increases in smooth muscle cell DNA synthesis. This facile genome editing system paves the way for a new generation of studies designed to test the importance of individual regulatory elements in living animals, including regulatory variants in conserved sequence blocks linked to human disease.

Junk DNA and the long non-coding RNA twist in cancer genetics

The central dogma of molecular biology states that the flow of genetic information moves from DNA to RNA to protein. However, in the last decade this dogma has been challenged by new findings on non-coding RNAs (ncRNAs) such as microRNAs (miRNAs). More recently, long non-coding RNAs (lncRNAs) have attracted much attention due to their large number and biological significance. Many lncRNAs have been identified as mapping to regulatory elements including gene promoters and enhancers, ultraconserved regions, and intergenic regions of protein-coding genes. Yet, the biological function and molecular mechanisms of lncRNA in human diseases in general and cancer in particular remain largely unknown. Data from the literature suggest that lncRNA, often via interaction with proteins, functions in specific genomic loci or use their own transcription loci for regulatory activity. In this review, we summarize recent findings supporting the importance of DNA loci in lncRNA function, and the underlying molecular mechanisms via cis or trans regulation, and discuss their implications in cancer. In addition, we use the 8q24 genomic locus, a region containing interactive SNPs, DNA regulatory elements and lncRNAs, as an example to illustrate how single nucleotide polymorphism (SNP) located within lncRNAs may be functionally associated with the individual’s susceptibility to cancer.

A NOTCH1-driven MYC enhancer promotes T cell development, transformation and acute lymphoblastic leukemia

Efforts to identify and annotate cancer driver genetic lesions have been focused primarily on the analysis of protein-coding genes; however, most genetic abnormalities found in human cancer are located in intergenic regions. Here we identify a new long range-acting MYC enhancer controlled by NOTCH1 that is targeted by recurrent chromosomal duplications in human T cell acute lymphoblastic leukemia (T-ALL). This highly conserved regulatory element, hereby named N-Me for NOTCH MYC enhancer, is located within a broad super-enhancer region +1.47 Mb from the MYC transcription initiating site, interacts with the MYC proximal promoter and induces orientation-independent MYC expression in reporter assays. Moreover, analysis of N-Me knockout mice demonstrates a selective and essential role of this regulatory element during thymocyte development and in NOTCH1-induced T-ALL. Together these results identify N-Me as a long-range oncogenic enhancer implicated directly in the pathogenesis of human leukemia and highlight the importance of the NOTCH1-MYC regulatory axis in T cell transformation and as a therapeutic target in T-ALL.

Expression QTL-based analyses reveal candidate causal genes and loci across five tumor types

The majority of trait-associated loci discovered through genome-wide association studies are located outside of known protein coding regions. Consequently, it is difficult to ascertain the mechanism underlying these variants and to pinpoint the causal alleles. Expression quantitative trait loci (eQTLs) provide an organizing principle to address both of these issues. eQTLs are genetic loci that correlate with RNA transcript levels. Large-scale data sets such as the Cancer Genome Atlas (TCGA) provide an ideal opportunity to systematically evaluate eQTLs as they have generated multiple data types on hundreds of samples. We evaluated the determinants of gene expression (germline variants and somatic copy number and methylation) and performed cis-eQTL analyses for mRNA expression and miRNA expression in five tumor types (breast, colon, kidney, lung and prostate). We next tested 149 known cancer risk loci for eQTL effects, and observed that 42 (28.2%) were significantly associated with at least one transcript. Lastly, we described a fine-mapping strategy for these 42 eQTL target-gene associations based on an integrated strategy that combines the eQTL level of significance and the regulatory potential as measured by DNaseI hypersensitivity. For each of the risk loci, our analyses suggested 1 to 81 candidate causal variants that may be prioritized for downstream functional analysis. In summary, our study provided a comprehensive landscape of the genetic determinants of gene expression in different tumor types and ranked the genes and loci for further functional assessment of known cancer risk loci.

Functional annotation of colon cancer risk SNPs

Colorectal cancer (CRC) is a leading cause of cancer-related deaths in the United States. Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with increased risk for CRC. A molecular understanding of the functional consequences of this genetic variation has been complicated because each GWAS SNP is a surrogate for hundreds of other SNPs, most of which are located in non-coding regions. Here we use genomic and epigenomic information to test the hypothesis that the GWAS SNPs and/or correlated SNPs are in elements that regulate gene expression, and identify 23 promoters and 28 enhancers. Using gene expression data from normal and tumour cells, we identify 66 putative target genes of the risk-associated enhancers (10 of which were also identified by promoter SNPs). Employing CRISPR nucleases, we delete one risk-associated enhancer and identify genes showing altered expression. We suggest that similar studies be performed to characterize all CRC risk-associated enhancers.

Previous studies identified genetic variants associated with colorectal cancer (CRC), but the functional consequences of these genetic risk factors remain poorly understood. Here, the authors report that CRC risk variants reside in promoters and enhancers and could increase colon cancer risk through gene expression regulation.

Genome-wide interrogation of longitudinal FEV1 in children with asthma

RATIONALE

Most genomic studies of lung function have used phenotypic data derived from a single time-point (e.g., presence/absence of disease) without considering the dynamic progression of a chronic disease.

OBJECTIVES

To characterize lung function change over time in subjects with asthma and identify genetic contributors to a longitudinal phenotype.

METHODS

We present a method that models longitudinal FEV1 data, collected from 1,041 children with asthma who participated in the Childhood Asthma Management Program. This longitudinal progression model was built using population-based nonlinear mixed-effects modeling with an exponential structure and the determinants of age and height.

MEASUREMENTS AND MAIN RESULTS

We found ethnicity was a key covariate for FEV1 level. Budesonide-treated children with asthma had a slight but significant effect on FEV1 when compared with those treated with placebo or nedocromil (P < 0.001). A genome-wide association study identified seven single-nucleotide polymorphisms nominally associated with longitudinal lung function phenotypes in 581 white Childhood Asthma Management Program subjects (P < 10(-4) in the placebo ["discovery"] and P < 0.05 in the nedocromil treatment ["replication"] group). Using ChIP-seq and RNA-seq data, we found that some of the associated variants were in strong enhancer regions in human lung fibroblasts and may affect gene expression in human lung tissue. Genetic mapping restricted to genome-wide enhancer single-nucleotide polymorphisms in lung fibroblasts revealed a highly significant variant (rs6763931; P = 4 × 10(-6); false discovery rate < 0.05).

CONCLUSIONS

This study offers a strategy to explore the genetic determinants of longitudinal phenotypes, provide a comprehensive picture of disease pathophysiology, and suggest potential treatment targets.

Expression QTL-based analyses reveal the mechanisms underlying colorectal cancer predisposition

Genome-wide association studies have identified many risk loci associated with colorectal cancer. Strategies integrating biological data sets with GWAS results will provide insights into the roles of risk single-nucleotide polymorphisms. We performed expression quantitative trait locus-based analyses using the information provided in The Cancer Genome Atlas. Analysis of the cis-expression quantitative trait loci (eQTLs) of 18 previously reported colorectal cancer risk loci resulted in the discovery of five variants that were significantly associated with gene expressions. Analysis of the trans-eQTLs identified three risk loci that affect the expression levels of a neighboring transcription factor, MYC. These findings provide a more comprehensive picture of gene expression determinants in colorectal cancer and insights into the underlying biology of risk loci.

MYC association with cancer risk and a new model of MYC-mediated repression

MYC is one of the most frequently mutated and overexpressed genes in human cancer but the regulation of MYC expression and the ability of MYC protein to repress cellular genes (including itself) have remained mysterious. Recent genome-wide association studies show that many genetic polymorphisms associated with disease risk map to distal regulatory elements that regulate the MYC promoter through large chromatin loops. Cancer risk-associated single-nucleotide polymorphisms (SNPs) contain more potent enhancer activity, promoting higher MYC levels and a greater risk of disease. The MYC promoter is also subject to complex regulatory circuits and limits its own expression by a feedback loop. A model for MYC autoregulation is discussed which involves a signaling pathway between the PTEN (phosphatase and tensin homolog) tumor suppressor and repressive histone modifications laid down by the EZH2 methyltransferase.

Enhancer malfunction in cancer

Why certain point mutations in a general transcription factor are associated with specific forms of cancer has been a major question in cancer biology. Enhancers are DNA regulatory elements that are key regulators of tissue-specific gene expression. Recent studies suggest that enhancer malfunction through point mutations in either regulatory elements or factors modulating enhancer-promoter communication could be the cause of tissue-specific cancer development. In this Perspective, we will discuss recent findings in the identification of cancer-related enhancer mutations and the role of Drosophila Trr and its human homologs, the MLL3 and MLL4/COMPASS-like complexes, as enhancer histone H3 lysine 4 (H3K4) monomethyltransferases functioning in enhancer-promoter communication. Recent genome-wide studies in the cataloging of somatic mutations in cancer have identified mutations in intergenic sequences encoding regulatory elements-and in MLL3 and MLL4 in both hematological malignancies and solid tumors. We propose that cancer-associated mutations in MLL3 and MLL4 exert their properties through the malfunction of Trr/MLL3/MLL4-dependent enhancers.

A comprehensive resequence-analysis of 250 kb region of 8q24.21 in men of African ancestry

BACKGROUND

Genome-wide association studies (GWAS) have identified that a ∼1 M region centromeric to the MYC oncogene on chromosome 8q24.21 harbors at least five independent loci associated with prostate cancer risk and additional loci associated with cancers of breast, colon, bladder, and chronic lymphocytic leukemia (CLL). Because GWAS identify genetic markers that may be indirectly associated with disease, fine-mapping based on sequence analysis provides important insights into patterns of linkage disequilibrium (LD) and is critical in defining the optimal variants to nominate for biological follow-up.

METHODS

To catalog variation in individuals of African ancestry, we resequenced a region (250 kb; chr8:128,050, 768–128, 300,801, hg19) containing several prostate cancer susceptibility loci as well as a locus associated with CLL. Our samples included 78 individuals from Ghana and 47 of African-Americans from Johns Hopkins University.

RESULTS

After quality control metrics were applied to next-generation sequence data, 1,838 SNPs were identified. Of these, 285 were novel and not yet reported in any public database. Using genotypes derived from sequencing, we refined the LD and recombination hotspots within the region and determined a set of tag SNPs to be used in future fine-mapping studies. Based on LD, we annotated putative risk loci and their surrogates using ENCODE data, which should help guide laboratory studies.

CONCLUSIONS

In comparison to the 1000 Genome Project data, we have identified additional variants that could be important in establishing priorities for future functional work designed to explain the biological basis of associations between SNPs and both prostate cancer and CLL.

Pleiotropic effects of genetic risk variants for other cancers on colorectal cancer risk: PAGE, GECCO and CCFR consortia

OBJECTIVE

Genome-wide association studies have identified a large number of single nucleotide polymorphisms (SNPs) associated with a wide array of cancer sites. Several of these variants demonstrate associations with multiple cancers, suggesting pleiotropic effects and shared biological mechanisms across some cancers. We hypothesised that SNPs previously associated with other cancers may additionally be associated with colorectal cancer. In a large-scale study, we examined 171 SNPs previously associated with 18 different cancers for their associations with colorectal cancer.

DESIGN

We examined 13 338 colorectal cancer cases and 40 967 controls from three consortia: Population Architecture using Genomics and Epidemiology (PAGE), Genetic Epidemiology of Colorectal Cancer (GECCO), and the Colon Cancer Family Registry (CCFR). Study-specific logistic regression results, adjusted for age, sex, principal components of genetic ancestry, and/or study specific factors (as relevant) were combined using fixed-effect meta-analyses to evaluate the association between each SNP and colorectal cancer risk. A Bonferroni-corrected p value of 2.92×10(-4) was used to determine statistical significance of the associations.

RESULTS

Two correlated SNPs--rs10090154 and rs4242382--in Region 1 of chromosome 8q24, a prostate cancer susceptibility region, demonstrated statistically significant associations with colorectal cancer risk. The most significant association was observed with rs4242382 (meta-analysis OR=1.12; 95% CI 1.07 to 1.18; p=1.74×10(-5)), which also demonstrated similar associations across racial/ethnic populations and anatomical sub-sites.

CONCLUSIONS

This is the first study to clearly demonstrate Region 1 of chromosome 8q24 as a susceptibility locus for colorectal cancer; thus, adding colorectal cancer to the list of cancer sites linked to this particular multicancer risk region at 8q24.

Dissection of thousands of cell type-specific enhancers identifies dinucleotide repeat motifs as general enhancer features

Gene expression is determined by genomic elements called enhancers, which contain short motifs bound by different transcription factors (TFs). However, how enhancer sequences and TF motifs relate to enhancer activity is unknown, and general sequence requirements for enhancers or comprehensive sets of important enhancer sequence elements have remained elusive. Here, we computationally dissect thousands of functional enhancer sequences from three different Drosophila cell lines. We find that the enhancers display distinct cis-regulatory sequence signatures, which are predictive of the enhancers’ cell type-specific or broad activities. These signatures contain transcription factor motifs and a novel class of enhancer sequence elements, dinucleotide repeat motifs (DRMs). DRMs are highly enriched in enhancers, particularly in enhancers that are broadly active across different cell types. We experimentally validate the importance of the identified TF motifs and DRMs for enhancer function and show that they can be sufficient to create an active enhancer de novo from a nonfunctional sequence. The function of DRMs as a novel class of general enhancer features that are also enriched in human regulatory regions might explain their implication in several diseases and provides important insights into gene regulation.

Frequent mutation of rs13281615 and its association with PVT1 expression and cell proliferation in breast cancer

The q24 band of chromosome 8 (8q24) is frequently amplified in human cancers including breast cancer, and several SNPs (single nucleotide polymorphisms) at 8q24, including rs13281615, have been identified for their association with cancer risks. These SNPs are in a "gene desert" region, and their functions in cancer development remain to be illustrated, although several of the SNPs appear to influence the genes in the "desert" in a long-range manner, including the v-myc avian myelocytomatosis viral oncogene homolog (MYC) and the non-protein coding plasmacytoma variant translocation 1 (PVT1), both of which have been implicated in human cancers. In the current study, we examined rs13281615 for its potential role in breast cancer using normal and cancer tissues from 121 Chinese women with breast cancer. In addition to confirming the association of the GG genotype of rs13281615 with breast cancer risk, we found that germline GG genotype was significantly associated with estrogen receptor (ER) positivity, higher tumor grade and higher proliferation index. We also found frequent somatic mutations (22/121 or 18.2%) of this SNP in breast cancer. Interestingly, the majority of the mutations (17/22 or 77%) involved a G→A change, resulting in a decrease in the number of cancers with the GG risk genotype and subsequent loss of GG association with higher tumor grade and proliferation index in cancers. Furthermore, PVT1 expression was increased in cancers, and the increase was associated with the GG genotype of rs13281615. These results suggest that the GG genotype of SNP rs13281615 plays a role in breast cancer likely by influencing PVT1 expression, and that during oncogenesis, "protective" mutations could occur.

Characterization of the transcriptional machinery bound across the widely presumed type 2 diabetes causal variant, rs7903146, within TCF7L2

Resolving the underlying functional mechanism to a given genetic association has proven extremely challenging. However, the strongest associated type 2 diabetes (T2D) locus reported to date, TCF7L2, presents an opportunity for translational analyses, as many studies in multiple ethnicities strongly point to SNP rs7903146 in intron 3 as being the causal variant within this gene. We carried out oligo pull-down combined with mass spectrophotometry (MS) to elucidate the specific transcriptional machinery across this SNP using protein extracts from HCT116 cells. We observed that poly (ADP-ribose) polymerase 1 (PARP-1) is by far the most abundant binding factor. Pursuing the possibility of a feedback mechanism, we observed that PARP-1, along with the next most abundant binding proteins, DNA topoisomerase I and ATP-dependent RNA helicase A, dimerize with the TCF7L2 protein and with each other. We uncovered further evidence of a feedback mechanism using a luciferase reporter approach, including observing expression differences between alleles for rs7903146. We also found that there was an allelic difference in the MS results for proteins with less abundant binding, namely X-ray repair cross-complementing 5 and RPA/p70. Our results point to a protein complex binding across rs7903146 within TCF7L2 and suggests a possible mechanism by which this locus confers its T2D risk.

Human colorectal cancer-specific CCAT1-L lncRNA regulates long-range chromatin interactions at the MYC locus

The human 8q24 gene desert contains multiple enhancers that form tissue-specific long-range chromatin loops with the MYC oncogene, but how chromatin looping at the MYC locus is regulated remains poorly understood. Here we demonstrate that a long noncoding RNA (lncRNA), CCAT1-L, is transcribed specifically in human colorectal cancers from a locus 515 kb upstream of MYC. This lncRNA plays a role in MYC transcriptional regulation and promotes long-range chromatin looping. Importantly, the CCAT1-L locus is located within a strong super-enhancer and is spatially close to MYC. Knockdown of CCAT1-L reduced long-range interactions between the MYC promoter and its enhancers. In addition, CCAT1-L interacts with CTCF and modulates chromatin conformation at these loop regions. These results reveal an important role of a previously unannotated lncRNA in gene regulation at the MYC locus.

Transcriptional enhancers: from properties to genome-wide predictions

Cellular development, morphology and function are governed by precise patterns of gene expression. These are established by the coordinated action of genomic regulatory elements known as enhancers or cis-regulatory modules. More than 30 years after the initial discovery of enhancers, many of their properties have been elucidated; however, despite major efforts, we only have an incomplete picture of enhancers in animal genomes. In this Review, we discuss how properties of enhancer sequences and chromatin are used to predict enhancers in genome-wide studies. We also cover recently developed high-throughput methods that allow the direct testing and identification of enhancers on the basis of their activity. Finally, we discuss recent technological advances and current challenges in the field of regulatory genomics.

Long-range interaction and correlation between MYC enhancer and oncogenic long noncoding RNA CARLo-5

The mechanism by which the 8q24 MYC enhancer region, including cancer-associated variant rs6983267, increases cancer risk is unknown due to the lack of protein-coding genes at 8q24.21. Here we report the identification of long noncoding RNAs named cancer-associated region long noncoding RNAs (CARLos) in the 8q24 region. The expression of one of the long noncoding RNAs, CARLo-5, is significantly correlated with the rs6983267 allele associated with increased cancer susceptibility. We also found the MYC enhancer region physically interacts with the active regulatory region of the CARLo-5 promoter, suggesting long-range interaction of MYC enhancer with the CARLo-5 promoter regulates CARLo-5 expression. Finally, we demonstrate that CARLo-5 has a function in cell-cycle regulation and tumor development. Overall, our data provide a key of the mystery of the 8q24 gene desert.

The myc 3' wnt-responsive element suppresses colonic tumorigenesis

Mutations in components of the Wnt/β-catenin signaling pathway are commonly found in colorectal cancers, and these mutations cause aberrant expression of genes controlled by Wnt-responsive DNA elements (WREs). While the c-Myc proto-oncogene (Myc) is required for intestinal phenotypes associated with pathogenic Wnt/β-catenin signaling in vivo, the WREs that control Myc expression in this setting have yet to be fully described. Previously, we demonstrated that the Myc 3' WRE was required for intestinal homeostasis and intestinal repair in response to damage. Here, we tested the role of the Myc 3' WRE in intestinal tumorigenesis using two independent mouse models. In comparison to Apc(Min/+) mice, Apc(Min/+) Myc 3' WRE(-/-) mice contained 25% fewer tumors in the small intestine. Deletion of the Myc 3' WRE(-/-) in the Apc(Min/+) background resulted in 4-fold more colonic tumors. In a model of colitis-associated colorectal cancer, the Myc 3' WRE suppressed colonic tumorigenesis, most notably within the cecum. Using chromatin immunoprecipitation and transcript analysis of purified colonic crypts, we found that the Myc 3' WRE is required for the transcriptional regulation of Myc expression in vivo. These results emphasize the critical role of the Myc 3' WRE in maintaining homeostatic Myc expression.

Allelic imbalance at an 8q24 oncogenic SNP is involved in activating MYC in human colorectal cancer

BACKGROUND

The rs6983267 at 8q24.21 has been established as a significant cancer-related single nucleotide polymorphism (SNP). The risk allele showed similarity to the binding site of transcription factor TCF4/LEF1 that activates transcription of MYC. However, little is known about the role of this SNP in increasing MYC activity in colorectal cancers (CRCs).

METHODS

The genotypes of rs6983267 in peripheral blood and primary cancers, MYC activity and copy number (CN) alteration were examined in 107 CRCs. Next, we plotted the number of cancers cell lines exhibiting specific G/T genotypes in 746 cancer cell lines of the Sanger Institute database. Then we validated the relationship between the 8q24 SNP status and clinicopathologic parameters in 68 CRCs with loss of heterozygosity (LOH).

RESULTS

The MYC module activity was activated by either transcription in the risk allele (G) or by amplification in the non-risk allele (T). Then, we confirmed that the CN amplification dominantly occurred in the non-risk allele, whereas CN neutral LOH, which indicated uniparental disomy (UPD) was more frequently observed for the risk allele. Finally, we confirmed that risk allele dominant cases, either by amplification or by UPD, indicated a more malignant clinical phenotype than non-risk allele dominant cases.

CONCLUSIONS

The development of CRC requires MYC activation through retention of the risk allele, or amplification of the non-risk allele at the oncogenic SNP in the site of primary tumor.

An overview of MYC and its interactome

This review is intended to provide a broad outline of the biological and molecular functions of MYC as well as of the larger protein network within which MYC operates. We present a view of MYC as a sensor that integrates multiple cellular signals to mediate a broad transcriptional response controlling many aspects of cell behavior. We also describe the larger transcriptional network linked to MYC with emphasis on the MXD family of MYC antagonists. Last, we discuss evidence that the network has evolved for millions of years, dating back to the emergence of animals.

Targeted knock-in of the polymorphism rs61764370 does not affect KRAS expression but reduces let-7 levels

Understanding the role of single-nucleotide polymorphisms (SNPs) in the pathological process represents a unique experimental challenge especially when the variants occur outside of coding regions. The noncoding SNP rs61764370 located in the 3'-untranslated region of Kirsten rat sarcoma viral oncogene homolog (KRAS) has been implicated as a risk factor for the development of cancer and the response to targeted therapies. This cancer-associated variant is thought to affect the binding of the microRNA let-7, which allegedly modulates KRAS expression. Using site-specific homologous recombination, we inserted the rs61764370:T>G KRAS gene variant in the colorectal cancer cell line SW48 (SW48 +SNP) and assessed the cellular and biochemical phenotype. We observed a significant increase in cellular proliferation, as well as a reduction in the levels of the microRNA let-7a, let-7b, and let-7c. Transcriptional and biochemical analysis showed no concomitant change in the KRAS protein expression or modulation of the downstream mitogen activated kinase or PI3K/AKT signaling. These results suggest that the cancer-associated rs61764370 variant exerts a biological effect not through transcriptional modulation of KRAS but rather by tuning the expression of the microRNA let-7.

From GWAS to function: genetic variation in sodium channel gene enhancer influences electrical patterning

The electrical activity of the heart depends on the correct interplay between key transcription factors and cis-regulatory elements, which together regulate the proper heterogeneous expression of genes encoding for ion channels and other proteins. Genome-wide association studies of ECG parameters implicated genetic variants in the genes for these factors and ion channels modulating conduction and depolarization. Here, we review recent insights into the regulation of localized expression of ion channel genes and the mechanism by which a single-nucleotide polymorphism (SNP) associated with alterations in cardiac conduction patterns in humans affects the transcriptional regulation of the sodium channel genes, SCN5A and SCN10A. The identification of regulatory elements of electrical activity genes helps to explain the impact of genetic variants in non-coding regulatory DNA sequences on regulation of cardiac conduction and the predisposition for cardiac arrhythmias.