Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Gender Differences in Genetic Risk Profiles for Cardiovascular Disease

  • Kaisa Silander ,

    kaisa.silander@ktl.fi

    Affiliations Department of Molecular Medicine, National Public Health Institute, Helsinki, Finland, FIMM, Institute for Molecular Medicine Finland, Helsinki, Finland

  • Mervi Alanne ,

    Contributed equally to this work with: Mervi Alanne, Kati Kristiansson

    Affiliations Department of Molecular Medicine, National Public Health Institute, Helsinki, Finland, FIMM, Institute for Molecular Medicine Finland, Helsinki, Finland

  • Kati Kristiansson ,

    Contributed equally to this work with: Mervi Alanne, Kati Kristiansson

    Affiliations Department of Molecular Medicine, National Public Health Institute, Helsinki, Finland, FIMM, Institute for Molecular Medicine Finland, Helsinki, Finland

  • Olli Saarela,

    Affiliation Department of Health Promotion and Chronic Disease Prevention, National Public Health Institute, Helsinki, Finland

  • Samuli Ripatti,

    Affiliations Department of Molecular Medicine, National Public Health Institute, Helsinki, Finland, FIMM, Institute for Molecular Medicine Finland, Helsinki, Finland, Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden

  • Kirsi Auro,

    Affiliations Department of Molecular Medicine, National Public Health Institute, Helsinki, Finland, FIMM, Institute for Molecular Medicine Finland, Helsinki, Finland

  • Juha Karvanen,

    Affiliation Department of Health Promotion and Chronic Disease Prevention, National Public Health Institute, Helsinki, Finland

  • Sangita Kulathinal,

    Affiliation Indic Society for Education and Development, Nashik, India

  • Matti Niemelä,

    Affiliation Department of Health Promotion and Chronic Disease Prevention, National Public Health Institute, Helsinki, Finland

  • Pekka Ellonen,

    Affiliations Department of Molecular Medicine, National Public Health Institute, Helsinki, Finland, FIMM, Institute for Molecular Medicine Finland, Helsinki, Finland

  • Erkki Vartiainen,

    Affiliation Department of Health Promotion and Chronic Disease Prevention, National Public Health Institute, Helsinki, Finland

  • Pekka Jousilahti,

    Affiliation Department of Health Promotion and Chronic Disease Prevention, National Public Health Institute, Helsinki, Finland

  • Janna Saarela,

    Affiliations Department of Molecular Medicine, National Public Health Institute, Helsinki, Finland, FIMM, Institute for Molecular Medicine Finland, Helsinki, Finland

  • Kari Kuulasmaa,

    Affiliation Department of Health Promotion and Chronic Disease Prevention, National Public Health Institute, Helsinki, Finland

  • Alun Evans,

    Affiliation The Queen's University of Belfast, Belfast, United Kingdom

  • Markus Perola,

    Affiliations Department of Molecular Medicine, National Public Health Institute, Helsinki, Finland, FIMM, Institute for Molecular Medicine Finland, Helsinki, Finland, Department of Medical Genetics, University of Helsinki, Helsinki, Finland

  • Veikko Salomaa,

    Affiliation Department of Health Promotion and Chronic Disease Prevention, National Public Health Institute, Helsinki, Finland

  •  [ ... ],
  • Leena Peltonen

    Affiliations Department of Molecular Medicine, National Public Health Institute, Helsinki, Finland, FIMM, Institute for Molecular Medicine Finland, Helsinki, Finland, Department of Medical Genetics, University of Helsinki, Helsinki, Finland, Department of Human Genetics, Wellcome Trust Sanger Institute, Cambridge, United Kingdom, The Broad Institute of MIT and Harvard, Boston, Massachusetts, United States of America

  • [ view all ]
  • [ view less ]

Abstract

Background

Cardiovascular disease (CVD) incidence, complications and burden differ markedly between women and men. Although there is variation in the distribution of lifestyle factors between the genders, they do not fully explain the differences in CVD incidence and suggest the existence of gender-specific genetic risk factors. We aimed to estimate whether the genetic risk profiles of coronary heart disease (CHD), ischemic stroke and the composite end-point of CVD differ between the genders.

Methodology/Principal Findings

We studied in two Finnish population cohorts, using the case-cohort design the association between common variation in 46 candidate genes and CHD, ischemic stroke, CVD, and CVD-related quantitative risk factors. We analyzed men and women jointly and also conducted genotype-gender interaction analysis. Several allelic variants conferred disease risk for men and women jointly, including rs1801020 in coagulation factor XII (HR = 1.31 (1.08–1.60) for CVD, uncorrected p = 0.006 multiplicative model). Variant rs11673407 in the fucosyltransferase 3 gene was strongly associated with waist/hip ratio (uncorrected p = 0.00005) in joint analysis. In interaction analysis we found statistical evidence of variant-gender interaction conferring risk of CHD and CVD: rs3742264 in the carboxypeptidase B2 gene, p(interaction) = 0.009 for CHD, and rs2774279 in the upstream stimulatory factor 1 gene, p(interaction) = 0.007 for CHD and CVD, showed strong association in women but not in men, while rs2069840 in interleukin 6 gene, p(interaction) = 0.004 for CVD, showed strong association in men but not in women (uncorrected p-values). Also, two variants in the selenoprotein S gene conferred risk for ischemic stroke in women, p(interaction) = 0.003 and 0.007. Importantly, we identified a larger number of gender-specific effects for women than for men.

Conclusions/Significance

A false discovery rate analysis suggests that we may expect half of the reported findings for combined gender analysis to be true positives, while at least third of the reported genotype-gender interaction results are true positives. The asymmetry in positive findings between the genders could imply that genetic risk loci for CVD are more readily detectable in women, while for men they are more confounded by environmental/lifestyle risk factors. The possible differences in genetic risk profiles between the genders should be addressed in more detail in genetic studies of CVD, and more focus on female CVD risk is also warranted in genome-wide association studies.

Introduction

According to world statistics for 2006, cardiovascular diseases (CVD) are responsible for 30% of all deaths globally, and are the leading cause of death amongst non-communicable diseases. Cardiovascular diseases are also responsible for 10% of the global burden of disease [1]. Differences in CVD incidence, complications and burden exist between men and women. Women are afflicted with cardiovascular disease at an older age than men, and many risk variables for coronary heart disease (CHD) and stroke have different distributions in men and women [2][7]. However, the differences in lifestyle factors do not fully explain the differences in CVD incidence between the genders [2]. Genetic factors also contribute to CHD and stroke susceptibility [8][10]. A recent large population-based prospective study suggested that heritability of ischemic stroke was greater in women than men [11]. Some of the traditional CVD risk factors also have high heritability [10], some of which show gender differences [12]. A large scale study of CVD traits in a Sardinian population showed that for several traits in which heritability estimates differed by gender, for example weight and hip circumference, the heritability was larger among women [13]. The evidence for gender differences in trait heritabilities implies possible gender-gene interaction in the etiology of these traits [12].

The effect of genetic variables on CHD and ischemic stroke has been studied for several decades, yet there are only a few consistent risk factors identified to date [10], [12], [14][19]. These genetic studies include few large scale candidate gene studies, as well as numerous smaller studies, and very recently several genome-wide association studies. Most of the large scale candidate gene studies published so far on CHD or stroke have performed combined analyses of both genders, using gender as a covariate [20][25]. In a Japanese case-control study of myocardial infarction, men and women were analyzed separately, and the significant results obtained for men and women were for different variants [26], indicating different genetic risk factors. In a large-scale genetic association study of the metabolic syndrome among CHD patients, McCarthy and colleagues identified several variants which displayed significant genotype-gender interaction [27]. In recent genome-wide association studies of CHD [14][17], [19] and ischemic stroke [28], the association results were reported for the combined study sample of both genders.

We estimated the effect of genetic variation on CHD, ischemic stroke and the composite end-point of CVD in two prospectively followed population cohorts. Our study had a case-cohort design on the FINRISK-92 and -97 cohorts participating in the MORGAM Project [29]. We selected 46 genes for study as putatively involved in cardiovascular pathobiology, based on their function, previous association with cardiovascular disease, and/or relevant phenotype in animal models. These genes represent a selected array of pathways, including lipid and energy metabolism, inflammation, coagulation, and thrombosis. We assessed the risk associated with common variation in each gene and CHD, ischemic stroke, and CVD while the cohort setting allowed us to control for classic CVD risk factors. We also assessed whether the variants affect relevant quantitative traits that are related to CVD risk: lipid and C-reactive protein (CRP) levels, blood pressure, body mass index (BMI), and waist/hip ratio (WHR). Our previous analysis of candidate genes like upstream stimulatory factor 1 (USF1) and Selenoprotein S (SEPS1, SELS, or SELENOS) mainly showed genetic effects in women [30], [31]. In this study, we therefore proceeded with a formal genotype-gender interaction analysis for all variants, and show that for several of the associated variants, there is evidence for statistical interaction between gender and genotype.

Materials and Methods

FINRISK cohort description

FINRISK surveys are carried out every 5 years to assess the prevalence and risk factors of CVD in Finland [32]. Baseline information on all randomly sampled individuals includes anthropometric measurements, serum lipids, blood pressure and questionnaire data on CVD risk factors. Information on fatal and non-fatal coronary and stroke events and all-cause mortality during the follow-up period is obtained from national registers. We utilized the FINRISK-92 cohort (n = 5999) and FINRISK-97 cohort (n = 8141), which have been followed up for 10 and 7 years, respectively. On these large cohorts, we conducted a case-cohort study, as previously described in detail [29], [31], [33][35]. The cohorts constituted respondents to surveys of independent random samples of the same geographically defined population. The resulting few overlaps were identified on the basis of personal ID codes, unique to every resident of Finland, and removed from the FINRISK-97 case-cohort set to ensure there was no overlap between the sets used for the analyses.

We initially studied the FINRISK-92 case-cohort set, which consisted of a total of 190 incident CHD cases, 66 incident ischemic stroke cases, 219 individuals with a history of either CHD or stroke event, 276 individuals who died during the follow-up, and a random sample (sub-cohort) of 398 individuals from the cohort. We also analyzed a second case-cohort set selected from the FINRISK-97 cohort, for genes associated with risk for CHD, ischemic stroke, the composite end-point of CVD or all-cause mortality, or strongly associated with quantitative traits in the FINRISK-92 case-cohort set. This sample included 210 incident CHD cases, 84 incident stroke cases, 436 individuals with a history of either CHD or stroke event, 352 individuals who died during the follow-up, and 407 sub-cohort individuals. The sub-cohort was a sex- and geographic-region stratified random sample, drawn from each of the original cohorts with unequal sampling probabilities so that the age distribution was similar to the cases. The selection procedure for the cases and the sub-cohort, and the exact diagnostic criteria used for CHD and ischemic stroke have been described in detail previously [29], [31], [33][35]. The case-cohort sets included in this study are described in Tables 1 and 2. All participants gave informed consent. In 1992 it was not yet customary to ask for a written consent, thus only oral informed consent exists for that survey. In 1997 a written informed consent was obtained from all survey participants. The law about the National Public Health Institute of Finland gives the Institute a possibility to also use the samples from the 1992 survey for public health research. The study was approved by the Ethics Committee of the National Public Health Institute of Finland and conformed to the principles expressed in the Declaration of Helsinki.

thumbnail
Table 1. Number of individuals in each of the case status categories in the FINRISK-92 and FINRISK-97 case-cohort study.

https://doi.org/10.1371/journal.pone.0003615.t001

thumbnail
Table 2. Baseline characteristics according to case status, mean (std).

https://doi.org/10.1371/journal.pone.0003615.t002

Quality control of DNA samples

We implemented several quality control measures to minimize errors associated with DNA sample handling and DNA quality, and excluded a total of 19 samples chosen as cases or in the sub-cohort. These 19 individuals are not included in Table 1. A gender-specific PCR test identified a total of 9 samples (0.4%) that had a different gender than expected, and they were subsequently excluded from the study. We also verified that the DNA sample was of good quality by testing five highly polymorphic microsatellite markers for each sample. In these analyses, one sample was found to be contaminated and was excluded. DNA samples with low DNA yield (<7.5 µg of genomic DNA) as measured by fluorescent label PicoGreen (Invitrogen, Carlsbad, CA, USA) were subjected to whole genome amplification before genotyping, followed by additional quality control checks [36]. A total of five samples were excluded due to biased whole genome amplification, and a further 4 samples were excluded due to extremely low quantities of DNA which was insufficient for whole genome amplification.

Variant selection

For each gene, we aimed to genotype a set of variants that would capture the common variation present in the gene, as well as variants that have been previously associated with CVD or related traits. For the majority of the genes, haplotype-tagging single nucleotide polymorphism (SNP) variants were selected from the SeattleSNPs database (http://pga.gs.washington.edu/). The SeattleSNPs project has resequenced the genes using 24 Centre d'Etude du Polymorphisme Humain DNA samples, and tag SNPs have been selected using LDSelect, an algorithm that is based on the linkage disequilibrium (LD) statistic r2 [37]. We selected tag SNPs from each multi-SNP bin with a frequency >10%. For genes that were not included in the SeattleSNPs sequencing project, we selected variants from public databases (Celera, dbSNP), at approximately 5 kb distance from one another, giving priority to variants with known frequency information. Once HapMap phase I data were available, we selected additional variants to better capture the common variation in these genes. More detailed information about gene cladistics, sequence and haplotype structure information was available for apolipoprotein E (APOE), lactase (LCT), and lipin 1 (LPIN1)-genes, and here variant selection was based on previously published sequencing and haplotype analysis [38][41]. A full list of the variants selected for study and successfully genotyped (see below) is provided in Table S1.

Variant genotyping

Variant genotyping was done using several genotyping platforms (Table S1). Approximately 5.5% of the genotypes were created with an in-house developed method of allele-specific primer extension on microarrays, as previously described [36]. Approximately 93.0% of the genotypes were produced with the MassARRAY System (Sequenom, San Diego, CA, USA), either with the homogeneous Mass Extension (hME) reaction or iPLEX reaction, using the protocols recommended by the manufacturer with these modifications: hME reactions were carried out with 5–7.5 ng of DNA and for the majority of the variants, the hME extension reaction was run using TERMIPol DNA polymerase (Solis Biodyne OÜ, Tartu, Estonia) [42] instead of ThermoSequenase (GE Healthcare Life Sciences, Chalfont St. Giles, UK). The two APOE variants that define the epsilon genotypes (rs429358 and rs7412) were genotyped on the MassARRAY with a modified protocol as previously described [43] (full protocol available from authors upon request). Three of the variants were genotyped with other platforms: rs4340 was genotyped by a PCR assay followed by separation on 2% agarose gel with ethidium bromide staining and rs28665122 and rs3216183 were genotyped with TaqMan (Applied Biosystems, Foster City, CA) [30], [44]. For 100 samples where inadequate amount of genomic DNA was available, the DNA was amplified with GenomiPhi DNA amplification kit (GE Healthcare Life Sciences), as previously described [36].

Before genotyping the FINRISK case-cohort samples we genotyped all variants on 60 anonymous Finnish trio samples and 180 unrelated control samples. The FINRISK samples were genotyped in plates containing 2% negative control samples, 2% known duplicate samples, and 5% blind duplicate samples to allow assessment of genotyping quality. The disease status of each individual genotyped was unknown to the genotyping laboratory and samples from cases and sub-cohort individuals were distributed on the plates independently of the disease status. All genotypes were manually reviewed for various quality control aspects as previously described [36], [42], [45]. The genotyping success rate for each variant included in the analysis was >90%, with an average genotyping success of 95.3%. Among the 27,522 successful blind duplicate genotypic pairs, we detected 37 genotypic inconsistencies (99.87% concordance between genotypes). All variants included in analyses were in Hardy-Weinberg equilibrium (HWE) in the sub-cohort sample (p>0.01). A single Mendelian error was identified for 3 variants among the 60 trio samples (rs1926446, rs3212478, and rs1081106). However, since the genotypes for these variants were in HWE and no errors were detected among known and blind duplicates, these variants were included in the analysis.

Statistical analysis

Genotype frequencies in sub-cohort individuals were tested for deviation from HWE using Pearson's chi-square test statistics with 1 degree of freedom for bi-allelic variants and 3 for three-allelic variants, applying a threshold of p<0.01. For variants in which one of the genotype groups had less than 5 individuals, HWE was calculated using an exact test. Allele segregation within trio families was analyzed with the PedCheck program [46]. Pair-wise LD between the variants in each gene, haplotype frequencies, and haplotype tags were assessed with Haploview software version 3.32 [47]. For variants in high LD with each other (r2>0.95), only one of the results is shown.

Time-to-event analysis was used to assess whether any of the tested allelic variants have effect on the incidence of CHD, ischemic stroke, or CVD. The effects under recessive, dominant and multiplicative models of individual variants were tested using the proportional hazards regression model where the case-cohort design was taken into account by applying a modification of the Prentice weighting [48], with the non-case sub-cohort members and sub-cohort cases before events weighted with the inverses of their individual inclusion probabilities to account for the over-sampling of cases[34]. Estimation of model parameters and standard errors was carried out in R statistical environment, using the coxph function of the package survival and its robust variance estimator. We adjusted for classic CVD risk factors: smoking, high density lipoprotein-cholesterol (HDL-C), non-HDL-cholesterol, history of diabetes, BMI, and hypertension, as well as geographic region (western Finland, northern Finland, and eastern Finland), and cohort (and gender for combined analysis in women and men). Age was used in the models as the time scale. We fitted two types of models. In the first model, men and women from both cohorts were analyzed jointly, as described above. In the second model, we carried out a test for genotype-gender interaction, defined as a departure from multiplicative, dominant or recessive model, using similar regression models and testing the null hypothesis of equality of genotype effect parameters between men and women. We report results in which the variant genotype specific p-value is ≤0.01 for either men or women. We verified that these results do not stem from a single cohort by testing the null hypothesis of equality of genotype effect parameters between FINRISK-92 and FINRISK-97 cohorts, using a similar regression model. For variants that conferred a risk at p<0.05 for CHD, we also studied the association in prevalent CHD cases (documented or self-reported myocardial infarction or unstable angina pectoris at baseline), using healthy sub-cohort subjects as controls. The analysis of prevalent cases was carried out using logistic regression, again with inverse sampling probability weighting, and using age, cohort and geographic region, and gender as covariates for the combined analysis of men and women. Analysis of haplotype effects was done for two variants of the F12 gene that were not in very high LD with each other and were both associated at p<0.01 with CHD and CVD. Haplotype analysis was done with an additive model, in which the common haplotype (containing the ‘non-risk’ alleles) was used as reference, and modeling an additive effect for the other haplotypes, in a weighted Cox proportional hazards model, applying the same weighting scheme and covariates that were used for single variant analysis, and using the PHREG procedure implemented in SAS version 9.1.3 SP4. Haplotype uncertainty was taken into account using multiple imputations, where a sample of haplotypes was obtained using Phase 2.1.1 software and the analysis was repeated for each sampled haplotype pair.

Additionally, we tested whether allelic variants were associated with quantitative traits measured at baseline in sub-cohort individuals without a history of CVD. The lipid variables studied were: serum total cholesterol, HDL-C, triglycerides, and low density lipoprotein-cholesterol (LDL-C). LDL-C was calculated from measured values of total cholesterol, HDL-C and triglycerides using Friedewald's formula and excluding individuals with triglyceride value >4.0 mmol/l. Additional variables studied were mean blood pressure (average of systolic and diastolic blood pressure, each value based on two subsequent measurements), high sensitivity CRP, BMI, and WHR. Association of the variants with baseline measurements was tested using standard linear regression, employing additive, dominant, and recessive models, while adjusting for cohort, age, geographic region, and gender. Tests for genotype-gender interaction, defined as a departure from additive, dominant or recessive model, were carried out using similar regression models and testing the null hypothesis of equality of genotype effect parameters between men and women. Individuals using lipid lowering medication were excluded from the analyses of lipid variables, and individuals using drugs for hypertension were excluded from the analysis of blood pressure. We used logarithmic transformation for CRP and triglycerides. We verified that the results reported do not stem from a single cohort by testing the null hypothesis of equality of genotype effect parameters between FINRISK-92 and FINRISK-97 cohorts, using a similar regression model.

For genes in which two or more variants (not in perfect LD) were associated at p<0.01 with a given quantitative trait, we also performed haplotype analysis to discern which allelic haplotype might be contributing to variation in the trait. Haplotype tagging variants were identified with the Haploview software version 3.32 using default settings. Analyses with the haplotype-tagging variants were performed with the haplo.stats package of the R statistical software [49], using the function haplo.glm with an additive model, and adjusting for age, cohort, geographic region and gender. The haplo.glm function estimates haplotype frequencies with the EM algorithm and calculates for each haplotype linear regression coefficient and p-value, comparing each haplotype to a base haplotype, defined as the most common haplotype. Rare haplotypes (frequency <0.05) were combined with the base haplotype for this analysis. The global p-value for haplotype effect coefficients was calculated for the null hypothesis of no effect for any haplotype.

For the initial analyses of the FINRISK-92 case-cohort alone, time-to-event analyses and quantitative trait analyses were done as previously described [30], [31], [33], analyzing women and men both separately and together. We did not perform formal gender-genotype interaction analysis or haplotype analysis at this stage.

In reporting the findings, we used a cut-off value of 0.01 for the p-values and reported uncorrected p-values. The cut-off value of 0.01 corresponds to posterior odds 6:1 of a finding being a true signal when we expect to see two signals among the 27 independent genes and our power is 70% (see The Wellcome Trust Case-control Consortium's 2007 paper for details) [19]. The effect of multiple testing was addressed with standard Q-Q-plots for the individual test statistics and with false discovery rate (FDR) analysis [50], [51]. The tail-area FDR statistic for a group of tests can be interpreted as the expected proportion of null results given the observed test statistics. The analysis was carried out using the R package “fdrtool” [52]. The method used for power simulations is described in more detail elsewhere [34]. The reported results are for both cohorts combined, for tests of the null hypothesis of no genotype effects (or no genotype-gender interaction) at 1% significance level. While simulating genotype-gender interaction we assumed no genotype effects for men while varying the effect for women.

Results

Study outline

The case-cohort sets from the FINRISK-92 (10 year follow up, 57,858 person-years) and FINRISK-97 (7 year follow up, 54,577 person-years) population cohorts [31] are presented in Tables 1 and 2. The list of genes and the number of variants successfully genotyped for each gene are presented in Table 3, and detailed information on all variants is presented in Table S1. In addition to known CVD candidate genes, we explored the effect of variation in the LCT gene on CVD risk and CVD related quantitative traits, because of previous findings of reduced triglyceride and cholesterol values in individuals with lactose malabsorption [53], [54]. We also studied one novel gene, apolipoprotein B mRNA editing enzyme (APOBEC2), which is located directly under a linkage peak (lod score of 4.44) for total cholesterol in our linkage study of 5775 individuals from twin families from the GenomEUtwin (www.genomeutwin.org). Individual results of the analysis for several of the genes have already been published: USF1, thrombomodulin (THBD), SEPS1, coagulation factor V (F5), protein C (PROC), and intercellular adhesion molecule 1 (ICAM1) [30], [31], [33], [44]. We include these genes here to provide a more complete picture of the observed difference in genetic susceptibility between men and women, and because formal genotype-gender interaction analysis was not reported for any of the genes in our previous publications.

The study outline is presented in Figure 1. Initially, we studied the 46 genes in the FINRISK-92 case-cohort set. We selected for further study in the FINRISK-97 sample 27 genes in which one or more variants showed an association with CHD, ischemic stroke, CVD, total mortality, or any of the quantitative traits in the FINRISK-92 cohort, either in women or men separately, or in combined analyses. The selection criterion was 60% FDR. A total of 172 variants were thus typed also in the FINRISK-97 case-cohort samples, as indicated in Table 3 and Table S1, and analyzed using the combined FINRISK-92 and FINRISK-97 case-cohort sets. Power simulations are presented in Figures S1 and S2. For time-to-event analysis (Figure S1), the combination of the two cohorts has a 88% power to detect a dominant gene main effect on CVD risk of 1.8 in men at p = 0.01, a 39% power to detect a similar effect in women, and 96% power to detect this effect size when analyzing women and men together, given a risk allele frequency of 0.2 assuming a proportional hazards model. For a higher allele frequency the power is somewhat higher. For gene-gender interaction analysis, our study sample has power to detect only large differences in risk effects at p = 0.01, for example 38% power to detect a difference of HR = 1.0 versus HR = 1.8 for allele frequency of 0.4. For quantitative traits (Figure S2), combining both cohorts provides a power of 75% for detecting a 0.3 standard deviation difference at allele frequency of 0.2 in men at p = 0.01, while the power is much lower for the smaller study sample of women. For gene-gender interaction analyses the power is >85% only for large differences in the effects, for example no effect in men and a coefficient of 0.6 in women.

Time-to-event analysis results

Analysis of both genders jointly.

Time-to-event analysis was used to assess the association between variants and CHD, ischemic stroke and the composite end point of CVD. Results with p≤0.01 from combined analysis of both cohorts and both genders are shown in Table 4. The estimated FDR for the set of all association tests (including tests for quantitative traits) with p≤0.01 is 53%. These analyses identified variants in angiotensin II receptor type 1 (AGTR1), APOE, carboxypeptidase B2 (CPB2), and coagulation factor XII (F12) as conferring risk of CHD. The two variants of the F12 gene also conferred risk of CVD, as did one variant of fibrinogen alpha chain (FGA) gene. Haplotype analysis for the two F12 variants, rs4976691 and rs1801020, in which carriers of the specific ‘risk’ haplotypes (CA, CG, or GA for rs4976691 and rs1801020, respectively) were compared to individuals homozygous for the non-risk haplotype GG did not reveal stronger association with CHD or CVD than analysis of single variants. For ischemic stroke, only one SEPS1 variant, rs7178239, was associated at p≤0.01 in the combined analysis of both genders, but only the women contributed to this effect (see below). The most consistent result was for CHD association with the F12 variant rs1801020 (men and women combined, p = 0.005 for additive model), which also conferred risk at the p<0.05 level for CHD in both women and men when analyzed separately. The rest of the variants showed association at p<0.05 level in only one gender. We tested whether the results were driven by only one of the cohorts by assessing genotype-cohort interaction, and observed no interaction at p<0.05, suggesting that the results are similar in both cohorts. Variant rs440446 of APOE showed association at p<0.05 also in both cohorts separately, while the rest of the variants showed association at p<0.05 in one cohort only, though a similar trend was observed in the other cohort.

thumbnail
Table 4. Results with p≤0.01for coronary heart disease and cardiovascular disease for the variants studied, analyzing men and women together.

https://doi.org/10.1371/journal.pone.0003615.t004

Gender-genotype interactions.

We performed gender-genotype interaction analysis to identify variants that showed different genetic effects in women and men. This test is sensitive to both effect direction and effect size. The variants that gave interaction p-value≤0.01 and were associated with CHD, ischemic stroke or the composite end-point of CVD at p≤0.01 in either women or men in combined analysis of both cohorts are presented in Figure 2 and Table S2. The estimated FDR for the set of all interaction tests with p≤0.01 is 70%, but by using the additional criteria of association p-value≤0.01 in at least one of the genders, the actual FDR is likely to be smaller. The gender-genotype interaction analysis supports our previous findings for USF1 and SEPS1 variants in which the disease risk was limited to women [30], [31], providing a gender-genotype interaction p-values<0.01 for the USF1 variant rs2774279 and for two SEPS1 variants, rs4965814 and rs9874. For the USF1 variant rs2774279, the results were also at p<0.05 for women in each cohort separately. Furthermore, for rs2774279 we also found evidence for association when analyzing prevalent female CHD cases in both cohorts combined (odds ratio of 1.58, 95% CI 1.04–2.40, p = 0.03). We identified variants in additional genes which showed gender-genotype interaction: CPB2 and coagulation factor XIII, A1 polypeptide (F13A1) conferred gender-specific risk in women for CHD, another variant in CPB2 conferred risk for CVD, and F5 for ischemic stroke; and for men, interleukin 6 (IL6) for CVD. The data obtained with F5 variant rs970741 is based on relatively small groups, with only 12 women incident stroke cases carrying the protective allele, and the result should be interpreted with caution. Genotype-cohort interaction analysis showed that none of the gender-specific results emerge from a strong effect in only one of the cohorts but rather both cohorts contribute to the result. For purpose of future meta-analyses, we provide data for all variants analyzed in both cohorts showing genotype-specific hazard ratios for men and women separately and number of individuals and person years in each genotype group (Tables S3a–c).

thumbnail
Figure 2. Gender-specific association between variants and coronary heart disease, ischemic stroke, and cardiovascular disease.

Results for gender-genotype interaction at p<0.05, and association in either women or men were at p≤0.01 (uncorrected p-values). Allele information: allele 1/allele 2, the minor allele is underlined. Multiplicative model: 11>12>22, dominant model: 11+12 vs 22, recessive model: 11 vs 12+22. Variants showing high pair-wise LD: CPB2 rs3581419 and rs3742264 (r2 = 0.827), SEPS1 rs496581 and rs7178239 (r2>0.7), and SEPS1 rs9874 and rs7178239 (r2>0.7). Detailed information is found in Table S2.

https://doi.org/10.1371/journal.pone.0003615.g002

We tested which of the variants conferring a CHD risk at p <0.05 were also associated with CHD in the prevalent cases. In addition to USF1 variant rs2774279, also the T allele of variant rs2073658 of USF1 conferred risk in both incident and prevalent female cases (HR = 1.62, 95% CI 1.04–2.52, p = 0.03 for incident cases, and odds ratio = 1.87, 95% CI 1.26–2.76, p = 0.002 for prevalent cases, additive model, T risk allele). A variant in the APOBEC2 gene, rs2395754, was associated with CHD in both prevalent and incident male cases (HR = 1.45, 95% CI 1.04–2.02, p = 0.03 for incident cases, and odds ratio = 1.43, 95% CI 1.06–1.94, p = 0.02 for prevalent cases, C allele homozygotes compared to T allele carriers).

Quantitative trait analysis results

We tested whether any of the 172 variants was associated with the CVD-related quantitative traits: total cholesterol, HDL-C, LDL-C, triglycerides, CRP, BMI, WHR and mean blood pressure, analyzing the sub-cohort individuals without a history of CVD at baseline examination. The results showing association in the combined data analyses of both genders at significance level of p<0.01 are shown in Table 5. The estimated FDR for the set of all association tests (including tests for time-to-event responses) with p≤0.01 is 53%. We identified 3 variants displaying effect differences between the cohorts using genotype-cohort interaction analysis (interaction p-value <0.05), and they were removed. APOE variant rs440446, conferring risk for CHD in time-to-event analysis (Table 4), was associated with triglyceride values, and FGA variant rs2070018 was associated with mean blood pressure, with heterozygotes having the highest blood pressure values. None of the other variants associated with CHD, ischemic stroke, or CVD at p≤0.01 in women and men combined, was associated at p<0.01 with the quantitative traits tested here. However, we identified several interesting associations with each of the traits studied, as discussed below.

thumbnail
Table 5. Results with p<0.01for associations between variants and quantitative traits as measured at baseline examination in sub-cohort subjects free of CVD at baseline, women and men combined.

https://doi.org/10.1371/journal.pone.0003615.t005

The strongest association identified for quantitative traits in the combined analysis of women and men was for fucosyltransferase 3 (FUT3) variant rs11673407 and WHR. For men the additive model gave a p-value = 0.00006; for women the association was weaker, but in the same direction (p = 0.07). Haplotype analysis for WHR in men using the FUT3 variants rs874232, rs778986, and rs11673407 identified haplotype CAG as the only one associated with WHR, compared to base haplotype TAA (p = 0.00008) (Table S4a), suggesting that the true causal variant is not one of these 3 variants. Another strong association was found for a rare synonymous CRP variant, rs1800947, and CRP levels in men (p = 0.0001, recessive model).

The LCT variants were associated with total cholesterol and LDL-C in the combined data: The lactase non-persistence genotype (defined as minor allele homozygotes for variant rs4988235) was associated with higher cholesterol values. Similarly to FUT3 variant, the association was stronger for men (for total cholesterol, p = 0.003 and p = 0.005 for variants rs4988235 and rs6719488, respectively, and for LDL-C p = 0.002, and p = 0.0005 for variants rs4988235 and rs6719488, respectively), and in females the association was weaker but in the same direction. Haplotype analysis using the 3 haplotype-tagging variants rs2304371, rs6719488, and rs4988235 for men implied that haplotype GGG, tagged by the G allele of variant rs2304371 was the one associated with both traits, p = 0.003 for total cholesterol and p = 0.005 for LDL-C (compared to base haplotype ATA) (Table S4b). Sub-cohort men homozygotes for the G allele of rs2304371 have the highest LDL-C values, 4.02 mmol/l (n = 14), compared to 3.74 for GA genotype (n = 109) and 3.55 for AA genotype (n = 243), p = 0.014 for the additive model. Variants rs6719488 and rs2304371 are located in the LCT gene itself, while the lactase non-persistence variant is located at 14 kb distance upstream of the LCT gene. The LCT locus on chromosome 2q21.3 is known for being strongly selected during human evolution, with the lactase persistence allele varying in frequency in different populations and even between geographic regions [55]. We observed no differences in allele frequencies of the lactase persistence genotype in the geographic regions studied here (G allele frequency 0.46 in Western Finland and 0.44 in Eastern Finland).

Variants that showed different effects on CVD-related quantitative traits in women and men are shown in Table 6, using an interaction p-value cut off ≤0.01 and an association cut off p<0.01 in either women or men in combined analysis of both cohorts. The estimated FDR for the set of all interaction tests with p≤0.01 is 70%, but the additional criteria of association p-value<0.01 in at least one of the genders makes the actual FDR smaller than the upper limit of 70%. As for the disease risk, also here variants in different genes were associated with the traits in women and men. In women, variants in the fibrinogen genes (FGA and FGG) were associated with HDL-C. Interestingly, none of the genes that are in lipid pathways were associated with lipid variables in women at p<0.01. For weight-related variables, variants that showed gender-specific effect were identified only in women. USF1 variant rs2774279, which was associated with CHD and CVD risk, was also associated with BMI in women, though risk allele carriers had lower BMI. Women with the risk allele also had lower values of CRP. Three variants in ICAM1 gene associated with WHR in women. Haplotype analysis did not reveal any ICAM1 haplotypes associated more strongly with the trait than single alleles. The largest number of gender-genotype interactions was identified for CRP levels in females.

thumbnail
Table 6. Results with p<0.01 for associations between variant and quantitative traits as measured at baseline examination, for variants showing evidence of gender-genotype interaction (p≤0.01), women and men of the sub-cohort free of CVD at baseline.

https://doi.org/10.1371/journal.pone.0003615.t006

For men, the APOBEC2 variant rs2395754, which associated with CHD in both incident and prevalent cases, was also associated with cholesterol variables. Men carrying the risk allele had higher levels of LDL-C, p = 0.001. In men also a variant in the serpin peptidase inhibitor, clade E member 1 gene was associated with mean blood pressure. In addition to these findings, very few male-specific results at p<0.01 were identified, as shown in Table 6. The strongest associations with lipids for men were for variants that also showed the same trend in women, as discussed above.

Discussion

The hormonal environment as well as tissue specific gene expression is known to differ significantly between the genders in vertebrates. For many human diseases, gender-dependent differences in the progression and extent of disease have been explained by sex hormones. These hormones may differentially affect gene expression in somatic tissues, thus leading to the gender specific susceptibility to disease [56]. Also for cardiovascular disease, critical determinants of gender differences are sex steroid hormones and their receptors [57]. They interact with and activate, together with other proteins, genes that are possibly involved in CVD pathogenesis in the endothelial and smooth muscle cells [2], [57]. Sex steroid hormones are also expressed in the liver and regulate lipid levels, mostly through hepatic effects on lipoprotein metabolism [57].

Although women and men differ in various aspects related to CHD and ischemic stroke [2][4], [6], [7], the difference in genetic effects on disease and its risk factors between women and men remains largely unexplored territory [12]. Recent genome-wide association studies also do not address this issue [14][17], [19], [58]. In this candidate gene study we explored the genetic risk profiles for CHD, ischemic stroke and the composite end point of CVD in men and women, as well as the effect of the specific genetic variants on CVD-related quantitative risk factors. Our case-cohort study was based on two prospective cohorts from the relatively homogeneous Finnish population, the sub-cohort representing a random subsample of the original cohort. Detailed information on CVD risk factors recorded before the occurrence of CVD events allowed us to control for confounding factors, such as smoking, lipid levels, blood pressure and obesity, while the inclusion of two separate cohorts allowed for the verification of results. We identified variants in several genes as conferring disease risk for both men and women jointly, while other variants showed evidence for a gender-specific effect. We also identified variants that were associated with quantitative CVD risk factors in both men and women combined, and other variants that showed evidence for gender-genotype interaction. A recent review of gender differences in genetic effects has suggested three criteria for appropriately documented gender differences: (1) The genetic effect is based on the same genetic contrast in both genders; (2) Different genetic subsets in the 2 genders are not compared; and (3) Evidence for a nominally statistically significant gender-gene interaction exists [59]. Our study fulfils all these criteria for the genetic variants showing different effects in men and women. However, studies that replicate these results in larger study samples would be required to confirm or refute the gender-specific associations presented here. With pooling of information across the latest genome-wide association studies [14], [16], [17], [19], [28], there is ample opportunity to test for the presence of gender-genotype interactions behind CHD and ischemic stroke at a genomic level.

In this study, we identified variants in CPB2, F13A1 and LPIN1 as contributing to female-specific risk for CHD and/or CVD, in addition to a variant in USF1 which we have previously reported [31]. Other variants of USF1 have also been reported as showing significant gender-genotype interaction for triglycerides and BMI in familial combined hyperlipidemia families [60]. For ischemic stroke, we identified a variant in F5 as conferring gender-specific risk, in addition to our previously reported association between SEPS1 variants and ischemic stroke in women [30]. Importantly, we identified a larger number of gender-specific effects for women than for men. For men, only one variant in IL6 gene was associated with CVD at p≤0.01 and interaction p<0.01. The asymmetry in positive results is similar to a previous large scale candidate gene study of the metabolic syndrome, in which genetic effects were stronger in women [27]. This is also consistent with the larger heritability estimates for stroke and several CVD-related traits in women [11], [13]. These results suggest that genetic effects on CVD risk may be more readily detectable in women, while for men the genetic effects are more confounded by environmental/lifestyle risk factors.

The most consistent result we identified when analyzing women and men jointly was for a variant in the F12 gene, rs1801020. This promoter variant is located in the untranslated exon 1 of the gene, and the T allele was found to be less common in patients with acute coronary syndrome compared to patients with stable coronary artery disease [61]. In our study sample, in which the A ( = T) allele was associated with risk of CHD and CVD, the study setting was very different, and therefore the results are not readily comparable. Variants in the F12 gene were not present in the Affymetrix 500K and Illumina 300K chips that have been used for the recent genome-wide association studies. The strongest association for quantitative trait variable was between WHR and an intronic variant of the FUT3 gene, rs11673407. The associated variant is not one of the four variants previously associated with Lewis blood phenotype [62] (rs778986 studied here) and which have been reported to be associated with several CVD-related risk factors [63].

Two of the genes we selected to this study, LCT and APOBEC2, have not been previously associated with molecular pathogenesis of cardiovascular disease. We found association between LCT variants and both total and LDL cholesterol. Haplotype analysis implied that the associated variants are in the LCT gene itself, and not necessarily related to the lactase persistence variant upstream of the gene. The C allele of the exonic variant rs2304371, which was associated with highest cholesterol values, is the ancestral allele, present in other mammals and located in a highly conserved region. We also found that a variant in APOBEC2 conferred risk of CHD in men and was associated with higher levels of LDL-C. APOBEC2 belongs to the cytidine deaminase superfamily, and is closely related to APOBEC1 [64]. APOBEC1 mediates the editing of apolipoprotein B mRNA [65]. APOBEC2 is expressed exclusively in heart and skeletal muscle [64], and its function is still largely unknown.

To summarize, we have identified several variants of relevant candidate genes that may confer risk of CHD, ischemic stroke or CVD and/or associate with quantitative CVD-risk factors in a gender-specific manner, and other variants which probably confer risk in both women and men. The identified disease associations and quantitative trait associations had uncorrected p-values≤0.01 for both genders combined and on the basis of the FDR analysis we expect that half of the findings are true positives. For interaction analysis, we may expect that at least third of the reported results are true positives. However, the FDR analysis for the interaction analysis is conservative, because it does not account for the additional criteria we used of association p-value<0.01 for the trait itself in either men or women. Thus, we are convinced that some of the results represent a real effect of variants on disease/trait, but obviously require replication in other studies. In addition, our study had low power to detect genetic effects with HR<1.8 or coefficient<0.3, thus some of the variants we have studied that show no genetic effect might represent false negative results. The possible differences in genetic risk profiles between the genders should be addressed in more detail in genetic studies of CVD, and more focus on female CVD risk is warranted also in genome-wide association studies.

Supporting Information

Table S1.

The genetic variants analyzed in the current study

https://doi.org/10.1371/journal.pone.0003615.s001

(0.09 MB XLS)

Table S2.

Gender-specific results (p≤0.01, uncorrected), in which there was gender-genotype interaction (p≤0.01, uncorrected) in coronary heart disease, ischemic stroke and cardiovascular disease

https://doi.org/10.1371/journal.pone.0003615.s002

(0.06 MB DOC)

Table S3.

Genotype association results for women and men (uncorrected p-values) for coronary heart disease (a), ischemic stroke (b) and cardiovascular disease (c)

https://doi.org/10.1371/journal.pone.0003615.s003

(0.27 MB XLS)

Table S4.

(a) Haplotype analysis for body mass index (BMI) and waist/hip ratio (WHR) for FUT3 variants, sub-cohort men free of CVD at baseline; (b)Haplotype analysis for total cholesterol and LDL cholesterol for LCT variants, sub-cohort men free of CVD at baseline.

https://doi.org/10.1371/journal.pone.0003615.s004

(0.05 MB DOC)

Figure S1.

Power simulations for time-to-event analysis for risk allele frequencies of 0.2 and 0.4, combining both cohorts, using p-value cut-off of 0.01 and assuming for interaction analysis (IA) no effect for men while testing different effect values for women. The lines connect different value points and are not interpolations. HR = hazard ratio.

https://doi.org/10.1371/journal.pone.0003615.s005

(7.13 MB TIF)

Figure S2.

Power simulations for quantitative trait analysis using BMI as an example, testing risk allele frequencies of 0.2 and 0.4, combining both cohorts, using p-value cut-off of 0.01 and assuming for interaction analysis (IA) no effect for men while testing different effect values for women. The lines connect different value points and are not interpolations. Regression coefficients are given in standard deviation scale. BMI = body mass index.

https://doi.org/10.1371/journal.pone.0003615.s006

(7.13 MB TIF)

Acknowledgments

We warmly thank all the Finnish volunteers who participated in the FINRISK-92 and FINRISK-97 studies. KTL Analytical Biochemistry laboratory is acknowledged for the laboratory measurements. Thanks are due to Siv Knaappila, Minna Suvela, Arja Tapio, and Anne Nyberg for their fastidious laboratory work. Minttu Jussila and Elina Mäkinen are acknowledged for DNA aliquotting. Reetta Rautiainen, Tuuli Kauttu, and Aimee Trudeau are acknowledged for partially genotyping the AGTR1, APOA5, and LCT variants, respectively. Zygimantas Cepaitis, Ari Haukijärvi and Dr. Bijoy Joseph from the MORGAM Data Center and Juha Knuuttila and Juri Ahokas from the Genome Informatics Unit at the National Public Health Institute are acknowledged for data handling.

Author Contributions

Conceived and designed the experiments: KS MA KKristiansson KA JK SK KKuulasmaa AE MP VS LP. Performed the experiments: KS MA KKristiansson KA PE JS. Analyzed the data: KS MA KKristiansson OS SR JK SK KKuulasmaa MP LP. Contributed reagents/materials/analysis tools: MN EV PJ VS. Wrote the paper: KS MA KKristiansson OS SR KA JK SK MN PE EV PJ JS KKuulasmaa AE MP VS LP. Responsible for the DNA sample management: MA. Participated in SNP selection: MA KKristiansson KA MP. Participated in genotyping: MA KKristiansson KA PE JS. Participated in data analyses: KS MA KKristiansson SR MP. Responsible for the statistical analysis of the data: OS. Participated in data management, case-cohort selection, and statistical analysis: JK SK. Responsible for phenotype data quality assessment: MN. Participated in recruitment of the FINRISK cohorts: EV PJ VS. Participated in whole genome amplification of DNA samples: JS. Responsible for data management and case-cohort selection: KKuulasmaa. Participated in study design: KS KKuulasmaa AE MP VS. MORGAM Project Coordinator: AE. Participated in the assessment of the clinical diagnoses: VS. Has had an overall responsibility of the study design and data analyses: LP. Responsible for SNP selection: KS. Responsible for DNA sample quality, whole genome amplification of DNA, genotyping, and writing the manuscript: KS. Holds overall responsibility of the data presented: KS.

References

  1. 1. World Health Organization (2006) World health statistics 2006. Geneva: World Health Organization.
  2. 2. Barrett-Connor E (1997) Sex differences in coronary heart disease. Why are women so superior? The 1995 Ancel Keys Lecture. Circulation 95: 252–264.
  3. 3. Jousilahti P, Vartiainen E, Tuomilehto J, Puska P (1999) Sex, age, cardiovascular risk factors, and coronary heart disease: a prospective follow-up study of 14 786 middle-aged men and women in Finland. Circulation 99: 1165–1172.
  4. 4. Pilote L, Dasgupta K, Guru V, Humphries KH, McGrath J, et al. (2007) A comprehensive view of sex-specific issues related to cardiovascular disease. Cmaj 176: S1–44.
  5. 5. Tunstall-Pedoe H (1998) Myth and paradox of coronary risk and the menopause. Lancet 351: 1425–1427.
  6. 6. Tunstall-Pedoe H, Kuulasmaa K, Amouyel P, Arveiler D, Rajakangas AM, et al. (1994) Myocardial infarction and coronary deaths in the World Health Organization MONICA Project. Registration procedures, event rates, and case-fatality rates in 38 populations from 21 countries in four continents. Circulation 90: 583–612.
  7. 7. Yusuf S, Hawken S, Ounpuu S, Dans T, Avezum A, et al. (2004) Effect of potentially modifiable risk factors associated with myocardial infarction in 52 countries (the INTERHEART study): case-control study. Lancet 364: 937–952.
  8. 8. Bak S, Gaist D, Sindrup SH, Skytthe A, Christensen K (2002) Genetic liability in stroke: a long-term follow-up study of Danish twins. Stroke 33: 769–774.
  9. 9. Evans A, Van Baal GC, McCarron P, DeLange M, Soerensen TI, et al. (2003) The genetics of coronary heart disease: the contribution of twin studies. Twin Res 6: 432–441.
  10. 10. Lusis AJ, Mar R, Pajukanta P (2004) Genetics of atherosclerosis. Annu Rev Genomics Hum Genet 5: 189–218.
  11. 11. Touze E, Rothwell PM (2007) Heritability of ischaemic stroke in women compared with men: a genetic epidemiological study. Lancet Neurol 6: 125–133.
  12. 12. McCarthy JJ (2007) Gene by sex interaction in the etiology of coronary heart disease and the preceding metabolic syndrome. Nutr Metab Cardiovasc Dis 17: 153–161.
  13. 13. Pilia G, Chen WM, Scuteri A, Orru M, Albai G, et al. (2006) Heritability of cardiovascular and personality traits in 6,148 Sardinians. PLoS Genet 2: e132.
  14. 14. Helgadottir A, Thorleifsson G, Manolescu A, Gretarsdottir S, Blondal T, et al. (2007) A common variant on chromosome 9p21 affects the risk of myocardial infarction. Science 316: 1491–1493.
  15. 15. Larson MG, Atwood LD, Benjamin EJ, Cupples LA, D'Agostino RB Sr, et al. (2007) Framingham Heart Study 100K project: genome-wide associations for cardiovascular disease outcomes. BMC Med Genet 8: Suppl 1S5.
  16. 16. McPherson R, Pertsemlidis A, Kavaslar N, Stewart A, Roberts R, et al. (2007) A common allele on chromosome 9 associated with coronary heart disease. Science 316: 1488–1491.
  17. 17. Samani NJ, Erdmann J, Hall AS, Hengstenberg C, Mangino M, et al. (2007) Genomewide association analysis of coronary artery disease. N Engl J Med 357: 443–453.
  18. 18. Wang MM (2006) Genetics of ischemic stroke: future clinical applications. Semin Neurol 26: 523–530.
  19. 19. Wellcome Trust Case Control Consortium (2007) Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 447: 661–678.
  20. 20. McCarthy JJ, Parker A, Salem R, Moliterno DJ, Wang Q, et al. (2004) Large scale association analysis for identification of genes underlying premature coronary heart disease: cumulative perspective from analysis of 111 candidate genes. J Med Genet 41: 334–341.
  21. 21. Ozaki K, Ohnishi Y, Iida A, Sekine A, Yamada R, et al. (2002) Functional SNPs in the lymphotoxin-alpha gene that are associated with susceptibility to myocardial infarction. Nat Genet 32: 650–654.
  22. 22. Pare G, Serre D, Brisson D, Anand SS, Montpetit A, et al. (2007) Genetic analysis of 103 candidate genes for coronary artery disease and associated phenotypes in a founder population reveals a new association between endothelin-1 and high-density lipoprotein cholesterol. Am J Hum Genet 80: 673–682.
  23. 23. Topol EJ, McCarthy J, Gabriel S, Moliterno DJ, Rogers WJ, et al. (2001) Single nucleotide polymorphisms in multiple novel thrombospondin genes may be associated with familial premature myocardial infarction. Circulation 104: 2641–2644.
  24. 24. Yamada Y, Matsuo H, Segawa T, Watanabe S, Kato K, et al. (2006) Assessment of genetic risk for myocardial infarction. Thromb Haemost 96: 220–227.
  25. 25. Yamada Y, Metoki N, Yoshida H, Satoh K, Ichihara S, et al. (2006) Genetic risk for ischemic and hemorrhagic stroke. Arterioscler Thromb Vasc Biol 26: 1920–1925.
  26. 26. Yamada Y, Izawa H, Ichihara S, Takatsu F, Ishihara H, et al. (2002) Prediction of the risk of myocardial infarction from polymorphisms in candidate genes. N Engl J Med 347: 1916–1923.
  27. 27. McCarthy JJ, Meyer J, Moliterno DJ, Newby LK, Rogers WJ, et al. (2003) Evidence for substantial effect modification by gender in a large-scale genetic association study of the metabolic syndrome among coronary heart disease patients. Hum Genet 114: 87–98.
  28. 28. Matarin M, Brown WM, Scholz S, Simon-Sanchez J, Fung HC, et al. (2007) A genome-wide genotyping study in patients with ischaemic stroke: initial analysis and data release. Lancet Neurol 6: 414–420.
  29. 29. Evans A, Salomaa V, Kulathinal S, Asplund K, Cambien F, et al. (2005) MORGAM (an international pooling of cardiovascular cohorts). Int J Epidemiol 34: 21–27.
  30. 30. Alanne M, Kristiansson K, Auro K, Silander K, Kuulasmaa K, et al. (2007) Variation in the selenoprotein S gene locus is associated with coronary heart disease and ischemic stroke in two independent Finnish cohorts. Hum Genet 122: 355–365.
  31. 31. Komulainen K, Alanne M, Auro K, Kilpikari R, Pajukanta P, et al. (2006) Risk alleles of USF1 gene predict cardiovascular disease of women in two prospective studies. PLoS Genet 2: e69.
  32. 32. Vartiainen E, Jousilahti P, Alfthan G, Sundvall J, Pietinen P, et al. (2000) Cardiovascular risk factor changes in Finland, 1972–1997. Int J Epidemiol 29: 49–56.
  33. 33. Auro K, Komulainen K, Alanne M, Silander K, Peltonen L, et al. (2006) Thrombomodulin gene polymorphisms and haplotypes and the risk of cardiovascular events: a prospective follow-up study. Arterioscler Thromb Vasc Biol 26: 942–947.
  34. 34. Kulathinal S, Karvanen J, Saarela O, Kuulasmaa K, for the MORGAM Project (2007) Case-cohort design in practice–experiences from the MORGAM Project. Epidemiol Perspect Innov 4: 15.
  35. 35. Kulathinal S, NM , Kuulasmaa K, contributors from Participating Centresfor the MORGAM Project (2005) Description of MORGAM Cohorts. Available from: http://www.ktl.fi/publications/morgam/cohorts/index.html (access date 18.12.2007).
  36. 36. Silander K, Komulainen K, Ellonen P, Jussila M, Alanne M, et al. (2005) Evaluating whole genome amplification via multiply-primed rolling circle amplification for SNP genotyping of samples with low DNA yield. Twin Res Hum Genet 8: 368–375.
  37. 37. Carlson CS, Eberle MA, Rieder MJ, Yi Q, Kruglyak L, et al. (2004) Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am J Hum Genet 74: 106–120.
  38. 38. Enattah NS, Sahi T, Savilahti E, Terwilliger JD, Peltonen L, et al. (2002) Identification of a variant associated with adult-type hypolactasia. Nat Genet 30: 233–237.
  39. 39. Fullerton SM, Clark AG, Weiss KM, Nickerson DA, Taylor SL, et al. (2000) Apolipoprotein E variation at the sequence haplotype level: implications for the origin and maintenance of a major human polymorphism. Am J Hum Genet 67: 881–900.
  40. 40. Nickerson DA, Taylor SL, Fullerton SM, Weiss KM, Clark AG, et al. (2000) Sequence diversity and large-scale typing of SNPs in the human apolipoprotein E gene. Genome Res 10: 1532–1545.
  41. 41. Suviolahti E, Reue K, Cantor RM, Phan J, Gentile M, et al. (2006) Cross-species analyses implicate Lipin 1 involvement in human glucose metabolism. Hum Mol Genet 15: 377–386.
  42. 42. Tikka-Kleemola P, Hamalainen E, Tuomainen K, Suvela M, Artma A, et al. (2007) The enhancement of homogenous mass extension reaction: comparison of two enzymes. Mol Cell Probes 21: 216–221.
  43. 43. Janis MT, Siggins S, Tahvanainen E, Vikstedt R, Silander K, et al. (2004) Active and low-active forms of serum phospholipid transfer protein in a normal Finnish population sample. J Lipid Res 45: 2303–2309.
  44. 44. Auro K, Alanne M, Kristiansson K, Silander K, Kuulasmaa K, et al. (2007) Combined effects of thrombosis pathway gene variants predict cardiovascular events. PLoS Genet 3: e120.
  45. 45. Silander K, Axelsson T, Widen E, Dahlgren A, Palotie A, et al. (2003) Analysis of genetic variation in the GenomEUtwin project. Twin Res 6: 391–398.
  46. 46. O'Connell JR, Weeks DE (1998) PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet 63: 259–266.
  47. 47. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–265.
  48. 48. Prentice RL (1986) A case-cohort design for epidemiologic cohort studies and disease prevention trials. Biometrika 73: 1–11.
  49. 49. Sinnwell JP, SD , Yu Z (2007) The haplo.stats Package, version 1.3.1.
  50. 50. Efron B, Tibshirani R, Storey JD, Tusher V (2001) Empirical Bayes analysis of a microarray experiment. Journal of the American Statistical Association 96: 1151–1160.
  51. 51. Efron B (2004) Large-scale simultaneous hypothesis testing: The choice of a null hypothesis. Journal of the American Statistical Association 99: 96–104.
  52. 52. Strimmer K (2007) The fdrtool Package, version 1.1.4.
  53. 53. Russo F, De Carne M, Buonsante A, Cisternino AM, Centonze S, et al. (1997) Hypolactasia and metabolic changes in post-menopausal women. Maturitas 26: 193–202.
  54. 54. Sahi T, Jussila J, Penttila IM, Sarna S, Isokoski M (1977) Serum lipids and proteins in lactose malabsorption. Am J Clin Nutr 30: 476–481.
  55. 55. Enattah NS, Trudeau A, Pimenoff V, Maiuri L, Auricchio S, et al. (2007) Evidence of still-ongoing convergence evolution of the lactase persistence T-13910 alleles in humans. Am J Hum Genet 81: 615–625.
  56. 56. Rinn JL, Snyder M (2005) Sexual dimorphism in mammalian gene expression. Trends Genet 21: 298–305.
  57. 57. Mendelsohn ME, Karas RH (2005) Molecular and cellular basis of cardiovascular gender differences. Science 308: 1583–1587.
  58. 58. Schunkert H, Gotz A, Braund P, McGinnis R, Tregouet DA, et al. (2008) Repeated Replication and a Prospective Meta-Analysis of the Association Between Chromosome 9p21.3 and Coronary Artery Disease. Circulation.
  59. 59. Patsopoulos NA, Tatsioni A, Ioannidis JP (2007) Claims of sex differences: an empirical assessment in genetic associations. Jama 298: 880–893.
  60. 60. Lee JC, Weissglas-Volkov D, Kyttälä M, Sinsheimer JS, Jokiaho A, et al. (2007) USF1 Contributes to High Serum Lipid Levels in Dutch FCHL Families and U.S. Whites With Coronary Artery Disease. Arterioscler Thromb Vasc Biol. In press.
  61. 61. Endler G, Mannhalter C, Sunder-Plassmann H, Lalouschek W, Kapiotis S, et al. (2001) Homozygosity for the C→T polymorphism at nucleotide 46 in the 5' untranslated region of the factor XII gene protects from development of acute coronary syndrome. Br J Haematol 115: 1007–1009.
  62. 62. Salomaa V, Pankow J, Heiss G, Cakir B, Eckfeldt JH, et al. (2000) Genetic background of Lewis negative blood group phenotype and its association with atherosclerotic disease in the NHLBI family heart study. J Intern Med 247: 689–698.
  63. 63. Cakir B, Heiss G, Pankow JS, Salomaa V, Sharrett AR, et al. (2004) Association of the Lewis genotype with cardiovascular risk factors and subclinical carotid atherosclerosis: the Atherosclerosis Risk in Communities (ARIC) study. J Intern Med 255: 40–51.
  64. 64. Liao W, Hong SH, Chan BH, Rudolph FB, Clark SC, et al. (1999) APOBEC-2, a cardiac- and skeletal muscle-specific member of the cytidine deaminase supergene family. Biochem Biophys Res Commun 260: 398–404.
  65. 65. Lau PP, Zhu HJ, Baldini A, Charnsangavej C, Chan L (1994) Dimeric structure of a human apolipoprotein B mRNA editing protein and cloning and chromosomal localization of its gene. Proc Natl Acad Sci U S A 91: 8522–8526.