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 2757: Evaluation of EYA4 as a candidate risk locus in familial lung cancer families linked to 6q
Cancer Research | August 5, 2015 | Volume: 75 | Issue: 15_Supplement | Pages: 2757-2757
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
Trans‐Ethnic Meta‐Analysis Identifies Common and Rare Variants Associated with Hepatocyte Growth Factor Levels in the Multi‐Ethnic Study of Atherosclerosis (MESA)
Annals of Human Genetics | May 22, 2015 | Volume: 79 | Issue: 4 | Pages: 264-274
SummaryHepatocyte growth factor (HGF) is a mesenchyme‐derived pleiotropic factor that regulates cell growth, motility, mitogenesis, and morphogenesis in a variety of cells, and increased serum levels of HGF have been linked to a number of clinical and subclinical cardiovascular disease phenotypes. However, little is currently known regarding which genetic factors influence HGF levels, despite evidence of substantial genetic contributions to HGF variation. Based upon ethnicity‐stratified single‐variant association analysis and trans‐ethnic meta‐analysis of 6201 participants of the Multi‐Ethnic Study of Atherosclerosis (MESA), we discovered five statistically significant common and low‐frequency variants: HGF missense polymorphism rs5745687 (p.E299K) as well as four variants (rs16844364, rs4690098, rs114303452, rs3748034) within or in proximity to HGFAC. We also identified two significant ethnicity‐specific gene‐level associations (A1BG in African Americans; FASN in Chinese Americans) based upon low‐frequency/rare variants, while meta‐analysis of gene‐level results identified a significant association for HGFAC. However, identified single‐variant associations explained modest proportions of the total trait variation and were not significantly associated with coronary artery calcium or coronary heart disease. Our findings indicate that genetic factors influencing circulating HGF levels may be complex and ethnically diverse.
Using Item Response Theory to Model Multiple Phenotypes and Their Joint Heritability in Family Data
ABSTRACTMany important complex diseases are composed of a series of phenotypes, which makes the disease diagnosis and its genetic dissection difficult. The standard procedures to determine heritability in such complex diseases are either applied for single phenotype analyses or to compare findings across phenotypes or multidimensional reduction procedures, such as principal components analysis using all phenotypes. However each method has its own problems and the challenges are even more complex for extended family data and categorical phenotypes. In this paper, we propose a methodology to determine a scale for complex outcomes involving multiple categorical phenotypes in extended pedigrees using item response theory (IRT) models that take all categorical phenotypes into account, allowing informative comparison among individuals. An advantage of the IRT framework is that a straightforward joint heritability parameter can be estimated for categorical phenotypes. Furthermore, our methodology allows many possible extensions such as the inclusion of covariates and multiple variance components. We use Markov Chain Monte Carlo algorithm for the parameter estimation and validate our method through simulated data. As an application we consider the metabolic syndrome as the multiple phenotype disease using data from the Baependi Heart Study consisting of 1,696 individuals in 95 families. We adjust IRT models without covariates and include age and age squared as covariates. The results showed that adjusting for covariates yields a higher joint heritability () than without co variates () indicating that the covariates absorbed some of the error variance.