Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Epidemiologic Factors in Relation to Colorectal Cancer Risk and Survival by Genotoxic Colibactin Mutational Signature

BACKGROUND

The genotoxin colibactin causes a tumor single-base substitution (SBS) mutational signature, SBS88. It is unknown whether epidemiologic factors' association with colorectal cancer risk and survival differs by SBS88.

METHODS

Within the Genetic Epidemiology of Colorectal Cancer Consortium and Colon Cancer Family Registry, we measured SBS88 in 4,308 microsatellite stable/microsatellite instability low tumors. Associations of epidemiologic factors with colorectal cancer risk by SBS88 were assessed using multinomial regression (N = 4,308 cases, 14,192 controls; cohort-only cases N = 1,911), and with colorectal cancer-specific survival using Cox proportional hazards regression (N = 3,465 cases).

RESULTS

392 (9%) tumors were SBS88 positive. Among all cases, the highest quartile of fruit intake was associated with lower risk of SBS88-positive colorectal cancer than SBS88-negative colorectal cancer [odds ratio (OR) = 0.53, 95% confidence interval (CI) 0.37-0.76; OR = 0.75, 95% CI 0.66-0.85, respectively, Pheterogeneity = 0.047]. Among cohort studies, associations of body mass index (BMI), alcohol, and fruit intake with colorectal cancer risk differed by SBS88. BMI ≥30 kg/m2 was associated with worse colorectal cancer-specific survival among those SBS88-positive [hazard ratio (HR) = 3.40, 95% CI 1.47-7.84], but not among those SBS88-negative (HR = 0.97, 95% CI 0.78-1.21, Pheterogeneity = 0.066).

CONCLUSIONS

Most epidemiologic factors did not differ by SBS88 for colorectal cancer risk or survival. Higher BMI may be associated with worse colorectal cancer-specific survival among those SBS88-positive; however, validation is needed in samples with whole-genome or whole-exome sequencing available.

IMPACT

This study highlights the importance of identification of tumor phenotypes related to colorectal cancer and understanding potential heterogeneity for risk and survival.

Genetic risk impacts the association of menopausal hormone therapy with colorectal cancer risk

BACKGROUND

Menopausal hormone therapy (MHT), a common treatment to relieve symptoms of menopause, is associated with a lower risk of colorectal cancer (CRC). To inform CRC risk prediction and MHT risk-benefit assessment, we aimed to evaluate the joint association of a polygenic risk score (PRS) for CRC and MHT on CRC risk.

METHODS

We used data from 28,486 postmenopausal women (11,519 cases and 16,967 controls) of European descent. A PRS based on 141 CRC-associated genetic variants was modeled as a categorical variable in quartiles. Multiplicative interaction between PRS and MHT use was evaluated using logistic regression. Additive interaction was measured using the relative excess risk due to interaction (RERI). 30-year cumulative risks of CRC for 50-year-old women according to MHT use and PRS were calculated.

RESULTS

The reduction in odds ratios by MHT use was larger in women within the highest quartile of PRS compared to that in women within the lowest quartile of PRS (p-value = 2.7 × 10-8). At the highest quartile of PRS, the 30-year CRC risk was statistically significantly lower for women taking any MHT than for women not taking any MHT, 3.7% (3.3%-4.0%) vs 6.1% (5.7%-6.5%) (difference 2.4%, P-value = 1.83 × 10-14); these differences were also statistically significant but smaller in magnitude in the lowest PRS quartile, 1.6% (1.4%-1.8%) vs 2.2% (1.9%-2.4%) (difference 0.6%, P-value = 1.01 × 10-3), indicating 4 times greater reduction in absolute risk associated with any MHT use in the highest compared to the lowest quartile of genetic CRC risk.

CONCLUSIONS

MHT use has a greater impact on the reduction of CRC risk for women at higher genetic risk. These findings have implications for the development of risk prediction models for CRC and potentially for the consideration of genetic information in the risk-benefit assessment of MHT use.

Body size and risk of colorectal cancer molecular defined subtypes and pathways: Mendelian randomization analyses

BACKGROUND

Obesity has been positively associated with most molecular subtypes of colorectal cancer (CRC); however, the magnitude and the causality of these associations is uncertain.

METHODS

We used Mendelian randomization (MR) to examine potential causal relationships between body size traits (body mass index [BMI], waist circumference, and body fat percentage) with risks of Jass classification types and individual subtypes of CRC (microsatellite instability [MSI] status, CpG island methylator phenotype [CIMP] status, BRAF and KRAS mutations). Summary data on tumour markers were obtained from two genetic consortia (CCFR, GECCO).

FINDINGS

A 1-standard deviation (SD:5.1 kg/m2) increment in BMI levels was found to increase risks of Jass type 1MSI-high,CIMP-high,BRAF-mutated,KRAS-wildtype (odds ratio [OR]: 2.14, 95% confidence interval [CI]: 1.46, 3.13; p-value = 9 × 10-5) and Jass type 2non-MSI-high,CIMP-high,BRAF-mutated,KRAS-wildtype CRC (OR: 2.20, 95% CI: 1.26, 3.86; p-value = 0.005). The magnitude of these associations was stronger compared with Jass type 4non-MSI-high,CIMP-low/negative,BRAF-wildtype,KRAS-wildtype CRC (p-differences: 0.03 and 0.04, respectively). A 1-SD (SD:13.4 cm) increment in waist circumference increased risk of Jass type 3non-MSI-high,CIMP-low/negative,BRAF-wildtype,KRAS-mutated (OR 1.73, 95% CI: 1.34, 2.25; p-value = 9 × 10-5) that was stronger compared with Jass type 4 CRC (p-difference: 0.03). A higher body fat percentage (SD:8.5%) increased risk of Jass type 1 CRC (OR: 2.59, 95% CI: 1.49, 4.48; p-value = 0.001), which was greater than Jass type 4 CRC (p-difference: 0.03).

INTERPRETATION

Body size was more strongly linked to the serrated (Jass types 1 and 2) and alternate (Jass type 3) pathways of colorectal carcinogenesis in comparison to the traditional pathway (Jass type 4).

FUNDING

Cancer Research UK, National Institute for Health Research, Medical Research Council, National Institutes of Health, National Cancer Institute, American Institute for Cancer Research, Brigham and Women's Hospital, Prevent Cancer Foundation, Victorian Cancer Agency, Swedish Research Council, Swedish Cancer Society, Region Västerbotten, Knut and Alice Wallenberg Foundation, Lion's Cancer Research Foundation, Insamlingsstiftelsen, Umeå University. Full funding details are provided in acknowledgements.

Genome-Wide Gene-Environment Interaction Analyses to Understand the Relationship between Red Meat and Processed Meat Intake and Colorectal Cancer Risk

BACKGROUND

High red meat and/or processed meat consumption are established colorectal cancer risk factors. We conducted a genome-wide gene-environment (GxE) interaction analysis to identify genetic variants that may modify these associations.

METHODS

A pooled sample of 29,842 colorectal cancer cases and 39,635 controls of European ancestry from 27 studies were included. Quantiles for red meat and processed meat intake were constructed from harmonized questionnaire data. Genotyping arrays were imputed to the Haplotype Reference Consortium. Two-step EDGE and joint tests of GxE interaction were utilized in our genome-wide scan.

RESULTS

Meta-analyses confirmed positive associations between increased consumption of red meat and processed meat with colorectal cancer risk [per quartile red meat OR = 1.30; 95% confidence interval (CI) = 1.21-1.41; processed meat OR = 1.40; 95% CI = 1.20-1.63]. Two significant genome-wide GxE interactions for red meat consumption were found. Joint GxE tests revealed the rs4871179 SNP in chromosome 8 (downstream of HAS2); greater than median of consumption ORs = 1.38 (95% CI = 1.29-1.46), 1.20 (95% CI = 1.12-1.27), and 1.07 (95% CI = 0.95-1.19) for CC, CG, and GG, respectively. The two-step EDGE method identified the rs35352860 SNP in chromosome 18 (SMAD7 intron); greater than median of consumption ORs = 1.18 (95% CI = 1.11-1.24), 1.35 (95% CI = 1.26-1.44), and 1.46 (95% CI = 1.26-1.69) for CC, CT, and TT, respectively.

CONCLUSIONS

We propose two novel biomarkers that support the role of meat consumption with an increased risk of colorectal cancer.

IMPACT

The reported GxE interactions may explain the increased risk of colorectal cancer in certain population subgroups.

Impact of weight loss on cancer-related proteins in serum: results from a cluster randomised controlled trial of individuals with type 2 diabetes

BACKGROUND

Type 2 diabetes is associated with higher risk of several cancer types. However, the biological intermediates driving this relationship are not fully understood. As novel interventions for treating and managing type 2 diabetes become increasingly available, whether they also disrupt the pathways leading to increased cancer risk is currently unknown. We investigated the effect of a type 2 diabetes intervention, in the form of intentional weight loss, on circulating proteins associated with cancer risk to gain insight into potential mechanisms linking type 2 diabetes and adiposity with cancer development.

METHODS

Fasting serum samples from participants with diabetes enrolled in the Diabetes Remission Clinical Trial (DiRECT) receiving the Counterweight-Plus weight-loss programme (intervention, N = 117, mean weight-loss 10 kg, 46% diabetes remission) or best-practice care by guidelines (control, N = 143, mean weight-loss 1 kg, 4% diabetes remission) were subject to proteomic analysis using the Olink Oncology-II platform (48% of participants were female; 52% male). To identify proteins which may be altered by the weight-loss intervention, the difference in protein levels between groups at baseline and 1 year was examined using linear regression. Mendelian randomization (MR) was performed to extend these results to evaluate cancer risk and elucidate possible biological mechanisms linking type 2 diabetes and cancer development. MR analyses were conducted using independent datasets, including large cancer meta-analyses, UK Biobank, and FinnGen, to estimate potential causal relationships between proteins modified during intentional weight loss and the risk of colorectal, breast, endometrial, gallbladder, liver, and pancreatic cancers.

FINDINGS

Nine proteins were modified by the intervention: glycoprotein Nmb; furin; Wnt inhibitory factor 1; toll-like receptor 3; pancreatic prohormone; erb-b2 receptor tyrosine kinase 2; hepatocyte growth factor; endothelial cell specific molecule 1 and Ret proto-oncogene (Holm corrected P-value <0.05). Mendelian randomization analyses indicated a causal relationship between predicted circulating furin and glycoprotein Nmb on breast cancer risk (odds ratio (OR) = 0.81, 95% confidence interval (CI) = 0.67-0.99, P-value = 0.03; and OR = 0.88, 95% CI = 0.78-0.99, P-value = 0.04 respectively), though these results were not supported in sensitivity analyses examining violations of MR assumptions.

INTERPRETATION

Intentional weight loss among individuals with recently diagnosed diabetes may modify levels of cancer-related proteins in serum. Further evaluation of the proteins identified in this analysis could reveal molecular pathways that mediate the effect of adiposity and type 2 diabetes on cancer risk.

FUNDING

The main sources of funding for this work were Diabetes UK, Cancer Research UK, World Cancer Research Fund, and Wellcome.

Selection, optimization and validation of ten chronic disease polygenic risk scores for clinical implementation in diverse US populations

Polygenic risk scores (PRSs) have improved in predictive performance, but several challenges remain to be addressed before PRSs can be implemented in the clinic, including reduced predictive performance of PRSs in diverse populations, and the interpretation and communication of genetic results to both providers and patients. To address these challenges, the National Human Genome Research Institute-funded Electronic Medical Records and Genomics (eMERGE) Network has developed a framework and pipeline for return of a PRS-based genome-informed risk assessment to 25,000 diverse adults and children as part of a clinical study. From an initial list of 23 conditions, ten were selected for implementation based on PRS performance, medical actionability and potential clinical utility, including cardiometabolic diseases and cancer. Standardized metrics were considered in the selection process, with additional consideration given to strength of evidence in African and Hispanic populations. We then developed a pipeline for clinical PRS implementation (score transfer to a clinical laboratory, validation and verification of score performance), and used genetic ancestry to calibrate PRS mean and variance, utilizing genetically diverse data from 13,475 participants of the All of Us Research Program cohort to train and test model parameters. Finally, we created a framework for regulatory compliance and developed a PRS clinical report for return to providers and for inclusion in an additional genome-informed risk assessment. The initial experience from eMERGE can inform the approach needed to implement PRS-based testing in diverse clinical settings.

Genotoxic colibactin mutational signature in colorectal cancer is associated with clinicopathological features, specific genomic alterations and better survival

Background and Aims

The microbiome has long been suspected of a role in colorectal cancer (CRC) tumorigenesis. The mutational signature SBS88 mechanistically links CRC development with the strain of Escherichia coli harboring the pks island that produces the genotoxin colibactin, but the genomic, pathological and survival characteristics associated with SBS88-positive tumors are unknown.

Methods

SBS88-positive CRCs were identified from targeted sequencing data from 5,292 CRCs from 17 studies and tested for their association with clinico-pathological features, oncogenic pathways, genomic characteristics and survival.

Results

In total, 7.5% (398/5,292) of the CRCs were SBS88-positive, of which 98.7% (392/398) were microsatellite stable/microsatellite instability low (MSS/MSI-L), compared with 80% (3916/4894) of SBS88 negative tumors (p=1.5x10-28). Analysis of MSS/MSI-L CRCs demonstrated that SBS88 positive CRCs were associated with the distal colon (OR=1.84, 95% CI=1.40-2.42, p=1x10-5) and rectum (OR=1.90, 95% CI=1.44-2.51, p=6x10-6) tumor sites compared with the proximal colon. The top seven recurrent somatic mutations associated with SBS88-positive CRCs demonstrated mutational contexts associated with colibactin-induced DNA damage, the strongest of which was the APC:c.835-8A>G mutation (OR=65.5, 95%CI=39.0-110.0, p=3x10-80). Large copy number alterations (CNAs) including CNA loss on 14q and gains on 13q, 16q and 20p were significantly enriched in SBS88-positive CRCs. SBS88-positive CRCs were associated with better CRC-specific survival (p=0.007; hazard ratio of 0.69, 95% CI=0.52-0.90) when stratified by age, sex, study, and by stage.

Conclusion

SBS88-positivity, a biomarker of colibactin-induced DNA damage, can identify a novel subtype of CRC characterized by recurrent somatic mutations, copy number alterations and better survival. These findings provide new insights for treatment and prevention strategies for this subtype of CRC.

FORGEdb: a tool for identifying candidate functional variants and uncovering target genes and mechanisms for complex diseases

The majority of disease-associated variants identified through genome-wide association studies are located outside of protein-coding regions. Prioritizing candidate regulatory variants and gene targets to identify potential biological mechanisms for further functional experiments can be challenging. To address this challenge, we developed FORGEdb ( https://forgedb.cancer.gov/ ; https://forge2.altiusinstitute.org/files/forgedb.html ; and https://doi.org/10.5281/zenodo.10067458 ), a standalone and web-based tool that integrates multiple datasets, delivering information on associated regulatory elements, transcription factor binding sites, and target genes for over 37 million variants. FORGEdb scores provide researchers with a quantitative assessment of the relative importance of each variant for targeted functional experiments.

Obesity and prostate cancer screening, incidence, and mortality in the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial

BACKGROUND

Though obesity, measured by body mass index (BMI), is an established risk factor for several cancer sites, there is conflicting evidence on whether obesity increases prostate cancer risk or mortality and, if it does, whether it increases risk directly or indirectly by affecting prostate cancer screening efficacy.

METHODS

We examined associations between BMI and prostate cancer screening outcomes, incidence, and mortality in men randomly assigned to the intervention arm of the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial (n = 36 756) between 1993 and 2001. Participants received annual screening with the prostate-specific antigen test and digital rectal exam. Associations between baseline BMI and screening outcomes were assessed via multinomial logistic regression, and associations with prostate cancer incidence and mortality were assessed via Cox proportional hazards regression.

RESULTS

Individuals with higher BMI were less likely to screen positive via the prostate-specific antigen test and/or digital rectal exam and more likely to have an inadequate screen (all Ptrend < .01). Higher BMI was inversely associated with prostate cancer incidence (per 5 kg/m2 BMI increase: hazard ratio [HR] = 0.94, 95% confidence interval [CI] = 0.91 to 0.97), including incidence of early stage (HR = 0.94, 95% CI = 0.90 to 0.97) and advanced-stage (HR = 0.91, 95% CI = 0.82 to 1.02) disease, but positively associated with prostate cancer mortality (HR = 1.21, 95% CI = 1.06 to 1.37). The association with mortality was not modified by screening outcome (Pinteraction = .13).

CONCLUSIONS

Within this screened population, individuals with higher BMI had lower risk of prostate cancer diagnosis but higher risk of prostate cancer mortality. As higher BMI was not positively associated with advanced-stage prostate cancer risk, the increased mortality is unlikely to be due to delayed prostate cancer detection.

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.

Glyphosate Use and Mosaic Loss of Chromosome Y among Male Farmers in the Agricultural Health Study

BACKGROUND

Glyphosate is the most commonly used herbicide worldwide and has been implicated in the development of certain hematologic cancers. Although mechanistic studies in human cells and animals support the genotoxic effects of glyphosate, evidence in human populations is scarce.

OBJECTIVES

We evaluated the association between lifetime occupational glyphosate use and mosaic loss of chromosome Y (mLOY) as a marker of genotoxicity among male farmers.

METHODS

We analyzed blood-derived DNA from 1,606 farmers ≥ 50 years of age in the Biomarkers of Exposure and Effect in Agriculture study, a subcohort of the Agricultural Health Study. mLOY was detected using genotyping array intensity data in the pseudoautosomal region of the sex chromosomes. Cumulative lifetime glyphosate use was assessed using self-reported pesticide exposure histories. Using multivariable logistic regression, we estimated odds ratios (ORs) and 95% confidence intervals (CIs) for the associations between glyphosate use and any detectable mLOY (overall mLOY) or mLOY affecting ≥ 10 % of cells (expanded mLOY).

RESULTS

Overall, mLOY was detected in 21.4% of farmers, and 9.8% of all farmers had expanded mLOY. Increasing total lifetime days of glyphosate use was associated with expanded mLOY [highest vs. lowest quartile; OR = 1.75 (95% CI: 1.00, 3.07), p trend = 0.03 ] but not with overall mLOY; the associations with expanded mLOY were most apparent among older (≥ 70 years of age) men [OR = 2.30 (95% CI: 1.13, 4.67), p trend = 0.01 ], never smokers [OR = 2.32 (95% CI: 1.04, 5.21), p trend = 0.04 ], and nonobese men [OR = 2.04 (95% CI: 0.99, 4.19), p trend = 0.03 ]. Similar patterns of associations were observed for intensity-weighted lifetime days of glyphosate use.

DISCUSSION

High lifetime glyphosate use could be associated with mLOY affecting a larger fraction of cells, suggesting glyphosate could confer genotoxic or selective effects relevant for clonal expansion. As the first study to investigate this association, our findings contribute novel evidence regarding the carcinogenic potential of glyphosate and require replication in future studies. https://doi.org/10.1289/EHP12834.

Joint associations between established genetic susceptibility loci, pesticide exposures, and risk of prostate cancer

More than 200 genetic variants have been independently associated with prostate cancer risk. Studies among farmers have also observed increased prostate cancer risk associated with exposure to specific organophosphate (fonofos, terbufos, malathion, dimethoate) and organochlorine (aldrin, chlordane) insecticides. We examined the joint associations between these pesticides, established prostate cancer loci, and prostate cancer risk among 1,162 cases (588 aggressive) and 2,206 frequency-matched controls nested in the Agricultural Health Study cohort. History of lifetime pesticide use was combined with a polygenic risk score (PRS) generated using 256 established prostate cancer risk variants. Logistic regression models estimated the joint associations of the pesticides, the PRS, and the 256 individual genetic variants with risk of total and aggressive prostate cancer. Likelihood ratio tests assessed multiplicative interaction. We observed interaction between ever use of fonofos and the PRS in relation to total and aggressive prostate cancer risk. Compared to the reference group (never use, PRS < median), men with ever use of fonofos and PRS > median had elevated risks of total (OR 1.35 [1.06-1.73], p-interaction = 0.03) and aggressive (OR 1.49 [1.09-2.04], p-interaction = 0.19) prostate cancer. There was also suggestion of interaction between pesticides and individual genetic variants occurring in regions associated with DNA damage response (CDH3, EMSY genes) and with variants related to altered androgen receptor-driven transcriptional programs critical for prostate cancer. Our study provides evidence that men with greater genetic susceptibility to prostate cancer may be at higher risk if they are also exposed to pesticides and suggests potential mechanisms by which pesticides may increase prostate cancer risk.

Genome-wide interaction analysis of folate for colorectal cancer risk

BACKGROUND

Epidemiological and experimental evidence suggests that higher folate intake is associated with decreased colorectal cancer (CRC) risk; however, the mechanisms underlying this relationship are not fully understood. Genetic variation that may have a direct or indirect impact on folate metabolism can provide insights into folate's role in CRC.

OBJECTIVES

Our aim was to perform a genome-wide interaction analysis to identify genetic variants that may modify the association of folate on CRC risk.

METHODS

We applied traditional case-control logistic regression, joint 3-degree of freedom, and a 2-step weighted hypothesis approach to test the interactions of common variants (allele frequency >1%) across the genome and dietary folate, folic acid supplement use, and total folate in relation to risk of CRC in 30,550 cases and 42,336 controls from 51 studies from 3 genetic consortia (CCFR, CORECT, GECCO).

RESULTS

Inverse associations of dietary, total folate, and folic acid supplement with CRC were found (odds ratio [OR]: 0.93; 95% confidence interval [CI]: 0.90, 0.96; and 0.91; 95% CI: 0.89, 0.94 per quartile higher intake, and 0.82 (95% CI: 0.78, 0.88) for users compared with nonusers, respectively). Interactions (P-interaction < 5×10-8) of folic acid supplement and variants in the 3p25.2 locus (in the region of Synapsin II [SYN2]/tissue inhibitor of metalloproteinase 4 [TIMP4]) were found using traditional interaction analysis, with variant rs150924902 (located upstream to SYN2) showing the strongest interaction. In stratified analyses by rs150924902 genotypes, folate supplementation was associated with decreased CRC risk among those carrying the TT genotype (OR: 0.82; 95% CI: 0.79, 0.86) but increased CRC risk among those carrying the TA genotype (OR: 1.63; 95% CI: 1.29, 2.05), suggesting a qualitative interaction (P-interaction = 1.4×10-8). No interactions were observed for dietary and total folate.

CONCLUSIONS

Variation in 3p25.2 locus may modify the association of folate supplement with CRC risk. Experimental studies and studies incorporating other relevant omics data are warranted to validate this finding.

Germline Sequencing Analysis to Inform Clinical Gene Panel Testing for Aggressive Prostate Cancer

Importance

Germline gene panel testing is recommended for men with advanced prostate cancer (PCa) or a family history of cancer. While evidence is limited for some genes currently included in panel testing, gene panels are also likely to be incomplete and missing genes that influence PCa risk and aggressive disease.

Objective

To identify genes associated with aggressive PCa.

Design, Setting, and Participants

A 2-stage exome sequencing case-only genetic association study was conducted including men of European ancestry from 18 international studies. Data analysis was performed from January 2021 to March 2023. Participants were 9185 men with aggressive PCa (including 6033 who died of PCa and 2397 with confirmed metastasis) and 8361 men with nonaggressive PCa.

Exposure

Sequencing data were evaluated exome-wide and in a focused investigation of 29 DNA repair pathway and cancer susceptibility genes, many of which are included on gene panels.

Main Outcomes and Measures

The primary study outcomes were aggressive (category T4 or both T3 and Gleason score ≥8 tumors, metastatic PCa, or PCa death) vs nonaggressive PCa (category T1 or T2 and Gleason score ≤6 tumors without known recurrence), and metastatic vs nonaggressive PCa.

Results

A total of 17 546 men of European ancestry were included in the analyses; mean (SD) age at diagnosis was 65.1 (9.2) years in patients with aggressive PCa and 63.7 (8.0) years in those with nonaggressive disease. The strongest evidence of association with aggressive or metastatic PCa was noted for rare deleterious variants in known PCa risk genes BRCA2 and ATM (P ≤ 1.9 × 10-6), followed by NBN (P = 1.7 × 10-4). This study found nominal evidence (P < .05) of association with rare deleterious variants in MSH2, XRCC2, and MRE11A. Five other genes had evidence of greater risk (OR≥2) but carrier frequency differences between aggressive and nonaggressive PCa were not statistically significant: TP53, RAD51D, BARD1, GEN1, and SLX4. Deleterious variants in these 11 candidate genes were carried by 2.3% of patients with nonaggressive, 5.6% with aggressive, and 7.0% with metastatic PCa.

Conclusions and Relevance

The findings of this study provide further support for DNA repair and cancer susceptibility genes to better inform disease management in men with PCa and for extending testing to men with nonaggressive disease, as men carrying deleterious alleles in these genes are likely to develop more advanced disease.

Antinuclear Antibodies Are Associated with an Increased Risk of Diffuse Large B-Cell Lymphoma

Immune dysregulation is thought to increase the risk of non-Hodgkin lymphoma (NHL), but the evidence varies by subtype. We evaluated whether antinuclear antibodies (ANA), double-stranded DNA antibodies (anti-dsDNA), and extractable nuclear antigen antibodies (anti-ENA) were associated with the risk of common NHL subtypes in a nested case-control study. The autoantibodies were tested in serum collected years prior to NHL diagnosis in 832 cases and 809 controls from the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. Logistic regression was used to determine odds ratios (ORs) and 95% confidence intervals (95% CI) for the association with NHL risk. No association was observed between ANA positivity and NHL risk overall (OR: 1.18, 95% CI: 0.88-1.58); however, ANA positivity was associated with an increased risk of diffuse large B-cell lymphoma (DLBCL) (OR: 1.83, 95% CI: 1.15-2.91), with 19.7% of cases and 12.2% of controls testing positive. The presence of either anti-ENA or anti-dsDNA was associated with an increased risk of NHL (OR: 2.93, 95% CI: 1.18-7.28), particularly DLBCL (OR: 3.51, 95% CI: 1.02-12.0) and marginal zone lymphoma (OR: 8.86, 95% CI: 1.26-62.0). Our study demonstrates that autoantibodies are associated with an elevated risk of DLBCL, providing support for autoimmunity as a risk factor.

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.

Heterogeneous genetic architectures and evolutionary genomics of prostate cancer in Sub-Saharan Africa

Men of African descent have the highest prostate cancer (CaP) incidence and mortality rates, yet the genetic basis of CaP in African men has been understudied. We used genomic data from 3,963 CaP cases and 3,509 controls recruited in Ghana, Nigeria, Senegal, South Africa, and Uganda, to infer ancestry-specific genetic architectures and fine-mapped disease associations. Fifteen independent associations at 8q24.21, 6q22.1, and 11q13.3 reached genome-wide significance, including four novel associations. Intriguingly, multiple lead SNPs are private alleles, a pattern arising from recent mutations and the out-of-Africa bottleneck. These African-specific alleles contribute to haplotypes with odds ratios above 2.4. We found that the genetic architecture of CaP differs across Africa, with effect size differences contributing more to this heterogeneity than allele frequency differences. Population genetic analyses reveal that African CaP associations are largely governed by neutral evolution. Collectively, our findings emphasize the utility of conducting genetic studies that use diverse populations.

Combining Asian and European genome-wide association studies of colorectal cancer improves risk prediction across racial and ethnic populations

Polygenic risk scores (PRS) have great potential to guide precision colorectal cancer (CRC) prevention by identifying those at higher risk to undertake targeted screening. However, current PRS using European ancestry data have sub-optimal performance in non-European ancestry populations, limiting their utility among these populations. Towards addressing this deficiency, we expand PRS development for CRC by incorporating Asian ancestry data (21,731 cases; 47,444 controls) into European ancestry training datasets (78,473 cases; 107,143 controls). The AUC estimates (95% CI) of PRS are 0.63(0.62-0.64), 0.59(0.57-0.61), 0.62(0.60-0.63), and 0.65(0.63-0.66) in independent datasets including 1681-3651 cases and 8696-115,105 controls of Asian, Black/African American, Latinx/Hispanic, and non-Hispanic White, respectively. They are significantly better than the European-centric PRS in all four major US racial and ethnic groups (p-values < 0.05). Further inclusion of non-European ancestry populations, especially Black/African American and Latinx/Hispanic, is needed to improve the risk prediction and enhance equity in applying PRS in clinical practice.

A metabolomic investigation of serum perfluorooctane sulfonate and perfluorooctanoate

BACKGROUND

Exposures to perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA), environmentally persistent chemicals detectable in the blood of most Americans, have been associated with several health outcomes. To offer insight into their possible biologic effects, we evaluated the metabolomic correlates of circulating PFOS and PFOA among 3,647 participants in eight nested case-control serum metabolomic profiling studies from the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial.

METHODS

Metabolomic profiling was conducted by Metabolon Inc., using ultra high-performance liquid chromatography/tandem accurate mass spectrometry. We conducted study-specific multivariable linear regression analyses estimating the associations of metabolite levels with levels of PFOS or PFOA. For metabolites measured in at least 3 of 8 nested case-control studies, random effects meta-analysis was used to summarize study-specific results (1,038 metabolites in PFOS analyses and 1,100 in PFOA analyses).

RESULTS

The meta-analysis identified 51 and 38 metabolites associated with PFOS and PFOA, respectively, at a Bonferroni-corrected significance level (4.8x10-5 and 4.6x10-5, respectively). For both PFOS and PFOA, the most common types of associated metabolites were lipids (sphingolipids, fatty acid metabolites) and xenobiotics (xanthine metabolites, chemicals). Positive associations were commonly observed with lipid metabolites sphingomyelin (d18:1/18:0) (P = 2.0x10-10 and 2.0x10-8, respectively), 3-carboxy-4-methyl-5-pentyl-2-furanpropionate (P = 2.7x10-15, 1.1x10-17), and lignoceroylcarnitine (C24) (P = 2.6x10-8, 6.2x10-6). The strongest positive associations were observed for chemicals 3,5-dichloro-2,6-dihydroxybenzoic acid (P = 3.0x10-112 and 6.8x10-13, respectively) and 3-bromo-5-chloro-2,6-dihydroxybenzoic acid (P = 1.6x10-14, 2.3x10-6). Other metabolites positively associated with PFOS included D-glucose (carbohydrate), carotene diol (vitamin A metabolism), and L-alpha-aminobutyric acid (glutathione metabolism), while uric acid (purine metabolite) was positively associated with PFOA. PFOS associations were consistent even after adjusting for PFOA as a covariate, while PFOA associations were greatly attenuated with PFOS adjustment.

CONCLUSIONS

In this large metabolomic study, we observed robust positive associations with PFOS for several molecules. Further investigation of these metabolites may offer insight into PFOS-related biologic effects.

Genetic Susceptibility to Nonalcoholic Fatty Liver Disease and Risk for Pancreatic Cancer: Mendelian Randomization

BACKGROUND

There are conflicting data on whether nonalcoholic fatty liver disease (NAFLD) is associated with susceptibility to pancreatic cancer. Using Mendelian randomization (MR), we investigated the relationship between genetic predisposition to NAFLD and risk for pancreatic cancer.

METHODS

Data from genome-wide association studies (GWAS) within the Pancreatic Cancer Cohort Consortium (PanScan; cases n = 5,090, controls n = 8,733) and the Pancreatic Cancer Case Control Consortium (PanC4; cases n = 4,163, controls n = 3,792) were analyzed. We used data on 68 genetic variants with four different MR methods [inverse variance weighting (IVW), MR-Egger, simple median, and penalized weighted median] separately to predict genetic heritability of NAFLD. We then assessed the relationship between each of the four MR methods and pancreatic cancer risk, using logistic regression to calculate ORs and 95% confidence intervals (CI), adjusting for PC risk factors, including obesity and diabetes.

RESULTS

No association was found between genetically predicted NAFLD and pancreatic cancer risk in the PanScan or PanC4 samples [e.g., PanScan, IVW OR, 1.04; 95% confidence interval (CI), 0.88-1.22; MR-Egger OR, 0.89; 95% CI, 0.65-1.21; PanC4, IVW OR, 1.07; 95% CI, 0.90-1.27; MR-Egger OR, 0.93; 95% CI, 0.67-1.28]. None of the four MR methods indicated an association between genetically predicted NAFLD and pancreatic cancer risk in either sample.

CONCLUSIONS

Genetic predisposition to NAFLD is not associated with pancreatic cancer risk.

IMPACT

Given the close relationship between NAFLD and metabolic conditions, it is plausible that any association between NAFLD and pancreatic cancer might reflect host metabolic perturbations (e.g., obesity, diabetes, or metabolic syndrome) and does not necessarily reflect a causal relationship between NAFLD and pancreatic cancer.

Genetically inferred birthweight, height, and puberty timing and risk of osteosarcoma

INTRODUCTION

Several studies have linked increased risk of osteosarcoma with tall stature, high birthweight, and early puberty, although evidence is inconsistent. We used genetic risk scores (GRS) based on established genetic loci for these traits and evaluated associations between genetically inferred birthweight, height, and puberty timing with osteosarcoma.

METHODS

Using genotype data from two genome-wide association studies, totaling 1039 cases and 2923 controls of European ancestry, association analyses were conducted using logistic regression for each study and meta-analyzed to estimate pooled odds ratios (ORs) and 95% confidence intervals (CIs). Subgroup analyses were conducted by case diagnosis age, metastasis status, tumor location, tumor histology, and presence of a known pathogenic variant in a cancer susceptibility gene.

RESULTS

Genetically inferred higher birthweight was associated with an increased risk of osteosarcoma (OR =1.59, 95% CI 1.07-2.38, P = 0.02). This association was strongest in cases without metastatic disease (OR =2.46, 95% CI 1.44-4.19, P = 9.5 ×10-04). Although there was no overall association between osteosarcoma and genetically inferred taller stature (OR=1.06, 95% CI 0.96-1.17, P = 0.28), the GRS for taller stature was associated with an increased risk of osteosarcoma in 154 cases with a known pathogenic cancer susceptibility gene variant (OR=1.29, 95% CI 1.03-1.63, P = 0.03). There were no significant associations between the GRS for puberty timing and osteosarcoma.

CONCLUSION

A genetic propensity to higher birthweight was associated with increased osteosarcoma risk, suggesting that shared genetic factors or biological pathways that affect birthweight may contribute to osteosarcoma pathogenesis.

Probing the diabetes and colorectal cancer relationship using gene - environment interaction analyses

BACKGROUND

Diabetes is an established risk factor for colorectal cancer. However, the mechanisms underlying this relationship still require investigation and it is not known if the association is modified by genetic variants. To address these questions, we undertook a genome-wide gene-environment interaction analysis.

METHODS

We used data from 3 genetic consortia (CCFR, CORECT, GECCO; 31,318 colorectal cancer cases/41,499 controls) and undertook genome-wide gene-environment interaction analyses with colorectal cancer risk, including interaction tests of genetics(G)xdiabetes (1-degree of freedom; d.f.) and joint testing of Gxdiabetes, G-colorectal cancer association (2-d.f. joint test) and G-diabetes correlation (3-d.f. joint test).

RESULTS

Based on the joint tests, we found that the association of diabetes with colorectal cancer risk is modified by loci on chromosomes 8q24.11 (rs3802177, SLC30A8 - ORAA: 1.62, 95% CI: 1.34-1.96; ORAG: 1.41, 95% CI: 1.30-1.54; ORGG: 1.22, 95% CI: 1.13-1.31; p-value3-d.f.: 5.46 × 10-11) and 13q14.13 (rs9526201, LRCH1 - ORGG: 2.11, 95% CI: 1.56-2.83; ORGA: 1.52, 95% CI: 1.38-1.68; ORAA: 1.13, 95% CI: 1.06-1.21; p-value2-d.f.: 7.84 × 10-09).

DISCUSSION

These results suggest that variation in genes related to insulin signaling (SLC30A8) and immune function (LRCH1) may modify the association of diabetes with colorectal cancer risk and provide novel insights into the biology underlying the diabetes and colorectal cancer relationship.

A Genetic Locus within the FMN1/GREM1 Gene Region Interacts with Body Mass Index in Colorectal Cancer Risk

Colorectal cancer risk can be impacted by genetic, environmental, and lifestyle factors, including diet and obesity. Gene-environment interactions (G × E) can provide biological insights into the effects of obesity on colorectal cancer risk. Here, we assessed potential genome-wide G × E interactions between body mass index (BMI) and common SNPs for colorectal cancer risk using data from 36,415 colorectal cancer cases and 48,451 controls from three international colorectal cancer consortia (CCFR, CORECT, and GECCO). The G × E tests included the conventional logistic regression using multiplicative terms (one degree of freedom, 1DF test), the two-step EDGE method, and the joint 3DF test, each of which is powerful for detecting G × E interactions under specific conditions. BMI was associated with higher colorectal cancer risk. The two-step approach revealed a statistically significant G×BMI interaction located within the Formin 1/Gremlin 1 (FMN1/GREM1) gene region (rs58349661). This SNP was also identified by the 3DF test, with a suggestive statistical significance in the 1DF test. Among participants with the CC genotype of rs58349661, overweight and obesity categories were associated with higher colorectal cancer risk, whereas null associations were observed across BMI categories in those with the TT genotype. Using data from three large international consortia, this study discovered a locus in the FMN1/GREM1 gene region that interacts with BMI on the association with colorectal cancer risk. Further studies should examine the potential mechanisms through which this locus modifies the etiologic link between obesity and colorectal cancer.

SIGNIFICANCE

This gene-environment interaction analysis revealed a genetic locus in FMN1/GREM1 that interacts with body mass index in colorectal cancer risk, suggesting potential implications for precision prevention strategies.

Evaluating approaches for constructing polygenic risk scores for prostate cancer in men of African and European ancestry

Genome-wide polygenic risk scores (GW-PRSs) have been reported to have better predictive ability than PRSs based on genome-wide significance thresholds across numerous traits. We compared the predictive ability of several GW-PRS approaches to a recently developed PRS of 269 established prostate cancer-risk variants from multi-ancestry GWASs and fine-mapping studies (PRS269). GW-PRS models were trained with a large and diverse prostate cancer GWAS of 107,247 cases and 127,006 controls that we previously used to develop the multi-ancestry PRS269. Resulting models were independently tested in 1,586 cases and 1,047 controls of African ancestry from the California Uganda Study and 8,046 cases and 191,825 controls of European ancestry from the UK Biobank and further validated in 13,643 cases and 210,214 controls of European ancestry and 6,353 cases and 53,362 controls of African ancestry from the Million Veteran Program. In the testing data, the best performing GW-PRS approach had AUCs of 0.656 (95% CI = 0.635-0.677) in African and 0.844 (95% CI = 0.840-0.848) in European ancestry men and corresponding prostate cancer ORs of 1.83 (95% CI = 1.67-2.00) and 2.19 (95% CI = 2.14-2.25), respectively, for each SD unit increase in the GW-PRS. Compared to the GW-PRS, in African and European ancestry men, the PRS269 had larger or similar AUCs (AUC = 0.679, 95% CI = 0.659-0.700 and AUC = 0.845, 95% CI = 0.841-0.849, respectively) and comparable prostate cancer ORs (OR = 2.05, 95% CI = 1.87-2.26 and OR = 2.21, 95% CI = 2.16-2.26, respectively). Findings were similar in the validation studies. This investigation suggests that current GW-PRS approaches may not improve the ability to predict prostate cancer risk compared to the PRS269 developed from multi-ancestry GWASs and fine-mapping.

A prospective nested case-control study of serum concentrations of per- and polyfluoroalkyl substances and aggressive prostate cancer risk

Per- and polyfluoroalkyl substances (PFAS) are environmentally persistent organic pollutants detectable in the serum of most U.S. adults. Some studies of highly-exposed individuals have suggested an association between PFAS and prostate cancer, but evidence from population-based studies is limited. We investigated the association between pre-diagnostic serum PFAS concentrations and aggressive prostate cancer risk in a large prospective study. We measured pre-diagnostic serum concentrations of eight PFAS, including perfluorooctanoate (PFOA), for 750 aggressive prostate cancer cases and 750 individually matched controls within the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. We assessed the reproducibility of PFAS concentrations in serial samples collected up to six years apart among 60 controls using intraclass correlation coefficients (ICCs). Conditional logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs) for the association with prostate cancer, adjusting for other PFAS and potential confounders. Concentrations of most PFAS were consistent (ICC>0.7) across the serial samples over time. We observed an inverse association between PFOA and aggressive prostate cancer (OR_continuous = 0.79, 95% CI = 0.63, 0.99), but the association was limited to cases diagnosed ≤3 years after blood collection and became statistically non-significant for cases diagnosed with later follow-up (>3 years, OR_continuous = 0.90, 95% CI = 0.79, 1.03). Other PFAS were not associated with aggressive prostate cancer risk. Although we cannot rule out an increased risk at higher levels, our findings from a population with PFAS serum concentrations comparable to the general population do not support an association with increased risk of aggressive prostate cancer.

Evidence of Novel Susceptibility Variants for Prostate Cancer and a Multiancestry Polygenic Risk Score Associated with Aggressive Disease in Men of African Ancestry

BACKGROUND

Genetic factors play an important role in prostate cancer (PCa) susceptibility.

OBJECTIVE

To discover common genetic variants contributing to the risk of PCa in men of African ancestry.

DESIGN, SETTING, AND PARTICIPANTS

We conducted a meta-analysis of ten genome-wide association studies consisting of 19378 cases and 61620 controls of African ancestry.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Common genotyped and imputed variants were tested for their association with PCa risk. Novel susceptibility loci were identified and incorporated into a multiancestry polygenic risk score (PRS). The PRS was evaluated for associations with PCa risk and disease aggressiveness.

RESULTS AND LIMITATIONS

Nine novel susceptibility loci for PCa were identified, of which seven were only found or substantially more common in men of African ancestry, including an African-specific stop-gain variant in the prostate-specific gene anoctamin 7 (ANO7). A multiancestry PRS of 278 risk variants conferred strong associations with PCa risk in African ancestry studies (odds ratios [ORs] >3 and >5 for men in the top PRS decile and percentile, respectively). More importantly, compared with men in the 40-60% PRS category, men in the top PRS decile had a significantly higher risk of aggressive PCa (OR = 1.23, 95% confidence interval = 1.10-1.38, p = 4.4 × 10-4).

CONCLUSIONS

This study demonstrates the importance of large-scale genetic studies in men of African ancestry for a better understanding of PCa susceptibility in this high-risk population and suggests a potential clinical utility of PRS in differentiating between the risks of developing aggressive and nonaggressive disease in men of African ancestry.

PATIENT SUMMARY

In this large genetic study in men of African ancestry, we discovered nine novel prostate cancer (PCa) risk variants. We also showed that a multiancestry polygenic risk score was effective in stratifying PCa risk, and was able to differentiate risk of aggressive and nonaggressive disease.

Improve the model of disease subtype heterogeneity by leveraging external summary data

Researchers are often interested in understanding the disease subtype heterogeneity by testing whether a risk exposure has the same level of effect on different disease subtypes. The polytomous logistic regression (PLR) model provides a flexible tool for such an evaluation. Disease subtype heterogeneity can also be investigated with a case-only study that uses a case-case comparison procedure to directly assess the difference between risk effects on two disease subtypes. Motivated by a large consortium project on the genetic basis of non-Hodgkin lymphoma (NHL) subtypes, we develop PolyGIM, a procedure to fit the PLR model by integrating individual-level data with summary data extracted from multiple studies under different designs. The summary data consist of coefficient estimates from working logistic regression models established by external studies. Examples of the working model include the case-case comparison model and the case-control comparison model, which compares the control group with a subtype group or a broad disease group formed by merging several subtypes. PolyGIM efficiently evaluates risk effects and provides a powerful test for disease subtype heterogeneity in situations when only summary data, instead of individual-level data, is available from external studies due to various informatics and privacy constraints. We investigate the theoretic properties of PolyGIM and use simulation studies to demonstrate its advantages. Using data from eight genome-wide association studies within the NHL consortium, we apply it to study the effect of the polygenic risk score defined by a lymphoid malignancy on the risks of four NHL subtypes. These results show that PolyGIM can be a valuable tool for pooling data from multiple sources for a more coherent evaluation of disease subtype heterogeneity.

Increase in power by obtaining 10 or more controls per case when type-1 error is small in large-scale association studies

BACKGROUND

The rule of thumb that there is little gain in statistical power by obtaining more than 4 controls per case, is based on type-1 error α = 0.05. However, association studies that evaluate thousands or millions of associations use smaller α and may have access to plentiful controls. We investigate power gains, and reductions in p-values, when increasing well beyond 4 controls per case, for small α.

METHODS

We calculate the power, the median expected p-value, and the minimum detectable odds-ratio (OR), as a function of the number of controls/case, as α decreases.

RESULTS

As α decreases, at each ratio of controls per case, the increase in power is larger than for α = 0.05. For α between 10-6 and 10-9 (typical for thousands or millions of associations), increasing from 4 controls per case to 10-50 controls per case increases power. For example, a study with power = 0.2 (α = 5 × 10-8) with 1 control/case has power = 0.65 with 4 controls/case, but with 10 controls/case has power = 0.78, and with 50 controls/case has power = 0.84. For situations where obtaining more than 4 controls per case provides small increases in power beyond 0.9 (at small α), the expected p-value can decrease by orders-of-magnitude below α. Increasing from 1 to 4 controls/case reduces the minimum detectable OR toward the null by 20.9%, and from 4 to 50 controls/case reduces by an additional 9.7%, a result which applies regardless of α and hence also applies to "regular" α = 0.05 epidemiology.

CONCLUSIONS

At small α, versus 4 controls/case, recruiting 10 or more controls/cases can increase power, reduce the expected p-value by 1-2 orders of magnitude, and meaningfully reduce the minimum detectable OR. These benefits of increasing the controls/case ratio increase as the number of cases increases, although the amount of benefit depends on exposure frequencies and true OR. Provided that controls are comparable to cases, our findings suggest greater sharing of comparable controls in large-scale association studies.

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.

Coffee Intake, Caffeine Metabolism Genotype, and Survival Among Men with Prostate Cancer

BACKGROUND

Coffee intake may lower prostate cancer risk and progression, but postdiagnosis outcomes by caffeine metabolism genotype are not well characterized.

OBJECTIVE

To evaluate associations between coffee intake, caffeine metabolism genotype, and survival in a large, multicenter study of men with prostate cancer.

DESIGN, SETTING, AND PARTICIPANTS

Data from The PRACTICAL Consortium database for 5727 men with prostate cancer from seven US, Australian, and European studies were included. The cases included had data available for the CYP1A2 -163C>A rs762551 single-nucleotide variant associated with caffeine metabolism, coffee intake, and >6 mo of follow-up.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Multivariable-adjusted Cox proportional hazards models across pooled patient-level data were used to compare the effect of coffee intake (categorized as low [reference], high, or none/very low) in relation to overall survival (OS) and prostate cancer-specific survival (PCSS), with stratified analyses conducted by clinical disease risk and genotype.

RESULTS AND LIMITATIONS

High coffee intake appeared to be associated with longer PCSS (hazard ratio [HR] 0.85, 95% confidence interval [CI] 0.68-1.08; p = 0.18) and OS (HR 0.90, 95% CI 0.77-1.07; p = 0.24), although results were not statistically significant. In the group with clinically localized disease, high coffee intake was associated with longer PCSS (HR 0.66, 95% CI 0.44-0.98; p = 0.040), with comparable results for the group with advanced disease (HR 0.92, 95% CI 0.69-1.23; p = 0.6). High coffee intake was associated with longer PCSS among men with the CYP1A2 AA (HR 0.67, 95% CI 0.49-0.93; p = 0.017) but not the AC/CC genotype (p = 0.8); an interaction was detected (p = 0.042). No associations with OS were observed in subgroup analyses (p > 0.05). Limitations include the nominal statistical significance and residual confounding.

CONCLUSIONS

Coffee intake was associated with longer PCSS among men with a CYP1A2 -163AA (*1F/*1F) genotype, a finding that will require further replication.

PATIENT SUMMARY

It is likely that coffee intake is associated with longer prostate cancer-specific survival in certain groups, but more research is needed to fully understand which men may benefit and why.

Genome-wide association study of lung adenocarcinoma in East Asia and comparison with a European population

Lung adenocarcinoma is the most common type of lung cancer. Known risk variants explain only a small fraction of lung adenocarcinoma heritability. Here, we conducted a two-stage genome-wide association study of lung adenocarcinoma of East Asian ancestry (21,658 cases and 150,676 controls; 54.5% never-smokers) and identified 12 novel susceptibility variants, bringing the total number to 28 at 25 independent loci. Transcriptome-wide association analyses together with colocalization studies using a Taiwanese lung expression quantitative trait loci dataset (n = 115) identified novel candidate genes, including FADS1 at 11q12 and ELF5 at 11p13. In a multi-ancestry meta-analysis of East Asian and European studies, four loci were identified at 2p11, 4q32, 16q23, and 18q12. At the same time, most of our findings in East Asian populations showed no evidence of association in European populations. In our studies drawn from East Asian populations, a polygenic risk score based on the 25 loci had a stronger association in never-smokers vs. individuals with a history of smoking (Pinteraction = 0.0058). These findings provide new insights into the etiology of lung adenocarcinoma in individuals from East Asian populations, which could be important in developing translational applications.

Evaluating Approaches for Constructing Polygenic Risk Scores for Prostate Cancer in Men of African and European Ancestry

Genome-wide polygenic risk scores (GW-PRS) have been reported to have better predictive ability than PRS based on genome-wide significance thresholds across numerous traits. We compared the predictive ability of several GW-PRS approaches to a recently developed PRS of 269 established prostate cancer risk variants from multi-ancestry GWAS and fine-mapping studies (PRS 269 ). GW-PRS models were trained using a large and diverse prostate cancer GWAS of 107,247 cases and 127,006 controls used to develop the multi-ancestry PRS 269 . Resulting models were independently tested in 1,586 cases and 1,047 controls of African ancestry from the California/Uganda Study and 8,046 cases and 191,825 controls of European ancestry from the UK Biobank and further validated in 13,643 cases and 210,214 controls of European ancestry and 6,353 cases and 53,362 controls of African ancestry from the Million Veteran Program. In the testing data, the best performing GW-PRS approach had AUCs of 0.656 (95% CI=0.635-0.677) in African and 0.844 (95% CI=0.840-0.848) in European ancestry men and corresponding prostate cancer OR of 1.83 (95% CI=1.67-2.00) and 2.19 (95% CI=2.14-2.25), respectively, for each SD unit increase in the GW-PRS. However, compared to the GW-PRS, in African and European ancestry men, the PRS 269 had larger or similar AUCs (AUC=0.679, 95% CI=0.659-0.700 and AUC=0.845, 95% CI=0.841-0.849, respectively) and comparable prostate cancer OR (OR=2.05, 95% CI=1.87-2.26 and OR=2.21, 95% CI=2.16-2.26, respectively). Findings were similar in the validation data. This investigation suggests that current GW-PRS approaches may not improve the ability to predict prostate cancer risk compared to the multi-ancestry PRS 269 constructed with fine-mapping.

Polymorphisms within Autophagy-Related Genes as Susceptibility Biomarkers for Multiple Myeloma: A Meta-Analysis of Three Large Cohorts and Functional Characterization

Multiple myeloma (MM) arises following malignant proliferation of plasma cells in the bone marrow, that secrete high amounts of specific monoclonal immunoglobulins or light chains, resulting in the massive production of unfolded or misfolded proteins. Autophagy can have a dual role in tumorigenesis, by eliminating these abnormal proteins to avoid cancer development, but also ensuring MM cell survival and promoting resistance to treatments. To date no studies have determined the impact of genetic variation in autophagy-related genes on MM risk. We performed meta-analysis of germline genetic data on 234 autophagy-related genes from three independent study populations including 13,387 subjects of European ancestry (6863 MM patients and 6524 controls) and examined correlations of statistically significant single nucleotide polymorphisms (SNPs; p < 1 × 10^-9) with immune responses in whole blood, peripheral blood mononuclear cells (PBMCs), and monocyte-derived macrophages (MDM) from a large population of healthy donors from the Human Functional Genomic Project (HFGP). We identified SNPs in six loci, CD46, IKBKE, PARK2, ULK4, ATG5, and CDKN2A associated with MM risk (p = 4.47 × 10^-4-5.79 × 10^-14). Mechanistically, we found that the ULK4_rs6599175 SNP correlated with circulating concentrations of vitamin D3 (p = 4.0 × 10^-4), whereas the IKBKE_rs17433804 SNP correlated with the number of transitional CD24⁺CD38⁺ B cells (p = 4.8 × 10^-4) and circulating serum concentrations of Monocyte Chemoattractant Protein (MCP)-2 (p = 3.6 × 10^-4). We also found that the CD46_rs1142469 SNP correlated with numbers of CD19⁺ B cells, CD19⁺CD3^- B cells, CD5⁺IgD^- cells, IgM^- cells, IgD^-IgM^- cells, and CD4^-CD8^- PBMCs (p = 4.9 × 10^-4-8.6 × 10^-4) and circulating concentrations of interleukin (IL)-20 (p = 0.00082). Finally, we observed that the CDKN2A_rs2811710 SNP correlated with levels of CD4⁺EMCD45RO⁺CD27^- cells (p = 9.3 × 10^-4). These results suggest that genetic variants within these six loci influence MM risk through the modulation of specific subsets of immune cells, as well as vitamin D3^-, MCP-2^-, and IL20-dependent pathways.

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.

Do GWAS-Identified Risk Variants for Chronic Lymphocytic Leukemia Influence Overall Patient Survival and Disease Progression?

Chronic lymphocytic leukemia (CLL) is the most common leukemia among adults worldwide. Although genome-wide association studies (GWAS) have uncovered the germline genetic component underlying CLL susceptibility, the potential use of GWAS-identified risk variants to predict disease progression and patient survival remains unexplored. Here, we evaluated whether 41 GWAS-identified risk variants for CLL could influence overall survival (OS) and disease progression, defined as time to first treatment (TTFT) in a cohort of 1039 CLL cases ascertained through the CRuCIAL consortium. Although this is the largest study assessing the effect of GWAS-identified susceptibility variants for CLL on OS, we only found a weak association of ten single nucleotide polymorphisms (SNPs) with OS (p < 0.05) that did not remain significant after correction for multiple testing. In line with these results, polygenic risk scores (PRSs) built with these SNPs in the CRuCIAL cohort showed a modest association with OS and a low capacity to predict patient survival, with an area under the receiver operating characteristic curve (AUROC) of 0.57. Similarly, seven SNPs were associated with TTFT (p < 0.05); however, these did not reach the multiple testing significance threshold, and the meta-analysis with previous published data did not confirm any of the associations. As expected, PRSs built with these SNPs showed reduced accuracy in prediction of disease progression (AUROC = 0.62). These results suggest that susceptibility variants for CLL do not impact overall survival and disease progression in CLL patients.

Abstract 3508: Association between clonal hematopoiesis and risk of prostate cancer in a large sample of African ancestry men

Clonal hematopoiesis of indeterminate potential (CHIP) has been associated with inflammation, which is a risk factor for cancer, including prostate cancer. We previously reported weak evidence of an association between CHIP and prostate cancer risk in men of European ancestry. However, little is known for African ancestry populations. We investigated the association of age-related CHIP with overall and aggressive prostate cancer risk in a large whole-exome sequencing study of 12,049 African ancestry men, including 7,176 prostate cancer cases (of which 3,283 had aggressive disease and 1,074 had metastatic disease) and 4,873 controls. Somatic variant calling was carried out using GATK Mutect2, and only variants with minor allele frequencies (MAF) &lt;0.1% and a variant allelic fraction (VAF) &gt;5% were included. Variants with a MAF ≥0.1% in gnomAD were excluded. CHIP variants were identified from a list of pre-specified mutations in 74 genes. Associations were tested using regression models adjusting for age, sub-study, and top 10 principal components, with statistical significance tested by the likelihood ratio test and applying a Bonferroni correction to account for multiple testing. In total, 998 variants in 57 CHIP genes were identified. Consistent with previous results, we observed a strong association between CHIP and age at blood draw. CHIP genes in aggregate were not statistically significantly associated with risks of total (OR=1.12, 95% CI=0.97-1.28), aggressive (OR=1.14, 95% CI=0.92-1.43) or metastatic (OR=1.17, 95% CI=0.91-1.49) prostate cancer. We observed that carriers of variants in DNMT3A, which is the gene that harbors the most CHIP driver mutations, had a nominally elevated risk of prostate cancer compared to non-carriers (OR=1.35, 95% CI=1.08-1.68, p=0.007). Additionally, carriers of variants in EZH2, which is implicated in cancer progression, showed a suggestive association with aggressive prostate cancer (OR=7.33, 95% CI=1.01-53.21, p=0.029). After adjusting for age at blood draw, CHIP genes in aggregate were not associated with age at prostate cancer diagnosis. However, we found that EZH2 variants carriers were diagnosed 12.9 years earlier on average than non-carriers (95% CI=6.1-19.7, adjusted p=0.01). A prostate cancer polygenic risk score (PRS) constructed using 269 risk variants was not associated with CHIP carrier status in aggregate (OR=0.99, 95% CI=0.92-1.06, p=0.70) or with any individual gene (all adjusted p&gt;0.05). In summary, overall CHIP is not likely to be a risk factor of prostate cancer or aggressive disease in men of African ancestry. However, our results do confirm the association of CHIP in DNMT3A with prostate cancer risk as reported in previous studies in men of European ancestry. Future work will be needed to evaluate the biological causality of DNMT3A- and EZH2- related CHIP on prostate cancer.

Citation Format: Anqi Wang, Yili Xu, Xin Sheng, Raymond Hughley, Ben Adusei, Mohamed Jalloh, Serigne Magueye Gueye, Andrew A Adjei, James Mensah, Pedro W. Fernandez, Akin Olupelumi Adebiyi, Oseremen Inokhoife Aisuodionoe-Shadrach, Lindsay Petersen, Maureen Joffe, Jeannette T. Bensen, James L. Mohler, Jack A. Taylor, Eboneé N. Butler, Sue A. Ingles, Benjamin A. Rybicki, Janet L. Stanford, Wei Zheng, Sonja I. Berndt, Chad D. Huff, Joseph Lachance, Luc Multigner, Caroline Andrews, Timothy R. Rebbeck, Laurent Brureau, Stephen J. Chanock, Adam de Smith, Fei Chen, Burcu F. Darst, David V. Conti, Christopher A. Haiman. Association between clonal hematopoiesis and risk of prostate cancer in a large sample of African ancestry men [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3508.

Biochemical activity induced by a germline variation in KLK3 (PSA) associates with cellular function and clinical outcome in prostate cancer

Genetic variation at the 19q13.3 KLK locus is linked with prostate cancer susceptibility. The non-synonymous KLK3 SNP, rs17632542 (c.536T>C; Ile163Thr-substitution in PSA) is associated with reduced prostate cancer risk, however, the functional relevance is unknown. Here, we identify that the SNP variant-induced change in PSA biochemical activity as a previously undescribed function mediating prostate cancer pathogenesis. The 'Thr' PSA variant led to small subcutaneous tumours, supporting reduced prostate cancer risk. However, 'Thr' PSA also displayed higher metastatic potential with pronounced osteolytic activity in an experimental metastasis in-vivo model. Biochemical characterization of this PSA variant demonstrated markedly reduced proteolytic activity that correlated with differences in in-vivo tumour burden. The SNP is associated with increased risk for aggressive disease and prostate cancer-specific mortality in three independent cohorts, highlighting its critical function in mediating metastasis. Carriers of this SNP allele had reduced serum total PSA and a higher free/total PSA ratio that could contribute to late biopsy decisions and delay in diagnosis. Our results provide a molecular explanation for the prominent 19q13.3 KLK locus, rs17632542 SNP, association with a spectrum of prostate cancer clinical outcomes.

Genome-wide Interaction Study with Smoking for Colorectal Cancer Risk Identifies Novel Genetic Loci Related to Tumor Suppression, Inflammation, and Immune Response

BACKGROUND

Tobacco smoking is an established risk factor for colorectal cancer. However, genetically defined population subgroups may have increased susceptibility to smoking-related effects on colorectal cancer.

METHODS

A genome-wide interaction scan was performed including 33,756 colorectal cancer cases and 44,346 controls from three genetic consortia.

RESULTS

Evidence of an interaction was observed between smoking status (ever vs. never smokers) and a locus on 3p12.1 (rs9880919, P = 4.58 × 10-8), with higher associated risk in subjects carrying the GG genotype [OR, 1.25; 95% confidence interval (CI), 1.20-1.30] compared with the other genotypes (OR <1.17 for GA and AA). Among ever smokers, we observed interactions between smoking intensity (increase in 10 cigarettes smoked per day) and two loci on 6p21.33 (rs4151657, P = 1.72 × 10-8) and 8q24.23 (rs7005722, P = 2.88 × 10-8). Subjects carrying the rs4151657 TT genotype showed higher risk (OR, 1.12; 95% CI, 1.09-1.16) compared with the other genotypes (OR <1.06 for TC and CC). Similarly, higher risk was observed among subjects carrying the rs7005722 AA genotype (OR, 1.17; 95% CI, 1.07-1.28) compared with the other genotypes (OR <1.13 for AC and CC). Functional annotation revealed that SNPs in 3p12.1 and 6p21.33 loci were located in regulatory regions, and were associated with expression levels of nearby genes. Genetic models predicting gene expression revealed that smoking parameters were associated with lower colorectal cancer risk with higher expression levels of CADM2 (3p12.1) and ATF6B (6p21.33).

CONCLUSIONS

Our study identified novel genetic loci that may modulate the risk for colorectal cancer of smoking status and intensity, linked to tumor suppression and immune response.

IMPACT

These findings can guide potential prevention treatments.

Prognostic role of detailed colorectal location and tumor molecular features: analyses of 13,101 colorectal cancer patients including 2994 early-onset cases

BACKGROUND

The pathogenic effect of colorectal tumor molecular features may be influenced by several factors, including those related to microbiota, inflammation, metabolism, and epigenetics, which may change along colorectal segments. We hypothesized that the prognostic association of colon cancer location might differ by tumor molecular characteristics.

METHODS

Utilizing a consortium dataset of 13,101 colorectal cancer cases, including 2994 early-onset cases, we conducted survival analyses of detailed tumor location stratified by statuses of microsatellite instability (MSI), CpG island methylator phenotype (CIMP), and KRAS and BRAF oncogenic mutation.

RESULTS

There was a statistically significant trend for better colon cancer-specific survival in relation to tumor location from the cecum to sigmoid colon (Ptrend = 0.002), excluding the rectum. The prognostic association of colon location differed by MSI status (Pinteraction = 0.001). Non-MSI-high tumors exhibited the cecum-to-sigmoid trend for better colon cancer-specific survival [Ptrend < 0.001; multivariable hazard ratio (HR) for the sigmoid colon (vs. cecum), 0.80; 95% confidence interval (CI) 0.70-0.92], whereas MSI-high tumors demonstrated a suggestive cecum-to-sigmoid trend for worse survival (Ptrend = 0.020; the corresponding HR, 2.13; 95% CI 1.15-3.92). The prognostic association of colon tumor location also differed by CIMP status (Pinteraction = 0.003) but not significantly by age, stage, or other features. Furthermore, MSI-high status was a favorable prognostic indicator in all stages.

CONCLUSIONS

Both detailed colonic location and tumor molecular features need to be accounted for colon cancer prognostication to advance precision medicine. Our study indicates the important role of large-scale studies to robustly examine detailed colonic subsites in molecular oncology research.

Serum perfluorooctane sulfonate and perfluorooctanoate and risk of postmenopausal breast cancer according to hormone receptor status: An analysis in the Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial

Per- and polyfluoroalkyl substances (PFAS) are highly persistent endocrine-disrupting chemicals that may contribute to breast cancer development; however, epidemiologic evidence is limited. We investigated associations between prediagnostic serum levels of perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA) and postmenopausal breast cancer risk, overall and by hormone receptor status, in a nested case-control study of 621 cases and 621 matched controls in the Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial. PFOS and PFOA levels were determined based on serum metabolomic profiling performed using ultraperformance liquid chromatography-tandem mass spectrometry. We used multivariable conditional logistic regression to estimate odds ratios (ORs) and 95% confidence intervals (CIs) for the association between each PFAS and breast cancer risk, overall, by estrogen receptor (ER) or progesterone receptor (PR) status, and by joint ER/PR status. We found little evidence of association between PFOS or PFOA and breast cancer risk overall. However, in subtype-specific analyses, we observed statistically significant increased risks of ER+, PR+, and ER+/PR+ tumors for the third vs lowest quartile of serum PFOS (ORs [95% CIs] = 1.59 [1.01-2.50], 2.34 [1.29-4.23], and 2.19 [1.21-3.98], respectively) and elevated but nonstatistically significant ORs for the fourth quartile. Conversely, for PFOA, modest positive associations with ER-, PR-, ER+/PR-, and ER-/PR- tumors were generally seen in the upper quartiles. Our findings contribute evidence supporting positive associations between serum PFOS and hormone receptor-positive tumors, and possibly between PFOA and receptor-negative tumors. Future prospective studies incorporating tumor hormone receptor status are needed to better understand the role of PFAS in breast cancer etiology.

Body mass index and molecular subtypes of colorectal cancer

BACKGROUND

Obesity is an established risk factor for colorectal cancer (CRC), but the evidence for the association is inconsistent across molecular subtypes of the disease.

METHODS

We pooled data on body mass index (BMI), tumor microsatellite instability status, CpG island methylator phenotype status, BRAF and KRAS mutations, and Jass classification types for 11 872 CRC cases and 11 013 controls from 11 observational studies. We used multinomial logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CI) adjusted for covariables.

RESULTS

Higher BMI was associated with increased CRC risk (OR per 5 kg/m2 = 1.18, 95% CI = 1.15 to 1.22). The positive association was stronger for men than women but similar across tumor subtypes defined by individual molecular markers. In analyses by Jass type, higher BMI was associated with elevated CRC risk for types 1-4 cases but not for type 5 CRC cases (considered familial-like/Lynch syndrome microsatellite instability-H, CpG island methylator phenotype-low or negative, BRAF-wild type, KRAS-wild type, OR = 1.04, 95% CI = 0.90 to 1.20). This pattern of associations for BMI and Jass types was consistent by sex and design of contributing studies (cohort or case-control).

CONCLUSIONS

In contrast to previous reports with fewer study participants, we found limited evidence of heterogeneity for the association between BMI and CRC risk according to molecular subtype, suggesting that obesity influences nearly all major pathways involved in colorectal carcinogenesis. The null association observed for the Jass type 5 suggests that BMI is not a risk factor for the development of CRC for individuals with Lynch syndrome.

Combining Asian-European Genome-Wide Association Studies of Colorectal Cancer Improves Risk Prediction Across Race and Ethnicity

Polygenic risk scores (PRS) have great potential to guide precision colorectal cancer (CRC) prevention by identifying those at higher risk to undertake targeted screening. However, current PRS using European ancestry data have sub-optimal performance in non-European ancestry populations, limiting their utility among these populations. Towards addressing this deficiency, we expanded PRS development for CRC by incorporating Asian ancestry data (21,731 cases; 47,444 controls) into European ancestry training datasets (78,473 cases; 107,143 controls). The AUC estimates (95% CI) of PRS were 0.63(0.62-0.64), 0.59(0.57-0.61), 0.62(0.60-0.63), and 0.65(0.63-0.66) in independent datasets including 1,681-3,651 cases and 8,696-115,105 controls of Asian, Black/African American, Latinx/Hispanic, and non-Hispanic White, respectively. They were significantly better than the European-centric PRS in all four major US racial and ethnic groups (p-values<0.05). Further inclusion of non-European ancestry populations, especially Black/African American and Latinx/Hispanic, is needed to improve the risk prediction and enhance equity in applying PRS in clinical practice.

A pleiotropic variant in DNAJB4 is associated with multiple myeloma risk

Pleiotropy, which consists of a single gene or allelic variant affecting multiple unrelated traits, is common across cancers, with evidence for genome-wide significant loci shared across cancer and noncancer traits. This feature is particularly relevant in multiple myeloma (MM) because several susceptibility loci that have been identified to date are pleiotropic. Therefore, the aim of this study was to identify novel pleiotropic variants involved in MM risk using 28 684 independent single nucleotide polymorphisms (SNPs) from GWAS Catalog that reached a significant association (P < 5 × 10^-8 ) with their respective trait. The selected SNPs were analyzed in 2434 MM cases and 3446 controls from the International Lymphoma Epidemiology Consortium (InterLymph). The 10 SNPs showing the strongest associations with MM risk in InterLymph were selected for replication in an independent set of 1955 MM cases and 1549 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and 418 MM cases and 147 282 controls from the FinnGen project. The combined analysis of the three studies identified an association between DNAJB4-rs34517439-A and an increased risk of developing MM (OR = 1.22, 95%CI 1.13-1.32, P = 4.81 × 10^-7 ). rs34517439-A is associated with a modified expression of the FUBP1 gene, which encodes a multifunctional DNA and RNA-binding protein that it was observed to influence the regulation of various genes involved in cell cycle regulation, among which various oncogenes and oncosuppressors. In conclusion, with a pleiotropic scan approach we identified DNAJB4-rs34517439 as a potentially novel MM risk locus.

Clonal hematopoiesis and risk of prostate cancer in large samples of European ancestry men

Little is known regarding the potential relationship between clonal hematopoiesis (CH) of indeterminate potential (CHIP), which is the expansion of hematopoietic stem cells with somatic mutations, and risk of prostate cancer, the fifth leading cause of cancer death of men worldwide. We evaluated the association of age-related CHIP with overall and aggressive prostate cancer risk in two large whole-exome sequencing studies of 75 047 European ancestry men, including 7663 prostate cancer cases, 2770 of which had aggressive disease, and 3266 men carrying CHIP variants. We found that CHIP, defined by over 50 CHIP genes individually and in aggregate, was not significantly associated with overall (aggregate HR = 0.93, 95% CI = 0.76-1.13, P = 0.46) or aggressive (aggregate OR = 1.14, 95% CI = 0.92-1.41, P = 0.22) prostate cancer risk. CHIP was weakly associated with genetic risk of overall prostate cancer, measured using a polygenic risk score (OR = 1.05 per unit increase, 95% CI = 1.01-1.10, P = 0.01). CHIP was not significantly associated with carrying pathogenic/likely pathogenic/deleterious variants in DNA repair genes, which have previously been found to be associated with aggressive prostate cancer. While findings from this study suggest that CHIP is likely not a risk factor for prostate cancer, it will be important to investigate other types of CH in association with prostate cancer risk.

Abstract PR008: Evaluation of intra-tumoral pathogenic bacteria pks+ E. coli, enterotoxigenic B. fragilis and fusobacterium nucleatum in 3695 colorectal cancer cases

Background: The microbiome in colorectal cancer (CRC) tumors is now recognized as a potent facilitator of tumorigenesis. Identification of the intra-tumoral microbiome may provide an avenue to enhance prediction of adverse outcomes and ultimately lead to improved survival of CRC, the second leading cause of cancer death. Methods: We performed targeted sequencing of tumor and matched normal DNA samples for 3695 colorectal cancer cases within the Genetics and Colorectal Cancer Consortium (GECCO) and the Colon Cancer Family Registry (CCFR). Tumor DNA was extracted primarily from formalin fixed paraffin embedded (FFPE) or from fresh frozen tissue. FFPE slides were dissected to ensure &gt;70% tumor content. Normal DNA was primarily extracted from blood and otherwise normal surrounding colorectal tissue. We designed probes for bacterial pathogens including Bacteroides fragilis (ETBF), polyketide synthase positive E. coli (pks+ EC), enterotoxigenic and Fusobacterium nucleatum (Fn). Logistic regression was used to assess tumor-associated bacteria and various clinical parameters, tumor molecular characteristics, and tumor stage adjusted for age, sex, and study and a Cox promotional hazard to assess association with CRC-specific survival adjusted for age at diagnosis, sex, mutational burden and microsatellite instability (MSI) and stratified baseline hazards by study. Survival data were available for 2357 patients. Results: The prevalence of patients with bacteria in tumors varied for ETBF (5.2%), pks+ EC (12.9%), and Fn (10.4%). Men had a lower probability of tumors with Fn (odds ratio (OR)=0.8, 95% confidence interval (CI)=0.64-0.99,p=0.04) but higher odds of tumors with pks+ EC (OR=1.38, CI=1.14-1.68,p=0.001). Patients with proximal tumors had higher odds of carrying Fn (OR=1.76, CI=1.38-2.26, p&lt;0.001) as did patients with MSI high tumors (OR 2.79, CI=2.08-3.74) p&lt;0.001). Patients with MSI-high tumors had lower odds of carrying pks+ EC (OR=0.63, CI=0.44-0.88, p=0.008) when compared to patients with microsatellite stable tumors. Among non-hypermutated tumors, we observed that tumors positive for pks+ EC were associated with better CRC-specific survival (OR=0.73, CI=0.57-0.93,p=0.01 whereas in patients with ETBF positive tumors had shorter survival (OR=1.99,CI=1,18-3.36,p=0.009) as did patients with tumors positive for Fn (OR=1.33, CI=1.04-1.72, p=0.03). Conclusions: In this largest investigation of intra-tumoral bacteria in CRC, we demonstrate that inclusion of pathogenic intra-tumoral microbiome may enhance prediction of survival and can identify a subset of patients with CRC that may benefit from targeted therapies that modulate the tumor microbiome.

Citation Format: Meredith A. J. Hullar, Keith R. Curtis, Tabitha A. Harrison, Yi Lin, Robert S. Steinfelder, Sonja I. Berndt, Daniel D. Buchanan, Andrew T. Chan, David Drew, Jane C. Figueiredo, Amy J. French, Ellen L. Goode, Mark A. Jenkins, Yohannes A. Melaku, Victor Moreno, Shuji Ogino, Stephen N. Thibideau, Amanda I. Phipps, Ulrike Peters. Evaluation of intra-tumoral pathogenic bacteria pks+ E. coli, enterotoxigenic B. fragilis and fusobacterium nucleatum in 3695 colorectal cancer cases [abstract]. In: Proceedings of the AACR Special Conference on Colorectal Cancer; 2022 Oct 1-4; Portland, OR. Philadelphia (PA): AACR; Cancer Res 2022;82(23 Suppl_1):Abstract nr PR008.

Distinct germline genetic susceptibility profiles identified for common non-Hodgkin lymphoma subtypes

Lymphoma risk is elevated for relatives with common non-Hodgkin lymphoma (NHL) subtypes, suggesting shared genetic susceptibility across subtypes. To evaluate the extent of mutual heritability among NHL subtypes and discover novel loci shared among subtypes, we analyzed data from eight genome-wide association studies within the InterLymph Consortium, including 10,629 cases and 9505 controls. We utilized Association analysis based on SubSETs (ASSET) to discover loci for subsets of NHL subtypes and evaluated shared heritability across the genome using Genome-wide Complex Trait Analysis (GCTA) and polygenic risk scores. We discovered 17 genome-wide significant loci (P < 5 × 10-8) for subsets of NHL subtypes, including a novel locus at 10q23.33 (HHEX) (P = 3.27 × 10-9). Most subset associations were driven primarily by only one subtype. Genome-wide genetic correlations between pairs of subtypes varied broadly from 0.20 to 0.86, suggesting substantial heterogeneity in the extent of shared heritability among subtypes. Polygenic risk score analyses of established loci for different lymphoid malignancies identified strong associations with some NHL subtypes (P < 5 × 10-8), but weak or null associations with others. Although our analyses suggest partially shared heritability and biological pathways, they reveal substantial heterogeneity among NHL subtypes with each having its own distinct germline genetic architecture.

A saturated map of common genetic variants associated with human height

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.