The full spectrum of SLC22 OCT1 mutations illuminates the bridge between drug transporter biophysics and pharmacogenomics

Mutations in transporters can impact an individual's response to drugs and cause many diseases. Few variants in transporters have been evaluated for their functional impact. Here, we combine saturation mutagenesis and multi-phenotypic screening to dissect the impact of 11,213 missense single-amino-acid deletions, and synonymous variants across the 554 residues of OCT1, a key liver xenobiotic transporter. By quantifying in parallel expression and substrate uptake, we find that most variants exert their primary effect on protein abundance, a phenotype not commonly measured alongside function. Using our mutagenesis results combined with structure prediction and molecular dynamic simulations, we develop accurate structure-function models of the entire transport cycle, providing biophysical characterization of all known and possible human OCT1 polymorphisms. This work provides a complete functional map of OCT1 variants along with a framework for integrating functional genomics, biophysical modeling, and human genetics to predict variant effects on disease and drug efficacy.

Association of immunoglobulin E levels with glioma risk and survival

Background

Previous epidemiological studies have reported an association of serum immunoglobulin E (IgE) levels with reduced glioma risk, but the association between IgE and glioma prognosis is not well characterized. This study aimed to examine how sex, tumor subtype, and IgE class modulate the association of serum IgE levels with glioma risk and survival.

Methods

We conducted a case-control study using participants from the University of California, San Francisco Adult Glioma Study (1997-2010). Serum IgE levels for total, respiratory and food allergy were measured in adults diagnosed with glioma (n=1,696) and cancer-free controls (n=1,135) matched based on age, sex, and race/ethnicity. Logistic regression was adjusted for patient demographics to assess the association between IgE levels and glioma risk. Multivariable Cox regression adjusted for patient-specific and tumor-specific factors compared survival between the elevated and normal IgE groups.

Results

Elevated total IgE was associated with reduced risk of IDH wildtype (OR=0.65, 95% CI: 0.54-0.78) and IDH mutant glioma (OR=0.65, 95% CI: 0.50-0.85). In multivariable Cox regression, elevated respiratory IgE was associated with improved survival for IDH wildtype glioma (HR=0.78, 95% CI: 0.67-0.91). The reduction in mortality risk was more pronounced in females (HR=0.71, 95% CI: 0.53-0.96) than in males (HR=0.80, 95% CI: 0.66-0.97), with improvements in median survival of 6.2 months (P<.001) and 1.6 months (P=0.003), respectively.

Conclusion

Elevated serum IgE was associated with improved prognosis for IDH wildtype glioma, with a more pronounced protective effect in females. These results suggest a possible sexual dimorphism and antitumor activity of IgE-mediated immune responses.

Improved multi-ancestry fine-mapping identifies cis -regulatory variants underlying molecular traits and disease risk

Multi-ancestry statistical fine-mapping of cis -molecular quantitative trait loci ( cis -molQTL) aims to improve the precision of distinguishing causal cis -molQTLs from tagging variants. However, existing approaches fail to reflect shared genetic architectures. To solve this limitation, we present the Sum of Shared Single Effects (SuShiE) model, which leverages LD heterogeneity to improve fine-mapping precision, infer cross-ancestry effect size correlations, and estimate ancestry-specific expression prediction weights. We apply SuShiE to mRNA expression measured in PBMCs (n=956) and LCLs (n=814) together with plasma protein levels (n=854) from individuals of diverse ancestries in the TOPMed MESA and GENOA studies. We find SuShiE fine-maps cis -molQTLs for 16 % more genes compared with baselines while prioritizing fewer variants with greater functional enrichment. SuShiE infers highly consistent cis -molQTL architectures across ancestries on average; however, we also find evidence of heterogeneity at genes with predicted loss-of-function intolerance, suggesting that environmental interactions may partially explain differences in cis -molQTL effect sizes across ancestries. Lastly, we leverage estimated cis -molQTL effect-sizes to perform individual-level TWAS and PWAS on six white blood cell-related traits in AOU Biobank individuals (n=86k), and identify 44 more genes compared with baselines, further highlighting its benefits in identifying genes relevant for complex disease risk. Overall, SuShiE provides new insights into the cis -genetic architecture of molecular traits.

A noncoding regulatory variant in IKZF1 increases acute lymphoblastic leukemia risk in Hispanic/Latino children

Hispanic/Latino children have the highest risk of acute lymphoblastic leukemia (ALL) in the US compared to other racial/ethnic groups, yet the basis of this remains incompletely understood. Through genetic fine-mapping analyses, we identified a new independent childhood ALL risk signal near IKZF1 in self-reported Hispanic/Latino individuals, but not in non-Hispanic White individuals, with an effect size of ∼1.44 (95% confidence interval = 1.33-1.55) and a risk allele frequency of ∼18% in Hispanic/Latino populations and <0.5% in European populations. This risk allele was positively associated with Indigenous American ancestry, showed evidence of selection in human history, and was associated with reduced IKZF1 expression. We identified a putative causal variant in a downstream enhancer that is most active in pro-B cells and interacts with the IKZF1 promoter. This variant disrupts IKZF1 autoregulation at this enhancer and results in reduced enhancer activity in B cell progenitors. Our study reveals a genetic basis for the increased ALL risk in Hispanic/Latino children.

Polygenic risk score for ulcerative colitis predicts immune checkpoint inhibitor-mediated colitis

Immune checkpoint inhibitor-mediated colitis (IMC) is a common adverse event of treatment with immune checkpoint inhibitors (ICI). We hypothesize that genetic susceptibility to Crohn's disease (CD) and ulcerative colitis (UC) predisposes to IMC. In this study, we first develop a polygenic risk scores for CD (PRSCD) and UC (PRSUC) in cancer-free individuals and then test these PRSs on IMC in a cohort of 1316 patients with ICI-treated non-small cell lung cancer and perform a replication in 873 ICI-treated pan-cancer patients. In a meta-analysis, the PRSUC predicts all-grade IMC (ORmeta=1.35 per standard deviation [SD], 95% CI = 1.12-1.64, P = 2×10-03) and severe IMC (ORmeta=1.49 per SD, 95% CI = 1.18-1.88, P = 9×10-04). PRSCD is not associated with IMC. Furthermore, PRSUC predicts severe IMC among patients treated with combination ICIs (ORmeta=2.20 per SD, 95% CI = 1.07-4.53, P = 0.03). Overall, PRSUC can identify patients receiving ICI at risk of developing IMC and may be useful to monitor patients and improve patient outcomes.

Genome-wide Polygenic Risk Scores Predict Risk of Glioma and Molecular Subtypes

Background

Polygenic risk scores (PRS) aggregate the contribution of many risk variants to provide a personalized genetic susceptibility profile. Since sample sizes of glioma genome-wide association studies (GWAS) remain modest, there is a need to find efficient ways of capturing genetic risk factors using available germline data.

Methods

We developed a novel PRS (PRS-CS) that uses continuous shrinkage priors to model the joint effects of over 1 million polymorphisms on disease risk and compared it to an approach (PRS-CT) that selects a limited set of independent variants that reach genome-wide significance (P<5×10-8). PRS models were trained using GWAS results stratified by histological (10,346 cases, 14,687 controls) and molecular subtype (2,632 cases, 2,445 controls), and validated in two independent cohorts.

Results

PRS-CS was consistently more predictive than PRS-CT across glioma subtypes with an average increase in explained variance (R2) of 21%. Improvements were particularly pronounced for glioblastoma tumors, with PRS-CS yielding larger effect sizes (odds ratio (OR)=1.93, P=2.0×10-54 vs. OR=1.83, P=9.4×10-50) and higher explained variance (R2=2.82% vs. R2=2.56%). Individuals in the 95th percentile of the PRS-CS distribution had a 3-fold higher lifetime absolute risk of IDH mutant (0.63%) and IDH wildtype (0.76%) glioma relative to individuals with average PRS. PRS-CS also showed high classification accuracy for IDH mutation status among cases (AUC=0.895).

Conclusions

Our novel genome-wide PRS may improve the identification of high-risk individuals and help distinguish between prognostic glioma subtypes, increasing the potential clinical utility of germline genetics in glioma patient management.

Principles and methods for transferring polygenic risk scores across global populations

Polygenic risk scores (PRSs) summarize the genetic predisposition of a complex human trait or disease and may become a valuable tool for advancing precision medicine. However, PRSs that are developed in populations of predominantly European genetic ancestries can increase health disparities due to poor predictive performance in individuals of diverse and complex genetic ancestries. We describe genetic and modifiable risk factors that limit the transferability of PRSs across populations and review the strengths and weaknesses of existing PRS construction methods for diverse ancestries. Developing PRSs that benefit global populations in research and clinical settings provides an opportunity for innovation and is essential for health equity.

Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

Clinical consequences of a genetic predisposition toward higher benign prostate-specific antigen levels

BACKGROUND

Prostate-specific antigen (PSA) levels are influenced by genetic variation unrelated to prostate cancer risk. Whether a genetic predisposition to a higher PSA level predisposes to a diagnostic work-up for prostate cancer is not known.

METHODS

Participants were 3110 men of African and European ancestries ages 45-70, without prostate cancer and with a baseline PSA < 4 ng/mL, undergoing routine clinical PSA screening. The exposure was a polygenic score (PGS) comprising 111 single nucleotide polymorphisms associated with PSA level, but not prostate cancer. We tested whether the PGS was associated with a: 1) PSA value > 4 ng/mL, 2) International Classification of Diseases (ICD) code for an elevated PSA, 3) encounter with a urologist, or 4) prostate biopsy. Multivariable Cox proportional hazards models were adjusted for age and genetic principal components. Analyses were stratified by age (45-59 years, and 60-70 years old). Association estimates are per standard deviation change in the PGS.

FINDINGS

The median age was 56.6 years, and 2118 (68%) participants were 45-59 years. The median (IQR) baseline PSA level was 1.0 (0.6-1.7) ng/mL. Among men ages 45-59, the PGS was associated with a PSA > 4 (hazard ratio [HR] = 1.35 [95% CI, 1.17-1.57], p = 4.5 × 10-5), an ICD code for elevated PSA (HR = 1.30 [1.12-1.52], p = 8.0 × 10-4), a urological evaluation (HR = 1.34 [1.14-1.57], p = 4.8 × 10-4), and undergoing a prostate biopsy (HR = 1.35 [1.11-1.64], p = 0.002). Among men ages 60-70, association effect sizes were smaller and not significant.

INTERPRETATION

A predisposition toward higher PSA levels was associated with clinical evaluations of an elevated PSA among men ages 45-59 years.

FUNDING

National Institutes of Health (NIH).

Genome-wide association study of prostate-specific antigen levels in 392,522 men identifies new loci and improves cross-ancestry prediction

We conducted a multi-ancestry genome-wide association study of prostate-specific antigen (PSA) levels in 296,754 men (211,342 European ancestry; 58,236 African ancestry; 23,546 Hispanic/Latino; 3,630 Asian ancestry; 96.5% of participants were from the Million Veteran Program). We identified 318 independent genome-wide significant (p≤5e-8) variants, 184 of which were novel. Most demonstrated evidence of replication in an independent cohort (n=95,768). Meta-analyzing discovery and replication (n=392,522) identified 447 variants, of which a further 111 were novel. Out-of-sample variance in PSA explained by our new polygenic risk score reached 16.9% (95% CI=16.1%-17.8%) in European ancestry, 9.5% (95% CI=7.0%-12.2%) in African ancestry, 18.6% (95% CI=15.8%-21.4%) in Hispanic/Latino, and 15.3% (95% CI=12.7%-18.1%) in Asian ancestry, and lower for higher age. Our study highlights how including proportionally more participants from underrepresented populations improves genetic prediction of PSA levels, with potential to personalize prostate cancer screening.

Genetic predisposition to altered blood cell homeostasis is associated with glioma risk and survival

Glioma is a highly fatal brain tumor comprised of molecular subtypes with distinct clinical trajectories. Observational studies have suggested that variability in immune response may play a role in glioma etiology. However, their findings have been inconsistent and susceptible to reverse causation due to treatment effects and the immunosuppressive nature of glioma. We applied genetic variants associated (p<5×10-8) with blood cell traits to a meta-analysis of 3418 glioma cases and 8156 controls. Genetically predicted increase in the platelet to lymphocyte ratio (PLR) was associated with an increased risk of glioma (odds ratio (OR)=1.25, p=0.005), especially in IDH-mutant (IDHmut OR=1.38, p=0.007) and IDHmut 1p/19q non-codeleted (IDHmut-noncodel OR=1.53, p=0.004) tumors. However, reduced glioma risk was observed for higher counts of lymphocytes (IDHmut-noncodel OR=0.70, p=0.004) and neutrophils (IDHmut OR=0.69, p=0.019; IDHmut-noncodel OR=0.60, p=0.009), which may reflect genetic predisposition to enhanced immune-surveillance. In contrast to susceptibility, there was no association with survival in IDHmut-noncodel; however, in IDHmut 1p/19q co-deleted tumors, we observed higher mortality with increasing genetically predicted counts of lymphocytes (hazard ratio (HR)=1.65, 95% CI: 1.24-2.20), neutrophils (HR=1.49, 1.13-1.97), and eosinophils (HR=1.59, 1.18-2.14). Polygenic scores for blood cell traits were also associated with tumor immune microenvironment features, with heterogeneity by IDH status observed for 17 signatures related to interferon signaling, PD-1 expression, and T-cell/Cytotoxic responses. In summary, we identified novel, immune-mediated susceptibility mechanisms for glioma with potential disease management implications.

Polygenic risk score for ulcerative colitis predicts immune checkpoint inhibitor-mediated colitis

Immune checkpoint inhibitors (ICIs) are a remarkable advancement in cancer therapeutics; however, a substantial proportion of patients develop severe immune-related adverse events (irAEs). Understanding and predicting irAEs is a key to advancing precision immuno-oncology. Immune checkpoint inhibitor-mediated colitis (IMC) is a significant complication from ICI and can have life-threatening consequences. Based on clinical presentation, IMC mimics inflammatory bowel disease, however the link is poorly understood. We hypothesized that genetic susceptibility to Crohn's disease (CD) and ulcerative colitis (UC) may predispose to IMC. We developed and validated polygenic risk scores for CD (PRSCD) and UC (PRSUC) in cancer-free individuals and assessed the role of each of these PRSs on IMC in a cohort of 1,316 patients with non-small cell lung cancer who received ICIs. Prevalence of all-grade IMC in our cohort was 4% (55 cases), and for severe IMC, 2.5% (32 cases). The PRSUC predicted the development of all-grade IMC (HR=1.34 per standard deviation [SD], 95% CI=1.02-1.76, P=0.04) and severe IMC (HR=1.62 per SD, 95% CI=1.12-2.35, P=0.01). PRSCD was not associated with IMC or severe IMC. The association between PRSUC and IMC (all-grade and severe) was consistent in an independent pan-cancer cohort of patients treated with ICIs. Furthermore, PRSUC predicted severe IMC among patients treated with combination ICIs (OR = 2.20 per SD, 95% CI = 1.07-4.53, P=0.03). This is the first study to demonstrate the potential clinical utility of a PRS for ulcerative colitis in identifying patients receiving ICI at high risk of developing IMC, where risk reduction and close monitoring strategies could help improve overall patient outcomes.

The full spectrum of OCT1 (SLC22A1) mutations bridges transporter biophysics to drug pharmacogenomics

Membrane transporters play a fundamental role in the tissue distribution of endogenous compounds and xenobiotics and are major determinants of efficacy and side effects profiles. Polymorphisms within these drug transporters result in inter-individual variation in drug response, with some patients not responding to the recommended dosage of drug whereas others experience catastrophic side effects. For example, variants within the major hepatic Human organic cation transporter OCT1 (SLC22A1) can change endogenous organic cations and many prescription drug levels. To understand how variants mechanistically impact drug uptake, we systematically study how all known and possible single missense and single amino acid deletion variants impact expression and substrate uptake of OCT1. We find that human variants primarily disrupt function via folding rather than substrate uptake. Our study revealed that the major determinants of folding reside in the first 300 amino acids, including the first 6 transmembrane domains and the extracellular domain (ECD) with a stabilizing and highly conserved stabilizing helical motif making key interactions between the ECD and transmembrane domains. Using the functional data combined with computational approaches, we determine and validate a structure-function model of OCT1s conformational ensemble without experimental structures. Using this model and molecular dynamic simulations of key mutants, we determine biophysical mechanisms for how specific human variants alter transport phenotypes. We identify differences in frequencies of reduced function alleles across populations with East Asians vs European populations having the lowest and highest frequency of reduced function variants, respectively. Mining human population databases reveals that reduced function alleles of OCT1 identified in this study associate significantly with high LDL cholesterol levels. Our general approach broadly applied could transform the landscape of precision medicine by producing a mechanistic basis for understanding the effects of human mutations on disease and drug response.

Genetically adjusted PSA levels for prostate cancer screening

Prostate-specific antigen (PSA) screening for prostate cancer remains controversial because it increases overdiagnosis and overtreatment of clinically insignificant tumors. Accounting for genetic determinants of constitutive, non-cancer-related PSA variation has potential to improve screening utility. In this study, we discovered 128 genome-wide significant associations (P < 5 × 10-8) in a multi-ancestry meta-analysis of 95,768 men and developed a PSA polygenic score (PGSPSA) that explains 9.61% of constitutive PSA variation. We found that, in men of European ancestry, using PGS-adjusted PSA would avoid up to 31% of negative prostate biopsies but also result in 12% fewer biopsies in patients with prostate cancer, mostly with Gleason score <7 tumors. Genetically adjusted PSA was more predictive of aggressive prostate cancer (odds ratio (OR) = 3.44, P = 6.2 × 10-14, area under the curve (AUC) = 0.755) than unadjusted PSA (OR = 3.31, P = 1.1 × 10-12, AUC = 0.738) in 106 cases and 23,667 controls. Compared to a prostate cancer PGS alone (AUC = 0.712), including genetically adjusted PSA improved detection of aggressive disease (AUC = 0.786, P = 7.2 × 10-4). Our findings highlight the potential utility of incorporating PGS for personalized biomarkers in prostate cancer screening.

Gene expression in African Americans, Puerto Ricans and Mexican Americans reveals ancestry-specific patterns of genetic architecture

We explored ancestry-related differences in the genetic architecture of whole-blood gene expression using whole-genome and RNA sequencing data from 2,733 African Americans, Puerto Ricans and Mexican Americans. We found that heritability of gene expression significantly increased with greater proportions of African genetic ancestry and decreased with higher proportions of Indigenous American ancestry, reflecting the relationship between heterozygosity and genetic variance. Among heritable protein-coding genes, the prevalence of ancestry-specific expression quantitative trait loci (anc-eQTLs) was 30% in African ancestry and 8% for Indigenous American ancestry segments. Most anc-eQTLs (89%) were driven by population differences in allele frequency. Transcriptome-wide association analyses of multi-ancestry summary statistics for 28 traits identified 79% more gene-trait associations using transcriptome prediction models trained in our admixed population than models trained using data from the Genotype-Tissue Expression project. Our study highlights the importance of measuring gene expression across large and ancestrally diverse populations for enabling new discoveries and reducing disparities.

Development and testing of a polygenic risk score for breast cancer aggressiveness

Aggressive breast cancers portend a poor prognosis, but current polygenic risk scores (PRSs) for breast cancer do not reliably predict aggressive cancers. Aggressiveness can be effectively recapitulated using tumor gene expression profiling. Thus, we sought to develop a PRS for the risk of recurrence score weighted on proliferation (ROR-P), an established prognostic signature. Using 2363 breast cancers with tumor gene expression data and single nucleotide polymorphism (SNP) genotypes, we examined the associations between ROR-P and known breast cancer susceptibility SNPs using linear regression models. We constructed PRSs based on varying p-value thresholds and selected the optimal PRS based on model r2 in 5-fold cross-validation. We then used Cox proportional hazards regression to test the ROR-P PRS's association with breast cancer-specific survival in two independent cohorts totaling 10,196 breast cancers and 785 events. In meta-analysis of these cohorts, higher ROR-P PRS was associated with worse survival, HR per SD = 1.13 (95% CI 1.06-1.21, p = 4.0 × 10-4). The ROR-P PRS had a similar magnitude of effect on survival as a comparator PRS for estrogen receptor (ER)-negative versus positive cancer risk (PRSER-/ER+). Furthermore, its effect was minimally attenuated when adjusted for PRSER-/ER+, suggesting that the ROR-P PRS provides additional prognostic information beyond ER status. In summary, we used integrated analysis of germline SNP and tumor gene expression data to construct a PRS associated with aggressive tumor biology and worse survival. These findings could potentially enhance risk stratification for breast cancer screening and prevention.

Transcriptome-Wide Association Analysis Identifies Novel Candidate Susceptibility Genes for Prostate-Specific Antigen Levels in Men Without Prostate Cancer

Deciphering the genetic basis of prostate-specific antigen (PSA) levels may improve their utility to screen for prostate cancer (PCa). We thus conducted a transcriptome-wide association study (TWAS) of PSA levels using genome-wide summary statistics from 95,768 PCa-free men, the MetaXcan framework, and gene prediction models trained in Genotype-Tissue Expression (GTEx) project data. Tissue-specific analyses identified 41 statistically significant (p < 0.05/12,192 = 4.10e-6) associations in whole blood and 39 statistically significant (p < 0.05/13,844 = 3.61e-6) associations in prostate tissue, with 18 genes associated in both tissues. Cross-tissue analyses that combined associations across 45 tissues identified 155 genes that were statistically significantly (p < 0.05/22,249 = 2.25e-6) associated with PSA levels. Based on conditional analyses that assessed whether TWAS associations were attributable to a lead GWAS variant, we found 20 novel genes (11 single-tissue, 9 cross-tissue) that were associated with PSA levels in the TWAS. Of these novel genes, five showed evidence of colocalization (colocalization probability > 0.5): EXOSC9, CCNA2, HIST1H2BN, RP11-182L21.6, and RP11-327J17.2. Six of the 20 novel genes are not known to impact PCa risk. These findings yield new hypotheses for genetic factors underlying PSA levels that should be further explored toward improving our understanding of PSA biology.

Common genetic variations in telomere length genes and lung cancer: a Mendelian randomisation study and its novel application in lung tumour transcriptome

Background

Genome-wide association studies (GWASs) have identified genetic susceptibility variants for both leukocyte telomere length (LTL) and lung cancer susceptibility. Our study aims to explore the shared genetic basis between these traits and investigate their impact on somatic environment of lung tumours.

Methods

We performed genetic correlation, Mendelian randomisation (MR), and colocalisation analyses using the largest available GWASs summary statistics of LTL (N=464,716) and lung cancer (N=29,239 cases and 56,450 controls). Principal components analysis based on RNA-sequencing data was used to summarise gene expression profile in lung adenocarcinoma cases from TCGA (N=343).

Results

Although there was no genome-wide genetic correlation between LTL and lung cancer risk, longer LTL conferred an increased risk of lung cancer regardless of smoking status in the MR analyses, particularly for lung adenocarcinoma. Of the 144 LTL genetic instruments, 12 colocalised with lung adenocarcinoma risk and revealed novel susceptibility loci, including MPHOSPH6, PRPF6, and POLI. The polygenic risk score for LTL was associated with a specific gene expression profile (PC2) in lung adenocarcinoma tumours. The aspect of PC2 associated with longer LTL was also associated with being female, never smokers, and earlier tumour stages. PC2 was strongly associated with cell proliferation score and genomic features related to genome stability, including copy number changes and telomerase activity.

Conclusions

This study identified an association between longer genetically predicted LTL and lung cancer and sheds light on the potential molecular mechanisms related to LTL in lung adenocarcinomas.

Funding

Institut National du Cancer (GeniLuc2017-1-TABAC-03-CIRC-1-TABAC17-022), INTEGRAL/NIH (5U19CA203654-03), CRUK (C18281/A29019), and Agence Nationale pour la Recherche (ANR-10-INBS-09).

A Polygenic Risk Score for Prostate Cancer Risk Prediction

This retrospective cohort study compares 2 risk calculator systems that compute the probabilities of finding high-grade or any cancer on biopsy results in men undergoing a first prostate biopsy.

Abstract PR008: Development and testing of a polygenic risk score for breast cancer. Aggressiveness

Background: Aggressive breast cancers have increased proliferation or metastatic potential and portend a poor prognosis. The ability to identify women at elevated risk of aggressive cancers could have major implications for screening and prevention, yet there are no available tools for predicting aggressive cancer risk. We sought to construct a polygenic risk score (PRS) for aggressive breast cancers by leveraging the associations of single nucleotide polymorphisms (SNPs) with tumor gene expression. We used as our measure of aggressiveness the risk of recurrence score weighted on proliferation (ROR-P), a validated tumor prognostic signature. We hypothesized that known breast cancer susceptibility SNPs would have differential associations with ROR-P, which could then be used to construct a PRS for ROR-P. Methods: We developed our PRS in a case-only analysis of 3 studies containing SNP genotypes and tumor gene expression: The Cancer Genome Atlas, METABRIC, and the I-SPY 2 TRIAL (total n=2,363). We used linear regression models to evaluate individual SNP associations with ROR-P, adjusted for genetic ancestry and study. We then constructed PRS using varying p-value thresholds and used cross-validation to identify the PRS with highest model r2. To assess whether the ROR-P PRS was associated with poor prognosis, we performed survival analysis in two longitudinal cohorts of breast cancer patients: the UK Biobank (women with incident invasive cancers only) and the Pathways Study. These studies included 10,196 total cases with 785 deaths. We built Cox proportional hazards models to evaluate the association between the ROR-P PRS (adjusted for genetic ancestry) and breast cancer-specific survival (BCSS) in both studies. We then performed meta-analysis of the Cox model results. We also constructed joint models containing the ROR-P PRS and a PRS representing the case-case risk of ER-negative vs. ER-positive cancer, PRSER-/ER+. Results: We tested the associations between 226 breast cancer susceptibility SNPs and ROR-P. The best-performing PRS contained 76 SNPs and had a cross-validated r2 of 0.051. In the UK Biobank and Pathways Study, higher ROR-P PRS was associated with worse BCSS, with nearly identical effects observed in each study, HR per standard deviation of 1.13 (95% CI 1.05-1.21, p=9.0x10-4) in meta-analysis. The ROR-P PRS’s effect was minimally attenuated when adjusted for PRSER-/ER+, suggesting that the ROR-P PRS was providing additional prognostic information beyond ER status. Conclusions: We used breast cancer susceptibility SNPs to construct a PRS for ROR-P, a prognostic signature recapitulating aggressiveness, and found the ROR-P PRS to be associated with worse BCSS. Our findings represent an improvement on current PRS for overall breast cancer risk, which preferentially predict cancers with favorable prognosis. Given that aggressive cancers are more likely to present as advanced cancers even among women undergoing routine screening, our findings could potentially identify women who may benefit from more intensive screening.

Citation Format: Yiwey Shieh, Jacquelyn Roger, Christina Yau, Denise Wolf, Gillian Hirst, Lamorna Swigart, Scott Huntsman, Donglei Hu, Jovia Nierenberg, Pooja Middha, Rachel Heise, Linda Kachuri, Qianqian Zhu, Song Yao, Christine Ambrosone, Marilyn Kwan, Bette Caan, John Witte, Lawrence Kushi, Laura van ’T. Veer, Laura Esserman, Elad Ziv. Development and testing of a polygenic risk score for breast cancer. Aggressiveness. [abstract]. In: Proceedings of the AACR Special Conference: Precision Prevention, Early Detection, and Interception of Cancer; 2022 Nov 17-19; Austin, TX. Philadelphia (PA): AACR; Can Prev Res 2023;16(1 Suppl): Abstract nr PR008.

Inverted genomic regions between reference genome builds in humans impact imputation accuracy and decrease the power of association testing

Over the last two decades, the human reference genome has undergone multiple updates as we complete a linear representation of our genome. Two versions of human references are currently used in the biomedical literature, GRCh37/hg19 and GRCh38. Conversions between these versions are critical for quality control, imputation, and association analysis. In the present study, we show that single-nucleotide variants (SNVs) in regions inverted between different builds of the reference genome are often mishandled bioinformatically. Depending on the array type, SNVs are found in approximately 2-5 Mb of the genome that are inverted between reference builds. Coordinate conversions of these variants are mishandled by both the TOPMed imputation server as well as routine in-house quality control pipelines, leading to underrecognized downstream analytical consequences. Specifically, we observe that undetected allelic conversion errors for palindromic (i.e., A/T or C/G) variants in these inverted regions would destabilize the local haplotype structure, leading to loss of imputation accuracy and power in association analyses. Though only a small proportion of the genome is affected, these regions include important disease susceptibility variants that would be affected. For example, the p value of a known locus associated with prostate cancer on chromosome 10 (chr10) would drop from 2.86 × 10^-7 to 0.0011 in a case-control analysis of 20,286 Africans and African Americans (10,643 cases and 9,643 controls). We devise a straight-forward heuristic based on the popular tool, liftOver, that can easily detect and correct these variants in the inverted regions between genome builds to locally improve imputation accuracy.

Assessment of genetic susceptibility to multiple primary cancers through whole-exome sequencing in two large multi-ancestry studies

BACKGROUND

Up to one of every six individuals diagnosed with one cancer will be diagnosed with a second primary cancer in their lifetime. Genetic factors contributing to the development of multiple primary cancers, beyond known cancer syndromes, have been underexplored.

METHODS

To characterize genetic susceptibility to multiple cancers, we conducted a pan-cancer, whole-exome sequencing study of individuals drawn from two large multi-ancestry populations (6429 cases, 165,853 controls). We created two groupings of individuals diagnosed with multiple primary cancers: (1) an overall combined set with at least two cancers across any of 36 organ sites and (2) cancer-specific sets defined by an index cancer at one of 16 organ sites with at least 50 cases from each study population. We then investigated whether variants identified from exome sequencing were associated with these sets of multiple cancer cases in comparison to individuals with one and, separately, no cancers.

RESULTS

We identified 22 variant-phenotype associations, 10 of which have not been previously discovered and were significantly overrepresented among individuals with multiple cancers, compared to those with a single cancer.

CONCLUSIONS

Overall, we describe variants and genes that may play a fundamental role in the development of multiple primary cancers and improve our understanding of shared mechanisms underlying carcinogenesis.

Genetic Analysis of Lung Cancer and the Germline Impact on Somatic Mutation Burden

BACKGROUND

Germline genetic variation contributes to lung cancer (LC) susceptibility. Previous genome-wide association studies (GWAS) have implicated susceptibility loci involved in smoking behaviors and DNA repair genes, but further work is required to identify susceptibility variants.

METHODS

To identify LC susceptibility loci, a family history-based genome-wide association by proxy (GWAx) of LC (48 843 European proxy LC patients, 195 387 controls) was combined with a previous LC GWAS (29 266 patients, 56 450 controls) by meta-analysis. Colocalization was used to explore candidate genes and overlap with existing traits at discovered susceptibility loci. Polygenic risk scores (PRS) were tested within an independent validation cohort (1 666 LC patients vs 6 664 controls) using variants selected from the LC susceptibility loci and a novel selection approach using published GWAS summary statistics. Finally, the effects of the LC PRS on somatic mutational burden were explored in patients whose tumor resections have been profiled by exome (n = 685) and genome sequencing (n = 61). Statistical tests were 2-sided.

RESULTS

The GWAx-GWAS meta-analysis identified 8 novel LC loci. Colocalization implicated DNA repair genes (CHEK1), metabolic genes (CYP1A1), and smoking propensity genes (CHRNA4 and CHRNB2). PRS analysis demonstrated that these variants, as well as subgenome-wide significant variants related to expression quantitative trait loci and/or smoking propensity, assisted in LC genetic risk prediction (odds ratio = 1.37, 95% confidence interval = 1.29 to 1.45; P < .001). Patients with higher genetic PRS loads of smoking-related variants tended to have higher mutation burdens in their lung tumors.

CONCLUSIONS

This study has expanded the number of LC susceptibility loci and provided insights into the molecular mechanisms by which these susceptibility variants contribute to LC development.

Cross-ancestry genome-wide meta-analysis of 61,047 cases and 947,237 controls identifies new susceptibility loci contributing to lung cancer

To identify new susceptibility loci to lung cancer among diverse populations, we performed cross-ancestry genome-wide association studies in European, East Asian and African populations and discovered five loci that have not been previously reported. We replicated 26 signals and identified 10 new lead associations from previously reported loci. Rare-variant associations tended to be specific to populations, but even common-variant associations influencing smoking behavior, such as those with CHRNA5 and CYP2A6, showed population specificity. Fine-mapping and expression quantitative trait locus colocalization nominated several candidate variants and susceptibility genes such as IRF4 and FUBP1. DNA damage assays of prioritized genes in lung fibroblasts indicated that a subset of these genes, including the pleiotropic gene IRF4, potentially exert effects by promoting endogenous DNA damage.

Abstract 1441: Genetic determinants of PSA levels improve prostate cancer screening

Prostate-specific antigen (PSA) screening for prostate cancer (PCa) remains controversial due to poor sensitivity and specificity that lead to overdiagnosis and overtreatment. The aim of our study is to characterize genetic determinants of PSA levels in cancer-free men in order to personalize PCa screening. We hypothesize that test accuracy may be improved by accounting for PSA variation that is due to genetic factors and does not reflect PCa.

We conducted the largest ever genome-wide association study (GWAS) of PSA in men without PCa (N=95,768; 85,924 predominantly European ancestry) using data from the UK Biobank, BioVU, PLCO, and Kaiser Permanente cohorts. Our GWAS discovered 129 PSA-associated variants (P&lt;5×10-8), 82 of which were novel. A polygenic score (PGSPSA) comprised of these 129 variants was successfully validated in two cancer prevention trials: PCPT (n=5737) and SELECT (n=22,247). PGSPSA explained 7.3% (p=7.0×10-98) of variation in baseline PSA in PCPT and 8.7% (p=7.0×10-476) in SELECT. Importantly, PGSPSA was not associated with PCa status in PCPT (OR=0.98, p=0.71) or SELECT (OR=1.04, p=0.98), which confirms that it reflects benign PSA variation.

Potential clinical utility of PSA genetic adjustment was explored by examining reclassification at thresholds used for biopsy referrals in a real-world setting at Kaiser Permanente. We estimated that correction using PGSPSA would have avoided 21.2% of negative biopsies in non-cases. Reclassification below the biopsy referral threshold was also more common in cases, particularly with low-grade disease with Gleason score &lt;7 (7.3% below vs. 2.6% above). Overall, genetic correction of PSA appeared to improve the accuracy of referral decisions, with a Net Reclassification Index of 0.148.

Next, we evaluated genetically adjusted PSA in the context of detection of aggressive PCa, defined as Gleason score ≥7, PSA ≥10 ng/mL, T3-T4 stage, and/or distant or nodal metastases. Genetically adjusted baseline PSA was more robustly associated with aggressive PCa than observed PSA and yielded a higher area under the curve (AUC) in PCPT (OR=3.03, p=3.5×10-7; AUC: 0.72 vs. 0.68) and SELECT (OR=3.37, p=3.5×10-11; AUC: 0.78 vs. 0.74) when added to a baseline model with age and trial arm. Furthermore, genetically adjusted PSA provides complementary information to PCa risk variants. In PCPT, a logistic regression model that included genetically corrected PSA and the 269-variant PGSPCa achieved a significantly higher AUC than PGSPCa-269 alone for aggressive PCa (AUC: 0.73 vs. 0.65, p=3.3×10-4) and overall PCa (AUC=0.69 vs. 0.66, p=3.3×10-6).

Our work provides evidence that accounting for genetic determinants of PSA has the potential to reduce unnecessary testing and overdiagnosis of low-risk PCa, as well as increase detection of aggressive disease. Larger and more diverse study populations are required to fully characterize the genetic basis of PSA variation and optimize its clinical utility.

Citation Format: Linda Kachuri, Rebecca E. Graff, Sonja I. Berndt, Mitchell Machiela, Neal D. Freedman, Stephen J. Chanock, John P. Shelley, Kerry Schaffer, Jonathan D. Mosley, Phyllis J. Goodman, Cathee Till, Ian Thompson, Robert J. Klein, Stephen K. Van Den Eeden, Thomas J. Hoffmann, John S. Witte. Genetic determinants of PSA levels improve prostate cancer screening [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1441.

Abstract 1446: Genetic risk factors for the development of multiple primary cancers

Introduction: Individuals with one cancer are at greater risk of new primary cancers than the general population. While several hereditary cancer syndromes are known, genetic risk factors for multiple primary cancers in an individual are not well understood. Identification of susceptibility variants to multiple primary tumors could enhance screening for subsequent cancers among those at highest risk.

Methods: We conducted a pan-cancer genome-wide association study (GWAS) of multiple primary cancers among participants from 2 prospective cohorts: Kaiser Permanente and the UK Biobank. The primary GWAS within cohorts used logistic regression to estimate associations for diagnosis with ≥2 invasive or in situ primary cancers other than non-melanoma skin (N=11,773, 8,928 invasive only) compared to cancer-free controls (N=420,101). Case-case analyses were conducted to distinguish associations with multiple cancers from single-cancer (N=90,576) susceptibility signals. Regression models were adjusted for age, sex, first 10 genetic ancestry principal components, and array. Cohort-specific GWAS results were meta-analyzed. We highlight genome-wide significant (p&lt;5×10-8) results with consistent effect direction across the 2 studies.

Results: We identified 8 variants associated with multiple primary cancers.

Discussion: To our knowledge, rs192703567 has no previous cancer associations. Three identified variants are in known cancer predisposition genes (rs2293607 in TERC, rs6983267 in CASC8, rs35850753 in TP53). Three variants (rs34379047, rs612611, and rs9419958) are previously associated with multiple cancers, 2 (rs2293607 and rs6983267) with cancer pleiotropy, and 2 (rs283732 and rs35850753) with individual cancers. Most variants from the cancer-free control analyses had consistent effects in the single-cancer case-case analyses, suggesting pleiotropic mechanisms. Our preliminary findings offer insight into genetic risk factors associated with developing multiple primary cancers.

Multiple vs. none Multiple invasive vs. none Multiple vs. single Multiple invasive vs. single Chr Position rsID A1 A2 Gene OR P OR P OR P OR P 3 169482335 rs2293607 T C TERC 1.11 1.0×10-9 1.10 2.2×10-7 1.06 2.0×10-3 1.05 0.01 8 128281644 rs283732 C T Intergenic 1.09 2.3×10-7 1.11 9.2×10-9 1.06 3.8×10-4 1.09 3.9×10-5 8 128413305 rs6983267 G T CASC8 POU5F1B CCAT2 PCAT1 1.08 4.1×10-7 1.09 2.9×10-8 1.00 0.83 1.02 0.34 10 105644473 rs34379047 T A OBFC1 1.16 8.5×10-12 1.19 6.0×10-11 1.19 2.3×10-4 1.10 4.0×10-4 10 105675946 rs9419958 T C STN1 1.14 1.8×10-11 1.16 4.5×10-11 1.08 3.1×10-4 1.09 2.8×10-4 11 69307463 rs612611 G A Intergenic 1.11 3.7×10-8 1.10 5.3×10-6 1.06 9.6×10-4 1.05 0.02 17 7578671 rs35850753 T C TP53 1.27 4.9×10-7 1.34 1.6×10-8 1.16 2.8×10-3 1.24 8.8×10-5 22 40738280 rs192703567 C T Intergenic 1.37 5.2×10-8 1.37 1.8×10-6 1.43 4.9×10-9 1.43 2.8×10-7 Chr=Chromosome; A1=Effect allele; A2=Other allele; OR=Odds ratio.

Citation Format: Jovia L. Nierenberg, Linda Kachuri, Taylor B. Cavazos, Rebecca E. Graff, Thomas J. Hoffmann, Jie Zhang, Stacey Alexeeff, Laurel Habel, Douglas Corley, Stephen Van Den Eeden, Elad Ziv, Lori C. Sakoda, John S. Witte. Genetic risk factors for the development of multiple primary cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1446.

Abstract 3631: Transcriptome-wide association study identifies novel genes associated with breast cancer susceptibility in Latinas

Background: Genetic susceptibility to breast cancer has been studied extensively in European ancestry populations, but few studies have addressed genetic susceptibility in non-European women. Latinas are a genetically diverse group with contributions from European, African, and Indigenous American ancestries. Genome-wide association studies (GWAS) have identified unique variants in this population, particularly at the 6q25 locus. We conducted a transcriptome-wide association study (TWAS) to identify novel genes associated with risk of breast cancer in Latinas.

Methods: We used individual level GWAS data from 2,396 Latina cases and 6,505 Latina controls from the studies in Northern California (San Francisco Bay Area Breast Cancer Study, Northern California Breast Cancer Family Registry and Kaiser Permanente Genetic Epidemiology Research on Aging Cohort), Southern California (Multi-ethnic Cohort) and Mexico (CAMA study). We analyzed the association between genetically predicted whole blood (WB) gene expression and breast cancer risk using newly developed TWAS models based on 784 Mexican American individuals. We also conducted parallel analyses using breast mammary tissue (BT) TWAS models from GTEx v8. All analyses were adjusted for age, ancestry, and study. Associations with false discovery rate (FDR) probability &lt;0.05 were considered statistically significant.

Results: At FDR&lt;0.05, we identified 20 genes from BT and 39 genes from WB. Seven of the genes were significantly associated in both the WB models and GTEx BT models. Increased expression of MIB2 (pFDR = 4.74x10-17 (WB) and 1.22x10-4 (BT)), NBPF26 (pFDR = 1.30x10-7 (WB) and 7.43x10-8 (BT)), SLC35E2B (pFDR = 1.12x10-4 (WB) and 5.47x10-5 (BT)), and FAM30A (pFDR = 1.18x10-10 (WB) and 9.27x10-3 (BT)) was associated with increased risk of breast cancer risk, whereas increased expression of SLC35E2A (pFDR = 8.17x10-6 (WB) and 1.55x10-3 (BT)) and HCP5B (pFDR = 1.84x10-3 (WB) and 2.19x10-3 (BT)) was associated with decreased breast cancer risk. Increased expression of PDGFA was associated with increased risk (pFDR = 1.30x10-7) in GTEx BT reference models but decreased risk (pFDR = 2.34x10-10) in the ancestry-specific WB model.

Conclusion: Our study is the first TWAS investigating the relationship between genetically predicted gene expression and breast cancer risk in Latinas. By leveraging gene expression prediction models that capture eQTLs that are more common in populations with Indigenous American ancestry, we have identified some novel genes associated with breast cancer risk in Latinas. Of these, MIB2 is a strong candidate for a mechanistic role in breast carcinogenesis in Latinas. MIB2 is involved in Notch signaling which plays an important role in breast carcinogenesis via its mismatched receptor-ligand interaction. Our study highlights the role of ancestry-based prediction models in TWAS analyses.

Citation Format: Pooja Middha Kapoor, Angel C. Mak, Linda Kachuri, Donglei Hu, Scott Huntsman, Lawrence H. Kushi, Christopher Haiman, Esther M. John, Gabriela Torres-Mejia, Esteban G. Burchard, Susan L. Neuhausen, Laura Fejerman, Elad Ziv. Transcriptome-wide association study identifies novel genes associated with breast cancer susceptibility in Latinas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3631.

The immunogenetics of viral antigen response is associated with subtype-specific glioma risk and survival

Glioma is a highly fatal cancer with prognostically significant molecular subtypes and few known risk factors. Multiple studies have implicated infections in glioma susceptibility, but evidence remains inconsistent. Genetic variants in the human leukocyte antigen (HLA) region modulate host response to infection and have been linked to glioma risk. In this study, we leveraged genetic predictors of antibody response to 12 viral antigens to investigate the relationship with glioma risk and survival. Genetic reactivity scores (GRSs) for each antigen were derived from genome-wide-significant (p < 5 × 10-8) variants associated with immunoglobulin G antibody response in the UK Biobank cohort. We conducted parallel analyses of glioma risk and survival for each GRS and HLA alleles imputed at two-field resolution by using data from 3,418 glioma-affected individuals subtyped by somatic mutations and 8,156 controls. Genetic reactivity scores to Epstein-Barr virus (EBV) ZEBRA and EBNA antigens and Merkel cell polyomavirus (MCV) VP1 antigen were associated with glioma risk and survival (Bonferroni-corrected p < 0.01). GRSZEBRA and GRSMCV were associated in opposite directions with risk of IDH wild-type gliomas (ORZEBRA = 0.91, p = 0.0099/ORMCV = 1.11, p = 0.0054). GRSEBNA was associated with both increased risk for IDH mutated gliomas (OR = 1.09, p = 0.040) and improved survival (HR = 0.86, p = 0.010). HLA-DQA1∗03:01 was significantly associated with decreased risk of glioma overall (OR = 0.85, p = 3.96 × 10-4) after multiple testing adjustment. This systematic investigation of the role of genetic determinants of viral antigen reactivity in glioma risk and survival provides insight into complex immunogenomic mechanisms of glioma pathogenesis. These results may inform applications of antiviral-based therapies in glioma treatment.

Investigating the effect of sexual behaviour on oropharyngeal cancer risk: a methodological assessment of Mendelian randomization

BACKGROUND

Human papilloma virus infection is known to influence oropharyngeal cancer (OPC) risk, likely via sexual transmission. However, sexual behaviour has been correlated with other risk factors including smoking and alcohol, meaning independent effects are difficult to establish. We aimed to evaluate the causal effect of sexual behaviour on the risk of OPC using Mendelian randomization (MR).

METHODS

Genetic variants robustly associated with age at first sex (AFS) and the number of sexual partners (NSP) were used to perform both univariable and multivariable MR analyses with summary data on 2641 OPC cases and 6585 controls, obtained from the largest available genome-wide association studies (GWAS). Given the potential for genetic pleiotropy, we performed a number of sensitivity analyses: (i) MR methods to account for horizontal pleiotropy, (ii) MR of sexual behaviours on positive (cervical cancer and seropositivity for Chlamydia trachomatis) and negative control outcomes (lung and oral cancer), (iii) Causal Analysis Using Summary Effect estimates (CAUSE), to account for correlated and uncorrelated horizontal pleiotropic effects, (iv) multivariable MR analysis to account for the effects of smoking, alcohol, risk tolerance and educational attainment.

RESULTS

In univariable MR, we found evidence supportive of an effect of both later AFS (IVW OR = 0.4, 95%CI (0.3, 0.7), per standard deviation (SD), p = < 0.001) and increasing NSP (IVW OR = 2.2, 95%CI (1.3, 3.8) per SD, p = < 0.001) on OPC risk. These effects were largely robust to sensitivity analyses accounting for horizontal pleiotropy. However, negative control analysis suggested potential violation of the core MR assumptions and subsequent CAUSE analysis implicated pleiotropy of the genetic instruments used to proxy sexual behaviours. Finally, there was some attenuation of the univariable MR results in the multivariable models (AFS IVW OR = 0.7, 95%CI (0.4, 1.2), p = 0.21; NSP IVW OR = 0.9, 95%CI (0.5 1.7), p = 0.76).

CONCLUSIONS

Despite using genetic variants strongly related sexual behaviour traits in large-scale GWAS, we found evidence for correlated pleiotropy. This emphasizes a need for multivariable approaches and the triangulation of evidence when performing MR of complex behavioural traits.

Genetic Factors Associated with Prostate Cancer Conversion from Active Surveillance to Treatment

Men diagnosed with low-risk prostate cancer (PC) are increasingly electing active surveillance (AS) as their initial management strategy. While this may reduce the side effects of treatment for prostate cancer, many men on AS eventually convert to active treatment. PC is one of the most heritable cancers, and genetic factors that predispose to aggressive tumors may help distinguish men who are more likely to discontinue AS. To investigate this, we undertook a multi-institutional genome-wide association study (GWAS) of 5,222 PC patients and 1,139 other patients from replication cohorts, all of whom initially elected AS and were followed over time for the potential outcome of conversion from AS to active treatment. In the GWAS we detected 18 variants associated with conversion, 15 of which were not previously associated with PC risk. With a transcriptome-wide association study (TWAS), we found two genes associated with conversion (MAST3, p = 6.9×10-7 and GAB2, p = 2.0×10-6). Moreover, increasing values of a previously validated 269-variant genetic risk score (GRS) for PC was positively associated with conversion (e.g., comparing the highest to the two middle deciles gave a hazard ratio [HR] = 1.13; 95% Confidence Interval [CI]= 0.94-1.36); whereas, decreasing values of a 36-variant GRS for prostate-specific antigen (PSA) levels were positively associated with conversion (e.g., comparing the lowest to the two middle deciles gave a HR = 1.25; 95% CI, 1.04-1.50). These results suggest that germline genetics may help inform and individualize the decision of AS-or the intensity of monitoring on AS-versus treatment for the initial management of patients with low-risk PC.

Genetic determinants of blood-cell traits influence susceptibility to childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Despite overlap between genetic risk loci for ALL and hematologic traits, the etiological relevance of dysregulated blood-cell homeostasis remains unclear. We investigated this question in a genome-wide association study (GWAS) of childhood ALL (2,666 affected individuals, 60,272 control individuals) and a multi-trait GWAS of nine blood-cell indices in the UK Biobank. We identified 3,000 blood-cell-trait-associated (p < 5.0 × 10-8) variants, explaining 4.0% to 23.9% of trait variation and including 115 loci associated with blood-cell ratios (LMR, lymphocyte-to-monocyte ratio; NLR, neutrophil-to-lymphocyte ratio; PLR, platelet-to-lymphocyte ratio). ALL susceptibility was genetically correlated with lymphocyte counts (rg = 0.088, p = 4.0 × 10-4) and PLR (rg = -0.072, p = 0.0017). In Mendelian randomization analyses, genetically predicted increase in lymphocyte counts was associated with increased ALL risk (odds ratio [OR] = 1.16, p = 0.031) and strengthened after accounting for other cell types (OR = 1.43, p = 8.8 × 10-4). We observed positive associations with increasing LMR (OR = 1.22, p = 0.0017) and inverse effects for NLR (OR = 0.67, p = 3.1 × 10-4) and PLR (OR = 0.80, p = 0.002). Our study shows that a genetically induced shift toward higher lymphocyte counts, overall and in relation to monocytes, neutrophils, and platelets, confers an increased susceptibility to childhood ALL.

A Large-Scale Association Study Detects Novel Rare Variants, Risk Genes, Functional Elements, and Polygenic Architecture of Prostate Cancer Susceptibility

To identify rare variants associated with prostate cancer susceptibility and better characterize the mechanisms and cumulative disease risk associated with common risk variants, we conducted an integrated study of prostate cancer genetic etiology in two cohorts using custom genotyping microarrays, large imputation reference panels, and functional annotation approaches. Specifically, 11,984 men (6,196 prostate cancer cases and 5,788 controls) of European ancestry from Northern California Kaiser Permanente were genotyped and meta-analyzed with 196,269 men of European ancestry (7,917 prostate cancer cases and 188,352 controls) from the UK Biobank. Three novel loci, including two rare variants (European ancestry minor allele frequency < 0.01, at 3p21.31 and 8p12), were significant genome wide in a meta-analysis. Gene-based rare variant tests implicated a known prostate cancer gene (HOXB13), as well as a novel candidate gene (ILDR1), which encodes a receptor highly expressed in prostate tissue and is related to the B7/CD28 family of T-cell immune checkpoint markers. Haplotypic patterns of long-range linkage disequilibrium were observed for rare genetic variants at HOXB13 and other loci, reflecting their evolutionary history. In addition, a polygenic risk score (PRS) of 188 prostate cancer variants was strongly associated with risk (90th vs. 40th-60th percentile OR = 2.62, P = 2.55 × 10-191). Many of the 188 variants exhibited functional signatures of gene expression regulation or transcription factor binding, including a 6-fold difference in log-probability of androgen receptor binding at the variant rs2680708 (17q22). Rare variant and PRS associations, with concomitant functional interpretation of risk mechanisms, can help clarify the full genetic architecture of prostate cancer and other complex traits. SIGNIFICANCE: This study maps the biological relationships between diverse risk factors for prostate cancer, integrating different functional datasets to interpret and model genome-wide data from over 200,000 men with and without prostate cancer.See related commentary by Lachance, p. 1637.

Cross-cancer evaluation of polygenic risk scores for 16 cancer types in two large cohorts

Even distinct cancer types share biological hallmarks. Here, we investigate polygenic risk score (PRS)-specific pleiotropy across 16 cancers in European ancestry individuals from the Genetic Epidemiology Research on Adult Health and Aging cohort (16,012 cases, 50,552 controls) and UK Biobank (48,969 cases, 359,802 controls). Within cohorts, each PRS is evaluated in multivariable logistic regression models against all other cancer types. Results are then meta-analyzed across cohorts. Ten positive and one inverse cross-cancer associations are found after multiple testing correction. Two pairs show bidirectional associations; the melanoma PRS is positively associated with oral cavity/pharyngeal cancer and vice versa, whereas the lung cancer PRS is positively associated with oral cavity/pharyngeal cancer, and the oral cavity/pharyngeal cancer PRS is inversely associated with lung cancer. Overall, we validate known, and uncover previously unreported, patterns of pleiotropy that have the potential to inform investigations of risk prediction, shared etiology, and precision cancer prevention strategies.

Pan-cancer analysis demonstrates that integrating polygenic risk scores with modifiable risk factors improves risk prediction

Cancer risk is determined by a complex interplay of environmental and heritable factors. Polygenic risk scores (PRS) provide a personalized genetic susceptibility profile that may be leveraged for disease prediction. Using data from the UK Biobank (413,753 individuals; 22,755 incident cancer cases), we quantify the added predictive value of integrating cancer-specific PRS with family history and modifiable risk factors for 16 cancers. We show that incorporating PRS measurably improves prediction accuracy for most cancers, but the magnitude of this improvement varies substantially. We also demonstrate that stratifying on levels of PRS identifies significantly divergent 5-year risk trajectories after accounting for family history and modifiable risk factors. At the population level, the top 20% of the PRS distribution accounts for 4.0% to 30.3% of incident cancer cases, exceeding the impact of many lifestyle-related factors. In summary, this study illustrates the potential for improving cancer risk assessment by integrating genetic risk scores.

The landscape of host genetic factors involved in immune response to common viral infections

BACKGROUND

Humans and viruses have co-evolved for millennia resulting in a complex host genetic architecture. Understanding the genetic mechanisms of immune response to viral infection provides insight into disease etiology and therapeutic opportunities.

METHODS

We conducted a comprehensive study including genome-wide and transcriptome-wide association analyses to identify genetic loci associated with immunoglobulin G antibody response to 28 antigens for 16 viruses using serological data from 7924 European ancestry participants in the UK Biobank cohort.

RESULTS

Signals in human leukocyte antigen (HLA) class II region dominated the landscape of viral antibody response, with 40 independent loci and 14 independent classical alleles, 7 of which exhibited pleiotropic effects across viral families. We identified specific amino acid (AA) residues that are associated with seroreactivity, the strongest associations presented in a range of AA positions within DRβ1 at positions 11, 13, 71, and 74 for Epstein-Barr virus (EBV), Varicella zoster virus (VZV), human herpesvirus 7, (HHV7), and Merkel cell polyomavirus (MCV). Genome-wide association analyses discovered 7 novel genetic loci outside the HLA associated with viral antibody response (P < 5.0 × 10-8), including FUT2 (19q13.33) for human polyomavirus BK (BKV), STING1 (5q31.2) for MCV, and CXCR5 (11q23.3) and TBKBP1 (17q21.32) for HHV7. Transcriptome-wide association analyses identified 114 genes associated with response to viral infection, 12 outside of the HLA region, including ECSCR: P = 5.0 × 10-15 (MCV), NTN5: P = 1.1 × 10-9 (BKV), and P2RY13: P = 1.1 × 10-8 EBV nuclear antigen. We also demonstrated pleiotropy between viral response genes and complex diseases, from autoimmune disorders to cancer to neurodegenerative and psychiatric conditions.

CONCLUSIONS

Our study confirms the importance of the HLA region in host response to viral infection and elucidates novel genetic determinants beyond the HLA that contribute to host-virus interaction.

Investigation of Leukocyte Telomere Length and Genetic Variants in Chromosome 5p15.33 as Prognostic Markers in Lung Cancer

BACKGROUND

Lung cancer remains the leading cause of cancer mortality with relatively few prognostic biomarkers. We investigated associations with overall survival for telomere length (TL) and genetic variation in chromosome 5p15.33, an established telomere maintenance locus.

METHODS

Leukocyte TL was measured after diagnosis in 807 patients with non-small cell lung cancer (NSCLC) from the Princess Margaret Cancer Center in Toronto and assessed prospectively in 767 NSCLC cases from the Copenhagen City Heart Study and the Copenhagen General Population Study. Associations with all-cause mortality were tested for 723 variants in 5p15.33, genotyped in 4,672 NSCLC cases.

RESULTS

Short telomeres (≤10th percentile) were associated with poor prognosis for adenocarcinoma in both populations: TL measured 6 months after diagnosis [HR = 1.65; 95% confidence intervals (CI), 1.04-2.64] and for those diagnosed within 5 years after blood sampling (HR = 2.42; 95% CI, 1.37-4.28). Short TL was associated with mortality in never smokers with NSCLC (HR = 10.29; 95% CI, 1.86-56.86) and adenocarcinoma (HR = 11.31; 95% CI, 1.96-65.24). Analyses in 5p15.33 identified statistically significant prognostic associations for rs56266421-G in LPCAT1 (HR = 1.86; 95% CI, 1.38-2.52; P = 4.5 × 10^-5) in stage I-IIIA NSCLC, and for the SLC6A3 gene with OS in females with NSCLC (P = 1.6 × 10^-3).

CONCLUSIONS

Our findings support the potential clinical utility of TL, particularly for adenocarcinoma patients, while associations in chromosome 5p15.33 warrant further exploration.

IMPACT

This is the largest lung cancer study of leukocyte TL and OS, and the first to examine the impact of the timing of TL measurement. Our findings suggest that extremely short telomeres are indicative of poor prognosis in NSCLC.

The landscape of host genetic factors involved in immune response to common viral infections

INTRODUCTION

Humans and viruses have co-evolved for millennia resulting in a complex host genetic architecture. Understanding the genetic mechanisms of immune response to viral infection provides insight into disease etiology and therapeutic opportunities.

METHODS

We conducted a comprehensive study including genome-wide and transcriptome-wide association analyses to identify genetic loci associated with immunoglobulin G antibody response to 28 antigens for 16 viruses using serological data from 7924 European ancestry participants in the UK Biobank cohort.

RESULTS

Signals in human leukocyte antigen (HLA) class II region dominated the landscape of viral antibody response, with 40 independent loci and 14 independent classical alleles, 7 of which exhibited pleiotropic effects across viral families. We identified specific amino acid (AA) residues that are associated with seroreactivity, the strongest associations presented in a range of AA positions within DRβ1 at positions 11, 13, 71, and 74 for Epstein-Barr Virus (EBV), Varicella Zoster Virus (VZV), Human Herpes virus 7, (HHV7) and Merkel cell polyomavirus (MCV). Genome-wide association analyses discovered 7 novel genetic loci outside the HLA associated with viral antibody response (P<5.×10-8), including FUT2 (19q13.33) for human polyomavirus BK (BKV), STING1 (5q31.2) for MCV, as well as CXCR5 (11q23.3) and TBKBP1 (17q21.32) for HHV7. Transcriptome-wide association analyses identified 114 genes associated with response to viral infection, 12 outside of the HLA region, including ECSCR: P=5.0×10-15 (MCV), NTN5: P=1.1×10-9 (BKV), and P2RY13: P=1.1×10-8 EBV nuclear antigen. We also demonstrated pleiotropy between viral response genes and complex diseases; from autoimmune disorders to cancer to neurodegenerative and psychiatric conditions.

CONCLUSIONS

Our study confirms the importance of the HLA region in host response to viral infection and elucidates novel genetic determinants beyond the HLA that contribute to host-virus interaction.

Pan-cancer study detects genetic risk variants and shared genetic basis in two large cohorts

Deciphering the shared genetic basis of distinct cancers has the potential to elucidate carcinogenic mechanisms and inform broadly applicable risk assessment efforts. Here, we undertake genome-wide association studies (GWAS) and comprehensive evaluations of heritability and pleiotropy across 18 cancer types in two large, population-based cohorts: the UK Biobank (408,786 European ancestry individuals; 48,961 cancer cases) and the Kaiser Permanente Genetic Epidemiology Research on Adult Health and Aging cohorts (66,526 European ancestry individuals; 16,001 cancer cases). The GWAS detect 21 genome-wide significant associations independent of previously reported results. Investigations of pleiotropy identify 12 cancer pairs exhibiting either positive or negative genetic correlations; 25 pleiotropic loci; and 100 independent pleiotropic variants, many of which are regulatory elements and/or influence cross-tissue gene expression. Our findings demonstrate widespread pleiotropy and offer further insight into the complex genetic architecture of cross-cancer susceptibility.

Insecticide use and risk of non-Hodgkin lymphoma subtypes: A subset meta-analysis of the North American Pooled Project

Insecticide use has been linked to increased risk of non-Hodgkin lymphoma (NHL), however, findings of epidemiologic studies have been inconsistent, particularly for NHL subtypes. We analyzed 1690 NHL cases and 5131 controls in the North American Pooled Project (NAPP) to investigate self-reported insecticide use and risk of NHL overall and by subtypes: follicular lymphoma (FL), diffuse large B-cell lymphoma (DLBCL) and small lymphocytic lymphoma (SLL). Odds ratios (OR) and 95% confidence intervals for each insecticide were estimated using logistic regression. Subtype-specific associations were evaluated using ASSET (Association analysis for SubSETs). Increased risks of multiple NHL subtypes were observed for lindane (OR = 1.60, 1.20-2.10: FL, DLCBL, SLL), chlordane (OR = 1.59, 1.17-2.16: FL, SLL) and DDT (OR = 1.36, 1.06-1.73: DLBCL, SLL). Positive trends were observed, within the subsets with identified associations, for increasing categories of exposure duration for lindane (P_trend = 1.7 × 10^-4 ), chlordane (P_trend = 1.0 × 10^-3 ) and DDT (P_trend = 4.2 × 10^-3 ), however, the exposure-response relationship was nonlinear. Ever use of pyrethrum was associated with an increased risk of FL (OR = 3.65, 1.45-9.15), and the relationship with duration of use appeared monotonic (OR for >10 years: OR = 5.38, 1.75-16.53; P_trend = 3.6 × 10^-3 ). Our analysis identified several novel associations between insecticide use and specific NHL subtypes, suggesting possible etiologic heterogeneity in the context of pesticide exposure.

Silica and asbestos exposure at work and the risk of bladder cancer in Canadian men: a population-based case-control study

Background

Silica and asbestos are recognized lung carcinogens. However, their role in carcinogenesis at other organs is less clear. Clearance of inhaled silica particles and asbestos fibers from the lungs may lead to translocation to sites such as the bladder where they may initiate carcinogenesis. We used data from a Canadian population-based case-control study to evaluate the associations between these workplace exposures and bladder cancer.

Methods

Data from a population-based case-control study were used to characterize associations between workplace exposure to silica and asbestos and bladder cancer among men. Bladder cancer cases (N = 658) and age-frequency matched controls (N = 1360) were recruited within the National Enhanced Cancer Surveillance System from eight Canadian provinces (1994–97). Exposure concentration, frequency and reliability for silica and asbestos were assigned to each job, based on lifetime occupational histories, using a combination of job-exposure profiles and expert review. Exposure was modeled as ever/never, highest attained concentration, duration (years), highest attained frequency (% worktime) and cumulative exposure. Odds ratios (OR) and their 95% confidence intervals (CI) were estimated using adjusted logistic regression.

Results

A modest (approximately 20%) increase in bladder cancer risk was found for ever having been exposed to silica, highest attained concentration and frequency of exposure but this increase was not statistically significant. Relative to unexposed, the odds of bladder cancer were 1.41 (95%CI: 1.01–1.98) times higher among men exposed to silica at work for ≥27 years. For asbestos, relative to unexposed, an increased risk of bladder cancer was observed for those first exposed ≥20 years ago (OR:2.04, 95%CI:1.25–3.34), those with a frequency of exposure of 5–30% of worktime (OR:1.45, 95%CI:1.06–1.98), and for those with < 10 years of exposure at low concentrations (OR:1.75, 95%CI:1.10–2.77) and the lower tertile of cumulative exposure (OR:1.69, 95%CI:1.07–2.65). However, no clear exposure-response relationships emerged.

Conclusions

Our results indicate a slight increase in risk of bladder cancer with exposure to silica and asbestos, suggesting that the effects of these agents are broader than currently recognized. The findings from this study inform evidence-based action to enhance cancer prevention efforts, particularly for workers in industries with regular exposure.

Immune-mediated genetic pathways resulting in pulmonary function impairment increase lung cancer susceptibility

Impaired lung function is often caused by cigarette smoking, making it challenging to disentangle its role in lung cancer susceptibility. Investigation of the shared genetic basis of these phenotypes in the UK Biobank and International Lung Cancer Consortium (29,266 cases, 56,450 controls) shows that lung cancer is genetically correlated with reduced forced expiratory volume in one second (FEV1: rg = 0.098, p = 2.3 × 10-8) and the ratio of FEV1 to forced vital capacity (FEV1/FVC: rg = 0.137, p = 2.0 × 10-12). Mendelian randomization analyses demonstrate that reduced FEV1 increases squamous cell carcinoma risk (odds ratio (OR) = 1.51, 95% confidence intervals: 1.21-1.88), while reduced FEV1/FVC increases the risk of adenocarcinoma (OR = 1.17, 1.01-1.35) and lung cancer in never smokers (OR = 1.56, 1.05-2.30). These findings support a causal role of pulmonary impairment in lung cancer etiology. Integrative analyses reveal that pulmonary function instruments, including 73 novel variants, influence lung tissue gene expression and implicate immune-related pathways in mediating the observed effects on lung carcinogenesis.

Occupational Exposure to Diesel and Gasoline Engine Exhausts and the Risk of Kidney Cancer in Canadian Men

Introduction

Kidney cancer is the fifth most common incident cancer in Canadian men. Diesel and gasoline exhausts are common workplace exposures that have been examined as risk factors for non-lung cancer sites, including the kidney, but limitations in exposure assessment methods have contributed to inconsistent findings. The objective of this study was to assess the relationship between occupational gasoline and diesel engine exhausts and the risk of kidney cancer in men.

Methods

The National Enhanced Cancer Surveillance System (NECSS) is a Canadian population-based case–control study conducted in 1994–1997. Incident kidney cancer cases were identified using provincial registries, while the control series was identified through random-digit dialing, or provincial administrative databases. Self-reported questionnaires were used to obtain information on lifetime occupational history and cancer risk factors. Two hygienists, blinded to case status, coded occupational histories for diesel and gasoline exhaust exposures using concentration, frequency, duration, and reliability. Logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs) separately by exhaust type. The separate and combined impacts of both engine exhausts were also examined. ORs were adjusted for age, province, body mass index, occupational secondhand smoke exposure, and education.

Results

Of the kidney cancer cases (n = 712), 372 (52%) had exposure to both exhausts at some point, and 984 (40%) of the controls (n = 2457) were ever exposed. Workers who had ever been exposed to engine exhausts were more likely to have kidney cancer than those who were never exposed (OR diesel = 1.23, 95% CI = 0.99–1.53; OR gasoline = 1.51, 95% CI = 1.23–1.86). Exposure to gasoline exhaust was consistently associated with kidney cancer in a dose–response manner (P value for trends in highest attained and cumulative exposure both <0.0001). Those men with high cumulative exposure to both gasoline and diesel exhaust had a 76% increased odds of kidney cancer (95% CI = 1.27–2.43).

Conclusions

This study provides evidence that occupational gasoline, and to a lesser extent, diesel exhaust exposure may increase the risk of kidney cancer.

Author Correction: Association of imputed prostate cancer transcriptome with disease risk reveals novel mechanisms

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Exposure to crystalline silica in Canadian workplaces and the risk of kidney cancer

OBJECTIVES

The causes of kidney cancer are not well understood though occupational exposures are thought to play a role. Crystalline silica is a known human carcinogen, and despite previous links with kidney disease, there have been few studies investigating its association with kidney cancer. We addressed this research gap using a population-based case-control study of Canadian men.

METHODS

Questionnaire data were obtained from individuals with histologically confirmed kidney cancer, and population-based controls recruited from eight Canadian provinces (1994-1997). An industrial hygienist characterised participants' lifetime occupational exposure, and their confidence in the assessment (possibly, probably or definitely exposed) to silica on three dimensions (intensity, frequency and duration), and cumulative exposure was estimated. Logistic regression was used to estimate ORs and 95% CIs, adjusting for potential confounders.

RESULTS

Nearly half of the 689 kidney cancer cases (49%) and 2369 controls (44%) had ever been occupationally exposed to crystalline silica. In a fully adjusted model, workers ever-exposed to silica had a slightly increased risk of kidney cancer relative to those who were unexposed (OR 1.10, 95% CI 0.92 to 1.32). Odds were modestly (and generally not statistically significantly) increased for models with duration of exposure and cumulative exposure, though exposure-response relationships were not evident.

CONCLUSIONS

Our findings do not provide evidence that occupational exposure to crystalline silica increases risk of kidney cancer in men.

Non-Hodgkin lymphoma risk and organophosphate and carbamate insecticide use in the north American pooled project

Organophosphates and carbamates have been among the most commonly used insecticides, with both agricultural and residential uses. Previous studies have suggested associations of non-Hodgkin lymphoma (NHL) with some of these chemicals; however, many studies have been limited in their ability to evaluate associations with lymphoma subtypes. We evaluated the use of eleven organophosphate and two carbamate insecticides in association with NHL in the North American Pooled Project, which includes data from case-control studies in the United States and Canada (1690 cases/5131 controls). We used unconditional logistic regression adjusting for potential confounders, including use of other pesticides, to estimate odds ratios (OR) and 95% confidence intervals (CI) for associations between these chemicals and NHL overall, and NHL subtypes, i.e., follicular (FL), diffuse large B-cell (DLBCL), small lymphocytic lymphoma (SLL) and others. Ever use of malathion was associated with increased risk of NHL overall (OR = 1.43; 95% CI: 1.14-1.81) compared with never users. Categories using tertiles of duration (<4 yrs., 4-12 yrs., and >12 yrs) also showed a significant exposure-response for increasing years of use of malathion and risk of NHL (OR_<4vsUnex = 1.33 (0.88, 2.03), OR_4-12vsUnex = 1.42 (1.02, 1.96), OR_>12vsUnex = 1.55 (1.05, 2.28, p-trend < 0.01)). In addition, malathion use was statistically significantly associated with FL (OR = 1.58; 95% CI: 1.11-2.27) and DLBCL (OR = 1.61; 95% CI: 1.16-2.22) while there were no apparent associations with SLL or other subtypes, the p-value for heterogeneity across subtypes, however, was not significant. These results support previous studies suggesting an association between insecticide use and NHL overall, and provide new information on associations with NHL subtypes.

Mendelian Randomization and mediation analysis of leukocyte telomere length and risk of lung and head and neck cancers

BACKGROUND

Evidence from observational studies of telomere length (TL) has been conflicting regarding its direction of association with cancer risk. We investigated the causal relevance of TL for lung and head and neck cancers using Mendelian Randomization (MR) and mediation analyses.

METHODS

We developed a novel genetic instrument for TL in chromosome 5p15.33, using variants identified through deep-sequencing, that were genotyped in 2051 cancer-free subjects. Next, we conducted an MR analysis of lung (16 396 cases, 13 013 controls) and head and neck cancer (4415 cases, 5013 controls) using eight genetic instruments for TL. Lastly, the 5p15.33 instrument and distinct 5p15.33 lung cancer risk loci were evaluated using two-sample mediation analysis, to quantify their direct and indirect, telomere-mediated, effects.

RESULTS

The multi-allelic 5p15.33 instrument explained 1.49-2.00% of TL variation in our data (p = 2.6 × 10-9). The MR analysis estimated that a 1000 base-pair increase in TL increases risk of lung cancer [odds ratio (OR) = 1.41, 95% confidence interval (CI): 1.20-1.65] and lung adenocarcinoma (OR = 1.92, 95% CI: 1.51-2.22), but not squamous lung carcinoma (OR = 1.04, 95% CI: 0.83-1.29) or head and neck cancers (OR = 0.90, 95% CI: 0.70-1.05). Mediation analysis of the 5p15.33 instrument indicated an absence of direct effects on lung cancer risk (OR = 1.00, 95% CI: 0.95-1.04). Analysis of distinct 5p15.33 susceptibility variants estimated that TL mediates up to 40% of the observed associations with lung cancer risk.

CONCLUSIONS

Our findings support a causal role for long telomeres in lung cancer aetiology, particularly for adenocarcinoma, and demonstrate that telomere maintenance partially mediates the lung cancer susceptibility conferred by 5p15.33 loci.

Elevated Platelet Count Appears to Be Causally Associated with Increased Risk of Lung Cancer: A Mendelian Randomization Analysis

BACKGROUND

Platelets are a critical element in coagulation and inflammation, and activated platelets are linked to cancer risk through diverse mechanisms. However, a causal relationship between platelets and risk of lung cancer remains unclear.

METHODS

We performed single and combined multiple instrumental variable Mendelian randomization analysis by an inverse-weighted method, in addition to a series of sensitivity analyses. Summary data for associations between SNPs and platelet count are from a recent publication that included 48,666 Caucasian Europeans, and the International Lung Cancer Consortium and Transdisciplinary Research in Cancer of the Lung data consisting of 29,266 cases and 56,450 controls to analyze associations between candidate SNPs and lung cancer risk.

RESULTS

Multiple instrumental variable analysis incorporating six SNPs showed a 62% increased risk of overall non-small cell lung cancer [NSCLC; OR, 1.62; 95% confidence interval (CI), 1.15-2.27; P = 0.005] and a 200% increased risk for small-cell lung cancer (OR, 3.00; 95% CI, 1.27-7.06; P = 0.01). Results showed only a trending association with NSCLC histologic subtypes, which may be due to insufficient sample size and/or weak effect size. A series of sensitivity analysis retained these findings.

CONCLUSIONS

Our findings suggest a causal relationship between elevated platelet count and increased risk of lung cancer and provide evidence of possible antiplatelet interventions for lung cancer prevention.

IMPACT

These findings provide a better understanding of lung cancer etiology and potential evidence for antiplatelet interventions for lung cancer prevention.