Figures
Abstract
Aims
Genetic factors explain a proportion of the inter-individual variation in the risk for atherosclerotic events, but the genetic basis of atherosclerosis and atherothrombosis in families with Mendelian forms of premature atherosclerosis is incompletely understood. We set out to unravel the molecular pathology in a large kindred with an autosomal dominant inherited form of premature atherosclerosis.
Methods and Results
Parametric linkage analysis was performed in a pedigree comprising 4 generations, of which a total of 11 members suffered from premature vascular events. A parametric LOD-score of 3.31 was observed for a 4.4 Mb interval on chromosome 12. Upon sequencing, a non-synonymous variant in KERA (c.920C>G; p.Ser307Cys) was identified. The variant was absent from nearly 28,000 individuals, including 2,571 patients with premature atherosclerosis. KERA, a proteoglycan protein, was expressed in lipid-rich areas of human atherosclerotic lesions, but not in healthy arterial specimens. Moreover, KERA expression in plaques was significantly associated with plaque size in a carotid-collar Apoe−/− mice (r2 = 0.69; p<0.0001).
Citation: Maiwald S, Sivapalaratnam S, Motazacker MM, van Capelleveen JC, Bot I, de Jager SC, et al. (2014) Mutation in KERA Identified by Linkage Analysis and Targeted Resequencing in a Pedigree with Premature Atherosclerosis. PLoS ONE 9(5): e98289. https://doi.org/10.1371/journal.pone.0098289
Editor: Florian Kronenberg, Innsbruck Medical University, Austria
Received: March 12, 2014; Accepted: April 30, 2014; Published: May 30, 2014
Copyright: © 2014 Maiwald et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. Data are with the corresponding author.
Funding: SM has been supported by the Marie-Curie NetSim ITN grant [grant number EC-215820] from the European Commission 7th Framework Programme and Ipse Movet. SS is supported by a grant from ATHEROS. MMM is supported by Fondation LeDucq (Transatlantic Network, 2009-2014). MvE is an Established Investigator [2007T056] of the Dutch Heart Foundation. JJPK is a recipient of the Lifetime Achievement Award of the Dutch Heart Foundation [2010T082]. Research in WHO group is supported by programme grants from the National Institute for Health Research (NIHR) (RP-PG-0310-1002), from the British Heart Foundation [RG/09/12/28096], from the European Commission and the Medical Research Council. The project made use of NHS Blood and Transplant donors from the Cambridge BioResource (http://www.cambridgebioresource.org.uk/). This local resource for genotype-phenotype association studies is supported by a grant from the NIHR to the Cambridge Biomedical Research Centre. GKH is holder of a Veni grant [91612122] from the Netherlands Organisation for Scientific Research (NWO). This work is supported by CardioVascular Research Initiative [CVON2011-19; Genius] and the European Union [Resolve: FP7-305707’ TransCard: FP7-603091-2]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interest exist.
Introduction
In both cardiovascular disease (CVD) and stroke, atherosclerosis is the underlying pathology. Genetic factors explain a proportion of the observed inter-individual variability in atherosclerosis progression, which is exemplified by the observed 30–60% heritability in twin studies [1], and the finding that a positive family history for premature atherosclerosis is an independent risk factor [2]. Both common and rare genetic variants contribute to the heritability [3]. A recent meta-analysis of Genome Wide Association Studies (GWAS) of nearly 64,000 cases with CVD has identified 46 common single nucleotide polymorphisms (SNPs) of small effect size, which account for about 10.6% of the estimated heritability [4]. The remaining heritability is assumed to be explained by a combination of common variants with effect sizes so small that they remained undetected in the recent GWAS meta-analysis, by rare variants with an intermediate effect, and by pedigree-specific mutations with a large effect. The latter have been identified in several pedigrees with Mendelian forms of atherosclerosis [5]–[11]. A well-known example of such a monogenic dominant disorder, that underlies atherosclerosis, is Familial Hypercholesterolemia (FH), caused by loss of function (LOF) causing mutations in the genes encoding for the Low-density lipoprotein receptor (LDLR) or Apolipoprotein B (APOB), or gain of function (GOF) mutations in Proprotein convertase subtilisin/kexin type 9 (PCSK9). Carriers of mutations in these genes are characterized by high plasma levels of LDL-cholesterol (LDL-c) and early onset atherosclerosis [12].
The molecular basis of premature atherosclerosis in the absence of high LDL-c levels is largely unknown, but recently rare and putative causal variants in Myocyte enhancer factor–2 (MEF2A) and Low-density lipoprotein receptor-related protein 6 (LRP6) have been identified in pedigrees with Mendelian forms of atherosclerosis where FH as a causal factor was ruled out [5], [10], [11]. The predictive power of a family based approach has recently been documented by Erdmann and co-workers [13], who identified 2 novel private mutations in Guanylate Cyclase 1 soluble alpha 3 (GUCY1A3) and Chaperone Containing TCP1 subunit 7 (CCT7), in a pedigree with a mendelian form of CVD. The identification of such mutations in novel genes provides new and pivotal information about the pathobiology of premature atherosclerosis, and may ultimately lead to new pharmacological interventions to address the burden of atherosclerosis.
The aim of the current study was to identify the molecular defect in a large non-FH pedigree with an autosomal dominant form of atherosclerosis. We identified a non-synonymous (ns) mutation in the Keratocan (KERA) gene, which encodes the extracellular proteoglycan KERA. Additional genetic, histological and animal studies were performed to further establish the role of this variant in atherosclerosis.
Methods
Recruitment of the Pedigree with Early Onset Atherosclerosis
A male subject who suffered from an acute myocardial infarction (AMI) at the age of 49 years was referred to the outpatient clinic of the Academic Medical Center (AMC) in Amsterdam (Figure 1B; index case III:8). An autosomal dominant form of inheritance of early onset atherosclerosis in the pedigree was identified. Premature atherosclerosis was defined as a documented atherosclerotic event, either CVD or stroke, before the age of 55 (male) and 65 (female). Blood was obtained after an overnight fast. Plasma was isolated by centrifugation at 3000 rpm, for 20 minutes at 4°C and was stored at −80°C for further analyses. Plasma cholesterol, LDL-c, high-density lipoprotein cholesterol (HDL-c) and triglycerides (TG) were analysed using commercially available assays (Randox, Antrim, United Kingdom and Wako, Neuss, Germany) on a Cobas-Mira autoanalyzer (Roche, Basel, Switzerland). Hypertension was defined as a systolic blood pressure >140 mmHg and/or diastolic blood pressure >90 mmHg or the use of anti-hypertensive lowering drugs. Diabetes mellitus was defined as fasting plasma glucose ≥7.0 mmol/l or use of anti-diabetic medication. The study complies with the Declaration of Helsinki and the Institutional Review Board of the AMC (Medische Ethische Toetsings Commissie, METC) of the University of Amsterdam approved the study. All participants provided written informed consent.
Genetic Studies (Figure 2)
Genomic DNA was extracted from whole blood on an AutopureLS apparatus according to the manufacturer’s protocol (Gentra Systems, Minneapolis, MN, USA). Human CytoSNP-12 BeadChip kits were used for genome wide single nucleotide polymorphism (SNP) genotyping (Illumina, San Diego, CA, USA) in 12 relatives (Figure 1B; 9 affected and 3 unaffected). A Nimblegen (Madison, WI, USA) custom sequence capture array comprising 395K probes was designed to enrich for the genomic region that was identified by linkage analysis and used to sequence the DNA region with an Illumina GAII platform. Confirmation of the identified mutations and analysis of co-segregation of the variant in the pedigree was by Sanger sequencing as previously described [14]. The following primer pairs were used: KERA: exon 1 forward 5′-AAG ATT ACC AGC CAA TAC AAT GC-3′, reverse 5′-TGA TGG GAG ACC CTC ATC TG-3′ exon 2 forward: 5′-GCC ACT AAG CCC TCC ATA GG-3′; reverse (2) 5′-AGC AAT GGG GAA TAT GAC TTG-3′. After establishing the segregation in the core pedigree, the family was further expanded (Figure S1 in File S1).
A, A significant LOD-score of 3.31 was obtained on chromosome 12 after parametric linkage analysis by Allegro software using genotypes from 12 family members. B, Haplotype analyses of the linkage interval showing two recombination events in III-3 and III-8; top to bottom SNP order: rs11104542, rs7974491, rs10506959, rs1688545, rs1347846, rs704106, rs1948839, rs704144, rs790455, rs7980716, rs10777477, rs1493848, rs12308959, rs3847810, rs12426730. The linkage interval is defined by rs1688545 and rs1493848 and harbours 21 genes (Tables S2 and S3) Additional sequencing of the region led to the identification of the KERA p.Ser307Cys variant. This was the only rare non-synonymous variant. All affected relatives in this pedigree are heterozygous carriers of the KERA mutation. In the pedigree we show the type of event and the age at which the event occurred for each relative. III-8 is the index case. AMI = acute myocardial infarction; TIA = transient ischemic attack; PTCA = percutaneous transluminal coronary angioplasty; CVA = cerebrovascular accident; AP = angina pectoris; ACS = acute coronary syndrome; CABG = coronary artery bypass graft; C, DNA Sanger sequencing chromatogram showing the heterozygote c.920C>G; p.Ser307Cys KERA mutation.
Validation Cohorts
The KERA variant was genotyped in:
- Premature Atherosclerosis (PAS) Cohort: this cohort comprises 935 patients with early symptomatic atherosclerosis before the age of 51 years. Atherosclerosis is defined as myocardial infarction, coronary revascularization, or evidence of at least 70% stenosis in a major epicardial artery. [15] Patients were recruited at the cardiology and vascular outpatient clinic of the AMC. [16] To identify possible further cases with mutation or rare variants in the coding fraction of KERA the DNA samples of 296 randomly chosen PAS cases were sequenced.
- Sanquin Blood Bank common Controls: DNA samples from 1,440 healthy volunteers were recruited from a large cohort of healthy blood donors, who were free of clinical CVD, at one of the collection clinics of the Sanquin Blood Bank covering the northwest section of the Netherlands, which geographically overlaps the PAS case cohort [16].
- Cambridge BioResource Collection: DNA samples of 8,946 healthy volunteers were enrolled by NHS Blood and Transplant Unit in a resource for genotype-phenotype association studies [17]. In addition, genotyping results from 16,515 samples were retrieved from the UK10K (http://www.uk10k.org) and the NHLBI Exome Sequencing (ESP) projects [18].
Human Plaque Quantification
Specimens of cornea (n = 2), tonsil (n = 2), mammary artery (n = 2), atherosclerotic- and non-diseased arteries (n = 9) were collected from patients at autopsy. The tissues from autopsy material for histological verification of protein expression were obtained from the Department of Pathology (Prof Dr AC van der Wal) within the AMC, Amsterdam. They were anonymously provided to us according to the GPC guidelines. No METC conformation was required. For details of the staining procedures see Supplemental Methods in File S1.
Molecular Dynamics Computer Simulations
Molecular dynamics (MD) computer stimulations were performed to assess the effect of the p.Ser307Cys mutation on the structure of the KERA protein (see Supplemental Methods in File S1).
Animal Experiments
The animal protocol was approved by the Ethics Committee for animal Experiments of the Leiden University (Leiden, The Netherlands) and carried out in compliance with the Dutch government guidelines. The animals were bred in house in the Gorleaus Laboratories of the Leiden/Amsterdam Center for Drug Research Leiden, the Netherlands. All surgery was performed under sodium pentobarbital anesthesia, and all efforts were made to minimize suffering. Male Apoe−/− mice, aged 10–12 weeks, were fed a Western type diet containing 0.25% cholesterol and 15% cacao butter (Special Diet Services, Sussex, UK) starting two weeks prior to collar placement surgery. Mice (n = 3–4 per group) were anaesthetized at t = 0 or 2, 4, 6, 8 and 10 weeks after collar placement and in-situ perfusion-fixation was performed, after which carotid arteries were sectioned and lesions were analyzed. [19] Morphometric analysis using Leica Qwin image analysis software (Leica Microsystems, Rijswijk, The Netherlands) was performed on hematoxylin and eosin stained 5 µm sections of the carotid arteries at the site of maximal stenosis as previously described. [19] Cryosections were stained for KERA protein using a rabbit anti KERA polyclonal antibody (clone H-50, Sc-66941) (Supplemental Methods in File S1).
Statistical Analysis
Results are expressed as mean ± standard deviation unless otherwise stated. Two-sided P values of <0.05 were considered statistically significant. All statistical analysis were performed with SPSS version 18.0. Multipoint parametric linkage analysis assuming a fully penetrant autosomal dominant model, based on the clinical segregation of the disease phenotype in the pedigree and minor allele frequency (MAF) of <0.001 was conducted using Allegro software [20]. The identified variants were compared with the dbSNP (build 129) database and available data from the 1000 Genomes pilot project [21]. Only those variants that were not present in the database were selected for further testing.
Results
Identification of the Mutation in KERA
The index case III.8 (Figure 1B) suffered from an AMI at the age of 49 years, and apart from obesity, no other atherosclerosis risk factors were identified (Tables 1 and 2). Ten relatives did suffer from premature atherosclerosis, of whom 2 died prior to the conduct of this study. In addition to the nine living affected cases, we identified three unaffected relatives, who were above the predefined age for PAS and had not suffered from an atherosclerotic event. (Figure 1B, Tables 1 and 2). Genotyping was performed in these 12 individuals and multi-point linkage analysis resulted in a parametric LOD-score of 3.31 in a region located on chromosome 12q21.33-q22 (Figure 1A). The interval of 4.4 Mb was located between SNPs rs1688545 and rs1493848 (Figure 1B) and harboured 21 annotated genes (Table S1 in File S1 and Table S2 in File S1). Sequencing of the interval in the index case III.8 revealed 1 non-synonymous variant, which was absent from the 1000 Genomes project. This novel variant in exon 3 (NM_007035.3: c.920C>G; p.Ser307Cys, Figure 1C) was identified in a single KERA haplotype, which showed complete co-segregation with premature atherosclerosis in the core pedigree. Further expansion of the pedigree led to the identification of 10 additional KERA mutation carriers of whom 1 brother, aged 73 years, had no signs of atherosclerotic disease. Two female carriers of the mutation are still under the PAS age (Figure S1in File S1 and Table S3 in File S1). The KERA p.Ser307Cys variant was absent from 27,901 other DNA samples, including samples from 935 cases from the PAS cohort (Table S4 in File S1) and 1,636 cases of premature MI collection of the NHLBI Exome Sequencing project. The coding sequence and intron-exon boundaries of KERA were additionally sequenced in 296 PAS subjects and no additional rare variants with MAFs<0.5% were identified.
Autosomal recessive LOF mutations in KERA cause cornea plana type 2 (CP2), an ophthalmologic disorder characterised by corneal flattening, but split lamp examination in 2 carriers of the KERA p.Ser307Cys mutation did not reveal CP2 characteristics. Additionally, in 9 CP2 cases and their 9 first-degree relatives (age>55 and >65 years for males and females) no evidence of clinical events related with atherosclerosis were detected.
The Effect of the Mutation on the KERA Protein Structure; In silico Analysis (Figure 3)
The cysteine residues at 303 and 343 in KERA are highly conserved among a large range of animal species. The mutation, which introduces a novel cysteine at residue 307 leads to a substantial change in this conserved region. This was studied in more detail by molecular dynamics (MD) computer simulations started from a homology model based on the crystal structure of decorin. [22] The Ser307Cys mutation is localised in the C-terminal part of KERA and is flanked by cysteine residues at positions 303 and 343, which are assumed to be covalently bonded (Figure 3A). Whilst the horseshoe fold, typical of leucine-rich repeat domain containing proteins, was maintained during simulation, the mutation at residue 307 may form a preferred disulphide bond with Cys303, resulting in Cys343 being unpaired and to become more solvent exposed than in the wild-type KERA protein (Figures 3A–D).
Molecular dynamic simulations were performed as described in File S1. The residue Cys303 is highlighted in yellow, Cys343 in green, Ser307 in cyan and Cys307 residue in blue. A, Structure of wild-type KERA. C Structure of the KERA mutant p.Ser307Cys. B, A detailed view of the C-terminal part of wild type KERA highlighting the Cys303–Cys343 disulphide bond. D, Possible structural effects of the substitution of a serine for a cysteine at residue 307 showing a favourable Cys303–Cys307 disulphide bond. Consequently, Cys343 is available for binding with cysteine residues of other proteins.
KERA is Associated with Atherosclerotic Burden in Humans
KERA was not expressed in tonsil tissue (Figure 4A) and mammary artery segments (Figure 4C). KERA was nicely expressed in cornea (Figure 4B). Interestingly KERA was present in a human artery segment obtained from a patient who is a carrier of the KERA mutation (Figure 4D). Additionally we have tested arterial segments from 9 patients at different stages of atherosclerosis progression (Figures S2B and S2C in File S1) and found that KERA was not expressed in 2 non-diseased arterial segments (Figure S2A in File S1) but was expressed in atherosclerotic plaque regions. KERA was localised within the extracellular matrix, in close vicinity to the lipid core of the atherosclerotic plaque. Domains rich in KERA showed a strong staining by antibodies against CD3 and Chemokine (C-X-C motif) ligand 1 (CXCL1), indicating the presence of lymphoid T helper (Th1) cells in addition to the myeloid cells (Figures 5A–C).
A, KERA is not expressed in tonsil tissue; B: KERA is expressed in corneal tissue C: KERA is not expressed in mammary artery. D: KERA is expressed in arterial segments obtained from a patient with the KERA mutation.
A: triple stain with KERA (blue), CXCL1 (red) and CD3 positive cells (brown). B: Spectral Imaging of triple staining. C: The yellow staining demonstrates that these cells are positive for all three components.
KERA Expression is Associated with Atherosclerotic Burden in Apoe−/− Mice
Next we tested whether KERA protein expression was associated with atherosclerotic burden in an established Apoe−/− mouse model for atherosclerosis, in which atherosclerosis was induced by perivascular collar placement. At different time points after collar placement the degree of atherosclerosis and KERA protein expression was analysed in carotid artery. KERA was detected in both early and advanced atherosclerotic lesions (Figures 6A and B). Interestingly, the extent of atherosclerotic lesion formation, quantified by intima area, was over time significantly correlated with KERA protein expression (r2 = 0.69, P<0.0001) (Figures 6C and D).
A–B, Early (week = 2, A) and advanced (week = 8, B) atherosclerotic tissue from murine carotid arteries were stained for KERA (brown) and hematoxylin (blue). While present mainly near endothelial cells in early lesion, KERA is predominantly present in the matrix of the plaque at more advanced lesions. C–D, KERA expression overtime in Apoe−/− mice with collar placement show significant correlation with plaque size (r2 = 0.69, P<0.0001).
Discussion
We identified a novel p.Ser307Cys mutation in the extracellular matrix protein KERA in a large pedigree with a Mendelian form of premature atherosclerosis by linkage analysis combined with next generation sequencing. The variant was absent from nearly 28,000 DNA samples. KERA encodes for a keratan sulfate proteoglycan expressed in the cornea, trachea and at lower levels in the intestine, skeletal muscle, ovary and lung, However no expression in the vasculature has been reported [23]. Interestingly, the protein was absent from healthy artery segments, but heavily abundant within the lipid core of atherosclerotic lesions, which emphasizes that KERA might be a novel actor in the pathobiology underlying atherosclerosis. This notion is further strengthened by the fact that in Apoe−/− mice, induction of atherosclerosis in the carotid arteries by cuff placements significantly correlated with KERA expression in the plaques.
The mutation introduces a cysteine at residue 307 in the carboxy-terminal of the gene and is flanked by similar residues at positions 303 and 343, which are assumed to be disulphide bonded in wild-type KERA. Molecular dynamic simulations suggested that the new Cys307 favours form a disulphide bond with Cys303, which may lead to improved stability of the protein fold and structure around this carboxy-terminal domain. In addition, during simulation Cys343 was found to be more solvent exposed in the mutant protein thus enabling novel protein-protein interactions. Collectively, the results from the molecular dynamic simulations favour a gain of function (GOF) for the p.Ser307Cys mutation.
The concept of GOF is further substantiated by studies in patients with LOF mutations in KERA. Thus far, autosomal recessive LOF mutations in KERA cause cornea plane type 2 (CP2) (OMIM 121400; 217300), a rare disorder characterized by excessively flat and thin corneas. Worldwide about 100 cases have been described [24]–[28]. An increased risk for atherosclerosis has not been reported in these individuals, but this might have been overlooked because previous studies may have focussed solely on the ophthalmological consequences of KERA mutations. However, in the 9 identified CP2 patients in the Netherlands, and in their first degree relatives, no symptoms of premature atherosclerosis were observed, suggesting that both partial or complete LOF for KERA does not confer a substantial increase in the risk of atherosclerosis. No corneal abnormalities were identified in 2 affected family members (IV:I and IV:2, Figure 2B) carrying the novel p.Ser305Cys mutation. The absence of opthalmic effects in the assumed GOF mutation carriers is in line with previous studies, showing that KERA mutations are not identified in subjects suffering from cornea plane type 1 (CP1), the autosomal dominant form of cornea plana [29].
The clinical diversity between heterozygous carriers of our novel GOF KERA variant, characterized by atherosclerosis risk, and homozygous carriers of LOF KERA variants, resulting in CP2, might be related to the putative GOF effect of the novel mutation. The role of KERA in atherosclerosis has not been investigated so far, but it is interesting to note that KERA has been shown to play a role in neutrophil migration [30]. In mice, neutrophil migration is orchestrated by a chemical gradient of a wide range of chemokines including CXCL1 in the vessel wall [31]. KERA is one of the regulators of the CXCL1 gradient [32]. Interestingly, endothelial CXCL1 has recently been shown to play a crucial role in hyperlipidemia-induced arterial leukocyte arrest [33]. The receptor for CXCL1 CXCR2 is present on myeloid cells like neutrophils and monocyte/macrophages, which are directly involved in all aspects of atherosclerosis [30]. We hypothesize that the GOF mutation observed in the large PAS pedigree results in an augmented binding of KERA to CXCL1, which may lead to increased neutrophil migration into the vessel wall. This notion is supported by recent observations in mice, where both arterial CXCL1 and leukocyte-specific CXCR2 expression are central to macrophage accumulation in established fatty streak lesions [34]. Concomitantly in mice lacking Cxcl1, atherosclerosis is significantly reduced [35]. Notably, in human artery segments we observe a co-localisation of KERA with CXCL1 and the lymphoid T cell marker CD3. Further functional studies are required to confirm our proposed model.
A number of considerations have to be taken into account while interpreting the data. Although the results from this study are suggestive of a role for KERA in atherosclerosis, a direct causative role has not been established thus far.
A specific concern is that neither we, nor others did observe possible KERA GOF mutations in other cases with premature atherosclerosis. However, the identification of this extremely rare variant, which is absent from nearly 56,000 alleles and the confirmation of the presence of KERA in the plaque in mice and men does suggest that KERA may be an active player in an, as of yet, not fully elucidated novel pathway in atherosclerosis. Further studies are warranted to confirm our findings and to establish whether KERA might be an attractive target for therapy.
Supporting Information
File S1.
This file contains a more extended method sextoin is presented including detailed methodologic information and Table S1–Table S4, Figure S1–Figure S2, and a Reference list (References S1). Table S1, Annotated Genes in linkage interval on chromosome 12. Table S2, SNPs in linkage interval on chromosome 12. Table S3, Demographics of additional relatives of the extensive pedigree. Table S4, Age, BMI and plasma lipids for participants in the PAS cohort (N = 935). Figure S1, Extended pedigree. Figure S2, KERA expression in atherosclerotic plaque segments. References S1.
https://doi.org/10.1371/journal.pone.0098289.s001
(DOCX)
Acknowledgments
We would like to acknowledge all participants. We thank Christy Holtkamp for the genetic fieldwork and Jitske Kuipers, Caroline Aalbers, Sylvia Nurnberg and Jorge Peter for technical and laboratory assistance. Also we would like to acknowledge the Cornea Workgroup of the Dutch Ophthalmic Society for providing us with cornea tissue material. We gratefully acknowledge the participation of all NIHR Cambridge BioResource volunteers. We thank the Cambridge BioResource staff for their help with volunteer recruitment.
Author Contributions
Conceived and designed the experiments: MDT WHO GKH GDT ME JK DPG. Performed the experiments: SM SS MMM JCC IB SCJ JJ JS CAA CRV HK ACW CML. Analyzed the data: SM MMM GDT ME JK. Contributed reagents/materials/analysis tools: SM MMM DPG GDT. Wrote the paper: SM MM GKH GDT MDT WHO JJPK JK.
References
- 1. Lee K, Sung J, Lee SC, Park SW, Kim YS, et al. (2012) Segment-specific carotid intima-media thickness and cardiovascular risk factors in Koreans: the Healthy Twin Study. Eur J Prev Cardiol 19: 1161–1172.
- 2. Sivapalaratnam S, Boekholdt SM, Trip MD, Sandhu MS, Luben R, et al. (2010) Family history of premature coronary heart disease and risk prediction in the EPIC-Norfolk prospective population study. Heart 96: 1985–1989.
- 3. Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, et al. (2009) Finding the missing heritability of complex diseases. Nature 461: 747–753.
- 4. Deloukas P, Kanoni S, Willenborg C, Farrall M, Assimes TL, et al. (2013) Large-scale association analysis identifies new risk loci for coronary artery disease. Nat Genet 45: 25–33.
- 5. Mani A, Radhakrishnan J, Wang H, Mani A, Mani MA, et al. (2007) LRP6 Mutation in a Family with Early Coronary Disease and Metabolic Risk Factors. Science 315: 1278–1282.
- 6. Wang Q, Rao S, Shen GQ, Li L, Moliterno DJ, et al. (2004) Premature myocardial infarction novel susceptibility locus on chromosome 1P34–36 identified by genomewide linkage analysis. Am J Hum Genet 74: 262–271.
- 7. Guella I, Rimoldi V, Asselta R, Ardissino D, Francolini M, et al. (2009) Association and functional analyses of MEF2A as a susceptibility gene for premature myocardial infarction and coronary artery disease. Circ Cardiovasc Genet 2: 165–172.
- 8. Keramati AR, Singh R, Lin A, Faramarzi S, Ye ZJ, et al. (2011) Wild-type LRP6 inhibits, whereas atherosclerosis-linked LRP6R611C increases PDGF-dependent vascular smooth muscle cell proliferation. Proc Natl Acad Sci USA 108: 1914–1918.
- 9. Lieb W, Mayer B, Konig IR, Borwitzky I, Gotz A, et al. (2008) Lack of association between the MEF2A gene and myocardial infarction. Circulation 117: 185–191.
- 10. Liu Y, Niu W, Wu Z, Su X, Chen Q, et al. (2012) Variants in exon 11 of MEF2A gene and coronary artery disease: evidence from a case-control study, systematic review, and meta-analysis. PLoS One 7: e31406.
- 11. Wang L, Fan C, Topol SE, Topol EJ, Wang Q (2003) Mutation of MEF2A in an Inherited Disorder with Features of Coronary Artery Disease. Science 302: 1578–1581.
- 12. Huijgen R, Vissers MN, Defesche JC, Lansberg PJ, Kastelein JJ, et al. (2008) Familial hypercholesterolemia: current treatment and advances in management. Expert Rev Cardiovasc Ther 6: 567–581.
- 13. Erdmann J, Stark K, Esslinger UB, Rumpf PM, Koesling D, et al. (2013) Dysfunctional nitric oxide signalling increases risk of myocardial infarction. Nature 504: 432–436.
- 14. Surendran RP, Visser ME, Heemelaar S, Wang J, Peter J, et al. (2012) Mutations in LPL, APOC2, APOA5, GPIHBP1 and LMF1 in patients with severe hypertriglyceridaemia. J Intern Med 272: 185–196.
- 15. Hauser ER, Mooser V, Crossman DC, Haines JL, Jones CH, et al. (2003) Design of the Genetics of Early Onset Cardiovascular Disease (GENECARD) study. Am Heart J 145: 602–613.
- 16. Vergeer M, Boekholdt SM, Sandhu MS, Ricketts SL, Wareham NJ, et al. (2010) Genetic variation at the phospholipid transfer protein locus affects its activity and high-density lipoprotein size and is a novel marker of cardiovascular disease susceptibility. Circulation 122: 470–477.
- 17. Dendrou CA, Plagnol V, Fung E, Yang JH, Downes K, et al. (2009) Cell-specific protein phenotypes for the autoimmune locus IL2RA using a genotype-selectable human bioresource. Nat Genet 41: 1011–1015.
- 18. Fu W, O’Connor TD, Jun G, Kang HM, Abecasis G, et al. (2013) Analysis of 6,515 exomes reveals the recent origin of most human protein-coding variants. Nature 493: 216–220.
- 19. von der Thusen JH, Van Berkel TJ, Biessen EA (2001) Induction of rapid atherogenesis by perivascular carotid collar placement in apolipoprotein E-deficient and low-density lipoprotein receptor-deficient mice. Circulation 103: 1164–1170.
- 20. Gudbjartsson DF, Bjornsdottir US, Halapi E, Helgadottir A, Sulem P, et al. (2009) Sequence variants affecting eosinophil numbers associate with asthma and myocardial infarction. Nat Genet 41: 342–347.
- 21. A map of human genome variation from population-scale sequencing. Nature 467: 1061–1073.
- 22. Scott PG, McEwan PA, Dodd CM, Bergmann EM, Bishop PN, et al. (2004) Crystal structure of the dimeric protein core of decorin, the archetypal small leucine-rich repeat proteoglycan. Proc Natl Acad Sci USA 101: 15633–15638.
- 23. Wentz-Hunter K, Cheng EL, Ueda J, Sugar J, Yue BY (2001) Keratocan expression is increased in the stroma of keratoconus corneas. Mol Med 7: 470–477.
- 24. Forsius H, Damsten M, Eriksson AW, Fellman J, Lindh S, et al. (1998) Autosomal recessive cornea plana. A clinical and genetic study of 78 cases in Finland. Acta Ophthalmol Scand 76: 196–203.
- 25. Khan A, Al-Saif A, Kambouris M (2004) A novel KERA mutation associated with autosomal recessive cornea plana. Ophthalmic Genet 25: 147–152.
- 26. Liskova P, Hysi PG, Williams D, Ainsworth JR, Shah S, et al. (2007) Study of p.N247S KERA mutation in a British family with cornea plana. Mol Vis 13: 1339–1347.
- 27. Khan AO, Aldahmesh M, Meyer B (2006) Recessive cornea plana in the Kingdom of Saudi Arabia. Ophthalmology 113: 1773–1778.
- 28.
Dudakova L, Palos M, Hardcastle AJ, Liskova P (2013) Corneal Endothelial Findings in a Czech Patient with Compound Heterozygous Mutations in KERA. Ophthalmic Genet: doi:https://doi.org/10.3109/13816810.2013.811272.
- 29. Aldave AJ, Sonmez B, Bourla N, Schultz G, Papp JC, et al. (2007) Autosomal dominant cornea plana is not associated with pathogenic mutations in DCN, DSPG3, FOXC1, KERA, LUM, or PITX2. Ophthalmic Genet 28: 57–67.
- 30. Carlson EC, Sun Y, Auletta J, Kao WW, Liu CY, et al. (2010) Regulation of corneal inflammation by neutrophil-dependent cleavage of keratan sulfate proteoglycans as a model for breakdown of the chemokine gradient. J Leukoc Biol 88: 517–522.
- 31. Zernecke A, Shagdarsuren E, Weber C (2008) Chemokines in atherosclerosis: an update. Arterioscler Thromb Vasc Biol 28: 1897–1908.
- 32. Carlson EC, Lin M, Liu CY, Kao WW, Perez VL, et al. (2007) Keratocan and lumican regulate neutrophil infiltration and corneal clarity in lipopolysaccharide-induced keratitis by direct interaction with CXCL1. J Biol Chem 282: 35502–35509.
- 33. Zhou Z, Subramanian P, Sevilmis G, Globke B, Soehnlein O, et al. (2011) Lipoprotein-derived lysophosphatidic acid promotes atherosclerosis by releasing CXCL1 from the endothelium. Cell Metab 13: 592–600.
- 34. Boisvert WA, Rose DM, Johnson KA, Fuentes ME, Lira SA, et al. (2006) Up-regulated expression of the CXCR2 ligand KC/GRO-alpha in atherosclerotic lesions plays a central role in macrophage accumulation and lesion progression. Am J Pathol 168: 1385–1395.
- 35. Boisvert WA, Santiago R, Curtiss LK, Terkeltaub RA (1998) A leukocyte homologue of the IL-8 receptor CXCR-2 mediates the accumulation of macrophages in atherosclerotic lesions of LDL receptor-deficient mice. J Clin Invest 101: 353–363.