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⁻⁶⁶. 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⁻⁴ 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.     

Genome-wide association of familial late-onset Alzheimer's disease replicates BIN1 and CLU and nominates CUGBP2 in interaction with APOE

Late-onset Alzheimer''s disease (LOAD) is the most common form of dementia in the elderly. The National Institute of Aging-Late Onset Alzheimer''s Disease Family Study and the National Cell Repository for Alzheimer''s Disease conducted a joint genome-wide association study (GWAS) of multiplex LOAD families (3,839 affected and unaffected individuals from 992 families plus additional unrelated neurologically evaluated normal subjects) using the 610 IlluminaQuad panel. This cohort represents the largest family-based GWAS of LOAD to date, with analyses limited here to the European-American subjects. SNPs near APOE gave highly significant results (e.g., rs2075650, p = 3.2×10⁻⁸¹), but no other genome-wide significant evidence for association was obtained in the full sample. Analyses that stratified on APOE genotypes identified SNPs on chromosome 10p14 in CUGBP2 with genome-wide significant evidence for association within APOE ε4 homozygotes (e.g., rs201119, p = 1.5×10⁻⁸). Association in this gene was replicated in an independent sample consisting of three cohorts. There was evidence of association for recently-reported LOAD risk loci, including BIN1 (rs7561528, p = 0.009 with, and p = 0.03 without, APOE adjustment) and CLU (rs11136000, p = 0.023 with, and p = 0.008 without, APOE adjustment), with weaker support for CR1. However, our results provide strong evidence that association with PICALM (rs3851179, p = 0.69 with, and p = 0.039 without, APOE adjustment) and EXOC3L2 is affected by correlation with APOE, and thus may represent spurious association. Our results indicate that genetic structure coupled with ascertainment bias resulting from the strong APOE association affect genome-wide results and interpretation of some recently reported associations. We show that a locus such as APOE, with large effects and strong association with disease, can lead to samples that require appropriate adjustment for this locus to avoid both false positive and false negative evidence of association. We suggest that similar adjustments may also be needed for many other large multi-site studies.     

Association of Genetic Variants Influencing Lipid Levels with Coronary Artery Disease in Japanese Individuals

Background/Objective

In Japanese populations, we performed a replication study of genetic loci previously identified in European-descent populations as being associated with lipid levels and risk of coronary artery disease (CAD).

Methods

We genotyped 48 single nucleotide polymorphisms (SNPs) from 22 candidate loci that had previously been identified by genome-wide association (GWA) meta-analyses for low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and/or triglycerides in Europeans. We selected 22 loci with 2 parallel tracks from 95 reported loci: 16 significant loci (p<1×10⁻³⁰ in Europeans) and 6 other loci including those with suggestive evidence of lipid associations in 1292 GWA-scanned Japanese samples. Genotyping was done in 4990 general population samples, and 1347 CAD cases and 1337 controls. For 9 SNPs, we further examined CAD associations in an additional panel of 3052 CAD cases and 6335 controls.

Principal Findings

Significant lipid associations (one-tailed p<0.05) were replicated for 18 of 22 loci in Japanese samples, with significant inter-ethnic heterogeneity at 4 loci–APOB, APOE-C1, CETP, and APOA5–and allelic heterogeneity. The strongest association was detected at APOE rs7412 for LDL-C (p = 1.3×10⁻⁴¹), CETP rs3764261 for HDL-C (p = 5.2×10⁻²⁴), and APOA5 rs662799 for triglycerides (p = 5.8×10⁻⁵⁴). CAD association was replicated and/or verified for 4 loci: SORT1 rs611917 (p = 1.7×10⁻⁸), APOA5 rs662799 (p = 0.0014), LDLR rs1433099 (p = 2.1×10⁻⁷), and APOE rs7412 (p = 6.1×10⁻¹³).

Conclusions

Our results confirm that most of the tested lipid loci are associated with lipid traits in the Japanese, further indicating that in genetic susceptibility to lipid levels and CAD, the related metabolic pathways are largely common across the populations, while causal variants at individual loci can be population-specific.     

Mitochondrial DNA haplogroup background affects LHON, but not suspected LHON, in Chinese patients

Recent studies have shown that mtDNA background could affect the clinical expression of Leber hereditary optic neuropathy (LHON). We analyzed the mitochondrial DNA (mtDNA) variation of 304 Chinese patients with m.11778G>A (sample #1) and of 843 suspected LHON patients who lack the three primary mutations (sample #2) to discern mtDNA haplogroup effect on disease onset. Haplogroup frequencies in the patient group was compared to frequencies in the general Han Chinese population (n = 1,689; sample #3). The overall matrilineal composition of the suspected LHON population resembles that of the general Han Chinese population, suggesting no association with mtDNA haplogroup. In contrast, analysis of these LHON patients confirms mtDNA haplogroup effect on LHON. Specifically, the LHON sample significantly differs from the general Han Chinese and suspected LHON populations by harboring an extremely lower frequency of haplogroup R9, in particular of its main sub-haplogroup F (#1 vs. #3, P-value = 1.46×10⁻¹⁷, OR = 0.051, 95% CI: 0.016–0.162; #1 vs. #2, P-value = 4.44×10⁻¹⁷, OR = 0.049, 95% CI: 0.015–0.154; in both cases, adjusted P-value <10⁻⁵) and higher frequencies of M7b (#1 vs. #3, adjusted P-value = 0.001 and #1 vs. #2, adjusted P-value = 0.004). Our result shows that mtDNA background affects LHON in Chinese patients with m.11778G>A but not suspected LHON. Haplogroup F has a protective effect against LHON, while M7b is a risk factor.     

Genetic loci associated with plasma phospholipid n-3 fatty acids: a meta-analysis of genome-wide association studies from the CHARGE Consortium

Long-chain n-3 polyunsaturated fatty acids (PUFAs) can derive from diet or from α-linolenic acid (ALA) by elongation and desaturation. We investigated the association of common genetic variation with plasma phospholipid levels of the four major n-3 PUFAs by performing genome-wide association studies in five population-based cohorts comprising 8,866 subjects of European ancestry. Minor alleles of SNPs in FADS1 and FADS2 (desaturases) were associated with higher levels of ALA (p = 3×10⁻⁶⁴) and lower levels of eicosapentaenoic acid (EPA, p = 5×10⁻⁵⁸) and docosapentaenoic acid (DPA, p = 4×10⁻¹⁵⁴). Minor alleles of SNPs in ELOVL2 (elongase) were associated with higher EPA (p = 2×10⁻¹²) and DPA (p = 1×10⁻⁴³) and lower docosahexaenoic acid (DHA, p = 1×10⁻¹⁵). In addition to genes in the n-3 pathway, we identified a novel association of DPA with several SNPs in GCKR (glucokinase regulator, p = 1×10⁻⁸). We observed a weaker association between ALA and EPA among carriers of the minor allele of a representative SNP in FADS2 (rs1535), suggesting a lower rate of ALA-to-EPA conversion in these subjects. In samples of African, Chinese, and Hispanic ancestry, associations of n-3 PUFAs were similar with a representative SNP in FADS1 but less consistent with a representative SNP in ELOVL2. Our findings show that common variation in n-3 metabolic pathway genes and in GCKR influences plasma phospholipid levels of n-3 PUFAs in populations of European ancestry and, for FADS1, in other ancestries.     

Variants from GIPR, TCF7L2, DGKB, MADD, CRY2, GLIS3, PROX1, SLC30A8 and IGF1 are associated with glucose metabolism in the Chinese

Background

Recent meta-analysis of genome-wide association studies in European descent samples identified novel loci influencing glucose and insulin related traits. In the current study, we aimed to evaluate the association between these loci and traits related to glucose metabolism in the Chinese.

Methods/Principal Findings

We genotyped seventeen single nucleotide polymorphisms (SNPs) from fifteen loci including GIPR, ADCY5, TCF7L2, VPS13C, DGKB, MADD, ADRA2A, FADS1, CRY2, SLC2A2, GLIS3, PROX1, C2CD4B, SLC30A8 and IGF1 in 6,822 Shanghai Chinese Hans comprising 3,410 type 2 diabetic patients and 3,412 normal glucose regulation subjects. MADD rs7944584 showed strong association to type 2 diabetes (p = 3.5×10⁻⁶, empirical p = 0.0002) which was not observed in the European descent populations. SNPs from GIPR, TCF7L2, CRY2, GLIS3 and SLC30A8 were also associated with type 2 diabetes (p = 0.0487∼2.0×10⁻⁸). Further adjusting age, gender and BMI as confounders found PROX1 rs340874 was associated with type 2 diabetes (p = 0.0391). SNPs from DGKB, MADD and SLC30A8 were associated with fasting glucose while PROX1 rs340874 was significantly associated with OGTT 2-h glucose (p = 0.0392∼0.0014, adjusted for age, gender and BMI), the glucose-raising allele also showed association to lower insulin secretion. IGF1 rs35767 showed significant association to both fasting and 2-h insulin levels as well as insulin secretion and sensitivity indices (p = 0.0160∼0.0035, adjusted for age, gender and BMI).

Conclusions/Significance

Our results indicated that SNPs from GIPR, TCF7L2, DGKB, MADD, CRY2, GLIS3, PROX1, SLC30A8 and IGF1 were associated with traits related to glucose metabolism in the Chinese population.     

Genetic Determinants of Long-Term Changes in Blood Lipid Concentrations: 10-Year Follow-Up of the GLACIER Study

Recent genome-wide meta-analyses identified 157 loci associated with cross-sectional lipid traits. Here we tested whether these loci associate (singly and in trait-specific genetic risk scores [GRS]) with longitudinal changes in total cholesterol (TC) and triglyceride (TG) levels in a population-based prospective cohort from Northern Sweden (the GLACIER Study). We sought replication in a southern Swedish cohort (the MDC Study; N = 2,943). GLACIER Study participants (N = 6,064) were genotyped with the MetaboChip array. Up to 3,495 participants had 10-yr follow-up data available in the GLACIER Study. The TC- and TG-specific GRSs were strongly associated with change in lipid levels (β = 0.02 mmol/l per effect allele per decade follow-up, P = 2.0×10⁻¹¹ for TC; β = 0.02 mmol/l per effect allele per decade follow-up, P = 5.0×10⁻⁵ for TG). In individual SNP analysis, one TC locus, apolipoprotein E (APOE) rs4420638 (β = 0.12 mmol/l per effect allele per decade follow-up, P = 2.0×10⁻⁵), and two TG loci, tribbles pseudokinase 1 (TRIB1) rs2954029 (β = 0.09 mmol/l per effect allele per decade follow-up, P = 5.1×10⁻⁴) and apolipoprotein A-I (APOA1) rs6589564 (β = 0.31 mmol/l per effect allele per decade follow-up, P = 1.4×10⁻⁸), remained significantly associated with longitudinal changes for the respective traits after correction for multiple testing. An additional 12 loci were nominally associated with TC or TG changes. In replication analyses, the APOE rs4420638, TRIB1 rs2954029, and APOA1 rs6589564 associations were confirmed (P≤0.001). In summary, trait-specific GRSs are robustly associated with 10-yr changes in lipid levels and three individual SNPs were strongly associated with 10-yr changes in lipid levels.     

Knowledge-driven analysis identifies a gene-gene interaction affecting high-density lipoprotein cholesterol levels in multi-ethnic populations

Total cholesterol, low-density lipoprotein cholesterol, triglyceride, and high-density lipoprotein cholesterol (HDL-C) levels are among the most important risk factors for coronary artery disease. We tested for gene–gene interactions affecting the level of these four lipids based on prior knowledge of established genome-wide association study (GWAS) hits, protein–protein interactions, and pathway information. Using genotype data from 9,713 European Americans from the Atherosclerosis Risk in Communities (ARIC) study, we identified an interaction between HMGCR and a locus near LIPC in their effect on HDL-C levels (Bonferroni corrected P _c = 0.002). Using an adaptive locus-based validation procedure, we successfully validated this gene–gene interaction in the European American cohorts from the Framingham Heart Study (P _c = 0.002) and the Multi-Ethnic Study of Atherosclerosis (MESA; P _c = 0.006). The interaction between these two loci is also significant in the African American sample from ARIC (P _c = 0.004) and in the Hispanic American sample from MESA (P _c = 0.04). Both HMGCR and LIPC are involved in the metabolism of lipids, and genome-wide association studies have previously identified LIPC as associated with levels of HDL-C. However, the effect on HDL-C of the novel gene–gene interaction reported here is twice as pronounced as that predicted by the sum of the marginal effects of the two loci. In conclusion, based on a knowledge-driven analysis of epistasis, together with a new locus-based validation method, we successfully identified and validated an interaction affecting a complex trait in multi-ethnic populations.     

Variants in GLIS3 and CRY2 are associated with type 2 diabetes and impaired fasting glucose in Chinese Hans

Background

Recent genome-wide association studies have identified a number of common variants associated with fasting glucose homeostasis and type 2 diabetes in populations of European origin. This is a replication study to examine whether such associations are also observed in Chinese Hans.

Methods

We genotyped nine variants in or near MADD, ADRA2A, CRY2, GLIS3, PROX1, FADS1, C2CD4B, IGF1 and IRS1 in a population-based cohort including 3,210 unrelated Chinese Hans from Beijing and Shanghai.

Results

We confirmed the associations of GLIS3-rs7034200 with fasting glucose (beta = 0.07 mmol/l, P = 0.03), beta cell function (HOMA-B) (beta = −3.03%, P = 0.009), and type 2 diabetes (OR [95%CI]  = 1.27 [1.09–1.49], P = 0.003) after adjustment for age, sex, region and BMI. The association for type 2 diabetes remained significant after adjusting for other diabetes related risk factors including family history of diabetes, lipid profile, medication information, hypertension and life style factors, while further adjustment for HOMA-B abolished the association. The A-allele of CRY2-rs11605924 was moderately associated with increased risk of combined IFG/type 2 diabetes (OR [95%CI]  = 1.15[1.01–1.30], P = 0.04). SNPs in or near MADD, ADRA2A, PROX1, FADS1, C2CD4B, IGF1, and IRS1 did not exhibit significant associations with type 2 diabetes or related glycemic traits (P≥0.10).

Conclusions

In conclusion, our results indicate the associations of GLIS3 locus with type 2 diabetes and impaired fasting glucose in Chinese Hans, partially mediated through impaired beta-cell function. In addition, we also found modest evidence for the association of CRY2-rs11605924 with combined IFG/type 2 diabetes.     

Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10⁻⁸–1.2×10⁻⁴³). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10⁻⁴). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10⁻³, n = 22,044), increased triglycerides (p = 2.6×10⁻¹⁴, n = 93,440), increased waist-to-hip ratio (p = 1.8×10⁻⁵, n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10⁻³, n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10⁻¹³, n = 96,748) and decreased BMI (p = 1.4×10⁻⁴, n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.     

European bone mineral density loci are also associated with BMD in East-Asian populations

Most genome-wide association (GWA) studies have focused on populations of European ancestry with limited assessment of the influence of the sequence variants on populations of other ethnicities. To determine whether markers that we have recently shown to associate with Bone Mineral Density (BMD) in Europeans also associate with BMD in East-Asians we analysed 50 markers from 23 genomic loci in samples from Korea (n = 1,397) and two Chinese Hong Kong sample sets (n = 3,869 and n = 785). Through this effort we identified fourteen loci that associated with BMD in East-Asian samples using a false discovery rate (FDR) of 0.05; 1p36 (ZBTB40, P = 4.3×10⁻⁹), 1p31 (GPR177, P = 0.00012), 3p22 (CTNNB1, P = 0.00013), 4q22 (MEPE, P = 0.0026), 5q14 (MEF2C, P = 1.3×10⁻⁵), 6q25 (ESR1, P = 0.0011), 7p14 (STARD3NL, P = 0.00025), 7q21 (FLJ42280, P = 0.00017), 8q24 (TNFRSF11B, P = 3.4×10⁻⁵), 11p15 (SOX6, P = 0.00033), 11q13 (LRP5, P = 0.0033), 13q14 (TNFSF11, P = 7.5×10⁻⁵), 16q24 (FOXL1, P = 0.0010) and 17q21 (SOST, P = 0.015). Our study marks an early effort towards the challenge of cataloguing bone density variants shared by many ethnicities by testing BMD variants that have been established in Europeans, in East-Asians.     

Replication of genetic associations with plasma lipoprotein traits in a multiethnic sample

Recent genome-wide association studies (GWAS) have reproducibly identified loci associated with plasma triglycerides (TG), HDL cholesterol, and LDL cholesterol. We sought to replicate these findings in a multiethnic population-based cohort using the curated single nucleotide polymorphism (SNP) set found on the new Illumina cardiovascular disease (CVD) beadchip, which contains approximately 50,000 SNPs densely mapping approximately 2,100 genes, selected based on their potential role in CVD. The sample consisted of individuals with European (n = 272), South Asian (n = 330), and Chinese (n = 304) ancestry. Identity by state clustering successfully classified individuals according to self-reported ethnicities. Associations between TG and APOA5, TG and LPL, HDL and CETP, and LDL and APOE were all identified (P < 2 × 10⁻⁶). In 13 loci, associations with the same SNP or a proxy SNP were identified in the same direction as previously reported (P < 0.05). Assessing the cumulative number of risk-associated alleles at multiple replicated SNPs increased the proportion of explained lipoprotein variance over and above traditional variables such as age, sex, body mass index, and ethnicity. The findings indicate the potential utility of the Illumina CVD beadchip, but they underscore the need to consider meta-analysis of results from commonly studied clinical or epidemiological samples.     

Associations between Apolipoprotein E Genotype, Diet, Body Mass Index, and Serum Lipids in Lithuanian Adult Population

Background

Apolipoprotein E (APOE) polymorphism is associated with lipid levels. Some studies have reported that blood lipid response to diet or obesity varies depending on APOE genotypes. The aim of this study was to assess the effect of APOE genotypes, the intake of saturated fatty acids (SFA), and obesity on serum lipid levels in Lithuanian adult population.

Methodology/Principal Findings

A cross-sectional health survey was carried out in five municipalities of Lithuania. The random sample was obtained from lists of 25–64 year-old inhabitants registered at primary health care centres. The data from 996 subjects (416 men and 580 women) were analysed in this study. Two single-nucleotide polymorphisms (rs429358 and rs7412) were assessed using a real-time polymerase chain reaction. 24-hour recall and food frequency questionnaire were used for evaluation of dietary habits. Serum lipids were determined using enzymatic methods.Men and women with the APOE2 genotype had the lowest level of total cholesterol (TC) (p = 0.002 for men, and p = 0.02 for women) and low-density lipoprotein cholesterol (LDL-C) (p<0.001). Multivariate linear regression analysis showed that age, genotype APOE2, SFA intake, and body mass index (BMI) were significant determinants of TC and LDL-C level (with p values ranging from 0.043 to 0.001). Our data did not reveal any statistically significant interactions between APOE genotype and SFA intake or between APOE genotype and BMI regarding TC and LDL-C level (all p>0.05). However, the predictive power of the regression model for LDL-C improved when gene-BMI interaction and gene-BMI interaction plus gene-nutrient interaction were added (p = 0.04 and p = 0.032 for R² change, respectively).

Conclusions/Significance

APOE genotypes, SFA intake, and obesity were found to be associated with blood lipid levels in Lithuanian adult population. Analysis of gene-diet and gene-obesity interactions did not confirm that the effects of diet and obesity on TC and LDL-C level significantly depended on APOE genotype.     

A genome-wide scan of Ashkenazi Jewish Crohn's disease suggests novel susceptibility loci

Crohn''s disease (CD) is a complex disorder resulting from the interaction of intestinal microbiota with the host immune system in genetically susceptible individuals. The largest meta-analysis of genome-wide association to date identified 71 CD–susceptibility loci in individuals of European ancestry. An important epidemiological feature of CD is that it is 2–4 times more prevalent among individuals of Ashkenazi Jewish (AJ) descent compared to non-Jewish Europeans (NJ). To explore genetic variation associated with CD in AJs, we conducted a genome-wide association study (GWAS) by combining raw genotype data across 10 AJ cohorts consisting of 907 cases and 2,345 controls in the discovery stage, followed up by a replication study in 971 cases and 2,124 controls. We confirmed genome-wide significant associations of 9 known CD loci in AJs and replicated 3 additional loci with strong signal (p<5×10⁻⁶). Novel signals detected among AJs were mapped to chromosomes 5q21.1 (rs7705924, combined p = 2×10⁻⁸; combined odds ratio OR = 1.48), 2p15 (rs6545946, p = 7×10⁻⁹; OR = 1.16), 8q21.11 (rs12677663, p = 2×10⁻⁸; OR = 1.15), 10q26.3 (rs10734105, p = 3×10⁻⁸; OR = 1.27), and 11q12.1 (rs11229030, p = 8×10⁻⁹; OR = 1.15), implicating biologically plausible candidate genes, including RPL7, CPAMD8, PRG2, and PRG3. In all, the 16 replicated and newly discovered loci, in addition to the three coding NOD2 variants, accounted for 11.2% of the total genetic variance for CD risk in the AJ population. This study demonstrates the complementary value of genetic studies in the Ashkenazim.     

Geographic differences in genetic susceptibility to IgA nephropathy: GWAS replication study and geospatial risk analysis

IgA nephropathy (IgAN), major cause of kidney failure worldwide, is common in Asians, moderately prevalent in Europeans, and rare in Africans. It is not known if these differences represent variation in genes, environment, or ascertainment. In a recent GWAS, we localized five IgAN susceptibility loci on Chr.6p21 (HLA-DQB1/DRB1, PSMB9/TAP1, and DPA1/DPB2 loci), Chr.1q32 (CFHR3/R1 locus), and Chr.22q12 (HORMAD2 locus). These IgAN loci are associated with risk of other immune-mediated disorders such as type I diabetes, multiple sclerosis, or inflammatory bowel disease. We tested association of these loci in eight new independent cohorts of Asian, European, and African-American ancestry (N = 4,789), followed by meta-analysis with risk-score modeling in 12 cohorts (N = 10,755) and geospatial analysis in 85 world populations. Four susceptibility loci robustly replicated and all five loci were genome-wide significant in the combined cohort (P = 5×10⁻³²–3×10⁻¹⁰), with heterogeneity detected only at the PSMB9/TAP1 locus (I² = 0.60). Conditional analyses identified two new independent risk alleles within the HLA-DQB1/DRB1 locus, defining multiple risk and protective haplotypes within this interval. We also detected a significant genetic interaction, whereby the odds ratio for the HORMAD2 protective allele was reversed in homozygotes for a CFHR3/R1 deletion (P = 2.5×10⁻⁴). A seven–SNP genetic risk score, which explained 4.7% of overall IgAN risk, increased sharply with Eastward and Northward distance from Africa (r = 0.30, P = 3×10⁻¹²⁸). This model paralleled the known East–West gradient in disease risk. Moreover, the prediction of a South–North axis was confirmed by registry data showing that the prevalence of IgAN–attributable kidney failure is increased in Northern Europe, similar to multiple sclerosis and type I diabetes. Variation at IgAN susceptibility loci correlates with differences in disease prevalence among world populations. These findings inform genetic, biological, and epidemiological investigations of IgAN and permit cross-comparison with other complex traits that share genetic risk loci and geographic patterns with IgAN.     

Knowledge-driven multi-locus analysis reveals gene-gene interactions influencing HDL cholesterol level in two independent EMR-linked biobanks

Genome-wide association studies (GWAS) are routinely being used to examine the genetic contribution to complex human traits, such as high-density lipoprotein cholesterol (HDL-C). Although HDL-C levels are highly heritable (h²∼0.7), the genetic determinants identified through GWAS contribute to a small fraction of the variance in this trait. Reasons for this discrepancy may include rare variants, structural variants, gene-environment (GxE) interactions, and gene-gene (GxG) interactions. Clinical practice-based biobanks now allow investigators to address these challenges by conducting GWAS in the context of comprehensive electronic medical records (EMRs). Here we apply an EMR-based phenotyping approach, within the context of routine care, to replicate several known associations between HDL-C and previously characterized genetic variants: CETP (rs3764261, p = 1.22e-25), LIPC (rs11855284, p = 3.92e-14), LPL (rs12678919, p = 1.99e-7), and the APOA1/C3/A4/A5 locus (rs964184, p = 1.06e-5), all adjusted for age, gender, body mass index (BMI), and smoking status. By using a novel approach which censors data based on relevant co-morbidities and lipid modifying medications to construct a more rigorous HDL-C phenotype, we identified an association between HDL-C and TRIB1, a gene which previously resisted identification in studies with larger sample sizes. Through the application of additional analytical strategies incorporating biological knowledge, we further identified 11 significant GxG interaction models in our discovery cohort, 8 of which show evidence of replication in a second biobank cohort. The strongest predictive model included a pairwise interaction between LPL (which modulates the incorporation of triglyceride into HDL) and ABCA1 (which modulates the incorporation of free cholesterol into HDL). These results demonstrate that gene-gene interactions modulate complex human traits, including HDL cholesterol.     

Associations between APOE variants and metabolic traits and the impact of psychological stress

Objective

In a previous study, we observed that associations between APOE rs439401 and metabolic traits were moderated by chronic stress. Thus, in a population of stressed and non-stressed Danish men, we examined whether associations between APOE rs439401 and a panel of metabolic quantitative traits, all metabolic traits which may lead to T2D and CVD were moderated by psychological stress.

Methods

Obese young men (n = 475, BMI≥31.0 kg/m²) and a randomly selected control group (n = 709) identified from a population of 141,800 men were re-examined in two surveys (S-46: mean age 46, S-49: mean age 49 years) where anthropometric and biochemical measures were available. Psychological stress factors were assessed by a self-administered 7-item questionnaire. Each item had the possible response categories “yes” and “no” and assessed familial problems and conflicts. Summing positive responses constituted a stress item score, which was then dichotomized into stressed and non-stressed. Logistic regression analysis, applying a recessive genetic model, was used to assess odds ratios (OR) of the associations between APOE rs439401 genotypes and adverse levels of metabolic traits.

Results

The APOE rs439401 TT-genotype associated positively with BMI (OR = 1.09 [1.01; 1.17]), waist circumference (OR = 1.09 [1.02; 1.17]) in stressed men at S-46. Positive associations were observed for fasting plasma glucose (OR = 1.42 [1.07; 1.87]), serum triglycerides (OR = 1.41 [1.05; 1.91]) and with fasting plasma insulin (OR = 1.48 [1.05; 2.08]) in stressed men at S-49. Rs439401 TT-genotype also associated positively with surrogate measures of insulin resistance (HOMA-IR; OR = 1.21 [1.03; 1.41]) and inversely with insulin sensitivity (Stumvoll index; OR = 0.90 [0.82; 0.99], BIGTT-S_I; OR = 0.60 [0.43; 0.85]) in stressed men. No significant associations were observed in non-stressed men, albeit the estimates showed similar but weaker trends as in stressed men.

Conclusion

The present results suggest that the APOE rs439401 TT-genotype is associated with an adverse metabolic profile in a population of psychologically stressed Danish men.     

Genetic association for renal traits among participants of African ancestry reveals new loci for renal function

Chronic kidney disease (CKD) is an increasing global public health concern, particularly among populations of African ancestry. We performed an interrogation of known renal loci, genome-wide association (GWA), and IBC candidate-gene SNP association analyses in African Americans from the CARe Renal Consortium. In up to 8,110 participants, we performed meta-analyses of GWA and IBC array data for estimated glomerular filtration rate (eGFR), CKD (eGFR <60 mL/min/1.73 m²), urinary albumin-to-creatinine ratio (UACR), and microalbuminuria (UACR >30 mg/g) and interrogated the 250 kb flanking region around 24 SNPs previously identified in European Ancestry renal GWAS analyses. Findings were replicated in up to 4,358 African Americans. To assess function, individually identified genes were knocked down in zebrafish embryos by morpholino antisense oligonucleotides. Expression of kidney-specific genes was assessed by in situ hybridization, and glomerular filtration was evaluated by dextran clearance. Overall, 23 of 24 previously identified SNPs had direction-consistent associations with eGFR in African Americans, 2 of which achieved nominal significance (UMOD, PIP5K1B). Interrogation of the flanking regions uncovered 24 new index SNPs in African Americans, 12 of which were replicated (UMOD, ANXA9, GCKR, TFDP2, DAB2, VEGFA, ATXN2, GATM, SLC22A2, TMEM60, SLC6A13, and BCAS3). In addition, we identified 3 suggestive loci at DOK6 (p-value = 5.3×10⁻⁷) and FNDC1 (p-value = 3.0×10⁻⁷) for UACR, and KCNQ1 with eGFR (p = 3.6×10⁻⁶). Morpholino knockdown of kcnq1 in the zebrafish resulted in abnormal kidney development and filtration capacity. We identified several SNPs in association with eGFR in African Ancestry individuals, as well as 3 suggestive loci for UACR and eGFR. Functional genetic studies support a role for kcnq1 in glomerular development in zebrafish.     

Association of variants at 1q32 and STAT3 with ankylosing spondylitis suggests genetic overlap with Crohn's disease

Ankylosing spondylitis (AS) is a common inflammatory arthritic condition. Overt inflammatory bowel disease (IBD) occurs in about 10% of AS patients, and in addition 70% of AS cases may have subclinical terminal ileitis. Spondyloarthritis is also common in IBD patients. We therefore tested Crohn''s disease susceptibility genes for association with AS, aiming to identify pleiotropic genetic associations with both diseases. Genotyping was carried out using Sequenom and Applied Biosystems TaqMan and OpenArray technologies on 53 markers selected from 30 Crohn''s disease associated genomic regions. We tested genotypes in a population of unrelated individual cases (n = 2,773) and controls (n = 2,215) of white European ancestry for association with AS. Statistical analysis was carried out using a Cochran-Armitage test for trend in PLINK. Strong association was detected at chr1q32 near KIF21B (rs11584383, P = 1.6×10⁻¹⁰, odds ratio (OR) = 0.74, 95% CI:0.68–0.82). Association with disease was also detected for 2 variants within STAT3 (rs6503695, P = 4.6×10⁻⁴. OR = 0.86 (95% CI:0.79–0.93); rs744166, P = 2.6×10⁻⁵, OR = 0.84 (95% CI:0.77–0.91)). Association was confirmed for IL23R (rs11465804, P = 1.2×10⁻⁵, OR = 0.65 (95% CI:0.54–0.79)), and further associations were detected for IL12B (rs10045431, P = 5.2×10⁻⁵, OR = 0.83 (95% CI:0.76–0.91)), CDKAL1 (rs6908425, P = 1.1×10⁻⁴, OR = 0.82 (95% CI:0.74–0.91)), LRRK2/MUC19 (rs11175593, P = 9.9×10⁻⁵, OR = 1.92 (95% CI: 1.38–2.67)), and chr13q14 (rs3764147, P = 5.9×10⁻⁴, OR = 1.19 (95% CI: 1.08–1.31)). Excluding cases with clinical IBD did not significantly affect these findings. This study identifies chr1q32 and STAT3 as ankylosing spondylitis susceptibility loci. It also further confirms association for IL23R and detects suggestive association with another 4 loci. STAT3 is a key signaling molecule within the Th17 lymphocyte differentiation pathway and further enhances the case for a major role of this T-lymphocyte subset in ankylosing spondylitis. Finally these findings suggest common aetiopathogenic pathways for AS and Crohn''s disease and further highlight the involvement of common risk variants across multiple diseases.     

A genome-wide meta-analysis of six type 1 diabetes cohorts identifies multiple associated loci

Diabetes impacts approximately 200 million people worldwide, of whom approximately 10% are affected by type 1 diabetes (T1D). The application of genome-wide association studies (GWAS) has robustly revealed dozens of genetic contributors to the pathogenesis of T1D, with the most recent meta-analysis identifying in excess of 40 loci. To identify additional genetic loci for T1D susceptibility, we examined associations in the largest meta-analysis to date between the disease and ∼2.54 million SNPs in a combined cohort of 9,934 cases and 16,956 controls. Targeted follow-up of 53 SNPs in 1,120 affected trios uncovered three new loci associated with T1D that reached genome-wide significance. The most significantly associated SNP (rs539514, P = 5.66×10⁻¹¹) resides in an intronic region of the LMO7 (LIM domain only 7) gene on 13q22. The second most significantly associated SNP (rs478222, P = 3.50×10⁻⁹) resides in an intronic region of the EFR3B (protein EFR3 homolog B) gene on 2p23; however, the region of linkage disequilibrium is approximately 800 kb and harbors additional multiple genes, including NCOA1, C2orf79, CENPO, ADCY3, DNAJC27, POMC, and DNMT3A. The third most significantly associated SNP (rs924043, P = 8.06×10⁻⁹) lies in an intergenic region on 6q27, where the region of association is approximately 900 kb and harbors multiple genes including WDR27, C6orf120, PHF10, TCTE3, C6orf208, LOC154449, DLL1, FAM120B, PSMB1, TBP, and PCD2. These latest associated regions add to the growing repertoire of gene networks predisposing to T1D.