- Research article
- Open Access
Coronary risk in relation to genetic variation in MEOX2 and TCF15 in a Flemish population
BMC Geneticsvolume 16, Article number: 116 (2015)
In mice MEOX2/TCF15 heterodimers are highly expressed in heart endothelial cells and are involved in the transcriptional regulation of lipid transport. In a general population, we investigated whether genetic variation in these genes predicted coronary heart disease (CHD).
In 2027 participants randomly recruited from a Flemish population (51.0 % women; mean age 43.6 years), we genotyped six SNPs in MEOX2 and four in TCF15. Over 15.2 years (median), CHD, myocardial infarction, coronary revascularisation and ischaemic cardiomyopathy occurred in 106, 53, 78 and 22 participants. For SNPs, we contrasted CHD risk in minor-allele heterozygotes and homozygotes (variant) vs. major-allele homozygotes (reference) and for haplotypes carriers (variant) vs. non-carriers. In multivariable-adjusted analyses with correction for multiple testing, CHD risk was associated with MEOX2 SNPs (P ≤ 0.049), but not with TCF15 SNPs (P ≥ 0.29). The MEOX2 GTCCGC haplotype (frequency 16.5 %) was associated with the sex- and age-standardised CHD incidence (5.26 vs. 3.03 events per 1000 person-years; P = 0.036); the multivariable-adjusted hazard ratio [HR] of CHD was 1.78 (95 % confidence interval, 1.25–2.56; P = 0.0054). For myocardial infarction, coronary revascularisation, and ischaemic cardiomyopathy, the corresponding HRs were 1.96 (1.16–3.31), 1.87 (1.20–2.91) and 3.16 (1.41–7.09), respectively. The MEOX2 GTCCGC haplotype significantly improved the prediction of CHD over and beyond traditional risk factors and was associated with similar population-attributable risk as smoking (18.7 % vs. 16.2 %).
Genetic variation in MEOX2, but not TCF15, is a strong predictor of CHD. Further experimental studies should elucidate the underlying molecular mechanisms.
Endothelial cells lining the microvasculature constitute the interface between the circulating blood and tissues . They differentiate to acquire the molecular, morphological and functional characteristics required for proper organ function . In the heart, endothelial cells play an active role in the transport of fatty acids, the principal energy source for the continuously beating muscle [1, 2]. Using microarray profiling on endothelial cells isolated from the heart, brain, and liver of mice, we recently identified a specific genetic signature for heart endothelial cells, including MEOX2/TCF15 heterodimers as novel transcriptional determinants . This signature was largely shared with skeletal muscle and adipose tissue endothelium and was enriched in genes encoding fatty acid transport-related proteins . Using gain- and loss-of-function approaches, we showed that MEOX2/TCF15 mediates fatty acid uptake in heart endothelial cells, in part, by driving endothelial CD36 and lipoprotein lipase (LPL) expression and thereby facilitating fatty acid transport across cardiac endothelial cells .
LPL is expressed at the luminal endothelial surface of arteries and capillaries and hydrolyses circulating lipoprotein triglycerides into free fatty acids and glycerol . Local  and systemic  dysregulation of lipid metabolism and endothelial dysfunction [6, 7] are hallmarks of coronary atherosclerosis and long precede clinically overt disease. These observations suggest that genetic predisposition plays an important role in the pathogenesis of coronary heart disease (CHD) [6, 7]. In view of our recent observations of enriched expression of MEOX2 and TCF15 in heart endothelial cells , we hypothesised that genetic variation in the genes encoding these transcription factors might be associated with coronary risk. To test this hypothesis, we analysed data accumulated since 1985 in a Flemish population study [8, 9].
The Flemish Study on Environment, Genes and Health Outcomes (FLEMENGHO) complies with the Helsinki declaration for research in human subjects and the Belgian legislation for the protection of privacy (http://www.privacycommission.be). The Ethics Committee of the University of Leuven approved the study. Recruitment for the FLEMENGHO study started in 1985 [8, 9]. From August 1985 to November 1990, a random sample of the households living in a geographically defined area of Northern Belgium was investigated with the goal to recruit an equal number of participants in each of six strata by sex and age (20–39, 40–59, and ≥60 years). All household members aged 20 years or older were invited, if the quota of their sex-age group had not yet been met. From June 1996 until January 2004 recruitment of families continued using the former participants (1985–1990) as index persons and including teenagers. The participants were repeatedly followed up. In all study phases, we used the same standardised methods to measure blood pressure and to administer questionnaires. The participation rate at enrolment was 78.0 %. At each contact, participants gave or renewed informed written consent.
Of 3343 enrolled participants, we excluded 1316 from analysis, because blood stored in the biobank was exhausted with no material left for genotyping (n = 521), because of DNA degradation (n = 314), because at enrolment they were less than 20 years old (n = 372), or because one or more of the six MEOX2 or four TCF15 SNPs were unavailable (n = 109). Thus, the number of participants statistically analysed totalled 2027.
Measurements at baseline
Trained nurses measured the participants’ anthropometric characteristics and blood pressure. Body mass index was weight in kilograms divided by the square of height in meters. Blood pressure was the average of five consecutive auscultatory readings obtained with a standard mercury sphygmomanometer after participants had rested in the sitting position for at least 5 min. Hypertension was a blood pressure of at least 140 mm Hg systolic or 90 mm Hg diastolic, or use of antihypertensive drugs. The nurses also administered a standardised questionnaire inquiring about each participant’s medical history, smoking and drinking habits, and intake of medications. Plasma glucose and serum total and high-density lipoprotein (HDL) cholesterol and serum creatinine were measured by automated methods in certified laboratories. Diabetes mellitus was a fasting or random plasma glucose level exceeding 7.0 or 11.1 mmol/L, or use of antidiabetic agents .
Ascertainment of coronary events
FLEMENGHO received ethical approval. The database was registered with the Privacy Commission. These legal requirements being fulfilled, we could ascertain the vital status of participants at annual intervals until 06 December 2012 via the Belgian Population Registry. In addition, we could obtain the International Classification of Disease codes for the immediate and underlying causes of death from the Flemish Registry of Death Certificates. For 1853 participants, we collected information on the incidence of non-fatal endpoints either via face-to-face follow-up visits with repeated administration of the same standardised questionnaire as used at baseline (n = 1521) or via a structured telephone interview (n = 332). Follow-up data were available from one visit in 360 participants, from two in 304, from three in 436, and from four or more in 421 participants.
Trained nurses used the International Classification of Diseases to code incident cases of CHD. Two investigators blinded with regard to the genotypic results adjudicated all coronary events against the medical records of general practitioners or hospitals. Coronary events included sudden death, fatal and non-fatal myocardial infarction, acute coronary syndrome requiring hospitalisation, ischaemic cardiomyopathy, and surgical or percutaneous coronary revascularisation. In the outcome analyses, we only considered the first event within each category.
Ethics approval and informed consent covered genotyping. After DNA extraction from peripheral blood , SNPs were genotyped using the TaqMan® OpenArray™ Genotyping System (Life Technologies, Foster City, CA). All DNA samples were loaded at 50 ng/mL and amplified according to the manufacturer’s instructions. For analysis of the genotypes, we used autocalling methods, as implemented in the TaqMan Genotyper software version 1.3 (Life Technologies). Next, genotype clusters were evaluated manually with the sample call rate set above 0.90. Sixteen duplicate samples gave 100 % reproducibility for all 64 SNPs on the custom made array, including the genes of interest in the current article .
MEOX2 (75601 base pairs) maps to chromosome 7 (p22.1–p21.3). To select the MEOX2 SNPs to genotype, we first reviewed all SNPs in this gene, including the flanking regions, as available in the Illumina 1 M Duo and OmniExpress arrays (San Diego, CA). We excluded SNPs with a minor allele frequency of less than 1 % and those that were in high linkage disequilibrium (r 2 ≥ 0.80). Next, based on the availability of SNPs on the TaqMan OpenArray Genotyping System, we selected 12 tagging SNPs (rs6946099, rs10777, rs7800473, rs13438001, rs12056299, rs7787043, rs758297, rs4532497, rs10263561, rs6959056, rs740566, rs1050290) that are in high linkage disequilibrium (r 2 ≥ 0.80) with 92 neighbouring SNPs (Additional file 1: Figure S1 and Table S1), but were not in high linkage disequilibrium (r 2 < 0.80) with one another. The 12 selected SNPs covered the entire gene with extension into the 3’ and 5’ flanking regions. We excluded six SNPs with a successful genotyping call rate of less than 0.98. Finally, we retained six MEOX2 SNPs (rs10777, rs12056299, rs7787043, rs4532497, rs6959056, and rs1050290) in the analysis (Additional file 1: Table S2) that are in linkage disequilibrium (r 2 > 0.80) with 23 other SNPs (Additional file 1: Table S1). TCF15 (6602 base pairs) maps to chromosome 20p13. We genotyped five SNPs covering the whole gene (rs282152, rs6116745, rs282162, rs3761308 and rs12624577), but excluded rs282152, because the SNP call rate was less than 0.98 (Additional file 1: Figure S2).
For database management and statistical analysis, we used SAS software, version 9.3 (SAS Institute, Cary, NC). For comparison of means and proportions, we applied the large sample z-test or ANOVA and Fisher’s exact, respectively. We tested Hardy-Weinberg equilibrium in unrelated founders, using the exact statistics available in the PROC ALLELE procedure of the SAS package. For analysis of single SNPs, we combined the least frequent homozygous group with heterozygous subjects. We tested linkage disequilibrium and reconstructed haplotypes using the SAS procedures PROC ALLELE and PROC HAPLOTYPE. To check for consistency, we repeated haplotype construction accounting for pedigree information using SHAPEIT version 2 (http://mathgen.stats.ox.ac.uk/genetics_software/shapeit/shapeit.html ).
We compared the incidence of coronary endpoints in relation to genetic variants, using (i) rates standardised by the direct method for sex and age (<40, 40–59, ≥60 years) and (ii) the cumulative incidence derived from Cox models adjusted for sex and age. Next, we assessed the prognostic value of the genetic variants in multivariable-adjusted Cox regression. We checked the proportional hazard assumption by applying a Kolmogorov-type supremum test as implemented in the ASSESS statement of the PROC PHREG procedure. To account for family clusters, we used the PROC SURVIVAL procedure of the SUDAAN 11.0.1 software (Research Triangle Institute, NC). In this procedure, non-independence among family members was taken into account by including family as a random effect along with other covariables as fixed effects. We analysed genotypes and haplotypes using major allele homozygotes and non-carriers as the reference group, respectively. We adjusted P values for the associations between outcomes and genetic variants, using the Benjamini and Hochberg false discovery rate  according to the number of SNPs retained in the analysis.
We computed the positive predictive value of the risk carrying MEOX2 haplotype GTCCGC as (R × D)/([G/100] × [R −1] + 1), where R is the multivariable-adjusted hazard ratio, D is the incidence of CHD in the whole population, and G is the prevalence of the GTCCGC haplotype . The attributable risk is given by ([R −1] × 100)/R and the population-attributable risk by ([G/100] × [R −1] × 100)/([G/100] × [R −1] + 1) . Finally, we assessed the power of the MEOX2 GTCCGC haplotype to predict CHD over and beyond classical risk factors, using the integrated discrimination improvement (IDI) and the net reclassification improvement (NRI), as described by Pencina and colleagues for survival data .
All 2027 participants were White Europeans, of whom 1034 (51.0 %) were women. The study population consisted of 332 singletons and 1695 related subjects, belonging to 49 single-generation families and 191 multi-generation pedigrees. Age averaged (±SD) 43.6 ± 14.3 years, blood pressure 125.0 ± 15.4 mm Hg systolic and 76.2 ± 9.5 mm Hg diastolic, body mass index 25.7 ± 4.3 kg/m2, and total cholesterol 5.49 ± 1.15 mmol/L. Among all participants, 486 (24.0 %) had hypertension, of whom 214 (44.0 %) were on antihypertensive drug treatment, 33 (1.6 %) had diabetes mellitus, and 41 (2.0 %) reported a history of CHD. Previous coronary complications included angiographically proven coronary stenosis, myocardial infarction, and coronary revascularisation in 8 (0.4 %), 11 (0.5 %) and 22 (1.1 %) patients, respectively. Of 1034 women and 993 men, 277 (26.8 %) women and 328 (33.0 %) men were smokers, and 168 (16.3 %) women and 418 (42.1 %) men reported intake of alcohol. In smokers, median tobacco use was 15 cigarettes per day (interquartile range, 10 to 20 cigarettes per day). In drinkers, the median alcohol consumption was 14 g per day (8 to 26 g per day).
Table 1 lists the baseline characteristics of participants according to CHD incidence. Most risk factors differed in the expected direction between cases and non-cases. However, compared with non-cases, the prevalence of smoking was not different in patients with incident CHD (29.6 % vs. 36.8 %; P = 0.13), while the prevalence of drinking was lower among cases (29.4 % vs. 19.8 %; P = 0.036). Heart rate at baseline was similar in participants without and with incident CHD (69.2 vs. 69.9 beats per minute; P = 0.43).
Incidence of events
Over a median follow-up of 15.2 years (5th to 95th percentile interval, 5.7 to 27.1 years), 106 new coronary events occurred, 24 fatal and 82 non-fatal. Coronary events comprised 12 fatal and 34 non-fatal myocardial infarcts and 7 sudden deaths. There were 78 patients who underwent surgical (n = 29) or percutaneous (n = 56) coronary revascularisation. Coronary events also included 5 fatal and 17 non-fatal cases of ischaemic cardiomyopathy.
Analyses of single SNPs in MEOX2 and TCF15
Additional file 1: Table S2 describes the position and the SNPs retained in the analysis and the allele and genotype frequencies in 825 unrelated founders. The six SNPs in MEOX2 and the four SNPs in TCF15 complied with Hardy-Weinberg equilibrium (0.30 ≤ P ≤ 0.80). In the whole study population (Additional file 1: Table S3), the frequencies of the minor alleles ranged from 21.1 to 43.0 % for MEOX2 and from 9.6 to 38.5 % for TCF15. The prevalence of minor allele homozygotes ranged from 4.4 to 17.8 % for MEOX2, and from 0.7 to 15.9 % for TCF15.
The sex- and age-standardised incidence rates of coronary events associated with the MEOX2 SNPs appear in Additional file 1: Table S4. Compared with major allele homozygotes, minor allele carriers experienced a higher CHD incidence except for rs6959056. The sex- and age-adjusted cumulative incidence of coronary events (Fig. 1) showed significant association (P ≤ 0.012) with the MEOX2 SNPs except for rs1050290 (P = 0.058). There were no differences in these estimates between homozygous and heterozygous minor allele carriers (0.23 ≤ P ≤ 0.98) except for rs12056299 (P = 0.014). For all coronary events combined, the sex- and age-standardised incidence rates (0.11 ≤ P ≤ 0.39) and the sex- and age-adjusted cumulative incidence (0.11 ≤ P ≤ 0.71) did not differ among minor allele carriers and major allele homozygotes of the four TCF15 SNPs.
Next, we accounted for family clusters and adjusted the hazard ratios for baseline characteristics, including sex, age, body mass index, systolic pressure, the total-to-HDL cholesterol ratio, smoking and drinking, and antihypertensive drug treatment. Compared with homozygotes of the major allele, for rs10777, rs12056299, rs7787043, rs4532497, and rs1050290, CHD risk was higher in minor allele carriers, whereas the opposite was the case for rs6959056 (Table 2). These findings remained consistent with correction for multiple testing (Table 2) and after excluding patients who had a history of CHD at baseline (Additional file 1: Table S5).
For the four TCF15 SNPs, the multivariable-adjusted hazard ratios modelling the CHD risk of minor allele carriers vs. major allele homozygotes did not reach significance (0.75 ≤ hazard ratio ≤ 1.45; 0.072 ≤ P ≤ 0.52). However, Additional file 1: Figure S3 shows interaction (P = 0.011) between rs12624577 in TCF15 and rs4532497 in MEOX2. Among C allele carriers of MEOX2 rs4532497, the hazard ratio expressing the risk of the C allele relative to the TT genotype in TCF15 was 2.44 (95 % confidence interval, 1.38–4.29; P = 0.0021), whereas among MEOX2 TT homozygotes the corresponding hazard ratio was 0.75 (0.37–1.51; P = 0.42).
Analysis of MEOX2 haplotypes
Using the expectation-maximisation algorithm as implemented in the PROC HAPLOTYPE procedure of the SAS software version 9.3, three haplotypes of MEOX2 SNPs had a frequency of over 10 % and were carried through in the analysis. With letters referring to the alleles in rs10777, rs12056299, rs7787043, rs4532497, rs6959056, and rs1050290 in MEOX2, these haplotypes were TCTTAT (27.5 %), TCTTGT (26.4 %), and GTCCGC (16.5 %). For all coronary events, rates standardised for sex and age (5.26 vs. 3.03 events per 1000 person-years; Additional file 1: Table S6) and cumulative incidence estimates adjusted for sex and age (Additional file 1: Figure S4), both before (P ≤ 0.012) and after adjustment for multiple testing (P ≤ 0.036), showed significant association with GTCCGC. In multivariable-adjusted analyses (Table 3), GTCCGC carriers, compared to non-carriers, had a 78 % higher CHD risk. The p-values before and after correction for multiple testing were 0.0018 and 0.0054, respectively. Myocardial infarction, coronary revascularisation and ischaemic cardiomyopathy all showed significant association with the GTCCGC haplotype. These findings were not materially different when we reconstructed haplotypes accounting for pedigree information (Additional file 1: Table S7) or after excluding patients who had a history of CHD at baseline (Additional file 1: Table S8).
With adjustments applied as before, the positive predictive value and the attributable and population-attributable CHD risks associated with GTCCGC were 7.5, 43.1 and 18.7 %, respectively. For smoking, analysed as reference, the corresponding estimates (unadjusted for genetic risk) were 7.2, 39.3 and 16.2 %, respectively. Table 4 shows that in all participants and in those without CHD at entry, adding the GTCCGC haplotype to the basic model including traditional risk factors improved (0.016 ≤ P ≤ 0.056) IDI and NRI.
To our knowledge, our study is the first to relate in a general population CHD incidence to genetic variation in MEOX2 and TCF15, two transcription factors that are highly expressed by cardiac endothelium and that in a heterodimeric fashion interfere with cardiac energy metabolism by driving endothelial CD36 and LPL expression, thereby facilitating fatty acid transport across the cardiac endothelium . The key finding of our current study was that the risk of advanced CHD was associated with genetic variation in MEOX2, as captured by six tagging SNPs. On the other hand, genetic variation in TCF15, coding for the heterodimeric partner of MEOX2, was not associated with the incidence of coronary events. Nonetheless, the CHD risk associated with MEOX2 rs4532497 was confined to TCF15 rs12624577 variant allele carriers, which might reflect the known heterodimeric action picked up in our experimental studies . Although our current study firmly established an association between CHD risk and genetic variation in MEOX2, the molecular mechanisms underlying this relation need further clarification in experimental studies back translating our epidemiological findings. For now, working hypotheses might be developed along two lines respectively involving disturbed lipid handling [5, 17–19], a key mechanism in atherosclerosis, or the involvement of MEOX2 in the angiogenic responses to stressors [20–24] or in the migration or proliferation of endothelial and vascular smooth muscle cells [25, 26].
A large-scale genome-wide association study identified 46 significant lead SNPs associated with CHD. Twelve showed a significant association with a lipid trait. Variation in the MEOX2 gene was not among these SNPs, but LPL (rs264) was . LPL catalyses the hydrolysis of triglycerides in plasma triglyceride-rich lipoproteins, chylomicrons and very low density lipoproteins at the capillary endothelial cell surface, providing free fatty acids and glycerol as energy source for tissues . Genetic variation in LPL is associated with the levels of circulating LPL activity , the plasma concentration of triglycerides [5, 18] and HDL cholesterol [5, 18] and in some , albeit not all , studies with the risk of CHD. Parenchymal cells in adipose, skeletal and cardiac muscle widely express LPL throughout the body . Cardiac endothelial cells have a particular expression profile including MEOX2, a gene that facilitates fatty acid transport into and through heart endothelial cells. As for genetic variation in LPL [5, 18], we hypothesised that dysregulation of lipid transport in cardiac endothelial cells might increase CHD risk. Important in this regard is that in our current study, in contrast to the studies on genetic variants of LPL [5, 18], we did not find any association of the circulating lipid levels with MEOX2 variants (data not shown), pointing to a local coronary rather than a systemic underlying mechanism.
Moving to the second hypothetical pathophysiological pathway, MEOX2 is also known as growth arrest-specific homeobox (GAX) , During embryonic development, the three muscle lineages express GAX . In adult life, vascular smooth muscle cells also express GAX . Mitogenic stimuli, such as platelet-derived growth factor and angiotensin II or injury of the endothelium , inhibit GAX expression, whereas growth arrest signals, such as serum deprivation of cultured cells, enhance its expression and negatively regulate the cell cycle . Observations in transfected cells also point to MEOX2 as a potentially important regulatory gene inhibiting not only the angiogenic response of endothelial cells to pro-angiogenic factors, but also their response to chronic inflammatory stimulation that normally activates NF-κB . Inflammatory pathways identified in a network analysis of 233 candidate genes play key roles in development of coronary atherosclerosis . These observations [22–24, 27] may offer an alternative explanation why in our current study coronary risk was associated with genetic variation in MEOX2. Chen and Gorski did an in silico search for micro-RNA binding sites in the GAX 5’UTR and identified consensus sites for multiple candidate micro-RNAs, of which only miR-130a was expressed in proliferating endothelial cells . miR-130a was largely responsible for the down-regulation of GAX expression in response to mitogens and pro-angiogenic factors and antagonised the antiangiogenic activity of GAX .
To our knowledge, previously published GWAS results did not demonstrate association between coronary heart disease and MEOX2. However, these GWAS studies relied on comparing cases and controls drawn from heterogeneous sources [27–29] or on a retrospective cross-sectional analysis of patients referred for coronary angiography . GWAS case–control studies offer the opportunity for searching for association between CHD and densely distributed SNPs across the whole genome in large numbers of patients and controls. Such studies require significance levels of 10-6 to 10-8. In contrast, our study was prospective and population-based and tested a prior hypothesis involving only six SNPs in MEOX2 and four in TCF15. We did therefore not rely on such extreme P-values, but applied the Benjamin-Hochberg approach for multiple testing. Admittedly, our sample size was smaller than in the GWAS studies. This is particularly relevant for the interaction between MEOX2 rs4532497 and TCF15 rs12624577 (Additional file 1: Figure S3), a finding, which although in line with our experimental findings  can only be considered as hypothesis generating. Future population-based research projects might address this issue.
We performed the annotation of the genomic context surrounding the SNPs retained for MEOX2, based on the data of the ENCODE project (http://www.genome.gov/encode). As shown in Additional file 1: Table S1, rs10777 and rs1050290 map into the 3’UTR and 5’UTR regions, respectively, whereas rs12056299, rs7787043, rs4532497 and rs6959056 are in the first intron of MEOX2. Moreover, according to Ensemble 75 annotation (http://www.ensembl.org/Homo_sapiens/Info/Annotation) and GENCODE 22 (http://www.gencodegenes.org/releases/22.html), rs4532497 and rs6959056 also map into the ENSG00000237070 (AC005550.3) antisense non-protein coding gene. In particular, rs4532497 maps into intron 1, whereas rs6959056 maps into exon 4 of ENSG00000237070. rs6959056 and rs1050290 fall into a promoter regulatory region (ENST00000622287) both in human umbilical vein endothelial cells and in human dermal fibroblasts, where MEOX2 is expressed. Further functional studies are needed to investigate the possible modulation of MEOX2 expression. Moreover, rs6959056 in exon 4 of the non-coding transcript ENST00000451240 could affect gene function, although no information on the ENSG00000237070 gene is currently available.
TCF15, also known as Paraxis, is a member of the Twist subfamily of basic helix-loop-helix transcription factors that regulate specification of mesodermal derivatives during vertebrate embryogenesis . TCF15 primes pluripotent cells for differentiation . During dermomyotome formation in Xenopus laevis, TCF15 directly activates the expression of MEOX2 . Our experimental studies demonstrated that the MEOX2/TCF15 heterodimer facilitates the transport of fatty acids across cardiac endothelial cells and that in mice haplodeficiency in these genes results in impaired contractility of cardiomyocytes and heart failure . In our population study, we therefore also searched for association between the incidence of well-documented heart failure and variation in the MEOX2 and TCF15 genes. Several reasons may explain why such associations were not detected. Indeed, heart failure is a heterogeneous disease caused by a multitude of instigators, including ischaemic or valvular heart disease, comorbidities, or risk factors such as hypertension. Moreover, the diagnosis of heart failure depends on the clinical interpretation of a combination of signs and symptoms, that are difficult to recognise [33, 34].
The present study must be interpreted within the context of some potential limitations. First, even though our analysis was hypothesis-driven based on published evidence from experimental studies , we adjusted significance levels for multiple testing according to the number of SNPs tested, using the Benjamini and Hochberg false discovery rate . Even applying the most stringent approach described by Bonferroni did not remove the significance for rs12056299 (P = 0.031), rs45324977 (P = 0.019) and haplotype GTCCGC (P = 0.0054) in relation to CHD. On the other hand, we did not consider applying a correction for multiple testing based on the four coronary endpoints. Indeed, such events are highly correlated. Multiple testing is therefore not indicated, because each new test does not provide an independent opportunity for a type-I error . Second, although an observational study cannot prove causation, the Bradford-Hill criteria  suggest that the association between coronary risk and genetic variation in MEOX2 might be causal, taking into account (i) the strength and consistency of the association across different SNPs; (ii) temporality, genetic variability preceding the event; (iii) plausibility based on the experimental studies ; and (iv) the analogy observed with genetic variability in LPL [5, 17–19]. Third, only few genes regulated by MEOX2 and TCF15 are currently known. In this regard, it is important to note that we noticed that of the genes overexpressed in cardiac endothelial cells, those most upregulated included genes involved in lipid homeostasis, including LPL . Finally, as is common in many population studies, follow-up was inconsistent, with varying numbers of follow-up visits across participants. In addition, participants without blood sample, compared with those included in the analyses (Additional file 1: Table S9), were slightly older (49.4 vs. 43.6 years) and had a higher systolic blood pressure (129.5 vs. 125.0 mm Hg), resulting in a higher prevalence of hypertension (40.0 vs. 24.0 %). We cannot ascertain whether these factors might have biased our analyses.
The clinical implications of our current findings can be gauged by the observation that the attributable and population-attributable CHD risks were similar for the MEOX2 GTCCGC carrying state and smoking. Several investigators proposed the use of genetic risk scores based on genome-wide association studies to stratify for the probability of CHD [37, 38]. In the Framingham Heart Study , a score consisting of 13 SNPs did not refine the prediction of CHD or cardiovascular disease, but led to modest improvements in risk reclassification. In contrast, in the Rotterdam Study , a score based on 152 SNPs was associated with incident CHD, but did not enhance risk prediction. SNP discovery based on prevalent rather than incident CHD might explain these discrepancies . In our current study, genetic variation in MEOX2 improved IDI and NRI over and beyond the basic model including traditional CHD risk factors.
Our current study based on a predefined hypothesis generated by data from our experimental studies , identified genetic variation in the transcription factor MEOX2 gene as a novel risk factor for CHD in a white population. However, further experimental studies are required to back-translate our epidemiological observations into underlying molecular mechanisms. Elucidation of these pathways might reveal new targets for the prevention and treatment of CHD.
Coronary heart disease
Cluster of differentiation 36
Encyclopaedia of DNA elements
Growth arrest-specific homeobox
Encyclopaedia of genes and genes variants
Flemish Study on Environment, Genes and Health Outcomes
Genome-wide association study
High density lipoprotein
Integrated discrimination improvement
Mesenchyme homeobox 2
Micro Ribonucleic acid
Nuclear factor of kappa light polypeptide gene enhancer in B cells
Net reclassification improvement
Single nucleotide polymorphism
Transcription factor 15
Aird WC. Phenotypic heterogeneity of the endothelium: I. Structure, function, and mechanisms. Circ Res. 2007;100:158–73.
Stanley WC, Recchia FA, Lopaschuk GD. Myocardial substrate metabolism in the normal and failing heart. Physiol Rev. 2005;85:1093–129.
Coppiello G, Collantes M, Sirerol-Piquer MS, Vandenwijngaert S, Schoors S, Swinnen M, et al. Meox2/Tcf15 heterodimers program the heart capillary endothelium for cardiac fatty acid uptake. Circulation. 2015;131:815–26.
Goldberg IJ. Lipoprotein lipase and lipolysis: central roles in lipoprotein metabolism and atherogenesis. J Lipid Res. 1996;37:693–707.
Tang W, Apostel G, Schreiner PJ, Jacobs Jr DR, Boerwinkle E, Fornage M. Associations of lipoprotein lipase gene polymorphisms with longitudinal plasma lipid trends in yound adults: The Coronary Artery Risk Development in Young Adults (CARDIA) Study. Circ Cardiovasc Genet. 2010;3:179–86.
Treasure CB, Manoukian SV, Klein JL, Vita JA, Nabel EG, Renwick GH, et al. Epicardial coronary artery responses to acetyl choline are impaired in hypertensive patients. Circ Res. 1992;71:776–81.
Celermajer DS, Sorensen KE, Bull C, Robinson J, Deanfield JE. Endothelium-dependent dilation in the systemic arteries of asymptomatic subjects relates to coronary risk factors and their interaction. J Am Coll Cardiol. 1994;24:1468–74.
Staessen JA, Wang JG, Brand E, Barlassina C, Birkenhäger WH, Herrmann SM, et al. Effects of three candidate genes on prevalence and incidence of hypertension in a Caucasian population. J Hypertens. 2001;19:1349–58.
Stolarz-Skrzypek K, Kuznetsova T, Thijs L, Tikhonoff V, Seidlerová J, Richart T, et al. Fatal and nonfatal outcomes, incidence of hypertension and blood pressure changes in relation to urinary sodium excretion in White Europeans. JAMA. 2011;305:1777–85.
Expert Committee on the Diagnosis and Classification of Diabetes Mellitus. Report of the expert committee on the diagnosis and classification of diabetes mellitus. Diabet Care. 2003;26 Suppl 1:5–20.
Citterio L, Simonini M, Zagato L, Salvi E, Delli Carpini S, Lanzani C, et al. Genes involved in vasoconstriction and vasodilation affect salt-sensitive hypertension. PLoS One. 2011;6:e19620.
Liu YP, Gu YM, Thijs L, Knapen MHJ, Salvi E, Citterio L, et al. Inactive matrix gla protein is causally related to adverse health outcomes: a Mendelian randomization study in a Flemish population. Hypertension. 2015;65:463–70.
O'Connell J, Gurdasani D, Delaneau O, Pirastu N, Ulivi S, Cocca M, et al. A general approach for haplotype phasing across the full spectrum of relatedness. PLoS Genet. 2014;10:e1004234.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc B. 1995;57:289–300.
Holtzman NA, Marteau TM. Will genetics revolutionize medicine? N Engl J Med. 2000;343:141–4.
Pencina MJ, D'Agostino Sr RB, D'Agostino Jr RB, Vasan RS. Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond. Stat Med. 2008;27:157–72.
Nordestgaard BG, Albigaard S, Wittrup HH, Steffensen R, Jensen G, Tybjærg-Hansen A. Heterozygous lipoproteinase lipase deficiency. Circulation. 1997;96:1737–44.
Wittrup HH, Tybjærg-Hansen A, Nordestgaard BG. Lipoprotein lipase mutations, plasma lipids and lipoproteins, and the risk of ischemic heart disease. A meta-analysis. Circulation. 1999;99:2901–7.
Li Y, He PP, Zhang DW, Zheng XL, Cayabyab FS, Yin WD, et al. Lipoprotein lipase: from gene to atherosclerosis. Atherosclerosis. 2014;237:597–608.
Cantile M, Schiavo G, Terracciano L, Cillo C. Homeobox genes in normal and abnormal vasculogenesis. Nutr Metab Cardiovasc Dis. 2008;18:651–8.
Skopicki HA, Lyons GE, Schatterman G, Smith RC, Andrés V, Schirm S, et al. Embryonic expression of the Gax homeodomain protein in cardiac, smooth, and skeletal muscle. Circ Res. 1997;80:452–62.
Weir L, Chen D, Pastore C, Isner JM, Walsh K. Expression of gax, a growth arrest homeobox gene, is rapidly down-regulated in the rat carotid artery during the proliferative response to balloon injury. J Biol Chem. 1995;270:5457–61.
Smith RC, Branellec D, Gorski DH, Guo K, Perlman H, Dedieu JF, et al. p21CIP1-mediated inhibition of cell proliferation by overexpression of the gax homeodomain gene. Genes Dev. 1997;11:1674–89.
Chen Y, Rabson AB, Gorski DH. MEOX2 regulates nuclear factor-κB activity in vascular endothelial cells through interactions with p65 and IκBβ. Cardiovasc Res. 2010;87:723–31.
Tabas I, García-Cardeña G, Owens GK. Recent insights into the cellular biology of atherosclerosis. J Cell Biol. 2015;209:13–22.
Chen Y, Gorski DH. Regulation of angiogenesis through a microRNA (miR-130a) that down-regulates antiangiogenic homeobox genes GAX and HOXA5. Blood. 2008;1217–26.
The CARDIoGRAMplusCD4 Consortium. Large-scale association analysis identifies new risk loci for coronary heart disease. Nat Genet. 2013;45:25–33.
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:1393–400.
Mehta NN. Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease. Circ Cardiovasc Genet. 2011;4:327–9.
Saade S, Cazier JBG-SM, Youhanna S, Badro DA, Kamatani Y, Hager J, et al. Large scale association analysis identifies three susceptibility loci for coronary artery disease. PLoS One. 2011;6:e29427.
Della Gaspera B, Armand AS, Lecolle S, Charbonnier F, Chanoine C. Mef2d acts upstream of muscle identity genes and couples lateral myogenesis to dermomyotome formation in Xenopus laevis. PLoS One. 2012;7:e52359.
Davies OR, Lin CY, Radzisheuskaya A, Zhou X, Taube J, Blin G, et al. Tcf15 primes pluripotent cells for differentiation. Cell Rep. 2013;3:472–84.
Paulus WJ, Tschöpe C, Sanderson JE, Rusconi C, Flachskampf FA, Rademakers FE, et al. How to diagnose heart failure: a consensus statement on the diagnosis of heart failure with normal left ventricular ejection fraction by the Heart Failure and Echocardiography Associations of the European Society of Cardiology. Eur Heart J. 2007;28:2539–50.
Zile MR, Brutsaert DL. New concepts in diastolic dysfunction and diastolic heart failure: part I. Diagnosis, prognosis, and measurements of diastolic function. Circulation. 2002;105:1387–93.
Nyholt DR. Genetic case–control association studies - correcting for multiple testing. Hum Genet. 2001;109:564–5.
Hill AB. The environment and disease: association or causation? Proc R Soc Med. 1965;58:295–300.
Thanassoulis G, Peolso 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:113–21.
de Vries PS, Kavousi M, Lightart S, Uitterlinden AG, Hofman A, Franco OH, et al. Incremental predictive value of 152 single nucleotide polymorphisms in the 10-year risk prediction of incident coronary heart disease: the Rotterdam Study. Int J Epidemiol. 2015;44:682–8.
The authors gratefully acknowledge the clerical assistance of Annick De Soete and Renilde Wolfs and the contribution of Linda Custers, Marie-Jeanne Jehoul, Daisy Thijs and Hanne Truyens in data collection at the field centre. The European Union (HEALTH-2011.2.4.2-2-EU-MASCARA, HEALTH-F7-305507 HOMAGE, and the European Research Council Advanced Researcher Grant-2011-294713-EPLORE) and the Fonds voor Wetenschappelijk Onderzoek Vlaanderen, Ministry of the Flemish Community, Brussels, Belgium (G.0881.13 and G.088013) currently support the Studies Coordinating Centre in Leuven. Hypothesis-generating studies at the Centre for Molecular and Vascular Biology were supported by the European Research Council Starting Gran-2007-203291-IMAGINED and the Fonds voor Wetenschappelijk Onderzoek Vlaanderen, Ministry of the Flemish Community, Brussels, Belgium (G.0393.12). The funding source had no role in study design, data extraction, data analysis, data interpretation, or writing of the report. The corresponding author had full access to all the data in the study and had responsibility for the decision to submit for publication.
The authors declare that they have no competing interests.
JAS coordinated the Flemish Study on Environment, Genes and Health Outcomes and managed funding. LT and LJ coordinated the construction and updates of the master database and the management of the biobank. MB, LC, SDC, and ES did the genotyping. CB, DC and PM supervised DNA extraction and amplification, set up the genotyping procedures and managed quality control of genotyping. NC, Y-MG, TP, F-FW, Z-YZ collected phenotypic data and participated in the quality control of the database. AH and TK collected outcome data. W-YY and JAS did the statistical analysis with guidance provided by LT. W-YY and JAS wrote the first draft of the manuscript. XLA, GC, AL and PV assisted in translating the basic science data. All authors interpreted the results, commented on successive drafts of the manuscript and approved the final version.
Table S1. Common tagging SNPs in MEOX2. Table S2. MEOX2 and TCF15 SNPs and allele and genotype frequencies in unrelated founders. Table S3. MEOX2 and TCF15 allele and genotype frequencies in 2027 analysed participants. Table S4. Sex- and age-standardised CHD rates by MEOX2 SNPs. Table S5. Hazard ratios for CHD by MEOX2 SNPs in participants free of CHD at baseline. Table S6. Sex- and age-standardised CHD rates by MEOX2 haplotypes. Table S7. Hazard ratios for CHD by MEOX2 haplotypes reconstructed while accounting for pedigree information. Table S8. Hazard ratios for CHD by MEOX2 haplotypes in participants free of CHD at baseline. Table S9. Baseline characteristics of participants without blood left for genotyping compared with those included in the analyses. Figure S1. Plot of the MEOX2 gene and flanking regions on chromosome 7. Figure S2. Plot of the TCF15 gene and flanking regions on chromosome 20. Figure S3. Interaction between TCF15 rs12624577 and MEOX2 rs4532497. Figure S4. Incidence of coronary endpoints, myocardial infarction and coronary revascularisation in MEOX2 GTCCGC carriers and non-carriers.