Hsiung, Chao A.

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.

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.

Abstract To identify the genetic factors that influence overall survival in never smokers who have non–small cell lung carcinoma (NSCLC), we conducted a consistency meta-analysis study using genome-wide association approaches for overall survival in 327 never smoker patients with NSCLC from The University of Texas MD Anderson Cancer Center (Houston, TX) and 293 cases from the Mayo Clinic (Rochester, MN). We then conducted a two-pronged validation of the top 25 variants that included additional validation in 1,256 patients with NSCLC from Taiwan and assessment of expression quantitative trait loci (eQTL) and differential expression of genes surrounding the top loci in 70 tumors and matched normal tissues. A total of 94 loci were significant for overall survival in both MD Anderson and Mayo studies in the consistency meta-analysis phase, with the top 25 variants reaching a P value of 10−6. Two variants of these 25 were also significant in the Taiwanese population: rs6901416 [HR, 1.44; 95% confidence interval (CI), 1.01–2.06] and rs10766739 (HR, 1.23; 95%CI, 1.00–1.51). These loci resulted in a reduction of median survival time of at least eight and five months in three populations, respectively. An additional six variants (rs4237904, rs7976914, rs4970833, rs954785, rs485411, and rs10906104) were validated through eQTL analysis that identified significant correlations with expression levels of six genes (LEMD3, TMBIM, ATXN7L2, SHE, ITIH2, and NUDT5, respectively) in normal lung tissue. These genes were also significantly differentially expressed between the tumor and normal lung tissue. These findings identify several novel, candidate prognostic markers for NSCLC in never smokers, with eQTL analysis suggesting a potential biologic mechanism for a subset of these observed associations. Cancer Res; 73(13); 4028–38. ©2013 AACR.

Genome-wide variance components linkage analysis was performed on 4 latent factors underlying metabolic syndrome derived from 10 risk factors. The latent factors represent obesity and insulin, blood pressure, lipids and insulin, and central obesity. The metabolic syndrome factor scores were derived in 4 ethnic groups recruited in 3 Networks of the Family Blood Pressure Program: GENOA (blacks, Hispanics, and whites), HyperGEN (blacks and whites), SAPPHIRe (Asians). Heritabilities of metabolic syndrome factors ranged from 66% for obesity and insulin to 11% for blood pressure factor. We observed higher heritabilities for obesity and insulin, and lipids and insulin, whereas those for blood pressure and central obesity were smaller. Linkage analysis detected two major quantitative trait loci. One of them linked to the obesity and insulin factor with a lod score of 3.94 ( P =0.00001, marker GATA11A06, D18S53, 41.24 cM) at marker positions linkage (lod 4.71, at 46.84 cM at 1-cM-apart distances linkage), located on chromosome 18p11.21 in GENOA black. The other linked to the blood pressure factor with a lod score of 3.22 ( P =0.000059, marker GATA49C09, D17S1290, 82 cM) at marker positions linkage (lod 3.56, at 84.63 cM for 1 cM apart distances linkage) located on chromosome 17q23.1 in Hispanics. These quantitative trait loci, together with 4 additional ones with lod scores >2.5, and 30 additional ones with lod score >1.7, offer hope for dissecting the genetic architecture of metabolic syndrome with beneficial implications for molecular diagnosis, prognosis, and in potential medical intervention.