Liu, Geoffrey

Abstract Understanding tumor cell dynamics can improve prognosis and treatment but remains limited for lung adenocarcinoma in people who have never smoked (NS-LUAD). With RNA sequencing data from 684 NS-LUAD cases and validation in an independent dataset, we identified three subtypes with distinct phenotypic traits and cell compositions. Additional genomic and histologic data further characterized the subtypes. “Steady,” marked by low proliferation, high alveolar cell fraction, moderate-to-well differentiation, and fewer driver gene alterations, is linked to prolonged survival and low immune evasion. “Proliferative” shows high proliferation markers, TP53 mutations, and gene fusions. “Chaotic,” with high epithelial-to-mesenchymal transition markers, has the worst prognosis, even within stage I tumors. Lacking known molecular or histologic characteristics, this aggressive subtype is solely identified by transcriptomic data. A 60-gene signature recapitulates the classification and predicts survival even within subgroups based on tumor stage or known genomic features, emphasizing its potential for improving early-stage NS-LUAD prognostication in clinical settings. Significance: The transcriptome of 684 NS-LUAD identifies three subtypes with different cellular dynamics and genomic and morphologic features. A 60-gene signature accurately stratifies subjects for mortality risk, even in stage I, offering a potential clinically applicable tool for treatment decision-making in patients with NS-LUAD. See related commentary by Azizi et al., p. 423

Abstract Background Lung adenocarcinoma in never-smokers is a major public health burden, especially among East Asian women. Polygenic risk scores (PRSs) are promising for risk stratification but are primarily developed in European-ancestry populations. We aimed to develop and validate single- and multi-ancestry PRSs for East Asian never-smokers to improve lung adenocarcinoma risk prediction. Methods PRSs were developed using genome-wide association study summary statistics from East Asian (8002 cases; 20 782 controls) and European (2058 cases; 5575 controls) populations. Single-ancestry models included a PRS with 25 single nucleotide variants (PRS-25), PRS clumping and thresholding (CT), and LDpred2; multi-ancestry models included LDpred2 plus European PRS-128, PRS-CSx, and clumping and thresholding, super learning, empirical Bayes (CT-SLEB). Performance was evaluated in independent East Asian data from the Female Lung Cancer Consortium and externally validated in the Nanjing Lung Cancer Cohort. We assessed predictive accuracy via area under the curve operating procedure (AUC), with 10-year and (ages 30-80 years) absolute risks estimates. Results The best multi-ancestry PRS, using East Asian and European data via CT-SLEB, outperformed the best East Asian-only (LDpred2; AUC = 0.629, 95% confidence interval [CI] = 0.618 to 0.641), achieving an AUC of 0.640 (95% CI = 0.629 to 0.653) and odds ratio of 1.71 (95% CI = 1.61 to 1.82) per standard deviation increase. Nanjing Lung Cancer Cohort validation confirmed robust performance (AUC = 0.649, 95% CI = 0.623 to 0.676). The top 20% PRS group had a 3.92-fold higher lung adenocarcinoma risk than the bottom 20%. Further, the top 5% PRS group reached a 6.69% lifetime absolute risk. Notably, this group reached the average population 10-year lung adenocarcinoma risk at age 50 years (0.42%) by age 41 years, 9 years earlier. Conclusions Multi-ancestry PRS approaches enhance lung adenocarcinoma risk stratification in East Asian never-smokers, with consistent external validation, suggesting future clinical utility.

Background: Rare, functionally deleterious genetic variants have been implicated as having substantially larger effect sizes than common variants, potentially accounting for much of the missing heritability in lung cancer (LC). Over 25 years, the Genetic Epidemiology of LC Consortium (GELCC) has curated an invaluable resource of specimens and data from individuals with a strong family history of LC. Methods: Whole-genome and exome sequencing data were analyzed in 129 high-risk familial LC (FLC) families, defined as having two first-degree relatives or three or more second-degree relatives affected by LC (177 affected FLC and 309 unaffected relatives). The discovery focused on rare (allele frequency < 1% in gnomAD global population) and deleterious variations (missense, stop-gained, frameshift) enriched in FLC. Promising candidates were further external validated in 2,408 sporadic lung cancer (SLC) cases and 885 population controls from the International Lung Cancer Consortium. Results: We identified 168 variants with moderate-to-large effects on FLC risk, 100 of which were validated for SLC risk (Table 1 Table 1. Top 20 candidates Gene Rare (MAF < 1%), deleterious variants No. of variant carrier counts Odds Ratio, LC vs Control FLC | Relative | SLC | Control FLC+SLC SLC Oncogene ERBB3 / HER3 p.A1131T 3 | 1 | 2 | 0 35 15 JAK p.V651M 2 | 3 | 7 | 0 62 51 NOTCH3 p.R1560P 3 | 1 | 4 | 0 48 29 ROS1 p.N785X 6 | 4 | 10 | 0 110 74 USP6 p.R522 fs del 3 | 6 | 6 | 0 62 45 Tumor suppressor gene ATM p.L2332P 3 | 2 | 8 | 0 76 59 LZTS2 p.K458R 4 | 4 | 18 | 0 152 133 NPIPB13 p.V270 fs del 13 | 14 | 10 | 0 159 74 PARK2 / PRKN p.P153R 3 | 3 | 8 | 0 76 59 RB1 p.A525G 3 | 0 | 6 | 0 62 45 SLX4 p.S1716T 3 | 2 | 9 | 0 83 66 SYNE1 p.M6566I 2 | 0 | 7 | 0 62 52 TGFBRAP1 p.T301R 5 | 0 | 11 | 0 110 82 WNK1 p.K583 fs del 4 | 0 | 23 | 0 187 171 Other cancer-associated gene CFTR p.R74W 4 | 2 | 4 | 0 55 30 FLG p.R3404 fs del 7 | 2 | 4 | 0 69 30 LTN1 p.R349H 6 | 1 | 10 | 0 110 74 MLNR p.Q334 fs del 3 | 2 | 19 | 0 152 134 SLC39A11 p.R38Q 4 | 2 | 13 | 0 117 96 TTYH2 p.R493C 5 | 4 | 14 | 0 131 104 Dose effect of the top candidates No. candidates 0 variant allele 67 | 194 | 1911 | 834 reference reference 1 variant allele 50 | 66 | 363 | 50 3.5 (2.6 - 4.7) 3.2 (2.3 - 4.3) 2 variant alleles 15 | 18 | 62 | 1 32 (4.5 - 230) 27 (3.8 - 195) 3+ variant alleles 45 | 31 | 72 | 0 98 (6 - 1506) 62 (3.9 -1012) lists 20 top hits). Among these, 15 variants mapped to our previously identified lung cancer linkage locus at 6q23-25 (ROS1, PARK2, SYNE1). Others were enriched in oncogenes (ERBB3, JAK1, MET), tumor suppressor genes (ATM, BRCA2, RB1), and genes involved in the extracellular matrix (COL6A3, FLG, MUC5B). Importantly, individuals carrying more than two combinations of these variants exhibited exceptionally strong dose effects, with odds ratios (OR) exceeding 27. Conclusion: Our findings underscore the significant role of rare, high-penetrance variants in FLC etiology and highlight promising targets for early detection and personalized treatment. Future in-depth mechanistic studies are planned to evaluate the pathogenic effects of these specific alleles. Citation Format: Yanhong Liu, Claudio Pikielny, Xiangjun Xiao, Bo Peng, Yafang Li, Jinyoung Byun, Chao Cheng, Dakai Zhu, Spiridon Tsavachidis, Colette Gaba, Elena Kupert, Ellen L. Goode, Erin L. Crawford, Kristen Purrington, Marshall Anderson, Michael Cole, Paul Brennan, Geoffrey Liu, James McKay, John K. Field, Rayjean J. Hung, David C. Christiani, Diptasri Mandal, James C. Willey, Ann Schwartz, Joan Bailey-Wilson, Susan M. Pinney, Christopher I. Amos. High-penetrance rare variants underlying familial lung cancer risk: Insights from the Genetic Epidemiology of Lung Cancer Consortium [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 3815.

ABSTRACT Understanding lung cancer evolution can identify tools for intercepting its growth. In a landscape analysis of 1024 lung adenocarcinomas (LUAD) with deep whole-genome sequencing integrated with multiomic data, we identified 542 LUAD that displayed diverse clonal architecture. In this group, we observed an interplay between mobile elements, endogenous and exogenous mutational processes, distinct driver genes, and epidemiological features. Our results revealed divergent evolutionary trajectories based on tobacco smoking exposure, ancestry, and sex. LUAD from smokers showed an abundance of tobacco-related C:G>A:T driver mutations in KRAS plus short subclonal diversification. LUAD in never smokers showed early occurrence of copy number alterations and EGFR mutations associated with SBS5 and SBS40a mutational signatures. Tumors harboring EGFR mutations exhibited long latency, particularly in females of European-ancestry (EU_N). In EU_N, EGFR mutations preceded the occurrence of other driver genes, including TP53 and RBM10. Tumors from Asian never smokers showed a short clonal evolution and presented with heterogeneous repetitive patterns for the inferred mutational order. Importantly, we found that the mutational signature ID2 is a marker of a previously unrecognized mechanism for LUAD evolution. Tumors with ID2 showed short latency and high L1 retrotransposon activity linked to L1 promoter demethylation. These tumors exhibited an aggressive phenotype, characterized by increased genomic instability, elevated hypoxia scores, low burden of neoantigens, propensity to develop metastasis, and poor overall survival. Re-activated L1 retrotransposition-induced mutagenesis can contribute to the origin of the mutational signature ID2, including through the regulation of the transcriptional factor ZNF695 , a member of the KZFP family. The complex nature of LUAD evolution creates both challenges and opportunities for screening and treatment plans.

Accurate assessment of GFR is crucial to guiding drug eligibility, dosing of systemic therapy, and minimizing the risks of both undertreatment and toxicity in patients with cancer. Up to 32% of patients with cancer have baseline CKD, and both malignancy and treatment may cause kidney injury and subsequent CKD. To date, there has been lack of guidance to standardize approaches to GFR estimation in the cancer population. In this two-part statement from the American Society of Onco-Nephrology, we present key messages for estimation of GFR in patients with cancer, including the choice of GFR estimating equation, use of race and body surface area adjustment, and anticancer drug dose-adjustment in the setting of CKD. These key messages are based on a systematic review of studies assessing GFR estimating equations using serum creatinine and cystatin C in patients with cancer, against a measured GFR comparator. The preponderance of current data involving validated GFR estimating equations involves the CKD Epidemiology Collaboration (CKD-EPI) equations, with 2508 patients in whom CKD-EPI using serum creatinine and cystatin C was assessed (eight studies) and 15,349 in whom CKD-EPI with serum creatinine was assessed (22 studies). The former may have improved performance metrics and be less susceptible to shortfalls of eGFR using serum creatinine alone. Since included studies were moderate quality or lower, the American Society of Onco-Nephrology Position Committee rated the certainty of evidence as low. Additional studies are needed to assess the accuracy of other validated eGFR equations in patients with cancer. Given the importance of accurate and timely eGFR assessment, we advocate for the use of validated GFR estimating equations incorporating both serum creatinine and cystatin C in patients with cancer. Measurement of GFR via exogenous filtration markers should be considered in patients with cancer for whom eGFR results in borderline eligibility for therapies or clinical trials.

e24143 Background: Accurate assessment of glomerular filtration rate (GFR) is crucial to guiding drug eligibility, dose-adjusting systemic therapy, and minimizing the risks of both undertreatment and toxicity in patients with cancer. To date, there has been a lack of guidance to standardize approaches to GFR estimation in the cancer population based on the available evidence. We aim to present data on the first systematic review (SR) of GFR estimating equations in patients with cancer. Methods: We conducted a SR on behalf of the American Society of Onco-Nephrology (ASON) on the accuracy of equations to estimate GFR in adult patients with cancer, in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. We required that studies use a reference methodology as comparator involving measured GFR via urinary or plasma clearance of an exogenous filtration marker or carboplatin. Two search strategies focusing on studies using serum creatinine (Cr) and cystatin C (Cys) were developed. Risk of bias and quality assessment was performed for all studies. Results: A total of 39 studies (21,949 patients in total, 19,025 with solid cancers) were included. Studies assessing populations in Europe were most frequent (50%), followed by Asia (27%) and the remaining from Australia/New Zealand and North and South America. Risk of bias and quality assessment allocated 15 studies (38%) as low or very low quality, and 24 (62%) as low-moderate or moderate quality. The most commonly assessed GFR estimating equations included Cockroft-Gault (32 studies), Modification of Diet Renal Disease (MDRD) (24 studies), Chronic Kidney Disese Epidemiology Collaboration Study (CKD-EPI)Cr (22 studies), Wright (13 studies), CKD-EPICys (8 studies), CKD-EPICr-Cys (8 studies), and Jeliffe (8 studies). Given the importance of accurate and timely eGFR assessment and the evidence supporting the performance (i.e., accuracy) of the CKD-EPI incorporating both Cr and Cys, we suggest the use of this equation in patients with cancer (grade 2C). Where Cys is unavailable, CKD-EPICr is the next preferred equation given the supporting data for its performance and widespread availability. Since included studies were moderate quality or lower, the ASON rated the certainty of evidence as low. Also, additional data on clinical outcomes related to the selection of eGFR equations is needed. We suggest measurement of GFR via an exogenous filtration marker in patients with cancer for whom eGFR results in borderline eligibility for therapies or clinical trials (ungraded). Conclusions: Based on the results of this SR, we suggest use of contemporary methods for estimation of GFR (CKD-EPICr-Cys and CKD-EPICr-Cr) and avoidance of older and less accurate formulae such as Cockcroft-Gault.

Background and Aims: Despite the substantial impact of environmental factors, individuals with a family history of liver cancer have an increased risk for HCC. However, genetic factors have not been studied systematically by genome-wide approaches in large numbers of individuals from European descent populations (EDP). Approach and Results: We conducted a 2-stage genome-wide association study (GWAS) on HCC not affected by HBV infections. A total of 1872 HCC cases and 2907 controls were included in the discovery stage, and 1200 HCC cases and 1832 controls in the validation. We analyzed the discovery and validation samples separately and then conducted a meta-analysis. All analyses were conducted in the presence and absence of HCV. The liability-scale heritability was 24.4% for overall HCC. Five regions with significant ORs (95% CI) were identified for nonviral HCC: 3p22.1, MOBP, rs9842969, (0.51, [0.40–0.65]); 5p15.33, TERT, rs2242652, (0.70, (0.62–0.79]); 19q13.11, TM6SF2, rs58542926, (1.49, [1.29–1.72]); 19p13.11 MAU2, rs58489806, (1.53, (1.33–1.75]); and 22q13.31, PNPLA3, rs738409, (1.66, [1.51–1.83]). One region was identified for HCV-induced HCC: 6p21.31, human leukocyte antigen DQ beta 1, rs9275224, (0.79, [0.74–0.84]). A combination of homozygous variants of PNPLA3 and TERT showing a 6.5-fold higher risk for nonviral-related HCC compared to individuals lacking these genotypes. This observation suggests that gene-gene interactions may identify individuals at elevated risk for developing HCC. Conclusions: Our GWAS highlights novel genetic susceptibility of nonviral HCC among European descent populations from North America with substantial heritability. Selected genetic influences were observed for HCV-positive HCC. Our findings indicate the importance of genetic susceptibility to HCC development.

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