Veronique Vitart

Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts

Recent genome-wide association (GWA) studies of lipids have been conducted in samples ascertained for other phenotypes, particularly diabetes. Here we report the first GWA analysis of loci affecting total cholesterol (TC), low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol and triglycerides sampled randomly from 16 population-based cohorts and genotyped using mainly the Illumina HumanHap300-Duo platform. Our study included a total of 17,797–22,562 persons, aged 18–104 years and from geographic regions spanning from the Nordic countries to Southern Europe. We established 22 loci associated with serum lipid levels at a genome-wide significance level (P < 5 × 10−8), including 16 loci that were identified by previous GWA studies. The six newly identified loci in our cohort samples are ABCG5 (TC, P = 1.5 × 10−11; LDL, P = 2.6 × 10−10), TMEM57 (TC, P = 5.4 × 10−10), CTCF-PRMT8 region (HDL, P = 8.3 × 10−16), DNAH11 (LDL, P = 6.1 × 10−9), FADS3-FADS2 (TC, P = 1.5 × 10−10; LDL, P = 4.4 × 10−13) and MADD-FOLH1 region (HDL, P = 6 × 10−11). For three loci, effect sizes differed significantly by sex. Genetic risk scores based on lipid loci explain up to 4.8% of variation in lipids and were also associated with increased intima media thickness (P = 0.001) and coronary heart disease incidence (P = 0.04). The genetic risk score improves the screening of high-risk groups of dyslipidemia over classical risk factors.

SLC2A9 is a newly identified urate transporter influencing serum urate concentration, urate excretion and gout

Uric acid is the end product of purine metabolism in humans and great apes, which have lost hepatic uricase activity, leading to uniquely high serum uric acid concentrations (200–500 μM) compared with other mammals (3–120 μM). About 70% of daily urate disposal occurs via the kidneys, and in 5–25% of the human population, impaired renal excretion leads to hyperuricemia. About 10% of people with hyperuricemia develop gout, an inflammatory arthritis that results from deposition of monosodium urate crystals in the joint. We have identified genetic variants within a transporter gene, SLC2A9, that explain 1.7–5.3% of the variance in serum uric acid concentrations, following a genome-wide association scan in a Croatian population sample. SLC2A9 variants were also associated with low fractional excretion of uric acid and/or gout in UK, Croatian and German population samples. SLC2A9 is a known fructose transporter, and we now show that it has strong uric acid transport activity in Xenopus laevis oocytes.

Meta-Analysis of 28,141 Individuals Identifies Common Variants within Five New Loci That Influence Uric Acid Concentrations

Elevated serum uric acid levels cause gout and are a risk factor for cardiovascular disease and diabetes. To investigate the polygenetic basis of serum uric acid levels, we conducted a meta-analysis of genome-wide association scans from 14 studies totalling 28,141 participants of European descent, resulting in identification of 954 SNPs distributed across nine loci that exceeded the threshold of genome-wide significance, five of which are novel. Overall, the common variants associated with serum uric acid levels fall in the following nine regions: SLC2A9 (p = 5.2×10−201), ABCG2 (p = 3.1×10−26), SLC17A1 (p = 3.0×10−14), SLC22A11 (p = 6.7×10−14), SLC22A12 (p = 2.0×10−9), SLC16A9 (p = 1.1×10−8), GCKR (p = 1.4×10−9), LRRC16A (p = 8.5×10−9), and near PDZK1 (p = 2.7×10−9). Identified variants were analyzed for gender differences. We found that the minor allele for rs734553 in SLC2A9 has greater influence in lowering uric acid levels in women and the minor allele of rs2231142 in ABCG2 elevates uric acid levels more strongly in men compared to women. To further characterize the identified variants, we analyzed their association with a panel of metabolites. rs12356193 within SLC16A9 was associated with DL-carnitine (p = 4.0×10−26) and propionyl-L-carnitine (p = 5.0×10−8) concentrations, which in turn were associated with serum UA levels (p = 1.4×10−57 and p = 8.1×10−54, respectively), forming a triangle between SNP, metabolites, and UA levels. Taken together, these associations highlight additional pathways that are important in the regulation of serum uric acid levels and point toward novel potential targets for pharmacological intervention to prevent or treat hyperuricemia. In addition, these findings strongly support the hypothesis that transport proteins are key in regulating serum uric acid levels.

A Comprehensive Evaluation of Potential Lung Function Associated Genes in the SpiroMeta General Population Sample

Rationale Lung function measures are heritable traits that predict population morbidity and mortality and are essential for the diagnosis of chronic obstructive pulmonary disease (COPD). Variations in many genes have been reported to affect these traits, but attempts at replication have provided conflicting results. Recently, we undertook a meta-analysis of Genome Wide Association Study (GWAS) results for lung function measures in 20,288 individuals from the general population (the SpiroMeta consortium). Objectives To comprehensively analyse previously reported genetic associations with lung function measures, and to investigate whether single nucleotide polymorphisms (SNPs) in these genomic regions are associated with lung function in a large population sample. Methods We analysed association for SNPs tagging 130 genes and 48 intergenic regions (+/−10 kb), after conducting a systematic review of the literature in the PubMed database for genetic association studies reporting lung function associations. Results The analysis included 16,936 genotyped and imputed SNPs. No loci showed overall significant association for FEV1 or FEV1/FVC traits using a carefully defined significance threshold of 1.3×10−5. The most significant loci associated with FEV1 include SNPs tagging MACROD2 (P = 6.81×10−5), CNTN5 (P = 4.37×10−4), and TRPV4 (P = 1.58×10−3). Among ever-smokers, SERPINA1 showed the most significant association with FEV1 (P = 8.41×10−5), followed by PDE4D (P = 1.22×10−4). The strongest association with FEV1/FVC ratio was observed with ABCC1 (P = 4.38×10−4), and ESR1 (P = 5.42×10−4) among ever-smokers. Conclusions Polymorphisms spanning previously associated lung function genes did not show strong evidence for association with lung function measures in the SpiroMeta consortium population. Common SERPINA1 polymorphisms may affect FEV1 among smokers in the general population.

NRXN3 Is a Novel Locus for Waist Circumference: A Genome-Wide Association Study from the CHARGE Consortium

Central abdominal fat is a strong risk factor for diabetes and cardiovascular disease. To identify common variants influencing central abdominal fat, we conducted a two-stage genome-wide association analysis for waist circumference (WC). In total, three loci reached genome-wide significance. In stage 1, 31,373 individuals of Caucasian descent from eight cohort studies confirmed the role of FTO and MC4R and identified one novel locus associated with WC in the neurexin 3 gene [NRXN3 (rs10146997, p = 6.4×10−7)]. The association with NRXN3 was confirmed in stage 2 by combining stage 1 results with those from 38,641 participants in the GIANT consortium (p = 0.009 in GIANT only, p = 5.3×10−8 for combined analysis, n = 70,014). Mean WC increase per copy of the G allele was 0.0498 z-score units (0.65 cm). This SNP was also associated with body mass index (BMI) [p = 7.4×10−6, 0.024 z-score units (0.10 kg/m2) per copy of the G allele] and the risk of obesity (odds ratio 1.13, 95% CI 1.07–1.19; p = 3.2×10−5 per copy of the G allele). The NRXN3 gene has been previously implicated in addiction and reward behavior, lending further evidence that common forms of obesity may be a central nervous system-mediated disorder. Our findings establish that common variants in NRXN3 are associated with WC, BMI, and obesity.

Novel Associations of Multiple Genetic Loci With Plasma Levels of Factor VII, Factor VIII, and von Willebrand Factor The CHARGE (Cohorts for Heart and Aging Research in Genome Epidemiology) Consortium

A complex cascade of coagulation factors underlies hemostasis and prevents life-threatening blood loss from damaged blood vessels. The hemostatic factors VII and VIII, both produced in the liver, play central roles in the initiation and propagation, respectively, of fibrin formation. In the tissue-factor pathway, blood coagulation factor VII (FVII), once activated, serves as a catalyst for factor X (FX) activation, which converts prothrombin to thrombin. During propagation, activated factor VIII (FVIII) activates FX in the presence of activated factor IX. Von Willebrand factor (vWF), produced by endothelial cells and megakaryocytes, has multiple roles in hemostasis. Its primary role is to serve as an adhesion molecule that anchors platelets to exposed collagen after endothelial cell damage. The factor also acts as a carrier protein of FVIII, thereby prolonging the half-life of FVIII. Elevated circulating levels of FVIII and vWF are risk factors for venous thrombosis but the data supporting an association of FVII levels with arterial thrombosis are less consistent.1-5 Hemorrhagic complications are associated with deficiency in FVII and vWF (von Willebrand disease), as well as X-linked deficiency in FVIII (Hemophilia A).6-9 Plasma levels of these proteins are affected by environmental factors but they also are genetically influenced.10-13 Heritability estimates range from 0.53-0.63 for FVII, 0.40-0.61 for FVIII, and 0.31-0.75 for vWF.12, 13 To date, our understanding of genetic variation influencing plasma levels has been focused primarily on cis-acting variation in the genes encoding each protein product (F7, F8, and VWF, respectively). A large-scale genome-wide investigation of the genomic correlates of plasma levels has not been previously published. Using data from 23,608 adults, we investigated genome-wide associations between common genetic variation and plasma levels of FVII, FVIII, and vWF.Background— Plasma levels of coagulation factors VII (FVII), VIII (FVIII), and von Willebrand factor (vWF) influence risk of hemorrhage and thrombosis. We conducted genome-wide association studies to identify new loci associated with plasma levels. Methods and Results— The setting of the study included 5 community-based studies for discovery comprising 23 608 European-ancestry participants: Atherosclerosis Risk In Communities Study, Cardiovascular Health Study, British 1958 Birth Cohort, Framingham Heart Study, and Rotterdam Study. All subjects had genome-wide single-nucleotide polymorphism (SNP) scans and at least 1 phenotype measured: FVII activity/antigen, FVIII activity, and vWF antigen. Each study used its genotype data to impute to HapMap SNPs and independently conducted association analyses of hemostasis measures using an additive genetic model. Study findings were combined by meta-analysis. Replication was conducted in 7604 participants not in the discovery cohort. For FVII, 305 SNPs exceeded the genome-wide significance threshold of 5.0×10−8 and comprised 5 loci on 5 chromosomes: 2p23 (smallest P value 6.2×10−24), 4q25 (3.6×10−12), 11q12 (2.0×10−10), 13q34 (9.0×10−259), and 20q11.2 (5.7×10−37). Loci were within or near genes, including 4 new candidate genes and F7 (13q34). For vWF, 400 SNPs exceeded the threshold and marked 8 loci on 6 chromosomes: 6q24 (1.2×10−22), 8p21 (1.3×10−16), 9q34 (<5.0×10−324), 12p13 (1.7×10−32), 12q23 (7.3×10−10), 12q24.3 (3.8×10−11), 14q32 (2.3×10−10), and 19p13.2 (1.3×10−9). All loci were within genes, including 6 new candidate genes, as well as ABO (9q34) and VWF (12p13). For FVIII, 5 loci were identified and overlapped vWF findings. Nine of the 10 new findings were replicated. Conclusions— New genetic associations were discovered outside previously known biological pathways and may point to novel prevention and treatment targets of hemostasis disorders.

A General Approach for Haplotype Phasing across the Full Spectrum of Relatedness

Many existing cohorts contain a range of relatedness between genotyped individuals, either by design or by chance. Haplotype estimation in such cohorts is a central step in many downstream analyses. Using genotypes from six cohorts from isolated populations and two cohorts from non-isolated populations, we have investigated the performance of different phasing methods designed for nominally ‘unrelated’ individuals. We find that SHAPEIT2 produces much lower switch error rates in all cohorts compared to other methods, including those designed specifically for isolated populations. In particular, when large amounts of IBD sharing is present, SHAPEIT2 infers close to perfect haplotypes. Based on these results we have developed a general strategy for phasing cohorts with any level of implicit or explicit relatedness between individuals. First SHAPEIT2 is run ignoring all explicit family information. We then apply a novel HMM method (duoHMM) to combine the SHAPEIT2 haplotypes with any family information to infer the inheritance pattern of each meiosis at all sites across each chromosome. This allows the correction of switch errors, detection of recombination events and genotyping errors. We show that the method detects numbers of recombination events that align very well with expectations based on genetic maps, and that it infers far fewer spurious recombination events than Merlin. The method can also detect genotyping errors and infer recombination events in otherwise uninformative families, such as trios and duos. The detected recombination events can be used in association scans for recombination phenotypes. The method provides a simple and unified approach to haplotype estimation, that will be of interest to researchers in the fields of human, animal and plant genetics.

Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus

Central corneal thickness (CCT) is associated with eye conditions including keratoconus and glaucoma. We performed a meta-analysis on >20,000 individuals in European and Asian populations that identified 16 new loci associated with CCT at genome-wide significance (P < 5 × 10−8). We further showed that 2 CCT-associated loci, FOXO1 and FNDC3B, conferred relatively large risks for keratoconus in 2 cohorts with 874 cases and 6,085 controls (rs2721051 near FOXO1 had odds ratio (OR) = 1.62, 95% confidence interval (CI) = 1.4–1.88, P = 2.7 × 10−10, and rs4894535 in FNDC3B had OR = 1.47, 95% CI = 1.29–1.68, P = 4.9 × 10−9). FNDC3B was also associated with primary open-angle glaucoma (P = 5.6 × 10−4; tested in 3 cohorts with 2,979 cases and 7,399 controls). Further analyses implicate the collagen and extracellular matrix pathways in the regulation of CCT.

Genetic determinants of circulating sphingolipid concentrations in European populations

Sphingolipids have essential roles as structural components of cell membranes and in cell signalling, and disruption of their metabolism causes several diseases, with diverse neurological, psychiatric, and metabolic consequences. Increasingly, variants within a few of the genes that encode enzymes involved in sphingolipid metabolism are being associated with complex disease phenotypes. Direct experimental evidence supports a role of specific sphingolipid species in several common complex chronic disease processes including atherosclerotic plaque formation, myocardial infarction (MI), cardiomyopathy, pancreatic β-cell failure, insulin resistance, and type 2 diabetes mellitus. Therefore, sphingolipids represent novel and important intermediate phenotypes for genetic analysis, yet little is known about the major genetic variants that influence their circulating levels in the general population. We performed a genome-wide association study (GWAS) between 318,237 single-nucleotide polymorphisms (SNPs) and levels of circulating sphingomyelin (SM), dihydrosphingomyelin (Dih-SM), ceramide (Cer), and glucosylceramide (GluCer) single lipid species (33 traits); and 43 matched metabolite ratios measured in 4,400 subjects from five diverse European populations. Associated variants (32) in five genomic regions were identified with genome-wide significant corrected p-values ranging down to 9.08×10−66. The strongest associations were observed in or near 7 genes functionally involved in ceramide biosynthesis and trafficking: SPTLC3, LASS4, SGPP1, ATP10D, and FADS1–3. Variants in 3 loci (ATP10D, FADS3, and SPTLC3) associate with MI in a series of three German MI studies. An additional 70 variants across 23 candidate genes involved in sphingolipid-metabolizing pathways also demonstrate association (p = 10−4 or less). Circulating concentrations of several key components in sphingolipid metabolism are thus under strong genetic control, and variants in these loci can be tested for a role in the development of common cardiovascular, metabolic, neurological, and psychiatric diseases.

Genomics Meets Glycomics—The First GWAS Study of Human N-Glycome Identifies HNF1α as a Master Regulator of Plasma Protein Fucosylation

Over half of all proteins are glycosylated, and alterations in glycosylation have been observed in numerous physiological and pathological processes. Attached glycans significantly affect protein function; but, contrary to polypeptides, they are not directly encoded by genes, and the complex processes that regulate their assembly are poorly understood. A novel approach combining genome-wide association and high-throughput glycomics analysis of 2,705 individuals in three population cohorts showed that common variants in the Hepatocyte Nuclear Factor 1α (HNF1α) and fucosyltransferase genes FUT6 and FUT8 influence N-glycan levels in human plasma. We show that HNF1α and its downstream target HNF4α regulate the expression of key fucosyltransferase and fucose biosynthesis genes. Moreover, we show that HNF1α is both necessary and sufficient to drive the expression of these genes in hepatic cells. These results reveal a new role for HNF1α as a master transcriptional regulator of multiple stages in the fucosylation process. This mechanism has implications for the regulation of immunity, embryonic development, and protein folding, as well as for our understanding of the molecular mechanisms underlying cancer, coronary heart disease, and metabolic and inflammatory disorders.