Background Endothelial dysfunction is an early stage of atherosclerosis. Single‐nucleotide polymorphisms ( SNP s) have been associated with vascular dysfunction, cardiac events, and coronary artery remodeling. We aimed to detect SNP s associated with endothelial dysfunction and determine whether these associations are sex specific. Methods and Results Six hundred forty‐three subjects without significant obstructive coronary artery disease underwent invasive coronary endothelial function assessment. We collected data from 1536 SNP s that had previously been associated with vasoreactivity, angiogenesis, inflammation, artery calcification, atherosclerotic risk factors, insulin resistance, hormone levels, blood coagulability, or with coronary heart disease. Coronary vascular reactivity was assessed by the percent change in coronary artery diameter ≤ −20% after an intracoronary bolus injection of acetylcholine on invasive coronary physiology study. SNP s significantly associated with coronary epicardial endothelial dysfunction were ADORA 1 , KCNQ 1 , and DNAJC 4 in the whole cohort, LPA , MYBPH , ADORA 3 , and PON 1 in women and KIF 6 and NFKB 1 in men ( P <0.01). Conclusions We have identified several significant SNP s that are associated with an increased risk of coronary endothelial dysfunction. These associations appear to be sex specific and may explain gender‐related differences in development of atherosclerosis.
Bayesian Variable Selection in Multilevel Item Response Theory Models with Application in Genomics
ObjectiveAt the cellular level, how excess adiposity promotes atherogenesis is not fully understood. One pathway involves secretion of adipokines that stimulate endothelial dysfunction through increased expression of adhesion molecules. However, the relationship of adiposity to adhesion molecules that promote atherosclerosis is largely unknown.MethodsLinear regression models were used to assess the sex‐specific associations of soluble cellular adhesion molecules (sP‐ and sL‐selectin, sICAM‐1, sVCAM‐1, and sHGF) and adiposity in 5,974 adults examined as part of the Multi‐Ethnic Study of Atherosclerosis (MESA). Adiposity measures included body mass index (BMI), waist‐to‐hip‐ratio (WHR), and computed tomography measures of subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT).ResultsThe mean age was 64 years and 52% were female. In multivariable models adjusting for traditional cardiovascular risk factors, sHGF was positively associated with BMI, WHR, and VAT in both males and females, and sP‐selectin with WHR and VAT in males. sVCAM‐1 was inversely associated with VAT in females only.ConclusionsOur results showed the relation of adiposity to soluble cellular adhesion proteins was similar across adiposity measures and for both sexes. However, the relationship between adiposity and sVCAM‐1 and P‐selectin may be modified by sex and the measure used to assess adiposity.
Pharmacogenomics of estrogens on changes in carotid artery intima-medial thickness and coronary arterial calcification: Kronos Early Estrogen Prevention Study
Prior to the initiation of menopausal hormone treatment (MHT), genetic variations in the innate immunity pathway were found to be associated with carotid artery intima-medial thickness (CIMT) and coronary arterial calcification (CAC) in women ( n = 606) enrolled in the Kronos Early Estrogen Prevention Study (KEEPS). Whether MHT might affect these associations is unknown. The association of treatment outcomes with variation in the same 764 candidate genes was evaluated in the same KEEPS participants 4 yr after randomization to either oral conjugated equine estrogens (0.45 mg/day), transdermal 17β-estradiol (50 μg/day), each with progesterone (200 mg/day) for 12 days each month, or placebo pills and patch. Twenty SNPs within the innate immunity pathway most related with CIMT after 4 yr were not among those associated with CIMT prior to MHT. In 403 women who completed the study in their assigned treatment group, single nucleotide polymorphisms (SNPs) within the innate immunity pathway were found to alter the treatment effect on 4 yr change in CIMT (i.e., significant interaction between treatment and genetic variation in the innate immunity pathway; P < 0.001). No SNPs by treatment effects were observed with changes of CAC >5 Agatston units after 4 yr. Results of this study suggest that hormonal status may interact with genetic variants to influence cardiovascular phenotypes, specifically, the pharmacogenomic effects within the innate immunity pathway for CIMT.
Abstract LB-189: Genetic Epidemiology of Lung Cancer Consortium: genome-wide association study of familial lung cancer cases
Cancer Research | August 5, 2015 | Volume: 75 | Issue: 15_Supplement | Pages: LB-189-LB-189
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.
Biology‐Driven Gene‐Gene Interaction Analysis of Age‐Related Cataract in the eMERGE Network
ABSTRACTBioinformatics approaches to examine gene‐gene models provide a means to discover interactions between multiple genes that underlie complex disease. Extensive computational demands and adjusting for multiple testing make uncovering genetic interactions a challenge. Here, we address these issues using our knowledge‐driven filtering method, Biofilter, to identify putative single nucleotide polymorphism (SNP) interaction models for cataract susceptibility, thereby reducing the number of models for analysis. Models were evaluated in 3,377 European Americans (1,185 controls, 2,192 cases) from the Marshfield Clinic, a study site of the Electronic Medical Records and Genomics (eMERGE) Network, using logistic regression. All statistically significant models from the Marshfield Clinic were then evaluated in an independent dataset of 4,311 individuals (742 controls, 3,569 cases), using independent samples from additional study sites in the eMERGE Network: Mayo Clinic, Group Health/University of Washington, Vanderbilt University Medical Center, and Geisinger Health System. Eighty‐three SNP‐SNP models replicated in the independent dataset at likelihood ratio test P < 0.05. Among the most significant replicating models was rs12597188 (intron of CDH1)–rs11564445 (intron of CTNNB1). These genes are known to be involved in processes that include: cell‐to‐cell adhesion signaling, cell‐cell junction organization, and cell‐cell communication. Further Biofilter analysis of all replicating models revealed a number of common functions among the genes harboring the 83 replicating SNP‐SNP models, which included signal transduction and PI3K‐Akt signaling pathway. These findings demonstrate the utility of Biofilter as a biology‐driven method, applicable for any genome‐wide association study dataset.
A comprehensive assessment of environmental exposures among 1000 North American patients with primary sclerosing cholangitis, with and without inflammatory bowel disease
SummaryBackgroundThe relationships between primary sclerosing cholangitis (PSC) and the environment are largely unknown.AimTo validate associations reported in previous studies and to identify novel environmental exposures among PSC patients.MethodsWe performed a multicenter, case–control analysis utilising self‐administered questionnaires. Responses between cases (n = 1000) and controls (n = 663) were compared using multivariable logistic regression adjusted for age and gender. The model was further stratified based on inflammatory bowel disease (IBD) status (with IBD n = 741 without IBD n = 259).ResultsSmoking was associated with PSC only when IBD was present (OR, 0.5; 95% CI 0.4–0.7) but not among those PSC patients without IBD (OR, 0.9; 95% CI 0.7–1.2). Compared to controls, women with PSC (irrespective of the presence of IBD) were less likely to have received hormone replacement therapy (HRT; OR, 0.5; 95% CI 0.4–0.7) and were more likely to have recurrent urinary tract infections (OR, 1.6; 95% CI 1.2–2.3). PSC patients regardless of gender or IBD status were less likely to eat fish (OR, 0.4; 95% CI 0.3–0.6) and grilled/barbecued meat (OR, 0.8; 95% CI 0.7–0.9). In contrast, PSC patients with and without IBD were more likely to consume steak/burgers that were more well done (OR, 1.3; 95% CI 1.2–1.5).ConclusionsIBD (rather than PSC) is associated with smoking. Women with PSC are more likely to have recurrent urinary tract infections and less likely to receive HRT. Dietary intake and methods of food preparation differ in PSC patients when compared to controls.
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.