Blood lipid levels including low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and triglycerides (TG) are highly heritable. Genome-wide association is a promising approach to map genetic loci related to these heritable phenotypes.
In 1087 Framingham Heart Study Offspring cohort participants (mean age 47 years, 52% women), we conducted genome-wide analyses (Affymetrix 100K GeneChip) for fasting blood lipid traits. Total cholesterol, HDL-C, and TG were measured by standard enzymatic methods and LDL-C was calculated using the Friedewald formula. The long-term averages of up to seven measurements of LDL-C, HDL-C, and TG over a ~30 year span were the primary phenotypes. We used generalized estimating equations (GEE), family-based association tests (FBAT) and variance components linkage to investigate the relationships between SNPs (on autosomes, with minor allele frequency ≥10%, genotypic call rate ≥80%, and Hardy-Weinberg equilibrium p ≥ 0.001) and multivariable-adjusted residuals. We pursued a three-stage replication strategy of the GEE association results with 287 SNPs (P < 0.001 in Stage I) tested in Stage II (n ~1450 individuals) and 40 SNPs (P < 0.001 in joint analysis of Stages I and II) tested in Stage III (n~6650 individuals).
Long-term averages of LDL-C, HDL-C, and TG were highly heritable (h2 = 0.66, 0.69, 0.58, respectively; each P < 0.0001). Of 70,987 tests for each of the phenotypes, two SNPs had p < 10-5 in GEE results for LDL-C, four for HDL-C, and one for TG. For each multivariable-adjusted phenotype, the number of SNPs with association p < 10-4 ranged from 13 to 18 and with p < 10-3, from 94 to 149. Some results confirmed previously reported associations with candidate genes including variation in the lipoprotein lipase gene (LPL) and HDL-C and TG (rs7007797; P = 0.0005 for HDL-C and 0.002 for TG). The full set of GEE, FBAT and linkage results are posted at the database of Genotype and Phenotype (dbGaP). After three stages of replication, there was no convincing statistical evidence for association (i.e., combined P < 10-5 across all three stages) between any of the tested SNPs and lipid phenotypes.
Using a 100K genome-wide scan, we have generated a set of putative associations for common sequence variants and lipid phenotypes. Validation of selected hypotheses in additional samples did not identify any new loci underlying variability in blood lipids. Lack of replication may be due to inadequate statistical power to detect modest quantitative trait locus effects (i.e., <1% of trait variance explained) or reduced genomic coverage of the 100K array. GWAS in FHS using a denser genome-wide genotyping platform and a better-powered replication strategy may identify novel loci underlying blood lipids.
Increased circulating levels of hemostatic factors as well as anemia have been associated with increased risk of cardiovascular disease (CVD). Known associations between hemostatic factors and sequence variants at genes encoding these factors explain only a small proportion of total phenotypic variation. We sought to confirm known putative loci and identify novel loci that may influence either trait in genome-wide association and linkage analyses using the Affymetrix GeneChip 100K single nucleotide polymorphism (SNP) set.
Plasma levels of circulating hemostatic factors (fibrinogen, factor VII, plasminogen activator inhibitor-1, von Willebrand factor, tissue plasminogen activator, D-dimer) and hematological phenotypes (platelet aggregation, viscosity, hemoglobin, red blood cell count, mean corpuscular volume, mean corpuscular hemoglobin concentration) were obtained in approximately 1000 Framingham Heart Study (FHS) participants from 310 families. Population-based association analyses using the generalized estimating equations (GEE), family-based association test (FBAT), and multipoint variance components linkage analyses were performed on the multivariable adjusted residuals of hemostatic and hematological phenotypes.
In association analysis, the lowest GEE p-value for hemostatic factors was p = 4.5*10-16 for factor VII at SNP rs561241, a variant located near the F7 gene and in complete linkage disequilibrium (LD) (r2 = 1) with the Arg353Gln F7 SNP previously shown to account for 9% of total phenotypic variance. The lowest GEE p-value for hematological phenotypes was 7*10-8 at SNP rs2412522 on chromosome 4 for mean corpuscular hemoglobin concentration. We presented top 25 most significant GEE results with p-values in the range of 10-6 to 10-5 for hemostatic or hematological phenotypes. In relating 100K SNPs to known candidate genes, we identified two SNPs (rs1582055, rs4897475) in erythrocyte membrane protein band 4.1-like 2 (EPB41L2) associated with hematological phenotypes (GEE p < 10-3). In linkage analyses, the highest linkage LOD score for hemostatic factors was 3.3 for factor VII on chromosome 10 around 15 Mb, and for hematological phenotypes, LOD 3.4 for hemoglobin on chromosome 4 around 55 Mb. All GEE and FBAT association and variance components linkage results can be found at
Using genome-wide association methodology, we have successfully identified a SNP in complete LD with a sequence variant previously shown to be strongly associated with factor VII, providing proof of principle for this approach. Further study of additional strongly associated SNPs and linked regions may identify novel variants that influence the inter-individual variability in hemostatic factors and hematological phenotypes.
Systemic biomarkers provide insights into disease pathogenesis, diagnosis, and risk stratification. Many systemic biomarker concentrations are heritable phenotypes. Genome-wide association studies (GWAS) provide mechanisms to investigate the genetic contributions to biomarker variability unconstrained by current knowledge of physiological relations.
We examined the association of Affymetrix 100K GeneChip single nucleotide polymorphisms (SNPs) to 22 systemic biomarker concentrations in 4 biological domains: inflammation/oxidative stress; natriuretic peptides; liver function; and vitamins. Related members of the Framingham Offspring cohort (n = 1012; mean age 59 ± 10 years, 51% women) had both phenotype and genotype data (minimum-maximum per phenotype n = 507–1008). We used Generalized Estimating Equations (GEE), Family Based Association Tests (FBAT) and variance components linkage to relate SNPs to multivariable-adjusted biomarker residuals. Autosomal SNPs (n = 70,987) meeting the following criteria were studied: minor allele frequency ≥ 10%, call rate ≥ 80% and Hardy-Weinberg equilibrium p ≥ 0.001.
With GEE, 58 SNPs had p < 10-6: the top SNPs were rs2494250 (p = 1.00*10-14) and rs4128725 (p = 3.68*10-12) for monocyte chemoattractant protein-1 (MCP1), and rs2794520 (p = 2.83*10-8) and rs2808629 (p = 3.19*10-8) for C-reactive protein (CRP) averaged from 3 examinations (over about 20 years). With FBAT, 11 SNPs had p < 10-6: the top SNPs were the same for MCP1 (rs4128725, p = 3.28*10-8, and rs2494250, p = 3.55*10-8), and also included B-type natriuretic peptide (rs437021, p = 1.01*10-6) and Vitamin K percent undercarboxylated osteocalcin (rs2052028, p = 1.07*10-6). The peak LOD (logarithm of the odds) scores were for MCP1 (4.38, chromosome 1) and CRP (3.28, chromosome 1; previously described) concentrations; of note the 1.5 support interval included the MCP1 and CRP SNPs reported above (GEE model). Previous candidate SNP associations with circulating CRP concentrations were replicated at p < 0.05; the SNPs rs2794520 and rs2808629 are in linkage disequilibrium with previously reported SNPs. GEE, FBAT and linkage results are posted at .
The Framingham GWAS represents a resource to describe potentially novel genetic influences on systemic biomarker variability. The newly described associations will need to be replicated in other studies.
Osteoporosis is characterized by low bone mass and compromised bone structure, heritable traits that contribute to fracture risk. There have been no genome-wide association and linkage studies for these traits using high-density genotyping platforms.
We used the Affymetrix 100K SNP GeneChip marker set in the Framingham Heart Study (FHS) to examine genetic associations with ten primary quantitative traits: bone mineral density (BMD), calcaneal ultrasound, and geometric indices of the hip. To test associations with multivariable-adjusted residual trait values, we used additive generalized estimating equation (GEE) and family-based association tests (FBAT) models within each sex as well as sexes combined. We evaluated 70,987 autosomal SNPs with genotypic call rates ≥80%, HWE p ≥ 0.001, and MAF ≥10% in up to 1141 phenotyped individuals (495 men and 646 women, mean age 62.5 yrs). Variance component linkage analysis was performed using 11,200 markers.
Heritability estimates for all bone phenotypes were 30–66%. LOD scores ≥3.0 were found on chromosomes 15 (1.5 LOD confidence interval: 51,336,679–58,934,236 bp) and 22 (35,890,398–48,603,847 bp) for femoral shaft section modulus. The ten primary phenotypes had 12 associations with 100K SNPs in GEE models at p < 0.000001 and 2 associations in FBAT models at p < 0.000001. The 25 most significant p-values for GEE and FBAT were all less than 3.5 × 10-6 and 2.5 × 10-5, respectively. Of the 40 top SNPs with the greatest numbers of significantly associated BMD traits (including femoral neck, trochanter, and lumbar spine), one half to two-thirds were in or near genes that have not previously been studied for osteoporosis. Notably, pleiotropic associations between BMD and bone geometric traits were uncommon. Evidence for association (FBAT or GEE p < 0.05) was observed for several SNPs in candidate genes for osteoporosis, such as rs1801133 in MTHFR; rs1884052 and rs3778099 in ESR1; rs4988300 in LRP5; rs2189480 in VDR; rs2075555 in COLIA1; rs10519297 and rs2008691 in CYP19, as well as SNPs in PPARG (rs10510418 and rs2938392) and ANKH (rs2454873 and rs379016). All GEE, FBAT and linkage results are provided as an open-access results resource at .
The FHS 100K SNP project offers an unbiased genome-wide strategy to identify new candidate loci and to replicate previously suggested candidate genes for osteoporosis.
Intima-media thickness (IMT) of the common and internal carotid arteries is an established surrogate for atherosclerosis and predicts risk of stroke and myocardial infarction. Often IMT is measured as the average of these two arteries, yet they are believed to result from separate biological mechanisms. The aim of this study was to conduct a family-based genome-wide association study (GWAS) for IMT to identify polymorphisms influencing IMT and to determine if distinct carotid artery segments are influenced by different genetic components.
Methods and Results
IMT for the common and internal carotid arteries was determined through B-mode ultrasound in 772 Mexican Americans from the San Antonio Family Heart Study. A GWAS utilizing 931,219 single nucleotide polymorphisms (SNPs) was undertaken with six internal and common carotid artery IMT phenotypes utilizing an additive measured genotype model. The most robust association detected was for two SNPs (rs16983261, rs6113474, p=1.60e−7) in complete linkage disequilibrium on chromosome 20p11 for the internal carotid artery near wall, next to the gene PAX1. We also replicated previously reported GWAS regions on chromosomes 19q13 and 7q22. We found no overlapping associations between internal and common carotid artery phenotypes at p<5.0e0−6. The genetic correlation between the two carotid IMT arterial segments was 0.51.
This study represents the first large scale GWAS of carotid IMT in a non-European population and identified several novel loci. We do not detect any shared GWAS signals between common and internal carotid arterial segments but the moderate genetic correlation implies both common and unique genetic components.
intima-media thickness; carotid artery; GWAS; Hispanics
Brain magnetic resonance imaging (MRI) and cognitive tests can identify heritable endophenotypes associated with an increased risk of developing stroke, dementia and Alzheimer's disease (AD). We conducted a genome-wide association (GWA) and linkage analysis exploring the genetic basis of these endophenotypes in a community-based sample.
A total of 705 stroke- and dementia-free Framingham participants (age 62 +9 yrs, 50% male) who underwent volumetric brain MRI and cognitive testing (1999–2002) were genotyped. We used linear models adjusting for first degree relationships via generalized estimating equations (GEE) and family based association tests (FBAT) in additive models to relate qualifying single nucleotide polymorphisms (SNPs, 70,987 autosomal on Affymetrix 100K Human Gene Chip with minor allele frequency ≥ 0.10, genotypic call rate ≥ 0.80, and Hardy-Weinberg equilibrium p-value ≥ 0.001) to multivariable-adjusted residuals of 9 MRI measures including total cerebral brain (TCBV), lobar, ventricular and white matter hyperintensity (WMH) volumes, and 6 cognitive factors/tests assessing verbal and visuospatial memory, visual scanning and motor speed, reading, abstract reasoning and naming. We determined multipoint identity-by-descent utilizing 10,592 informative SNPs and 613 short tandem repeats and used variance component analyses to compute LOD scores.
The strongest gene-phenotype association in FBAT analyses was between SORL1 (rs1131497; p = 3.2 × 10-6) and abstract reasoning, and in GEE analyses between CDH4 (rs1970546; p = 3.7 × 10-8) and TCBV. SORL1 plays a role in amyloid precursor protein processing and has been associated with the risk of AD. Among the 50 strongest associations (25 each by GEE and FBAT) were other biologically interesting genes. Polymorphisms within 28 of 163 candidate genes for stroke, AD and memory impairment were associated with the endophenotypes studied at p < 0.001. We confirmed our previously reported linkage of WMH on chromosome 4 and describe linkage of reading performance to a marker on chromosome 18 (GATA11A06), previously linked to dyslexia (LOD scores = 2.2 and 5.1).
Our results suggest that genes associated with clinical neurological disease also have detectable effects on subclinical phenotypes. These hypothesis generating data illustrate the use of an unbiased approach to discover novel pathways that may be involved in brain aging, and could be used to replicate observations made in other studies.
Echocardiographic left ventricular (LV) measurements, exercise responses to standardized treadmill test (ETT) and brachial artery (BA) vascular function are heritable traits that are associated with cardiovascular disease risk. We conducted a genome-wide association study (GWAS) in the community-based Framingham Heart Study.
We estimated multivariable-adjusted residuals for quantitative echocardiography, ETT and BA function traits. Echocardiography residuals were averaged across 4 examinations and included LV mass, diastolic and systolic dimensions, wall thickness, fractional shortening, left atrial and aortic root size. ETT measures (single exam) included systolic blood pressure and heart rate responses during exercise stage 2, and at 3 minutes post-exercise. BA measures (single exam) included vessel diameter, flow-mediated dilation (FMD), and baseline and hyperemic flow responses. Generalized estimating equations (GEE), family-based association tests (FBAT) and variance-components linkage were used to relate multivariable-adjusted trait residuals to 70,987 SNPs (Human 100K GeneChip, Affymetrix) restricted to autosomal SNPs with minor allele frequency ≥0.10, genotype call rate ≥0.80, and Hardy-Weinberg equilibrium p ≥ 0.001.
We summarize results from 17 traits in up to 1238 related middle-aged to elderly men and women. Results of all association and linkage analyses are web-posted at . We confirmed modest-to-strong heritabilities (estimates 0.30–0.52) for several Echo, ETT and BA function traits. Overall, p < 10-5 in either GEE or FBAT models were observed for 21 SNPs (nine for echocardiography, eleven for ETT and one for BA function). The top SNPs associated were (GEE results): LV diastolic dimension, rs1379659 (SLIT2, p = 1.17*10-7); LV systolic dimension, rs10504543 (KCNB2, p = 5.18*10-6); LV mass, rs10498091 (p = 5.68*10-6); Left atrial size, rs1935881 (FAM5C, p = 6.56*10-6); exercise heart rate, rs6847149 (NOLA1, p = 2.74*10-6); exercise systolic blood pressure, rs2553268 (WRN, p = 6.3*10-6); BA baseline flow, rs3814219 (OBFC1, 9.48*10-7), and FMD, rs4148686 (CFTR, p = 1.13*10-5). Several SNPs are reasonable biological candidates, with some being related to multiple traits suggesting pleiotropy. The peak LOD score was for LV mass (4.38; chromosome 5); the 1.5 LOD support interval included NRG2.
In hypothesis-generating GWAS of echocardiography, ETT and BA vascular function in a moderate-sized community-based sample, we identified several SNPs that are candidates for replication attempts and we provide a web-based GWAS resource for the research community.
Type 2 diabetes mellitus (T2DM) is a major cardiovascular disease (CVD) risk factor. Identification of genetic risk factors for CVD is important to understand disease risk. Two recent genome-wide association study (GWAS) meta-analyses in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium detected CVD-associated loci.
Variants identified in CHARGE were tested for association with CVD phenotypes, including vascular calcification, and conventional CVD risk factors, in the Diabetes Heart Study (DHS) (n = 1208; >80% T2DM affected). This included 36 genotyped or imputed single nucleotide polymorphisms (SNPs) from DHS GWAS data. 28 coding SNPs from 14 top CHARGE genes were also identified from exome sequencing resources and genotyped, along with 209 coding variants from the Illumina HumanExome BeadChip genotype data in the DHS were also tested. Genetic risk scores (GRS) were calculated to evaluate the association of combinations of variants with CVD measures.
After correction for multiple comparisons, none of the CHARGE SNPs were associated with vascular calcification (p < 0.0014). Multiple SNPs showed nominal significance with calcification, including rs599839 (PSRC1, p = 0.008), rs646776 (CELSR2, p = 0.01), and rs17398575 (PIK3CG, p = 0.009). Additional COL4A2 and CXCL12 SNPs were nominally associated with all-cause or CVD-cause mortality. Three SNPs were significantly or nominally associated with serum lipids: rs3135506 (Ser19Trp, APOA5) with triglycerides (TG) (p = 5×10−5), LDL (p = 0.00070), and nominally with high density lipoprotein (HDL) (p = 0.0054); rs651821 (5′UTR, APOA5) with increased TGs (p = 0.0008); rs13832449 (splice donor, APOC3) associated with decreased TGs (p = 0.0015). Rs45456595 (CDKN2A, Gly63Arg), rs5128 (APOC3, 3′UTR), and rs72650673 (SH2B3, Glu400Lys) were nominally associated with history of CVD, subclinical CVD, or CVD risk factors (p < 0.010). From the exome chip, rs3750103 (CHN2, His204Arg/His68Arg) with carotid intima-medial thickness (IMT) (p = 3.9×10−5), and rs61937878 (HAL, Val549Met) with infra-renal abdominal aorta CP (AACP) (p = 7.1×10−5). The unweighted GRS containing coronary artery calcified plaque (CAC) SNPs was nominally associated with history of prior CVD (p = 0.033; OR = 1.09). The weighted GRS containing SNPs was associated with CAC and myocardial infarction (MI) was associated with history of MI (p = 0.026; OR = 1.15).
Genetic risk factors for subclinical CVD in the general population (CHARGE) were modestly associated with T2DM-related risk factors and CVD outcomes in the DHS.
Coronary artery calcified plaque; Type 2 diabetes mellitus; Cardiovascular disease; Genetic risk score
Family studies and heritability estimates provide evidence for a genetic contribution to variation in the human life span.
We conducted a genome wide association study (Affymetrix 100K SNP GeneChip) for longevity-related traits in a community-based sample. We report on 5 longevity and aging traits in up to 1345 Framingham Study participants from 330 families. Multivariable-adjusted residuals were computed using appropriate models (Cox proportional hazards, logistic, or linear regression) and the residuals from these models were used to test for association with qualifying SNPs (70, 987 autosomal SNPs with genotypic call rate ≥80%, minor allele frequency ≥10%, Hardy-Weinberg test p ≥ 0.001).
In family-based association test (FBAT) models, 8 SNPs in two regions approximately 500 kb apart on chromosome 1 (physical positions 73,091,610 and 73, 527,652) were associated with age at death (p-value < 10-5). The two sets of SNPs were in high linkage disequilibrium (minimum r2 = 0.58). The top 30 SNPs for generalized estimating equation (GEE) tests of association with age at death included rs10507486 (p = 0.0001) and rs4943794 (p = 0.0002), SNPs intronic to FOXO1A, a gene implicated in lifespan extension in animal models. FBAT models identified 7 SNPs and GEE models identified 9 SNPs associated with both age at death and morbidity-free survival at age 65 including rs2374983 near PON1. In the analysis of selected candidate genes, SNP associations (FBAT or GEE p-value < 0.01) were identified for age at death in or near the following genes: FOXO1A, GAPDH, KL, LEPR, PON1, PSEN1, SOD2, and WRN. Top ranked SNP associations in the GEE model for age at natural menopause included rs6910534 (p = 0.00003) near FOXO3a and rs3751591 (p = 0.00006) in CYP19A1. Results of all longevity phenotype-genotype associations for all autosomal SNPs are web posted at .
Longevity and aging traits are associated with SNPs on the Affymetrix 100K GeneChip. None of the associations achieved genome-wide significance. These data generate hypotheses and serve as a resource for replication as more genes and biologic pathways are proposed as contributing to longevity and healthy aging.
Breast and prostate cancer are two commonly diagnosed cancers in the United States. Prior work suggests that cancer causing genes and cancer susceptibility genes can be identified.
We conducted a genome-wide association study (Affymetrix 100K SNP GeneChip) of cancer in the community-based Framingham Heart Study. We report on 2 cancer traits – prostate cancer and breast cancer – in up to 1335 participants from 330 families (54% women, mean entry age 33 years). Multivariable-adjusted residuals, computed using Cox proportional hazards models, were tested for association with qualifying SNPs (70, 987 autosomal SNPs with genotypic call rate ≥80%, minor allele frequency ≥10%, Hardy-Weinberg test p ≥ 0.001) using generalized estimating equations (GEE) models and family based association tests (FBAT).
There were 58 women with breast cancer and 59 men with prostate cancer. No SNP associations attained genome-wide significance. The top SNP associations in GEE models for each trait were as follows: breast cancer, rs2075555, p = 8.0 × 10-8 in COL1A1; and prostate cancer, rs9311171, p = 1.75 × 10-6 in CTDSPL. In analysis of selected candidate cancer susceptibility genes, two MSR1 SNPs (rs9325782, GEE p = 0.008 and rs2410373, FBAT p = 0.021) were associated with prostate cancer and three ERBB4 SNPs (rs905883 GEE p = 0.0002, rs7564590 GEE p = 0.003, rs7558615 GEE p = 0.0078) were associated with breast cancer. The previously reported risk SNP for prostate cancer, rs1447295, was not included on the 100K chip. Results of cancer phenotype-genotype associations for all autosomal SNPs are web posted at .
Although no association attained genome-wide significance, several interesting associations emerged for breast and prostate cancer. These findings can serve as a resource for replication in other populations to identify novel biologic pathways contributing to cancer susceptibility.
A number of genetic variants have been discovered by recent genome-wide association studies for their associations with clinical coronary heart disease (CHD). However, it is unclear whether these variants are also associated with the development of CHD as measured by subclinical atherosclerosis phenotypes, ankle brachial index (ABI), carotid artery intima-media thickness (cIMT) and carotid plaque.
Ten CHD risk single nucleotide polymorphisms (SNPs) were genotyped in individuals of European American (EA), African American (AA), American Indian (AI), and Mexican American (MA) ancestry in the Population Architecture using Genomics and Epidemiology (PAGE) study. In each individual study, we performed linear or logistic regression to examine population-specific associations between SNPs and ABI, common and internal cIMT, and plaque. The results from individual studies were meta-analyzed using a fixed effect inverse variance weighted model.
None of the ten SNPs was significantly associated with ABI and common or internal cIMT, after Bonferroni correction. In the sample of 13,337 EA, 3,809 AA, and 5,353 AI individuals with carotid plaque measurement, the GCKR SNP rs780094 was significantly associated with the presence of plaque in AI only (OR = 1.32, 95% confidence interval: 1.17, 1.49, P = 1.08 × 10−5), but not in the other populations (P = 0.90 in EA and P = 0.99 in AA). A 9p21 region SNP, rs1333049, was nominally associated with plaque in EA (OR = 1.07, P = 0.02) and in AI (OR = 1.10, P = 0.05).
We identified a significant association between rs780094 and plaque in AI populations, which needs to be replicated in future studies. There was little evidence that the index CHD risk variants identified through genome-wide association studies in EA influence the development of CHD through subclinical atherosclerosis as assessed by cIMT and ABI across ancestries.
ankle brachial index; carotid artery intima-media thickness; carotid plaque; coronary heart disease; genetic association study; multiethnic populations; subclinical atherosclerosis
Circulating adiponectin has been associated with both clinical and subclinical cardiovascular disease (CVD). Variants of the adiponectin gene (ADIPOQ) are associated with clinical CVD, but little is known about associations with subclinical CVD. We studied the association of 11 ADIPOQ SNPs with common and internal carotid intima media thickness (cIMT), presence of coronary artery calcification (CAC), and CAC scores (in those with CAC) in 2847 participants in the Multi-Ethnic Study of Atherosclerosis (MESA). Participants were Caucasian (n=712), African-American (n=712), Chinese (n=718), and Hispanic (n=705). All models were adjusted for age, sex, and field site, and stratified by race/ethnic group. African-Americans with genotypes AG/GG of rs2241767 had 36% greater (95% CI (16%, 59%), p=0.0001) CAC prevalence; they also had a larger common cIMT (p=0.0043). Also in African-Americans, genotypes AG/AA of rs1063537 were associated with a 35% (95% CI (14%, 59%), p=0.0005) greater CAC prevalence. Hispanics with the AA genotype of rs11711353 had a 37% (95% CI (14%, 66%), p=0.0011), greater CAC prevalence compared to those with the GG genotype. Additional adjustment for ancestry in African-American and Hispanic participants did not change the results. No single SNP was associated with subclinical CVD phenotypes in Chinese or Caucasian participants. There appears to be an association between ADIPOQ SNPs and subclinical CVD in African-American and Hispanics. Replication as well as assessment of other ADIPOQ SNPs appears warranted.
The role of host genetics in the development of subclinical atherosclerosis in the context of HIV infected persons who are being treated with highly active antiretroviral therapy (HAART) is not well understood.
The present genome-wide association study (GWAS) is based on 177 HIV-positive Caucasian males receiving HAART who participated in the Fat Redistribution and Metabolic Change in HIV Infection (FRAM) Study. Common and internal carotid intima-media thicknesses (cIMT) measured by B-mode ultrasound were used as a subclinical measure of atherosclerosis. Single nucleotide polymorphisms (SNPs) were assayed using the Illumina HumanCNV370-quad beadchip. Copy Number Variants (CNV) were inferred using a hidden Markov Model (PennCNV). Regression analyses were used to assess the association of common and internal cIMT with individual SNPs and CNVs, adjusting for age, duration of antiretroviral treatment, and principal components to account for potential population stratification.
Two SNPs in tight linkage disequilibrium, rs2229116 (a missense, nonsynonymous polymorphism (IIe to Val)) and rs7177922, located in the Ryanodine receptor (RYR3) gene on chromosome 15 were significantly associated with common cIMT (p-value<1.61×10−7). The RYR gene family has been known to play a role in the etiology of cardiovascular disease and has been shown to be regulated by HIV TAT protein.
These results suggest that in the context of HIV infection and HAART, a functional SNP in a biologically plausible candidate gene, RYR3, is associated with increased common carotid IMT, which is a surrogate for atherosclerosis.
HIV; HAART; atherosclerosis; GWAS; intima-media thickness
Susceptibility to type 2 diabetes may be conferred by genetic variants having modest effects on risk. Genome-wide fixed marker arrays offer a novel approach to detect these variants.
We used the Affymetrix 100K SNP array in 1,087 Framingham Offspring Study family members to examine genetic associations with three diabetes-related quantitative glucose traits (fasting plasma glucose (FPG), hemoglobin A1c, 28-yr time-averaged FPG (tFPG)), three insulin traits (fasting insulin, HOMA-insulin resistance, and 0–120 min insulin sensitivity index); and with risk for diabetes. We used additive generalized estimating equations (GEE) and family-based association test (FBAT) models to test associations of SNP genotypes with sex-age-age2-adjusted residual trait values, and Cox survival models to test incident diabetes.
We found 415 SNPs associated (at p < 0.001) with at least one of the six quantitative traits in GEE, 242 in FBAT (18 overlapped with GEE for 639 non-overlapping SNPs), and 128 associated with incident diabetes (31 overlapped with the 639) giving 736 non-overlapping SNPs. Of these 736 SNPs, 439 were within 60 kb of a known gene. Additionally, 53 SNPs (of which 42 had r2 < 0.80 with each other) had p < 0.01 for incident diabetes AND (all 3 glucose traits OR all 3 insulin traits, OR 2 glucose traits and 2 insulin traits); of these, 36 overlapped with the 736 other SNPs. Of 100K SNPs, one (rs7100927) was in moderate LD (r2 = 0.50) with TCF7L2 (rs7903146), and was associated with risk of diabetes (Cox p-value 0.007, additive hazard ratio for diabetes = 1.56) and with tFPG (GEE p-value 0.03). There were no common (MAF > 1%) 100K SNPs in LD (r2 > 0.05) with ABCC8 A1369S (rs757110), KCNJ11 E23K (rs5219), or SNPs in CAPN10 or HNFa. PPARG P12A (rs1801282) was not significantly associated with diabetes or related traits.
Framingham 100K SNP data is a resource for association tests of known and novel genes with diabetes and related traits posted at . Framingham 100K data replicate the TCF7L2 association with diabetes.
Coronary artery calcium (CAC) and carotid intima-media thickness (IMT) are noninvasive measures of atherosclerosis that consensus panels have recommended as possible additions to risk factor assessment for predicting the probability of cardiovascular disease (CVD) occurrence.
To assess whether maximum carotid IMT or CAC (Agatston Score) is the better predictor of incident CVD.
Design, Setting, Patients
Prospective cohort study of 45–84 year-olds initially free of CVD (n = 6,698) in four ethnic groups, with standardized carotid IMT and CAC measures at baseline, in six field centers of the Multi-Ethnic Study of Atherosclerosis (MESA).
Main Outcome Measure(s)
Incident CVD events (coronary heart disease, stroke, and fatal CVD) over a maximum of 5.3 years of follow-up.
There were 222 CVD events during follow-up. CAC was associated more strongly than carotid IMT with risk of incident CVD. After adjustment for each other and traditional CVD risk factors, the hazard of CVD increased 2.1-fold (95% CI 1.8–2.5) for each standard deviation greater level of log-transformed CAC, versus 1.3-fold (95% CI 1.1–1.4) for each standard deviation greater maximum IMT. For coronary heart disease, the hazard ratios per standard deviation increment were 2.5-fold (95% CI 2.1–3.1) for CAC and 1.2-fold (95% CI 1.0–1.4) for IMT. An ROC analysis also suggested that CAC predicted incident CVD better than IMT did.
Although whether and how to clinically use bio-imaging tests of subclinical atherosclerosis remains a topic of debate, this study found that CAC predicts subsequent CVD events better than does carotid IMT.
Genome-wide association studies (GWASs) have identified a large number of variants (SNPs) associating with an increased risk of coronary artery disease (CAD). Recently, the CARDIoGRAM consortium published a GWAS based on the largest study population so far. They successfully replicated twelve already known associations and discovered thirteen new SNPs associating with CAD. We examined whether the genetic profiling of these variants improves prediction of subclinical atherosclerosis – i.e., carotid intima-media thickness (CIMT) and carotid artery elasticity (CAE) – beyond classical risk factors.
Subjects and Methods
We genotyped 24 variants found in a population of European ancestry and measured CIMT and CAE in 2001 and 2007 from 2,081, and 2,015 subjects (aged 30–45 years in 2007) respectively, participating in the Cardiovascular Risk in Young Finns Study (YFS). The Bogalusa Heart Study (BHS; n = 1179) was used as a replication cohort (mean age of 37.5). For additional replication, a sub-sample of 5 SNPs was genotyped for 1,291 individuals aged 46–76 years participating in the Health 2000 population survey. We tested the impact of genetic risk score (GRS24SNP/CAD) calculated as a weighted (by allelic odds ratios for CAD) sum of CAD risk alleles from the studied 24 variants on CIMT, CAE, the incidence of carotid atherosclerosis and the progression of CIMT and CAE during a 6-year follow-up.
CIMT or CAE did not significantly associate with GRS24SNP/CAD before or after adjusting for classical CAD risk factors (p>0.05 for all) in YFS or in the BHS. CIMT and CAE associated with only one SNP each in the YFS. The findings were not replicated in the replication cohorts. In the meta-analysis CIMT or CAE did not associate with any of the SNPs.
Genetic profiling, by using known CAD risk variants, should not improve risk stratification for subclinical atherosclerosis beyond conventional risk factors among healthy young adults.
Polymorphisms within the ICAM1 structural gene have been shown to influence circulating levels of soluble intercellular adhesion molecule -1 (sICAM-1) but their relation to atherosclerosis has not been clearly established. We sought to determine whether ICAM1 SNPs are associated with circulating sICAM-1 concentration, coronary artery calcium (CAC), and common and internal carotid intima medial thickness (IMT).
Methods and Results
3,550 black and white Coronary Artery Risk Development in Young Adults (CARDIA) Study subjects who participated in the year 15 and/or 20 examinations and were part of the Young Adult Longitudinal Study of Antioxidants (YALTA) ancillary study were included in this analysis. In whites, rs5498 was significantly associated with sICAM-1 (p < 0.001) and each G-allele of rs5498 was associated with 5% higher sICAM-1 concentration. In blacks, each C-allele of rs5490 was associated with 6 % higher sICAM-1 level; this SNP was in strong linkage disequilibrium with rs5491, a functional variant. Subclinical measurements of atherosclerosis in either year 15 or year 20 were not significantly related to ICAM1 SNPs.
In CARDIA, ICAM1 DNA segment variants were associated with sICAM-1 protein level including the novel finding that levels differ by the functional variant rs5491. However, ICAM1 SNPs were not strongly related to either IMT or CAC. Our findings in CARDIA suggest that ICAM1 variants are not major early contributors to subclinical atherosclerosis.
cell adhesion molecules; atherosclerosis; coronary calcium; genetics; inflammation
Little is known regarding the association of scavenger receptor class B type I (SCARB1) single nucleotide polymorphisms (SNPs) and subclinical atherosclerosis (SCA), particularly in subjects of different racial/ethnic backgrounds. We examined this relationship in the Multi-Ethnic Study of Atherosclerosis (MESA).
Methods and Results
Forty-three SCARB1 tagging SNPs were genotyped. Baseline examinations included fasting lipids and SCA phenotypes (coronary artery calcium [CAC], and common and internal carotid artery thickness [CCIMT and ICIMT]). Examining SNP associations with different SCA phenotypes across multiple racial/ethnic groups with adjustment for multiple covariates, we found the C allele of SNP rs10846744 was associated with higher CCIMT in African American (P=0.03), Chinese (P=0.02), European American (P=0.05), and Hispanic participants (P=0.03), and was strongly associated in pooled analyses (P=0.0002). The results also showed that the association of this SNP with CCIMT was independent of lipids and other well-established cardiovascular risk factors. Stratifying by sex, there appeared to be a strong association of rs10846744 with CCIMT in females, but no genotype-sex interactions were observed.
Variation in SCARB1 at rs10846744 was significantly associated with CCIMT across racial/ethnic groups in MESA.
genetics; atherosclerosis; cholesterol; lipids; prospective cohort study; genetic association
Genome-wide association (GWA) studies that use population-based association approaches may identify spurious associations in the presence of population admixture. In this paper, we propose a novel three-stage approach that is computationally efficient and robust to population admixture and more powerful than the family-based association test (FBAT) for GWA studies with family data.
We propose a three-stage approach for GWA studies with family data. The first stage is to perform linear regression ignoring phenotypic correlations among family members. SNPs with a first stage p-value below a liberal cut-off (e.g. 0.1) are then analyzed in the second stage that employs a linear mixed effects (LME) model that accounts for within family correlations. Next, SNPs that reach genome-wide significance (e.g. 10-6 for 34,625 genotyped SNPs in this paper) are analyzed in the third stage using FBAT, with correction of multiple testing only for SNPs that enter the third stage. Simulations are performed to evaluate type I error and power of the proposed method compared to LME adjusting for 10 principal components (PC) of the genotype data. We also apply the three-stage approach to the GWA analyses of uric acid in Framingham Heart Study's SNP Health Association Resource (SHARe) project.
Our simulations show that whether or not population admixture is present, the three-stage approach has no inflated type I error. In terms of power, using LME adjusting PC is only slightly more powerful than the three-stage approach. When applied to the GWA analyses of uric acid in the SHARe project of FHS, the three-stage approach successfully identified and confirmed three SNPs previously reported as genome-wide significant signals.
For GWA analyses of quantitative traits with family data, our three-stage approach provides another appealing solution to population admixture, in addition to LME adjusting for genetic PC.
Elevated fasting glucose level is associated with increased carotid intima-media thickness (IMT), a measure of subclinical atherosclerosis. It is unclear if this association is causal. Using the principle of Mendelian randomization, we sought to explore the causal association between circulating glucose and IMT by examining the association of a genetic risk score with IMT.
RESEARCH DESIGN AND METHODS
The sample was drawn from the Atherosclerosis Risk in Communities (ARIC) study and included 7,260 nondiabetic Caucasian individuals with IMT measurements and relevant genotyping. Components of the fasting glucose genetic risk score (FGGRS) were selected from a fasting glucose genome-wide association study in ARIC. The score was created by combining five single nucleotide polymorphisms (SNPs) (rs780094 [GCKR], rs560887 [G6PC2], rs4607517 [GCK], rs13266634 [SLC30A8], and rs10830963 [MTNR1B]) and weighting each SNP by its strength of association with fasting glucose. IMT was measured through bilateral carotid ultrasound. Mean IMT was regressed on the FGGRS and on the component SNPs, individually.
The FGGRS was significantly associated (P = 0.009) with mean IMT. The difference in IMT predicted by a 1 SD increment in the FGGRS (0.0048 mm) was not clinically relevant but was larger than would have been predicted based on observed associations between the FFGRS, fasting glucose, and IMT. Additional adjustment for baseline measured glucose in regression models attenuated the association by about one third.
The significant association of the FGGRS with IMT suggests a possible causal association of elevated fasting glucose with atherosclerosis, although it may be that these loci influence IMT through nonglucose pathways.
Carotid and coronary atherosclerosis are associated to each other in imaging and autopsy studies. We evaluated whether carotid artery plaque seen on carotid ultrasound can predict incident coronary artery calcification (CAC).
Materials and Methods
We repeated Agatston calcium score measurements in 5445 participants of the Multi-Ethnic Study of Atherosclerosis (MESA) (mean age 57.9 years; 62.9% female). Internal carotid artery lesions were graded as 0%, 1-24%, >25% diameter narrowing and intima-media thickness (IMT) was measured. Plaque was present for any stenosis > 0%. CAC progression was evaluated with multivariable relative risk regression in cases with CAC = 0 at baseline and with multivariable linear regression for CAC > 0 adjusting for cardiovascular risk factors, body mass index, ethnicity, and common carotid IMT.
CAC was positive at baseline in 2708/5445 (49.7%) participants and became positive in 458/2837 (16.1%) at mean interval of 2.4 years between repeat examinations. Plaque and ICA IMT were both strongly associated with presence of CAC. After statistical adjustment, presence of carotid artery plaque significantly predicted incident CAC with a relative risk(RR) of 1.37 (95% Confidence Intervals: 1.12, 1.67). Incident CAC was associated with ICA IMT with an RR of 1.13 (95% Confidence Intervals: 1.03, 1.25) for each mm increase. Progression of CAC was also significantly associated (p < 0.001) with plaque and ICA IMT.
In individuals free of cardiovascular disease, subjective and quantitative measures of carotid artery plaques by ultrasound imaging are associated with CAC incidence and progression.
Chromosome 9p21 has recently been shown to be a risk region for a broad range of vascular diseases. Since carotid intima-media thickness (IMT) and plaque are independent predictors for vascular diseases, the association between 9p21 and these two phenotypes was investigated.
Carotid segment-specific IMT and plaques were obtained in 1083 stroke- and myocardial infarction-free volunteers. We tested the genotypes and haplotypes of key single nucleotide polymorphisms (SNPs) on chromosome 9p21 for the associations with carotid IMT and plaque. Multivariate permutation analyses demonstrated that carriers of the T allele of SNP rs1333040 were significantly associated with thicker common carotid artery (CCA) IMT (p = 0.021) and internal carotid artery (ICA) IMT (p = 0.033). The risk G allele of SNP rs2383207 was associated with ICA IMT (p = 0.007). Carriers of the C allele of SNP rs1333049 were found to be significantly associated with thicker ICA IMT (p = 0.010) and the greater risk for the presence of carotid plaque (OR = 1.57 for heterozygous carriers; OR = 1.75 for homozygous carriers). Haplotype analysis showed a global p value of 0.031 for ICA IMT and 0.115 for the presence of carotid plaque. Comparing with the other haplotypes, the risk TGC haplotype yielded an adjusted p value of 0.011 and 0.017 for thicker ICA IMT and the presence of carotid plaque respectively. Further analyzing the data separated by sex, the results were significant only in men but not in women.
Chromosome 9p21 had a significant association with carotid atherosclerosis, especially ICA IMT. Furthermore, such genetic effect was in a gender-specific manner in the Han Chinese population.
Background and Purpose
Carotid intima-media thickness (IMT) is a surrogate marker of subclinical atherosclerosis and a strong predictor of stroke and myocardial infarction. The object of this study was to determine the association between carotid IMT and 702 single nucleotide polymorphisms in 145 genes.
B-mode carotid ultrasound was performed among 408 Hispanics from the Northern Manhattan Study. The common carotid artery IMT and bifurcation IMT were phenotypes of interest. Genetic effects were evaluated by the multivariate regression model adjusting for traditional vascular risk factors. For each individual, we calculated a gene risk score (GRS) defined as the total number of the significant single nucleotide polymorphisms in different genes. Subjects were then divided into 3 GRS categories using the 2 cutoff points: mean GRS ±1 SD.
We identified 6 significant single nucleotide polymorphisms in 6 genes for common carotid artery IMT and 7 single nucleotide polymorphisms in 7 genes for bifurcation IMT using the probability value of 0.005 as the significant level. There were no common significant genes for both phenotypes. The most significant genes were the tissue plasminogen activator (P=0.0005 for common carotid artery IMT) and matrix metallopeptidase-12 genes (P=0.0004 for bifurcation IMT). Haplotype analysis did not yield a more significant result. Subjects with GRS ≥9 had significantly increased IMT than those with GRS ≤5 (P<0.001). GRS was an independent predictor of both common carotid artery IMT (P=2.3×10−9) and bifurcation MT (P=7.2×10−8).
Multiple genes contributed to the variation in carotid IMT. IMT in different carotid segments may be regulated by different sets of susceptibility genes.
atherosclerosis; carotid intima-media thickness; genetics; polymorphism
Background and Purpose
Higher plasma concentrations of the endogenous nitric oxides synthase (NOS) inhibitor asymmetric dimethylarginine (ADMA) are associated with increased risk of cardio- and cerebrovascular events and death, presumably by promoting endothelial dysfunction and subclinical atherosclerosis. We hypothesized that plasma ADMA concentrations are positively related to common carotid artery intimal media thickness (CCA-IMT) and to internal carotid (ICA)/bulb-IMT.
We investigated the cross-sectional relations of plasma ADMA with CCA-IMT and ICA/bulb-IMT in 2958 Framingham Heart Study participants (mean age 58 years, 55% women).
In unadjusted analyses, ADMA was positively related to both CCA-IMT (β per SD increment 0.012, p<0.001) and ICA/bulb IMT (β per SD increment 0.059, p<0.001). In multivariable analyses (adjusting for age, sex, systolic blood pressure, antihypertensive treatment, smoking status, diabetes, body mass index (BMI), Total to HDL cholesterol ratio, log C-reactive protein, and serum creatinine), plasma ADMA was not associated with CCA-IMT (p=0.991), but remained significantly and positively related to ICA/bulb IMT (β per SD increment 0.0246, p=0.002).
In our large community-based sample, we observed that higher plasma ADMA concentrations were associated with greater ICA/bulb-IMT but not with CCA-IMT. These data are consistent with the notion that ADMA promotes subclinical atherosclerosis in a site-specific manner, with a greater proatherogenic influence at known vulnerable sites in the arterial tree.
Carotid Intimal Medial Thickness; Endothelium; Epidemiology; Risk Factors; Nitric Oxide
Arterial stiffness leads to left ventricular (LV) mass through non-atherosclerotic pathways in mice. In humans, a high ankle brachial index (ABI) indicates stiff peripheral arteries, and is associated with cardiovascular disease (CVD) events. Whether high ABI is associated with LV mass in humans, and whether this may reflect consequences of arterial stiffness, atherosclerosis, or both is unknown.
Among 4,972 MESA participants without clinical CVD, we used linear regression to evaluate the association of low (< 0.90) and high (>1.40 or incompressible) ABI with LV mass by cardiac magnetic resonance imaging (MRI). Intermediate ABIs served as the reference category. To determine the effect of subclinical atherosclerosis, models were adjusted for common and internal carotid intima media thickness (cIMT) and log-transformed coronary artery calcification (Ln[CAC+1]).
Compared to subjects with intermediate ABI, LV mass was higher with either low (2.70g/m2 higher, 95% CI 0.65–4.75) or high ABI (6.84 g/m2 higher, 95% CI 3.2–10.47) after adjustment for traditional CVD risk factors, kidney function, and CRP. However, further adjustment for cIMT and CAC substantially attenuated the association of low ABI with LVMI (1.24 g/m2 higher, 95% CI −0.84–3.33), whereas the association of high ABI was minimally altered (6.01 g/m2 higher, 95% CI 2.36–9.67).
High ABI is associated with greater LV mass; an association that is not attenuated with adjustment for subclinical atherosclerosis in non-peripheral arterial beds. High ABI may lead to greater LV mass through non-atherosclerotic pathways.
vascular stiffness; medial arterial calcification; left ventricular mass; heart failure; cardiovascular disease