Li, Yafang

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 The roles of sex hormones such as estradiol, testosterone, and sex hormone-binding globulin (SHBG) in the etiology of lung and colorectal cancers in women, among the most common cancers after breast cancer, are unclear. This Mendelian randomization (MR) study evaluated such potential causal associations in women of European ancestry. We used summary statistics data from genome-wide association studies (GWASs) on sex hormones and from the Trøndelag Health (HUNT) Study and large consortia on cancers. There was suggestive evidence of genetically predicted 1-standard deviation increase in total testosterone levels being associated with a lower risk of lung non-adenocarcinoma (hazard ratio (HR) 0.60, 95% CI 0.37–0.98) in the HUNT Study. However, this was not confirmed by using data from a larger consortium. In general, we did not find convincing evidence to support a causal role of sex hormones on risk of lung and colorectal cancers in women of European ancestry.

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 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 Although lung cancer is known to be caused by environmental factors, it has also been shown to have genetic components, and the genetic etiology of lung cancer remains understudied. We previously identified a lung cancer risk locus on 6q23-25 using microsatellite data in families with a history of lung cancer. To further elucidate that signal, we performed targeted sequencing on nine of our most strongly linked families. Two-point linkage analysis of the sequencing data revealed that the signal was heterogeneous and that different families likely had different risk variants. Three specific haplotypes were shared by some of the families: 6q25.3-26 in families 42 and 44, 6q25.2-25.3 in families 47 and 59, and 6q24.2-25.1 in families 30, 33, and 35. Region-based logarithm of the odds scores and expression data identified the likely candidate genes for each haplotype overlap: ARID1B at 6q25.3, MAP3K4 at 6q26, and UTRN (6q24.1) and PHACTR2 (6q24.2). Further annotation was used to zero in on potential risk variants in those genes. All four genes are good candidate genes for lung cancer risk, having been linked to either lung cancer specifically or other cancers. However, this is the first time any of these genes has been implicated in germline risk. Functional analysis of these four genes is planned for future work. Significance: This study identifies four genes associated with lung cancer risk, which could help guide future lung cancer prevention and treatment approaches.

Abstract Background: Lung cancer kills more people than any other cancer in the United States. In addition to environmental factors, lung cancer has genetic risk factors as well, though the genetic etiology is still not well understood. We have performed whole exome sequencing on 262 individuals from 28 extended families with a family history of lung cancer. Methods: Parametric genetic linkage analysis was performed on these samples using two distinct analyses—the lung cancer only (LCO) analysis, where only patients with lung cancer were coded as affected, and the all aggregated cancers (AAC) analysis, where other cancers seen in the pedigree were coded as affected. Results: The AAC analysis yielded a genome-wide significant result at rs61943670 in POLR3B at 12q23.3. POLR3B has been implicated somatically in lung cancer, but this germline finding is novel and is a significant expression quantitative trait locus in lung tissue. Interesting genome-wide suggestive haplotypes were also found within individual families, particularly near SSPO at 7p36.1 in one family and a large linked haplotype spanning 4q21.3-28.3 in a different family. The 4q haplotype contains potential causal rare variants in DSPP at 4q22.1 and PTPN13 at 4q21.3. Conclusions: Regions on 12q, 7p, and 4q are linked to increased cancer risk in highly aggregated lung cancer families, 12q across families and 7p and 4q within a single family. POLR3B, SSPO, DSPP, and PTPN13 are currently the best candidate genes. Impact: Functional work on these genes is planned for future studies and if confirmed would lead to potential biomarkers for risk in cancer.

Abstract Lung cancer (LC) kills more people in the United States each year than any other cancer. While it is well known that a variety of environmental factors (particularly tobacco smoke) strongly increase the risk of LC, there are multiple associated genetic variants with small contributions to risk. High aggregation of LC within rare individual families suggest that there are high-risk genetic variants as well. However, these genetic risk factors for LC are under studied due to the rapid fatality of LC. We studied 28 highly aggregated extended high-risk familial lung cancer (HRFLC) families collected from eight different sites across the US. Whole exome sequencing was performed on 290 individuals from these families to identify potential risk variants for HRFLC using genetic linkage analysis. Quality control was performed on the sequence data, filtering on parameters such as read depth, genotype quality, missingness, and Mendelian inconsistencies. Identity-by-descent (IBD) values were also calculated to verify correct familial relationships. Quality control procedures left approximately 400,000 SNVs and indels for analysis.Parametric two-point linkage analysis was performed assuming an autosomal dominant mode of inheritance. Disease allele frequency was set to 1% with a penetrance of 80% for carriers and 1% phenocopy rate. While we did not identify any genome-wide significant variants across the 28 families, multiple suggestive variants were identified. The largest cluster of suggestive variants was located at 14q32 in the CATSPERB gene. Given the likely locus heterogeneity in LC (combined with the lack of power in some families), it is not surprising that none of the variants were significant across the families; looking at the individual family results proved more informative. Long haplotypes linked to LC risk were identified in multiple families. These long runs of positive linkages, which have little to no negative linkage evidence across them, are characteristic of true linkage signals in these types of analysis. Two of the most interesting linked regions were at 7p36.1 and 4q21.23-28.23. The 7p signal (observed in a single family) was genome-wide suggestive and located within the SSPO gene. SSPO has been implicated in breast and skin cancer (melanoma); it is a novel lung cancer signal. The 4q linkage (again observed in a single family) covers a large chunk of 4q and contains multiple potential candidate genes, however, the best candidate gene is PTPN13, a gene implicated in lung cancer but never in familial lung cancer. We are currently evaluating the individual family results of several other pedigrees and plan to perform additional analyses to confirm these linkages. Citation Format: Anthony Musolf, Haiming Sun, Bilal A. Moiz, Claudio W. Pikielny, Mariza de Andrade, Diptasri Mandal, Colette Gaba, Ping Yang, Yafang Li, Ming You, Richard K. Wilson, Elena Y. Kupert, Marshall W. Anderson, Ann G. Schwartz, Susan M. Pinney, Christopher I. Amos, Ramaswamy Govindan, Joan E. Bailey-Wilson. Familial lung cancer exhibits multiple novel linked haplotypes within pedigrees [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4176.

Abstract Lung cancer is the deadliest cancer in the United States, contributing approximately 25% of all cancer deaths. Lung cancer risk is well-documented to increase in response to environmental factors, particularly tobacco smoking. As one might expect for a complex trait, there is also a significant genetic risk component in lung cancer, though it is not very well studied. We obtained whole exome sequencing (WES) genotype data from Washington University on 204 subjects from 25 extended families. These individuals were recruited from families with a history of lung cancer and previous analyses showed these 25 families be informative. The purpose of this study is to identify potential risk variants for lung cancer by performing genetic linkage analysis. Quality control was performed on the sequence data, filtering on parameters such as depth (less than 10), genotype quality (less than 10), missingness, and Mendelian inconsistencies. Identity-by-descent (IBD) values were also calculated to verify correct familial relationships. Quality control procedures left approximately 500,000 SNVs and indels for analysis. We performed two-point parametric linkage analysis assuming an autosomal dominant mode of inheritance with a disease allele frequency of 1%, a 10% penetrance for carriers and a 1% penetrance for non-carriers. Two discrete sets of linkage analyses were performed. One was a variant-based analysis, which evaluated linkage between the phenotype and individual SNVs or indels. The second was gene-based analysis, which created a multi-allelic pseudomarker corresponding to a gene from haplotypes of rare variants (minor allele frequency &lt;= 0.01) located within that particular gene. Two-point linkage analysis was then performed on the pseudomarkers. While the variant-based analysis did not identify any genome-wide significant results, several were identified by the gene-based analysis. The highest HLOD scores were both greater than the genome-wide significance level of HLOD = 3.3 and were located on two regions of chromosome 1q: 1q42.13-43 and 1q21.2-21.1. In both regions, the significant HLOD scores clustered around known cancer genes. At the 1q42.13-43 region, signals centered on the cancer-implicated genes OBSCN and RYR2 (also a known mesothelioma gene), while at 1q21.2-21.1 the signals centered on five genes in the neuroblastoma breakpoint family (NBPF) genes, a cluster of recently duplicated genes. NBPF genes have previously been implicated in a variety of different cancers including lung cancer. We are currently performing additional analyses to corroborate the significant results, as well as examining the individual families to determine which families are driving the significant signals. Citation Format: Anthony M. Musolf, Haiming Sun, Bilal A. Moiz, Diptasri Mandal, Mariza de Andrade, Colette Gaba, Ping Yang, Yafang Li, Ming You, Elena Y. Kupert, Marshall W. Anderson, Ann G. Schwartz, Susan M. Pinney, Christopher I. Amos, Joan E. Bailey-Wilson. Whole exome sequencing identifies significantly linked regions on multiple chromosomes in families with a history of lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3276.

Abstract Lung cancer is the leading cancer-related cause of death in the United States. We had previously performed multipoint linkage analysis on families with a strong family history of lung cancer and found significant linkage to the 6q25 region. In order to find the source of that linkage signal, we performed targeted sequencing 6q23-6q27 on 75 individuals from the 9 most highly linked families. We performed two types of parametric linkage analysis, using an autosomal dominant mode of inheritance with 10% penetrance for carriers and a 1% phenocopy rate using a disease allele frequency of 1%. The first type was a two-point analysis using an Elston-Stewart algorithm. This approach did not lead to any genome-wide significant results but demonstrated significant heterogeneity throughout the nine families.In an effort to recover the power from the original multipoint analysis, we performed a regional linkage analysis using SEQLinkage and MERLIN. SEQLinkage built regional haplotypes that corresponded to a gene or an intergenic region - the regions were based on a customized map of our design. The regional markers were multiallelic which allowed for greater information content and were similar to the microsatellites that were used in the original multipoint analysis. The regional markers were then analyzed in a two-point linkage analysis via MERLIN. This allowed us to identify two genome-wide significant signals at PACRG-AS1 at 6q26 (HLOD = 3.4) and SAMD5 at 6q24 (3.3). The PACRG-AS1 is novel, though it is associated with the known lung cancer gene PARK2. SAMD5 may be involved in lung cancer cell proliferation. The heterogeneity of the two signals was particularly interesting. The PACRG-AS1 was driven primarily by two families. The SAMD5 marker was not being strongly driven by any of the families and appears to be a small cumulative effect across the nine families. In addition, several families have large, but non-significant LOD scores at other loci across the region. This further reinforces the locus heterogeneity within the region, and it is likely that our two significant signals here are not the only variants affecting the phenotype. We further attempted to localize the signals by running SEQLinkage using a custom map where genes are broken into exons and introns. Though this resulted in no significant markers, the highest signal was located in the intronic region of SAMD5 (HLOD = 3) and several other suggestive signals were localized to intronic regions of good candidate genes like SASH1 and ARID1B. Examination of predicted effects of the candidate regulatory variants using eQTL databases is ongoing. This is a region full of promising candidates, and it is likely that the two significant signals found here are just part of many that could be affecting lung cancer risk. We plan to do further analysis within the individual families to elucidate any more genes affecting this signal. Citation Format: Anthony M. Musolf, Claire L. Simpson, Bilal A. Moiz, Candace Middlebrooks, Mariza de Andrade, Diptasri Mandal, Colette Gaba, Ping Yang, Yafang Li, Ming You, Elena Y. Kupert, Marshall W. Anderson, Ann G. Schwartz, Susan M. Pinney, Christopher I. Amos, Joan E. Bailey-Wilson. A study in locus heterogeneity: Targeted sequencing analysis of 6q reveals multiple significant loci as the source of a previous linkage peak in familial lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 4290. doi:10.1158/1538-7445.AM2017-4290