Lung Cancer in Ever- and Never-Smokers: Findings from Multi-Population GWAS Studies

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.

Impact of individual level uncertainty of lung cancer polygenic risk score (PRS) on risk stratification

BACKGROUND

Although polygenic risk score (PRS) has emerged as a promising tool for predicting cancer risk from genome-wide association studies (GWAS), the individual-level accuracy of lung cancer PRS and the extent to which its impact on subsequent clinical applications remains largely unexplored.

METHODS

Lung cancer PRSs and confidence/credible interval (CI) were constructed using two statistical approaches for each individual: (1) the weighted sum of 16 GWAS-derived significant SNP loci and the CI through the bootstrapping method (PRS-16-CV) and (2) LDpred2 and the CI through posteriors sampling (PRS-Bayes), among 17,166 lung cancer cases and 12,894 controls with European ancestry from the International Lung Cancer Consortium. Individuals were classified into different genetic risk subgroups based on the relationship between their own PRS mean/PRS CI and the population level threshold.

RESULTS

Considerable variances in PRS point estimates at the individual level were observed for both methods, with an average standard deviation (s.d.) of 0.12 for PRS-16-CV and a much larger s.d. of 0.88 for PRS-Bayes. Using PRS-16-CV, only 25.0% of individuals with PRS point estimates in the lowest decile of PRS and 16.8% in the highest decile have their entire 95% CI fully contained in the lowest and highest decile, respectively, while PRS-Bayes was unable to find any eligible individuals. Only 19% of the individuals were concordantly identified as having high genetic risk (> 90th percentile) using the two PRS estimators. An increased relative risk of lung cancer comparing the highest PRS percentile to the lowest was observed when taking the CI into account (OR = 2.73, 95% CI: 2.12-3.50, P-value = 4.13 × 10-15) compared to using PRS-16-CV mean (OR = 2.23, 95% CI: 1.99-2.49, P-value = 5.70 × 10-46). Improved risk prediction performance with higher AUC was consistently observed in individuals identified by PRS-16-CV CI, and the best performance was achieved by incorporating age, gender, and detailed smoking pack-years (AUC: 0.73, 95% CI = 0.72-0.74).

CONCLUSIONS

Lung cancer PRS estimates using different methods have modest correlations at the individual level, highlighting the importance of considering individual-level uncertainty when evaluating the practical utility of PRS.

Diagnostic accuracy of HPV16 early antigen serology for HPV-driven oropharyngeal cancer is independent of age and sex

A growing proportion of head and neck cancer (HNC), especially oropharyngeal cancer (OPC), is caused by human papillomavirus (HPV). There are several markers for HPV-driven HNC, one being HPV early antigen serology. We aimed to investigate the diagnostic accuracy of HPV serology and its performance across patient characteristics. Data from the VOYAGER consortium was used, which comprises five studies on HNC from North America and Europe. Diagnostic accuracy, that is, sensitivity, specificity, Cohen's kappa and correctly classified proportions of HPV16 E6 serology, was assessed for OPC and other HNC using p16INK4a immunohistochemistry (p16), HPV in situ hybridization (ISH) and HPV PCR as reference methods. Stratified analyses were performed for variables including age, sex, smoking and alcohol use, to test the robustness of diagnostic accuracy. A risk-factor analysis based on serology was conducted, comparing HPV-driven to non-HPV-driven OPC. Overall, HPV serology had a sensitivity of 86.8% (95% CI 85.1-88.3) and specificity of 91.2% (95% CI 88.6-93.4) for HPV-driven OPC using p16 as a reference method. In stratified analyses, diagnostic accuracy remained consistent across sex and different age groups. Sensitivity was lower for heavy smokers (77.7%), OPC without lymph node involvement (74.4%) and the ARCAGE study (66.7%), while specificity decreased for cases with <10 pack-years (72.1%). The risk-factor model included study, year of diagnosis, age, sex, BMI, alcohol use, pack-years, TNM-T and TNM-N stage. HPV serology is a robust biomarker for HPV-driven OPC, and its diagnostic accuracy is independent of age and sex. Future research is suggested on the influence of smoking on HPV antibody levels.

A risk prediction model for head and neck cancers incorporating lifestyle factors, HPV serology and genetic markers

Head and neck cancer is often diagnosed late and prognosis for most head and neck cancer patients remains poor. To aid early detection, we developed a risk prediction model based on demographic and lifestyle risk factors, human papillomavirus (HPV) serological markers and genetic markers. A total of 10 126 head and neck cancer cases and 5254 controls from five North American and European studies were included. HPV serostatus was determined by antibodies for HPV16 early oncoproteins (E6, E7) and regulatory early proteins (E1, E2, E4). The data were split into a training set (70%) for model development and a hold-out testing set (30%) for model performance evaluation, including discriminative ability and calibration. The risk models including demographic, lifestyle risk factors and polygenic risk score showed a reasonable predictive accuracy for head and neck cancer overall. A risk model that also included HPV serology showed substantially improved predictive accuracy for oropharyngeal cancer (AUC = 0.94, 95% CI = 0.92-0.95 in men and AUC = 0.92, 95% CI = 0.88-0.95 in women). The 5-year absolute risk estimates showed distinct trajectories by risk factor profiles. Based on the UK Biobank cohort, the risks of developing oropharyngeal cancer among 60 years old and HPV16 seropositive in the next 5 years ranged from 5.8% to 14.9% with an average of 8.1% for men, 1.3% to 4.4% with an average of 2.2% for women. Absolute risk was generally higher among individuals with heavy smoking, heavy drinking, HPV seropositivity and those with higher polygenic risk score. These risk models may be helpful for identifying people at high risk of developing head and neck cancer.

Investigation of Shared Genetic Risk Factors Between Parkinson's Disease and Cancers

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 (Gcorr = 0.16 [0.04; 0.28]) and reported an additional significant positive correlation of PD with prostate cancer (Gcorr = 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.

Association of germline TYK2 variation with lung cancer and non-Hodgkin lymphoma risk

Deucravacitinib, a novel, selective inhibitor of TYK2 is currently under review at the FDA and EMA for treatment of moderate-to-severe plaque psoriasis. It is unclear whether recent safety concerns (ie, elevated rates of lung cancer and lymphoma) related to similar medications (ie, other JAK inhibitors) are shared with this novel TYK2 inhibitor. We used a partial loss-of-function variant in TYK2 (rs34536443), previously shown to protect against psoriasis and other autoimmune diseases, to evaluate the potential effect of therapeutic TYK2 inhibition on risk of lung cancer and non-Hodgkin lymphoma. Summary genetic association data on lung cancer risk were obtained from a GWAS meta-analysis of 29 266 cases and 56 450 controls in the Integrative Analysis of Lung Cancer Risk and Aetiology (INTEGRAL) consortium. Summary genetic association data on non-Hodgkin lymphoma risk were obtained from a GWAS meta-analysis of 8489 cases and 374 506 controls in the UK Biobank and InterLymph consortium. In the primary analysis, each copy of the minor allele of rs34536443, representing partial TYK2 inhibition, was associated with an increased risk of lung cancer (OR 1.15, 95% CI 1.09-1.23, P = 2.29 × 10-6 ) and non-Hodgkin lymphoma (OR 1.18, 95% CI 1.05-1.33, P = 5.25 × 10-3 ). Our analyses using an established partial loss-of-function mutation to mimic TYK2 inhibition provide genetic evidence that therapeutic TYK2 inhibition may increase risk of lung cancer and non-Hodgkin lymphoma. These findings, consistent with recent reports from postmarketing trials of similar JAK inhibitors, could have important implications for future safety assessment of deucravacitinib and other TYK2 inhibitors in development.

Circulating Isovalerylcarnitine and Lung Cancer Risk: Evidence from Mendelian Randomization and Prediagnostic Blood Measurements

BACKGROUND

Tobacco exposure causes 8 of 10 lung cancers, and identifying additional risk factors is challenging due to confounding introduced by smoking in traditional observational studies.

MATERIALS AND METHODS

We used Mendelian randomization (MR) to screen 207 metabolites for their role in lung cancer predisposition using independent genome-wide association studies (GWAS) of blood metabolite levels (n = 7,824) and lung cancer risk (n = 29,266 cases/56,450 controls). A nested case-control study (656 cases and 1,296 matched controls) was subsequently performed using prediagnostic blood samples to validate MR association with lung cancer incidence data from population-based cohorts (EPIC and NSHDS).

RESULTS

An MR-based scan of 207 circulating metabolites for lung cancer risk identified that blood isovalerylcarnitine (IVC) was associated with a decreased odds of lung cancer after accounting for multiple testing (log10-OR = 0.43; 95% CI, 0.29-0.63). Molar measurement of IVC in prediagnostic blood found similar results (log10-OR = 0.39; 95% CI, 0.21-0.72). Results were consistent across lung cancer subtypes.

CONCLUSIONS

Independent lines of evidence support an inverse association of elevated circulating IVC with lung cancer risk through a novel methodologic approach that integrates genetic and traditional epidemiology to efficiently identify novel cancer biomarkers.

IMPACT

Our results find compelling evidence in favor of a protective role for a circulating metabolite, IVC, in lung cancer etiology. From the treatment of a Mendelian disease, isovaleric acidemia, we know that circulating IVC is modifiable through a restricted protein diet or glycine and L-carnatine supplementation. IVC may represent a modifiable and inversely associated biomarker for lung cancer.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

A Large-Scale Genome-Wide Gene-Gene Interaction Study of Lung Cancer Susceptibility in Europeans With a Trans-Ethnic Validation in Asians

INTRODUCTION

Although genome-wide association studies have been conducted to investigate genetic variation of lung tumorigenesis, little is known about gene-gene (G × G) interactions that may influence the risk of non-small cell lung cancer (NSCLC).

METHODS

Leveraging a total of 445,221 European-descent participants from the International Lung Cancer Consortium OncoArray project, Transdisciplinary Research in Cancer of the Lung and UK Biobank, we performed a large-scale genome-wide G × G interaction study on European NSCLC risk by a series of analyses. First, we used BiForce to evaluate and rank more than 58 billion G × G interactions from 340,958 single-nucleotide polymorphisms (SNPs). Then, the top interactions were further tested by demographically adjusted logistic regression models. Finally, we used the selected interactions to build lung cancer screening models of NSCLC, separately, for never and ever smokers.

RESULTS

With the Bonferroni correction, we identified eight statistically significant pairs of SNPs, which predominantly appeared in the 6p21.32 and 5p15.33 regions (e.g., rs521828C6orf10 and rs204999PRRT1, ORinteraction = 1.17, p = 6.57 × 10-13; rs3135369BTNL2 and rs2858859HLA-DQA1, ORinteraction = 1.17, p = 2.43 × 10-13; rs2858859HLA-DQA1 and rs9275572HLA-DQA2, ORinteraction = 1.15, p = 2.84 × 10-13; rs2853668TERT and rs62329694CLPTM1L, ORinteraction = 0.73, p = 2.70 × 10-13). Notably, even with much genetic heterogeneity across ethnicities, three pairs of SNPs in the 6p21.32 region identified from the European-ancestry population remained significant among an Asian population from the Nanjing Medical University Global Screening Array project (rs521828C6orf10 and rs204999PRRT1, ORinteraction = 1.13, p = 0.008; rs3135369BTNL2 and rs2858859HLA-DQA1, ORinteraction = 1.11, p = 5.23 × 10-4; rs3135369BTNL2 and rs9271300HLA-DQA1, ORinteraction = 0.89, p = 0.006). The interaction-empowered polygenetic risk score that integrated classical polygenetic risk score and G × G information score was remarkable in lung cancer risk stratification.

CONCLUSIONS

Important G × G interactions were identified and enriched in the 5p15.33 and 6p21.32 regions, which may enhance lung cancer screening models.

High-risk subtypes of chronic lymphocytic leukemia are detectable as early as 16 years prior to diagnosis

Chronic lymphocytic leukemia (CLL) is preceded by monoclonal B-cell lymphocytosis (MBL), a CLL precursor state with a prevalence of up to 12% in aged individuals; however, the duration of MBL and the mechanisms of its evolution to CLL remain largely unknown. In this study, we sequenced the B-cell receptor (BcR) immunoglobulin heavy chain (IGH) gene repertoire of 124 patients with CLL and 118 matched controls in blood samples taken up to 22 years prior to diagnosis. Significant skewing in the BcR IGH gene repertoire was detected in the majority of patients, even before the occurrence of lymphocytosis and irrespective of the clonotypic IGH variable gene somatic hypermutation status. Furthermore, we identified dominant clonotypes belonging to major stereotyped subsets associated with poor prognosis up to 16 years before diagnosis in 14 patients with CLL. In 22 patients with longitudinal samples, the skewing of the BcR IGH gene repertoire increased significantly over time to diagnosis or remained stable at high levels. For 14 of 16 patients with available samples at diagnosis, the CLL clonotype was already present in the prediagnostic samples. Overall, our data indicate that the preclinical phase of CLL could be longer than previously thought, even in adverse-prognostic cases.

Studying the Utility of Using Genetics to Predict Smoking-Related Outcomes in a Population-Based Study and a Selected Cohort

INTRODUCTION

The purpose of this study is to examine the predictive utility of polygenic risk scores (PRSs) for smoking behaviors.

AIMS AND METHODS

Using summary statistics from the Sequencing Consortium of Alcohol and Nicotine use consortium, we generated PRSs of ever smoking, age of smoking initiation, cigarettes smoked per day, and smoking cessation for participants in the population-based Atherosclerosis Risk in Communities (ARIC) study (N = 8638), and the Collaborative Genetic Study of Nicotine Dependence (COGEND) (N = 1935). The outcomes were ever smoking, age of smoking initiation, heaviness of smoking, and smoking cessation.

RESULTS

In the European ancestry cohorts, each PRS was significantly associated with the corresponding smoking behavior outcome. In the ARIC cohort, the PRS z-score for ever smoking predicted smoking (odds ratio [OR]: 1.37; 95% confidence interval [CI]: 1.31, 1.43); the PRS z-score for age of smoking initiation was associated with age of smoking initiation (OR: 0.87; 95% CI: 0.82, 0.92); the PRS z-score for cigarettes per day was associated with heavier smoking (OR: 1.17; 95% CI: 1.11, 1.25); and the PRS z-score for smoking cessation predicted successful cessation (OR: 1.24; 95% CI: 1.17, 1.32). In the African ancestry cohort, the PRSs did not predict smoking behaviors.

CONCLUSIONS

Smoking-related PRSs were associated with smoking-related behaviors in European ancestry populations. This improvement in prediction is greatest in the lowest and highest genetic risk categories. The lack of prediction in African ancestry populations highlights the urgent need to increase diversity in research so that scientific advances can be applied to populations other than those of European ancestry.

IMPLICATIONS

This study shows that including both genetic ancestry and PRSs in a single model increases the ability to predict smoking behaviors compared with the model including only demographic characteristics. This finding is observed for every smoking-related outcome. Even though adding genetics is more predictive, the demographics alone confer substantial and meaningful predictive power. However, with increasing work in PRSs, the predictive ability will continue to improve.

Germline determinants of humoral immune response to HPV-16 protect against oropharyngeal cancer

Although several oropharyngeal cancer (OPC) susceptibility loci have been identified, most previous studies lacked detailed information on human papillomavirus (HPV) status. We conduct a genome-wide analysis by HPV16 serology status in 4,002 oral cancer cases (OPC and oral cavity cancer (OCC)) and 5,256 controls. We detect four susceptibility loci pointing to a distinct genetic predisposition by HPV status. Our most notable finding in the HLA region, that is now confirmed to be specific of HPV(+)OPC risk, reveal two independent loci with strong protective effects, one refining the previously reported HLA class II haplotype association. Antibody levels against HPV16 viral proteins strongly implicate the protective HLA variants as major determinants of humoral response against L1 capsid protein or E6 oncoprotein suggesting a natural immune response against HPV(+)OPC promoted by HLA variants. This indicates that therapeutic vaccines that target E6 and attenuate viral response after established HPV infections might protect against HPV(+)OPC.

The shared genetic architecture between epidemiological and behavioral traits with lung cancer

The complex polygenic nature of lung cancer is not fully characterized. Our study seeks to identify novel phenotypes associated with lung cancer using cross-trait linkage disequilibrium score regression (LDSR). We measured pairwise genetic correlation (rg) and SNP heritability (h2) between 347 traits and lung cancer risk using genome-wide association study summary statistics from the UKBB and OncoArray consortium. Further, we conducted analysis after removing genomic regions previously associated with smoking behaviors to mitigate potential confounding effects. We found significant negative genetic correlations between lung cancer risk and dietary behaviors, fitness metrics, educational attainment, and other psychosocial traits. Alcohol taken with meals (rg = - 0.41, h2 = 0.10, p = 1.33 × 10-16), increased fluid intelligence scores (rg = - 0.25, h2 = 0.22, p = 4.54 × 10-8), and the age at which full time education was completed (rg = - 0.45, h2 = 0.11, p = 1.24 × 10-20) demonstrated negative genetic correlation with lung cancer susceptibility. The body mass index was positively correlated with lung cancer risk (rg = 0.20, h2 = 0.25, p = 2.61 × 10-9). This analysis reveals shared genetic architecture between several traits and lung cancer predisposition. Future work should test for causal relationships and investigate common underlying genetic mechanisms across these genetically correlated traits.

Large-scale cross-cancer fine-mapping of the 5p15.33 region reveals multiple independent signals

Genome-wide association studies (GWASs) have identified thousands of cancer risk loci revealing many risk regions shared across multiple cancers. Characterizing the cross-cancer shared genetic basis can increase our understanding of global mechanisms of cancer development. In this study, we collected GWAS summary statistics based on up to 375,468 cancer cases and 530,521 controls for fourteen types of cancer, including breast (overall, estrogen receptor [ER]-positive, and ER-negative), colorectal, endometrial, esophageal, glioma, head/neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancer, to characterize the shared genetic basis of cancer risk. We identified thirteen pairs of cancers with statistically significant local genetic correlations across eight distinct genomic regions. Specifically, the 5p15.33 region, harboring the TERT and CLPTM1L genes, showed statistically significant local genetic correlations for multiple cancer pairs. We conducted a cross-cancer fine-mapping of the 5p15.33 region based on eight cancers that showed genome-wide significant associations in this region (ER-negative breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, and prostate cancer). We used an iterative analysis pipeline implementing a subset-based meta-analysis approach based on cancer-specific conditional analyses and identified ten independent cross-cancer associations within this region. For each signal, we conducted cross-cancer fine-mapping to prioritize the most plausible causal variants. Our findings provide a more in-depth understanding of the shared inherited basis across human cancers and expand our knowledge of the 5p15.33 region in carcinogenesis.

Sequencing at lymphoid neoplasm susceptibility loci maps six myeloma risk genes

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.

Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology

Bipolar disorder is a heritable mental illness with complex etiology. We performed a genome-wide association study of 41,917 bipolar disorder cases and 371,549 controls of European ancestry, which identified 64 associated genomic loci. Bipolar disorder risk alleles were enriched in genes in synaptic signaling pathways and brain-expressed genes, particularly those with high specificity of expression in neurons of the prefrontal cortex and hippocampus. Significant signal enrichment was found in genes encoding targets of antipsychotics, calcium channel blockers, antiepileptics and anesthetics. Integrating expression quantitative trait locus data implicated 15 genes robustly linked to bipolar disorder via gene expression, encoding druggable targets such as HTR6, MCHR1, DCLK3 and FURIN. Analyses of bipolar disorder subtypes indicated high but imperfect genetic correlation between bipolar disorder type I and II and identified additional associated loci. Together, these results advance our understanding of the biological etiology of bipolar disorder, identify novel therapeutic leads and prioritize genes for functional follow-up studies.

The Shared Genetic Architectures Between Lung Cancer and Multiple Polygenic Phenotypes in Genome-Wide Association Studies

BACKGROUND

Prior genome-wide association studies have identified numerous lung cancer risk loci and reveal substantial etiologic heterogeneity across histologic subtypes. Analyzing the shared genetic architecture underlying variation in complex traits can elucidate common genetic etiologies across phenotypes. Exploring pairwise genetic correlations between lung cancer and other polygenic traits can reveal the common genetic etiology of correlated phenotypes.

METHODS

Using cross-trait linkage disequilibrium score regression, we estimated the pairwise genetic correlation and heritability between lung cancer and multiple traits using publicly available summary statistics. Identified genetic relationships were also examined after excluding genomic regions known to be associated with smoking behaviors, a major risk factor for lung cancer.

RESULTS

We observed several traits showing moderate single nucleotide polymorphism-based heritability and significant genetic correlations with lung cancer. We observed highly significant correlations between the genetic architectures of lung cancer and emphysema/chronic bronchitis across all histologic subtypes, as well as among lung cancer occurring among smokers. Our analyses revealed highly significant positive correlations between lung cancer and paternal history of lung cancer. We also observed a strong negative correlation with parental longevity. We observed consistent directions in genetic patterns after excluding genomic regions associated with smoking behaviors.

CONCLUSIONS

This study identifies numerous phenotypic traits that share genomic architecture with lung carcinogenesis and are not fully accounted for by known smoking-associated genomic loci.

IMPACT

These findings provide new insights into the etiology of lung cancer by identifying traits that are genetically correlated with increased risk of lung cancer.

The PI3K/mTOR Pathway Is Targeted by Rare Germline Variants in Patients with Both Melanoma and Renal Cell Carcinoma

Background: Malignant melanoma and RCC have different embryonic origins, no common lifestyle risk factors but intriguingly share biological properties such as immune regulation and radioresistance. An excess risk of malignant melanoma is observed in RCC patients and vice versa. This bidirectional association is poorly understood, and hypothetic genetic co-susceptibility remains largely unexplored. Results: We hereby provide a clinical and genetic description of a series of 125 cases affected by both malignant melanoma and RCC. Clinical germline mutation testing identified a pathogenic variant in a melanoma and/or RCC predisposing gene in 17/125 cases (13.6%). This included mutually exclusive variants in MITF (p.E318K locus, N = 9 cases), BAP1 (N = 3), CDKN2A (N = 2), FLCN (N = 2), and PTEN (N = 1). A subset of 46 early-onset cases, without underlying germline variation, was whole-exome sequenced. In this series, thirteen genes were significantly enriched in mostly exclusive rare variants predicted to be deleterious, compared to 19,751 controls of similar ancestry. The observed variation mainly consisted of novel or low-frequency variants (<0.01%) within genes displaying strong evolutionary mutational constraints along the PI3K/mTOR pathway, including PIK3CD, NFRKB, EP300, MTOR, and related epigenetic modifier SETD2. The screening of independently processed germline exomes from The Cancer Genome Atlas confirmed an association with melanoma and RCC but not with cancers of established differing etiology such as lung cancers. Conclusions: Our study highlights that an exome-wide case-control enrichment approach may better characterize the rare variant-based missing heritability of multiple primary cancers. In our series, the co-occurrence of malignant melanoma and RCC was associated with germline variation in the PI3K/mTOR signaling cascade, with potential relevance for early diagnostic and clinical management.

Sexual dimorphism in cancer: insights from transcriptional signatures in kidney tissue and renal cell carcinoma

Sexual dimorphism in cancer incidence and outcome is widespread. Understanding the underlying mechanisms is fundamental to improve cancer prevention and clinical management. Sex disparities are particularly striking in kidney cancer: across diverse populations, men consistently show unexplained 2-fold increased incidence and worse prognosis. We have characterized genome-wide expression and regulatory networks of 609 renal tumors and 256 non-tumor renal tissues. Normal kidney displayed sex-specific transcriptional signatures, including higher expression of X-linked tumor suppressor genes in women. Sex-dependent genotype-phenotype associations unraveled women-specific immune regulation. Sex differences were markedly expanded in tumors, with male-biased expression of key genes implicated in metabolism, non-malignant diseases with male predominance and carcinogenesis, including markers of tumor infiltrating leukocytes. Analysis of sex-dependent RCC progression and survival uncovered prognostic markers involved in immune response and oxygen homeostasis. In summary, human kidney tissues display remarkable sexual dimorphism at the molecular level. Sex-specific transcriptional signatures further shape renal cancer, with relevance for clinical management.

Assessing Lung Cancer Absolute Risk Trajectory Based on a Polygenic Risk Model

Lung cancer is the leading cause of cancer-related death globally. An improved risk stratification strategy can increase efficiency of low-dose CT (LDCT) screening. Here we assessed whether individual's genetic background has clinical utility for risk stratification in the context of LDCT screening. On the basis of 13,119 patients with lung cancer and 10,008 controls with European ancestry in the International Lung Cancer Consortium, we constructed a polygenic risk score (PRS) via 10-fold cross-validation with regularized penalized regression. The performance of risk model integrating PRS, including calibration and ability to discriminate, was assessed using UK Biobank data (N = 335,931). Absolute risk was estimated on the basis of age-specific lung cancer incidence and all-cause mortality as competing risk. To evaluate its potential clinical utility, the PRS distribution was simulated in the National Lung Screening Trial (N = 50,772 participants). The lung cancer ORs for individuals at the top decile of the PRS distribution versus those at bottom 10% was 2.39 [95% confidence interval (CI) = 1.92-3.00; P = 1.80 × 10-14] in the validation set (P trend = 5.26 × 10-20). The OR per SD of PRS increase was 1.26 (95% CI = 1.20-1.32; P = 9.69 × 10-23) for overall lung cancer risk in the validation set. When considering absolute risks, individuals at different PRS deciles showed differential trajectories of 5-year and cumulative absolute risk. The age reaching the LDCT screening recommendation threshold can vary by 4 to 8 years, depending on the individual's genetic background, smoking status, and family history. Collectively, these results suggest that individual's genetic background may inform the optimal lung cancer LDCT screening strategy. SIGNIFICANCE: Three large-scale datasets reveal that, after accounting for risk factors, an individual's genetics can affect their lung cancer risk trajectory, thus may inform the optimal timing for LDCT screening.

Genome-wide association meta-analysis identifies pleiotropic risk loci for aerodigestive squamous cell cancers

Squamous cell carcinomas (SqCC) of the aerodigestive tract have similar etiological risk factors. Although genetic risk variants for individual cancers have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. To identify novel and pleotropic SqCC risk variants, we performed a meta-analysis of GWAS data on lung SqCC (LuSqCC), oro/pharyngeal SqCC (OSqCC), laryngeal SqCC (LaSqCC) and esophageal SqCC (ESqCC) cancers, totaling 13,887 cases and 61,961 controls of European ancestry. We identified one novel genome-wide significant (Pmeta<5x10-8) aerodigestive SqCC susceptibility loci in the 2q33.1 region (rs56321285, TMEM273). Additionally, three previously unknown loci reached suggestive significance (Pmeta<5x10-7): 1q32.1 (rs12133735, near MDM4), 5q31.2 (rs13181561, TMEM173) and 19p13.11 (rs61494113, ABHD8). Multiple previously identified loci for aerodigestive SqCC also showed evidence of pleiotropy in at least another SqCC site, these include: 4q23 (ADH1B), 6p21.33 (STK19), 6p21.32 (HLA-DQB1), 9p21.33 (CDKN2B-AS1) and 13q13.1(BRCA2). Gene-based association and gene set enrichment identified a set of 48 SqCC-related genes rel to DNA damage and epigenetic regulation pathways. Our study highlights the importance of cross-cancer analyses to identify pleiotropic risk loci of histology-related cancers arising at distinct anatomical sites.

Patterns of Human Leukocyte Antigen Class I and Class II Associations and Cancer

Human leukocyte antigen (HLA) gene variation is associated with risk of cancers, particularly those with infectious etiology or hematopoietic origin, given its role in immune presentation. Previous studies focused primarily on HLA allele/haplotype-specific associations. To answer whether associations are driven by HLA class I (essential for T-cell cytotoxicity) or class II (important for T-cell helper responses) genes, we analyzed GWAS from 24 case-control studies and consortia comprising 27 cancers (totaling >71,000 individuals). Associations for most cancers with infectious etiology or of hematopoietic origin were driven by multiple HLA regions, suggesting that both cytotoxic and helper T-cell responses are important. Notable exceptions were observed for nasopharyngeal carcinoma, an EBV-associated cancer, and CLL/SLL forms of non-Hodgkin lymphomas; these cancers were associated with HLA class I region only and HLA class II region only, implying the importance of cytotoxic T-cell responses for the former and CD4+ T-cell helper responses for the latter. Our findings suggest that increased understanding of the pattern of HLA associations for individual cancers could lead to better insights into specific mechanisms involved in cancer pathogenesis. SIGNIFICANCE: GWAS of >71,000 individuals across 27 cancer types suggest that patterns of HLA Class I and Class II associations may provide etiologic insights for cancer.

Development of Sensitive Droplet Digital PCR Assays for Detecting Urinary TERT Promoter Mutations as Non-Invasive Biomarkers for Detection of Urothelial Cancer

Somatic mutations in the telomerase reverse transcriptase (TERT) promoter regions are frequent events in urothelial cancer (UC) and their detection in urine (supernatant cell-free DNA or DNA from exfoliated cells) could serve as putative non-invasive biomarkers for UC detection and monitoring. However, detecting these tumor-borne mutations in urine requires highly sensitive methods, capable of measuring low-level mutations. In this study, we developed sensitive droplet digital PCR (ddPCR) assays for detecting TERT promoter mutations (C228T, C228A, CC242-243TT, and C250T). We tested the C228T and C250T ddPCR assays on all samples with sufficient quantity of urinary DNA (urine supernatant cell-free DNA (US cfDNA) or urine pellet cellular DNA (UP cellDNA)) from the DIAGURO (n = 89/93 cases and n = 92/94 controls) and from the IPO-PORTO (n = 49/50 cases and n = 50/50 controls) series that were previously screened with the UroMuTERT assay and compared the performance of the two approaches. In the DIAGURO series, the sensitivity and specificity of the ddPCR assays for detecting UC using either US cfDNA or UP cellDNA were 86.8% and 92.4%. The sensitivity was slightly higher than that of the UroMuTERT assay in the IPO-PORTO series (67.4% vs. 65.3%, respectively), but not in the DIAGURO series (86.8% vs. 90.7%). The specificity was 100% in the IPO-PORTO controls for both the UroMuTERT and ddPCR assays, whereas in the DIAGURO series, the specificity dropped for ddPCR (92.4% versus 95.6%). Overall, an almost perfect agreement between the two methods was observed for both US cfDNA (n = 164; kappa coefficient of 0.91) and UP cellDNA (n = 280; kappa coefficient of 0.94). In a large independent series of serial urine samples from DIAGURO follow-up BC cases (n = 394), the agreement between ddPCR and UroMuTERT was (i) strong (kappa coefficient of 0.87), regardless of urine DNA types (kappa coefficient 0.89 for US cfDNA and 0.85 for UP cellDNA), (ii) the highest for samples with mutant allelic fractions (MAFs) > 2% (kappa coefficient of 0.99) and (iii) only minimal for the samples with the lowest MAFs (< 0.5%; kappa coefficient 0.32). Altogether, our results indicate that the two methods (ddPCR and UroMuTERT) for detecting urinary TERT promoter mutations are comparable and that the discrepancies relate to the detection of low-allelic fraction mutations. The simplicity of the ddPCR assays makes them suitable for implementation in clinical settings.

A molecular map of lung neuroendocrine neoplasms

BACKGROUND

Lung neuroendocrine neoplasms (LNENs) are rare solid cancers, with most genomic studies including a limited number of samples. Recently, generating the first multi-omic dataset for atypical pulmonary carcinoids and the first methylation dataset for large-cell neuroendocrine carcinomas led us to the discovery of clinically relevant molecular groups, as well as a new entity of pulmonary carcinoids (supra-carcinoids).

RESULTS

To promote the integration of LNENs molecular data, we provide here detailed information on data generation and quality control for whole-genome/exome sequencing, RNA sequencing, and EPIC 850K methylation arrays for a total of 84 patients with LNENs. We integrate the transcriptomic data with other previously published data and generate the first comprehensive molecular map of LNENs using the Uniform Manifold Approximation and Projection (UMAP) dimension reduction technique. We show that this map captures the main biological findings of previous studies and can be used as reference to integrate datasets for which RNA sequencing is available. The generated map can be interactively explored and interrogated on the UCSC TumorMap portal (https://tumormap.ucsc.edu/?p=RCG_lungNENomics/LNEN). The data, source code, and compute environments used to generate and evaluate the map as well as the raw data are available, respectively, in a Nextjournal interactive notebook (https://nextjournal.com/rarecancersgenomics/a-molecular-map-of-lung-neuroendocrine-neoplasms/) and at the EMBL-EBI European Genome-phenome Archive and Gene Expression Omnibus data repositories.

CONCLUSIONS

We provide data and all resources needed to integrate them with future LNENs transcriptomic studies, allowing meaningful conclusions to be drawn that will eventually lead to a better understanding of this rare understudied disease.

Dissecting the genetic overlap of smoking behaviors, lung cancer, and chronic obstructive pulmonary disease: A focus on nicotinic receptors and nicotine metabolizing enzyme

Smoking is a major contributor to lung cancer and chronic obstructive pulmonary disease (COPD). Two of the strongest genetic associations of smoking-related phenotypes are the chromosomal regions 15q25.1, encompassing the nicotinic acetylcholine receptor subunit genes CHRNA5-CHRNA3-CHRNB4, and 19q13.2, encompassing the nicotine metabolizing gene CYP2A6. In this study, we examined genetic relations between cigarettes smoked per day, smoking cessation, lung cancer, and COPD. Data consisted of genome-wide association study summary results. Genetic correlations were estimated using linkage disequilibrium score regression software. For each pair of outcomes, z-score-z-score (ZZ) plots were generated. Overall, heavier smoking and decreased smoking cessation showed positive genetic associations with increased lung cancer and COPD risk. The chromosomal region 19q13.2, however, showed a different correlational pattern. For example, the effect allele-C of the sentinel SNP (rs56113850) within CYP2A6 was associated with an increased risk of heavier smoking (z-score = 19.2; p = 1.10 × 10-81 ), lung cancer (z-score = 8.91; p = 5.02 × 10-19 ), and COPD (z-score = 4.04; p = 5.40 × 10-5 ). Surprisingly, this allele-C (rs56113850) was associated with increased smoking cessation (z-score = -8.17; p = 2.52 × 10-26 ). This inverse relationship highlights the need for additional investigation to determine how CYP2A6 variation could increase smoking cessation while also increasing the risk of lung cancer and COPD likely through increased cigarettes smoked per day.

Needlestack: an ultra-sensitive variant caller for multi-sample next generation sequencing data

The emergence of next-generation sequencing (NGS) has revolutionized the way of reaching a genome sequence, with the promise of potentially providing a comprehensive characterization of DNA variations. Nevertheless, detecting somatic mutations is still a difficult problem, in particular when trying to identify low abundance mutations, such as subclonal mutations, tumour-derived alterations in body fluids or somatic mutations from histological normal tissue. The main challenge is to precisely distinguish between sequencing artefacts and true mutations, particularly when the latter are so rare they reach similar abundance levels as artefacts. Here, we present needlestack, a highly sensitive variant caller, which directly learns from the data the level of systematic sequencing errors to accurately call mutations. Needlestack is based on the idea that the sequencing error rate can be dynamically estimated from analysing multiple samples together. We show that the sequencing error rate varies across alterations, illustrating the need to precisely estimate it. We evaluate the performance of needlestack for various types of variations, and we show that needlestack is robust among positions and outperforms existing state-of-the-art method for low abundance mutations. Needlestack, along with its source code is freely available on the GitHub platform: https://github.com/IARCbioinfo/needlestack.

Protein-altering germline mutations implicate novel genes related to lung cancer development

Few germline mutations are known to affect lung cancer risk. We performed analyses of rare variants from 39,146 individuals of European ancestry and investigated gene expression levels in 7,773 samples. We find a large-effect association with an ATM L2307F (rs56009889) mutation in adenocarcinoma for discovery (adjusted Odds Ratio = 8.82, P = 1.18 × 10-15) and replication (adjusted OR = 2.93, P = 2.22 × 10-3) that is more pronounced in females (adjusted OR = 6.81 and 3.19 and for discovery and replication). We observe an excess loss of heterozygosity in lung tumors among ATM L2307F allele carriers. L2307F is more frequent (4%) among Ashkenazi Jewish populations. We also observe an association in discovery (adjusted OR = 2.61, P = 7.98 × 10-22) and replication datasets (adjusted OR = 1.55, P = 0.06) with a loss-of-function mutation, Q4X (rs150665432) of an uncharacterized gene, KIAA0930. Our findings implicate germline genetic variants in ATM with lung cancer susceptibility and suggest KIAA0930 as a novel candidate gene for lung cancer risk.

Circulating tumour-derived KRAS mutations in pancreatic cancer cases are predominantly carried by very short fragments of cell-free DNA

BACKGROUND

The DNA released into the bloodstream by malignant tumours· called circulating tumour DNA (ctDNA), is often a small fraction of total cell-free DNA shed predominantly by hematopoietic cells and is therefore challenging to detect. Understanding the biological properties of ctDNA is key to the investigation of its clinical relevance as a non-invasive marker for cancer detection and monitoring.

METHODS

We selected 40 plasma DNA samples of pancreatic cancer cases previously reported to carry a KRAS mutation at the 'hotspot' codon 12 and re-screened the cell-free DNA using a 4-size amplicons strategy (57 bp, 79 bp, 167 bp and 218 bp) combined with ultra-deep sequencing in order to investigate whether amplicon lengths could impact on the capacity of detection of ctDNA, which in turn could provide inference of ctDNA and non-malignant cell-free DNA size distribution.

FINDINGS

Higher KRAS amplicon size (167 bp and 218 bp) was associated with lower detectable cell-free DNA mutant allelic fractions (p < 0·0001), with up to 4·6-fold (95% CI: 2·6-8·1) difference on average when comparing the 218bp- and the 57bp-amplicons. The proportion of cases with detectable KRAS mutations was also hampered with increased amplicon lengths, with only half of the cases having detectable ctDNA using the 218 bp assay relative to those detected with amplicons less than 80 bp.

INTERPRETATION

Tumour-derived mutations are carried by shorter cell-free DNA fragments than fragments of wild-type allele. Targeting short amplicons increases the sensitivity of cell-free DNA assays for pancreatic cancer and should be taken into account for optimized assay design and for evaluating their clinical performance.

FUNDING

IARC; MH CZ - DRO; MH SK; exchange program between IARC and Sao Paulo medical Sciences; French Cancer League.

Lipid Trait Variants and the Risk of Non-Hodgkin Lymphoma Subtypes: A Mendelian Randomization Study

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

Transcriptome-wide association study reveals candidate causal genes for lung cancer

We have recently completed the largest GWAS on lung cancer including 29,266 cases and 56,450 controls of European descent. The goal of our study has been to integrate the complete GWAS results with a large-scale expression quantitative trait loci (eQTL) mapping study in human lung tissues (n = 1,038) to identify candidate causal genes for lung cancer. We performed transcriptome-wide association study (TWAS) for lung cancer overall, by histology (adenocarcinoma, squamous cell carcinoma and small cell lung cancer) and smoking subgroups (never- and ever-smokers). We performed replication analysis using lung data from the Genotype-Tissue Expression (GTEx) project. DNA damage assays were performed in human lung fibroblasts for selected TWAS genes. As expected, the main TWAS signal for all histological subtypes and ever-smokers was on chromosome 15q25. The gene most strongly associated with lung cancer at this locus using the TWAS approach was IREB2 (pTWAS = 1.09E-99), where lower predicted expression increased lung cancer risk. A new lung adenocarcinoma susceptibility locus was revealed on 9p13.3 and associated with higher predicted expression of AQP3 (pTWAS = 3.72E-6). Among the 45 previously described lung cancer GWAS loci, we mapped candidate target gene for 17 of them. The association AQP3-adenocarcinoma on 9p13.3 was replicated using GTEx (pTWAS = 6.55E-5). Consistent with the effect of risk alleles on gene expression levels, IREB2 knockdown and AQP3 overproduction promote endogenous DNA damage. These findings indicate genes whose expression in lung tissue directly influences lung cancer risk.

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

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

Sex specific associations in genome wide association analysis of renal cell carcinoma

Renal cell carcinoma (RCC) has an undisputed genetic component and a stable 2:1 male to female sex ratio in its incidence across populations, suggesting possible sexual dimorphism in its genetic susceptibility. We conducted the first sex-specific genome-wide association analysis of RCC for men (3227 cases, 4916 controls) and women (1992 cases, 3095 controls) of European ancestry from two RCC genome-wide scans and replicated the top findings using an additional series of men (2261 cases, 5852 controls) and women (1399 cases, 1575 controls) from two independent cohorts of European origin. Our study confirmed sex-specific associations for two known RCC risk loci at 14q24.2 (DPF3) and 2p21(EPAS1). We also identified two additional suggestive male-specific loci at 6q24.3 (SAMD5, male odds ratio (ORmale) = 0.83 [95% CI = 0.78-0.89], Pmale = 1.71 × 10-8 compared with female odds ratio (ORfemale) = 0.98 [95% CI = 0.90-1.07], Pfemale = 0.68) and 12q23.3 (intergenic, ORmale = 0.75 [95% CI = 0.68-0.83], Pmale = 1.59 × 10-8 compared with ORfemale = 0.93 [95% CI = 0.82-1.06], Pfemale = 0.21) that attained genome-wide significance in the joint meta-analysis. Herein, we provide evidence of sex-specific associations in RCC genetic susceptibility and advocate the necessity of larger genetic and genomic studies to unravel the endogenous causes of sex bias in sexually dimorphic traits and diseases like RCC.

Genetic overlap between autoimmune diseases and non-Hodgkin lymphoma subtypes

Epidemiologic studies show an increased risk of non-Hodgkin lymphoma (NHL) in patients with autoimmune disease (AD), due to a combination of shared environmental factors and/or genetic factors, or a causative cascade: chronic inflammation/antigen-stimulation in one disease leads to another. Here we assess shared genetic risk in genome-wide-association-studies (GWAS). Secondary analysis of GWAS of NHL subtypes (chronic lymphocytic leukemia, diffuse large B-cell lymphoma, follicular lymphoma, and marginal zone lymphoma) and ADs (rheumatoid arthritis, systemic lupus erythematosus, and multiple sclerosis). Shared genetic risk was assessed by (a) description of regional genetic of overlap, (b) polygenic risk score (PRS), (c)"diseasome", (d)meta-analysis. Descriptive analysis revealed few shared genetic factors between each AD and each NHL subtype. The PRS of ADs were not increased in NHL patients (nor vice versa). In the diseasome, NHLs shared more genetic etiology with ADs than solid cancers (p = .0041). A meta-analysis (combing AD with NHL) implicated genes of apoptosis and telomere length. This GWAS-based analysis four NHL subtypes and three ADs revealed few weakly-associated shared loci, explaining little total risk. This suggests common genetic variation, as assessed by GWAS in these sample sizes, may not be the primary explanation for the link between these ADs and NHLs.

Redefining malignant pleural mesothelioma types as a continuum uncovers immune-vascular interactions

BACKGROUND

Malignant Pleural Mesothelioma (MPM) is an aggressive disease related to asbestos exposure, with no effective therapeutic options.

METHODS

We undertook unsupervised analyses of RNA-sequencing data of 284 MPMs, with no assumption of discreteness. Using immunohistochemistry, we performed an orthogonal validation on a subset of 103 samples and a biological replication in an independent series of 77 samples.

FINDINGS

A continuum of molecular profiles explained the prognosis of the disease better than any discrete model. The immune and vascular pathways were the major sources of molecular variation, with strong differences in the expression of immune checkpoints and pro-angiogenic genes; the extrema of this continuum had specific molecular profiles: a "hot" bad-prognosis profile, with high lymphocyte infiltration and high expression of immune checkpoints and pro-angiogenic genes; a "cold" bad-prognosis profile, with low lymphocyte infiltration and high expression of pro-angiogenic genes; and a "VEGFR2+/VISTA+" better-prognosis profile, with high expression of immune checkpoint VISTA and pro-angiogenic gene VEGFR2. We validated the gene expression levels at the protein level for a subset of five selected genes belonging to the immune and vascular pathways (CD8A, PDL1, VEGFR3, VEGFR2, and VISTA), in the validation series, and replicated the molecular profiles as well as their prognostic value in the replication series.

INTERPRETATION

The prognosis of MPM is best explained by a continuous model, which extremes show specific expression patterns of genes involved in angiogenesis and immune response.

Integrative and comparative genomic analyses identify clinically relevant pulmonary carcinoid groups and unveil the supra-carcinoids

The worldwide incidence of pulmonary carcinoids is increasing, but little is known about their molecular characteristics. Through machine learning and multi-omics factor analysis, we compare and contrast the genomic profiles of 116 pulmonary carcinoids (including 35 atypical), 75 large-cell neuroendocrine carcinomas (LCNEC), and 66 small-cell lung cancers. Here we report that the integrative analyses on 257 lung neuroendocrine neoplasms stratify atypical carcinoids into two prognostic groups with a 10-year overall survival of 88% and 27%, respectively. We identify therapeutically relevant molecular groups of pulmonary carcinoids, suggesting DLL3 and the immune system as candidate therapeutic targets; we confirm the value of OTP expression levels for the prognosis and diagnosis of these diseases, and we unveil the group of supra-carcinoids. This group comprises samples with carcinoid-like morphology yet the molecular and clinical features of the deadly LCNEC, further supporting the previously proposed molecular link between the low- and high-grade lung neuroendocrine neoplasms.

Safety and Efficacy of Belimumab Plus Standard Therapy for Up to Thirteen Years in Patients With Systemic Lupus Erythematosus

Objective

To investigate the long‐term safety and efficacy of intravenous (IV) belimumab plus standard of care (SOC) therapy for systemic lupus erythematosus (SLE) in patients with active, autoantibody‐positive SLE.

Methods

The study was designed as a multicenter, open‐label, continuation study of IV belimumab given every 4 weeks in conjunction with SOC therapy in patients with SLE who completed a phase II, double‐blind study. Adverse events (AEs) and laboratory data were monitored from the first belimumab dose (in either study) until 24 weeks after the final dose. Efficacy assessments included SLE Responder Index (SRI) and flare index scores (each assessed at 16‐week intervals) and glucocorticoid use (assessed at 4‐week intervals).

Results

Of the 476 patients in the parent study, 298 (62.6%) entered the continuation study, of whom 96 (32.2%) remained in the study. Patients received belimumab for up to 13 years (median duration of exposure 3,334.0 days [range 260–4,332 days], total belimumab exposure 2,294 patient‐years, median number of infusions 115.5 [range 7–155]). The percentage of patients with AEs each year remained stable or decreased. Normal serum IgG levels were maintained in the majority of patients over the study, and the rate of infections remained stable. The percentage of patients who achieved an SRI response increased from 32.8% (year 1) to 75.6% of those remaining on treatment at year 12. The glucocorticoid dose was decreased in patients who had been receiving >7.5 mg/day at baseline.

Conclusion

This study is the longest to date to assess belimumab treatment in patients with SLE in clinical trials. Belimumab was well tolerated with no new safety concerns, and efficacy was maintained in patients who continued the study. For patients who initially exhibited a satisfactory response to belimumab, the treatment continues to be well tolerated and provides long‐term disease control.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Genome-wide association study identifies 30 loci associated with bipolar disorder

Bipolar disorder is a highly heritable psychiatric disorder. We performed a genome-wide association study (GWAS) including 20,352 cases and 31,358 controls of European descent, with follow-up analysis of 822 variants with P < 1 × 10-4 in an additional 9,412 cases and 137,760 controls. Eight of the 19 variants that were genome-wide significant (P < 5 × 10-8) in the discovery GWAS were not genome-wide significant in the combined analysis, consistent with small effect sizes and limited power but also with genetic heterogeneity. In the combined analysis, 30 loci were genome-wide significant, including 20 newly identified loci. The significant loci contain genes encoding ion channels, neurotransmitter transporters and synaptic components. Pathway analysis revealed nine significantly enriched gene sets, including regulation of insulin secretion and endocannabinoid signaling. Bipolar I disorder is strongly genetically correlated with schizophrenia, driven by psychosis, whereas bipolar II disorder is more strongly correlated with major depressive disorder. These findings address key clinical questions and provide potential biological mechanisms for bipolar disorder.

Genetic interaction analysis among oncogenesis-related genes revealed novel genes and networks in lung cancer development

The development of cancer is driven by the accumulation of many oncogenesis-related genetic alterations and tumorigenesis is triggered by complex networks of involved genes rather than independent actions. To explore the epistasis existing among oncogenesis-related genes in lung cancer development, we conducted pairwise genetic interaction analyses among 35,031 SNPs from 2027 oncogenesis-related genes. The genotypes from three independent genome-wide association studies including a total of 24,037 lung cancer patients and 20,401 healthy controls with Caucasian ancestry were analyzed in the study. Using a two-stage study design including discovery and replication studies, and stringent Bonferroni correction for multiple statistical analysis, we identified significant genetic interactions between SNPs in RGL1:RAD51B (OR=0.44, p value=3.27x10-11 in overall lung cancer and OR=0.41, p value=9.71x10-11 in non-small cell lung cancer), SYNE1:RNF43 (OR=0.73, p value=1.01x10-12 in adenocarcinoma) and FHIT:TSPAN8 (OR=1.82, p value=7.62x10-11 in squamous cell carcinoma) in our analysis. None of these genes have been identified from previous main effect association studies in lung cancer. Further eQTL gene expression analysis in lung tissues provided information supporting the functional role of the identified epistasis in lung tumorigenesis. Gene set enrichment analysis revealed potential pathways and gene networks underlying molecular mechanisms in overall lung cancer as well as histology subtypes development. Our results provide evidence that genetic interactions between oncogenesis-related genes play an important role in lung tumorigenesis and epistasis analysis, combined with functional annotation, provides a valuable tool for uncovering functional novel susceptibility genes that contribute to lung cancer development by interacting with other modifier genes.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

IMPACT

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

The influence of obesity-related factors in the etiology of renal cell carcinoma-A mendelian randomization study

BACKGROUND

Several obesity-related factors have been associated with renal cell carcinoma (RCC), but it is unclear which individual factors directly influence risk. We addressed this question using genetic markers as proxies for putative risk factors and evaluated their relation to RCC risk in a mendelian randomization (MR) framework. This methodology limits bias due to confounding and is not affected by reverse causation.

METHODS AND FINDINGS

Genetic markers associated with obesity measures, blood pressure, lipids, type 2 diabetes, insulin, and glucose were initially identified as instrumental variables, and their association with RCC risk was subsequently evaluated in a genome-wide association study (GWAS) of 10,784 RCC patients and 20,406 control participants in a 2-sample MR framework. The effect on RCC risk was estimated by calculating odds ratios (ORSD) for a standard deviation (SD) increment in each risk factor. The MR analysis indicated that higher body mass index increases the risk of RCC (ORSD: 1.56, 95% confidence interval [CI] 1.44-1.70), with comparable results for waist-to-hip ratio (ORSD: 1.63, 95% CI 1.40-1.90) and body fat percentage (ORSD: 1.66, 95% CI 1.44-1.90). This analysis further indicated that higher fasting insulin (ORSD: 1.82, 95% CI 1.30-2.55) and diastolic blood pressure (DBP; ORSD: 1.28, 95% CI 1.11-1.47), but not systolic blood pressure (ORSD: 0.98, 95% CI 0.84-1.14), increase the risk for RCC. No association with RCC risk was seen for lipids, overall type 2 diabetes, or fasting glucose.

CONCLUSIONS

This study provides novel evidence for an etiological role of insulin in RCC, as well as confirmatory evidence that obesity and DBP influence RCC risk.

Multi-Omics Analysis Reveals a HIF Network and Hub Gene EPAS1 Associated with Lung Adenocarcinoma

Recent technological advancements have permitted high-throughput measurement of the human genome, epigenome, metabolome, transcriptome, and proteome at the population level. We hypothesized that subsets of genes identified from omic studies might have closely related biological functions and thus might interact directly at the network level. Therefore, we conducted an integrative analysis of multi-omic datasets of non-small cell lung cancer (NSCLC) to search for association patterns beyond the genome and transcriptome. A large, complex, and robust gene network containing well-known lung cancer-related genes, including EGFR and TERT, was identified from combined gene lists for lung adenocarcinoma. Members of the hypoxia-inducible factor (HIF) gene family were at the center of this network. Subsequent sequencing of network hub genes within a subset of samples from the Transdisciplinary Research in Cancer of the Lung-International Lung Cancer Consortium (TRICL-ILCCO) consortium revealed a SNP (rs12614710) in EPAS1 associated with NSCLC that reached genome-wide significance (OR = 1.50; 95% CI: 1.31-1.72; p = 7.75 × 10^-9). Using imputed data, we found that this SNP remained significant in the entire TRICL-ILCCO consortium (p = .03). Additional functional studies are warranted to better understand interrelationships among genetic polymorphisms, DNA methylation status, and EPAS1 expression.

Germline Lysine-Specific Demethylase 1 (LSD1/KDM1A) Mutations Confer Susceptibility to Multiple Myeloma

Given the frequent and largely incurable occurrence of multiple myeloma, identification of germline genetic mutations that predispose cells to multiple myeloma may provide insight into disease etiology and the developmental mechanisms of its cell of origin, the plasma cell (PC). Here, we identified familial and early-onset multiple myeloma kindreds with truncating mutations in lysine-specific demethylase 1 (LSD1/KDM1A), an epigenetic transcriptional repressor that primarily demethylates histone H3 on lysine 4 and regulates hematopoietic stem cell self-renewal. In addition, we found higher rates of germline truncating and predicted deleterious missense KDM1A mutations in patients with multiple myeloma unselected for family history compared with controls. Both monoclonal gammopathy of undetermined significance (MGUS) and multiple myeloma cells have significantly lower KDM1A transcript levels compared with normal PCs. Transcriptome analysis of multiple myeloma cells from KDM1A mutation carriers shows enrichment of pathways and MYC target genes previously associated with myeloma pathogenesis. In mice, antigen challenge followed by pharmacologic inhibition of KDM1A promoted PC expansion, enhanced secondary immune response, elicited appearance of serum paraprotein, and mediated upregulation of MYC transcriptional targets. These changes are consistent with the development of MGUS. Collectively, our findings show that KDM1A is the first autosomal-dominant multiple myeloma germline predisposition gene providing new insights into its mechanistic roles as a tumor suppressor during post-germinal center B-cell differentiation.Significance: KDM1A is the first germline autosomal dominant predisposition gene identified in multiple myeloma and provides new insights into multiple myeloma etiology and the mechanistic role of KDM1A as a tumor suppressor during post-germinal center B-cell differentiation. Cancer Res; 78(10); 2747-59. ©2018 AACR.

Genetic Contributions to The Association Between Adult Height and Head and Neck Cancer: A Mendelian Randomization Analysis

With the aim to dissect the effect of adult height on head and neck cancer (HNC), we use the Mendelian randomization (MR) approach to test the association between genetic instruments for height and the risk of HNC. 599 single nucleotide polymorphisms (SNPs) were identified as genetic instruments for height, accounting for 16% of the phenotypic variation. Genetic data concerning HNC cases and controls were obtained from a genome-wide association study. Summary statistics for genetic association were used in complementary MR approaches: the weighted genetic risk score (GRS) and the inverse-variance weighted (IVW). MR-Egger regression was used for sensitivity analysis and pleiotropy evaluation. From the GRS analysis, one standard deviation (SD) higher height (6.9 cm; due to genetic predisposition across 599 SNPs) raised the risk for HNC (Odds ratio (OR), 1.14; 95% Confidence Interval (95%CI), 0.99-1.32). The association analyses with potential confounders revealed that the GRS was associated with tobacco smoking (OR = 0.80, 95% CI (0.69-0.93)). MR-Egger regression did not provide evidence of overall directional pleiotropy. Our study indicates that height is potentially associated with HNC risk. However, the reported risk could be underestimated since, at the genetic level, height emerged to be inversely associated with smoking.

Genome-wide interaction study of smoking behavior and non-small cell lung cancer risk in Caucasian population

Non-small cell lung cancer is the most common type of lung cancer. Both environmental and genetic risk factors contribute to lung carcinogenesis. We conducted a genome-wide interaction analysis between single nucleotide polymorphisms (SNPs) and smoking status (never- versus ever-smokers) in a European-descent population. We adopted a two-step analysis strategy in the discovery stage: we first conducted a case-only interaction analysis to assess the relationship between SNPs and smoking behavior using 13336 non-small cell lung cancer cases. Candidate SNPs with P-value <0.001 were further analyzed using a standard case-control interaction analysis including 13970 controls. The significant SNPs with P-value <3.5 × 10-5 (correcting for multiple tests) from the case-control analysis in the discovery stage were further validated using an independent replication dataset comprising 5377 controls and 3054 non-small cell lung cancer cases. We further stratified the analysis by histological subtypes. Two novel SNPs, rs6441286 and rs17723637, were identified for overall lung cancer risk. The interaction odds ratio and meta-analysis P-value for these two SNPs were 1.24 with 6.96 × 10-7 and 1.37 with 3.49 × 10-7, respectively. In addition, interaction of smoking with rs4751674 was identified in squamous cell lung carcinoma with an odds ratio of 0.58 and P-value of 8.12 × 10-7. This study is by far the largest genome-wide SNP-smoking interaction analysis reported for lung cancer. The three identified novel SNPs provide potential candidate biomarkers for lung cancer risk screening and intervention. The results from our study reinforce that gene-smoking interactions play important roles in the etiology of lung cancer and account for part of the missing heritability of this disease.

Gene set enrichment analysis and expression pattern exploration implicate an involvement of neurodevelopmental processes in bipolar disorder

BACKGROUND

Bipolar disorder (BD) is a common and highly heritable disorder of mood. Genome-wide association studies (GWAS) have identified several independent susceptibility loci. In order to extract more biological information from GWAS data, multi-locus approaches represent powerful tools since they utilize knowledge about biological processes to integrate functional sets of genes at strongly to moderately associated loci.

METHODS

We conducted gene set enrichment analyses (GSEA) using 2.3 million single-nucleotide polymorphisms, 397 Reactome pathways and 24,025 patients with BD and controls. RNA expression of implicated individual genes and gene sets were examined in post-mortem brains across lifespan.

RESULTS

Two pathways showed a significant enrichment after correction for multiple comparisons in the GSEA: GRB2 events in ERBB2 signaling, for which 6 of 21 genes were BD associated (PFDR = 0.0377), and NCAM signaling for neurite out-growth, for which 11 out of 62 genes were BD associated (PFDR = 0.0451). Most pathway genes showed peaks of RNA co-expression during fetal development and infancy and mapped to neocortical areas and parts of the limbic system.

LIMITATIONS

Pathway associations were technically reproduced by two methods, although they were not formally replicated in independent samples. Gene expression was explored in controls but not in patients.

CONCLUSIONS

Pathway analysis in large GWAS data of BD and follow-up of gene expression patterns in healthy brains provide support for an involvement of neurodevelopmental processes in the etiology of this neuropsychiatric disease. Future studies are required to further evaluate the relevance of the implicated genes on pathway functioning and clinical aspects of BD.

Effect of alcohol dehydrogenase-1B and -7 polymorphisms on blood ethanol and acetaldehyde concentrations in healthy subjects with a history of moderate alcohol consumption

This study aims to evaluate the effect of ADH1B and ADH7 genotypes on blood acetaldehyde and ethanol levels after alcohol ingestion, and to measure the genotoxic effect of smoking and ethanol on the buccal cells, also controlling for ADH variants. We recruited healthy Italian subjects with at least a moderate history of alcohol consumption. All subjects were given an alcoholic drink of 0.4 g ethanol /kg of body weight. Blood venous samples were collected at baseline, and 30, 60, 90, and 120 minutes after ingestion. Buccal cells were collected before ethanol ingestion. Sixty subjects were enrolled in the study. Individuals with the ADH1B GG genotype had median ethanol levels of 5.0mM (IQR 3.4-7.2), and those with the ADH1B GT/TT genotype had 4.7mM (IQR 4.2-4.8). Corresponding acetaldehyde levels were 1.5μM (IQR 0.7-2.6) for ADH1B GG genotype and 1.6μM (IQR 1.5-1.7) for ADH1B CG/GG genotype. Individuals with the ADH7 CC genotype had median ethanol levels of 5.0mM (IQR 3.3-7.2), while 5.0mM (IQR 4.7-5.6) was in those with the ADH7 CG/GG genotype. Corresponding acetaldehyde levels were 1.5 μM (IQR 0.7-2.6) for ADH7 CC genotype and 1.5 μM (IQR 1.4-1.6) for ADH7 CG/GG genotypes. A non-significant increase in the frequency of karyolitic and pyknotic cells was found in the group of heavy drinkers and current smokers, when compared to the moderate drinkers and the non-smokers. Our study does not support the hypothesis that ADH1B and ADH7 genotypes affect blood ethanol and acetaldehyde concentration.