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.

METHODS: In total, DNA samples were obtained from 14,525 case subjects with invasive EOC and from 23,447 controls from 43 sites in the Ovarian Cancer Association Consortium (OCAC). Two hundred seventy nine SNPs, representing 131 genes, were genotyped using an Illumina Infinium iSelect BeadChip as part of the Collaborative Oncological Gene-environment Study (COGS). SNP analyses were conducted using unconditional logistic regression under a log-additive model, and the FDR q<0.2 was applied to adjust for multiple comparisons.

RESULTS: The most significant evidence of an association for all invasive cancers combined and for the serous subtype was observed for SNP rs17216603 in the iron transporter gene HEPH (invasive: OR = 0.85, P = 0.00026; serous: OR = 0.81, P = 0.00020); this SNP was also associated with the borderline/low malignant potential (LMP) tumors (P = 0.021). Other genes significantly associated with EOC histological subtypes (p<0.05) included the UGT1A (endometrioid), SLC25A45 (mucinous), SLC39A11 (low malignant potential), and SERPINA7 (clear cell carcinoma). In addition, 1785 SNPs in six genes (HEPH, MGST1, SERPINA, SLC25A45, SLC39A11 and UGT1A) were imputed from the 1000 Genomes Project and examined for association with INV EOC in white-European subjects. The most significant imputed SNP was rs117729793 in SLC39A11 (per allele, OR = 2.55, 95% CI = 1.5-4.35, p = 5.66x10-4).

METHODS: In a population of 15,596 epithelial ovarian cancer (EOC) cases and 23,236 controls, we measured genetic associations of 1,351 SNPs in Treg cell pathway genes with odds of ovarian cancer and tested pathway and gene-level associations, overall and by histotype, for the 25 genes, using the admixture likelihood (AML) method. The most significant single SNP associations were tested for correlation with expression levels in 44 ovarian cancer patients.

RESULTS: The most significant global associations for all genes in the pathway were seen in endometrioid ( p = 0.082) and clear cell ( p = 0.083), with the most significant gene level association seen with TGFBR2 ( p = 0.001) and clear cell EOC. Gene associations with histotypes at p < 0.05 included: IL12 ( p = 0.005 and p = 0.008, serous and high-grade serous, respectively), IL8RA ( p = 0.035, endometrioid and mucinous), LGALS1 ( p = 0.03, mucinous), STAT5B ( p = 0.022, clear cell), TGFBR1 ( p = 0.021 endometrioid) and TGFBR2 ( p = 0.017 and p = 0.025, endometrioid and mucinous, respectively).

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.