Lessons from functional analysis of genome-wide association studies

Association between novel genetic variants of Notch signaling pathway genes and survival of hepatitis B virus-related hepatocellular carcinoma

BACKGROUND

Although the Notch pathway plays an important role in formation and progression of hepatocellular carcinoma (HCC), few studies have reported the associations between functional genetic variants and the survival of hepatitis B virus (HBV)-related HCC.

METHODS

In the present study, we performed multivariable Cox proportional hazard regression analysis to evaluate associations between 36,101 SNPs in 264 Notch pathway-related genes and overall survival (OS) of 866 patients with HBV-related HCC.

RESULTS

It was found that three independent SNPs (NEURL1B rs4868192, CNTN1 rs444927 and FCER2 rs1990975) were significantly associated with the HBV-related HCC OS. The number of protective genotypes (NPGs) were significantly associated with better survival in a dose-response manner (ptrend <0.001). Compared with the model with sole clinical factors, the addition of protective genotypes to the predict models significantly increased the AUC, i.e., from 72.72% to 75.13% (p = 0.002) and from 72.04% to 74.76 (p = 0.004) for 3-year and 5-year OS, respectively. The expression quantitative trait loci (eQTL) analysis further revealed that the rs4868192 C allele was associated with lower mRNA expression levels of NEURL1B in the whole blood (p = 1.71 × 10-3), while the rs1990975 T allele was correlated with higher mRNA expression levels of FCER2 in the whole blood and normal liver tissues (p = 3.51 × 10-5 and 0.033, respectively).

CONCLUSIONS

Three potentially functional SNPs of NEURL1B, CNTN1 and FCER2 may serve as potential prognostic biomarkers for HBV-related HCC.

The PENGUIN approach to reconstruct protein interactions at enhancer-promoter regions and its application to prostate cancer

We introduce Promoter-Enhancer-Guided Interaction Networks (PENGUIN), a method for studying protein-protein interaction (PPI) networks within enhancer-promoter interactions. PENGUIN integrates H3K27ac-HiChIP data with tissue-specific PPIs to define enhancer-promoter PPI networks (EPINs). We validated PENGUIN using cancer (LNCaP) and benign (LHSAR) prostate cell lines. Our analysis detected EPIN clusters enriched with the architectural protein CTCF, a regulator of enhancer-promoter interactions. CTCF presence was coupled with the prevalence of prostate cancer (PrCa) single nucleotide polymorphisms (SNPs) within the same EPIN clusters, suggesting functional implications in PrCa. Within the EPINs displaying enrichments in both CTCF and PrCa SNPs, we also show enrichment in oncogenes. We substantiated our identified SNPs through CRISPR/Cas9 knockout and RNAi screens experiments. Here we show that PENGUIN provides insights into the intricate interplay between enhancer-promoter interactions and PPI networks, which are crucial for identifying key genes and potential intervention targets. A dedicated server is available at https://penguin.life.bsc.es/ .

Exploring the Association Between PRC2 Genes Variants and Lung Cancer Risk in Chinese Han Population

Background

Genetic susceptibilities play a large role in the pathogenesis of lung cancer (LC). The polycomb repressive complex 2 (PRC2) is a conserved chromatin-associated complex that represses gene expression and is crucial for proper organismal development and gene expression patterns. Despite PRC2 dysregulation has been observed in various human cancers, the relationship between PRC2 genes variants and lung cancer risk remains largely unexplored.

Methods

To investigate the association between single nucleotide polymorphisms (SNPs) in PRC2 genes and the risk of developing LC, we genotyped blood genomic DNA from 270 LC patients and 452 healthy individuals of Chinese Han ethnicity using the TaqMan™ genotyping technique.

Results

We found that rs17171119T>G(adjusted odds ratio (OR) = 0.662, 95% CI: 0.467-0.938, P < 0.05), rs10898459 T>C(adjusted OR = 0.615, 95% CI: 0.4-0.947, P < 0.05), and rs1136258 C>T(adjusted OR = 0.273, 95% CI: 0.186-0.401, P < 0.001) were significantly associated with a reduced risk of LC. Stratified analysis revealed a protective effect of rs17171119 in both male and female patients, specifically those with lung adenocarcinoma (LUAD). Additionally, rs1391221 showed a protective effect in both the LUAD and lung squamous cell carcinoma (LUSC) groups, while rs1136258 exhibited a protective effect in both females and males, as well as in both LUAD and LUSC groups. Furthermore, analysis of The Cancer Genome Atlas (TCGA) dataset revealed expression levels of EED and RBBP4 in both LUAD and LUSC.

Conclusion

This study provides evidence that allelic variants in EZH2, EED, and RBBP4 may act as protective factors against LC development and could serve as genetic markers associated with susceptibility to LC.

Genetic variants of SOS2, MAP2K1 and RASGRF2 in the RAS pathway genes predict survival of HBV-related hepatocellular carcinoma patients

The RAS pathway participates in the cascade of proliferation and cell division process, and the activated RAS pathway can lead to tumorigenesis including hepatocellular carcinoma (HCC). However, few studies have explored the effects of genetic variants in the RAS pathway-related genes on the survival of patients with HBV-related HCC. In the present study, we assessed the associations between 11,658 single-nucleotide polymorphisms (SNPs) in 62 RAS pathway genes and the overall survival (OS) of 866 HBV-related HCC individuals, which were randomly split (1:1) into discovery and validation datasets. As a result, three potentially functional SNPs were identified, based on multivariable cox proportional hazards regression analyses, in SOS Ras/Rho guanine nucleotide exchange factor 2 (SOS2, rs4632055 A > G), Ras protein-specific guanine nucleotide releasing factor 2 (RASGRF2, rs26418A > G) and mitogen-activated protein kinase 1 (MAP2K1,rs57120695 C > T), which were significantly and independently associated with OS of HBV-related HCC patients [adjusted hazards ratios (HRs) of 1.42, 1.32 and 1.50, respectively; 95% confidence intervals (CI), 1.14 to 1.76, 1.15 to 1.53 and 1.15 to 1.97, respectively; P = 0.001, < 0.001 and 0.003, respectively]. Additionally, the joint effects as the unfavorable genotypes of these three SNPs showed a significant association with the poor survival of HCC (trend test P < 0.001). The expression quantitative trait loci (eQTL) analysis further revealed that the rs4632055 G allele and the rs26418 A allele were associated with lower mRNA expression levels of SOS2 and RASGRF2, respectively. Collectively, these potentially functional SNPs of RASGRF2, SOS2 and M2PAK1 may become potential prognostic biomarkers for HBV-related HCC after hepatectomy.

Regulation of the Late Onset alzheimer's Disease Associated HLA-DQA1/DRB1 Expression

(Genome-wide Association Studies) GWAS have identified ∼42 late-onset Alzheimer's disease (LOAD)-associated loci, each of which contains multiple single nucleotide polymorphisms (SNPs) in linkage disequilibrium (LD) and most of these SNPs are in the non-coding region of human genome. However, how these SNPs regulate risk gene expression remains unknown. In this work, by using a set of novel techniques, we identified 6 functional SNPs (fSNPs) rs9271198, rs9271200, rs9281945, rs9271243, and rs9271247 on the LOAD-associated HLA-DRB1/DQA1 locus and 42 proteins specifically binding to five of these 6 fSNPs. As a proof of evidence, we verified the allele-specific binding of GATA2 and GATA3, ELAVL1 and HNRNPA0, ILF2 and ILF3, NFIB and NFIC, as well as CUX1 to these five fSNPs, respectively. Moreover, we demonstrate that all these nine proteins regulate the expression of both HLA-DQA1 and HLA-DRB1 in human microglial cells. The contribution of HLA class II to the susceptibility of LOAD is discussed.

Post-GWAS functional studies reveal an RA-associated CD40-induced NF-kB signal transduction and transcriptional regulation network targeted by class II HDAC inhibitors

Currently, it remains difficult to identify which single nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWAS) are functional and how various functional SNPs (fSNPs) interact and contribute to disease susceptibility. GWAS have identified a CD40 locus that is associated with rheumatoid arthritis (RA). We previously used two techniques developed in our laboratory, single nucleotide polymorphism-next-generation sequencing (SNP-seq) and flanking restriction enhanced DNA pulldown-mass spectrometry (FREP-MS), to determine that the RA risk gene RBPJ regulates CD40 expression via a fSNP at the RA-associated CD40 locus. In the present work, by applying the same approach, we report the identification of six proteins that regulate RBPJ expression via binding to two fSNPs on the RA-associated RBPJ locus. Using these findings, together with the published data, we constructed an RA-associated signal transduction and transcriptional regulation network (STTRN) that functionally connects multiple RA-associated risk genes via transcriptional regulation networks (TRNs) linked by CD40-induced nuclear factor kappa B (NF-kB) signaling. Remarkably, this STTRN provides insight into the potential mechanism of action for the histone deacetylase inhibitor givinostat, an approved therapy for systemic juvenile idiopathic arthritis. Thus, the generation of disease-associated STTRNs based on post-GWAS functional studies is demonstrated as a novel and effective approach to apply GWAS for mechanistic studies and target identification.

Gene regulation by long non-coding RNAs and its biological functions

Evidence accumulated over the past decade shows that long non-coding RNAs (lncRNAs) are widely expressed and have key roles in gene regulation. Recent studies have begun to unravel how the biogenesis of lncRNAs is distinct from that of mRNAs and is linked with their specific subcellular localizations and functions. Depending on their localization and their specific interactions with DNA, RNA and proteins, lncRNAs can modulate chromatin function, regulate the assembly and function of membraneless nuclear bodies, alter the stability and translation of cytoplasmic mRNAs and interfere with signalling pathways. Many of these functions ultimately affect gene expression in diverse biological and physiopathological contexts, such as in neuronal disorders, immune responses and cancer. Tissue-specific and condition-specific expression patterns suggest that lncRNAs are potential biomarkers and provide a rationale to target them clinically. In this Review, we discuss the mechanisms of lncRNA biogenesis, localization and functions in transcriptional, post-transcriptional and other modes of gene regulation, and their potential therapeutic applications.

Mapping of Diabetes Susceptibility Loci in a Domestic Cat Breed with an Unusually High Incidence of Diabetes Mellitus

Genetic variants that are associated with susceptibility to type 2 diabetes (T2D) are important for identification of individuals at risk and can provide insights into the molecular basis of disease. Analysis of T2D in domestic animals provides both the opportunity to improve veterinary management and breeding programs as well as to identify novel T2D risk genes. Australian-bred Burmese (ABB) cats have a 4-fold increased incidence of type 2 diabetes (T2D) compared to Burmese cats bred in the United States. This is likely attributable to a genetic founder effect. We investigated this by performing a genome-wide association scan on ABB cats. Four SNPs were associated with the ABB T2D phenotype with p values <0.005. All exons and splice junctions of candidate genes near significant single-nucleotide polymorphisms (SNPs) were sequenced, including the genes DGKG, IFG2BP2, SLC8A1, E2F6, ETV5, TRA2B and LIPH. Six candidate polymorphisms were followed up in a larger cohort of ABB cats with or without T2D and also in Burmese cats bred in America, which exhibit low T2D incidence. The original SNPs were confirmed in this cohort as associated with the T2D phenotype, although no novel coding SNPs in any of the seven candidate genes showed association with T2D. The identification of genetic markers associated with T2D susceptibility in ABB cats will enable preventative health strategies and guide breeding programs to reduce the prevalence of T2D in these cats.

A sequential methodology for the rapid identification and characterization of breast cancer-associated functional SNPs

GWAS cannot identify functional SNPs (fSNP) from disease-associated SNPs in linkage disequilibrium (LD). Here, we report developing three sequential methodologies including Reel-seq (Regulatory element-sequencing) to identify fSNPs in a high-throughput fashion, SDCP-MS (SNP-specific DNA competition pulldown-mass spectrometry) to identify fSNP-bound proteins and AIDP-Wb (allele-imbalanced DNA pulldown-Western blot) to detect allele-specific protein:fSNP binding. We first apply Reel-seq to screen a library containing 4316 breast cancer-associated SNPs and identify 521 candidate fSNPs. As proof of principle, we verify candidate fSNPs on three well-characterized loci: FGFR2, MAP3K1 and BABAM1. Next, using SDCP-MS and AIDP-Wb, we rapidly identify multiple regulatory factors that specifically bind in an allele-imbalanced manner to the fSNPs on the FGFR2 locus. We finally demonstrate that the factors identified by SDCP-MS can regulate risk gene expression. These data suggest that the sequential application of Reel-seq, SDCP-MS, and AIDP-Wb can greatly help to translate large sets of GWAS data into biologically relevant information.

The Oncogenic Potential of the Centromeric Border Protein FAM84B of the 8q24.21 Gene Desert

FAM84B is a risk gene in breast and prostate cancers. Its upregulation is associated with poor prognosis of prostate cancer, breast cancer, and esophageal squamous cell carcinoma. FAM84B facilitates cancer cell proliferation and invasion in vitro, and xenograft growth in vivo. The FAM84B and Myc genes border a 1.2 Mb gene desert at 8q24.21. Co-amplification of both occurs in 20 cancer types. Mice deficient of a 430 Kb fragment within the 1.2 Mb gene desert have downregulated FAM84B and Myc expressions concurrent with reduced breast cancer growth. Intriguingly, Myc works in partnership with other oncogenes, including Ras. FAM84B shares similarities with the H-Ras-like suppressor (HRASLS) family over their typical LRAT (lecithin:retinal acyltransferase) domain. This domain contains a catalytic triad, H23, H35, and C113, which constitutes the phospholipase A_1/2 and O-acyltransferase activities of HRASLS1-5. These enzymatic activities underlie their suppression of Ras. FAM84B conserves H23 and H35 but not C113 with both histidine residues residing within a highly conserved motif that FAM84B shares with HRASLS1-5. Deletion of this motif abolishes FAM84B oncogenic activities. These properties suggest a collaboration of FAM84B with Myc, consistent with the role of the gene desert in strengthening Myc functions. Here, we will discuss recent research on FAM84B-derived oncogenic potential.

Long non-coding RNA polymorphisms on 8q24 are associated with the prognosis of gastric cancer in a Chinese population

Background

Gastric cancer (GC) remains the third leading cause of cancer death in China. Although genome-wide association studies have identified the association between several single nucleotide polymorphisms (SNPs) on 8q24 and the risk of GC, the role of these SNPs in the prognosis of GC in Chinese populations has not yet been fully evaluated. Therefore, this study was conducted to explore the association between long non-coding RNA (lncRNA) polymorphisms on 8q24 and the prognosis of GC.

Methods

We genotyped 726 surgically resected GC patients to explore the association between eight SNPs in the lncRNAs CCAT1 (rs10087719, rs7816475), PCAT1 (rs1026411), PRNCR1 (rs12682421, rs13252298), and CASC8 (rs1562430, rs4871789, rs6983267) transcribed from the 8q24 locus and the prognosis of GC in a Chinese population.

Results

We found that the patients carrying rs12682421 AA genotypes survived for a shorter time than those with the GG/GA genotype (HR = 1.39, 95% confidence interval (CI) [1.09-1.78]). Compared with the CC/CT genotype, the TT genotype of rs1562430 was associated with an increased risk of death (HR = 1.38, 95% CI [1.06-1.80]). Furthermore, the results also identified the rs1026411 SNP as an independent prognostic factor for poor survival in GC patients. Patients carrying AA/AG variant genotypes had a 36% increased risk of death compared to those carrying the GG genotype (HR = 1.36, 95% CI [1.06-1.74]). These findings suggested that the rs12682421, rs1026411 and rs1562430 SNPs may contribute to the survival of GC and be prognostic markers for GC.

Integrative comparison of the genomic and transcriptomic landscape between prostate cancer patients of predominantly African or European genetic ancestry

Men of predominantly African Ancestry (AA) have higher prostate cancer (CaP) incidence and worse survival than men of predominantly European Ancestry (EA). While socioeconomic factors drive this disparity, genomic factors may also contribute to differences in the incidence and mortality rates. To compare the prevalence of prostate tumor genomic alterations and transcriptomic profiles by patient genetic ancestry, we evaluated genomic profiles from The Cancer Genome Atlas (TCGA) CaP cohort (n = 498). Patient global and local genetic ancestry were estimated by computational algorithms using genotyping data; 414 (83.1%) were EA, 61 (12.2%) were AA, 11 (2.2%) were East Asian Ancestry (EAA), 10 (2.0%) were Native American (NA), and 2 (0.4%) were other ancestry. Genetic ancestry was highly concordant with self-identified race/ethnicity. Subsequent analyses were limited to 61 AA and 414 EA cases. Significant differences were observed by ancestry in the frequency of SPOP mutations (20.3% AA vs. 10.0% EA; p = 5.6×10-03), TMPRSS2-ERG fusions (29.3% AA vs. 39.6% EA; p = 4.4×10-02), and PTEN deletions/losses (11.5% AA vs. 30.2% EA; p = 3.5×10-03). Differentially expressed genes (DEGs) between AAs and EAs showed significant enrichment for prostate eQTL target genes (p = 8.09×10-48). Enrichment of highly expressed DEGs for immune-related pathways was observed in AAs, and for PTEN/PI3K signaling in EAs. Nearly one-third of DEGs (31.3%) were long non-coding RNAs (DE-lncRNAs). The proportion of DE-lncRNAs with higher expression in AAs greatly exceeded that with lower expression in AAs (p = 1.2×10-125). Both ChIP-seq and RNA-seq data suggested a stronger regulatory role for AR signaling pathways in DE-lncRNAs vs. non-DE-lncRNAs. CaP-related oncogenic lncRNAs, such as PVT1, PCAT1 and PCAT10/CTBP1-AS, were found to be more highly expressed in AAs. We report substantial heterogeneity in the prostate tumor genome and transcriptome between EA and AA. These differences may be biological contributors to racial disparities in CaP incidence and outcomes.

Fruit and vegetable intake modifies the associations between suppressor of cytokine signaling 3 genetic variants and type 2 diabetes

PURPOSE

Type 2 diabetes is a complex disease determined by variable genes and environmental factors. The study was designed to investigate the effect of interactions of four polymorphisms of suppressor of cytokine signaling 3 (SOCS3) with fruit and vegetable (F&V) intake on type 2 diabetes in a rural population of China.

METHODS

A total of 4411 participants from the rural areas of Henan, China were included in the study. Multivariate logistic regression and restricted cubic splines were used to estimate the associations between polymorphisms and risk allele score of SOCS3 and type 2 diabetes in different groups. Haplotype analysis was conducted to examine the effects of linkage inheritance at these four loci on type 2 diabetes.

RESULTS

Three of the four polymorphisms showed significant associations with type 2 diabetes in the less F&V intake group after adjusting the covariates, the odds ratios (ORs) and corresponding 95% confidence intervals (95% CIs) were 1.24 (1.08-1.41) for rs4969168, 1.16 (1.02-1.32) for rs9892622, and 1.21 (1.06-1.39) for rs9914220. No significant association was detected in the more F&V intake group. The obvious dose-response relationship between the risk allele score and type 2 diabetes was also noted only in the less F&V intake group.

CONCLUSIONS

Variants of SOCS3 gene were associated with type 2 diabetes and the associations could be modified by the F&V intake.

Variants in oxidative stress-related genes affect the chemosensitivity through Nrf2-mediated signaling pathway in biliary tract cancer

Background

Oxidative stress and their effectors play critical roles in carcinogenesis and chemoresistance. However, the role of oxidative stress-related genes variants in biliary tract cancer (BTC) chemoresistance remains unknown. In this work, we aim to investigate oxidative stress-dependent molecular mechanisms underlying chemoresistance, and find potential biomarkers to predict chemotherapy response for BTC.

Methods

Sixty-six SNPs in 21 oxidative stress-related genes were genotyped and analyzed in 367 BTC patients. Immunoblot, immunohistochemical, immunofluorescent, quantitative PCR, chromatin immunoprecipitation analysis and study of animal xenograft models were performed to discover oxidative stress-related susceptibility genes underlying chemoresistance mechanism of BTC.

Findings

We found that 3 functional polymorphisms (CAT_rs769217, GPX4_rs4807542, and GSR_rs3779647), which were shown to affect their respective gene expression levels, modified the effect of chemotherapy on overall survival (OS). We then demonstrated that knockdown of GPX4, CAT, or GSR induced chemoresistance through elevation of ROS level and activation of Nrf2-ABCG2 pathway in BTC cell lines. Moreover, the association between Nrf2 expression and BTC prognosis is only found in patients who received chemotherapy. Knockdown of Nrf2 enhanced chemosensitivity or even eliminated postoperative recurrence in BTC xenograft mouse models. Importantly, upon chemotherapy treatment patients harboring high oxidative stress-related score received higher survival benefit from adjuvant chemotherapy compared with patients with low oxidative stress-related score.

Interpretation

The result of our study suggests, for the first time, that the oxidative stress-related score calculated by combining variations in CAT, GPX4, and GSR or Nrf2 expression could be used for predicting the chemosensitivity of BTC patients.

Fund

This work was supported by the National Science Foundation of China, Foundation of Shanghai Shen Kang Hospital Development Center, and Shanghai Outstanding Academic Leaders Plan.

HOTAIR lncRNA SNPs rs920778 and rs1899663 are associated with smoking, male gender, and squamous cell carcinoma in a Chinese lung cancer population

The abnormal expression of the long noncoding RNA (lncRNA) HOX transcript intergenic antisense RNA (HOTAIR) plays an important role in the development of various cancers; however, single nucleotide polymorphisms (SNPs) in HOTAIR and their association with primary lung cancer susceptibility have not yet been reported. Here, we performed a case-control study including 262 primary lung cancer patients and 451 cancer-free control individuals to investigate the association between four haplotype-tagging SNPs (rs920778, rs12826786, rs4759314, and rs1899663) in the HOTAIR lncRNA and the risk of developing primary lung cancer. We found a significant association between the SNPs rs920778 and rs1899663 in the HOTAIR and primary lung cancer susceptibility (P < 0.05). Moreover, homozygous C/T (C/T + TT) for rs920778 (C > T) sites was significantly associated with gender, smoking history, and pathological type. In addition, linkage disequilibrium and haplotype analysis of HOTAIR gene polymorphisms for susceptibility to lung cancer revealed a high degree of linkage disequilibrium between the rs920778 and rs1899663 loci (D' = 0.86, r2 = 0.52). The population of rs920778, rs1899663, and rs4759314 had a significantly increased risk of lung cancer (P < 0.001). In summary, the present study provides persuasive evidence that SNP rs920778 is closely correlated with susceptibility to primary lung cancer. Future studies are warranted to validate and expand these findings, and to further dissect the importance of these SNPs in the development of primary lung cancer.

A Genome-Wide Association Study Reveals Candidate Genes Related to Salt Tolerance in Rice (Oryza sativa) at the Germination Stage

Salt toxicity is the major factor limiting crop productivity in saline soils. In this paper, 295 accessions including a heuristic core set (137 accessions) and 158 bred varieties were re-sequenced and ~1.65 million SNPs/indels were used to perform a genome-wide association study (GWAS) of salt-tolerance-related phenotypes in rice during the germination stage. A total of 12 associated peaks distributed on seven chromosomes using a compressed mixed linear model were detected. Determined by linkage disequilibrium (LD) blocks analysis, we finally obtained a total of 79 candidate genes. By detecting the highly associated variations located inside the genic region that overlapped with the results of LD block analysis, we characterized 17 genes that may contribute to salt tolerance during the seed germination stage. At the same time, we conducted a haplotype analysis of the genes with functional variations together with phenotypic correlation and orthologous sequence analyses. Among these genes, OsMADS31, which is a MADS-box family transcription factor, had a down-regulated expression under the salt condition and it was predicted to be involved in the salt tolerance at the rice germination stage. Our study revealed some novel candidate genes and their substantial natural variations in the rice genome at the germination stage. The GWAS in rice at the germination stage would provide important resources for molecular breeding and functional analysis of the salt tolerance during rice germination.

Common Genetic Variation and Susceptibility to Ovarian Cancer: Current Insights and Future Directions

In this review, we summarize current progress in the genetic epidemiology of epithelial ovarian cancer (EOC), focusing exclusively on elucidating the role of common germline genetic variation in conferring susceptibility to EOC. We provide an overview of the more than 30 EOC risk loci identified to date by genome-wide association studies (GWAS) and describe the contribution of large-scale, cross-cancer type, custom genotyping projects, such as the OncoArray and the Collaborative Oncological Gene-Environment Study, to locus discovery and replication. We discuss the histotype-specific nature of these EOC risk loci, pleiotropy, or overlapping genetic effects between EOC and other hormone-related cancer types, and the application of findings to polygenic risk prediction for EOC. The second part of the article offers a concise review of primarily laboratory-based studies that have led to the identification of several putative EOC susceptibility genes using common variants at the known EOC risk loci as starting points. More global biological insights emerging from network- and pathway-based analyses of GWAS for EOC susceptibility are also highlighted. Finally, we delve into potential future directions, including the need to identify EOC risk loci in non-European populations and the next generation of GWAS functional studies that are likely to involve genome editing to establish the cell type-specific carcinogenic effects of EOC risk variants Cancer Epidemiol Biomarkers Prev; 27(4); 395-404. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in Cancer."

Rs2853677 modulates Snail1 binding to the TERT enhancer and affects lung adenocarcinoma susceptibility

Genome wide association studies (GWAS) have shown that SNPs in non-coding regions are associated with inherited susceptibility to cancer. The effect of one single SNP, however, is weak. To identify potential co-factors of SNPs, we investigated the underlying mechanism by which SNPs affect lung cancer susceptibility. We found that rs2853677 is located within the Snail1 binding site in a TERT enhancer. This enhancer increases TERT transcription when juxtaposed to the TERT promoter. The binding of Snail1 to the enhancer disrupts enhancer-promoter colocalization and silences TERT transcription. The high risk variant of rs2853677 disrupts the Snail1 binding site and derepresses TERT expression in response to Snail1 upregulation, thus increasing lung adenocarcinoma susceptibility. Our data suggest that Snail1 may be a co-factor of rs2853677 for predicting lung adenocarcinoma susceptibility and prognosis.

Inherited pancreatic cancer

Pancreatic cancers arise through a series of genetic events both inherited and acquired. Inherited genetic changes, both high penetrance and low penetrance, are an important component of pancreatic cancer risk, and may be used to characterize populations who will benefit from early detection. Furthermore, pancreatic cancer patients with inherited mutations may be particularly sensitive to certain targeted agents, providing an opportunity to personalized treatment. Family history of pancreatic cancer is one of the strongest risk factors for the disease, and is associated with an increased risk of caners at other sites, including but not limited to breast, ovarian and colorectal cancer. The goal of this chapter is to discuss the importance of family history of pancreatic cancer, and the known genes that account for a portion of the familial clustering of pancreatic cancer.

Genome-wide imaging association study implicates functional activity and glial homeostasis of the caudate in smoking addiction

Background

Nearly 6 million deaths and over a half trillion dollars in healthcare costs worldwide are attributed to tobacco smoking each year. Extensive research efforts have been pursued to elucidate the molecular underpinnings of smoking addiction and facilitate cessation. In this study, we genotyped and obtained both resting state and task-based functional magnetic resonance imaging from 64 non-smokers and 42 smokers. Smokers were imaged after having smoked normally (“sated”) and after having not smoked for at least 12 h (“abstinent”).

Results

While abstinent smokers did not differ from non-smokers with respect to pairwise resting state functional connectivities (RSFCs) between 12 brain regions of interest, RSFCs involving the caudate and putamen of sated smokers significantly differed from those of non-smokers (P < 0.01). Further analyses of caudate and putamen activity during elicited experiences of reward and disappointment show that caudate activity during reward (CR) correlated with smoking status (P = 0.015). Moreover, abstinent smokers with lower CR experienced greater withdrawal symptoms (P = 0.024), which suggests CR may be related to smoking urges. Associations between genetic variants and CR, adjusted for smoking status, were identified by genome-wide association study (GWAS). Genes containing or exhibiting caudate-specific expression regulation by these variants were enriched within Gene Ontology terms that describe cytoskeleton functions, synaptic organization, and injury response (P < 0.001, FDR < 0.05).

Conclusions

By integrating genomic and imaging data, novel insights into potential mechanisms of caudate activation and homeostasis are revealed that may guide new directions of research toward improving our understanding of addiction pathology.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4124-5) contains supplementary material, which is available to authorized users.

Altered TERT promoter and other genomic regulatory elements: occurrence and impact

Study of genetic alterations, inherited or acquired, that increase the risk or drive cancers and many other diseases had remained mostly confined to coding sequences of the human genome. Data from genome wide associations studies, development of the Encyclopedia of DNA Elements (ENCODE), and a spurt in detection of driver somatic mutations have shifted focus towards noncoding regions of the human genome. The majority of genetic variants robustly associated with cancers and other syndromes identified through genome wide studies are located within noncoding regulatory regions of the genome. Genome wide techniques have put an emphasis on the role of three-dimensional chromosomal structures and cis-acting elements in regulations of different genes. The variants within noncoding genomic regions can potentially alter a number of regulatory elements including promoters, enhancers, insulators, noncoding long RNAs and others that affect cancers and various diseases through altered expression of critical genes. With effect of genetic alterations within regulatory elements dependent on other partner molecules like transcription factors and histone marks, an understanding of such modifications can potentially identify extended therapeutic targets. That concept has been augmented by the detection of driver somatic noncoding mutations within the promoter region of the telomerase reverse transcriptase (TERT) gene in different cancers. The acquired somatic noncoding mutations within different regulatory elements are now being reported in different cancers with an increased regularity. In this review we discuss the occurrence and impact of germline and somatic alterations within the TERT promoter and other genomic regulatory elements.

Enrichment of putative PAX8 target genes at serous epithelial ovarian cancer susceptibility loci

BACKGROUND

Genome-wide association studies (GWAS) have identified 18 loci associated with serous ovarian cancer (SOC) susceptibility but the biological mechanisms driving these findings remain poorly characterised. Germline cancer risk loci may be enriched for target genes of transcription factors (TFs) critical to somatic tumorigenesis.

METHODS

All 615 TF-target sets from the Molecular Signatures Database were evaluated using gene set enrichment analysis (GSEA) and three GWAS for SOC risk: discovery (2196 cases/4396 controls), replication (7035 cases/21 693 controls; independent from discovery), and combined (9627 cases/30 845 controls; including additional individuals).

RESULTS

The PAX8-target gene set was ranked 1/615 in the discovery (PGSEA<0.001; FDR=0.21), 7/615 in the replication (PGSEA=0.004; FDR=0.37), and 1/615 in the combined (PGSEA<0.001; FDR=0.21) studies. Adding other genes reported to interact with PAX8 in the literature to the PAX8-target set and applying an alternative to GSEA, interval enrichment, further confirmed this association (P=0.006). Fifteen of the 157 genes from this expanded PAX8 pathway were near eight loci associated with SOC risk at P<10-5 (including six with P<5 × 10-8). The pathway was also associated with differential gene expression after shRNA-mediated silencing of PAX8 in HeyA8 (PGSEA=0.025) and IGROV1 (PGSEA=0.004) SOC cells and several PAX8 targets near SOC risk loci demonstrated in vitro transcriptomic perturbation.

CONCLUSIONS

Putative PAX8 target genes are enriched for common SOC risk variants. This finding from our agnostic evaluation is of particular interest given that PAX8 is well-established as a specific marker for the cell of origin of SOC.

Computational discovery of transcription factors associated with drug response

This study integrates gene expression, genotype and drug response data in lymphoblastoid cell lines with transcription factor (TF)-binding sites from ENCODE (Encyclopedia of Genomic Elements) in a novel methodology that elucidates regulatory contexts associated with cytotoxicity. The method, GENMi (Gene Expression iN the Middle), postulates that single-nucleotide polymorphisms within TF-binding sites putatively modulate its regulatory activity, and the resulting variation in gene expression leads to variation in drug response. Analysis of 161 TFs and 24 treatments revealed 334 significantly associated TF-treatment pairs. Investigation of 20 selected pairs yielded literature support for 13 of these associations, often from studies where perturbation of the TF expression changes drug response. Experimental validation of significant GENMi associations in taxanes and anthracyclines across two triple-negative breast cancer cell lines corroborates our findings. The method is shown to be more sensitive than an alternative, genome-wide association study-based approach that does not use gene expression. These results demonstrate the utility of GENMi in identifying TFs that influence drug response and provide a number of candidates for further testing.

Genomic Alterations of Non-Coding Regions Underlie Human Cancer: Lessons from T-ALL

It has been appreciated for decades that somatic genomic alterations that change coding sequences of proto-oncogenes, translocate enhancers/promoters near proto-oncogenes, or create fusion oncogenes can drive cancer by inducing oncogenic activities. An explosion of genome-wide technologies over the past decade has fueled discoveries of the roles of three-dimensional chromosome structure and powerful cis-acting elements (super-enhancers) in regulating gene transcription. In recent years, studies of human T cell acute lymphoblastic leukemia (T-ALL) using genome-wide technologies have provided paradigms for how non-coding genomic region alterations can disrupt 3D chromosome architecture or establish super-enhancers to activate oncogenic transcription of proto-oncogenes. These studies raise important issues to consider with the objective of leveraging basic knowledge into new diagnostic and therapeutic opportunities for cancer patients.

Three new pancreatic cancer susceptibility signals identified on chromosomes 1q32.1, 5p15.33 and 8q24.21

Genome-wide association studies (GWAS) have identified common pancreatic cancer susceptibility variants at 13 chromosomal loci in individuals of European descent. To identify new susceptibility variants, we performed imputation based on 1000 Genomes (1000G) Project data and association analysis using 5,107 case and 8,845 control subjects from 27 cohort and case-control studies that participated in the PanScan I-III GWAS. This analysis, in combination with a two-staged replication in an additional 6,076 case and 7,555 control subjects from the PANcreatic Disease ReseArch (PANDoRA) and Pancreatic Cancer Case-Control (PanC4) Consortia uncovered 3 new pancreatic cancer risk signals marked by single nucleotide polymorphisms (SNPs) rs2816938 at chromosome 1q32.1 (per allele odds ratio (OR) = 1.20, P = 4.88x10 -15), rs10094872 at 8q24.21 (OR = 1.15, P = 3.22x10 -9) and rs35226131 at 5p15.33 (OR = 0.71, P = 1.70x10 -8). These SNPs represent independent risk variants at previously identified pancreatic cancer risk loci on chr1q32.1 ( NR5A2), chr8q24.21 ( MYC) and chr5p15.33 ( CLPTM1L- TERT) as per analyses conditioned on previously reported susceptibility variants. We assessed expression of candidate genes at the three risk loci in histologically normal ( n = 10) and tumor ( n = 8) derived pancreatic tissue samples and observed a marked reduction of NR5A2 expression (chr1q32.1) in the tumors (fold change -7.6, P = 5.7x10 -8). This finding was validated in a second set of paired ( n = 20) histologically normal and tumor derived pancreatic tissue samples (average fold change for three NR5A2 isoforms -31.3 to -95.7, P = 7.5x10 -4-2.0x10 -3). Our study has identified new susceptibility variants independently conferring pancreatic cancer risk that merit functional follow-up to identify target genes and explain the underlying biology.

The search for cis-regulatory driver mutations in cancer genomes

With the advent of high-throughput and relatively inexpensive whole-genome sequencing technology, the focus of cancer research has begun to shift toward analyses of somatic mutations in non-coding cis-regulatory elements of the cancer genome. Cis-regulatory elements play an important role in gene regulation, with mutations in these elements potentially resulting in changes to the expression of linked genes. The recent discoveries of recurrent TERT promoter mutations in melanoma, and recurrent mutations that create a super-enhancer regulating TAL1 expression in T-cell acute lymphoblastic leukaemia (T-ALL), have sparked significant interest in the search for other somatic cis-regulatory mutations driving cancer development. In this review, we look more closely at the TERT promoter and TAL1 enhancer alterations and use these examples to ask whether other cis-regulatory mutations may play a role in cancer susceptibility. In doing so, we make observations from the data emerging from recent research in this field, and describe the experimental and analytical approaches which could be adopted in the hope of better uncovering the true functional significance of somatic cis-regulatory mutations in cancer.

The role of enhancers in cancer

Enhancer elements function as the logic gates of the genetic regulatory circuitry. One of their most important functions is the integration of extracellular signals with intracellular cell fate information to generate cell type-specific transcriptional responses. Mutations occurring in cancer often misregulate enhancers that normally control the signal-dependent expression of growth-related genes. This misregulation can result from trans-acting mechanisms, such as activation of the transcription factors or epigenetic regulators that control enhancer activity, or can be caused in cis by direct mutations that alter the activity of the enhancer or its target gene specificity. These processes can generate tumour type-specific super-enhancers and establish a 'locked' gene regulatory state that drives the uncontrolled proliferation of cancer cells. Here, we review the role of enhancers in cancer, and their potential as therapeutic targets.

Mismatch discrimination of lipidated DNA and LNA-probes (LiNAs) in hybridization-controlled liposome assembly

Assays for mismatch discrimination and detection of single nucleotide variations by hybridization-controlled assembly of liposomes, which do not require tedious surface chemistry, are versatile for both DNA and RNA targets. We report herein a comprehensive study on different DNA and LNA (locked nucleic acids) probe designs, including membrane-anchoring requirements, studies on different probes and target lengths (including overhangs), DNA and RNA targets (including sequences associated with pathogens) for lipidated nucleic acids (LiNAs). Advantages and limitations of the liposome assembly based assay in the context of mismatch discrimination and SNP detection are presented. The advantages of membrane-anchored LiNA-probes compared to chemically attached probes on solid nanoparticles (e.g. gold nanoparticles) are described. Key functionalities such as non-covalent attachment of LiNA probes without the need for long spacers and the inherent mobility of membrane-anchored probes in lipid-bilayer membranes will be described for several different probe designs.

Shared Risk Factors in Cardiovascular Disease and Cancer

Cardiovascular disease (CVD) and cancer are the 2 leading causes of death worldwide. Although commonly thought of as 2 separate disease entities, CVD and cancer possess various similarities and possible interactions, including a number of similar risk factors (eg, obesity, diabetes mellitus), suggesting a shared biology for which there is emerging evidence. Although chronic inflammation is an indispensable feature of the pathogenesis and progression of both CVD and cancer, additional mechanisms can be found at their intersection. Therapeutic advances, despite improving longevity, have increased the overlap between these diseases, with millions of cancer survivors now at risk of developing CVD. Cardiac risk factors have a major impact on subsequent treatment-related cardiotoxicity. In this review, we explore the risk factors common to both CVD and cancer, highlighting the major epidemiological studies and potential biological mechanisms that account for them.

The emerging role of lncRNAs in cancer

It is increasingly evident that many of the genomic mutations in cancer reside inside regions that do not encode proteins. However, these regions are often transcribed into long noncoding RNAs (lncRNAs). The recent application of next-generation sequencing to a growing number of cancer transcriptomes has indeed revealed thousands of lncRNAs whose aberrant expression is associated with different cancer types. Among the few that have been functionally characterized, several have been linked to malignant transformation. Notably, these lncRNAs have key roles in gene regulation and thus affect various aspects of cellular homeostasis, including proliferation, survival, migration or genomic stability. This review aims to summarize current knowledge of lncRNAs from the cancer perspective. It discusses the strategies that led to the identification of cancer-related lncRNAs and the methodologies and challenges involving the study of these molecules, as well as the imminent applications of these findings to the clinic.

A functional variant in TP63 at 3q28 associated with bladder cancer risk by creating an miR-140-5p binding site

The first genome-wide association study (GWAS) for bladder cancer has identified a susceptibility locus at 3q28 in the European ancestry. However, the causal variant at 3q28 remains unknown. We conducted a three-stage fine mapping study to identify potential causal variants in the region. A total of 41 single nucleotide polymorphisms (SNPs) across 120 kb at 3q28 were tested for association with bladder cancer risk among 3,094 bladder cancer cases and 3,738 controls. Resequencing and functional assays were further evaluated. The SNP rs35592567 in the 3'-UTR of TP63 was consistently associated with bladder cancer risk in all three stages. In the combined analysis, the T allele of rs35592567 was significantly associated with a decreased risk for bladder cancer (OR = 0.82, 95% CI = 0.75-0.90, P = 9.797 × 10(-6) in the additive model). Biochemical assays revealed that the T allele reduced the post-transcriptional levels of TP63 mediated by interfering binding efficiency of miR-140-5p. In addition, overexpression of miR-140-5p inhibited bladder cancer cell proliferation and attenuated cell migration, invasion and G1 cell-cycle arrest. Together, these results suggest that rs35592567 in TP63 may be a novel causal variant contributing to the susceptibility to bladder cancer, which provides additional insight into the pathogenesis of bladder carcinogenesis.

Association Analysis in Rice: From Application to Utilization

Association analysis based on linkage disequilibrium (LD) is an efficient way to dissect complex traits and to identify gene functions in rice. Although association analysis is an effective way to construct fine maps for quantitative traits, there are a few issues which need to be addressed. In this review, we will first summarize type, structure, and LD level of populations used for association analysis of rice, and then discuss the genotyping methods and statistical approaches used for association analysis in rice. Moreover, we will review current shortcomings and benefits of association analysis as well as specific types of future research to overcome these shortcomings. Furthermore, we will analyze the reasons for the underutilization of the results within association analysis in rice breeding.

Genetic predisposition to neuroblastoma mediated by a LMO1 super-enhancer polymorphism

Neuroblastoma is a paediatric malignancy that typically arises in early childhood, and is derived from the developing sympathetic nervous system. Clinical phenotypes range from localized tumours with excellent outcomes to widely metastatic disease in which long-term survival is approximately 40% despite intensive therapy. A previous genome-wide association study identified common polymorphisms at the LMO1 gene locus that are highly associated with neuroblastoma susceptibility and oncogenic addiction to LMO1 in the tumour cells. Here we investigate the causal DNA variant at this locus and the mechanism by which it leads to neuroblastoma tumorigenesis. We first imputed all possible genotypes across the LMO1 locus and then mapped highly associated single nucleotide polymorphism (SNPs) to areas of chromatin accessibility, evolutionary conservation and transcription factor binding sites. We show that SNP rs2168101 G>T is the most highly associated variant (combined P = 7.47 × 10(-29), odds ratio 0.65, 95% confidence interval 0.60-0.70), and resides in a super-enhancer defined by extensive acetylation of histone H3 lysine 27 within the first intron of LMO1. The ancestral G allele that is associated with tumour formation resides in a conserved GATA transcription factor binding motif. We show that the newly evolved protective TATA allele is associated with decreased total LMO1 expression (P = 0.028) in neuroblastoma primary tumours, and ablates GATA3 binding (P < 0.0001). We demonstrate allelic imbalance favouring the G-containing strand in tumours heterozygous for this SNP, as demonstrated both by RNA sequencing (P < 0.0001) and reporter assays (P = 0.002). These findings indicate that a recently evolved polymorphism within a super-enhancer element in the first intron of LMO1 influences neuroblastoma susceptibility through differential GATA transcription factor binding and direct modulation of LMO1 expression in cis, and this leads to an oncogenic dependency in tumour cells.

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.

BIDCHIPS: bias decomposition and removal from ChIP-seq data clarifies true binding signal and its functional correlates

Background

Unraveling transcriptional regulatory networks is a central problem in molecular biology and, in this quest, chromatin immunoprecipitation and sequencing (ChIP-seq) technology has given us the unprecedented ability to identify sites of protein-DNA binding and histone modification genome wide. However, multiple systemic and procedural biases hinder harnessing the full potential of this technology. Previous studies have addressed this problem, but a thorough characterization of different, interacting biases on ChIP-seq signals is still lacking.

Results

Here, we present a novel framework where the genome-wide ChIP-seq signal is viewed as being quantifiably influenced by different, measurable sources of bias, which can then be computationally subtracted away. We use a compendium of 123 human ENCODE ChIP-seq datasets to build regression models that tell us how much of a ChIP-seq signal can be attributed to mappability, GC-content, chromatin accessibility, and factors represented in input DNA and IgG controls. When we use the model to separate out these non-binding influences from the ChIP-seq signal, we obtain a purified signal that associates better to TF-DNA-binding motifs than do other measures of peak significance. We also carry out a multiscale analysis that reveals how ChIP-seq signal biases differ across different scales. Finally, we investigate previously reported associations between gene expression and ChIP-seq signals at transcription start sites. We show that our model can be used to discriminate ChIP-seq signals that are truly related to gene expression from those that are merely correlated by virtue of bias—in particular, chromatin accessibility bias, which shows up in ChIP-seq signals and also relates to gene expression.

Conclusions

Our study provides new insights into the behavior of ChIP-seq signal biases and proposes a novel mitigation framework that improves results compared to existing techniques. With ChIP-seq now being the central technology for studying transcriptional regulation, it is most crucial to accurately characterize, quantify, and adjust for the genome-wide effects of biases affecting ChIP-seq. Our study also emphasizes that properly accounting for confounders in ChIP-seq data is of paramount importance for obtaining biologically accurate insights into the workings of the complex regulatory mechanisms in living organisms. R and MATLAB packages implementing the framework can be obtained from http://www.perkinslab.ca/Software.html.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0028-2) contains supplementary material, which is available to authorized users.

Genetic Variation Determines PPARγ Function and Anti-diabetic Drug Response In Vivo

SNPs affecting disease risk often reside in non-coding genomic regions. Here, we show that SNPs are highly enriched at mouse strain-selective adipose tissue binding sites for PPARγ, a nuclear receptor for anti-diabetic drugs. Many such SNPs alter binding motifs for PPARγ or cooperating factors and functionally regulate nearby genes whose expression is strain selective and imbalanced in heterozygous F1 mice. Moreover, genetically determined binding of PPARγ accounts for mouse strain-specific transcriptional effects of TZD drugs, providing proof of concept for personalized medicine related to nuclear receptor genomic occupancy. In human fat, motif-altering SNPs cause differential PPARγ binding, provide a molecular mechanism for some expression quantitative trait loci, and are risk factors for dysmetabolic traits in genome-wide association studies. One PPARγ motif-altering SNP is associated with HDL levels and other metabolic syndrome parameters. Thus, natural genetic variation in PPARγ genomic occupancy determines individual disease risk and drug response.

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.

Genetic dissection of ozone tolerance in rice (Oryza sativa L.) by a genome-wide association study

Tropospheric ozone causes various negative effects on plants and affects the yield and quality of agricultural crops. Here, we report a genome-wide association study (GWAS) in rice (Oryza sativa L.) to determine candidate loci associated with ozone tolerance. A diversity panel consisting of 328 accessions representing all subgroups of O. sativa was exposed to ozone stress at 60 nl l(-1) for 7h every day throughout the growth season, or to control conditions. Averaged over all genotypes, ozone significantly affected biomass-related traits (plant height -1.0%, shoot dry weight -15.9%, tiller number -8.3%, grain weight -9.3%, total panicle weight -19.7%, single panicle weight -5.5%) and biochemical/physiological traits (symptom formation, SPAD value -4.4%, foliar lignin content +3.4%). A wide range of genotypic variance in response to ozone stress were observed in all phenotypes. Association mapping based on more than 30 000 single-nucleotide polymorphism (SNP) markers yielded 16 significant markers throughout the genome by applying a significance threshold of P<0.0001. Furthermore, by determining linkage disequilibrium blocks associated with significant SNPs, we gained a total of 195 candidate genes for these traits. The following sequence analysis revealed a number of novel polymorphisms in two candidate genes for the formation of visible leaf symptoms, a RING and an EREBP gene, both of which are involved in cell death and stress defence reactions. This study demonstrated substantial natural variation of responses to ozone in rice and the possibility of using GWAS in elucidating the genetic factors underlying ozone tolerance.

Cancer by super-enhancer

Tiny changes in our genomes can enhance oncogene expression and contribute to tumorigenesis [Also see Report by Mansour et al. ]

Single nucleotide polymorphisms of the FTO gene and cancer risk: an overview

The FTO (fat mass and obesity-associated) gene has a strong linkage disequilibrium block, within which SNPs have been identified that are involved in the development of obesity. Recently some of these variants have also been associated with cancer. However, identification of the possible mechanisms that could explain these associations has proven to be elusive. It has been found that FTO polymorphisms can regulate the expression of genes at large kilobases of distance as well as the expression of the FTO gene itself, and regions for transcription factor binding. To date it has been observed that variants rs9939609, rs17817449, rs8050136, rs1477196, rs6499640, rs16953002, rs11075995 and rs1121980 are associated with the risk of developing cancer. Some studies have produced negative results when comparing the same polymorphisms, but make a simple association between polymorphic variants and cancer, have proved difficult because this relation is by nature multifactorial. A certain degree of variation resulting from the improper design of studies or processing of data can lead to erroneous conclusions. However, it is now unquestionable that certain FTO polymorphisms regulate genetic expression related to cancer susceptibility, although this field is just beginning to be understood.

Genome-wide identification of long noncoding RNAs in rat models of cardiovascular and renal disease

Long noncoding RNAs (lncRNAs) are an emerging class of genomic regulatory molecules reported in various species. In the rat, which is one of the major mammalian model organisms, discovery of lncRNAs on a genome-wide scale is lagging. Renal lncRNA sequencing and lncRNA transcriptome analysis were conducted in 3 rat strains that are widely used in cardiovascular and renal research: the Dahl salt-sensitive rat, the spontaneously hypertensive rat, and the Dahl salt-resistant rat. Through the RNA sequencing approach, 3273 transcripts were identified as rat lncRNAs. A majority of lncRNAs were without predicted target genes. Differential expression of 273 and 749 lncRNAs was detected between Dahl salt-sensitive versus Dahl salt-resistant and Dahl salt-sensitive versus spontaneously hypertensive rat comparisons, respectively. To couple the observed differential expression of lncRNAs with the status of mRNAs, an mRNA transcriptome analysis was conducted. Several cis mRNA genes were coregulated with lncRNAs. Of these, the protein expression status of 4 target genes, Asb3, Chac2, Pex11b, and Sp5, were differentially expressed between the relevant strain comparisons, thereby suggesting that the differentially expressed lncRNAs associated with these genes are candidate genetic determinants of blood pressure. This study serves as a first-generation catalog of rat lncRNAs and illustrates the prioritization of lncRNAs as candidates for complex polygenic traits.

Frequent chromosomal aberrations and candidate genes in head and neck squamous cell carcinoma

The knowledge of the biology of head and neck squamous cell carcinoma (HNSCC) has had relatively little impact on the improvement in oncologic outcome up to date. However, the identification of oncogenes and tumor suppressor genes (TSGs) involved in cancer progression contributes to the understanding of the molecular pathways involved in oncogenesis and could contribute to individual risk assessment and provide tools for improvement of treatment and targets for therapy based on the alterations in these pathways. The aim of this article is to review the chromosomal aberrations commonly found in HNSCC, to identify the genes in these chromosomal regions suggested to act as (candidate) oncogenes or TSGs, and to discuss the molecular mechanisms modulating their expression.

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.

In pursuit of design principles of regulatory sequences

Instructions for when, where and to what level each gene should be expressed are encoded within regulatory sequences. The importance of motifs recognized by DNA-binding regulators has long been known, but their extensive characterization afforded by recent technologies only partly accounts for how regulatory instructions are encoded in the genome. Here, we review recent advances in our understanding of regulatory sequences that influence transcription and go beyond the description of motifs. We discuss how understanding different aspects of the sequence-encoded regulation can help to unravel the genotype-phenotype relationship, which would lead to a more accurate and mechanistic interpretation of personal genome sequences.

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.

Myc induced replicative stress response: How to cope with it and exploit it

Myc is a cellular oncogene frequently deregulated in cancer that has the ability to stimulate cellular growth by promoting a number of proliferative and pro-survival pathways. Here we will focus on how Myc controls a number of diverse cellular processes that converge to ensure processivity and robustness of DNA synthesis, thus preventing the inherent replicative stress responses usually evoked by oncogenic lesions. While these processes provide cancer cells with a long-term proliferative advantage, they also represent cancer liabilities that can be exploited to devise innovative therapeutic approaches to target Myc overexpressing tumors. This article is part of a Special Issue entitled: Myc proteins in cell biology and pathology.