Candidate pathway analysis of surfactant proteins identifies CTSH and SFTA2 that influences lung cancer risk.

Pulmonary surfactant is a lipoprotein synthesized and secreted by alveolar type II cells in lung. We evaluated the associations between 200,139 single nucleotide polymorphisms (SNPs) of 40 surfactant-related genes and lung cancer risk using genotyped data from two independent lung cancer genome-wide association studies. Discovery data included 18,082 cases and 13,780 controls of European ancestry. Replication data included 1,914 cases and 3,065 controls of European descent. Using multivariate logistic regression, we found novel SNPs in surfactant-related genes CTSH [rs34577742 C > T, odds ratio (OR) = 0.90, 95% confidence interval (CI) = 0.89-0.93, P = 7.64 × 10-9] and SFTA2 (rs3095153 G > A, OR = 1.16, 95% CI = 1.10-1.21, P = 1.27 × 10-9) associated with overall lung cancer in the discovery data and validated in an independent replication data-CTSH (rs34577742 C > T, OR = 0.88, 95% CI = 0.80-0.96, P = 5.76 × 10-3) and SFTA2 (rs3095153 G > A, OR = 1.14, 95% CI = 1.01-1.28, P = 3.25 × 10-2). Among ever smokers, we found SNPs in CTSH (rs34577742 C > T, OR = 0.89, 95% CI = 0.85-0.92, P = 1.94 × 10-7) and SFTA2 (rs3095152 G > A, OR = 1.20, 95% CI = 1.14-1.27, P = 4.25 × 10-11) associated with overall lung cancer in the discovery data and validated in the replication data-CTSH (rs34577742 C > T, OR = 0.88, 95% CI = 0.79-0.97, P = 1.64 × 10-2) and SFTA2 (rs3095152 G > A, OR = 1.15, 95% CI = 1.01-1.30, P = 3.81 × 10-2). Subsequent transcriptome-wide association study using expression weights from a lung expression quantitative trait loci study revealed genes most strongly associated with lung cancer are CTSH (PTWAS = 2.44 × 10-4) and SFTA2 (PTWAS = 2.32 × 10-6).

Genome-wide interaction analysis identified low-frequency variants with sex disparity in lung cancer risk.

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.

Prevalent occupational exposures and risk of lung cancer among women: Results from the application of the Canadian Job-Exposure Matrix (CANJEM) to a combined set of ten case-control studies.

BACKGROUND: Worldwide, lung cancer is the second leading cause of cancer death in women. The present study explored associations between occupational exposures that are prevalent among women, and lung cancer. METHODS: Data from 10 case-control studies of lung cancer from Europe, Canada, and New Zealand conducted between 1988 and 2008 were combined. Lifetime occupational history and information on nonoccupational factors including smoking were available for 3040 incident lung cancer cases and 4187 controls. We linked each reported job to the Canadian Job-Exposure Matrix (CANJEM), which provided estimates of probability, intensity, and frequency of exposure to each selected agent in each job. For this analysis, we selected 15 agents (cleaning agents, biocides, cotton dust, synthetic fibers, formaldehyde, cooking fumes, organic solvents, cellulose, polycyclic aromatic hydrocarbons from petroleum, ammonia, metallic dust, alkanes C18+, iron compounds, isopropanol, and calcium carbonate) that had lifetime exposure prevalence of at least 5% in the combined study population. For each agent, we estimated lung cancer risk in each study center for ever-exposure, by duration of exposure, and by cumulative exposure, using separate logistic regression models adjusted for smoking and other covariates. We then estimated the meta-odds ratios using random-effects meta-analysis. RESULTS AND CONCLUSIONS: None of the agents assessed showed consistent and compelling associations with lung cancer among women. The following agents showed elevated odds ratio in some analyses: metallic dust, iron compounds, isopropanol, and organic solvents. Future research into occupational lung cancer risk factors among women should prioritize these agents.

Integration of multiomic annotation data to prioritize and characterize inflammation and immune-related risk variants in squamous cell lung cancer.

Clinical trial results have recently demonstrated that inhibiting inflammation by targeting the interleukin-1ß pathway can offer a significant reduction in lung cancer incidence and mortality, highlighting a pressing and unmet need to understand the benefits of inflammation-focused lung cancer therapies at the genetic level. While numerous genome-wide association studies (GWAS) have explored the genetic etiology of lung cancer, there remains a large gap between the type of information that may be gleaned from an association study and the depth of understanding necessary to explain and drive translational findings. Thus, in this study we jointly model and integrate extensive multiomics data sources, utilizing a total of 40 genome-wide functional annotations that augment previously published results from the International Lung Cancer Consortium (ILCCO) GWAS, to prioritize and characterize single nucleotide polymorphisms (SNPs) that increase risk of squamous cell lung cancer through the inflammatory and immune responses. Our work bridges the gap between correlative analysis and translational follow-up research, refining GWAS association measures in an interpretable and systematic manner. In particular, reanalysis of the ILCCO data highlights the impact of highly associated SNPs from nuclear factor-κB signaling pathway genes as well as major histocompatibility complex mediated variation in immune responses. One consequence of prioritizing likely functional SNPs is the pruning of variants that might be selected for follow-up work by over an order of magnitude, from potentially tens of thousands to hundreds. The strategies we introduce provide informative and interpretable approaches for incorporating extensive genome-wide annotation data in analysis of genetic association studies.

A Likelihood Ratio Test for Gene-Environment Interaction Based on the Trend Effect of Genotype Under an Additive Risk Model Using the Gene-Environment Independence Assumption.

Several statistical methods have been proposed for testing gene-environment (G-E) interactions under additive risk models using data from genome-wide association studies. However, these approaches have strong assumptions from underlying genetic models, such as dominant or recessive effects that are known to be less robust when the true genetic model is unknown. We aimed to develop a robust trend test employing a likelihood ratio test for detecting G-E interaction under an additive risk model, while incorporating the G-E independence assumption to increase power. We used a constrained likelihood to impose 2 sets of constraints for: 1) the linear trend effect of genotype and 2) the additive joint effects of gene and environment. To incorporate the G-E independence assumption, a retrospective likelihood was used versus a standard prospective likelihood. Numerical investigation suggests that the proposed tests are more powerful than tests assuming dominant, recessive, or general models under various parameter settings and under both likelihoods. Incorporation of the independence assumption enhances efficiency by 2.5-fold. We applied the proposed methods to examine the gene-smoking interaction for lung cancer and gene-apolipoprotein E $\varepsilon$4 interaction for Alzheimer disease, which identified 2 interactions between apolipoprotein E $\varepsilon$4 and loci membrane-spanning 4-domains subfamily A (MS4A) and bridging integrator 1 (BIN1) genes at genome-wide significance that were replicated using independent data.

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.

Genome-Wide Meta-Analyses of FTND and TTFC Phenotypes.

INTRODUCTION: FTND (FagerstrÓ§m test for nicotine dependence) and TTFC (time to smoke first cigarette in the morning) are common measures of nicotine dependence (ND). However, genome-wide meta-analysis for these phenotypes has not been reported. METHODS: Genome-wide meta-analyses for FTND (N = 19,431) and TTFC (N = 18,567) phenotypes were conducted for adult smokers of European ancestry from 14 independent cohorts. RESULTS: We found that SORBS2 on 4q35 (p = 4.05 × 10-8), BG182718 on 11q22 (p = 1.02 × 10-8), and AA333164 on 14q21 (p = 4.11 × 10-9) were associated with TTFC phenotype. We attempted replication of leading candidates with independent samples (FTND, N = 7010 and TTFC, N = 10 061), however, due to limited power of the replication samples, the replication of these new loci did not reach significance. In gene-based analyses, COPB2 was found associated with FTND phenotype, and TFCP2L1, RELN, and INO80C were associated with TTFC phenotype. In pathway and network analyses, we found that the interconnected interactions among the endocytosis, regulation of actin cytoskeleton, axon guidance, MAPK signaling, and chemokine signaling pathways were involved in ND. CONCLUSIONS: Our analyses identified several promising candidates for both FTND and TTFC phenotypes, and further verification of these candidates was necessary. Candidates supported by both FTND and TTFC (CHRNA4, THSD7B, RBFOX1, and ZNF804A) were associated with addiction to alcohol, cocaine, and heroin, and were associated with autism and schizophrenia. We also identified novel pathways involved in cigarette smoking. The pathway interactions highlighted the importance of receptor recycling and internalization in ND. IMPLICATIONS: Understanding the genetic architecture of cigarette smoking and ND is critical to develop effective prevention and treatment. Our study identified novel candidates and biological pathways involved in FTND and TTFC phenotypes, and this will facilitate further investigation of these candidates and pathways.

Meta-analysis of genome-wide association studies identifies multiple lung cancer susceptibility loci in never-smoking Asian women.

Genome-wide association studies (GWAS) of lung cancer in Asian never-smoking women have previously identified six susceptibility loci associated with lung cancer risk. To further discover new susceptibility loci, we imputed data from four GWAS of Asian non-smoking female lung cancer (6877 cases and 6277 controls) using the 1000 Genomes Project (Phase 1 Release 3) data as the reference and genotyped additional samples (5878 cases and 7046 controls) for possible replication. In our meta-analysis, three new loci achieved genome-wide significance, marked by single nucleotide polymorphism (SNP) rs7741164 at 6p21.1 (per-allele odds ratio (OR) = 1.17; P = 5.8 × 10(-13)), rs72658409 at 9p21.3 (per-allele OR = 0.77; P = 1.41 × 10(-10)) and rs11610143 at 12q13.13 (per-allele OR = 0.89; P = 4.96 × 10(-9)). These findings identified new genetic susceptibility alleles for lung cancer in never-smoking women in Asia and merit follow-up to understand their biological underpinnings.

Germline mutations in Protection of Telomeres 1 in two families with Hodgkin lymphoma.

In a previous whole exome sequencing of patients from 41 families with Hodgkin lymphoma, we identified two families with distinct heterozygous rare coding variants in POT1 (D224N and Y36H), both in a highly conserved region of the gene. POT1 D224N mutant did not bind to a single-stranded telomere oligonucleotide in vitro suggesting the mutation perturbs POT1's ability to bind to the telomeric G-rich overhang. Human HT1080 cells expressing POT1 D224N and lymphoblastoid cells carrying Y36H both showed increased telomere length and fragility in comparison to wild type cells. This strongly suggests that mutant POT1 causes chromosome instability and may play a role in lymphomagenesis in these families.

Functional variants in DCAF4 associated with lung cancer risk in European populations.

Cullin-RING ubiquitin ligases (CRLs) responsible for substrate specificity of ubiquitination play a key role in cell-cycle control and DNA damage response. In this study, we assessed associations between 16 599 SNPs in 115 CRL genes and lung cancer risk by using summary data of six published genome-wide association studies (GWASs) of 12 160 cases and 16 838 cases of European ancestry. As a result, we identified three independent SNPs in DCAF4 (rs117781739, rs12587742 and rs2240980) associated with lung cancer risk (odds ratio = 0.91, 1.09 and 1.09, respectively; 95% confidence interval = 0.88-0.95, 1.05-1.14 and 1.05-1.13, respectively; and P = 3.99 × 10-6, 4.97 × 10-5 and 1.44 × 10-5, respectively) after multiple comparison correction by a false discovery rate <0.05. Since SNP rs12587742 is located within the promoter region and one CpG island of DCAF4, we further performed in silico functional analyses and found that the rs12587742 variant A allele was associated with an increased mRNA expression (P = 2.20 × 10-16, 1.79 × 10-13 and 0.001 in blood cells, normal lung tissues and tumor tissues of lung squamous carcinoma, respectively) and a decreased methylation status (P = 2.48 × 10-9 and 0.032 in adipose and lung tumor tissues, respectively). Moreover, evidence from differential expression analyses further supported an oncogenic effect of DCAF4 on lung cancer, with higher mRNA levels in both lung squamous carcinoma and adenocarcinoma (P = 4.48 × 10-11 and 1.22 × 10-9, respectively) than in adjacent normal tissues. Taken together, our results suggest that rs12587742 is associated with an increased lung cancer risk, possibly by up-regulating mRNA expression and decreasing methylation status of DCAF4.

HPV-associated lung cancers: an international pooled analysis.

Human papillomavirus (HPV) is the etiologic risk factor for cervical cancer. Some studies have suggested an association with a subset of lung tumors, but the etiologic link has not been firmly established. We performed an international pooled analysis of cross-sectional studies (27 datasets, n = 3249 patients) to evaluate HPV DNA prevalence in lung cancer and to investigate viral presence according to clinical and demographic characteristics. HPV16/18 were the most commonly detected, but with substantial variation in viral prevalence between geographic regions. The highest prevalence of HPV16/18 was observed in South and Central America, followed by Asia, North America and Europe (adjusted prevalence rates = 22, 5, 4 and 3%, respectively). Higher HPV16 prevalence was noted in each geographic region compared with HPV18, except in North America. HPV16/18-positive lung cancer was less likely observed among White race (adjusted odds ratio [OR] = 0.33, 95% confidence interval [CI] = 0.12-0.90), whereas no associations were observed with gender, smoking history, age, histology or stage. Comparisons between tumor and normal lung tissue show that HPV was more likely to be present in lung cancer rather than normal lung tissues (OR = 3.86, 95% CI = 2.87-5.19). Among a subset of patients with HPV16-positive tumors, integration was primarily among female patients (93%, 13/14), while the physical status in male cases (N = 14) was inconsistent. Our findings confirm that HPV DNA is present in a small fraction of lung tumors, with large geographic variations. Further comprehensive analysis is needed to assess whether this association reflects a causal relationship.

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.

Mosaic Chromosomal Alterations Are Associated With Increased Lung Cancer Risk: Insight From the INTEGRAL-ILCCO Cohort Analysis.

INTRODUCTION: Mosaic chromosomal alterations (mCAs) detected in white blood cells represent a type of clonal hematopoiesis (CH) that is understudied compared with CH-related somatic mutations. A few recent studies indicated their potential link with nonhematological cancers, especially lung cancer. METHODS: In this study, we investigated the association between mCAs and lung cancer using the high-density genotyping data from the OncoArray study of INTEGRAL-ILCCO, the largest single genetic study of lung cancer with 18,221 lung cancer cases and 14,825 cancer-free controls. RESULTS: We identified a comprehensive list of autosomal mCAs, ChrX mCAs, and mosaic ChrY (mChrY) losses from these samples. Autosomal mCAs were detected in 4.3% of subjects, in addition to ChrX mCAs in 3.6% of females and mChrY losses in 9.6% of males. Multivariable logistic regression analysis indicated that the presence of autosomal mCAs in white blood cells was associated with an increased lung cancer risk after adjusting for key confounding factors, including age, sex, smoking status, and race. This association was mainly driven by a specific type of mCAs: copy-neutral loss of heterozygosity on autosomal chromosomes. The association between autosome copy-neutral loss of heterozygosity and increased risk of lung cancer was further confirmed in two major histologic subtypes, lung adenocarcinoma and squamous cell carcinoma. In addition, we observed a significant increase of ChrX mCAs and mChrY losses in smokers compared with nonsmokers and racial differences in certain types of mCA events. CONCLUSIONS: Our study established a link between mCAs in white blood cells and increased risk of lung cancer.

Association of a novel functional promoter variant (rs2075533 C>T) in the apoptosis gene TNFSF8 with risk of lung cancer--a finding from Texas lung cancer genome-wide association study.

Published genome-wide association studies (GWASs) have identified few variants in the known biological pathways involved in lung cancer etiology. To mine the possibly hidden causal single nucleotide polymorphisms (SNPs), we explored all SNPs in the extrinsic apoptosis pathway from our published GWAS dataset for 1154 lung cancer cases and 1137 cancer-free controls. In an initial association analysis of 611 tagSNPs in 41 apoptosis-related genes, we identified only 10 tagSNPs associated with lung cancer risk with a P value<10(-2), including four tagSNPs in DAPK1 and three tagSNPs in TNFSF8. Unlike DAPK1 SNPs, TNFSF8 rs2181033 tagged other four predicted functional but untyped SNPs (rs776576, rs776577, rs31813148 and rs2075533) in the promoter region. Therefore, we further tested binding affinity of these four SNPs by performing the electrophoretic mobility shift assay. We found that only rs2075533T allele modified levels of nuclear proteins bound to DNA, leading to significantly decreased expression of luciferase reporter constructs by 5- to -10-fold in H1299, HeLa and HCT116 cell lines compared with the C allele. We also performed a replication study of the untyped rs2075533 in an independent Texas population but did not confirm the protective effect. We further performed a mini meta-analysis for SNPs of TNFSF8 obtained from other four published lung cancer GWASs with 12 214 cases and 47 721 controls, and we found that only rs3181366 (r2=0.69 with the untyped rs2075533) was associated to lung cancer risk (P=0.008). Our findings suggest a possible role of novel TNFSF8 variants in susceptibility to lung cancer.

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.

Expanding the genetic architecture of nicotine dependence and its shared genetics with multiple traits.

Cigarette smoking is the leading cause of preventable morbidity and mortality. Genetic variation contributes to initiation, regular smoking, nicotine dependence, and cessation. We present a Fagerström Test for Nicotine Dependence (FTND)-based genome-wide association study in 58,000 European or African ancestry smokers. We observe five genome-wide significant loci, including previously unreported loci MAGI2/GNAI1 (rs2714700) and TENM2 (rs1862416), and extend loci reported for other smoking traits to nicotine dependence. Using the heaviness of smoking index from UK Biobank (N = 33,791), rs2714700 is consistently associated; rs1862416 is not associated, likely reflecting nicotine dependence features not captured by the heaviness of smoking index. Both variants influence nearby gene expression (rs2714700/MAGI2-AS3 in hippocampus; rs1862416/TENM2 in lung), and expression of genes spanning nicotine dependence-associated variants is enriched in cerebellum. Nicotine dependence (SNP-based heritability = 8.6%) is genetically correlated with 18 other smoking traits (rg = 0.40-1.09) and co-morbidities. Our results highlight nicotine dependence-specific loci, emphasizing the FTND as a composite phenotype that expands genetic knowledge of smoking.

Frequency of Pathogenic Germline Variants in Cancer-Susceptibility Genes in Patients With Osteosarcoma.

Importance: Osteosarcoma, the most common malignant bone tumor in children and adolescents, occurs in a high number of cancer predisposition syndromes that are defined by highly penetrant germline mutations. The germline genetic susceptibility to osteosarcoma outside of familial cancer syndromes remains unclear. Objective: To investigate the germline genetic architecture of 1244 patients with osteosarcoma. Design, Setting, and Participants: Whole-exome sequencing (n = 1104) or targeted sequencing (n = 140) of the DNA of 1244 patients with osteosarcoma from 10 participating international centers or studies was conducted from April 21, 2014, to September 1, 2017. The results were compared with the DNA of 1062 individuals without cancer assembled internally from 4 participating studies who underwent comparable whole-exome sequencing and 27 173 individuals of non-Finnish European ancestry who were identified through the Exome Aggregation Consortium (ExAC) database. In the analysis, 238 high-interest cancer-susceptibility genes were assessed followed by testing of the mutational burden across 736 additional candidate genes. Principal component analyses were used to identify 732 European patients with osteosarcoma and 994 European individuals without cancer, with outliers removed for patient-control group comparisons. Patients were subsequently compared with individuals in the ExAC group. All data were analyzed from June 1, 2017, to July 1, 2019. Main Outcomes and Measures: The frequency of rare pathogenic or likely pathogenic genetic variants. Results: Among 1244 patients with osteosarcoma (mean [SD] age at diagnosis, 16 [8.9] years [range, 2-80 years]; 684 patients [55.0%] were male), an analysis restricted to individuals with European ancestry indicated a significantly higher pathogenic or likely pathogenic variant burden in 238 high-interest cancer-susceptibility genes among patients with osteosarcoma compared with the control group (732 vs 994, respectively; P = 1.3 × 10-18). A pathogenic or likely pathogenic cancer-susceptibility gene variant was identified in 281 of 1004 patients with osteosarcoma (28.0%), of which nearly three-quarters had a variant that mapped to an autosomal-dominant gene or a known osteosarcoma-associated cancer predisposition syndrome gene. The frequency of a pathogenic or likely pathogenic cancer-susceptibility gene variant was 128 of 1062 individuals (12.1%) in the control group and 2527 of 27 173 individuals (9.3%) in the ExAC group. A higher than expected frequency of pathogenic or likely pathogenic variants was observed in genes not previously linked to osteosarcoma (eg, CDKN2A, MEN1, VHL, POT1, APC, MSH2, and ATRX) and in the Li-Fraumeni syndrome-associated gene, TP53. Conclusions and Relevance: In this study, approximately one-fourth of patients with osteosarcoma unselected for family history had a highly penetrant germline mutation requiring additional follow-up analysis and possible genetic counseling with cascade testing.

On combining family and case-control studies.

Studies to detect genetic association with disease can be family-based, often using families with multiple affected members, or population based, as in population-based case-control studies. If data on both study types are available from the same population, it is useful to combine them to improve power to detect genetic associations. Two aspects of the data need to be accommodated, the sampling scheme and potential residual correlations among family members. We propose two approaches for combining data from a case-control study and a family study that collected families with multiple cases. In the first approach, we view a family as the sampling unit and specify the joint likelihood for the family members using a two-level mixed effects model to account for random familial effects and for residual genetic correlations among family members. The ascertainment of the families is accommodated by conditioning on the ascertainment event. The individuals in the case-control study are treated as families of size one, and their unconditional likelihood is combined with the conditional likelihood for the families. This approach yields subject specific maximum likelihood estimates of covariate effects. In the second approach, we view an individual as the sampling unit. The sampling scheme is accommodated using two-phase sampling techniques, marginal covariate effects are estimated, and correlations among family members are accounted for in the variance calculations. The models are compared in simulations. Data from a case-control and a family study from north-eastern Italy on melanoma and a low-risk melanoma-susceptibility gene, MC1R, are used to illustrate the approaches.

Alcohol consumption and lung cancer risk: A pooled analysis from the International Lung Cancer Consortium and the SYNERGY study.

BACKGROUND: There is inadequate evidence to determine whether there is an effect of alcohol consumption on lung cancer risk. We conducted a pooled analysis of data from the International Lung Cancer Consortium and the SYNERGY study to investigate this possible association by type of beverage with adjustment for other potential confounders. METHODS: Twenty one case-control studies and one cohort study with alcohol-intake data obtained from questionnaires were included in this pooled analysis (19,149 cases and 362,340 controls). Adjusted odds ratios (OR) or hazard ratios (HR) with corresponding 95% confidence intervals (CI) were estimated for each measure of alcohol consumption. Effect estimates were combined using random or fixed-effects models where appropriate. Associations were examined for overall lung cancer and by histological type. RESULTS: We observed an inverse association between overall risk of lung cancer and consumption of alcoholic beverages compared to non-drinkers, but the association was not monotonic. The lowest risk was observed for persons who consumed 10-19.9 g/day ethanol (OR vs. non-drinkers = 0.78; 95% CI: 0.67, 0.91), where 1 drink is approximately 12-15 g. This J-shaped association was most prominent for squamous cell carcinoma (SCC). The association with all lung cancer varied little by type of alcoholic beverage, but there were notable differences for SCC. We observed an association with beer intake (OR for ≥20 g/day vs nondrinker = 1.42; 95% CI: 1.06, 1.90). CONCLUSIONS: Whether the non-monotonic associations we observed or the positive association between beer drinking and squamous cell carcinoma reflect real effects await future analyses and insights about possible biological mechanisms.

VTET: a variable threshold exact test for identifying disease-associated copy number variations enriched in short genomic regions.

Copy number variations (CNVs) constitute a major source of genetic variations in human populations and have been reported to be associated with complex diseases. Methods have been developed for detecting CNVs and testing CNV associations in genome-wide association studies (GWAS) based on SNP arrays. Commonly used two-step testing procedures work well only for long CNVs while direct CNV association testing methods work only for recurrent CNVs. Assuming that short CNVs disrupting any part of a given genomic region increase disease risk, we developed a variable threshold exact test (VTET) for testing disease associations of CNVs randomly distributed in the genome using intensity data from SNP arrays. By extensive simulations, we found that VTET outperformed two-step testing procedures based on existing CNV calling algorithms for short CNVs and that the performance of VTET was robust to the length of the genomic region. In addition, VTET had a comparable performance with CNVtools for testing the association of recurrent CNVs. Thus, we expect VTET to be useful for testing disease associations of both recurrent and randomly distributed CNVs using existing GWAS data. We applied VTET to a lung cancer GWAS and identified a genome-wide significant region on chromosome 18q22.3 for lung squamous cell carcinoma.