Citations

Abstract Pancreatic cancer disproportionately affects Black individuals in the United States, but they have limited representation in genetic studies of pancreatic ductal adenocarcinoma (PDAC). To address this gap, we performed admixture mapping and genome-wide association analysis (GWAS) in genetically inferred African ancestry individuals (1,030 cases and 889 controls). Admixture mapping identified three regions with a significantly higher proportion of African ancestry in cases compared to controls (5q33.3, 10p1, 22q12.3). GWAS identified a genome-wide significant association at 5p15.33 ( CLPTM1L , rs383009:T>C, T Allele Frequency=0.51, OR:1.45, P value=1.24×10 -8 ), a locus previously associated with PDAC. Known loci at 5p15.33, 7q32.3, 8q24.21 and 7q25.1 also replicated (P value <0.01). Multi-ancestral fine-mapping identified two potential causal SNPs (rs3830069 and rs2735940) at 5p15.33. Collectively these findings identified novel PDAC risk loci and expanded our understanding of this deadly cancer in underrepresented populations, emphasizing the multifactorial nature of PDAC risk including inherited genetic and non-genetic factors. Statement of Significance To understand how genetic variation contributes to PDAC risk in Black people in North American, we studied individuals of genetically-inferred African ancestry. We identified novel risk loci and differences in the contribution of known loci. This demonstrates that ancestry-informed genetic analyses improve our understanding of PDAC risk and enhances discovery.

Abstract Monoclonal gammopathy of undetermined significance (MGUS) is a necessary precursor condition to multiple myeloma (MM). Given the role of autophagy in modulating MM risk, we investigated whether genetic variation in autophagy‐related genes influences susceptibility to MGUS. We analyzed the association of 34,042 common autophagy‐related single nucleotide polymorphisms (SNPs) with MGUS across six independent cohorts, five from Europe and one from North America, comprising 2317 MGUS cases and 282,358 controls. We also assessed their impact on immune parameters, including absolute counts of 91 blood‐derived immune cell subsets and 103 circulating immunological proteins. Meta‐analysis revealed a genome‐wide significant association between the ULK4 rs6599175C allele and increased MGUS risk ( p = 3.35 × 10 −8 ). Carriers of this allele showed reduced counts of memory B cell subsets (IgM+CD38+CD27+ and IgD+IgM+CD27+; p = .0038 and p = .0056, respectively) and natural effector B cells (CD24+CD38+IgD+IgM+ cells; p = .0060). Although these associations were not statistically significant after multiple testing correction, they suggest a role of ULK4 in early B‐cell differentiation. Additionally, the CDKN2A rs2811710 variant showed a suggestive association with MGUS risk ( p = 2.17 × 10 −4 ), affecting transcription factor binding involved in B‐cell proliferation and differentiation, although it lacked association with immune markers. In conclusion, we confirm a genome‐wide significant association of the ULK4 locus and MGUS risk, supporting its role in early B‐cell differentiation, and identify CDKN2A as a candidate susceptibility locus warranting further investigation.

Abstract We investigated the influence of 55 583 autophagy-related single-nucleotide polymorphisms (SNPs) on chronic lymphocytic leukemia (CLL) risk across 4 independent populations comprising 5472 CLL cases and 726 465 controls. We also examined their impact on overall survival (OS), time to first treatment (TTFT), autophagy flux, and immune responses. A meta-analysis of the 4 populations identified, to our knowledge, for the first time, significant associations between CDKN2A (rs3731204) and BCL2 (rs4940571, rs12457371, and rs1026825) SNPs and CLL risk, with CDKN2A showing the strongest association (P = 1.57 × 10−12). We also validated previously reported associations for FAS, BCL2, and BAK1 SNPs with CLL risk (P = 4.73 × 10−21 to 3.39 × 10−9). The CDKN2Ars3731204 and FASrs1926194 SNPs associated with increased CDKN2A and ACTA2 messenger RNA expression levels in the whole blood and/or lymphocytes (P = 5.1 × 10−7, P = 1.58 × 10−21, and P = 7.8 × 10−41), although no significant effect on autophagy flux was observed. However, associations were found between CDKN2A, BCL2, and FAS SNPs and various T-cell subsets, cytokine production, and circulating concentrations of interferon gamma, tumor necrosis factor–related apoptosis-inducing ligand, CD40, chemokine ligand 20, and interleukin-2 receptor subunit β proteins (P ≤ .005). No significant association was detected between autophagy variants and OS or TTFT, suggesting that these variants drive disease initiation rather than progression. In conclusion, this study identified 4 novel associations for CLL and provided insights into the biological pathways that influence CLL development.

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 ULK4rs6599175 SNP correlated with circulating concentrations of vitamin D3 (p = 4.0 × 10−4), whereas the IKBKErs17433804 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 CD46rs1142469 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 CDKN2Ars2811710 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.

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.

AbstractPleiotropy, 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.

Abstract Background: Genome-wide association studies (GWAS) of multiple myeloma in populations of European ancestry (EA) identified and confirmed 24 susceptibility loci. For other cancers (e.g., colorectum and melanoma), risk loci have also been associated with patient survival. Methods: We explored the possible association of all the known risk variants and their polygenic risk score (PRS) with multiple myeloma overall survival (OS) in multiple populations of EA [the International Multiple Myeloma rESEarch (IMMEnSE) consortium, the International Lymphoma Epidemiology consortium, CoMMpass, and the German GWAS] for a total of 3,748 multiple myeloma cases. Cox proportional hazards regression was used to assess the association between each risk SNP with OS under the allelic and codominant models of inheritance. All analyses were adjusted for age, sex, country of origin (for IMMEnSE) or principal components (for the others) and disease stage (ISS). SNP associations were meta-analyzed. Results: SNP associations were meta-analyzed. From the meta-analysis, two multiple myeloma risk SNPs were associated with OS (P &lt; 0.05), specifically POT1-AS1-rs2170352 [HR = 1.37; 95% confidence interval (CI) = 1.09–1.73; P = 0.007] and TNFRSF13B-rs4273077 (HR = 1.19; 95% CI = 1.01–1.41; P = 0.04). The association between the combined 24 SNP MM-PRS and OS, however, was not significant. Conclusions: Overall, our results did not support an association between the majority of multiple myeloma risk SNPs and OS. Impact: This is the first study to investigate the association between multiple myeloma PRS and OS in multiple myeloma.

AbstractGene expression profiling can be used for predicting survival in multiple myeloma (MM) and identifying patients who will benefit from particular types of therapy. Some germline single nucleotide polymorphisms (SNPs) act as expression quantitative trait loci (eQTLs) showing strong associations with gene expression levels. We performed an association study to test whether eQTLs of genes reported to be associated with prognosis of MM patients are directly associated with measures of adverse outcome. Using the genotype‐tissue expression portal, we identified a total of 16 candidate genes with at least one eQTL SNP associated with their expression with P < 10−7 either in EBV‐transformed B‐lymphocytes or whole blood. We genotyped the resulting 22 SNPs in 1327 MM cases from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and examined their association with overall survival (OS) and progression‐free survival (PFS), adjusting for age, sex, country of origin and disease stage. Three polymorphisms in two genes (TBRG4‐rs1992292, TBRG4‐rs2287535 and ENTPD1‐rs2153913) showed associations with OS at P < .05, with the former two also associated with PFS. The associations of two polymorphisms in TBRG4 with OS were replicated in 1277 MM cases from the International Lymphoma Epidemiology (InterLymph) Consortium. A meta‐analysis of the data from IMMEnSE and InterLymph (2579 cases) showed that TBRG4‐rs1992292 is associated with OS (hazard ratio = 1.14, 95% confidence interval 1.04‐1.26, P = .007). In conclusion, we found biologically a plausible association between a SNP in TBRG4 and OS of MM patients.

Abstract Germline variation and smoking are independently associated with pancreatic ductal adenocarcinoma (PDAC). We conducted genome-wide smoking interaction analysis of PDAC using genotype data from four previous genome-wide association studies in individuals of European ancestry (7,937 cases and 11,774 controls). Examination of expression quantitative trait loci data from the Genotype-Tissue Expression Project followed by colocalization analysis was conducted to determine whether there was support for common SNP(s) underlying the observed associations. Statistical tests were two sided and P &lt; 5 × 10–8 was considered statistically significant. Genome-wide significant evidence of qualitative interaction was identified on chr2q21.3 in intron 5 of the transmembrane protein 163 (TMEM163) and upstream of the cyclin T2 (CCNT2). The most significant SNP using the Empirical Bayes method, in this region that included 45 significantly associated SNPs, was rs1818613 [per allele OR in never smokers 0.87, 95% confidence interval (CI), 0.82–0.93; former smokers 1.00, 95% CI, 0.91–1.07; current smokers 1.25, 95% CI 1.12–1.40, Pinteraction = 3.08 × 10–9). Examination of the Genotype-Tissue Expression Project data demonstrated an expression quantitative trait locus in this region for TMEM163 and CCNT2 in several tissue types. Colocalization analysis supported a shared SNP, rs842357, in high linkage disequilibrium with rs1818613 (r2 = 0. 94) driving both the observed interaction and the expression quantitative trait loci signals. Future studies are needed to confirm and understand the differential biologic mechanisms by smoking status that contribute to our PDAC findings. Significance: This large genome-wide interaction study identifies a susceptibility locus on 2q21.3 that significantly modified PDAC risk by smoking status, providing insight into smoking-associated PDAC, with implications for prevention.

Abstract Background: It is not known whether modifiable lifestyle factors that predict survival after invasive breast cancer differ by subtype. Methods: We analyzed data for 121,435 women diagnosed with breast cancer from 67 studies in the Breast Cancer Association Consortium with 16,890 deaths (8,554 breast cancer specific) over 10 years. Cox regression was used to estimate associations between risk factors and 10-year all-cause mortality and breast cancer–specific mortality overall, by estrogen receptor (ER) status, and by intrinsic-like subtype. Results: There was no evidence of heterogeneous associations between risk factors and mortality by subtype (Padj &gt; 0.30). The strongest associations were between all-cause mortality and BMI ≥30 versus 18.5–25 kg/m2 [HR (95% confidence interval (CI), 1.19 (1.06–1.34)]; current versus never smoking [1.37 (1.27–1.47)], high versus low physical activity [0.43 (0.21–0.86)], age ≥30 years versus &lt;20 years at first pregnancy [0.79 (0.72–0.86)]; &gt;0–&lt;5 years versus ≥10 years since last full-term birth [1.31 (1.11–1.55)]; ever versus never use of oral contraceptives [0.91 (0.87–0.96)]; ever versus never use of menopausal hormone therapy, including current estrogen–progestin therapy [0.61 (0.54–0.69)]. Similar associations with breast cancer mortality were weaker; for example, 1.11 (1.02–1.21) for current versus never smoking. Conclusions: We confirm associations between modifiable lifestyle factors and 10-year all-cause mortality. There was no strong evidence that associations differed by ER status or intrinsic-like subtype. Impact: Given the large dataset and lack of evidence that associations between modifiable risk factors and 10-year mortality differed by subtype, these associations could be cautiously used in prognostication models to inform patient-centered care.