Genome-Wide Association Study of Obsessive-Compulsive Symptoms including 33,943 individuals from the general population.

While 1-2% of individuals meet the criteria for a clinical diagnosis of obsessive-compulsive disorder (OCD), many more (~13-38%) experience subclinical obsessive-compulsive symptoms (OCS) during their life. To characterize the genetic underpinnings of OCS and its genetic relationship to OCD, we conducted the largest genome-wide association study (GWAS) meta-analysis of parent- or self-reported OCS to date (N = 33,943 with complete phenotypic and genome-wide data), combining the results from seven large-scale population-based cohorts from Sweden, the Netherlands, England, and Canada (including six twin cohorts and one cohort of unrelated individuals). We found no genome-wide significant associations at the single-nucleotide polymorphism (SNP) or gene-level, but a polygenic risk score (PRS) based on the OCD GWAS previously published by the Psychiatric Genetics Consortium (PGC-OCD) was significantly associated with OCS (Pfixed = 3.06 × 10-5). Also, one curated gene set (Mootha Gluconeogenesis) reached Bonferroni-corrected significance (Ngenes = 28, Beta = 0.79, SE = 0.16, Pbon = 0.008). Expression of genes in this set is high at sites of insulin mediated glucose disposal. Dysregulated insulin signaling in the etiology of OCS has been suggested by a previous study describing a genetic overlap of OCS with insulin signaling-related traits in children and adolescents. We report a SNP heritability of 4.1% (P = 0.0044) in the meta-analyzed GWAS, and heritability estimates based on the twin cohorts of 33-43%. Genetic correlation analysis showed that OCS were most strongly associated with OCD (rG = 0.72, p = 0.0007) among all tested psychiatric disorders (N = 11). Of all 97 tested phenotypes, 24 showed a significant genetic correlation with OCS, and 66 traits showed concordant directions of effect with OCS and OCD. OCS have a significant polygenic contribution and share genetic risk with diagnosed OCD, supporting the hypothesis that OCD represents the extreme end of widely distributed OCS in the population.

Trans-ancestry meta-analysis of genome wide association studies of inhibitory control.

Deficits in effective executive function, including inhibitory control are associated with risk for a number of psychiatric disorders and significantly impact everyday functioning. These complex traits have been proposed to serve as endophenotypes, however, their genetic architecture is not yet well understood. To identify the common genetic variation associated with inhibitory control in the general population we performed the first trans-ancestry genome wide association study (GWAS) combining data across 8 sites and four ancestries (N = 14,877) using cognitive traits derived from the stop-signal task, namely - go reaction time (GoRT), go reaction time variability (GoRT SD) and stop signal reaction time (SSRT). Although we did not identify genome wide significant associations for any of the three traits, GoRT SD and SSRT demonstrated significant and similar SNP heritability of 8.2%, indicative of an influence of genetic factors. Power analyses demonstrated that the number of common causal variants contributing to the heritability of these phenotypes is relatively high and larger sample sizes are necessary to robustly identify associations. In Europeans, the polygenic risk for ADHD was significantly associated with GoRT SD and the polygenic risk for schizophrenia was associated with GoRT, while in East Asians polygenic risk for schizophrenia was associated with SSRT. These results support the potential of executive function measures as endophenotypes of neuropsychiatric disorders. Together these findings provide the first evidence indicating the influence of common genetic variation in the genetic architecture of inhibitory control quantified using objective behavioural traits derived from the stop-signal task.

Meta-analysis of genome-wide association studies of hoarding symptoms in 27,537 individuals.

Hoarding Disorder (HD) is a mental disorder characterized by persistent difficulties discarding or parting with possessions, often resulting in cluttered living spaces, distress, and impairment. Its etiology is largely unknown, but twin studies suggest that it is moderately heritable. In this study, we pooled phenotypic and genomic data from seven international cohorts (N = 27,537 individuals) and conducted a genome wide association study (GWAS) meta-analysis of parent- or self-reported hoarding symptoms (HS). We followed up the results with gene-based and gene-set analyses, as well as leave-one-out HS polygenic risk score (PRS) analyses. To examine a possible genetic association between hoarding symptoms and other phenotypes we conducted cross-trait PRS analyses. Though we did not report any genome-wide significant SNPs, we report heritability estimates for the twin-cohorts between 26-48%, and a SNP-heritability of 11% for an unrelated sub-cohort. Cross-trait PRS analyses showed that the genetic risk for schizophrenia and autism spectrum disorder were significantly associated with hoarding symptoms. We also found suggestive evidence for an association with educational attainment. There were no significant associations with other phenotypes previously linked to HD, such as obsessive-compulsive disorder, depression, anxiety, or attention-deficit hyperactivity disorder. To conclude, we found that HS are heritable, confirming and extending previous twin studies but we had limited power to detect any genome-wide significant loci. Much larger samples will be needed to further extend these findings and reach a "gene discovery zone". To move the field forward, future research should not only include genetic analyses of quantitative hoarding traits in larger samples, but also in samples of individuals meeting strict diagnostic criteria for HD, and more ethnically diverse samples.

Genetic Regulation of DNA Methylation Yields Novel Discoveries in GWAS of Colorectal Cancer.

BACKGROUND: Colorectal cancer has a strong epigenetic component that is accompanied by frequent DNA methylation (DNAm) alterations in addition to heritable genetic risk. It is of interest to understand the interrelationship of germline genetics, DNAm, and colorectal cancer risk. METHODS: We performed a genome-wide methylation quantitative trait locus (meQTL) analysis in 1,355 people, assessing the pairwise associations between genetic variants and lymphocytes methylation data. In addition, we used penalized regression with cis-genetic variants ± 1 Mb of methylation to identify genome-wide heritable DNAm. We evaluated the association of genetically predicted methylation with colorectal cancer risk based on genome-wide association studies (GWAS) of over 125,000 cases and controls using the multivariate sMiST as well as univariately via examination of marginal association with colorectal cancer risk. RESULTS: Of the 142 known colorectal cancer GWAS loci, 47 were identified as meQTLs. We identified four novel colorectal cancer-associated loci (NID2, ATXN10, KLHDC10, and CEP41) that reside over 1 Mb outside of known colorectal cancer loci and 10 secondary signals within 1 Mb of known loci. CONCLUSIONS: Leveraging information of DNAm regulation into genetic association of colorectal cancer risk reveals novel pathways in colorectal cancer tumorigenesis. Our summary statistics-based framework sMiST provides a powerful approach by combining information from the effect through methylation and residual direct effects of the meQTLs on disease risk. Further validation and functional follow-up of these novel pathways are needed. IMPACT: Using genotype, DNAm, and GWAS, we identified four new colorectal cancer risk loci. We studied the landscape of genetic regulation of DNAm via single-SNP and multi-SNP meQTL analyses.

Genome-wide association study of pediatric obsessive-compulsive traits: shared genetic risk between traits and disorder.

Using a novel trait-based measure, we examined genetic variants associated with obsessive-compulsive (OC) traits and tested whether OC traits and obsessive-compulsive disorder (OCD) shared genetic risk. We conducted a genome-wide association analysis (GWAS) of OC traits using the Toronto Obsessive-Compulsive Scale (TOCS) in 5018 unrelated Caucasian children and adolescents from the community (Spit for Science sample). We tested the hypothesis that genetic variants associated with OC traits from the community would be associated with clinical OCD using a meta-analysis of all currently available OCD cases. Shared genetic risk was examined between OC traits and OCD in the respective samples using polygenic risk score and genetic correlation analyses. A locus tagged by rs7856850 in an intron of PTPRD (protein tyrosine phosphatase δ) was significantly associated with OC traits at the genome-wide significance level (p = 2.48 × 10-8). rs7856850 was also associated with OCD in a meta-analysis of OCD case/control genome-wide datasets (p = 0.0069). The direction of effect was the same as in the community sample. Polygenic risk scores from OC traits were significantly associated with OCD in case/control datasets and vice versa (p's < 0.01). OC traits were highly, but not significantly, genetically correlated with OCD (rg = 0.71, p = 0.062). We report the first validated genome-wide significant variant for OC traits in PTPRD, downstream of the most significant locus in a previous OCD GWAS. OC traits measured in the community sample shared genetic risk with OCD case/control status. Our results demonstrate the feasibility and power of using trait-based approaches in community samples for genetic discovery.

Obsessive-compulsive disorder and attention-deficit/hyperactivity disorder: distinct associations with DNA methylation and genetic variation.

BACKGROUND: A growing body of research has demonstrated associations between specific neurodevelopmental disorders and variation in DNA methylation (DNAm), implicating this molecular mark as a possible contributor to the molecular etiology of these disorders and/or as a novel disease biomarker. Furthermore, genetic risk variants of neurodevelopmental disorders have been found to be enriched at loci associated with DNAm patterns, referred to as methylation quantitative trait loci (mQTLs). METHODS: We conducted two epigenome-wide association studies in individuals with attention-deficit/hyperactivity disorder (ADHD) or obsessive-compulsive disorder (OCD) (aged 4-18 years) using DNA extracted from saliva. DNAm data generated on the Illumina Human Methylation 450 K array were used to examine the interaction between genetic variation and DNAm patterns associated with these disorders. RESULTS: Using linear regression followed by principal component analysis, individuals with the most endorsed symptoms of ADHD or OCD were found to have significantly more distinct DNAm patterns from controls, as compared to all cases. This suggested that the phenotypic heterogeneity of these disorders is reflected in altered DNAm at specific sites. Further investigations of the DNAm sites associated with each disorder revealed that despite little overlap of these DNAm sites across the two disorders, both disorders were significantly enriched for mQTLs within our sample. CONCLUSIONS: Our DNAm data provide insights into the regulatory changes associated with genetic variation, highlighting their potential utility both in directing GWAS and in elucidating the pathophysiology of neurodevelopmental disorders.

Landscape of somatic single nucleotide variants and indels in colorectal cancer and impact on survival.

Colorectal cancer (CRC) is a biologically heterogeneous disease. To characterize its mutational profile, we conduct targeted sequencing of 205 genes for 2,105 CRC cases with survival data. Our data shows several findings in addition to enhancing the existing knowledge of CRC. We identify PRKCI, SPZ1, MUTYH, MAP2K4, FETUB, and TGFBR2 as additional genes significantly mutated in CRC. We find that among hypermutated tumors, an increased mutation burden is associated with improved CRC-specific survival (HR = 0.42, 95% CI: 0.21-0.82). Mutations in TP53 are associated with poorer CRC-specific survival, which is most pronounced in cases carrying TP53 mutations with predicted 0% transcriptional activity (HR = 1.53, 95% CI: 1.21-1.94). Furthermore, we observe differences in mutational frequency of several genes and pathways by tumor location, stage, and sex. Overall, this large study provides deep insights into somatic mutations in CRC, and their potential relationships with survival and tumor features.

Shared genetic etiology between obsessive-compulsive disorder, obsessive-compulsive symptoms in the population, and insulin signaling.

Obsessive-compulsive symptoms (OCS) in the population have been linked to obsessive-compulsive disorder (OCD) in genetic and epidemiological studies. Insulin signaling has been implicated in OCD. We extend previous work by assessing genetic overlap between OCD, population-based OCS, and central nervous system (CNS) and peripheral insulin signaling. We conducted genome-wide association studies (GWASs) in the population-based Philadelphia Neurodevelopmental Cohort (PNC, 650 children and adolescents) of the total OCS score and six OCS factors from an exploratory factor analysis of 22 questions. Subsequently, we performed polygenic risk score (PRS)-based analysis to assess shared genetic etiologies between clinical OCD (using GWAS data from the Psychiatric Genomics Consortium), the total OCS score and OCS factors. We then performed gene-set analyses with a set of OCD-linked genes centered around CNS insulin-regulated synaptic function and PRS-based analyses for five peripheral insulin signaling-related traits. For validation purposes, we explored data from the independent Spit for Science population cohort (5,047 children and adolescents). In the PNC, we found a significant shared genetic etiology between OCD and 'guilty taboo thoughts'. In the Spit for Science cohort, we additionally observed genetic sharing between 'symmetry/counting/ordering' and 'contamination/cleaning'. The CNS insulin-linked gene-set also associated with 'symmetry/counting/ordering' in the PNC. Further, we identified genetic sharing between peripheral insulin signaling-related traits: type 2 diabetes with 'aggressive taboo thoughts', and levels of fasting insulin and 2 h glucose with OCD. In conclusion, OCD, OCS in the population and insulin-related traits share genetic risk factors, indicating a common etiological mechanism underlying somatic and psychiatric disorders.

Mendelian Randomization of Circulating Polyunsaturated Fatty Acids and Colorectal Cancer Risk.

BACKGROUND: Results from epidemiologic studies examining polyunsaturated fatty acids (PUFA) and colorectal cancer risk are inconsistent. Mendelian randomization may strengthen causal inference from observational studies. Given their shared metabolic pathway, examining the combined effects of aspirin/NSAID use with PUFAs could help elucidate an association between PUFAs and colorectal cancer risk. METHODS: Information was leveraged from genome-wide association studies (GWAS) regarding PUFA-associated SNPs to create weighted genetic scores (wGS) representing genetically predicted circulating blood PUFAs for 11,016 non-Hispanic white colorectal cancer cases and 13,732 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). Associations per SD increase in the wGS were estimated using unconditional logistic regression. Interactions between PUFA wGSs and aspirin/NSAID use on colorectal cancer risk were also examined. RESULTS: Modest colorectal cancer risk reductions were observed per SD increase in circulating linoleic acid [ORLA = 0.96; 95% confidence interval (CI) = 0.93-0.98; P = 5.2 × 10-4] and α-linolenic acid (ORALA = 0.95; 95% CI = 0.92-0.97; P = 5.4 × 10-5), whereas modest increased risks were observed for arachidonic (ORAA = 1.06; 95% CI = 1.03-1.08; P = 3.3 × 10-5), eicosapentaenoic (OREPA = 1.04; 95% CI = 1.01-1.07; P = 2.5 × 10-3), and docosapentaenoic acids (ORDPA = 1.03; 95% CI = 1.01-1.06; P = 1.2 × 10-2). Each of these effects was stronger among aspirin/NSAID nonusers in the stratified analyses. CONCLUSIONS: Our study suggests that higher circulating shorter-chain PUFAs (i.e., LA and ALA) were associated with reduced colorectal cancer risk, whereas longer-chain PUFAs (i.e., AA, EPA, and DPA) were associated with an increased colorectal cancer risk. IMPACT: The interaction of PUFAs with aspirin/NSAID use indicates a shared colorectal cancer inflammatory pathway. Future research should continue to improve PUFA genetic instruments to elucidate the independent effects of PUFAs on colorectal cancer.

Whole genomes define concordance of matched primary, xenograft, and organoid models of pancreas cancer.

Pancreatic ductal adenocarcinoma (PDAC) has the worst prognosis among solid malignancies and improved therapeutic strategies are needed to improve outcomes. Patient-derived xenografts (PDX) and patient-derived organoids (PDO) serve as promising tools to identify new drugs with therapeutic potential in PDAC. For these preclinical disease models to be effective, they should both recapitulate the molecular heterogeneity of PDAC and validate patient-specific therapeutic sensitivities. To date however, deep characterization of the molecular heterogeneity of PDAC PDX and PDO models and comparison with matched human tumour remains largely unaddressed at the whole genome level. We conducted a comprehensive assessment of the genetic landscape of 16 whole-genome pairs of tumours and matched PDX, from primary PDAC and liver metastasis, including a unique cohort of 5 'trios' of matched primary tumour, PDX, and PDO. We developed a pipeline to score concordance between PDAC models and their paired human tumours for genomic events, including mutations, structural variations, and copy number variations. Tumour-model comparisons of mutations displayed single-gene concordance across major PDAC driver genes, but relatively poor agreement across the greater mutational load. Genome-wide and chromosome-centric analysis of structural variation (SV) events highlights previously unrecognized concordance across chromosomes that demonstrate clustered SV events. We found that polyploidy presented a major challenge when assessing copy number changes; however, ploidy-corrected copy number states suggest good agreement between donor-model pairs. Collectively, our investigations highlight that while PDXs and PDOs may serve as tractable and transplantable systems for probing the molecular properties of PDAC, these models may best serve selective analyses across different levels of genomic complexity.

Integration of Genomic and Transcriptional Features in Pancreatic Cancer Reveals Increased Cell Cycle Progression in Metastases.

We integrated clinical, genomic, and transcriptomic data from 224 primaries and 95 metastases from 289 patients to characterize progression of pancreatic ductal adenocarcinoma (PDAC). Driver gene alterations and mutational and expression-based signatures were preserved, with truncations, inversions, and translocations most conserved. Cell cycle progression (CCP) increased with sequential inactivation of tumor suppressors, yet remained higher in metastases, perhaps driven by cell cycle regulatory gene variants. Half of the cases were hypoxic by expression markers, overlapping with molecular subtypes. Paired tumor heterogeneity showed cancer cell migration by Halstedian progression. Multiple PDACs arising synchronously and metachronously in the same pancreas were actually intra-parenchymal metastases, not independent primary tumors. Established clinical co-variates dominated survival analyses, although CCP and hypoxia may inform clinical practice.

Diagnostics for Pleiotropy in Mendelian Randomization Studies: Global and Individual Tests for Direct Effects.

Diagnosing pleiotropy is critical for assessing the validity of Mendelian randomization (MR) analyses. The popular MR-Egger method evaluates whether there is evidence of bias-generating pleiotropy among a set of candidate genetic instrumental variables. In this article, we propose a statistical method-global and individual tests for direct effects (GLIDE)-for systematically evaluating pleiotropy among the set of genetic variants (e.g., single nucleotide polymorphisms (SNPs)) used for MR. As a global test, simulation experiments suggest that GLIDE is nearly uniformly more powerful than the MR-Egger method. As a sensitivity analysis, GLIDE is capable of detecting outliers in individual variant-level pleiotropy, in order to obtain a refined set of genetic instrumental variables. We used GLIDE to analyze both body mass index and height for associations with colorectal cancer risk in data from the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colon Cancer Family Registry (multiple studies). Among the body mass index-associated SNPs and the height-associated SNPs, several individual variants showed evidence of pleiotropy. Removal of these potentially pleiotropic SNPs resulted in attenuation of respective estimates of the causal effects. In summary, the proposed GLIDE method is useful for sensitivity analyses and improves the validity of MR.

Mendelian randomisation study of age at menarche and age at menopause and the risk of colorectal cancer.

BACKGROUND: Substantial evidence supports an association between use of menopausal hormone therapy and decreased colorectal cancer (CRC) risk, indicating a role of exogenous sex hormones in CRC development. However, findings on endogenous oestrogen exposure and CRC are inconsistent. METHODS: We used a Mendelian randomisation approach to test for a causal effect of age at menarche and age at menopause as surrogates for endogenous oestrogen exposure on CRC risk. Weighted genetic risk scores based on 358 single-nucleotide polymorphisms associated with age at menarche and 51 single-nucleotide polymorphisms associated with age at menopause were used to estimate the association with CRC risk using logistic regression in 12,944 women diagnosed with CRC and 10,741 women without CRC from three consortia. Sensitivity analyses were conducted to address pleiotropy and possible confounding by body mass index. RESULTS: Genetic risk scores for age at menarche (odds ratio per year 0.98, 95% confidence interval: 0.95-1.02) and age at menopause (odds ratio 0.98, 95% confidence interval: 0.94-1.01) were not significantly associated with CRC risk. The sensitivity analyses yielded similar results. CONCLUSIONS: Our study does not support a causal relationship between genetic risk scores for age at menarche and age at menopause and CRC risk.

Exome-Wide Association Study of Pancreatic Cancer Risk.

We conducted a case-control exome-wide association study to discover germline variants in coding regions that affect risk for pancreatic cancer, combining data from 5 studies. We analyzed exome and genome sequencing data from 437 patients with pancreatic cancer (cases) and 1922 individuals not known to have cancer (controls). In the primary analysis, BRCA2 had the strongest enrichment for rare inactivating variants (17/437 cases vs 3/1922 controls) (P = 3.27x10-6; exome-wide statistical significance threshold P < 2.5x10-6). Cases had more rare inactivating variants in DNA repair genes than controls, even after excluding 13 genes known to predispose to pancreatic cancer (adjusted odds ratio, 1.35; P = .045). At the suggestive threshold (P < .001), 6 genes were enriched for rare damaging variants (UHMK1, AP1G2, DNTA, CHST6, FGFR3, and EPHA1) and 7 genes had associations with pancreatic cancer risk, based on the sequence-kernel association test. We confirmed variants in BRCA2 as the most common high-penetrant genetic factor associated with pancreatic cancer and we also identified candidate pancreatic cancer genes. Large collaborations and novel approaches are needed to overcome the genetic heterogeneity of pancreatic cancer predisposition.

Genomics-Driven Precision Medicine for Advanced Pancreatic Cancer: Early Results from the COMPASS Trial.

Purpose: To perform real-time whole genome sequencing (WGS) and RNA sequencing (RNASeq) of advanced pancreatic ductal adenocarcinoma (PDAC) to identify predictive mutational and transcriptional features for better treatment selection.Experimental Design: Patients with advanced PDAC were prospectively recruited prior to first-line combination chemotherapy. Fresh tumor tissue was acquired by image-guided percutaneous core biopsy for WGS and RNASeq. Laser capture microdissection was performed for all cases. Primary endpoint was feasibility to report WGS results prior to first disease assessment CT scan at 8 weeks. The main secondary endpoint was discovery of patient subsets with predictive mutational and transcriptional signatures.Results: Sixty-three patients underwent a tumor biopsy between December 2015 and June 2017. WGS and RNASeq were successful in 62 (98%) and 60 (95%), respectively. Genomic results were reported at a median of 35 days (range, 19-52 days) from biopsy, meeting the primary feasibility endpoint. Objective responses to first-line chemotherapy were significantly better in patients with the classical PDAC RNA subtype compared with those with the basal-like subtype (P = 0.004). The best progression-free survival was observed in those with classical subtype treated with m-FOLFIRINOX. GATA6 expression in tumor measured by RNA in situ hybridization was found to be a robust surrogate biomarker for differentiating classical and basal-like PDAC subtypes. Potentially actionable genetic alterations were found in 30% of patients.Conclusions: Prospective genomic profiling of advanced PDAC is feasible, and our early data indicate that chemotherapy response differs among patients with different genomic/transcriptomic subtypes. Clin Cancer Res; 24(6); 1344-54. ©2017 AACR.

Enrichment of colorectal cancer associations in functional regions: Insight for using epigenomics data in the analysis of whole genome sequence-imputed GWAS data.

BACKGROUND: The evaluation of less frequent genetic variants and their effect on complex disease pose new challenges for genomic research. To investigate whether epigenetic data can be used to inform aggregate rare-variant association methods (RVAM), we assessed whether variants more significantly associated with colorectal cancer (CRC) were preferentially located in non-coding regulatory regions, and whether enrichment was specific to colorectal tissues. METHODS: Active regulatory elements (ARE) were mapped using data from 127 tissues and cell-types from NIH Roadmap Epigenomics and Encyclopedia of DNA Elements (ENCODE) projects. We investigated whether CRC association p-values were more significant for common variants inside versus outside AREs, or 2) inside colorectal (CR) AREs versus AREs of other tissues and cell-types. We employed an integrative epigenomic RVAM for variants with allele frequency <1%. Gene sets were defined as ARE variants within 200 kilobases of a transcription start site (TSS) using either CR ARE or ARE from non-digestive tissues. CRC-set association p-values were used to evaluate enrichment of less frequent variant associations in CR ARE versus non-digestive ARE. RESULTS: ARE from 126/127 tissues and cell-types were significantly enriched for stronger CRC-variant associations. Strongest enrichment was observed for digestive tissues and immune cell types. CR-specific ARE were also enriched for stronger CRC-variant associations compared to ARE combined across non-digestive tissues (p-value = 9.6 × 10-4). Additionally, we found enrichment of stronger CRC association p-values for rare variant sets of CR ARE compared to non-digestive ARE (p-value = 0.029). CONCLUSIONS: Integrative epigenomic RVAM may enable discovery of less frequent variants associated with CRC, and ARE of digestive and immune tissues are most informative. Although distance-based aggregation of less frequent variants in CR ARE surrounding TSS showed modest enrichment, future association studies would likely benefit from joint analysis of transcriptomes and epigenomes to better link regulatory variation with target genes.

Blood triglyceride levels are associated with DNA methylation at the serine metabolism gene PHGDH.

Efficient interventions to reduce blood triglycerides are few; newer and more tolerable intervention targets are needed. Understanding the molecular mechanisms underlying blood triglyceride levels variation is key to identifying new therapies. To explore the role of epigenetic mechanisms on triglyceride levels, a blood methylome scan was conducted in 199 individuals from 5 French-Canadian families ascertained on venous thromboembolism, and findings were replicated in 324 French unrelated patients with venous thromboembolism. Genetic context and functional relevance were investigated. Two DNA methylation sites associated with triglyceride levels were identified. The first one, located in the ABCG1 gene, was recently reported, whereas the second one, located in the promoter of the PHGDH gene, is novel. The PHGDH methylation site, cg14476101, was found to be associated with variation in triglyceride levels in a threshold manner: cg14476101 was inversely associated with triglyceride levels only when triglyceride levels were above 1.12 mmol/L (discovery P-value = 8.4 × 10-6; replication P-value = 0.0091). Public databases findings supported a functional role of cg14476101 on PHGDH expression. PHGDH catalyses the first step in the serine biosynthesis pathway. These findings highlight the role of epigenetic regulation of the PHGDH gene in triglyceride metabolism, providing novel insights on putative intervention targets.

In silico pathway analysis and tissue specific cis-eQTL for colorectal cancer GWAS risk variants.

BACKGROUND: Genome-wide association studies have identified 55 genetic variants associated with colorectal cancer risk to date. However, potential causal genes and pathways regulated by these risk variants remain to be characterized. Therefore, we performed gene ontology enrichment and pathway analyses to determine if there was an enrichment of genes in proximity to the colorectal cancer risk variants that could further elucidate the probable causal genes and pathways involved in colorectal cancer biology. RESULTS: For the 65 unique genes that either contained, or were immediately neighboring up- and downstream, of these variants there was a significant enrichment for the KEGG pathway, Pathways in Cancer (p-value = 2.67 × 10-5) and an enrichment for multiple biological processes (FDR < 0.05), such as cell junction organization, tissue morphogenesis, regulation of SMAD protein phosphorylation, and odontogenesis identified through Gene Ontology analysis. To identify potential causal genes, we conducted a cis-expression quantitative trait loci (cis-eQTL) analysis using gene expression and genotype data from the Genotype-Tissue Expression (GTEx) Project portal in normal sigmoid (n = 124) and transverse (n = 169) colon tissue. In addition, we also did a cis-eQTL analysis on colorectal tumor tissue (n = 147) from The Cancer Genome Atlas (TCGA). We identified two risk alleles that were significant cis-eQTLs for FADS2 (rs1535) and COLCA1 and 2 (rs3802842) genes in the normal transverse colon tissue and two risk alleles that were significant cis-eQTLs for the CABLES2 (rs2427308) and LIPG (rs7229639) genes in the normal sigmoid colon tissue, but not tumor tissue. CONCLUSIONS: Our data reaffirm the potential to identify an enrichment for biological processes and candidate causal genes based on expression profiles correlated with genetic risk alleles of colorectal cancer, however, the identification of these significant cis-eQTLs is context and tissue specific.

The dynamic DNA methylation landscape of the mutL homolog 1 shore is altered by MLH1-93G>A polymorphism in normal tissues and colorectal cancer.

BACKGROUND: Colorectal cancers (CRCs) undergo distinct genetic and epigenetic alterations. Expression of mutL homolog 1 (MLH1), a mismatch repair gene that corrects DNA replication errors, is lost in up to 15% of sporadic tumours due to mutation or, more commonly, due to DNA methylation of its promoter CpG island. A single nucleotide polymorphism (SNP) in the CpG island of MLH1 (MLH1-93G>A or rs1800734) is associated with CpG island hypermethylation and decreased MLH1 expression in CRC tumours. Further, in peripheral blood mononuclear cell (PBMC) DNA of both CRC cases and non-cancer controls, the variant allele of rs1800734 is associated with hypomethylation at the MLH1 shore, a region upstream of its CpG island that is less dense in CpG sites. RESULTS: To determine whether this genotype-epigenotype association is present in other tissue types, including colorectal tumours, we assessed DNA methylation in matched normal colorectal tissue, tumour, and PBMC DNA from 349 population-based CRC cases recruited from the Ontario Familial Colorectal Cancer Registry. Using the semi-quantitative real-time PCR-based MethyLight assay, MLH1 shore methylation was significantly higher in tumour tissue than normal colon or PBMCs (P < 0.01). When shore methylation levels were stratified by SNP genotype, normal colorectal DNA and PBMC DNA were significantly hypomethylated in association with variant SNP genotype (P < 0.05). However, this association was lost in tumour DNA. Among distinct stages of CRC, metastatic stage IV CRC tumours incurred significant hypomethylation compared to stage I-III cases, irrespective of genotype status. Shore methylation of MLH1 was not associated with MSI status or promoter CpG island hypermethylation, regardless of genotype. To confirm these results, bisulfite sequencing was performed in matched tumour and normal colorectal specimens from six CRC cases, including two cases per genotype (wildtype, heterozygous, and homozygous variant). Bisulfite sequencing results corroborated the methylation patterns found by MethyLight, with significant hypomethylation in normal colorectal tissue of variant SNP allele carriers. CONCLUSIONS: These results indicate that the normal tissue types tested (colorectum and PBMC) experience dynamic genotype-associated epigenetic alterations at the MLH1 shore, whereas tumour DNA incurs aberrant hypermethylation compared to normal DNA.