Evaluation of 71 Coronary Artery Disease Risk Variants in a Multiethnic Cohort
Wangjing Ke1 
Kristin A. Rand2 
David V. Conti1 
Veronica W. Setiawan1 
Daniel O. Stram1 
Lynne Wilkens3 
Loic Le Marchand3 
Themistocles L. Assimes4,5 
Christopher A. Haiman1*- 1Department of Preventive Medicine, Keck School of Medicine of USC, Los Angeles, CA, United States
- 2Ancestry, San Francisco, CA, United States
- 3Epidemiology Program, University of Hawaii Cancer Center, Honolulu, HI, United States
- 4Department of Medicine, Stanford University School of Medicine, Stanford, CA, United States
- 5Cardiovascular Institute, Stanford University School of Medicine, Stanford, CA, United States
Background: Coronary heart disease (CHD) is the most common cause of death worldwide. Previous studies have identified numerous common CHD susceptibility loci, with the vast majority identified in populations of European ancestry. How well these findings transfer to other racial/ethnic populations remains unclear.
Methods and Results: We examined the generalizability of the associations with 71 known CHD loci in African American, Latino and Japanese men and women in the Multiethnic Cohort (6,035 cases and 11,251 controls). In the combined multiethnic sample, 78% of the loci demonstrated odds ratios that were directionally consistent with those previously reported (p = 2 × 10−6), with this fraction ranging from 59% in Japanese to 70% in Latinos. The number of nominally significant associations across all susceptibility regions ranged from only 1 in Japanese to 11 in African Americans with the most statistically significant association observed through locus fine-mapping noted for rs3832016 (OR = 1.16, p = 2.5×10−5) in the SORT1 region on chromosome 1p13. Lastly, we examined the cumulative predictive effect of CHD SNPs across populations with improved power by creating genetic risk scores (GRSs) that summarize an individual’s aggregated exposure to risk variants. We found the GRSs to be significantly associated with risk in African Americans (OR = 1.03 per allele; p = 4.1×10−5) and Latinos (OR = 1.03; p = 2.2 × 10−8), but not in Japanese (OR = 1.01; p = 0.11).
Conclusions: While a sizable fraction of the known CHD loci appear to generalize in these populations, larger fine-mapping studies will be needed to localize the functional alleles and better define their contribution to CHD risk in these populations.
Introduction
Coronary heart disease (CHD) is the most common, chronic, life-threatening illness in the United States, affecting more than 11 million people (1). A study with twins has estimated the genetic contribution to the variation in CHD mortality to be 0.57 and 0.38 in males and females, respectively (2). Genome-wide association studies (GWAS) have been conducted primarily in populations of European ancestry and have identified ~65 regions associated with CHD risk (3–11). Many of the CHD loci were identified in a large study of 22,233 case and 64,762 control of European ancestry in the CARDIoGRAMplusC4D consortium, which reported 46 genome-wide significant variants with odds ratios ranging from 1.01 to 2.08 and effect allele frequencies of 0.06–0.91 (9). More recently, 10 additional loci were reported from the same consortium in a genome-wide association study involving 61,289 cases and 126,310 controls subjects following imputation to the 1,000 Genomes Project reference panel (12). Genome-wide scans have also revealed 7 CHD risk loci in Asian populations (13–17). The known genetic risk variants for CHD are estimated to explain only 10–11% of the heritability of CHD (9, 12), suggesting that many additional genetic susceptibility loci remain to be discovered.
Several studies in Asian populations have reported successful replication of known CHD regions (17–20), with a reproducible disease association consistently noted with the 9p21 region. A limited number studies have been performed to investigate risk associated with CHD variants in minority groups such as African Americans or Latinos (21–28). In 2011, a GWAS in African Americans found a SNP, rs1859023, located at 7q21 near the PFTK1 gene to be significantly associated with CHD (22), however this finding has never been replicated in African Americans or any other racial/ethnic group. In a study of 8,090 African Americans (~700 CHD cases) that examined known CHD risk regions, only 9p21 was found to be associated with CHD (25). In a study of 8,201 African Americans (~550 CHD cases) (26), investigators found consistent direction of effects compared to studies of European ancestry for 23 of 44 (binomial p = 0.52) known loci with two nominally statistically significant (rs599839 at 1p13/SORT1 and rs579459 at 4p23/ABO). Genetic studies of CHD in Latino populations have been extremely limited. In a Costa Rican study that examined only 14 CHD SNPs in 1,898 cases with MI and 2,096 controls, 7 variants at 3 regions (SORT1, CXCL12, and 9p21) were found to be significantly associated with risk (29). Thus, additional studies are needed to understand the generalizability and relevance of the known CHD risk loci in populations of non-European ancestry.
In this context, the objective of this study was threefold. First, we wished to determine whether associations involving 71 known susceptibility variants of CHD from 65 independent regions generalize across African-American, Latino and Japanese men and women in the Multiethnic Cohort, a study that includes over 6,000 cases and 11,000 controls. Second, we evaluated common genetic variation across each susceptibility region in an attempt to identify variation that might better define the risk associations compared to the index variants in the multiethnic sample. Lastly, we constructed genetic risk scores (GRS) summarizing one’s degree of exposure to high risk alleles of CHD and evaluated to what degree this GRS contributes to population differences in CHD risk.
Methods
Study Population
The Multiethnic Cohort study (MEC) is a large prospective cohort study that was established between 1993 and 1996. The MEC includes primarily African Americans, Japanese American, Native Hawaiians, Latinos and European Americans living in Hawaii and California. Cohort members were recruited through Department of Motor Vehicle license files and supplemented by voter registration and Health Care Financing Administration (Medicare) files. Participating individuals were between 45 and 75 years of age, and completed a 26-page self-administered, detailed questionnaire at cohort entry (baseline data, 1993–1996). The questionnaire included basic demographic factors (including race/ethnicity and education), lifestyle factors (e.g., diet, medication use and smoking history), and chronic medical conditions. Follow-up questionnaires were also administered in years 1999 and 2003 which contained updates on participant’s CHD status and lifestyle factors.
Several nested case-control studies have been assembled in the MEC for GWAS of a number of cancer and non-cancer traits (30–32) including breast cancer, prostate cancer, and type-2 diabetes, mainly in populations of non-European ancestry. In the current study, we identified CHD cases and non-cases within these nested studies for the genetic analysis of CHD risk SNPs.
The MEC study obtained written informed consent from study participants for genetic analysis, approval from the Health Science Review Board (HSIRB) at the University of Southern California, and IRB certification permitting data sharing in accordance with the NIH Policy for Sharing of Data Obtained in NIH Supported or Conducted Genome-Wide Association Studies (GWAS). Genetic data for the MEC is available on dbGAP (phs000517.v3.p1, phs000851.v1.p1, phs000356.v2.p1, phs000306.v4.p1, phs000683.v1.p1)
CHD Case/Control Definitions
CHD cases were identified through linkage of the MEC to the California Hospital Discharge Data (1990–2012) (CHDD) and the Centers for Medicare and Medicaid Services (CMS) claim files (MedPAR, outpatient) (1999–2011). Hospital discharge information was not available for the subjects from Hawaii which included 76.6% of the Japanese men and women. A CHD case was defined as having ischemic heart disease under ICD-9 codes (DX 410–414), by the principal or first diagnosis code and the principal or first procedure code. We also included cases with a primary cause of death due to myocardial infarction (ICD-9 DX410, ICD-10 I21), or other CHD conditions (ICD-9 DX411–414, ICD-10 I20, I22–25). Both prevalent and incident CHD cases were included in this study. Of the 6,035 CHD cases identified, 1,146 were identified from their baseline questionnaires at the time of enrollment in the MEC study, and a majority of these prevalent cases (1,122, 97.9%) were also identified from CHDD or Medicare.
Controls in this study were subjects with no history of heart attack or angina based on the baseline questionnaire or all subsequent follow-up questionnaires. Those taking nitrates at blood draw in subsequent examinations were also excluded. Individuals with non-primary CHD diagnosis codes (i.e., 2–24) from the CHDD and Medicare data were excluded from being either a case or control. A total of 11,251 controls were selected, of which 8,307 had at least one previous Medicare or CHDD claim (and thus would have been identified as a case). A sensitivity analysis using controls with definite claim information was performed.
Genotyping and Quality Control
We utilized genetic data generated from case-control studies in the MEC of breast cancer, prostate cancer, and type 2 diabetes in African Americans (2,976 males and 3,539 females), Japanese Americans (2,530 males and 2,132 females), and Latinos (3,340 males and 2,769 females). Genotyping was conducted using the Illumina platform with different arrays, including the Human 1M-Duo v3.0 BeadChip (31, 32), HumanOmni2.5-Quad BeadChip (33), Human 660W-Quad BeadChip (34), and the Cardio-MetaboChip (35) (Table S1). We used the following exclusion criteria to remove samples whose genetic or phenotypic data were questionable: (1) unknown replicates across studies, (2) call rates < 95%, (3) samples with mismatched gender, such as male samples with >10% mean heterozygosity of SNPs on the X chromosome and/or <10% mean intensity of the Y chromosome; or female samples with <15% mean heterozygosity of SNPs on the X chromosome and/or similar mean allele intensities of SNPs on the X and Y chromosomes, (4) ancestry outliers (>4 standard deviations from the mean of the first or second principal component), and (5) first degree relatives.
A subset of 2,717 African Americans (879 CHD cases and 1,838 controls) and 1,184 Japanese Americans (302 CHD cases and 882 controls) genotyped with the Cardio-MetaboChip were missing data for 20 of the 71 SNPs; these subjects were excluded from the risk score analysis.
SNP Imputation and Principal Components Analysis
All samples except for the African-American and Japanese samples genotyped with the Metabochip were imputed using the software IMPUTE2, based on build 37 (hg19) coordinates, to the 1000 Genome Project data phase 1 v3. Principal components were calculated by study in smartpca from EIGENSOFT (36) using a random selection of 10,000 SNPs across the genome (MAF >5% and call rate >95%).
Statistical Analysis
The log-additive effect of each SNP on CHD risk was estimated in PLINK using unconditional logistic regression adjusted for age, sex, BMI and the first 10 principal components to account for potential population stratification (37). All analyses were stratified by ethnicity, disease status (i.e., breast cancer, prostate cancer, or type 2 diabetes disease status). METAL was used to combine the results within and across populations, which included 18 case-control strata in the overall meta-analysis of all populations. For SNPs that were imputed, all were imputed with an IMPUTE2 INFO score >0.8 in each study and population. SNPs rs11752643 and rs3782886 in African Americans, rs180803 in Latinos, and rs6544713, rs4252120, rs2023938, rs3918226, rs3184504, and rs9982601 in Japanese had a minor allele frequency less than 1% and were not included in the ethnic-specific analysis. The cross-ethnic meta-analysis was performed on SNPs observed in at least two ethnic groups.
In addition to testing of the index SNP, we examined regional replication of the signal through testing SNPs in linkage disequilibrium (LD) with the index SNP in European ancestry groups (r2 ≥0.4 in EUR 1000 Genomes Project). Haploview (38) was used to assess pairwise tag SNPs among bins of markers in the AFR population [tagging r2 ≥0.8 for SNPs with a MAF >1% based on 1000 Genomes Project data (39)]. For each region, an alpha threshold of significance was set at 0.05 divided by the number of tag SNPs in AFR. We considered evidence of replication to be present in a region when one or more SNPs in LD with the index SNP had a p-value that was lower than the region-defined alpha threshold. For imputed SNPs, only those imputed with high quality (IMPUTE2 INFO score >0.8) were included in the regional replication testing. The regional association plots were generated with the LocusZoom program (40).
We also examined the aggregate effect of the CHD risk loci. Three genetic risk scores (GRS) were calculated for each individual: (1) An unweighted GRS comprised of risk-associated alleles from the 71 CHD SNPs, (2) a modified unweighted GRS (I) that substitutes the index SNP with the lead SNP reaching region-wide significance within a specific race/ethnic group of each known CHD locus, and (3) a modified unweighted GRS (II) similar to I but substituting index SNPs with the leading SNPs in each region from our cross-ethnic meta-analysis. The risk alleles for the substitution SNPs were determined based on their observed effects in our study. As outlined above, subjects genotyped with the non-GWAS Metabochip were excluded from the risk score analysis because of missing data on 20 SNPs. The risk score distributions across ethnic groups were compared using a two-sided t-test. The association of genetic risk scores with CHD were evaluated within each ethnicity in a logistic regression model adjusted for age, sex, BMI, and the first 10 principal components. Of the 71 SNPs selected, only one pair (rs16986953 and rs2123536) from TTC32-WDR35 was correlated. Since the association between rs2123536 and CHD was only observed in a Chinese population (16), both SNPs were kept in the GRS analysis.
Within each population, statistical power for each SNP was calculated in the R package “gap,” (41) using the allele frequency in each racial/ethnic group, and the documented OR from the literature. The allele frequency for the multiethnic sample was weighted by the sample size of each ethnic group. The power for detecting rare and common alleles within each ethnic group was calculated using QUANTO (42).
Results
Descriptive characteristics of the 6,035 CHD cases and 11,251 controls stratified by sex and race/ethnicity are presented in Table 1. We analyzed a total of 2,376 African-American cases and 4,139 controls, 2,291 Latino cases and 3,818 controls, and 1,368 Japanese cases and 2,294 controls. In general, compared to controls, CHD cases were slightly older at cohort entrance, were heavier in all three ethnic groups and were more likely to have ever smoked than controls in all three ethnic groups (Table 1). The associations of BMI and smoking with CHD were similar when further stratified by prevalent conditions, including prostate cancer, breast cancer, and diabetes (Table S2).
Table 1. Descriptive Characteristics of CHD Cases and Controls.
We had a priori greater than 80% power to detect reported per allele effect sizes for 6 out of 71 SNPs in African Americans, 9 out of 71 SNPs in Latinos, and 9 out of 71 SNPs in Japanese Americans and 16 out of 71 SNPs when combining samples from all three ethnic groups (Figure S1). Given the sample size in each ethnic group, we had 28.5% power to detect an OR of 1.12 (mean OR from the selected index SNPs) for a rare (MAF = 0.05) allele and 71.6% power to detect the same OR for a common (MAF = 0.20) allele in African Americans; we had 27.3% power to detect OR of 1.12 for a rare allele and 69.3% power for a common allele in Latinos; and 20% power for a rare allele and 52.7% power for a common allele in Japanese Americans.
We examined evidence of replication for 71 CHD variants from 65 regions (Table S3). Among these variants, 69 in African Americans, 70 in Latinos, and 65 in Japanese Americans had a MAF >1% and were included in the analysis. Compared to the null expectation that one-half of the examined SNPs show consistent direction of effects as previously reported, 65.2% (45 of 69, binomial p = 0.008) SNPs in African Americans, 70.0% (49 of 70, binomial p = 5.5 × 10−4) in Latinos, 58.5% (38 of 65, binomial p = 0.11) in Japanese, and 77.5% (55 of 71, binomial p = 2.0 × 10−6) in the combined multiethnic sample had the same direction of association as previously reported. In African Americans, nominally statistically significant associations (p < 0.05) and consistent directional effects were observed for 11 index SNPs in PPAP2B, SORT1, IL6R, REST-NOA1, BTN2A1, SLC22A3-LPAL2-LPA, 9p21, CXCL12, SH2B3, and KCNE2. In Latinos, nominal evidence of association (p < 0.05) and consistent directional effects were observed with 8 index SNPs at SORT1, APOB, NOS3, LPL, ZHF259-APOA5-APOA1, MFGE8-ABHD2, FURIN-FES, and BCAS3. In Japanese, only 1 index SNP at 9p21 was nominally significant and directionally consistent. In the combined multiethnic sample, 10 index SNPs at PPAP2B, SORT1, IL6R, REST-NOA1, EDNRA, PHACTR1, BTN2A1, NOS3, 9p21, and CXCL12 were directionally consistent and nominally statistically significant.
We observed evidence of regional replication for 6 regions in African Americans, 3 in Latinos, 1 in Japanese Americans, and 10 in the combined sample when examining SNPs correlated with the index SNPs (Table S4; see Methods). The previously reported index SNP in four of the 10 regions was not significant at the 0.05 level, but correlated SNPs with p-values smaller than the region specific significance levels were detected in these four regions: SLC22A4-SLC22A5 and RAI1-PEMT-RASD1 in African Americans, TTC32-WDR35 in the multiethnic analysis, and APOE-APOC1 in Latinos and the multiethnic sample.
The most statistically significant association was observed at the SORT1 locus (Figure 1). Two index SNPs in complete LD (rs602633 and rs599839) were initially reported from GWAS in European ancestry populations. The index SNP rs602633 was associated with risk in African Americans (OR = 1.13; p = 0.004), Latinos (OR = 1.11; p = 0.04), and in the cross-ethnic meta-analysis (OR = 1.11; p = 7.8×10−4), but not in Japanese Americans (OR = 1.01, p = 0.88). The most significant association in the region was with variant rs3832016 (OR = 1.16; p = 2.5×10−5 in the multiethnic sample), an INDEL (−/T) in high LD with rs602633 in EUR (r2 = 0.96) and with a MAF of 0.35 in African Americans, 0.20 in Latinos, and 0.07 in Japanese Americans. A previous fine-mapping study of the SORT1 region at 1p13 implicated a nearby non-coding polymorphism (rs12740374) to be the likely functional variant and to affect lipoprotein metabolism (43). SNP rs12740374 is in high LD not only with the index SNP rs602633 (r2 = 0.90 in EUR) but also with rs3832016 (r2 = 0.94 in EUR). Variant rs12740374 was less strongly associated with risk in the current study (p = 0.008 in African Americans with MAF = 0.25, p = 0.08 in Latinos with MAF = 0.20, p = 0.93 in Japanese Americans with MAF = 0.07, and p = 0.003 in the combined multiethnic analysis).
Figure 1. Results for the SORT1 region on chromosome 1p13 from the multiethnic meta-analysis. The r2 shown is for the EUR group in the 1000 Genomes Project relative to the index SNP rs602633. The stronger signal, rs3832016, is also highlighted. This regional association plot was generated with the LocusZoom plot (40).
Other regions where evidence of regional replication was observed in African Americans include PPAP2B (rs72664341, p = 0.00018), SLC22A4-SLC22A5 (rs17689550, p = 0.006), SLC22A3-LPAL2-LPA (rs4709431, p = 0.0077), SH2B3 (rs10774625, p = 0.0047) and RAI1-PEMT-RASD1 (rs9899364, p = 4.5 × 10−4). In Latinos, evidence for regional replication was observed at MFGE8-ABHD2 (rs8037001, p = 0.0017), FURIN-FES (rs8182016, p = 1.1 × 10−4), and APOE-APOC1 (rs7412, p = 0.0043). In the Japanese, regional replication was only observed at 9p21 (rs10811656, p = 0.0015). Five of the 10 regions that replicated in the multiethnic analysis were also significant in ethnic-specific analyses, whereas the remaining 5 regions were detected with significant regional associations in one or more of the ethnic-specific populations alone (TTC32-WDR35, APOB, EDNRA, PHACTR1 and BCAS3; Table S4).
Genetic risk scores (GRSs) were used to compare the distribution of genetic risk between populations. Japanese Americans carried, on average, more risk alleles (70.26 ± 4.61, mean ± SD) in comparison to African Americans and Latinos (67.03 ± 4.73 and 68.37 ± 5.11, respectively) (Table 2; Table S5). The greater number of risk alleles resulted in the distribution of the GRS to be shifted to the right in Japanese Americans compared to African Americans and Latinos (Figure 2). The distribution of the GRS was slightly higher in cases than in controls for every group (two-sided t-test, AA p = 4.4 × 10−5, LA p = 4.5 × 10−7, and JA p = 0.28). Only minor changes in the distribution of the GRS were noted when we included regionally significant leading SNPs from each ancestry (modified risk score I), or from the cross-ethnic meta-analysis (modified risk score II) (Table S5). The average risk scores remained highest in Japanese Americans whereas differences between African Americans and Latinos were reduced, especially when comparing CHD cases from these two ethnic groups (modified risk score II, p = 0.14).
Figure 2. A comparison of the aggregate allele count risk score for cases and controls in each race/ethnic group.
Table 2. Associations of the genetic risk score with CHD by ethnicity.
The unweighted risk scores were statistically significantly associated with CHD risk in African Americans (per allele OR = 1.03, p = 4.1 × 10−5) and Latinos OR = 1.03, (p = 2.2 × 10−8), but only weakly associated with CHD risk in Japanese Americans (OR = 1.01, p = 0.11) (Table 2). When comparing individuals within GRSs in the top quartile to individuals in the bottom quartile, we found both African Americans (OR = 1.40) and Latinos (OR = 1.39) to have a statistically significant ~40% increase in risk (Table 2). The analogous risk was lower (~10%) and not significant in Japanese-Americans (OR = 1.09). Results were similar for the modified risk scores (Table 2).
To evaluate the effect of existing conditions on the results, we repeated the analysis excluding cancer or diabetes cases; the ORs were comparable to those observed in each ethnic group and in the entire sample (Table S6).
A sensitivity analysis was also performed on the selected index SNPs using controls refined to those with medical claims from Medicare or CHDD. Despite the loss of statistical power due to smaller sample size, the effect sizes were comparable to those observed when using the entire control sample (Table S7).
Discussion
We evaluated 71 SNPs associated with CHD risk within 65 risk regions in a large multi-ethnic sample of African Americans, Latinos, and Japanese Americans and found that a statistically significant proportion of SNPs exhibited consistent directions of effect beyond the 50% expected by chance. However, only a subset of 11, 8, and 1 of these SNPs were found to be nominally statistically significant in African Americans, Latinos, and Japanese Americans, respectively. Exploration of common genetic variation in these CHD-associated regions provided additional support for association at 10 regions, with different ethnic-specific or cross-ethnic leading SNPs. These replication results provide additional evidence for shared common genetic effects across ethnicities, with previous studies only replicating signals at 9p21 (21, 24–26), SORT1 (26), and ABO (26) in African Americans, and SORT1, CXCL12, and 9p21 in Latinos (29). This is the first report of BTN2A1, a region initially reported in Japanese, replicating in African Americans.
Japanese Americans had a higher GRS on average when compared to African Americans and Latinos. However, the GRS was more strongly associated with CHD in African Americans and Latinos compared to the Japanese Americans. We note that the genetic markers reported from previous discovery efforts are unlikely to be the functional alleles. The correlation between the index and functional SNPs may vary depending on the LD structure of each ancestral group, which may contribute to the difference in the ethnic-specific odds ratios. In addition to having limited statistical power to replicate associations with index SNPs within and across these populations, differences in LD may serve as an alternative explanation for the lack of replication. In an attempt to address such issues, we conducted regional association testing and constructed modified risk scores incorporating regional association results. When substituting the index SNPs with leading SNPs from the regional analyses, the differences in the modified risk score distributions and per-allele aggregate effects were only modified slightly, but differences were still noted, particularly between the Japanese and the other populations. The reasons for these differences are unclear. Our findings may reflect the severity of subclinical coronary atherosclerosis among Japanese participants in the MEC that is on average greater than the severity observed in Africans and Hispanics (44). Although our analyses are preliminary, we deem it unlikely that these known risk alleles are major contributors to race/ethnic differences in the incidence of CHD, as the incidence of CHD in Japanese is lower than that in the other two groups (45). Given that Japanese had a higher average GRS compared to other ethnic groups, but their population risk is lower, it is possible that functional variants within CHD susceptibility genes not included in our GRS disproportionally affect non-Japanese race/ethnic groups. Alternatively, environmental risk factors such as suboptimal diet and smoking may be less prevalent in Japanese and primarily responsible for the lower rates of CHD despite the higher genetic risk. It is difficult to directly compare the GRS distribution reported in this study to those in studies in European ancestry populations as the methods and number of selected SNPs vary (46–53). The vast majority of studies in European ancestry populations have observed statistically significant per allele relative risks of 1.02–1.12 and relative risks of 1.5–1.9 in comparing the highest versus lowest quintile or quartile of the GRS. Our findings in African Americans and Latinos are generally consistent with these reports albeit smaller effect sizes were noted, perhaps due to differences in LD between the index and functional SNPs.
Our study has a number of limitations. First, the information used to define CHD cases and controls was based on a combination of health care claims data as well as self-report on questionnaires. Some of the Japanese cases from Hawaii may have been missed due to the lack of CHDD records. Of the 1,089 Japanese participants whose CHDD records were available (in California), 426 CHD cases were identified, with 103 classified as cases based solely on CHDD records. Given the same ratio, about 338 Japanese CHD cases from Hawaii where CHDD was not available, may have been misclassified as controls. Assuming an equal distribution of genotypes in these missed cases compared to recognized cases, this misclassification would result in effects being biased towards the null and a reduced power to detect associations. Similar misclassification may apply as Medicare or CHDD data were not available for all controls. In the sensitivity analysis, limiting controls to those with claims data, fewer SNPs reached nominal statistical significance (0.05) however effect sizes were relatively comparable to those observed in the entire control sample. In an attempt to increase specificity when using Medicare and CHDD claims, we only included CHD cases identified from the primary and the first diagnosis codes, with individuals identified with CHD beyond the primary diagnosis excluded from being a case or a control. The validity of our case and control definitions is indirectly supported by the observed associations of case-control status with known risk factors and by the detection of more directionally consistent genetic associations than expected. Another limitation is the potential inclusion of inappropriately labeled CHD deaths, as CHD is often reported on death certificates when the cause of death is unclear. However, of the 1,005 CHD deaths, 718 also had prior claims data from CHDD or Medicare (Figure S2), suggesting a high consistency between mortality records and health care claims.
Another limitation of the study was the selection of CHD cases and controls among MEC participants conditional on three existing medical conditions. However, in sensitivity analyses limited to those without existing conditions, we observed robust consistency in terms of effect size and effect direction between this subset and the entire sample.
We present the largest replication study of established GWAS loci for CHD in Latinos and Japanese Americans conducted to date. However, our power was still limited to detect the originally reported effect sizes even in the combined multiethnic sample. We observed a higher GRS in Japanese Americans compared with African Americans or Latinos but the GRS was paradoxically not significantly associated with CHD risk in Japanese Americans despite observing strong associations in the other two groups. Substantially larger samples that include multiple racial/ethnic groups will help to identify the functional alleles in these regions and characterize their associations with CHD risk and contributions to CHD disparities among ethnically diverse populations.
Ethics Statement
The MEC study obtained informed consent from study participants and approval from the Health Science Review Board (HSIRB) at the University of Southern California and obtained IRB certification permitting data sharing in accordance with the NIH Policy for Sharing of Data Obtained in NIH Supported or Conducted Genome-Wide Association studies (GWAS).
Author Contributions
WK conducted the analysis and wrote the paper. KR contributed to the analysis. DC contributed to the analysis. VS contributed data for the manuscript. DS contributed to the analysis of the manuscript. LW contributed data for the manuscript. LL contributed data for the manuscript. TA contributed data and assisted with the writing of the manuscript. CH oversaw the project, contributed to the data, analysis and writing of the manuscript.
Funding
This study was funded by the National Cancer Institute Grant Number 2U01CA164973- 06 and NHGRI grant U01 HG007397
Conflict of Interest Statement
Author KR is currently employed by company Ancestry.com. All other authors declare no competing interests.
The reviewer IB and handling Editor declared their shared affiliation.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2018.00019/full#supplementary-material
Table S1. Genetic arrays by studies.
Table S2. Detailed demographic information on participants by ethnicity, gender, and existing conditions.
Table S3. Associations with CHD SNPs in the MEC populations.
Table S4. Results of regional replication by index or leading SNPs.
Table S5. Comparison of risk score distribution.
Table S6. Sensitivity analysis: comparing subpopulation with/without existing conditions.
Table S7. Sensitivity analysis: associations with CHD SNPs using controls with claims.
Figure S1. Statistical power for known CHD SNPs.
Figure S2.. CHD case sources venn diagram.
References
1. Selwyn AP, Braunwald E. Ischemic heart disease. In: Braunwald E, Fauci AS, Kasper DL, Hauser SL, Longo DL, Jameson JL et al. , editors. Harrison’s Principles of Internal Medicine. 15th ed. New York: McGraw-Hill (2001). p. 1399–410.
2. Zdravkovic S, Wienke A, Pedersen NL, Marenberg ME, Yashin AI, de Faire U. Heritability of death from coronary heart disease: a 36-year follow-up of 20 966 Swedish twins. J Intern Med (2002) 252(3):247–54. doi: 10.1046/j.1365-2796.2002.01029.x
3. Mcpherson R, Pertsemlidis A, Kavaslar N, Stewart A, Roberts R, Cox DR, et al. A common allele on chromosome 9 associated with coronary heart disease. Science (2007) 316(5830):1488–91. doi: 10.1126/science.1142447
4. Samani NJ, Erdmann J, Hall AS, Hengstenberg C, Mangino M, Mayer B, et al. Genomewide association analysis of coronary artery disease. N Engl J Med (2007) 357(5):443–53. doi: 10.1056/NEJMoa072366
5. Wild PS, Zeller T, Schillert A, Szymczak S, Sinning CR, Deiseroth A, et al. A genome-wide association study identifies LIPA as a susceptibility gene for coronary artery disease. Circ Cardiovasc Genet (2011) 4(4):403–12. doi: 10.1161/CIRCGENETICS.110.958728
6. Schunkert H, König IR, Kathiresan S, Reilly MP, Assimes TL, Holm H, et al. Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease. Nat Genet (2011) 43(4):333–8. doi: 10.1038/ng.784
7. Peden JF, Hopewell JC, Saleheen D, Chambers JC, Hager J, Coronary Artery Disease (C4D) Genetics Consortium. A genome-wide association study in Europeans and South Asians identifies five new loci for coronary artery disease. Nat Genet (2011) 43(4):339–44. doi: 10.1038/ng.782
8. Erdmann J, Willenborg C, Nahrstaedt J, Preuss M, König IR, Baumert J, et al. Genome-wide association study identifies a new locus for coronary artery disease on chromosome 10p11.23. Eur Heart J (2011) 32(2):158–68. doi: 10.1093/eurheartj/ehq405
9. Deloukas P, Kanoni S, Willenborg C, Farrall M, Assimes TL, Thompson JR, et al. Large-scale association analysis identifies new risk loci for coronary artery disease. Nat Genet (2013) 45(1):25–33. doi: 10.1038/ng.2480
10. Davies RW, Wells GA, Stewart AF, Erdmann J, Shah SH, Ferguson JF, et al. A genome-wide association study for coronary artery disease identifies a novel susceptibility locus in the major histocompatibility complex. Circ Cardiovasc Genet (2012) 5(2):217–25. doi: 10.1161/CIRCGENETICS.111.961243
11. Vargas JD, Manichaikul A, Wang XQ, Rich SS, Rotter JI, Post WS, et al. Common genetic variants and subclinical atherosclerosis: The Multi-Ethnic Study of Atherosclerosis (MESA). Atherosclerosis (2016) 245:230–6. doi: 10.1016/j.atherosclerosis.2015.11.034
12. Nikpay M, Goel A, Won HH, Hall LM, Willenborg C, Kanoni S, et al. A comprehensive 1,000 Genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Genet (2015) 47(10):1121–30. doi: 10.1038/ng.3396
13. Ozaki K, Sato H, Inoue K, Tsunoda T, Sakata Y, Mizuno H, et al. SNPs in BRAP associated with risk of myocardial infarction in Asian populations. Nat Genet (2009) 41(3):329–33. doi: 10.1038/ng.326
14. Aoki A, Ozaki K, Sato H, Takahashi A, Kubo M, Sakata Y, et al. SNPs on chromosome 5p15.3 associated with myocardial infarction in Japanese population. J Hum Genet (2011) 56(1):47–51. doi: 10.1038/jhg.2010.141
15. Wang F, Xu CQ, He Q, Cai JP, Li XC, Wang D, et al. Genome-wide association identifies a susceptibility locus for coronary artery disease in the Chinese Han population. Nat Genet (2011) 43(4):345–9. doi: 10.1038/ng.783
16. Lu X, Wang L, Chen S, He L, Yang X, Shi Y, et al. Genome-wide association study in Han Chinese identifies four new susceptibility loci for coronary artery disease. Nat Genet (2012) 44(8):890–4. doi: 10.1038/ng.2337
17. Takeuchi F, Yokota M, Yamamoto K, Nakashima E, Katsuya T, Asano H, et al. Genome-wide association study of coronary artery disease in the Japanese. Eur J Hum Genet (2012) 20(3):333–40. doi: 10.1038/ejhg.2011.184
18. Assimes TL, Lee IT, Juang JM, Guo X, Wang TD, Kim ET, et al. Genetics of coronary artery disease in Taiwan: a cardiometabochip study by the taichi consortium. PLoS ONE (2016) 11(3):e0138014. doi: 10.1371/journal.pone.0138014
19. Lee JY, Lee BS, Shin DJ, Woo Park K, Shin YA, Joong Kim K, et al. A genome-wide association study of a coronary artery disease risk variant. J Hum Genet (2013) 58(3):120–6. doi: 10.1038/jhg.2012.124
20. Shen GQ, Li L, Rao S, Abdullah KG, Ban JM, Lee BS, et al. Four SNPs on chromosome 9p21 in a South Korean population implicate a genetic locus that confers high cross-race risk for development of coronary artery disease. Arterioscler Thromb Vasc Biol (2008) 28(2):360–5. doi: 10.1161/ATVBAHA.107.157248
21. Assimes TL, Knowles JW, Basu A, Iribarren C, Southwick A, Tang H, et al. Susceptibility locus for clinical and subclinical coronary artery disease at chromosome 9p21 in the multi-ethnic ADVANCE study. Hum Mol Genet (2008) 17(15):2320–8. doi: 10.1093/hmg/ddn132
22. Barbalic M, Reiner AP, Wu C, Hixson JE, Franceschini N, Eaton CB, et al. Genome-wide association analysis of incident coronary heart disease (CHD) in African Americans: a short report. PLoS Genet (2011) 7(8):e1002199. doi: 10.1371/journal.pgen.1002199
23. Franceschini N, Carty C, Bůzková P, Reiner AP, Garrett T, Lin Y, et al. Association of genetic variants and incident coronary heart disease in multiethnic cohorts: the PAGE study. Circ Cardiovasc Genet (2011) 4(6):661–72. doi: 10.1161/CIRCGENETICS.111.960096
24. Kral BG, Mathias RA, Suktitipat B, Ruczinski I, Vaidya D, Yanek LR, et al. A common variant in the CDKN2B gene on chromosome 9p21 protects against coronary artery disease in Americans of African ancestry. J Hum Genet (2011) 56(3):224–9. doi: 10.1038/jhg.2010.171
25. Lettre G, Palmer CD, Young T, Ejebe KG, Allayee H, Benjamin EJ, et al. Genome-wide association study of coronary heart disease and its risk factors in 8,090 African Americans: the NHLBI CARe Project. PLoS Genet (2011) 7(2):e1001300. doi: 10.1371/journal.pgen.1001300
26. Franceschini N, Hu Y, Reiner AP, Buyske S, Nalls M, Yanek LR, et al. Prospective associations of coronary heart disease loci in African Americans using the MetaboChip: the PAGE study. PLoS ONE (2014) 9(12):e113203. doi: 10.1371/journal.pone.0113203
27. Manichaikul A, Wang XQ, Musani SK, Herrington DM, Post WS, Wilson JG, et al. Association of the lipoprotein receptor SCARB1 common missense variant rs4238001 with incident coronary heart disease. PLoS ONE (2015) 10(5):e0125497. doi: 10.1371/journal.pone.0125497
28. Zhang L, Buzkova P, Wassel CL, Roman MJ, North KE, Crawford DC, et al. Lack of associations of ten candidate coronary heart disease risk genetic variants and subclinical atherosclerosis in four US populations: the Population Architecture using Genomics and Epidemiology (PAGE) study. Atherosclerosis (2013) 228(2):390–9. doi: 10.1016/j.atherosclerosis.2013.02.038
29. Qi L, Campos H. Genetic predictors for cardiovascular disease in hispanics. Trends Cardiovasc Med (2011) 21(1):15–20. doi: 10.1016/j.tcm.2012.01.002
30. Han Y, Hazelett DJ, Wiklund F, Schumacher FR, Stram DO, Berndt SI, et al. Integration of multiethnic fine-mapping and genomic annotation to prioritize candidate functional SNPs at prostate cancer susceptibility regions. Hum Mol Genet (2015) 24(19):5603–18. doi: 10.1093/hmg/ddv269
31. Feng Y, Stram DO, Rhie SK, Millikan RC, Ambrosone CB, John EM, et al. A comprehensive examination of breast cancer risk loci in African American women. Hum Mol Genet (2014) 23(20):5518–26. doi: 10.1093/hmg/ddu252
32. Han Y, Signorello LB, Strom SS, Kittles RA, Rybicki BA, Stanford JL, et al. Generalizability of established prostate cancer risk variants in men of African ancestry. Int J Cancer (2015) 136(5):1210–7. doi: 10.1002/ijc.29066
33. Williams AL, Jacobs SB, Moreno-Macías H, Huerta-Chagoya A, Churchhouse C, Márquez-Luna Cet al, et al. Sequence variants in SLC16A11 are a common risk factor for type 2 diabetes in Mexico. Nature (2014) 506(7486):97–101. doi: 10.1038/nature12828
34. Siddiq A, Couch FJ, Chen GK, Lindström S, Eccles D, Millikan RC, et al. A meta-analysis of genome-wide association studies of breast cancer identifies two novel susceptibility loci at 6q14 and 20q11. Hum Mol Genet (2012) 21(24):5373–84. doi: 10.1093/hmg/dds381
35. Voight BF, Kang HM, Ding J, Palmer CD, Sidore C, Chines PS, et al. The metabochip, a custom genotyping array for genetic studies of metabolic, cardiovascular, and anthropometric traits. PLoS Genet (2012) 8(8):e1002793. doi: 10.1371/journal.pgen.1002793
36. Patterson N, Price AL, Reich D. Population structure and eigenanalysis. PLoS Genet (2006) 2(12):e190. doi: 10.1371/journal.pgen.0020190
37. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet (2006) 38(8):904–9. doi: 10.1038/ng1847
38. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics (2005) 21(2):263–5. doi: 10.1093/bioinformatics/bth457
39. Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, Handsaker RE, et al. An integrated map of genetic variation from 1,092 human genomes. Nature (2012) 491(7422):56–65. doi: 10.1038/nature11632
40. Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, et al. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics (2010) 26(18):2336–7. doi: 10.1093/bioinformatics/btq419
41. Zhao JH. Gap : genetic analysis package. J Stat Softw (2007) 23(8):1–18. doi: 10.18637/jss.v023.i08
42. Gauderman WJ, Morrison JM. (2006) QUANTO 1.1: A computer program for power and sample size calculations for genetic-epidemiology studies. http://hydra.usc.edu/gxe
43. Musunuru K, Strong A, Frank-Kamenetsky M, Lee NE, Ahfeldt T, Sachs KV, et al. From noncoding variant to phenotype via SORT1 at the 1p13 cholesterol locus. Nature (2010) 466(7307):714–9. doi: 10.1038/nature09266
44. Salfati E, Nandkeolyar S, Fortmann SP, Sidney S, Hlatky MA, Quertermous T, et al. Susceptibility Loci for clinical coronary artery disease and subclinical coronary atherosclerosis throughout the life-course. Circ Cardiovasc Genet (2015) 8(6):803–11. doi: 10.1161/CIRCGENETICS.114.001071
45. Holland AT, Wong EC, Lauderdale DS, Palaniappan LP. Spectrum of cardiovascular diseases in Asian-American racial/ethnic subgroups. Ann Epidemiol (2011) 21(8):608–14. doi: 10.1016/j.annepidem.2011.04.004
46. Vaarhorst AA, Lu Y, Heijmans BT, Dollé ME, Böhringer S, Putter H, et al. Literature-based genetic risk scores for coronary heart disease: the Cardiovascular Registry Maastricht (CAREMA) prospective cohort study. Circ Cardiovasc Genet (2012) 5(2):202–9. doi: 10.1161/CIRCGENETICS.111.960708
47. Krarup NT, Borglykke A, Allin KH, Sandholt CH, Justesen JM, Andersson EA, et al. A genetic risk score of 45 coronary artery disease risk variants associates with increased risk of myocardial infarction in 6041 Danish individuals. Atherosclerosis (2015) 240(2):305–10. doi: 10.1016/j.atherosclerosis.2015.03.022
48. Paynter NP, Chasman DI, Paré G, Buring JE, Cook NR, Miletich JP, et al. Association between a literature-based genetic risk score and cardiovascular events in women. JAMA (2010) 303(7):631–7. doi: 10.1001/jama.2010.119
49. Ripatti S, Tikkanen E, Orho-Melander M, Havulinna AS, Silander K, Sharma A, et al. A multilocus genetic risk score for coronary heart disease: case-control and prospective cohort analyses. Lancet (2010) 376(9750):1393–400. doi: 10.1016/S0140-6736(10)61267-6
50. Khera AV, Emdin CA, Drake I, Natarajan P, Bick AG, Cook NR, et al. Genetic risk, adherence to a healthy lifestyle, and coronary disease. N Engl J Med (2016) 375(24):2349–58. doi: 10.1056/NEJMoa1605086
51. Thanassoulis G, Peloso GM, Pencina MJ, Hoffmann U, Fox CS, Cupples LA, et al. A genetic risk score is associated with incident cardiovascular disease and coronary artery calcium: the Framingham heart study. Circ Cardiovasc Genet (2012) 5(1):113–21. doi: 10.1161/CIRCGENETICS.111.961342
52. Brautbar A, Pompeii LA, Dehghan A, Ngwa JS, Nambi V, Virani SS, et al. A genetic risk score based on direct associations with coronary heart disease improves coronary heart disease risk prediction in the Atherosclerosis Risk in Communities (ARIC), but not in the Rotterdam and Framingham offspring, studies. Atherosclerosis (2012) 223(2):421–6. doi: 10.1016/j.atherosclerosis.2012.05.035
Keywords: coronary heart disease, genome wide association study (GWAS), multi-ethnic, African Americans, Latino American, Japanese Americans, SORT1
Citation: Ke W, Rand KA, Conti DV, Setiawan VW, Stram DO, Wilkens L, Le Marchand L, Assimes TL and Haiman CA (2018). Evaluation of 71 Coronary Artery Disease Risk Variants in a Multiethnic Cohort. Front. Cardiovasc. Med. 5:19. doi: 10.3389/fcvm.2018.00019
Received: 06 November 2017; Accepted: 21 February 2018;
Published: 14 March 2018
Edited by:
Mete Civelek, University of Virginia, United StatesReviewed by:
Thorsten Kessler, Deutsches Herzzentrum München, GermanyCasey E Romanoski, University of Arizona, United States
Ingrid Braenne, University of Virginia, United States
Copyright © 2018 Ke, Rand, Conti, Setiawan, Stram, Wilkens, LeMarchand, Assimes and Haiman. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Christopher A. Haiman, Christopher.Haiman@med.usc.edu