Yang, Ping

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.

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 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

Abstract Recently, Wilson et al. (doi:10.1038/onc.2013.396) found that EYA4 was frequently inactivated biallelically in sporadic lung cancer, displays tumor suppressor gene-like properties and affects DNA repair, and that 5 of 17 common SNPs in this 0.3Mb region showed nominal association to familial lung cancer (FLC, p<0.05) although these associations are not genome-wide significant after correction for multiple testing.. However, this is an interesting candidate gene for FLC risk because it is located on chromosome 6q in a region where we previously showed evidence of linkage to FLC (PMID: 15272417, 20215501). Here, we sequenced 75 individuals from our 9 most strongly 6q-linked families, using Illumina technology and a custom Agilent kit to capture 37Mb of chromosome 6 from 130Mb to 167Mb. We closely examined a region around EYA4, between 133Mb and 135Mb. This region contains a number of other genes in addition to EYA4, including the Vanin cluster; ribosomal protein S12; the tumor suppressor gene TCF21 and its antisense RNA TARID which directs methylation and activation of TCF21; solute carriers SLC18B1 and SLC2A12 and the serum/glucocorticoid regulated kinase SGK1. To detect sharing of variants among affecteds in each family while allowing for phenocopies, singlepoint affected-only linkage analysis was performed on lung cancer affection status for all sequence variants using the Elston-Stewart algorithm implemented in the R package paramlink assuming penetrance of 0.01, 0.1 and 0.1 for dd/Dd/DD and a disease allele (D) frequency of 0.01, at a recombination fraction (theta) of 0. The highest total LOD score (summed over all 9 families) in the entire 37Mb targeted region was not near EYA4. In only one family was there any variant in EYA4 that appeared on the linked haplotype with close to the highest family-specific LOD score: Family 30 had a LOD of 0.6514 at several EYA4 variants (max LOD across region 0.6519). However, many variants across the region had similar LODs. In the remaining families, all EYA4 variants had LODs much lower than the maximum family-specific LOD score in the region. Thus, EYA4 variants are unlikely to explain our 6q linkage peak in these 6q-linked families. The second highest total LOD score for the entire 37Mb targeted region was found near EYA4 at 133.9Mb (LOD = 3.82), in the middle of the TARID antisense RNA gene. This variant appears to be entirely novel, as it is not present in any of the online variant databases, such as NHLBI's ESP or 1000 Genomes. In two families, the family LOD score at this locus was at or close to the highest in the region: in Family 44 the LOD score was 1.206 (max LOD across entire region 1.207) and in Family 35 the LOD score was 0.3958 (max LOD across region 0.3959). Therefore this variant in the TARID antisense RNA is a candidate in some of our families most strongly linked to chromosome 6q. Citation Format: Joan E. Bailey-Wilson, Claire L. Simpson, Susan M. Pinney, Mariza de Andrade, Colette Gaba, Ping Yang, Ming You, Elena Y. Kupert, Ann G. Schwartz, Diptasri Mandal, Christopher I. Amos, Marshall W. Anderson. Evaluation of EYA4 as a candidate risk locus in familial lung cancer families linked to 6q. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2757. doi:10.1158/1538-7445.AM2015-2757

Abstract Purpose: Lung cancer in never-smokers (LCINS) is increasingly recognized as a distinct disease from that in ever-smokers owing to substantial differences in etiology, clinical characteristics, and prognosis. Therefore, we aimed to identify prognostic markers specific for LCINS. Experimental Design: First, 11,930 single-nucleotide polymorphisms (SNP) in 904 inflammation-related genes were genotyped, and their associations with overall survival in 411 patients with LCINS at MD Anderson Cancer Center were analyzed. Next, validation of the top 27 SNPs in 311 patients with LCINS at Mayo Clinic was conducted. Results: Three SNPs (IL17RA:rs879576, BMP8A:rs698141, and STY:rs290229) were validated (P < 0.05), and two SNPs (CD74:rs1056400 and CD38:rs10805347) reached borderline significance (P = 0.08) in the Mayo Clinic population. We validated a survival-tree created in the MD Anderson population exploring gene–gene interactions in the Mayo Clinic population. This survival-tree stratified patients into subsets with significantly different risks of death: patients with the rs1056400_GG/rs698141_GA + AA genotype had significantly higher risk of death in both MD Anderson (HR:2.32, 95%CI: 1.58–3.41) and Mayo (HR:1.97, 95%CI: 1.11–3.50) populations compared with those with the rs1056400_GG/rs698141_GG or rs1056400_GA + AA genotype. We evaluated these five SNPs in 996 ever-smokers from MD Anderson and found no significant associations. Conclusions: Our study provides strong evidence that inflammation-related genetic variations can affect clinical outcomes in LCINS, which may lead to significant biologic insight into these outcomes. Clin Cancer Res; 18(21); 5983–91. ©2012 AACR.

AbstractLung cancer in lifetime never smokers is distinct from that in smokers, but the role of separate or overlapping carcinogenic pathways has not been explored. We therefore evaluated a comprehensive panel of 11,737 single-nucleotide polymorphisms (SNP) in inflammatory-pathway genes in a discovery phase (451 lung cancer cases, 508 controls from Texas). SNPs that were significant were evaluated in a second external population (303 cases, 311 controls from the Mayo Clinic). An intronic SNP in the ACVR1B gene, rs12809597, was replicated with significance and restricted to those reporting adult exposure to environmental tobacco smoke. Another promising candidate was an SNP in NR4A1, although the replication OR did not achieve statistical significance. ACVR1B belongs to the TGFR-β superfamily, contributing to resolution of inflammation and initiation of airway remodeling. An inflammatory microenvironment (second-hand smoking, asthma, or hay fever) is necessary for risk from these gene variants to be expressed. These findings require further replication, followed by targeted resequencing, and functional validation.Significance: Beyond passive smoking and family history of lung cancer, little is known about the etiology of lung cancer in lifetime never smokers that accounts for about 15% of all lung cancers in the United States. Our two-stage candidate pathway approach examined a targeted panel of inflammation genes and has identified novel structural variants that appear to contribute to risk in patients who report prior exposure to sidestream smoking. Cancer Discovery; 1(5): 420–9. ©2011 AACR.This article is highlighted in the In This Issue feature, p. 367