METHODS: We conducted a two-sample Mendelian randomization (MR) study to examine the genetically predicted effects of epigenetic age acceleration as measured by HannumAge (nine single-nucleotide polymorphisms (SNPs)), Horvath Intrinsic Age (24 SNPs), PhenoAge (11 SNPs), and GrimAge (4 SNPs) on multiple cancers (i.e. breast, prostate, colorectal, ovarian and lung cancer). We obtained genome-wide association data for biological ageing from a meta-analysis (N = 34,710), and for cancer from the UK Biobank (N cases = 2671-13,879; N controls = 173,493-372,016), FinnGen (N cases = 719-8401; N controls = 74,685-174,006) and several international cancer genetic consortia (N cases = 11,348-122,977; N controls = 15,861-105,974). Main analyses were performed using multiplicative random effects inverse variance weighted (IVW) MR. Individual study estimates were pooled using fixed effect meta-analysis. Sensitivity analyses included MR-Egger, weighted median, weighted mode and Causal Analysis using Summary Effect Estimates (CAUSE) methods, which are robust to some of the assumptions of the IVW approach.
RESULTS: Meta-analysed IVW MR findings suggested that higher GrimAge acceleration increased the risk of colorectal cancer (OR = 1.12 per year increase in GrimAge acceleration, 95% CI 1.04-1.20, p = 0.002). The direction of the genetically predicted effects was consistent across main and sensitivity MR analyses. Among subtypes, the genetically predicted effect of GrimAge acceleration was greater for colon cancer (IVW OR = 1.15, 95% CI 1.09-1.21, p = 0.006), than rectal cancer (IVW OR = 1.05, 95% CI 0.97-1.13, p = 0.24). Results were less consistent for associations between other epigenetic clocks and cancers.
CONCLUSIONS: GrimAge acceleration may increase the risk of colorectal cancer. Findings for other clocks and cancers were inconsistent. Further work is required to investigate the potential mechanisms underlying the results.
FUNDING: FMB was supported by a Wellcome Trust PhD studentship in Molecular, Genetic and Lifecourse Epidemiology (224982/Z/22/Z which is part of grant 218495/Z/19/Z). KKT was supported by a Cancer Research UK (C18281/A29019) programme grant (the Integrative Cancer Epidemiology Programme) and by the Hellenic Republic's Operational Programme 'Competitiveness, Entrepreneurship & Innovation' (OΠΣ 5047228). PH was supported by Cancer Research UK (C18281/A29019). RMM was supported by the NIHR Biomedical Research Centre at University Hospitals Bristol and Weston NHS Foundation Trust and the University of Bristol and by a Cancer Research UK (C18281/A29019) programme grant (the Integrative Cancer Epidemiology Programme). RMM is a National Institute for Health Research Senior Investigator (NIHR202411). The views expressed are those of the author(s) and not necessarily those of the NIHR or the Department of Health and Social Care. GDS and CLR were supported by the Medical Research Council (MC_UU_00011/1 and MC_UU_00011/5, respectively) and by a Cancer Research UK (C18281/A29019) programme grant (the Integrative Cancer Epidemiology Programme). REM was supported by an Alzheimer's Society project grant (AS-PG-19b-010) and NIH grant (U01 AG-18-018, PI: Steve Horvath). RCR is a de Pass Vice Chancellor's Research Fellow at the University of Bristol.
OBJECTIVE: We performed an analysis of genetic variants associated with leukocyte telomere length to assess the relationship between telomere length and RCC risk using Mendelian randomization, an approach unaffected by biases from temporal variability and reverse causation that might have affected earlier investigations.
DESIGN, SETTING, AND PARTICIPANTS: Genotypes from nine telomere length-associated variants for 10 784 cases and 20 406 cancer-free controls from six genome-wide association studies (GWAS) of RCC were aggregated into a weighted genetic risk score (GRS) predictive of leukocyte telomere length.
OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Odds ratios (ORs) relating the GRS and RCC risk were computed in individual GWAS datasets and combined by meta-analysis.
RESULTS AND LIMITATIONS: Longer genetically inferred telomere length was associated with an increased risk of RCC (OR=2.07 per predicted kilobase increase, 95% confidence interval [CI]:=1.70-2.53, p<0.0001). As a sensitivity analysis, we excluded two telomere length variants in linkage disequilibrium (R2>0.5) with GWAS-identified RCC risk variants (rs10936599 and rs9420907) from the telomere length GRS; despite this exclusion, a statistically significant association between the GRS and RCC risk persisted (OR=1.73, 95% CI=1.36-2.21, p<0.0001). Exploratory analyses for individual histologic subtypes suggested comparable associations with the telomere length GRS for clear cell (N=5573, OR=1.93, 95% CI=1.50-2.49, p<0.0001), papillary (N=573, OR=1.96, 95% CI=1.01-3.81, p=0.046), and chromophobe RCC (N=203, OR=2.37, 95% CI=0.78-7.17, p=0.13).
CONCLUSIONS: Our investigation adds to the growing body of evidence indicating some aspect of longer telomere length is important for RCC risk.
PATIENT SUMMARY: Telomeres are segments of DNA at chromosome ends that maintain chromosomal stability. Our study investigated the relationship between genetic variants associated with telomere length and renal cell carcinoma risk. We found evidence suggesting individuals with inherited predisposition to longer telomere length are at increased risk of developing renal cell carcinoma.
OBJECTIVE: To identify biological pathways that contribute to risk for bipolar disorder (BP) using genes with consistent evidence for association in multiple genome-wide association studies (GWAS).
DATA SOURCES: Four independent data sets with individual genome-wide data available in July 2011 along with all data sets contributed to the Psychiatric Genomics Consortium Bipolar Group by May 2012. A prior meta-analysis was used as a source for brain gene expression data.
STUDY SELECTION: The 4 published GWAS were included in the initial sample. All independent BP data sets providing genome-wide data in the Psychiatric Genomics Consortium were included as a replication sample.
DATA EXTRACTION AND SYNTHESIS: We identified 966 genes that contained 2 or more variants associated with BP at P
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.
METHODS: The case-control portion of the study was conducted in nine UK centers with men ages 50-69 years who underwent prostate-specific antigen screening for prostate cancer within the Prostate Testing for Cancer and Treatment (ProtecT) trial. Two data sources were used to appraise causality: a genome-wide association study (GWAS) of metabolites in 24,925 participants and a GWAS of prostate cancer in 44,825 cases and 27,904 controls within the Association Group to Investigate Cancer Associated Alterations in the Genome (PRACTICAL) consortium.
RESULTS: Thirty-five metabolites were strongly associated with prostate cancer (P < 0.0014, multiple-testing threshold). These fell into four classes: (i) lipids and lipoprotein subclass characteristics (total cholesterol and ratios, cholesterol esters and ratios, free cholesterol and ratios, phospholipids and ratios, and triglyceride ratios); (ii) fatty acids and ratios; (iii) amino acids; (iv) and fluid balance. Fourteen top metabolites were proxied by genetic variables, but MR indicated these were not causal.
CONCLUSIONS: We identified 35 circulating metabolites associated with prostate cancer presence, but found no evidence of causality for those 14 testable with MR. Thus, the 14 MR-tested metabolites are unlikely to be mechanistically important in prostate cancer risk.
IMPACT: The metabolome provides a promising set of biomarkers that may aid prostate cancer classification.
METHODS: Using LDScore regression, we explored the genetic correlation between endometrial cancer and ovarian cancer. To identify loci associated with the risk of both cancers, we implemented a pipeline of statistical genetic analyses (i.e., inverse-variance meta-analysis, colocalization, and M-values) and performed analyses stratified by subtype. Candidate target genes were then prioritized using functional genomic data.
RESULTS: Genetic correlation analysis revealed significant genetic correlation between the two cancers (rG = 0.43, P = 2.66 × 10-5). We found seven loci associated with risk for both cancers (P Bonferroni < 2.4 × 10-9). In addition, four novel subgenome-wide regions at 7p22.2, 7q22.1, 9p12, and 11q13.3 were identified (P < 5 × 10-7). Promoter-associated HiChIP chromatin loops from immortalized endometrium and ovarian cell lines and expression quantitative trait loci data highlighted candidate target genes for further investigation.
CONCLUSIONS: Using cross-cancer GWAS meta-analysis, we have identified several joint endometrial and ovarian cancer risk loci and candidate target genes for future functional analysis.
IMPACT: Our research highlights the shared genetic relationship between endometrial cancer and ovarian cancer. Further studies in larger sample sets are required to confirm our findings.
Materials and Methods: All the variants' information was retrieved from the Ensembl genome database, and then different variation prediction analyses were performed. UTRScan was used to predict UTR variants while MaxEntScan was used to predict splice site variants. Meta-analysis by PredictSNP2 was used to predict sSNPs. Parallel prediction analyses by five different software packages including SIFT, PolyPhen-2, Mutation Assessor, I-Mutant2.0 and SNPs&GO were used to predict the effects of nsSNPs. The level of evolutionary conservation of deleterious nsSNPs was further analyzed using ConSurf server. Mutant protein structures of deleterious nsSNPs were modelled and refined using SPARKS-X and ModRefiner for structural comparison.
Results: A total of 56 deleterious variants were identified in this study, including 12 UTR variants, 18 splice site variants, eight sSNPs and 18 nsSNPs. Among these 56 deleterious variants, seven variants were also identified in the Alzheimer's Disease Sequencing Project (ADSP), Alzheimer's Disease Neuroimaging Initiative (ADNI) and Mount Sinai Brain Bank (MSBB) studies.
Discussion: The 56 deleterious variants were predicted to affect the regulation of gene expression, or have functional impacts on these three endocytosis genes and their gene products. The deleterious variants in these genes are expected to affect their cellular function in endocytosis and may be implicated in the pathogenesis of AD as well. The biological consequences of these deleterious variants and their potential impacts on the disease risks could be further validated experimentally and may be useful for gene-disease association study.
METHODS: A discovery cohort of Malaysian Chinese descent (NPC patients, n = 140; Healthy controls, n = 256) were genotyped using Illumina® HumanOmniExpress BeadChip. PennCNV and cnvPartition calling algorithms were applied for CNV calling. Taqman CNV assays and digital PCR were used to validate CNV calls and replicate candidate copy number variant region (CNVR) associations in a follow-up Malaysian Chinese (NPC cases, n = 465; and Healthy controls, n = 677) and Malay cohort (NPC cases, n = 114; Healthy controls, n = 124).
RESULTS: Six putative CNVRs overlapping GRM5, MICA/HCP5/HCG26, LILRB3/LILRA6, DPY19L2, RNase3/RNase2 and GOLPH3 genes were jointly identified by PennCNV and cnvPartition. CNVs overlapping GRM5 and MICA/HCP5/HCG26 were subjected to further validation by Taqman CNV assays and digital PCR. Combined analysis in Malaysian Chinese cohort revealed a strong association at CNVR on chromosome 11q14.3 (Pcombined = 1.54x10-5; odds ratio (OR) = 7.27; 95% CI = 2.96-17.88) overlapping GRM5 and a suggestive association at CNVR on chromosome 6p21.3 (Pcombined = 1.29x10-3; OR = 4.21; 95% CI = 1.75-10.11) overlapping MICA/HCP5/HCG26 genes.
CONCLUSION: Our results demonstrated the association of CNVs towards NPC susceptibility, implicating a possible role of CNVs in NPC development.