Anderson, Marshall W.

Abstract Lung cancer (LC) is the leading cancer killer of Americans; an estimated 158,000 people will die in the U.S. from LC in 2016. While it is well-known that LC risk is affected by the environment, particularly tobacco smoking, there is a substantial genetic risk to LC also. We examined genotype data from Illumina HumanCore-12v1-0 array (297,000 SNPs) on 175 individuals in 25 extended families with a strong history of LC recruited by the Genetic Epidemiology of Lung Cancer Consortium. The purpose of this study is to determine which families are segregating a high-penetrance genetic risk haplotype so that the most informative families can be selected for DNA sequencing studies. Quality control was performed to remove SNPs and individuals with greater than 1% missingness as well as monomorphic SNPs. SNPs with Mendelian errors in more than one family were dropped. Identity-by-descent values were calculated to confirm correct familial relationships; 1 individual was dropped. After quality control we were left with approximately 245,000 SNPs for analysis. We performed three types of parametric linkage analyses using an autosomal dominant model with 40% penetrance in carriers and 1% penetrance in non-carriers. A disease allele frequency of 0.01 was used. Standard single variant two-point analysis between the disease and each marker was performed using TwoPointLods. Multipoint linkage analysis was performed using SimWalk2. We also performed regional-based linkage analyses using SEQLinkage and MERLIN. SEQLinkage builds a multiallelic regional marker (similar to a microsatellite) that corresponds to a gene or a portion of a gene. Two-point linkage analyses were then performed on the regional markers using MERLIN. We identified five loci that were genome-wide significant (HLOD score ≥ 3.3) from the regional-based linkage analyses on 18p11.23 (HLOD = 4.1), 2p22.2 (3.9), 14q13.1 (3.7), 16p13.1 (3.4), and 20q13.11 (3.4). It is particularly exciting that the scores centered on prospective cancer genes. Our highest score was centered on PTPRM, a protein tyrosine phosphatase on chromosome 18 that has been implicated as a LC oncogene. The signal on 20q13 also centered on another protein tyrosine phosphatase (PTPRT) that has been shown to be mutated in LC cells. The source of the signal on 16p13 was RNA binding protein RBFOX1, which has been shown to be deleted in LC cell lines and the source of the signal on 2p22 was LRP1B, a LDL receptor-related protein that is often inactivated in LC cells. The source of the final signal on 14q13 was NPAS3, a transcription factor that is a tumor suppressor in brain tumors. It should be noted that all previous evidence linking these genes to cancer was based on somatic mutations; this is the first time any of these genes has been shown to be significantly linked to germline disease risk in a family-based study. We plan to perform targeted sequencing on the linked regions to elucidate the exact causal variant. Citation Format: Joan E. Bailey-Wilson, Anthony M. Musolf, Claire L. Simpson, Mariza de Andrade, Diptasri Mandal, Colette G. Gaba, Ping Yang, Ming You, Elena Y. Kupert, Marshall W. Anderson, Ann G. Schwartz, Susan M. Pinney, Christopher I. Amos. Familial lung cancer is significantly linked to cancer-associated genes on five chromosomes [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 4268. doi:10.1158/1538-7445.AM2017-4268

Lung cancer is the deadliest cancer in the United States, killing roughly one of four cancer patients in 2016. While it is well-established that lung cancer is caused primarily by environmental effects (particularly tobacco smoking), there is evidence for genetic susceptibility. Lung cancer has been shown to aggregate in families, and segregation analyses have hypothesized a major susceptibility locus for the disease. Genetic association studies have provided strong evidence for common risk variants of small-to-moderate effect. Rare and highly penetrant alleles have been identified by linkage studies, including on 6q23–25. Though not common, some germline mutations have also been identified via sequencing studies. Ongoing genomics studies aim to identify additional high penetrance germline susceptibility alleles for this deadly disease.

Abstract Several papers have implicated mutations in the AT-rich interactive domain containing protein 1B (ARID1B) in carcinoma formation in multiple tissue types, including the lung. The ARID1B protein functions as part of SWI/SNF, a multi-subunit chromatin remodeling complex that repositions nucleosomes and may act as a tumor suppressor in cancer. ARID1B is located in a region in 6q that was previously shown to have strong evidence of linkage with lung cancer risk. As it is one of the more promising candidates within the target region, this study presents the viability of ARID1B as a candidate gene for familial lung cancer. Here, we chose individuals from families in our previous linkage study that showed the strongest linkage signal on 6q. 75 individuals from 9 families were sequenced for 37Mb of chromosome 6 from 130Mb to 167Mb using Illumina technology and a custom Agilent kit. Genotypes were imported into GoldenHelix SVS 7 for quality control and dropped if read depth < 10, quality score < 10 or a quality score to read depth ratio < 0.5. Variants with more than 1 Mendelian error were dropped. Variants with a single Mendelian inconsistency were dropped only for that family. Two-point parametric linkage analysis at theta = 0 was performed using R and an implementation of the Elston-Stewart algorithm in the paramlink package using a penetrance model of 0.01/0.1/0.1 (dd/Dd/DD) and a disease allele (D) frequency of 0.01. Family 102 showed strong evidence for linkage near ARID1B. The highest LOD score (1.34) out of all sequenced variants was located near this gene in a known conserved intergenic region between ARID1B and NOX3. Four other LOD scores > 0.7 were found in the intergenic region between these two genes. Two are transcription factor binding sites; two are conserved genomic elements. We also observed two linked intronic variants of ARID1B in this family; an enhancer (LOD = 0.89) and a transcribed region (LOD = 0.73). Two other families (47 and 59) displayed evidence of linkage in ARID1B. Here, the interest is not the value of the LOD scores (top scores were 0.64 and 0.54), but the centering of the linked haplotype around the gene. Family 47 has 60 variants with LODs between 0.64 and 0.61; 45 located within ARID1B or known regulatory regions. Family 59 has 61 variants with LODs between 0.54 and 0.51; 29 located within ARID1B or known regulatory regions The variants appear to be spread evenly throughout ARID1B, suggesting the presence of linked haplotypes. Analyses are currently underway to reconstruct possible haplotypes in these families. The results of this study show that variants in the ARID1B gene at least partially account for the linkage signal on 6q in three of the nine most strongly linked families in our cohort and thus we believe that ARID1B is a plausible candidate gene for further studies involving familial lung cancer. Replication studies using 27 new extended families are currently underway to confirm this initial result. Citation Format: Anthony M. Musolf, Claire L. Simpson, 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. Chromatin remodeling gene ARID1B is linked to familial lung cancer. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 2542.

Abstract Several studies have identified common genetic factors that influence susceptibility to lung cancer but no large studies have systematically studied familial lung cancers. The Genetic Epidemiology of Lung Cancer Consortium completed a genome-wide association study (GWAS) using samples from 757 cases and 796 controls, all of European American ancestry, and frequency matched on age and sex. Genotyping was performed by the Center for Inherited Disease Research (CIDR) using an Illumina Human OmniExpressExome-8v1 array. The most significantly associated SNPs were identified in the cases as a whole, and were also stratified by histology (adenocarcinoma and squamous cell carcinoma) and family history (0-1 and 2-5 affected family members). Analysis of all cases and controls identified a large number of SNPs around the CFTR gene that are strongly associated with lung cancer risk. Upon imputation with European Ancestry 1000 genomes data, the most significant SNP in this region is rs43034 (p = 3.89E-07, OR = 0.657) in the promoter region. In addition, an exonic missense SNP, rs213950 (p = 6.32E-07, OR = 0.663), and an intronic SNP with a possible splicing effect, rs213935 (p = 5.20E-07, OR = 0.660) were also identified as possible protective alleles. Analysis by histological subtype identified several rare and exonic SNPs that were specific by histology. We highlight two SNPs, rs114719990 (p = 1.194E-04, OR = 336.8, SETD2, c.3229A>G, p.Thr1077Ala) and rs199710487 (p = 1.124E-04, OR = 366, NTN5, c.656G>A, p.Arg219Gln), that were identified as strongly associated SNPs in cases with histologically confirmed adenocarcinoma, and two SNPs, rs144465058 (p = 1.029E-04, OR = 484.6, ANK3, c.6555G>T, p.Gln2185His) and rs140216112 (p = 1.087E-04, OR = 476, GRIK4, c.1013C>G, p.Ala338Gly), that were identified as top SNPs in cases with histologically confirmed squamous cell carcinoma. The SETD2 gene is thought to be a tumor suppressor gene and variants in this gene have previously been associated with leukemia and glioma (Zhu X et al. 2014, Fontebasso AM et al. 2013), and may now have implications in lung cancer, specifically adenocarcinoma development. In conclusion, we identified SNPs in the CFTR gene that may have a protective effect on lung cancer and several rare and exonic variants that were specific to histology subtype. SNPChrPositionGeneGenomic/ Nucleotide/ Amino acid changeOdds Ratiop-valueCFTR Protective Allelesrs430347117108449CFTRg.117108449T>C; 2.4kb upstream0.6573.89×10−7rs2139507117199533CFTRc.1408G>A; p.Val470Met0.6636.32×10−7rs2139357117178208CFTRc.869+1481A>G0.6605.20×10−7 Citation Format: Diana M. Toledo, Susan M. Pinney, Diptasri Mandal, Mariza de Andrade, Elena Kupert, Jennifer Franks, Colette Gaba, Claire L. Simpson, Ming You, Marshall W. Anderson, Joan E. Bailey-Wilson, Christopher I. Amos, Ann Schwartz. Genetic Epidemiology of Lung Cancer Consortium: genome-wide association study of familial lung cancer cases. [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 LB-189. doi:10.1158/1538-7445.AM2015-LB-189

Abstract Background: Recent genome-wide association studies (GWAS) have identified 28 susceptibility loci for Lung cancer (LC), 32 loci for smoking behavior, 10 loci for chronic obstructive pulmonary disease (COPD), 63 loci for pulmonary function and related phenotypes, totaling 107 GWAS susceptibility loci (as of November 2014). Given the common variants that have been found to be associated with LC in GWAS, targeted-region exome sequencing provides a cost-effective approach for further investigation of high-priority regions of the genome and has the great potential to identify rare causal variants in GWAS loci. Patients and Methods: Using an extreme phenotype approach, we selected 48 sporadic LC patients reporting heavy smoking histories with the presence of severe COPD in whom environmental factors are considered overwhelming, and 54 familial LC cases from families with at least three first-degree relatives with LC who are likely enriched for genomic signal, to search for the disease-causing rare mutations within the target 107 GWAS loci using exome sequencing data. Results: By focusing on exome profiles of these target 107 loci, we identified two disease-causing rare mutations on 6p21.32 TNXB p. Arg504His and 10q25.1 CCDC147 p.Arg696Cys. The homozygous Arg504His mutation presented in four familial and one sporadic LC cases. Notably, the minor allele frequency (MAF) of this variant in Caucasians from 1000 Genomes is very rare (<0.001). The heterozygous CCDC147 Arg696Cys variant were identified in two familial cases and one sporadic case. The MAF for this variant in the Caucasians from 1000 Genomes is 0.004. Both deleterious changes lead to truncation of the putative proteins, which disrupts the predicted domain and alters the cellular localization. We also observed several other suggestive rare mutations on G7C, TNS1, DBH1, and CHRNA5. Conclusion: Our target exome sequencing results demonstrate novel highly disruptive LC risk-conferring TNXB and CCDC147 mutations. Since the rare variants have appreciable effects on disease risk prediction, families and patients may warrant screening for mutations in TNXB and CCDC147 (if validated) to assess their potential LC risk. Citation Format: Yanhong Liu, Farrah Kheradmand, Michael Scheurer, Caleb Davis, David Wheeler, Edwin Silverman, Shete Sanjay, Spiridon Tsavachidis, Georgina Armstrong, Elena Kupert, Marshall Anderson, Yafang Li, Claudio Pikielny, Joan E. Bailey-Wilson, Ming You, Colette Gaba, Mariza DeAndrade, Diptasri Mandal, Claire Simpson, Susan Pinney, Christopher Amos, Margaret Spitz. Target exome sequencing for disease-causing rare mutations in familial and sporadic lung cancer. [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 4600. doi:10.1158/1538-7445.AM2015-4600

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 In order to select samples that are optimally informative for sequence analysis, studies for the Genetic Epidemiology Consortium, we implemented three strategies for analysis, taking into account the known family histories of lung cancer. Susceptibility for lung cancer arises from a complex interplay between largely unknown genetic factors and smoking behavior. In order to optimize our ability to identify genetic factors for lung cancer we have recruited individuals with extensive family histories of lung cancer and also collected information about their smoking behaviors. For this study, we have modeled family history to identify individuals to select for sequencing using three approaches. First, we used an expert opinion approach in which multiple investigators evaluated pedigree structures with family history to reach a consensus priority in selecting individuals. Second we developed an approximation that assumes a low penetrance among noncarriers and uses known cancer status in first degree relatives. Third, we used the cancer, smoking behaviors and family information along with a likelihood that models penetrance according to genetic susceptibility and an interaction with smoking to estimate the probabilities individuals carry susceptibility factors. In this analysis we compare results among these methods. We evaluate effectiveness of these methods using simulation-based approaches. Citation Format: Christopher Amos, Joan Bailey-Wilson, Susan Pinney, Maria de Andrada, Colette Gaba, Diptasri Mandal, Ming You, Margaret Spitz, Marshall Anderson. Optimal strategies for selecting samples for sequencing families with a common cancer: The Genetic Epidemiology of Lung Cancer Consortium. [abstract]. In: Proceedings of the Eleventh Annual AACR International Conference on Frontiers in Cancer Prevention Research; 2012 Oct 16-19; Anaheim, CA. Philadelphia (PA): AACR; Cancer Prev Res 2012;5(11 Suppl):Abstract nr A106.

Abstract Background: Genetic susceptibility for cancer can differ substantially among families. We use trait-related covariates to identify a genetically homogeneous subset of families with the best evidence for linkage in the presence of heterogeneity. Methods: We performed a genome-wide linkage screen in 93 families. Samples and data were collected by the familial lung cancer recruitment sites of the Genetic Epidemiology of Lung Cancer Consortium. We estimated linkage scores for each family by the Markov chain Monte Carlo procedure using SimWalk2 software. We used ordered subset analysis (OSA) to identify genetically homogenous families by ordering families based on a disease-associated covariate. We performed permutation tests to determine the relationship between the trait-related covariate and the evidence for linkage. Results: A genome-wide screen for lung cancer loci identified strong evidence for linkage to 6q23–25 and suggestive evidence for linkage to 12q24 using OSA, with peak logarithm of odds (LOD) scores of 4.19 and 2.79, respectively. We found other chromosomes also suggestive for linkages, including 5q31–q33, 14q11, and 16q24. Conclusions: Our OSA results support 6q as a lung cancer susceptibility locus and provide evidence for disease linkage on 12q24. This study further increased our understanding of the inheritability for lung cancer. Validation studies using larger sample size are needed to verify the presence of several other chromosomal regions suggestive of an increased risk for lung cancer and/or other cancers. Impact: OSA can reduce genetic heterogeneity in linkage study and may assist in revealing novel susceptibility loci. Cancer Epidemiol Biomarkers Prev; 19(12); 3157–66. ©2010 AACR.

Abstract A common variant on chromosomal region 15q24-25.1, marked by rs1051730, was found to be associated with lung cancer risk. Here, we attempted to confirm the second variant on 15q24-25.1 in several large sporadic lung cancer populations and determined what percentage of additional risk for lung cancer is due to the genetic effect of the second variant. SNPs rs1051730 and rs481134 were genotyped in 2,818 lung cancer cases and 2,766 controls from four populations. Joint analysis of these two variants (rs1051730 and rs481134) on 15q24-25.1 identified three major haplotypes (G_T, A_C, and G_C) and provided stronger evidence for association of 15q24-25.1 with lung cancer (P = 9.72 × 10−9). These two variants represent three levels of risk associated with lung cancer. The most common haplotype G_T is neutral; the haplotype A_C is associated with increased risk for lung cancer with 5.0% higher frequency in cases than in controls [P = 1.68 × 10−7; odds ratio (OR), 1.24; 95% confidence interval (95% CI), 1.14–1.35]; whereas the haplotype G_C is associated with reduced risk for lung cancer with 4.4% lower frequency in cases than in controls (P = 7.39 × 10−7; OR, 0.80; 95% CI, 0.73–0.87). We further showed that these two genetic variants on 15q24-25.1 independently influence lung cancer risk (rs1051730: P = 4.42 × 10−11; OR, 1.60; 95% CI, 1.46–1.74; rs481134: P = 7.01 × 10−4; OR, 0.81; 95% CI, 0.72–0.92). The second variant on 15q24-25.1, marked by rs481134, explains an additional 13.2% of population attributable risk for lung cancer. Cancer Res; 70(8); 3128–35. ©2010 AACR.

Abstract Cigarette smoking is the major cause for lung cancer, but genetic factors also affect susceptibility. We studied families that included multiple relatives affected by lung cancer. Results from linkage analysis showed strong evidence that a region of chromosome 6q affects lung cancer risk. To characterize the effects that this region of chromosome 6q region has on lung cancer risk, we identified a haplotype that segregated with lung cancer. We then performed Cox regression analysis to estimate the differential effects that smoking behaviors have on lung cancer risk according to whether each individual carried a risk-associated haplotype or could not be classified and was assigned unknown haplotypic status. We divided smoking exposures into never smokers, light smokers (<20 pack-years), moderate smokers (20 to <40 pack-years), and heavy smokers (≥40 pack-years). Comparing results according to smoking behavior stratified by carrier status, compared with never smokers, there was weakly increasing risk for increasing smoking behaviors, with the hazards ratios being 3.44, 4.91, and 5.18, respectively, for light, moderate, or heavy smokers, whereas among the individuals from families without the risk haplotype, the risks associated with smoking increased strongly with exposure, the hazards ratios being, respectively, 4.25, 9.17, and 11.89 for light, moderate, and heavy smokers. The never smoking carriers had a 4.71-fold higher risk than the never smoking individuals without known risk haplotypes. These results identify a region of chromosome 6q that increases risk for lung cancer and that confers particularly higher risks to never and light smokers. Cancer Res; 70(6); 2359–67