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 37: Highly aggregated lung cancer families reveal a heterogeneous cause for a previous linkage signal on 6q
Abstract More Americans die every year of lung cancer than any other cancer. While the environmental risks for lung cancer are well understood, the genetic risk factors for this disease are not, even though it has been shown that lung cancer risk aggregates in families. Segregation analyses have confirmed the likelihood of rare, highly penetrant variants affecting lung cancer. Family studies offer a unique opportunity to identify such rare risk variants. We have previously identified a genome-wide significant risk locus for lung cancer using multipoint linkage analysis in highly aggregated lung cancer families. In this study, we performed targeted sequencing on this 6q23-6q27 region on 75 individuals from the 9 most strongly linked families. Parametric two-point linkage analysis using an autosomal dominant model with 10% penetrance for carriers and a 1% phenocopy rate and disease allele frequency of 1% was performed for each family. While we did not find any genome-wide significant results in any of the individual families, we did observe several interesting potential rare risk variants. First, we observed that the significant signal previously observed on 6q is recapitulated in these data but the location of the peak LOD is highly heterogeneous across families and it is likely that each family has a different risk gene. While most of the highest LOD scores were in noncoding variants, we did observe exonic variants with the highest LOD scores in two families (44 and 42). Both variants were rare, with a minor allele frequency (MAF) of approximately 1% and were in LPA and GPR126 respectively. Both these genes have been implicated somatically in lung cancer, however this is the first time that they have been implicated in the germline. While the other families have their highest LOD scores in noncoding variants, some of these genes are also good potential candidate genes including PARK2 (family 30), GRM1 (family 47), and PDE10A (family 102). All three genes have been implicated in lung cancer in some capacity, and this is the first time that GRM1 and PDE10A have been implicated as germline mutations. Again, all these variants are rare (MAF ⇐ 0.01) in the general population, which makes sense for a potential highly penetrant risk variant. In this study, we have elucidated that a previously identified risk locus on 6q25 is heterogeneous in the nine most strongly linked families, with different families appearing to be carrying different risk variants. Many of the top variants in each family are rare variants that are in good potential causal genes that have been identified for the first time here as germline risk variants for lung cancer. Further functional annotation is underway for these variants and some additional linkage analyses using other penetrance matrices may also be performed. Citation Format: Anthony M. Musolf, Claire L. Simpson, Bilal A. Moiz, Mariza de Andrade, Diptasri Mandal, Colette Gaba, Ping Yang, Ming You, Elena Y. Kupert, Marshall W. Anderson, Ann G. Schwartz, Susan M. Pinney, Christopher I. Amos, Joan E. Bailey-Wilson. Highly aggregated lung cancer families reveal a heterogeneous cause for a previous linkage signal on 6q [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 37.
Whole Exome Sequencing of Highly Aggregated Lung Cancer Families Reveals Linked Loci for Increased Cancer Risk on Chromosomes 12q, 7p, and 4q
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 4176: Familial lung cancer exhibits multiple novel linked haplotypes within pedigrees
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 3276: Whole exome sequencing identifies significantly linked regions on multiple chromosomes in families with a history of lung cancer
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 <= 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 4290: 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 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