Genome-wide association study identifies novel loci associated with circulating phospho- and sphingolipid concentrations

Phospho- and sphingolipids are crucial cellular and intracellular compounds. These lipids are required for active transport, a number of enzymatic processes, membrane formation, and cell signalling. Disruption of their metabolism leads to several diseases, with diverse neurological, psychiatric, and metabolic consequences. A large number of phospholipid and sphingolipid species can be detected and measured in human plasma. We conducted a meta-analysis of five European family-based genome-wide association studies (N = 4034) on plasma levels of 24 sphingomyelins (SPM), 9 ceramides (CER), 57 phosphatidylcholines (PC), 20 lysophosphatidylcholines (LPC), 27 phosphatidylethanolamines (PE), and 16 PE-based plasmalogens (PLPE), as well as their proportions in each major class. This effort yielded 25 genome-wide significant loci for phospholipids (smallest P-value = 9.88×10⁻²⁰⁴) and 10 loci for sphingolipids (smallest P-value = 3.10×10⁻⁵⁷). After a correction for multiple comparisons (P-value<2.2×10⁻⁹), we observed four novel loci significantly associated with phospholipids (PAQR9, AGPAT1, PKD2L1, PDXDC1) and two with sphingolipids (PLD2 and APOE) explaining up to 3.1% of the variance. Further analysis of the top findings with respect to within class molar proportions uncovered three additional loci for phospholipids (PNLIPRP2, PCDH20, and ABDH3) suggesting their involvement in either fatty acid elongation/saturation processes or fatty acid specific turnover mechanisms. Among those, 14 loci (KCNH7, AGPAT1, PNLIPRP2, SYT9, FADS1-2-3, DLG2, APOA1, ELOVL2, CDK17, LIPC, PDXDC1, PLD2, LASS4, and APOE) mapped into the glycerophospholipid and 12 loci (ILKAP, ITGA9, AGPAT1, FADS1-2-3, APOA1, PCDH20, LIPC, PDXDC1, SGPP1, APOE, LASS4, and PLD2) to the sphingolipid pathways. In large meta-analyses, associations between FADS1-2-3 and carotid intima media thickness, AGPAT1 and type 2 diabetes, and APOA1 and coronary artery disease were observed. In conclusion, our study identified nine novel phospho- and sphingolipid loci, substantially increasing our knowledge of the genetic basis for these traits.

Phospho- and sphingolipids are integral to membrane formation and are involved in crucial cellular functions such as signalling, membrane fluidity, membrane protein trafficking, neurotransmission, and receptor trafficking. In addition to severe monogenic diseases resulting from defective phospho- and sphingolipid function and metabolism, the evidence suggests that variations in these lipid levels at the population level are involved in the determination of cardiovascular and neurologic traits and subsequent disease. We took advantage of modern laboratory methods, including microarray-based genotyping and electrospray ionization tandem mass spectrometry, to hunt for genetic variation influencing the levels of more than 350 phospho- and sphingolipid phenotypes. We identified nine novel loci, in addition to confirming a number of previously described loci. Several other genetic regions provided substantial evidence of their involvement in these traits. All of these loci are strong candidates for further research in the field of lipid biology and are likely to yield considerable insights into the complex metabolic pathways underlying circulating phospho- and sphingolipid levels. Understanding these mechanisms might help to illuminate factors leading to the development of common cardiovascular and neurological diseases and might provide molecular targets for the development of new therapies.

Linkage and genome-wide association analysis of obesity-related phenotypes: association of weight with the MGAT1 gene

As major risk-factors for diabetes and cardiovascular diseases, the genetic contribution to obesity-related traits has been of interest for decades. Recently, a limited number of common genetic variants, which have replicated in different populations, have been identified. One approach to increase the statistical power in genetic mapping studies is to focus on populations with increased levels of linkage disequilibrium (LD) and reduced genetic diversity. We have performed joint linkage and genome-wide association analyses for weight and BMI in 3,448 (linkage) and 3,925 (association) partly overlapping healthy individuals from five European populations. A total of four chromosomal regions (two for weight and two for BMI) showed suggestive linkage (lod >2.69) either in one of the populations or in the joint data. At the genome-wide level (nominal P < 1.6 x 10(-7), Bonferroni-adjusted P < 0.05) one single-nucleotide polymorphism (SNP) (rs12517906) (nominal P = 7.3 x 10(-8)) was associated with weight, whereas none with BMI. The SNP associated with weight is located close to MGAT1. The monoacylglycerol acyltransferase (MGAT) enzyme family is known to be involved in dietary fat absorption. There was no overlap between the linkage regions and the associated SNPs. Our results show that genetic effects influencing weight and BMI are shared across diverse European populations, even though some of these populations have experienced recent population bottlenecks and/or been affected by genetic drift. The analysis enabled us to identify a new candidate gene, MGAT1, associated with weight in women.

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 beta-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.08x10(-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.

Extended haplotypes in the growth hormone releasing hormone receptor gene (GHRHR) are associated with normal variation in height

Mutations in the gene for growth hormone releasing hormone receptor (GHRHR) cause isolated growth hormone deficiency (IGHD) but this gene has not been found to affect normal variation in height. We performed a whole genome linkage analysis for height in a population from northern Sweden and identified a region on chromosome 7 with a lod-score of 4.7. The GHRHR gene is located in this region and typing of tagSNPs identified a haplotype that is associated with height (p = 0.00077) in the original study population. Analysis of a sample from an independent population from the most northern part of Sweden also showed an association with height (p = 0.0039) but with another haplotype in the GHRHR gene. Both haplotypes span the 3′ part of the GHRHR gene, including the region in which most of the mutations in IGHD have been located. The effect size of these haplotypes are larger than that of any gene previously associated with height, which indicates that GHRHR might be one of the most important genes so far identified affecting normal variation in human height.

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 x 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 x 10(-11); LDL, P = 2.6 x 10(-10)), TMEM57 (TC, P = 5.4 x 10(-10)), CTCF-PRMT8 region (HDL, P = 8.3 x 10(-16)), DNAH11 (LDL, P = 6.1 x 10(-9)), FADS3-FADS2 (TC, P = 1.5 x 10(-10); LDL, P = 4.4 x 10(-13)) and MADD-FOLH1 region (HDL, P = 6 x 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.

Common variants in the JAZF1 gene associated with height identified by linkage and genome-wide association analysis

Genes for height have gained interest for decades, but only recently have candidate genes started to be identified. We have performed linkage analysis and genome-wide association for height in approximately 4000 individuals from five European populations. A total of five chromosomal regions showed suggestive linkage and in one of these regions, two SNPs (rs849140 and rs1635852) were associated with height (nominal P = 7.0 x 10(-8) and P = 9.6 x 10(-7), respectively). In total, five SNPs across the genome showed an association with height that reached the threshold of genome-wide significance (nominal P < 1.6 x 10(-7)). The association with height was replicated for two SNPs (rs1635852 and rs849140) using three independent studies (n = 31 077, n=1268 and n = 5746) with overall meta P-values of 9.4 x 10(-10) and 5.3 x 10(-8). These SNPs are located in the JAZF1 gene, which has recently been associated with type II diabetes, prostate and endometrial cancer. JAZF1 is a transcriptional repressor of NR2C2, which results in low IGF1 serum concentrations, perinatal and early postnatal hypoglycemia and growth retardation when knocked out in mice. Both the linkage and association analyses independently identified the JAZF1 region affecting human height. We have demonstrated, through replication in additional independent populations, the consistency of the effect of the JAZF1 SNPs on height. Since this gene also has a key function in the metabolism of growth, JAZF1 represents one of the strongest candidates influencing human height identified so far.

A susceptibility locus for human systemic lupus erythematosus (hSLE1) on chromosome 2q

To identify chromosomal regions containing susceptibility loci for systemic lupus erythematosus (SLE), we performed genome scans in families with multiple SLE patients from Iceland, a geographical and genetic isolate, and from Sweden. A number of chromosomal regions showed maximum lod scores (Z) indicating possible linkage to SLE in both the Icelandic and Swedish families. In the Icelandic families, five regions showed lod scores greater than 2.0, three of which (4p15-13, Z=3.20; 9p22, Z=2.27; 19q13, Z=2.06) are homologous to the murine regions containing the lmb2, sle2 and sle3 loci, respectively. The fourth region is located on 19p13 (D19S247, Z=2.58) and the fifth on 2q37 (D2S125, Z=2.06). Only two regions showed lod scores above 2.0 in the Swedish families: on chromosome 2q11 (D2S436, Z=2. 13) and 2q37 (D2S125, Z=2.18). The combination of both family sets gave a highly significant lod score at D2S125 of Z=4.24 in favor of linkage for 2q37. This region represents a new locus for SLE. Our results underscore the importance of studying well-defined populations for genetic analysis of complex diseases such as SLE.

Genetic analysis of the contribution of IL10 to systemic lupus erythematosus

OBJECTIVE

To study the contribution of the IL10 gene to the susceptibility to systemic lupus erythematosus (SLE).

METHODS

Analysis by fluorescent-semiautomated genotyping of a dinucleotide repeat located in the promoter region of the IL10 locus (microsatellite G).

RESULTS

No significant difference was found in the frequencies of the microsatellite alleles of 330 Mexican patients with SLE compared to 368 controls from the same population. Two-point linkage analyses were carried out using 13 Mexican, 13 Swedish, and 8 Icelandic families with 2 or more cases with SLE. No linkage was revealed between IL10 and SLE, using a variety of parameter settings.

CONCLUSION

Our results do not support that the IL10 gene contributes to the susceptibility to SLE in the populations we studied.

A comparison of genome-scans performed in multicase families with systemic lupus erythematosus from different population groups

Systemic lupus erythematosus is a disease of unknown etiology. Multiple genetic factors are believed to be involved in its pathogenesis. In addition, and due to genetic heterogeneity, these factors and/or their combinations may be different in different ethnic groups, while some might be shared between populations. We have performed genome scans in multicase families from three different population groups, two from Northern Europe, with a high degree of homogeneity, and the third from a recently admixed population of Mexican Mestizos. Although our family material is relatively small, the results presented here show that using family sets from well defined populations are sufficient to detect susceptibility loci for SLE. Our results also reveal the chromosomal regions most likely to contain susceptibility genes for SLE.

Differential permeabilization of membranes by saponin treatment of isolated rat hepatocytes. Release of secretory proteins

Monolayer cultures of rat hepatocytes were treated with increasing concentrations of saponin (prepared from Gypsophila plants) for 30 min at 6 degrees C. Differential permeabilization of the intracellular membranes could be demonstrated: at 0.040 mg of saponin/ml the plasma membrane was permeabilized, as assessed by the release of 50% of the total cellular amount of lactate dehydrogenase, and at 0.20 mg/ml the endoplasmic reticulum was permeabilized, as measured by the release of 50% of pulse-35S-labelled albumin. The Golgi complex was permeabilized at an intermediate saponin concentration, as indicated by the release of homogeneously 35S-labelled albumin; about half the intracellular albumin is located in this organelle. At 1.0 up to 5.0 mg of saponin/ml 90-95% of the radioactively labelled albumin was released. Even at 5.0 mg/ml less than 10% of the membrane of the endoplasmic reticulum was solubilized, as judged by the degree of release of a membrane-bound enzyme specific for this organelle. These results demonstrate the usefulness of saponin as a tool for investigating the interior of different intracellular compartments.