Genome Analyses of >200,000 Individuals Identify 58 Loci for Chronic Inflammation and Highlight Pathways that Link Inflammation and Complex Disorders.

C-reactive protein (CRP) is a sensitive biomarker of chronic low-grade inflammation and is associated with multiple complex diseases. The genetic determinants of chronic inflammation remain largely unknown, and the causal role of CRP in several clinical outcomes is debated. We performed two genome-wide association studies (GWASs), on HapMap and 1000 Genomes imputed data, of circulating amounts of CRP by using data from 88 studies comprising 204,402 European individuals. Additionally, we performed in silico functional analyses and Mendelian randomization analyses with several clinical outcomes. The GWAS meta-analyses of CRP revealed 58 distinct genetic loci (p < 5 × 10-8). After adjustment for body mass index in the regression analysis, the associations at all except three loci remained. The lead variants at the distinct loci explained up to 7.0% of the variance in circulating amounts of CRP. We identified 66 gene sets that were organized in two substantially correlated clusters, one mainly composed of immune pathways and the other characterized by metabolic pathways in the liver. Mendelian randomization analyses revealed a causal protective effect of CRP on schizophrenia and a risk-increasing effect on bipolar disorder. Our findings provide further insights into the biology of inflammation and could lead to interventions for treating inflammation and its clinical consequences.

Chromosome X-wide association study identifies Loci for fasting insulin and height and evidence for incomplete dosage compensation.

The X chromosome (chrX) represents one potential source for the "missing heritability" for complex phenotypes, which thus far has remained underanalyzed in genome-wide association studies (GWAS). Here we demonstrate the benefits of including chrX in GWAS by assessing the contribution of 404,862 chrX SNPs to levels of twelve commonly studied cardiometabolic and anthropometric traits in 19,697 Finnish and Swedish individuals with replication data on 5,032 additional Finns. By using a linear mixed model, we estimate that on average 2.6% of the additive genetic variance in these twelve traits is attributable to chrX, this being in proportion to the number of SNPs in the chromosome. In a chrX-wide association analysis, we identify three novel loci: two for height (rs182838724 near FGF16/ATRX/MAGT1, joint P-value = 2.71×10(-9), and rs1751138 near ITM2A, P-value = 3.03×10(-10)) and one for fasting insulin (rs139163435 in Xq23, P-value = 5.18×10(-9)). Further, we find that effect sizes for variants near ITM2A, a gene implicated in cartilage development, show evidence for a lack of dosage compensation. This observation is further supported by a sex-difference in ITM2A expression in whole blood (P-value = 0.00251), and is also in agreement with a previous report showing ITM2A escapes from X chromosome inactivation (XCI) in the majority of women. Hence, our results show one of the first links between phenotypic variation in a population sample and an XCI-escaping locus and pinpoint ITM2A as a potential contributor to the sexual dimorphism in height. In conclusion, our study provides a clear motivation for including chrX in large-scale genetic studies of complex diseases and traits.

Plasma ceramides predict cardiovascular death in patients with stable coronary artery disease and acute coronary syndromes beyond LDL-cholesterol.

AIMS: The aim was to study the prognostic value of plasma ceramides (Cer) as cardiovascular death (CV death) markers in three independent coronary artery disease (CAD) cohorts. METHODS AND RESULTS: Corogene study is a prospective Finnish cohort including stable CAD patients (n = 160). Multiple lipid biomarkers and C-reactive protein were measured in addition to plasma Cer(d18:1/16:0), Cer(d18:1/18:0), Cer(d18:1/24:0), and Cer(d18:1/24:1). Subsequently, the association between high-risk ceramides and CV mortality was investigated in the prospective Special Program University Medicine-Inflammation in Acute Coronary Syndromes (SPUM-ACS) cohort (n = 1637), conducted in four Swiss university hospitals. Finally, the results were validated in Bergen Coronary Angiography Cohort (BECAC), a prospective Norwegian cohort study of stable CAD patients. Ceramides, especially when used in ratios, were significantly associated with CV death in all studies, independent of other lipid markers and C-reactive protein. Adjusted odds ratios per standard deviation for the Cer(d18:1/16:0)/Cer(d18:1/24:0) ratio were 4.49 (95% CI, 2.24-8.98), 1.64 (1.29-2.08), and 1.77 (1.41-2.23) in the Corogene, SPUM-ACS, and BECAC studies, respectively. The Cer(d18:1/16:0)/Cer(d18:1/24:0) ratio improved the predictive value of the GRACE score (net reclassification improvement, NRI = 0.17 and ΔAUC = 0.09) in ACS and the predictive value of the Marschner score in stable CAD (NRI = 0.15 and ΔAUC = 0.02). CONCLUSIONS: Distinct plasma ceramide ratios are significant predictors of CV death both in patients with stable CAD and ACS, over and above currently used lipid markers. This may improve the identification of high-risk patients in need of more aggressive therapeutic interventions.

Genetic variation on the BAT1-NFKBIL1-LTA region of major histocompatibility complex class III associates with periodontitis.

Periodontitis is a chronic inflammatory disease with a multifactorial etiology. We investigated whether human major histocompatibility complex (MHC) polymorphisms (6p21.3) are associated with periodontal parameters. Parogene 1 population samples (n = 169) were analyzed with 13,245 single nucleotide polymorphisms (SNPs) of the MHC region. Eighteen selected SNPs (P ≤ 0.001) were replicated in Parogene 2 population samples (n = 339) and the Health 2000 Survey (n = 1,420). All subjects had a detailed clinical and radiographic oral health examination. Serum lymphotoxin-α (LTA) concentrations were measured in the Parogene populations, and the protein was detected in inflamed periodontal tissue. In the Parogene 1 population, 10 SNPs were associated with periodontal parameters. The strongest associations emerged from the parameters bleeding on probing (BOP) and a probing pocket depth (PPD) of ≥6 mm with the genes BAT1, NFKBIL1, and LTA. Six SNPs, rs11796, rs3130059, rs2239527, rs2071591, rs909253, and rs1041981 (r(2), ≥0.92), constituted a risk haplotype. In the Parogene 1 population, the haplotype had the strongest association with the parameter BOP, a PPD of ≥6 mm, and severe periodontitis with odds ratios (95% confidence intervals) of 2.63 (2.21 to 3.20), 2.90 (2.37 to 3.52), and 3.10 (1.63 to 5.98), respectively. These results were replicated in the other two populations. High serum LTA concentrations in the Parogene population were associated with the periodontitis risk alleles of the LTA SNPs (rs909253 and rs1041981) of the haplotype. In addition, the protein was expressed in inflamed gingival connective tissue. We identified a novel BAT1-NFKBIL1-LTA haplotype as a significant contributor to the risk of periodontitis. The genetic polymorphisms in the MHC class III region may be functionally important in periodontitis susceptibility.

Genome-Wide Association Study of Peripheral Artery Disease.

BACKGROUND: Peripheral artery disease (PAD) affects >200 million people worldwide and is associated with high mortality and morbidity. We sought to identify genomic variants associated with PAD overall and in the contexts of diabetes and smoking status. METHODS: We identified genetic variants associated with PAD and then meta-analyzed with published summary statistics from the Million Veterans Program and UK Biobank to replicate their findings. Next, we ran stratified genome-wide association analysis in ever smokers, never smokers, individuals with diabetes, and individuals with no history of diabetes and corresponding interaction analyses, to identify variants that modify the risk of PAD by diabetic or smoking status. RESULTS: We identified 5 genome-wide significant (Passociation ≤5×10-8) associations with PAD in 449 548 (Ncases=12 086) individuals of European ancestry near LPA (lipoprotein [a]), CDKN2BAS1 (CDKN2B antisense RNA 1), SH2B3 (SH2B adaptor protein 3) - PTPN11 (protein tyrosine phosphatase non-receptor type 11), HDAC9 (histone deacetylase 9), and CHRNA3 (cholinergic receptor nicotinic alpha 3 subunit) loci (which overlapped previously reported associations). Meta-analysis with variants previously associated with PAD showed that 18 of 19 published variants remained genome-wide significant. In individuals with diabetes, rs116405693 at the CCSER1 (coiled-coil serine rich protein 1) locus was associated with PAD (odds ratio [95% CI], 1.51 [1.32-1.74], Pdiabetes=2.5×10-9, Pinteractionwithdiabetes=5.3×10-7). Furthermore, in smokers, rs12910984 at the CHRNA3 locus was associated with PAD (odds ratio [95% CI], 1.15 [1.11-1.19], Psmokers=9.3×10-10, Pinteractionwithsmoking=3.9×10-5). CONCLUSIONS: Our analyses confirm the published genetic associations with PAD and identify novel variants that may influence susceptibility to PAD in the context of diabetes or smoking status.

Genetic Predisposition to Coronary Artery Disease in Type 2 Diabetes Mellitus.

BACKGROUND: Coronary artery disease (CAD) is accelerated in subjects with type 2 diabetes mellitus (T2D). METHODS: To test whether this reflects differential genetic influences on CAD risk in subjects with T2D, we performed a systematic assessment of genetic overlap between CAD and T2D in 66 643 subjects (27 708 with CAD and 24 259 with T2D). Variants showing apparent association with CAD in stratified analyses or evidence of interaction were evaluated in a further 117 787 subjects (16 694 with CAD and 11 537 with T2D). RESULTS: None of the previously characterized CAD loci was found to have specific effects on CAD in T2D individuals, and a genome-wide interaction analysis found no new variants for CAD that could be considered T2D specific. When we considered the overall genetic correlations between CAD and its risk factors, we found no substantial differences in these relationships by T2D background. CONCLUSIONS: This study found no evidence that the genetic architecture of CAD differs in those with T2D compared with those without T2D.

Genetic Variants Contributing to Circulating Matrix Metalloproteinase 8 Levels and Their Association With Cardiovascular Diseases: A Genome-Wide Analysis.

BACKGROUND: Matrix metalloproteinase 8 (MMP-8) is a proinflammatory enzyme expressed mainly by neutrophils. Elevated serum and plasma concentrations of MMP-8 are associated with the risk for and outcome of cardiovascular diseases (CVDs). The origin of circulating MMP-8 is not completely clear. METHODS AND RESULTS: We performed a genome-wide association study of serum MMP-8 levels in 2 populations comprising altogether 6049 individuals. Moreover, we studied whether MMP-8-associated variants are linked to increased risk of CVDs and overall mortality in >20 000 subjects. The strongest association with serum MMP-8 was found in locus 1q31.3, containing the gene for complement factor H (lead single nucleotide polymorphism: rs800292; P=2.4×10-35). In functional experiments, activation of the alternative pathway of complement in the carriers of rs800292 minor allele (Ile62 in factor H) led to decreased release of MMP-8 from neutrophils compared with the major allele (Val62 in factor H). Another association was detected in 1q21.3, containing genes S100A8, S100A9, and S100A12 (strongest association: rs1560833; P=5.3×10-15). The minor allele of rs1560833 was inversely associated with CVD (odds ratio [95% confidence interval]: 0.90 [0.82-0.99]; P=0.032) and the time to incident CVD event (hazard ratio [95% confidence interval]: 0.91 [0.84-0.99]; P=0.032) in men but not in women. CONCLUSIONS: According to our results, the activation of the alternative pathway of the complement system strongly contributes to serum MMP-8 concentration. Genetic polymorphism in S100A9-S100A12-S100A8 locus affects serum and plasma MMP-8 and shows a suggestive association with the risk of CVDs. Our results show that genetic variation determines a significant portion of circulating MMP-8 concentrations.

Relations between lipoprotein(a) concentrations, LPA genetic variants, and the risk of mortality in patients with established coronary heart disease: a molecular and genetic association study.

BACKGROUND: Lipoprotein(a) concentrations in plasma are associated with cardiovascular risk in the general population. Whether lipoprotein(a) concentrations or LPA genetic variants predict long-term mortality in patients with established coronary heart disease remains less clear. METHODS: We obtained data from 3313 patients with established coronary heart disease in the Ludwigshafen Risk and Cardiovascular Health (LURIC) study. We tested associations of tertiles of lipoprotein(a) concentration in plasma and two LPA single-nucleotide polymorphisms ([SNPs] rs10455872 and rs3798220) with all-cause mortality and cardiovascular mortality by Cox regression analysis and with severity of disease by generalised linear modelling, with and without adjustment for age, sex, diabetes diagnosis, systolic blood pressure, BMI, smoking status, estimated glomerular filtration rate, LDL-cholesterol concentration, and use of lipid-lowering therapy. Results for plasma lipoprotein(a) concentrations were validated in five independent studies involving 10 195 patients with established coronary heart disease. Results for genetic associations were replicated through large-scale collaborative analysis in the GENIUS-CHD consortium, comprising 106 353 patients with established coronary heart disease and 19 332 deaths in 22 studies or cohorts. FINDINGS: The median follow-up was 9·9 years. Increased severity of coronary heart disease was associated with lipoprotein(a) concentrations in plasma in the highest tertile (adjusted hazard radio [HR] 1·44, 95% CI 1·14-1·83) and the presence of either LPA SNP (1·88, 1·40-2·53). No associations were found in LURIC with all-cause mortality (highest tertile of lipoprotein(a) concentration in plasma 0·95, 0·81-1·11 and either LPA SNP 1·10, 0·92-1·31) or cardiovascular mortality (0·99, 0·81-1·2 and 1·13, 0·90-1·40, respectively) or in the validation studies. INTERPRETATION: In patients with prevalent coronary heart disease, lipoprotein(a) concentrations and genetic variants showed no associations with mortality. We conclude that these variables are not useful risk factors to measure to predict progression to death after coronary heart disease is established. FUNDING: Seventh Framework Programme for Research and Technical Development (AtheroRemo and RiskyCAD), INTERREG IV Oberrhein Programme, Deutsche Nierenstiftung, Else-Kroener Fresenius Foundation, Deutsche Stiftung für Herzforschung, Deutsche Forschungsgemeinschaft, Saarland University, German Federal Ministry of Education and Research, Willy Robert Pitzer Foundation, and Waldburg-Zeil Clinics Isny.

Associations of functional alanine-glyoxylate aminotransferase 2 gene variants with atrial fibrillation and ischemic stroke.

Asymmetric and symmetric dimethylarginines (ADMA and SDMA) impair nitric oxide bioavailability and have been implicated in the pathogenesis of atrial fibrillation (AF). Alanine-glyoxylate aminotransferase 2 (AGXT2) is the only enzyme capable of metabolizing both of the dimethylarginines. We hypothesized that two functional AGXT2 missense variants (rs37369, V140I; rs16899974, V498L) are associated with AF and its cardioembolic complications. Association analyses were conducted using 1,834 individulas with AF and 7,159 unaffected individuals from two coronary angiography cohorts and a cohort comprising patients undergoing clinical exercise testing. In coronary angiography patients without structural heart disease, the minor A allele of rs16899974 was associated with any AF (OR = 2.07, 95% CI 1.59-2.68), and with paroxysmal AF (OR = 1.98, 95% CI 1.44-2.74) and chronic AF (OR = 2.03, 95% CI 1.35-3.06) separately. We could not replicate the association with AF in the other two cohorts. However, the A allele of rs16899974 was nominally associated with ischemic stroke risk in the meta-analysis of WTCCC2 ischemic stroke cohorts (3,548 cases, 5,972 controls) and with earlier onset of first-ever ischemic stroke (360 cases) in the cohort of clinical exercise test patients. In conclusion, AGXT2 variations may be involved in the pathogenesis of AF and its age-related thromboembolic complications.

Novel 6p21.3 Risk Haplotype Predisposes to Acute Coronary Syndrome.

BACKGROUND: The HLA-DRB1*01 allele of the human leukocyte antigen has been associated with acute coronary syndrome. Genome-wide association studies have revealed associations with human leukocyte antigen and non-human leukocyte antigen genes of 3 major histocompatibility complex gene classes but not at allelic level. METHODS AND RESULTS: We conducted a large-scale genetic analysis on a case-control cohort comprising 5376 acute coronary syndrome cases and 4852 unrelated controls from 4 populations of 2 European countries. We analyzed the risk candidate allele of HLA-DRB1*01 by genomic real-time polymerase chain reaction together with high-density single nucleotide polymorphisms of the major histocompatibility complex to precisely identify risk loci for acute coronary syndrome with effective clinical implications. We found a risk haplotype for the disease containing single nucleotide polymorphisms from BTNL2 and HLA-DRA genes and the HLA-DRB1*01 allele. The association of the haplotype appeared in 3 of the 4 populations, and the direction of the effect was consistent in the fourth. Coronary samples from subjects homozygous for the disease-associated haplotype showed higher BTNL2 mRNA levels (r=0.760; P<0.00001).We localized, with immunofluorescence staining, BTNL2 in CD68-positive macrophages of the coronary artery plaques. In homozygous cases, BTNL2 blocking, in T-cell stimulation assays, enhanced CD4(+)FOXP3(+) regulatory T cell proliferation significantly (blocking versus nonblocking; P<0.05). CONCLUSIONS: In cases with the risk haplotype for acute coronary syndrome, these results suggest involvement of enhanced immune reactions. BTNL2 may have an inhibitory effect on FOXP3(+) T cell proliferation, especially in patients homozygous for the risk alleles. CLINICAL TRIAL REGISTRATION: https://www.clinicaltrials.gov; Unique Identifier: NCT00417534.

Diversity of extended HLA-DRB1 haplotypes in the Finnish population.

The Major Histocompatibility Complex (MHC, 6p21) codes for traditional HLA and other host response related genes. The polymorphic HLA-DRB1 gene in MHC Class II has been associated with several complex diseases. In this study we focus on MHC haplotype structures in the Finnish population. We explore the variability of extended HLA-DRB1 haplotypes in relation to the other traditional HLA genes and a selected group of MHC class III genes. A total of 150 healthy Finnish individuals were included in the study. Subjects were genotyped for HLA alleles (HLA-A, -B, -DRB1, -DQB1, and -DPB1). The polymorphism of TNF, LTA, C4, BTNL2 and HLA-DRA genes was studied with 74 SNPs (single nucleotide polymorphism). The C4A and C4B gene copy numbers and a 2-bp silencing insertion at exon 29 in C4A gene were analysed with quantitative genomic realtime-PCR. The allele frequencies for each locus were calculated and haplotypes were constructed using both the traditional HLA alleles and SNP blocks. The most frequent Finnish A∼B∼DR -haplotype, uncommon in elsewhere in Europe, was A*03∼B*35∼DRB1*01∶01. The second most common haplotype was a common European ancestral haplotype AH 8.1 (A*01∼B*08∼DRB1*03∶01). Extended haplotypes containing HLA-B, TNF block, C4 and HLA-DPB1 strongly increased the number of HLA-DRB1 haplotypes showing variability in the extended HLA-DRB1 haplotype structures. On the contrary, BTNL2 block and HLA-DQB1 were more conserved showing linkage with the HLA-DRB1 alleles. We show that the use of HLA-DRB1 haplotypes rather than single HLA-DRB1 alleles is advantageous when studying the polymorphisms and LD patters of the MHC region. For disease association studies the HLA-DRB1 haplotypes with various MHC markers allows us to cluster haplotypes with functionally important gene variants such as C4 deficiency and cytokines TNF and LTA, and provides hypotheses for further assessment. Our study corroborates the importance of studying population-specific MHC haplotypes.