8q24 prostate, breast, and colon cancer risk loci show tissue-specific long-range interaction with MYC

4C-seq revealed long-range interactions of a functional enhancer at the 8q24 prostate cancer risk locus

Genome-wide association studies (GWAS) have identified >100 independent susceptibility loci for prostate cancer, including the hot spot at 8q24. However, how genetic variants at this locus confer disease risk hasn’t been fully characterized. Using circularized chromosome conformation capture (4C) coupled with next-generation sequencing and an enhancer at 8q24 as “bait”, we identified genome-wide partners interacting with this enhancer in cell lines LNCaP and C4-2B. These 4C-identified regions are distributed in open nuclear compartments, featuring active histone marks (H3K4me1, H3K4me2 and H3K27Ac). Transcription factors NKX3-1, FOXA1 and AR (androgen receptor) tend to occupy these 4C regions. We identified genes located at the interacting regions, and found them linked to positive regulation of mesenchymal cell proliferation in LNCaP and C4-2B, and several pathways (TGF beta signaling pathway in LNCaP and p53 pathway in C4-2B). Common genes (e.g. MYC and POU5F1B) were identified in both prostate cancer cell lines. However, each cell line also had exclusive genes (e.g. ELAC2 and PTEN in LNCaP and BRCA2 and ZFHX3 in C4-2B). In addition, BCL-2 identified in C4-2B might contribute to the progression of androgen-refractory prostate cancer. Overall, our work reveals key genes and pathways involved in prostate cancer onset and progression.

Germline genetic profiling in prostate cancer: latest developments and potential clinical applications

Familial and twin studies have demonstrated a significant inherited component to prostate cancer predisposition. Genome wide association studies have shown that there are 100 single nucleotide polymorphisms which have been associated with the development of prostate cancer. This review aims to discuss the scientific methods used to identify these susceptibility loci. It will also examine the current clinical utility of these loci, which include the development of risk models as well as predicting treatment efficacy and toxicity. In order to refine the clinical utility of the susceptibility loci, international consortia have been developed to combine statistical power as well as skills and knowledge to further develop models that could be used to predict risk and treatment outcomes.

Molecular Discriminators of Racial Disparities in Prostate Cancer

Recent molecular characterization of prostate cancer (PCa) identified novel genetic aberrations and disease subtypes. The frequencies of molecular aberrations show racial disparity. Clinical strategies and targeted therapies embracing these racial differences are required. Here we discuss ethnic differences in genetic alterations and their impact on the susceptibility, progression, and treatment of prostate cancer.

Prostate Cancer Susceptibility in Men of African Ancestry at 8q24

The 8q24 region harbors multiple risk variants for distinct cancers, including >8 for prostate cancer. In this study, we conducted fine mapping of the 8q24 risk region (127.8-128.8Mb) in search of novel associations with common and rare variation in 4853 prostate cancer case patients and 4678 control subjects of African ancestry. All statistical tests were two-sided. We identified three independent associations at P values of less than 5.00×10(-8), all of which were replicated in studies from Ghana and Uganda (combined sample = 5869 case patients, 5615 control subjects; rs114798100: risk allele frequency [RAF] = 0.04, per-allele odds ratio [OR] = 2.31, 95% confidence interval [CI] = 2.04 to 2.61, P = 2.38×10(-40); rs72725879: RAF = 0.33, OR = 1.37, 95% CI = 1.30 to 1.45, P = 3.04×10(-27); and rs111906932: RAF = 0.03, OR = 1.79, 95% CI = 1.53 to 2.08, P = 1.39×10(-13)). Risk variants rs114798100 and rs111906923 are only found in men of African ancestry, with rs111906923 representing a novel association signal. The three variants are located within or near a number of prostate cancer-associated long noncoding RNAs (lncRNAs), including PRNCR1, PCAT1, and PCAT2. These findings highlight ancestry-specific risk variation and implicate prostate-specific lncRNAs at the 8q24 prostate cancer susceptibility region.

Effects on the transcriptome upon deletion of a distal element cannot be predicted by the size of the H3K27Ac peak in human cells

Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with increased risk for colorectal cancer (CRC). A molecular understanding of the functional consequences of this genetic variation is complicated because most GWAS SNPs are located in non-coding regions. We used epigenomic information to identify H3K27Ac peaks in HCT116 colon cancer cells that harbor SNPs associated with an increased risk for CRC. Employing CRISPR/Cas9 nucleases, we deleted a CRC risk-associated H3K27Ac peak from HCT116 cells and observed large-scale changes in gene expression, resulting in decreased expression of many nearby genes. As a comparison, we showed that deletion of a robust H3K27Ac peak not associated with CRC had minimal effects on the transcriptome. Interestingly, although there is no H3K27Ac peak in HEK293 cells in the E7 region, deletion of this region in HEK293 cells decreased expression of several of the same genes that were downregulated in HCT116 cells, including the MYC oncogene. Accordingly, deletion of E7 causes changes in cell culture assays in HCT116 and HEK293 cells. In summary, we show that effects on the transcriptome upon deletion of a distal regulatory element cannot be predicted by the size or presence of an H3K27Ac peak.

Prostate cancer

Much progress has been made in research for prostate cancer in the past decade. There is now greater understanding for the genetic basis of familial prostate cancer with identification of rare but high-risk mutations (eg, BRCA2, HOXB13) and low-risk but common alleles (77 identified so far by genome-wide association studies) that could lead to targeted screening of patients at risk. This is especially important because screening for prostate cancer based on prostate-specific antigen remains controversial due to the high rate of overdiagnosis and unnecessary prostate biopsies, despite evidence that it reduces mortality. Classification of prostate cancer into distinct molecular subtypes, including mutually exclusive ETS-gene-fusion-positive and SPINK1-overexpressing, CHD1-loss cancers, could allow stratification of patients for different management strategies. Presently, men with localised disease can have very different prognoses and treatment options, ranging from observation alone through to radical surgery, with few good-quality randomised trials to inform on the best approach for an individual patient. The survival of patients with metastatic prostate cancer progressing on androgen-deprivation therapy (castration-resistant prostate cancer) has improved substantially. In addition to docetaxel, which has been used for more than a decade, in the past 4 years five new drugs have shown efficacy with improvements in overall survival leading to licensing for the treatment of metastatic castration-resistant prostate cancer. Because of this rapid change in the therapeutic landscape, no robust data exist to inform on the selection of patients for a specific treatment for castration-resistant prostate cancer or the best sequence of administration. Moreover, the high cost of the newer drugs limits their widespread use in several countries. Data from continuing clinical and translational research are urgently needed to improve, and, crucially, to personalise management.

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 identification of trans-associations between prostate cancer GWAS SNPs and RNA expression differences in tumor-adjacent stroma

Here we tested the hypothesis that SNPs associated with prostate cancer risk, might differentially affect RNA expression in prostate cancer stroma. The most significant 35 SNP loci were selected from Genome Wide Association (GWA) studies of ~40,000 patients. We also selected 4030 transcripts previously associated with prostate cancer diagnosis and prognosis. eQTL analysis was carried out by a modified BAYES method to analyze the associations between the risk variants and expressed transcripts jointly in a single model. We observed 47 significant associations between eight risk variants and the expression patterns of 46 genes. This is the first study to identify associations between multiple SNPs and multiple in trans gene expression differences in cancer stroma. Potentially, a combination of SNPs and associated expression differences in prostate stroma may increase the power of risk assessment for individuals, and for cancer progression.

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.

Molecular Updates in Prostate Cancer

A wide array of molecular markers and genomic signatures, reviewed in this article, may soon be used as adjuncts to currently established screening strategies, prognostic parameters, and early detection markers. Markers of genetic susceptibility to PCA, recurrent epigenetic and genetic alterations, including ETS gene fusions, PTEN alterations, and urine-based early detection marker PCA3, are discussed. Impact of recent genome-wide assessment on our understanding of key pathways of PCA development and progression and their potential clinical implications are highlighted.

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.

Integration of multiethnic fine-mapping and genomic annotation to prioritize candidate functional SNPs at prostate cancer susceptibility regions

Interpretation of biological mechanisms underlying genetic risk associations for prostate cancer is complicated by the relatively large number of risk variants (n = 100) and the thousands of surrogate SNPs in linkage disequilibrium. Here, we combined three distinct approaches: multiethnic fine-mapping, putative functional annotation (based upon epigenetic data and genome-encoded features), and expression quantitative trait loci (eQTL) analyses, in an attempt to reduce this complexity. We examined 67 risk regions using genotyping and imputation-based fine-mapping in populations of European (cases/controls: 8600/6946), African (cases/controls: 5327/5136), Japanese (cases/controls: 2563/4391) and Latino (cases/controls: 1034/1046) ancestry. Markers at 55 regions passed a region-specific significance threshold (P-value cutoff range: 3.9 × 10(-4)-5.6 × 10(-3)) and in 30 regions we identified markers that were more significantly associated with risk than the previously reported variants in the multiethnic sample. Novel secondary signals (P < 5.0 × 10(-6)) were also detected in two regions (rs13062436/3q21 and rs17181170/3p12). Among 666 variants in the 55 regions with P-values within one order of magnitude of the most-associated marker, 193 variants (29%) in 48 regions overlapped with epigenetic or other putative functional marks. In 11 of the 55 regions, cis-eQTLs were detected with nearby genes. For 12 of the 55 regions (22%), the most significant region-specific, prostate-cancer associated variant represented the strongest candidate functional variant based on our annotations; the number of regions increased to 20 (36%) and 27 (49%) when examining the 2 and 3 most significantly associated variants in each region, respectively. These results have prioritized subsets of candidate variants for downstream functional evaluation.

Cis-eQTL analysis and functional validation of candidate susceptibility genes for high-grade serous ovarian cancer

Genome-wide association studies have reported 11 regions conferring risk of high-grade serous epithelial ovarian cancer (HGSOC). Expression quantitative trait locus (eQTL) analyses can identify candidate susceptibility genes at risk loci. Here we evaluate cis-eQTL associations at 47 regions associated with HGSOC risk (P≤10(-5)). For three cis-eQTL associations (P<1.4 × 10(-3), FDR<0.05) at 1p36 (CDC42), 1p34 (CDCA8) and 2q31 (HOXD9), we evaluate the functional role of each candidate by perturbing expression of each gene in HGSOC precursor cells. Overexpression of HOXD9 increases anchorage-independent growth, shortens population-doubling time and reduces contact inhibition. Chromosome conformation capture identifies an interaction between rs2857532 and the HOXD9 promoter, suggesting this SNP is a leading causal variant. Transcriptomic profiling after HOXD9 overexpression reveals enrichment of HGSOC risk variants within HOXD9 target genes (P=6 × 10(-10) for risk variants (P<10(-4)) within 10 kb of a HOXD9 target gene in ovarian cells), suggesting a broader role for this network in genetic susceptibility to HGSOC.

MYC, Metabolism, and Cancer

The MYC oncogene encodes a transcription factor, MYC, whose broad effects make its precise oncogenic role enigmatically elusive. The evidence to date suggests that MYC triggers selective gene expression amplification to promote cell growth and proliferation. Through its targets, MYC coordinates nutrient acquisition to produce ATP and key cellular building blocks that increase cell mass and trigger DNA replication and cell division. In cancer, genetic and epigenetic derangements silence checkpoints and unleash MYC's cell growth- and proliferation-promoting metabolic activities. Unbridled growth in response to deregulated MYC expression creates dependence on MYC-driven metabolic pathways, such that reliance on specific metabolic enzymes provides novel targets for cancer therapy.

SIGNIFICANCE

MYC's expression and activity are tightly regulated in normal cells by multiple mechanisms, including a dependence upon growth factor stimulation and replete nutrient status. In cancer, genetic deregulation of MYC expression and loss of checkpoint components, such as TP53, permit MYC to drive malignant transformation. However, because of the reliance of MYC-driven cancers on specific metabolic pathways, synthetic lethal interactions between MYC overexpression and specific enzyme inhibitors provide novel cancer therapeutic opportunities.

Analysis of Prostate Cancer Susceptibility Variants in South African Men: Replicating Associations on Chromosomes 8q24 and 10q11

Genome-wide association studies (GWAS) have implicated single nucleotide polymorphisms (SNPs) on chromosomes 2p15, 6q25, 7p15.2, 7q21, 8q24, 10q11, 10q26, 11q13, 17q12, 17q24, 19q13, and Xp11, with prostate cancer (PCa) susceptibility and/or tumour aggressiveness, in populations of African, European, and Asian ancestry. The objective of this study was to confirm these associations in South African Mixed Ancestry and White men. We evaluated 17 prioritised GWAS SNPs in South African cases (331 Mixed Ancestry and 155 White) and controls (178 Mixed Ancestry and 145 White). The replicated SNP associations for the different South African ethnic groups were rs7008482 (8q24) (p = 2.45 × 10⁻⁵), rs6983267 (8q24) (p = 4.48 × 10⁻⁷), and rs10993994 (10q11) (p = 1.40 × 10⁻³) in Mixed Ancestry men and rs10993994 (p = 1.56 × 10⁻⁹) in White men. No significant associations were observed for the analyses stratified by disease aggressiveness in the individual and the combined population group analysis. The present study demonstrates that a number of known PCa susceptibility variants may contribute to disease susceptibility in South African men. Larger genetic investigations extended to other South African population groups are warranted to confirm the role of these and other SNPs in disease susceptibility.

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.

The Evolution of Bony Vertebrate Enhancers at Odds with Their Coding Sequence Landscape

Enhancers lie at the heart of transcriptional and developmental gene regulation. Therefore, changes in enhancer sequences usually disrupt the target gene expression and result in disease phenotypes. Despite the well-established role of enhancers in development and disease, evolutionary sequence studies are lacking. The current study attempts to unravel the puzzle of bony vertebrates’ conserved noncoding elements (CNE) enhancer evolution. Bayesian phylogenetics of enhancer sequences spotlights promising interordinal relationships among placental mammals, proposing a closer relationship between humans and laurasiatherians while placing rodents at the basal position. Clock-based estimates of enhancer evolution provided a dynamic picture of interspecific rate changes across the bony vertebrate lineage. Moreover, coelacanth in the study augmented our appreciation of the vertebrate cis-regulatory evolution during water–land transition. Intriguingly, we observed a pronounced upsurge in enhancer evolution in land-dwelling vertebrates. These novel findings triggered us to further investigate the evolutionary trend of coding as well as CNE nonenhancer repertoires, to highlight the relative evolutionary dynamics of diverse genomic landscapes. Surprisingly, the evolutionary rates of enhancer sequences were clearly at odds with those of the coding and the CNE nonenhancer sequences during vertebrate adaptation to land, with land vertebrates exhibiting significantly reduced rates of coding sequence evolution in comparison to their fast evolving regulatory landscape. The observed variation in tetrapod cis-regulatory elements caused the fine-tuning of associated gene regulatory networks. Therefore, the increased evolutionary rate of tetrapods’ enhancer sequences might be responsible for the variation in developmental regulatory circuits during the process of vertebrate adaptation to land.

Detection of Enhancer-Associated Rearrangements Reveals Mechanisms of Oncogene Dysregulation in B-cell Lymphoma

B-cell lymphomas frequently contain genomic rearrangements that lead to oncogene activation by heterologous distal regulatory elements. We used a novel approach called "pinpointing enhancer-associated rearrangements by chromatin immunoprecipitation," or PEAR-ChIP, to simultaneously map enhancer activity and proximal rearrangements in lymphoma cell lines and patient biopsies. This method detects rearrangements involving known cancer genes, including CCND1, BCL2, MYC, PDCD1LG2, NOTCH1, CIITA, and SGK1, as well as novel enhancer duplication events of likely oncogenic significance. We identify lymphoma subtype-specific enhancers in the MYC locus that are silenced in lymphomas with MYC-activating rearrangements and are associated with germline polymorphisms that alter lymphoma risk. We show that BCL6-locus enhancers are acetylated by the BCL6-activating transcription factor MEF2B, and can undergo genomic duplication, or target the MYC promoter for activation in the context of a "pseudo-double-hit" t(3;8)(q27;q24) rearrangement linking the BCL6 and MYC loci. Our work provides novel insights regarding enhancer-driven oncogene activation in lymphoma.

SIGNIFICANCE

We demonstrate a novel approach for simultaneous detection of genomic rearrangements and enhancer activity in tumor biopsies. We identify novel mechanisms of enhancer-driven regulation of the oncogenes MYC and BCL6, and show that the BCL6 locus can serve as an enhancer donor in an "enhancer hijacking" translocation.

Association of three 8q24 polymorphisms with prostate cancer susceptibility: evidence from a meta-analysis with 50,854 subjects

The 8q24 polymorphisms have been implicated in various cancers. Three 8q24 polymorphisms (rs1447295 C>A, rs16901979 C>A, and rs6983267 T>G) have been extensively investigated for their association with prostate cancer (PCa) susceptibility, yet conclusions are contradictory. We conducted a comprehensive meta-analysis to reevaluate the associations between those polymorphisms and PCa susceptibility, according to the latest meta-analysis guidelines (PRISMA). Eligible publications were searched from MEDLINE, EMBASE and CBM. False positive report possibility analysis was performed. We totally collected 20184 cases and 20439 controls from 20 studies for the rs1447295 C>A, 1850 cases and 2090 controls from 7 studies for the rs16901979 C>A, and 12233 cases and 7582 controls from 17 studies for the rs6983267 T>G. Overall, each of studied 8q24 polymorphisms was significantly associated with PCa risk individually. Significant associations were also observed in stratified analysis by ethnicity, source of control, and quality score. Interestingly, the effect of rs1447295 on PCa risk was observed among Caucasians and Asians, but not Africa-Americans. The effect of rs16901979 was more prominent among Africa-Americans than Asians. Likewise, rs6983267 conferred a higher Pca risk among Caucasians than Asians. Collectively, these 8q24 variant(s) may modulate PCa risk in an ethnic-specific manner.

Variants on 8q24 and prostate cancer risk in Chinese population: a meta-analysis

Previous studies have identified 8q24 as an important region to prostate cancer (PCa) susceptibility. The aim of this study was to investigate the role of six genetic variants on 8q24 (rs1447295, A; rs6983267, G; rs6983561, C; rs7837688, T; rs10090154, T and rs16901979, A) on PCa risk in Chinese population. Online electronic databases were searched to retrieve related articles concerning the association between 8q24 variants and PCa risk in men of Chinese population published between 2000 and 2014. Odds ratio (ORs) with its 95% correspondence interval (CI) were employed to assess the strength of association. Total eleven case-control studies were screened out, including 2624 PCa patients and 2438 healthy controls. Our results showed that three risk alleles of rs1447295 A (OR=1.35, 95% CI=1.19-1.53, P<0.00001), rs6983561 C (C vs. A: OR=1.41, 95% CI=1.21-1.63, P<0.00001) and rs10090154 T (T vs. C: OR=1.48, 95% CI=1.22-1.80, P<0.00001) on8q24 were significantly associated with PCa risk in Chinese population. Furthermore, genotypes of rs1447295, AA+AC; rs6983561, CC+AC and CC; rs10090154, TT+TC; and rs16901979, AA were associated with PCa as well (P<0.01). No association was found between rs6983267, rs7837688 and PCa risk. In conclusions, variants including rs1447295, rs6983561, rs10090154 and rs16901979 on 8q24 might be associated with PCa risk in Chinese population, indicating these four variations may contribute risk to this disease. This meta-analysis was the first study to assess the role of 8q24 variants on PCa risk in Chinese population.

Genetic variation at 8q24, family history of cancer, and upper gastrointestinal cancers in a Chinese population

Genetic variation at 8q24 is associated with prostate, bladder, breast, colorectal, thyroid, lung, ovarian, UADT, liver and stomach cancers. However, a role for variation at 8q24 in familial clustering of upper gastrointestinal cancers has not been studied. In order to explore potential inherited susceptibility, we analyzed epidemiologic data from a population-based case-control study of upper gastrointestinal cancers from Taixing, China. The study population includes 204 liver, 206 stomach, and 218 esophageal cancer cases and 415 controls. Associations between 8q24 rs1447295, rs16901979, rs6983267 and these cancers were stratified by family history of cancer. Odds ratios and 95% confidence intervals were adjusted for potential confounders: age, sex, education, tobacco smoking, alcohol consumption, and BMI at interview. We also adjusted for hepatitis B and aflatoxin (liver cancer) and Helicobacter pylori (stomach cancer). In a dominant model, among those with a family history of cancer, rs1447295 was positively associated with liver cancer (OR(adj) 2.80; 95% CI 1.15-6.80). Heterogeneity was observed (P(heterogeneity) = 0.029) with rs6983267 and liver cancer, with positive association in the dominant model among those with a family history of cancer and positive association in the recessive model among those without a family history of cancer. When considered in a genetic risk score model, each additional 8q24 risk genotype increased the odds of liver cancer by two-fold among those with a family history of cancer (OR(adj) 2.00; 95% CI 1.15-3.47). These findings suggest that inherited susceptibility to liver cancer may exist in the Taixing population and that variation at 8q24 might be a genetic component of that inherited susceptibility.

Association of common variants on chromosome 8q24 with gastric cancer in Venezuelan patients

Gastric cancer remains one of the leading causes of death in the world, being Central and South America among the regions showing the highest incidence and mortality rates worldwide. Although several single nucleotide polymorphisms (SNPs) identified in the chromosomal region 8q24 by genome-wide association studies have been related with the risk of different kinds of cancers, their role in the susceptibility of gastric cancer in Latin American populations has not been evaluated yet. Hereby, we performed a case-control study to explore the associations between three SNPs at 8q24 and gastric cancer risk in Venezuelan patients. We analyzed rs1447295, rs4733616 and rs6983267 SNPs in 122 paraffin-embedded tumor samples from archival bank and 129 samples with chronic gastritis (obtained by upper endoscopy during the study) from the Central Hospital of Barquisimeto (Lara, Venezuela). Genotypes were determined by PCR-RFLP reactions designed in this study for efficient genotyping of formalin-fixed/paraffin-embedded tissues. No significant differences in genotype frequencies between case and control groups were found. However, carriers of the homozygous TT genotype of SNP rs4733616 had an increased risk of developing poorly differentiated gastric cancer according to the codominant (OR=3.59, P=0.035) and the recessive models (OR=4.32, P=0.014, best-fitting model of inheritance), adjusted by age and gender. Our study suggests that the SNP rs4733616 is associated with susceptibility to poorly differentiated gastric cancer in Venezuelans. Additional studies are needed to further interrogate the prognostic value of the rs4733616 marker in this high-risk population for gastric cancer.

Capture Hi-C identifies the chromatin interactome of colorectal cancer risk loci

Multiple regulatory elements distant from their targets on the linear genome can influence the expression of a single gene through chromatin looping. Chromosome conformation capture implemented in Hi-C allows for genome-wide agnostic characterization of chromatin contacts. However, detection of functional enhancer-promoter interactions is precluded by its effective resolution that is determined by both restriction fragmentation and sensitivity of the experiment. Here we develop a capture Hi-C (cHi-C) approach to allow an agnostic characterization of these physical interactions on a genome-wide scale. Single-nucleotide polymorphisms associated with complex diseases often reside within regulatory elements and exert effects through long-range regulation of gene expression. Applying this cHi-C approach to 14 colorectal cancer risk loci allows us to identify key long-range chromatin interactions in cis and trans involving these loci.

Genetic study of complex diseases in the post-GWAS era

Genome-wide association studies (GWASs) have identified thousands of genes and genetic variants (mainly SNPs) that contribute to complex diseases in humans. Functional characterization and mechanistic elucidation of these SNPs and genes action are the next major challenge. It has been well established that SNPs altering the amino acids of protein-coding genes can drastically impact protein function, and play an important role in molecular pathogenesis. Functions of regulatory SNPs can be complex and elusive, and involve gene expression regulation through the effect on RNA splicing, transcription factor binding, DNA methylation and miRNA recruitment. In the present review, we summarize the recent progress in our understanding of functional consequences of GWAS-associated non-coding regulatory SNPs, and discuss the application of systems genetics and network biology in the interpretation of GWAS findings.

In silico functional pathway annotation of 86 established prostate cancer risk variants

Heritability is one of the strongest risk factors of prostate cancer, emphasizing the importance of the genetic contribution towards prostate cancer risk. To date, 86 established prostate cancer risk variants have been identified by genome-wide association studies (GWAS). To determine if these risk variants are located near genes that interact together in biological networks or pathways contributing to prostate cancer initiation or progression, we generated gene sets based on proximity to the 86 prostate cancer risk variants. We took two approaches to generate gene lists. The first strategy included all immediate flanking genes, up- and downstream of the risk variant, regardless of distance from the index variant, and the second strategy included genes closest to the index GWAS marker and to variants in high LD (r2 ≥0.8 in Europeans) with the index variant, within a 100 kb window up- and downstream. Pathway mapping of the two gene sets supported the importance of the androgen receptor-mediated signaling in prostate cancer biology. In addition, the hedgehog and Wnt/β-catenin signaling pathways were identified in pathway mapping for the flanking gene set. We also used the HaploReg resource to examine the 86 risk loci and variants high LD (r2 ≥0.8) for functional elements. We found that there was a 12.8 fold (p = 2.9 x 10-4) enrichment for enhancer motifs in a stem cell line and a 4.4 fold (p = 1.1 x 10-3) enrichment of DNase hypersensitivity in a prostate adenocarcinoma cell line, indicating that the risk and correlated variants are enriched for transcriptional regulatory motifs. Our pathway-based functional annotation of the prostate cancer risk variants highlights the potential regulatory function that GWAS risk markers, and their highly correlated variants, exert on genes. Our study also shows that these genes may function cooperatively in key signaling pathways in prostate cancer biology.

A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping

We use in situ Hi-C to probe the 3D architecture of genomes, constructing haploid and diploid maps of nine cell types. The densest, in human lymphoblastoid cells, contains 4.9 billion contacts, achieving 1 kb resolution. We find that genomes are partitioned into contact domains (median length, 185 kb), which are associated with distinct patterns of histone marks and segregate into six subcompartments. We identify ∼10,000 loops. These loops frequently link promoters and enhancers, correlate with gene activation, and show conservation across cell types and species. Loop anchors typically occur at domain boundaries and bind CTCF. CTCF sites at loop anchors occur predominantly (>90%) in a convergent orientation, with the asymmetric motifs "facing" one another. The inactive X chromosome splits into two massive domains and contains large loops anchored at CTCF-binding repeats.

Unbiased analysis of potential targets of breast cancer susceptibility loci by Capture Hi-C

Genome-wide association studies have identified more than 70 common variants that are associated with breast cancer risk. Most of these variants map to non-protein-coding regions and several map to gene deserts, regions of several hundred kilobases lacking protein-coding genes. We hypothesized that gene deserts harbor long-range regulatory elements that can physically interact with target genes to influence their expression. To test this, we developed Capture Hi-C (CHi-C), which, by incorporating a sequence capture step into a Hi-C protocol, allows high-resolution analysis of targeted regions of the genome. We used CHi-C to investigate long-range interactions at three breast cancer gene deserts mapping to 2q35, 8q24.21, and 9q31.2. We identified interaction peaks between putative regulatory elements ("bait fragments") within the captured regions and "targets" that included both protein-coding genes and long noncoding (lnc) RNAs over distances of 6.6 kb to 2.6 Mb. Target protein-coding genes were IGFBP5, KLF4, NSMCE2, and MYC; and target lncRNAs included DIRC3, PVT1, and CCDC26. For one gene desert, we were able to define two SNPs (rs12613955 and rs4442975) that were highly correlated with the published risk variant and that mapped within the bait end of an interaction peak. In vivo ChIP-qPCR data show that one of these, rs4442975, affects the binding of FOXA1 and implicate this SNP as a putative functional variant.

Genetic variation in regulatory DNA elements: the case of OCA2 transcriptional regulation

Mutations within the OCA2 gene or the complete absence of the OCA2 protein leads to oculocutaneous albinism type 2. The OCA2 protein plays a central role in melanosome biogenesis, and it is a strong determinant of the eumelanin content in melanocytes. Transcript levels of the OCA2 gene are strongly correlated with pigmentation intensities. Recent studies demonstrated that the transcriptional level of OCA2 is to a large extent determined by the noncoding SNP rs12913832 located 21.5 kb upstream of the OCA2 gene promoter. In this review, we discuss current hypotheses and the available data on the mechanism of OCA2 transcriptional regulation and how this is influenced by genetic variation. Finally, we will explore how future epigenetic studies can be used to advance our insight into the functional biology that connects genetic variation to human pigmentation.

DNA methylation levels at chromosome 8q24 in peripheral blood are associated with 8q24 cancer susceptibility loci

Chromosome 8q24 has emerged as an important region for genetic susceptibility to various cancers, but little is known about the contribution of DNA methylation at 8q24. To evaluate variability in DNA methylation levels at 8q24 and the relationship with cancer susceptibility single nucleotide polymorphisms (SNPs) in this region, we quantified DNA methylation levels in peripheral blood at 145 CpG sites nearby 8q24 cancer susceptibility SNPs or MYC using pyrosequencing among 80 Caucasian men in the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. For the 60 CpG sites meeting quality control, which also demonstrated temporal stability over a 5-year period, we calculated pairwise Spearman correlations for DNA methylation levels at each CpG site with 42 8q24 cancer susceptibility SNPs. To account for multiple testing, we adjusted P values into q values reflecting the false discovery rate (FDR). In contrast to the MYC CpG sites, most sites nearby the SNPs demonstrated good reproducibility, high methylation levels, and moderate-high between-individual variation. We observed 10 statistically significant (FDR < 0.05) CpG site-SNP correlations. These included correlations between an intergenic CpG site at Chr8:128393157 and the prostate cancer SNP rs16902094 (ρ = -0.54; P = 9.7 × 10(-7); q = 0.002), a PRNCR1 CpG site at Chr8:128167809 and the prostate cancer SNP rs1456315 (ρ = 0.52; P = 1.4 × 10(-6); q = 0.002), and two POU5F1B CpG sites and several prostate/colorectal cancer SNPs (for Chr8:128498051 and rs6983267, ρ = 0.46; P = 2.0 × 10(-5); q = 0.01). This is the first report of correlations between blood DNA methylation levels and cancer susceptibility SNPs at 8q24, suggesting that DNA methylation at this important susceptibility locus may contribute to cancer risk.

Enhancer alterations in cancer: a source for a cell identity crisis

Enhancers are selectively utilized to orchestrate gene expression programs that first govern pluripotency and then proceed to highly specialized programs required for the process of cellular differentiation. Whereas gene-proximal promoters are typically active across numerous cell types, distal enhancer activation is cell-type-specific and central to cell fate determination, thereby accounting for cell identity. Recent studies have highlighted the diversity of enhancer usage, cataloguing millions of such elements in the human genome. The disruption of enhancer activity, through genetic or epigenetic alterations, can impact cell-type-specific functions, resulting in a wide range of pathologies. In cancer, these alterations can promote a 'cell identity crisis', in which enhancers associated with oncogenes and multipotentiality are activated, while those promoting cell fate commitment are inactivated. Overall, these alterations favor an undifferentiated cellular phenotype. Here, we review the current knowledge regarding the role of enhancers in normal cell function, and discuss how genetic and epigenetic changes in enhancer elements potentiate oncogenesis. In addition, we discuss how understanding the mechanisms regulating enhancer activity can inform therapeutic opportunities in cancer cells and highlight key challenges that remain in understanding enhancer biology as it relates to oncology.

Functional annotation of putative regulatory elements at cancer susceptibility Loci

Most cancer-associated genetic variants identified from genome-wide association studies (GWAS) do not obviously change protein structure, leading to the hypothesis that the associations are attributable to regulatory polymorphisms. Translating genetic associations into mechanistic insights can be facilitated by knowledge of the causal regulatory variant (or variants) responsible for the statistical signal. Experimental validation of candidate functional variants is onerous, making bioinformatic approaches necessary to prioritize candidates for laboratory analysis. Thus, a systematic approach for recognizing functional (and, therefore, likely causal) variants in noncoding regions is an important step toward interpreting cancer risk loci. This review provides a detailed introduction to current regulatory variant annotations, followed by an overview of how to leverage these resources to prioritize candidate functional polymorphisms in regulatory regions.

Genome-wide association study identifies multiple susceptibility loci for pancreatic cancer

We performed a multistage genome-wide association study including 7,683 individuals with pancreatic cancer and 14,397 controls of European descent. Four new loci reached genome-wide significance: rs6971499 at 7q32.3 (LINC-PINT, per-allele odds ratio (OR) = 0.79, 95% confidence interval (CI) 0.74-0.84, P = 3.0 × 10(-12)), rs7190458 at 16q23.1 (BCAR1/CTRB1/CTRB2, OR = 1.46, 95% CI 1.30-1.65, P = 1.1 × 10(-10)), rs9581943 at 13q12.2 (PDX1, OR = 1.15, 95% CI 1.10-1.20, P = 2.4 × 10(-9)) and rs16986825 at 22q12.1 (ZNRF3, OR = 1.18, 95% CI 1.12-1.25, P = 1.2 × 10(-8)). We identified an independent signal in exon 2 of TERT at the established region 5p15.33 (rs2736098, OR = 0.80, 95% CI 0.76-0.85, P = 9.8 × 10(-14)). We also identified a locus at 8q24.21 (rs1561927, P = 1.3 × 10(-7)) that approached genome-wide significance located 455 kb telomeric of PVT1. Our study identified multiple new susceptibility alleles for pancreatic cancer that are worthy of follow-up studies.

Prostate cancer risk locus at 8q24 as a regulatory hub by physical interactions with multiple genomic loci across the genome

Chromosome 8q24 locus contains regulatory variants that modulate genetic risk to various cancers including prostate cancer (PC). However, the biological mechanism underlying this regulation is not well understood. Here, we developed a chromosome conformation capture (3C)-based multi-target sequencing technology and systematically examined three PC risk regions at the 8q24 locus and their potential regulatory targets across human genome in six cell lines. We observed frequent physical contacts of this risk locus with multiple genomic regions, in particular, inter-chromosomal interaction with CD96 at 3q13 and intra-chromosomal interaction with MYC at 8q24. We identified at least five interaction hot spots within the predicted functional regulatory elements at the 8q24 risk locus. We also found intra-chromosomal interaction genes PVT1, FAM84B and GSDMC and inter-chromosomal interaction gene CXorf36 in most of the six cell lines. Other gene regions appeared to be cell line-specific, such as RRP12 in LNCaP, USP14 in DU-145 and SMIN3 in lymphoblastoid cell line. We further found that the 8q24 functional domains more likely interacted with genomic regions containing genes enriched in critical pathways such as Wnt signaling and promoter motifs such as E2F1 and TCF3. This result suggests that the risk locus may function as a regulatory hub by physical interactions with multiple genes important for prostate carcinogenesis. Further understanding genetic effect and biological mechanism of these chromatin interactions will shed light on the newly discovered regulatory role of the risk locus in PC etiology and progression.

Admixture mapping identifies a quantitative trait locus associated with FEV1/FVC in the COPDGene Study

African Americans are admixed with genetic contributions from European and African ancestral populations. Admixture mapping leverages this information to map genes influencing differential disease risk across populations. We performed admixture and association mapping in 3,300 African American current or former smokers from the COPDGene Study. We analyzed estimated local ancestry and SNP genotype information to identify regions associated with FEV1 /FVC, the ratio of forced expiratory volume in one second to forced vital capacity, measured by spirometry performed after bronchodilator administration. Global African ancestry inversely associated with FEV1 /FVC (P = 0.035). Genome-wide admixture analysis, controlling for age, gender, body mass index, current smoking status, pack-years smoked, and four principal components summarizing the genetic background of African Americans in the COPDGene Study, identified a region on chromosome 12q14.1 associated with FEV1 /FVC (P = 2.1 × 10(-6) ) when regressed on local ancestry. Allelic association in this region of chromosome 12 identified an intronic variant in FAM19A2 (rs348644) as associated with FEV1 /FVC (P = 1.76 × 10(-6) ). By combining admixture and association mapping, a marker on chromosome 12q14.1 was identified as being associated with reduced FEV1 /FVC ratio among African Americans in the COPDGene Study.

Prostate cancer risk alleles are associated with prostate cancer volume and prostate size

PURPOSE

Genome-wide association studies have identified an increasing number of single nucleotide polymorphisms associated with prostate cancer risk. Some of these genetic variants are also associated with serum prostate specific antigen levels and lower urinary tract symptoms, raising the question of whether they are truly prostate cancer biomarkers or simply lead to detection bias. Therefore, we determined whether single nucleotide polymorphisms associated with prostate cancer risk are more strongly associated with tumor or prostate volume.

MATERIALS AND METHODS

The genotypes of 38 validated prostate cancer risk single nucleotide polymorphisms were determined in 1,321 white men who underwent radical prostatectomy. Univariate and multivariate analyses were performed to compare the relationship of single nucleotide polymorphism frequency with total prostate and tumor volumes.

RESULTS

On multivariate analysis 2 single nucleotide polymorphisms on chromosome 8q24, rs16901979 (A) and rs6983267 (G), were significantly associated with increased tumor volume (p=0.01 and 0.02, respectively). In contrast, rs17632542 (T) near the PSA gene on 19q13 was associated with significantly lower tumor volume and rs10788160 (A) on 10q26 was associated with significantly larger prostate volume (p=0.02 and 0.01, respectively).

CONCLUSIONS

Analysis of 38 single nucleotide polymorphisms associated with prostate cancer risk revealed a significant association between several on chromosome 8q24 and increased tumor volume but not prostate volume. This suggests that they are bona fide markers of prostate cancer susceptibility and possibly more aggressive disease. Other prostate cancer risk alleles are associated with prostate specific antigen and increased prostate or decreased tumor volume, suggesting detection bias due to their phenotypic influence.

Validation of six genetic determinants of susceptibility to estrogen-induced mammary cancer in the rat and assessment of their relevance to breast cancer risk in humans

When treated with 17β-estradiol, female ACI rats (Rattus norvegicus) rapidly develop mammary cancers that share multiple phenotypes with luminal breast cancers. Seven distinct quantitative trait loci that harbor genetic determinants of susceptibility to 17β-estradiol−induced mammary cancer have been mapped in reciprocal intercrosses between susceptible ACI rats and resistant Brown Norway (BN) rats. A panel of unique congenic rat strains has now been generated and characterized to confirm the existence of these quantitative trait loci, designated Emca3 through Emca9, and to quantify their individual effects on susceptibility to 17β-estradiol−induced mammary cancer. Each congenic strain carries BN alleles spanning an individual Emca locus, introgressed onto the ACI genetic background. Data presented herein indicate that BN alleles at Emca3, Emca4, Emca5, Emca6, and Emca9 reduce susceptibility to 17β-estradiol−induced mammary cancer, whereas BN alleles at Emca7 increase susceptibility, thereby confirming the previous interval mapping data. All of these Emca loci are orthologous to regions of the human genome that have been demonstrated in genome-wide association studies to harbor genetic variants that influence breast cancer risk. Moreover, four of the Emca loci are orthologous to loci in humans that have been associated with mammographic breast density, a biomarker of breast cancer risk. This study further establishes the relevance of the ACI and derived congenic rat models of 17β-estradiol−induced mammary cancer for defining the genetic bases of breast cancer susceptibility and elucidating the mechanisms through which 17β-estradiol contributes to breast cancer development.

Long noncoding RNAs and prostate carcinogenesis: the missing 'linc'?

Long noncoding RNAs (lncRNAs) are rapidly becoming essential pieces in the cancer puzzle. Our understanding of their functional capabilities is in its infancy. One certain fact, however, is that their molecular interactions extend beyond chromatin complexes into diverse biological processes. In prostate cancer, aberrant expression of lncRNAs is associated with disease progression. Overexpression of oncogenic lncRNAs promotes tumor-cell proliferation and metastasis through chromatin looping and distal engagement with the androgen receptor, antisense gene regulation, alternative splicing, and impeding DNA repair. Several lncRNAs, such as prostate cancer antigen 3 (PCA3), prostate cancer gene expression marker 1 (PCGEM1), and prostate cancer associated ncRNA transcript 1 (PCAT1), are highly prostate-specific, posing as attractive biomarkers. Herein we review the mechanisms of action of lncRNAs in prostate carcinogenesis and their potential clinical utility for disease.

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.

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.

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.

An expressed retrogene of the master embryonic stem cell gene POU5F1 is associated with prostate cancer susceptibility

Genetic association studies of prostate and other cancers have identified a major risk locus at chromosome 8q24. Several independent risk variants at this locus alter transcriptional regulatory elements, but an affected gene and mechanism for cancer predisposition have remained elusive. The retrogene POU5F1B within the locus has a preserved open reading frame encoding a homolog of the master embryonic stem cell transcription factor Oct4. We find that 8q24 risk alleles are expression quantitative trait loci correlated with reduced expression of POU5F1B in prostate tissue and that predicted deleterious POU5F1B missense variants are also associated with risk of transformation. POU5F1 is known to be self-regulated by the encoded Oct4 transcription factor. We further observe that POU5F1 expression is directly correlated with POU5F1B expression. Our results suggest that a pathway critical to self-renewal of embryonic stem cells may also have a role in the origin of cancer.

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.

The missing story behind Genome Wide Association Studies: single nucleotide polymorphisms in gene deserts have a story to tell

Genome wide association studies are central to the evolution of personalized medicine. However, the propensity for single nucleotide polymorphisms (SNPs) to fall outside of genes means that understanding how these polymorphisms alter cellular function requires an expanded view of human genetics. Integrating the study of genome structure (chromosome conformation capture) into its function opens up new avenues of exploration. Changes in the epigenome associated with SNPs in gene deserts will allow us to define complex diseases in a much clearer manner, and usher in a new era of disease pathway exploration.