Liu, Geoffrey

Abstract Background: Over 20 susceptibility single-nucleotide polymorphisms (SNP) have been identified for esophageal adenocarcinoma (EAC) and its precursor, Barrett esophagus (BE), explaining a small portion of heritability. Methods: Using genetic data from 4,323 BE and 4,116 EAC patients aggregated by international consortia including the Barrett's and Esophageal Adenocarcinoma Consortium (BEACON), we conducted a comprehensive transcriptome-wide association study (TWAS) for BE/EAC, leveraging Genotype Tissue Expression (GTEx) gene-expression data from six tissue types of plausible relevance to EAC etiology: mucosa and muscularis from the esophagus, gastroesophageal (GE) junction, stomach, whole blood, and visceral adipose. Two analytical approaches were taken: standard TWAS using the predicted gene expression from local expression quantitative trait loci (eQTL), and set-based SKAT association using selected eQTLs that predict the gene expression. Results: Although the standard approach did not identify significant signals, the eQTL set–based approach identified eight novel associations, three of which were validated in independent external data (eQTL SNP sets for EXOC3, ZNF641, and HSP90AA1). Conclusions: This study identified novel genetic susceptibility loci for EAC and BE using an eQTL set–based genetic association approach. Impact: This study expanded the pool of genetic susceptibility loci for EAC and BE, suggesting the potential of the eQTL set–based genetic association approach as an alternative method for TWAS analysis.

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

Abstract Background: Somatic EGFR mutations define a subset of non–small cell lung cancers (NSCLC) that have clinical impact on NSCLC risk and outcome. However, EGFR-mutation-status is often missing in epidemiologic datasets. We developed and tested pragmatic approaches to account for EGFR-mutation-status based on variables commonly included in epidemiologic datasets and evaluated the clinical utility of these approaches. Methods: Through analysis of the International Lung Cancer Consortium (ILCCO) epidemiologic datasets, we developed a regression model for EGFR-status; we then applied a clinical-restriction approach using the optimal cut-point, and a second epidemiologic, multiple imputation approach to ILCCO survival analyses that did and did not account for EGFR-status. Results: Of 35,356 ILCCO patients with NSCLC, EGFR-mutation-status was available in 4,231 patients. A model regressing known EGFR-mutation-status on clinical and demographic variables achieved a concordance index of 0.75 (95% CI, 0.74–0.77) in the training and 0.77 (95% CI, 0.74–0.79) in the testing dataset. At an optimal cut-point of probability-score = 0.335, sensitivity = 69% and specificity = 72.5% for determining EGFR-wildtype status. In both restriction-based and imputation-based regression analyses of the individual roles of BMI on overall survival of patients with NSCLC, similar results were observed between overall and EGFR-mutation-negative cohort analyses of patients of all ancestries. However, our approach identified some differences: EGFR-mutated Asian patients did not incur a survival benefit from being obese, as observed in EGFR-wildtype Asian patients. Conclusions: We introduce a pragmatic method to evaluate the potential impact of EGFR-status on epidemiological analyses of NSCLC. Impact: The proposed method is generalizable in the common occurrence in which EGFR-status data are missing.

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

Abstract Genome-wide association studies (GWAS) of esophageal adenocarcinoma (EAC) and its precursor, Barrett’s esophagus (BE), have uncovered significant genetic components of risk, but most heritability remains unexplained. Targeted assessment of genetic variation in biologically relevant pathways using novel analytical approaches may identify missed susceptibility signals. Central obesity, a key BE/EAC risk factor, is linked to systemic inflammation, altered hormonal signaling and insulin-like growth factor (IGF) axis dysfunction. Here, we assessed IGF-related genetic variation and risk of BE and EAC. Principal component analysis was employed to evaluate pathway-level and gene-level associations with BE/EAC, using genotypes for 270 single-nucleotide polymorphisms (SNPs) in or near 12 IGF-related genes, ascertained from 3295 BE cases, 2515 EAC cases and 3207 controls in the Barrett’s and Esophageal Adenocarcinoma Consortium (BEACON) GWAS. Gene-level signals were assessed using Multi-marker Analysis of GenoMic Annotation (MAGMA) and SNP summary statistics from BEACON and an expanded GWAS meta-analysis (6167 BE cases, 4112 EAC cases, 17 159 controls). Global variation in the IGF pathway was associated with risk of BE (P = 0.0015). Gene-level associations with BE were observed for GHR (growth hormone receptor; P = 0.00046, false discovery rate q = 0.0056) and IGF1R (IGF1 receptor; P = 0.0090, q = 0.0542). These gene-level signals remained significant at q &lt; 0.1 when assessed using data from the largest available BE/EAC GWAS meta-analysis. No significant associations were observed for EAC. This study represents the most comprehensive evaluation to date of inherited genetic variation in the IGF pathway and BE/EAC risk, providing novel evidence that variation in two genes encoding cell-surface receptors, GHR and IGF1R, may influence risk of BE.

AbstractAlthough the development of effective combined chemoradiation regimens for esophageal cancers has resulted in statistically significant survival benefits, the majority of patients treated with curative intent develop locoregional and/or distant relapse. Further improvements in disease control and survival will require the development of individualized therapy based on the knowledge of host and tumor genomics and potentially harnessing the host immune system. Although there are a number of gene targets that are amplified and proteins that are overexpressed in esophageal cancers, attempts to target several of these have not proven successful in unselected patients. Herein, we review our current state of knowledge regarding the molecular pathways implicated in esophageal carcinoma, and the available agents for targeting these pathways that may rationally be combined with standard chemoradiation, with the hope that this commentary will guide future efforts of novel combinations of therapy.

Abstract Background: Lung cancer remains the leading cause of cancer death worldwide, with 15% to 20% occurring in never smokers. To assess genetic determinants for prognosis among never smokers, we conducted a genome-wide investigation in the International Lung Cancer Consortium (ILCCO). Methods: Genomic and clinical data from 1,569 never-smoking patients with lung cancer of European ancestry from 10 ILCCO studies were included. HRs and 95% confidence intervals of overall survival were estimated. We assessed whether the associations were mediated through mRNA expression–based 1,553 normal lung tissues from the lung expression quantitative trait loci (eQTL) dataset and Genotype-Tissue Expression (GTEx). For cross-ethnicity generalization, we assessed the associations in a Japanese study (N = 887). Results: One locus at 13q22.2 was associated with lung adenocarcinoma survival at genome-wide level, with carriers of rs12875562-T allele exhibiting poor prognosis [HR = 1.71 (1.41–2.07), P = 3.60 × 10−8], and altered mRNA expression of LMO7DN in lung tissue (GTEx, P = 9.40 × 10−7; Lung eQTL dataset, P = 0.003). Furthermore, 2 of 11 independent loci that reached the suggestive significance level (P &lt; 10−6) were significant eQTL affecting mRNA expression of nearby genes in lung tissues, including CAPZB at 1p36.13 and UBAC1 at 9q34.3. One locus encoding NWD2/KIAA1239 at 4p14 showed associations in both European [HR = 0.50 (0.38–0.66), P = 6.92 × 10−7] and Japanese populations [HR = 0.79 (0.67–0.94), P = 0.007]. Conclusions: Based on the largest genomic investigation on the lung cancer prognosis of never smokers to date, we observed that lung cancer prognosis is affected by inherited genetic variants. Impact: We identified one locus near LMO7DN at genome-wide level and several potential prognostic genes with cis-effect on mRNA expression. Further functional genomics work is required to understand their role in tumor progression.

1512 Background: Tobacco smoking profoundly impacts LC risk; however, data are limited as to what extent SC prior to diagnosis impacts LC overall survival (OS) and lung cancer specific survival (LCSS). LC screening offers a possible teachable moment, but there is uncertainty of SC benefits after a lifetime of smoking. We use the ILCCO database to answer if SC prior to LC dx is associated with better OS and LCSS, considering time since smoking cessation (TSSC). Methods: Using individual data, analysis was performed on 17 ILCCO studies with available TSSC to estimate survival using univariable analysis and models of stage-adjusted and cumulative smoking-adjusted multivariable analysis. Adjusted Hazard Ratios (aHR) from Cox models, cubic spline smooth curves and Kaplan-Meier curves were created. Sensitivity analysis was performed for TSSC and LCSS on 13 studies. Results: Of 34649 patients, 14322 (41%) were current smokers 14273 (41%) ex-smokers and 6054 (18%) never smokers at diagnosis. We confirmed that ex-smokers (aHR 0.88 CI 0.86-0.91) and never smokers (aHR 0.76 CI 0.73-0.8) improved OS compared to current smokers. Amongst ex-smokers, < 2y TSSC (aHR 0.88 CI 0.82-0.94), 2-5y TSSC (aHR 0.83 CI 0.77-0.90) and > 5y TSSC (aHR 0.8 CI 0.76-0.84) had improved OS compared to CS. Sensitivity analysis showed a trend towards improved LCSS survival for < 2y TSSC (aHR 0.95 CI 0.86-1.05) and 2-5y TSSC (aHR 0.93 CI 0.83-1.04), whereas > 5y TSSC significantly improved LCSS by 15% (aHR 0.85 CI 0.78-0.92). To mimic the LC screening participants, in analysis of > 30 pack-years individuals, associations were strikingly strong: < 2y TSSC had improved OS by 14% (aHR 0.86 CI 0.80-0.93); 2-5y TSSC by 17% (aHR 0.83 CI 0.76-0.90); and > 5 TSSC by 22% (aHR 0.78 CI 0.74-0.83), compared to current smokers; for < 30 packs-years, a trend towards better OS was observed for < 2y TSSC (aHR 0.95 CI 0.92-1.02) and 2-5y TSSC (aHR 0.86 CI 0.74-1.01), whereas > 5y TSSC improved OS by 23% (aHR 0.77 CI 0.72-0.82). Conclusions: Among ex-smokers, the risk of overall death was reduced by 12% on < 2y TSSC, 17% on 2-5y TSSC and 20% > 5y TSSC, whereas for LCSS, the benefit was significant only for > 5y TSCC, compared to current smokers at time of diagnosis. Here we demonstrate that convincing screening participants to quit smoking at any point of their trajectory, even just prior to dx such as < 2y TSSC, improved OS, and LCSS benefit was present beyond 5y of quitting. These relationships are independent of pack-years, age, across all stages and other prognostic variables.

571 Background: Tumor-induced generation of splenic erythroblast- like cells (Ter-cells) has been shown to promote tumor progression. These Ter-cells produce the glial-line derived neurotropic factor, ARTN. We investigated the association of pre-operative serum ARTN and the risk of recurrence in HCC patients (pts) undergoing curative resection. Methods: Blood samples were collected prior to surgery as part of an institutional molecular epidemiologic study. Serum ARTN concentration was measured using an enzyme-linked immunosorbent assay (ELISA). Demographics, pathological variables known to be associated with outcomes and clinical outcomes were collected. Uni- and multi-variate analysis were conducted. The optimal cutpoint method was used to define high and low ARTN concentrations. Cox models (hazard ratios, HR) were used to compare progression-free (PFS, primary endpoint), and overall survival (OS) of pts with high vs low serum ARTN. Results: Pre-operative blood samples were available for 58 pts. Median age was 63 years (range: 25-85 years); 74% were male and 50% were Asian. Etiology of liver disease was hepatitis B (43%) and hepatitis C (22%); 43% of tumors were ≤5cm, and vascular invasion was seen in 62%. A baseline alpha-fetoprotein (AFP) of > 100 mcg/L was observed in 36% pts. Median follow-up was 18.9 months. Median ARTN concentration was 0.322 ng/mL. The optimal ARTN concentration cut-off was 0.206 ng/mL. Median PFS for pts with high ( > 0.206 ng/mL) vs low (≤0.206 ng/mL) ARTN was 15.7 vs 55 months (p = 0.04) respectively. Three year PFS was 34% vs 55%, and three year OS 54% vs 91% for high vs low groups. Univariate analysis found that high ARTN concentration (HR 2.44, p = 0.05) and multifocal tumors were associated with a worse PFS. In a multivariate analysis adjusted for AFP > 100mcg/L, vascular invasion, hepatitis status and multifocal tumors, high ARTN remained a negative prognostic factor for PFS, aHR 3.32 (95% CI: 1.22-9.05, p = 0.02). Conclusions: A high pre-surgery serum ARTN is associated with earlier recurrence in HCC pts undergoing curative resection. ARTN should be further studied in HCC to determine its value as a prognostic marker.

1562 Background: In Caucasian-predominant populations, overweight or obese NSCLC pts (BMI≥25-kg/m2) have better prognosis while underweight (BMI≤18.5-kg/m2) pts have a worse prognosis. A large pooled sample allowed us to evaluate the role of ethnicity in this BMI-NSCLC OS relationship. Methods: Using individual data, survival analysis was performed on 15 ILCCO studies, assessing the interactions between ethnicity and BMI on overall survival (OS). Adjusted Hazard Ratios (aHR) from Cox models and adjusted penalized smoothing spline plots were generated. Results: Among 13416 (77%) Caucasian, 2975 (17%) Asian, and 935 (5%) Black NSCLC pts analyzed, we confirmed that for all pts, being underweight at NSCLC diagnosis was associated with worse OS (aHR=1.68, CI 1.5-1.8, p<0.001), while overweight/obese pts had improved survival (aHR 0.89, CI 0.8-0.9, p<0.001), when compared to pts with normal BMI. In general, Black pts had poorer OS than Caucasian pts (aHR 1.26 CI 1.1-1.4, p<0.001). However, the BMI-OS relationship differed according to ethnicity (BMI-ethnicity interaction, p=0.009): Caucasian underweight pts had poorer OS (aHR 1.67, CI 1.5-1.8, P<0.001) while overweight/obese pts had improved OS (aHR 0.89, CI 0.8-0.9, P<0.001). In Asian pts, these two associations were aHR 1.15, CI 0.8-1.5, P=0.33, and aHR 0.91, CI 0.7-1.1, P=0.32, respectively. In Black pts, these two associations were aHR 1.06, CI 0.7-1.5, p=0.73 and 0.75, CI 0.6-0.8 p<0.001, respectively. For overweight/obese patients, as BMI rises from 25-kg/m2 through 45-kg/m2, the prognosis worsens for Caucasian pts, remains stable for Asian pts, but improves for Black pts. Conclusions: InCaucasian pts, being underweight had a greater negative impact on OS than for Asian or Black pts, while being obese had a greater beneficial impact in Blacks than in other ethnic groups. These ethnic differences likely reflect genetically-informed muscle/adipose tissue distributions, where Black pts may have less sarcopenic obesity than other ethnicities. In future prognostic studies, BMI relationships must account for ethnic differences.