Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity.
Aims Darapladib, a potent inhibitor of lipoprotein-associated phospholipase A(2) (Lp-PLA(2)), has not reduced risk of cardiovascular disease outcomes in recent randomized trials. We aimed to test whether Lp-PLA(2) enzyme activity is causally relevant to coronary heart disease. Methods In 72,657 patients with coronary heart disease and 110,218 controls in 23 epidemiological studies, we genotyped five functional variants: four rare loss-of-function mutations (c.109+2T > C (rs142974898), Arg82His (rs144983904), Val279Phe (rs76863441), Gln287Ter (rs140020965)) and one common modest-impact variant (Val379Ala (rs1051931)) in PLA2G7, the gene encoding Lp-PLA(2). We supplemented de-novo genotyping with information on a further 45,823 coronary heart disease patients and 88,680 controls in publicly available databases and other previous studies. We conducted a systematic review of randomized trials to compare effects of darapladib treatment on soluble Lp-PLA(2) activity, conventional cardiovascular risk factors, and coronary heart disease risk with corresponding effects of Lp-PLA(2)-lowering alleles. Results Lp-PLA(2) activity was decreased by 64% ( p = 2.4 x 10(-25)) with carriage of any of the four loss-of-function variants, by 45% ( p < 10(-300)) for every allele inherited at Val279Phe, and by 2.7% ( p = 1.9 x 10(-12)) for every allele inherited at Val379Ala. Darapladib 160 mg once-daily reduced Lp-PLA(2) activity by 65% ( p < 10(-300)). Causal risk ratios for coronary heart disease per 65% lower Lp-PLA(2) activity were: 0.95 (0.88-1.03) with Val279Phe; 0.92 (0.74-1.16) with carriage of any loss-of-function variant; 1.01 (0.68-1.51) with Val379Ala; and 0.95 (0.89-1.02) with darapladib treatment. Conclusions In a large-scale human genetic study, none of a series of Lp-PLA(2)-lowering alleles was related to coronary heart disease risk, suggesting that Lp-PLA(2) is unlikely to be a causal risk factor.
Importance: The activity of lipoprotein lipase (LPL) is the rate-determining step in clearing triglyceride-rich lipoproteins from the circulation. Mutations that damage the LPL gene (LPL) lead to lifelong deficiency in enzymatic activity and can provide insight into the relationship of LPL to human disease. Objective: To determine whether rare and/or common variants in LPL are associated with early-onset coronary artery disease (CAD). Design, Setting, and Participants: In a cross-sectional study, LPL was sequenced in 10 CAD case-control cohorts of the multinational Myocardial Infarction Genetics Consortium and a nested CAD case-control cohort of the Geisinger Health System DiscovEHR cohort between 2010 and 2015. Common variants were genotyped in up to 305699 individuals of the Global Lipids Genetics Consortium and up to 120600 individuals of the CARDIoGRAM Exome Consortium between 2012 and 2014. Study-specific estimates were pooled via meta-analysis. Exposures: Rare damaging mutations in LPL included loss-of-function variants and missense variants annotated as pathogenic in a human genetics database or predicted to be damaging by computer prediction algorithms trained to identify mutations that impair protein function. Common variants in the LPL gene region included those independently associated with circulating triglyceride levels. Main Outcomes and Measures: Circulating lipid levels and CAD. Results: Among 46891 individuals with LPL gene sequencing data available, the mean (SD) age was 50 (12.6) years and 51% were female. A total of 188 participants (0.40%; 95% CI, 0.35%-0.46%) carried a damaging mutation in LPL, including 105 of 32646 control participants (0.32%) and 83 of 14245 participants with early-onset CAD (0.58%). Compared with 46703 noncarriers, the 188 heterozygous carriers of an LPL damaging mutation displayed higher plasma triglyceride levels (19.6 mg/dL; 95% CI, 4.6-34.6 mg/dL) and higher odds of CAD (odds ratio = 1.84; 95% CI, 1.35-2.51; P < .001). An analysis of 6 common LPL variants resulted in an odds ratio for CAD of 1.51 (95% CI, 1.39-1.64; P = 1.1 x 10-22) per 1-SD increase in triglycerides. Conclusions and Relevance: The presence of rare damaging mutations in LPL was significantly associated with higher triglyceride levels and presence of coronary artery disease. However, further research is needed to assess whether there are causal mechanisms by which heterozygous lipoprotein lipase deficiency could lead to coronary artery disease.
OBJECTIVE: Autosomal recessive hypercholesterolemia is a rare inherited disorder, characterized by extremely high total and low-density lipoprotein cholesterol levels, that has been previously linked to mutations in LDLRAP1. We identified a family with autosomal recessive hypercholesterolemia not explained by mutations in LDLRAP1 or other genes known to cause monogenic hypercholesterolemia. The aim of this study was to identify the molecular pathogenesis of autosomal recessive hypercholesterolemia in this family. APPROACH AND RESULTS: We used exome sequencing to assess all protein-coding regions of the genome in 3 family members and identified a homozygous exon 8 splice junction mutation (c.894G>A, also known as E8SJM) in LIPA that segregated with the diagnosis of hypercholesterolemia. Because homozygosity for mutations in LIPA is known to cause cholesterol ester storage disease, we performed directed follow-up phenotyping by noninvasively measuring hepatic cholesterol content. We observed abnormal hepatic accumulation of cholesterol in the homozygote individuals, supporting the diagnosis of cholesterol ester storage disease. Given previous suggestions of cardiovascular disease risk in heterozygous LIPA mutation carriers, we genotyped E8SJM in >27 000 individuals and found no association with plasma lipid levels or risk of myocardial infarction, confirming a true recessive mode of inheritance. CONCLUSIONS: By integrating observations from Mendelian and population genetics along with directed clinical phenotyping, we diagnosed clinically unapparent cholesterol ester storage disease in the affected individuals from this kindred and addressed an outstanding question about risk of cardiovascular disease in LIPA E8SJM heterozygous carriers.
OBJECTIVE: Diminished serum paraoxonase and arylesterase activities (measures of paraoxonase-1 [PON-1] function) in humans have been linked to heightened systemic oxidative stress and atherosclerosis risk. The clinical prognostic use of measuring distinct PON-1 activities has not been established, and the genetic determinants of PON-1 activities are not known. METHODS AND RESULTS: We established analytically robust high-throughput assays for serum paraoxonase and arylesterase activities and measured these in 3668 stable subjects undergoing elective coronary angiography without acute coronary syndrome and were prospectively followed for major adverse cardiovascular events (MACE= death, myocardial infarction, stroke) over 3 years. Low serum arylesterase and paraoxonase activities were both associated with increased risk for MACE, with arylesterase activity showing greatest prognostic value (quartile 4 versus quartile 1; hazard ratio 2.63; 95% CI, 1.97-3.50; P<0.01). Arylesterase remained significant after adjusting for traditional risk factors, C-reactive protein, and creatinine clearance (hazard ratio, 2.20; 95% CI, 1.60-3.02; P<0.01), predicted future development of MACE in both primary and secondary prevention populations, and reclassified risk categories incrementally to traditional clinical variables. A genome-wide association study identified distinct single nucleotide polymorphisms within the PON-1 gene that were highly significantly associated with serum paraoxonase (1.18x10(-303)) or arylesterase (4.99x10(-116)) activity but these variants were not associated with either 3-year MACE risk in an angiographic cohort (n=2136) or history of either coronary artery disease or myocardial infarction in the Coronary Artery Disease Genome-Wide Replication and Meta-Analysis consortium (n approximately 80 000 subjects). CONCLUSIONS: Diminished serum arylesterase activity, but not the genetic determinants of PON-1 functional measures, provides incremental prognostic value and clinical reclassification of stable subjects at risk of developing MACE.
Genome-wide association studies (GWAS) have successfully identified loci associated with quantitative traits, such as blood lipids. Deep resequencing studies are being utilized to catalogue the allelic spectrum at GWAS loci. The goal of these studies is to identify causative variants and missing heritability, including heritability due to low frequency and rare alleles with large phenotypic impact. Whereas rare variant efforts have primarily focused on nonsynonymous coding variants, we hypothesized that noncoding variants in these loci are also functionally important. Using the HDL-C gene LIPG as an example, we explored the effect of regulatory variants identified through resequencing of subjects at HDL-C extremes on gene expression, protein levels, and phenotype. Resequencing a portion of the LIPG promoter and 5' UTR in human subjects with extreme HDL-C, we identified several rare variants in individuals from both extremes. Luciferase reporter assays were used to measure the effect of these rare variants on LIPG expression. Variants conferring opposing effects on gene expression were enriched in opposite extremes of the phenotypic distribution. Minor alleles of a common regulatory haplotype and noncoding GWAS SNPs were associated with reduced plasma levels of the LIPG gene product endothelial lipase (EL), consistent with its role in HDL-C catabolism. Additionally, we found that a common nonfunctional coding variant associated with HDL-C (rs2000813) is in linkage disequilibrium with a 5' UTR variant (rs34474737) that decreases LIPG promoter activity. We attribute the gene regulatory role of rs34474737 to the observed association of the coding variant with plasma EL levels and HDL-C. Taken together, the findings show that both rare and common noncoding regulatory variants are important contributors to the allelic spectrum in complex trait loci.
BACKGROUND: eQTL analyses are important to improve the understanding of genetic association results. We performed a genome-wide association and global gene expression study to identify functionally relevant variants affecting the risk of coronary artery disease (CAD). METHODS AND RESULTS: In a genome-wide association analysis of 2078 CAD cases and 2953 control subjects, we identified 950 single-nucleotide polymorphisms (SNPs) that were associated with CAD at P<10(-3). Subsequent in silico and wet-laboratory replication stages and a final meta-analysis of 21 428 CAD cases and 38 361 control subjects revealed a novel association signal at chromosome 10q23.31 within the LIPA (lysosomal acid lipase A) gene (P=3.7x10(-8); odds ratio, 1.1; 95% confidence interval, 1.07 to 1.14). The association of this locus with global gene expression was assessed by genome-wide expression analyses in the monocyte transcriptome of 1494 individuals. The results showed a strong association of this locus with expression of the LIPA transcript (P=1.3x10(-96)). An assessment of LIPA SNPs and transcript with cardiovascular phenotypes revealed an association of LIPA transcript levels with impaired endothelial function (P=4.4x10(-3)). CONCLUSIONS: The use of data on genetic variants and the addition of data on global monocytic gene expression led to the identification of the novel functional CAD susceptibility locus LIPA, located on chromosome 10q23.31. The respective eSNPs associated with CAD strongly affect LIPA gene expression level, which was related to endothelial dysfunction, a precursor of CAD.
Elevated plasma concentrations of HDL cholesterol (HDL-C) are associated with protection from atherosclerotic cardiovascular disease. Animal models indicate that decreased expression of endothelial lipase (LIPG) is inversely associated with HDL-C levels, and genome-wide association studies have identified LIPG variants as being associated with HDL-C levels in humans. We hypothesized that loss-of-function mutations in LIPG may result in elevated HDL-C and therefore performed deep resequencing of LIPG exons in cases with elevated HDL-C levels and controls with decreased HDL-C levels. We identified a significant excess of nonsynonymous LIPG variants unique to cases with elevated HDL-C. In vitro lipase activity assays demonstrated that these variants significantly decreased endothelial lipase activity. In addition, a meta-analysis across 5 cohorts demonstrated that the low-frequency Asn396Ser variant is significantly associated with increased HDL-C, while the common Thr111Ile variant is not. Functional analysis confirmed that the Asn396Ser variant has significantly decreased lipase activity both in vitro and in vivo, while the Thr111Ile variant has normal lipase activity. Our results establish that loss-of-function mutations in LIPG lead to increased HDL-C levels and support the idea that inhibition of endothelial lipase may be an effective mechanism to raise HDL-C.
BACKGROUND: Common single-nucleotide polymorphisms (SNPs) that are associated with blood low-density lipoprotein (LDL) or high-density lipoprotein (HDL) cholesterol modestly affect lipid levels. We tested the hypothesis that a combination of such SNPs contributes to the risk of cardiovascular disease. METHODS: We studied SNPs at nine loci in 5414 subjects from the cardiovascular cohort of the Malmo Diet and Cancer Study. We first validated the association between SNPs and either LDL or HDL cholesterol and subsequently created a genotype score on the basis of the number of unfavorable alleles. We used Cox proportional-hazards models to determine the time to the first cardiovascular event in relation to the genotype score. RESULTS: All nine SNPs showed replication of an association with levels of either LDL or HDL cholesterol. With increasing genotype scores, the level of LDL cholesterol increased from 152 mg to 171 mg per deciliter (3.9 to 4.4 mmol per liter), whereas HDL cholesterol decreased from 60 mg to 51 mg per deciliter (1.6 to 1.3 mmol per liter). During follow-up (median, 10.6 years), 238 subjects had a first cardiovascular event. The genotype score was associated with incident cardiovascular disease in models adjusted for covariates including baseline lipid levels (P<0.001). The use of the genotype score did not improve the clinical risk prediction, as assessed by the C statistic. However, there was a significant improvement in risk classification with the use of models that included the genotype score, as compared with those that did not include the genotype score. CONCLUSIONS: A genotype score of nine validated SNPs that are associated with modulation in levels of LDL or HDL cholesterol was an independent risk factor for incident cardiovascular disease. The score did not improve risk discrimination but did modestly improve clinical risk reclassification for individual subjects beyond standard clinical factors.