McKay, James

Abstract Background: Clinical, molecular, and genetic epidemiology studies displayed remarkable differences between ever- and never-smoking lung cancer. Methods: We conducted a stratified multi-population (European, East Asian, and African descent) association study on 44,823 ever-smokers and 20,074 never-smokers to identify novel variants that were missed in the non-stratified analysis. Functional analysis including expression quantitative trait loci (eQTL) colocalization and DNA damage assays, and annotation studies were conducted to evaluate the functional roles of the variants. We further evaluated the impact of smoking quantity on lung cancer risk for the variants associated with ever-smoking lung cancer. Results: Five novel independent loci, GABRA4, intergenic region 12q24.33, LRRC4C, LINC01088, and LCNL1 were identified with the association at two or three populations (P < 5 × 10−8). Further functional analysis provided multiple lines of evidence suggesting the variants affect lung cancer risk through excessive DNA damage (GABRA4) or cis-regulation of gene expression (LCNL1). The risk of variants from 12 independent regions, including the well-known CHRNA5, associated with ever-smoking lung cancer was evaluated for never-smokers, light-smokers (packyear ≤ 20), and moderate-to-heavy-smokers (packyear > 20). Different risk patterns were observed for the variants among the different groups by smoking behavior. Conclusions: We identified novel variants associated with lung cancer in only ever- or never-smoking groups that were missed by prior main-effect association studies. Impact: Our study highlights the genetic heterogeneity between ever- and never-smoking lung cancer and provides etiologic insights into the complicated genetic architecture of this deadly cancer.

Abstract Background People with cancer experience high rates of venous thromboembolism (VTE). Risk of subsequent cancer is also increased in people experiencing their first VTE. The causal mechanisms underlying this association are not completely understood, and it is unknown whether VTE is itself a risk factor for cancer. Methods We used data from large genome-wide association study meta-analyses to perform bidirectional Mendelian randomization analyses to estimate causal associations between genetic liability to VTE and risk of 18 different cancers. Results We found no conclusive evidence that genetic liability to VTE was causally associated with an increased incidence of cancer, or vice versa. We observed an association between liability to VTE and pancreatic cancer risk [odds ratio for pancreatic cancer: 1.23 (95% confidence interval: 1.08–1.40) per log-odds increase in VTE risk, P = 0.002]. However, sensitivity analyses revealed this association was predominantly driven by a variant proxying non-O blood group, with inadequate evidence to suggest a causal relationship. Conclusions These findings do not support the hypothesis that genetic liability to VTE is a cause of cancer. Existing observational epidemiological associations between VTE and cancer are therefore more likely to be driven by pathophysiological changes which occur in the setting of active cancer and anti-cancer treatments. Further work is required to explore and synthesize evidence for these mechanisms.

Abstract Background Epidemiological studies that examined the association between Parkinson's disease (PD) and cancers led to inconsistent results, but they face a number of methodological difficulties. Objective We used results from genome‐wide association studies (GWASs) to study the genetic correlation between PD and different cancers to identify common genetic risk factors. Methods We used individual data for participants of European ancestry from the Courage‐PD (Comprehensive Unbiased Risk Factor Assessment for Genetics and Environment in Parkinson's Disease; PD, N = 16,519) and EPITHYR (differentiated thyroid cancer, N = 3527) consortia and summary statistics of GWASs from iPDGC (International Parkinson Disease Genomics Consortium; PD, N = 482,730), Melanoma Meta‐Analysis Consortium (MMAC), Breast Cancer Association Consortium (breast cancer), the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (prostate cancer), International Lung Cancer Consortium (lung cancer), and Ovarian Cancer Association Consortium (ovarian cancer) (N comprised between 36,017 and 228,951 for cancer GWASs). We estimated the genetic correlation between PD and cancers using linkage disequilibrium score regression. We studied the association between PD and polymorphisms associated with cancers, and vice versa, using cross‐phenotypes polygenic risk score (PRS) analyses. Results We confirmed a previously reported positive genetic correlation of PD with melanoma (G corr = 0.16 [0.04; 0.28]) and reported an additional significant positive correlation of PD with prostate cancer (G corr = 0.11 [0.03; 0.19]). There was a significant inverse association between the PRS for ovarian cancer and PD (odds ratio [OR] = 0.89 [0.84; 0.94]). Conversely, the PRS of PD was positively associated with breast cancer (OR = 1.08 [1.06; 1.10]) and inversely associated with ovarian cancer (OR = 0.95 [0.91; 0.99]). The association between PD and ovarian cancer was mostly driven by rs183211 located in an intron of the NSF gene (17q21.31). Conclusions We show evidence in favor of a contribution of pleiotropic genes to the association between PD and specific cancers. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Abstract Differences by sex in lung cancer incidence and mortality have been reported which cannot be fully explained by sex differences in smoking behavior, implying existence of genetic and molecular basis for sex disparity in lung cancer development. However, the information about sex dimorphism in lung cancer risk is quite limited despite the great success in lung cancer association studies. By adopting a stringent two-stage analysis strategy, we performed a genome-wide gene–sex interaction analysis using genotypes from a lung cancer cohort including ~ 47 000 individuals with European ancestry. Three low-frequency variants (minor allele frequency < 0.05), rs17662871 [odds ratio (OR) = 0.71, P = 4.29×10−8); rs79942605 (OR = 2.17, P = 2.81×10−8) and rs208908 (OR = 0.70, P = 4.54×10−8) were identified with different risk effect of lung cancer between men and women. Further expression quantitative trait loci and functional annotation analysis suggested rs208908 affects lung cancer risk through differential regulation of Coxsackie virus and adenovirus receptor gene expression in lung tissues between men and women. Our study is one of the first studies to provide novel insights about the genetic and molecular basis for sex disparity in lung cancer development.

Abstract There is growing evidence showing that sex disparity in lung cancer risk cannot be fully explained by sex differences in smoking behavior, implying existence of genetic and molecular basis for sex disparity in lung cancer development. However, the information about sex dimorphism in lung cancer risk is quite limited despite the great success in lung cancer association studies. By adopting a stringent two-stage analysis strategy, we performed a genome-wide gene-sex interaction analysis using genotype from a lung cancer cohort including ~ 47,000 individuals with European ancestry. Three low-frequency variants (minor allele frequency < 0.05), rs17662871 (OR=0.71, p=4.29x10-8); rs79942605 (OR=2.17, p=2.81x10-8); and rs208908 (OR=0.70, p=4.54x10-8), were identified with different risk effect of lung cancer between men and women. Further eQTL and functional annotation analysis suggested rs208908 affect lung cancer risk through differential regulation of CXADR (Coxsackie Virus And Adenovirus Receptor) gene expression in lung tissues between men and women. Our study is one of the first studies to provide novel insights about the genetic and molecular basis for sex disparity in lung cancer development.

Abstract Inherited genetic risk factors play a role in multiple myeloma (MM), yet considerable missing heritability exists. Rare risk variants at genome-wide association study (GWAS) loci are a new avenue to explore. Pleiotropy between lymphoid neoplasms (LNs) has been suggested in family history and genetic studies, but no studies have interrogated sequencing for pleiotropic genes or rare risk variants. Sequencing genetically enriched cases can help discover rarer variants. We analyzed exome sequencing in familial or early-onset MM cases to identify rare, functionally relevant variants near GWAS loci for a range of LNs. A total of 149 distinct and significant LN GWAS loci have been published. We identified six recurrent, rare, potentially deleterious variants within 5 kb of significant GWAS single nucleotide polymorphisms in 75 MM cases. Mutations were observed in BTNL2, EOMES, TNFRSF13B, IRF8, ACOXL and TSPAN32. All six genes replicated in an independent set of 255 early-onset MM or familial MM or precursor cases. Expansion of our analyses to the full length of these six genes resulted in a list of 39 rare and deleterious variants, seven of which segregated in MM families. Three genes also had significant rare variant burden in 733 sporadic MM cases compared with 935 control individuals: IRF8 (P = 1.0 × 10−6), EOMES (P = 6.0 × 10−6) and BTNL2 (P = 2.1 × 10−3). Together, our results implicate six genes in MM risk, provide support for genetic pleiotropy between LN subtypes and demonstrate the utility of sequencing genetically enriched cases to identify functionally relevant variants near GWAS loci.

Abstract Background: Lipid traits have been inconsistently linked to risk of non-Hodgkin lymphoma (NHL). We examined the association of genetically predicted lipid traits with risk of diffuse large B-cell lymphoma (DLBCL), chronic lymphocytic leukemia (CLL), follicular lymphoma (FL), and marginal zone lymphoma (MZL) using Mendelian randomization (MR) analysis. Methods: Genome-wide association study data from the InterLymph Consortium were available for 2,661 DLBCLs, 2,179 CLLs, 2,142 FLs, 824 MZLs, and 6,221 controls. SNPs associated (P &lt; 5 × 10−8) with high-density lipoprotein (HDL, n = 164), low-density lipoprotein (LDL, n = 137), total cholesterol (TC, n = 161), and triglycerides (TG, n = 123) were used as instrumental variables (IV), explaining 14.6%, 27.7%, 16.8%, and 12.8% of phenotypic variation, respectively. Associations between each lipid trait and NHL subtype were calculated using the MR inverse variance–weighted method, estimating odds ratios (OR) per standard deviation and 95% confidence intervals (CI). Results: HDL was positively associated with DLBCL (OR = 1.14; 95% CI, 1.00–1.30) and MZL (OR = 1.09; 95% CI, 1.01–1.18), while TG was inversely associated with MZL risk (OR = 0.90; 95% CI, 0.83–0.99), all at nominal significance (P &lt; 0.05). A positive trend was observed for HDL with FL risk (OR = 1.08; 95% CI, 0.99–1.19; P = 0.087). No associations were noteworthy after adjusting for multiple testing. Conclusions: We did not find evidence of a clear or strong association of these lipid traits with the most common NHL subtypes. While these IVs have been previously linked to other cancers, our findings do not support any causal associations with these NHL subtypes. Impact: Our results suggest that prior reported inverse associations of lipid traits are not likely to be causal and could represent reverse causality or confounding.

Abstract Background Recent genome wide association studies (GWAS) of Hodgkin lymphoma (HL) have identified several associations at both HLA and non-HLA loci. However, much of HL heritability remains unexplained. Methods To identify novel risk loci, we performed a meta-analysis of 3 HL GWAS including a total of 1,810 cases and 7,879 controls. Results were replicated in an independent set of 1,163 cases and 2,580 controls, for a total of 3,097 and 11,097 cases and controls combined, respectively. participants in discovery and replication stages were of European descent. quality control and imputation we conducted a meta-analysis addressing 1,004,829 variants (λ= 1.10, λ1000= 1.03). Associations between SNP genotypes and HL risk were evaluated under a log-additive model of inheritance adjusting for sex, study center and significant principal components to control for population stratification. We performed an analysis with all HL cases and then conducted stratified analyses by histological subtype (classical, nodular sclerosis and mixed cellularity), age at diagnosis (nodular sclerosis among those diagnosed at 15- 35 years in all studies, and those diagnosed at 35 and older years in the European Study only) and EBV tumor status (negative and positive). We then used a bioinformatic approach (FunciSNP) to identify potential functional variants associated with HD risk correlated with risk loci of interest. We extracted publically available ENCODE data on biofeatures to identify potential functional motifs associated with the index SNP or correlated SNPs. Finally, we measured expression levels of the two alternative mRNA transcripts in lymphoblastoid cell lines (LCLs) from 49 post-therapy HL patients and 25 unaffected controls. RT-PCR was carried out in triplicate. Relative expression levels were calculated relative to TBP as housekeeping gene. Linear models were used to assess correlation between genotype and TCF3expression levels. Results The meta-analysis identified a novel susceptibility variant at chromosome 19p13.3 (rs1860661) associated with HL risk (Odds Ratio [OR]= 0.78, P=2.0*10-8, I2=0%). variant is located in intron 2 of TCF3 (also known as E2A), a regulator of B- and T-cell lineage commitment. was also significantly associated with HL (OR= 0.85, P=0.002) in the replication series of 1,281 cases and 3,218 controls. the combined analysis consisting of the discovery and replication sets, rs1860661was strongly associated with HL (OR=0.81,=3.5*10-10), with no evidence of heterogeneity between contributing studies (Phom=0.41, I2=0%). The number of G alleles defined by rs1860661 was significantly associated with a reduced risk of each HL subtype except EBV positive HL. rs1860661 and two correlated SNPs, rs10413888 (r2=0.90) and rs8103453 (r2=0.89) identified by FunciSNP analysis map in or near marks of open chromatin and in DNAse hypersensitivity sites in TCF3 in CD20+ B cell lines., the protective minor alleles of these SNPs as defined by the G-G-G haplotype map to the binding sites for ZBTB7a (rs10413888 and rs1860661) and (rs8103453) transcription factors, likely improving the binding efficiency to the sites which may result in increased transcription rates of TCF3. TCF3 is encoded by two alternative transcripts (E12 and E47). Higher expression levels of TCF3-E47, whose transcription start site is located close to rs1860661, was associated with the rs1860661-G allele in controls (P=0.02), but not in HL patients (P=0.22). Conclusion/Discussion TCF3 is essential for the commitment of lymphoid progenitors to both B-cell and T-cell lineage development. A molecular and phenotypic hallmark of classical HL is the loss of the B-cell phenotype in HRS cells, including lack of demonstrable B-cell receptor and most B-cell specific markers such as CD19 or CD20. HRS cells have a low level of TCF3, particularly homodimers of the isoform E47, due to expression of the ABF-1 and ID2 inhibitors that bind to TCF3. Thus, higher TCF3 levels in HRS precursor cells may lead to enhanced retention of the B cell phenotype, thereby conferring a protective effect. These data suggest a link between the 19p13.3 locus including TCF3 and HL risk, indicating that TCF3 could be relevant to HL etiology and pathogenesis. Disclosures: Link: Genentech: Consultancy; Millenium: Consultancy; Pharmacyclics: Consultancy; Spectrum: Consultancy.

Abstract Introduction Multiple myeloma (MM) is a result of a malignant transformation of plasma cells that is preceded by the presence of an asymptomatic clonal plasma cell expansion, a condition referred to as monoclonal gammopathy of undetermined significance (MGUS). We and others have shown familial aggregation of MM and MGUS. Evidence from epidemiologic, family and genome-wide association studies (GWAS) suggests a genetic component underlying MM etiology. GWAS have successfully established 17 common genetic risk loci for MM to date and recently, rare inherited susceptibility variants in the LSD1 / KDM1A and USP45 genes were identified in familial MM / MGUS kindreds. Family-based approaches may be used to elucidate genetic variation contributing to familial MM. Genetic linkage analysis has historically been used to detect the chromosomal location of disease genes. The objective of this study was to conduct a linkage analysis of MM / MGUS families to identify genomic regions for MM / MGUS. Methods Linkage analysis was performed on whole exome sequencing (WES) data generated from germline DNA extracted from peripheral blood from MM / MGUS families ascertained at four sites, of which 79 were selected with 2 or more affected members with any combination of MM or MGUS. Whole exome capture was performed using Agilent SureSelect 38 Mb paired end sequencing and ran on Illumina HiSeq 2000s/2500s; standard analyses for alignment (GRCh37) and quality control were conducted. The Genome Analysis Toolkit's (GATK) HaplotypeCaller in per-sample mode was used to jointly call all the samples together. Quality control included removing variants that had <75% call rate, <8x coverage, or minor allele frequency <0.01 (1KGenomes). The jointly called WES data from the 79 families consisted of 230 individuals, including 119 MM, 93 MGUS and 18 unaffected relatives. Additional early onset (defined as age <60) MM cases (n=400), sporadic MM cases (n=879) and controls (n=2389) were also included in the joint calling to be used for follow-up of linkage regions. We included independent variants defined as linkage disequilibrium (LD) r2 <0.05 using both PLINK 1.07's LD based variant pruner and also MERLIN's pairwise r2 cluster marker-marker modeling to ensure variant independence. After quality control and LD filtering: 12,946 variants remained. We conducted multipoint linkage analysis using MERLIN (Multipoint Engine for Rapid Likelihood Inference using the Lander-Green approach) and considered both non-parametric and parametric (dominant and recessive) models to test for co-segregation of chromosomal regions with MM / MGUS. A logarithm of the odds (LOD) score greater than 3.3 was significant evidence for linkage. Within the same families, we interrogated any suggestive regions by performing gene-level association tests using the PEDGENE R package to calculate burden and kernel statistics, accounting for family relationships. Follow-up of these genes in the early onset and sporadic MM cases and controls and interrogation of single variants driving the associations of these genes is ongoing. Results Among the 79 families, 28 consisted of two or more MM, 41 had two or more MM and MGUS cases and 10 had two or more MGUS. Analyses of the 79 MM / MGUS families identified chromosome 6 at 123420001-149070000 BP region, with evidence for linkage by the non-parametric model (LOD score=3.3) and supportive evidence for linkage by the dominant parametric model (LOD=2.5) (Figure). The chromosome 6 region, which is outside the HLA-region, contains a total of 72 genes. Using gene-level association tests, we identified four genes within the region (MOXD1, TRDN, ADAT2, LTV1) that were associated within MM / MGUS in the family cases (p<0.05). Discussion We found evidence for a locus on chromosome 6 linked to MM / MGUS. Four genes in this region may be associated with MM / MGUS in these families. Follow-up of these genes/variants in high-risk individuals (early-onset / age <60), sporadic cases and controls will help elucidate genetic mechanisms contributing to familial MM / MGUS. In addition, we demonstrate the benefit of applying a linkage analysis framework to familial WES data for discovery of genomic regions influencing MM / MGUS. Disclosures Dumontet: Roche: Research Funding; Merck: Consultancy, Membership on an entity's Board of Directors or advisory committees; Janssen: Honoraria; Sanofi: Honoraria. Kumar:AbbVie: Membership on an entity's Board of Directors or advisory committees, Research Funding; Janssen: Membership on an entity's Board of Directors or advisory committees, Research Funding; Celgene: Membership on an entity's Board of Directors or advisory committees, Research Funding; AbbVie: Membership on an entity's Board of Directors or advisory committees, Research Funding; Oncopeptides: Membership on an entity's Board of Directors or advisory committees; Takeda: Membership on an entity's Board of Directors or advisory committees; Merck: Membership on an entity's Board of Directors or advisory committees, Research Funding; Celgene: Membership on an entity's Board of Directors or advisory committees, Research Funding; KITE: Membership on an entity's Board of Directors or advisory committees, Research Funding; Novartis: Research Funding; Roche: Research Funding; Janssen: Membership on an entity's Board of Directors or advisory committees, Research Funding; KITE: Membership on an entity's Board of Directors or advisory committees, Research Funding.

Abstract We previously identified five single nucleotide polymorphisms (SNPs) at four susceptibility loci for diffuse large B-cell lymphoma (DLBCL) in individuals of European ancestry through a large genome-wide association study (GWAS). To further elucidate genetic susceptibility to DLBCL, we sought to validate two loci at 3q13.33 and 3p24.1 that were suggestive in the original GWAS with additional genotyping. In the meta-analysis (5662 cases and 9237 controls) of the four original GWAS discovery scans and three replication studies, the 3q13.33 locus (rs9831894; minor allele frequency [MAF] = 0.40) was associated with DLBCL risk [odds ratio (OR) = 0.83, P = 3.62 × 10−13]. rs9831894 is in linkage disequilibrium (LD) with additional variants that are part of a super-enhancer that physically interacts with promoters of CD86 and ILDR1. In the meta-analysis (5510 cases and 12 817 controls) of the four GWAS discovery scans and four replication studies, the 3p24.1 locus (rs6773363; MAF = 0.45) was also associated with DLBCL risk (OR = 1.20, P = 2.31 × 10−12). This SNP is 29 426-bp upstream of the nearest gene EOMES and in LD with additional SNPs that are part of a highly lineage-specific and tumor-acquired super-enhancer that shows long-range interaction with AZI2 promoter. These loci provide additional evidence for the role of immune function in the etiology of DLBCL, the most common lymphoma subtype.