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.     

A genome-wide association scan on the levels of markers of inflammation in Sardinians reveals associations that underpin its complex regulation

Identifying the genes that influence levels of pro-inflammatory molecules can help to elucidate the mechanisms underlying this process. We first conducted a two-stage genome-wide association scan (GWAS) for the key inflammatory biomarkers Interleukin-6 (IL-6), the general measure of inflammation erythrocyte sedimentation rate (ESR), monocyte chemotactic protein-1 (MCP-1), and high-sensitivity C-reactive protein (hsCRP) in a large cohort of individuals from the founder population of Sardinia. By analysing 731,213 autosomal or X chromosome SNPs and an additional ∼1.9 million imputed variants in 4,694 individuals, we identified several SNPs associated with the selected quantitative trait loci (QTLs) and replicated all the top signals in an independent sample of 1,392 individuals from the same population. Next, to increase power to detect and resolve associations, we further genotyped the whole cohort (6,145 individuals) for 293,875 variants included on the ImmunoChip and MetaboChip custom arrays. Overall, our combined approach led to the identification of 9 genome-wide significant novel independent signals—5 of which were identified only with the custom arrays—and provided confirmatory evidence for an additional 7. Novel signals include: for IL-6, in the ABO gene (rs657152, p = 2.13×10⁻²⁹); for ESR, at the HBB (rs4910472, p = 2.31×10⁻¹¹) and UCN119B/SPPL3 (rs11829037, p = 8.91×10⁻¹⁰) loci; for MCP-1, near its receptor CCR2 (rs17141006, p = 7.53×10⁻¹³) and in CADM3 (rs3026968, p = 7.63×10⁻¹³); for hsCRP, within the CRP gene (rs3093077, p = 5.73×10⁻²¹), near DARC (rs3845624, p = 1.43×10⁻¹⁰), UNC119B/SPPL3 (rs11829037, p = 1.50×10⁻¹⁴), and ICOSLG/AIRE (rs113459440, p = 1.54×10⁻⁰⁸) loci. Confirmatory evidence was found for IL-6 in the IL-6R gene (rs4129267); for ESR at CR1 (rs12567990) and TMEM57 (rs10903129); for MCP-1 at DARC (rs12075); and for hsCRP at CRP (rs1205), HNF1A (rs225918), and APOC-I (rs4420638). Our results improve the current knowledge of genetic variants underlying inflammation and provide novel clues for the understanding of the molecular mechanisms regulating this complex process.     

The role of inflammatory pathway genetic variation on maternal metabolic phenotypes during pregnancy

Background

Since mediators of inflammation are associated with insulin resistance, and the risk of developing diabetes mellitus and gestational diabetes, we hypothesized that genetic variation in members of the inflammatory gene pathway impact glucose levels and related phenotypes in pregnancy. We evaluated this hypothesis by testing for association between genetic variants in 31 inflammatory pathway genes in the Hyperglycemia and Adverse Pregnancy Outcome (HAPO) cohort, a large multiethnic multicenter study designed to address the impact of glycemia less than overt diabetes on pregnancy outcome.

Results

Fasting, 1-hour, and 2-hour glucose, fasting and 1-hour C-peptide, and HbA1c levels were measured in blood samples obtained from HAPO participants during an oral glucose tolerance test at 24-32 weeks gestation. We tested for association between 458 SNPs mapping to 31 genes in the inflammatory pathway and metabolic phenotypes in 3836 European ancestry and 1713 Thai pregnant women. The strongest evidence for association was observed with TNF alpha and HbA1c (rs1052248; 0.04% increase per allele C; p-value = 4.4×10⁻⁵), RETN and fasting plasma glucose (rs1423096; 0.7 mg/dl decrease per allele A; p-value = 1.1×10⁻⁴), IL8 and 1 hr plasma glucose (rs2886920; 2.6 mg/dl decrease per allele T; p-value = 1.3×10⁻⁴), ADIPOR2 and fasting C-peptide (rs2041139; 0.55 ug/L decrease per allele A; p-value = 1.4×10⁻⁴), LEPR and 1-hour C-peptide (rs1171278; 0.62 ug/L decrease per allele T; p-value = 2.4×10⁻⁴), and IL6 and 1-hour plasma glucose (rs6954897; −2.29 mg/dl decrease per allele G, p-value = 4.3×10⁻⁴).

Conclusions

Based on the genes surveyed in this study the inflammatory pathway is unlikely to have a strong impact on maternal metabolic phenotypes in pregnancy although variation in individual members of the pathway (e.g. RETN, IL8, ADIPOR2, LEPR, IL6, and TNF alpha,) may contribute to metabolic phenotypes in pregnant women.     

Associations of CFH Polymorphisms and CFHR1-CFHR3 Deletion with Blood Pressure and Hypertension in Chinese Population

Dysregulation of the complement system has been linked to pathogenesis of hypertension. However, whether genetic changes of complement factor H (CFH) and its related genes are associated with hypertension is unknown. We genotyped three SNPs in the CFH gene cluster that are closely linked to age-related macular degeneration, namely rs1061170 (Y402H), rs2274700 (A473A) and rs7542235 (CFHR1–3Δ), and tested for their associations with blood pressure and hypertension risk in a population-based cohort including 3,210 unrelated Chinese Hans (50–70 years of age) from Beijing and Shanghai. We found that rs2274700 (A473A) and rs7542235 (CFHR1–3Δ) were both significantly associated with diastolic blood pressure (DBP) (β = 0.632–1.431, P≤0.038) and systolic blood pressure (SBP) (β = 1.567–4.445, P≤0.008), and rs2274700 (A473A) was associated with hypertension risk (OR [95%CI]: 1.175 [1.005–1.373], P = 0.048). Notably, the associations of rs2274700 (A473A) with DBP (P = 2.1×10⁻³), SBP (P = 8×10⁻⁵) and hypertension risk (P = 7.9×10⁻³) were significant only in the individuals with low CRP levels (<2.0 mg/l), but not in those with CRP levels ≥2.0 mg/l (P≥0.0807) (P for interaction ≤0.0467). However, no significant association between rs1061170 (Y402H) and blood pressure or hypertension risk was observed (P≥0.259). In conclusion, our results suggest that genetic variations in CFH and its related genes may contribute to hypertension risk in Chinese Hans.     

CFH, VEGF, and PEDF genotypes and the response to intravitreous injection of bevacizumab for the treatment of age-related macular degeneration

We determined whether there is an association between complement factor H (CFH), high-temperature requirement A-1 (HTRA1), vascular endothelial growth factor (VEGF), and pigment epithelium-derived factor (PEDF) genotypes and the response to treatment with a single intravitreous injection of bevacizumab for age-related macular degeneration (AMD). Eighty-three patients with exudative AMD treated by bevacizumab injection were genotyped for three single nucleotide polymorphisms (SNPs; rs800292, rs1061170, rs1410996) in the CFH gene, a rs11200638-SNP in the HTRA1 gene, three SNPs (rs699947, rs1570360, rs2010963) in the VEGF gene, and four SNPs (rs12150053, rs12948385, rs9913583, rs1136287) in the PEDF gene using a TaqMan assay. The CT genotype (heterozygous) of CFH-rs1061170 was more frequently represented in nonresponders in vision than TT genotypes (nonrisk allele homozygous) at the time points of 1 and 3 months, while there was no CC genotype (risk allele homozygous) in our study cohort (p = 7.66 × 10⁻³, 7.83 × 10⁻³, respectively). VEGF-rs699947 was also associated with vision changes at 1 month and PEDF-rs1136287 at 3 months (p = 5.11 × 10⁻³, 2.05 × 10⁻², respectively). These variants may be utilized for genetic biomarkers to estimate visual outcomes in the response to intravitreal bevacizumab treatment for AMD.     

Common variants in CRP and LEPR influence high sensitivity C-reactive protein levels in North Indians

Background

High sensitivity C-reactive protein (hsCRP) levels are shown to be influenced by genetic variants in Europeans; however, little is explored in Indian population.

Methods

Herein, we comprehensively evaluated association of all previously reported genetic determinants of hsCRP levels, including 18 cis (proximal to CRP gene) and 73 trans-acting (distal to CRP gene) variants in 4,200 North Indians of Indo-European ethnicity. First, we evaluated association of 91 variants from 12 candidate loci with hsCRP levels in 2,115 North Indians (1,042 non-diabetic subjects and 1,073 patients with type 2 diabetes). Then, cis and trans-acting variants contributing maximally to hsCRP level variation were further replicated in an independent 2,085 North Indians (1,047 patients with type 2 diabetes and 1,038 non-diabetic subjects).

Results

We found association of 12 variants from CRP, LEPR, IL1A, IL6, and IL6R with hsCRP levels in non-diabetic subjects. However, only rs3093059-CRP [β = 0.33, P = 9.6×10⁻⁵] and the haplotype harboring rs3093059 risk allele [β = 0.32 µg/mL, P = 1.4×10⁻⁴/P_perm = 9.0×10⁻⁴] retained significance after correcting for multiple testing. The cis-acting variant rs3093059-CRP had maximum contribution to the variance in hsCRP levels (1.14%). Among, trans-acting variants, rs1892534-LEPR was observed to contribute maximally to hsCRP level variance (0.59%). Associations of rs3093059-CRP and rs1892534-LEPR were confirmed by replication and attained higher significance after meta-analysis [β_meta = 0.26/0.22; P_meta = 4.3×10⁻⁷/7.4×10⁻³ and β_meta = −0.15/−0.12; P_meta = 2.0×10⁻⁶/1.6×10⁻⁶ for rs3093059 and rs1892534, respectively in non-diabetic subjects and all subjects taken together].

Conclusion

In conclusion, we identified rs3093059 in CRP and rs1892534 in LEPR as major cis and trans-acting contributor respectively, to the variance in hsCRP levels in North Indian population.     

A genome-wide association search for type 2 diabetes genes in African Americans

African Americans are disproportionately affected by type 2 diabetes (T2DM) yet few studies have examined T2DM using genome-wide association approaches in this ethnicity. The aim of this study was to identify genes associated with T2DM in the African American population. We performed a Genome Wide Association Study (GWAS) using the Affymetrix 6.0 array in 965 African-American cases with T2DM and end-stage renal disease (T2DM-ESRD) and 1029 population-based controls. The most significant SNPs (n = 550 independent loci) were genotyped in a replication cohort and 122 SNPs (n = 98 independent loci) were further tested through genotyping three additional validation cohorts followed by meta-analysis in all five cohorts totaling 3,132 cases and 3,317 controls. Twelve SNPs had evidence of association in the GWAS (P<0.0071), were directionally consistent in the Replication cohort and were associated with T2DM in subjects without nephropathy (P<0.05). Meta-analysis in all cases and controls revealed a single SNP reaching genome-wide significance (P<2.5×10⁻⁸). SNP rs7560163 (P = 7.0×10⁻⁹, OR (95% CI) = 0.75 (0.67–0.84)) is located intergenically between RND3 and RBM43. Four additional loci (rs7542900, rs4659485, rs2722769 and rs7107217) were associated with T2DM (P<0.05) and reached more nominal levels of significance (P<2.5×10⁻⁵) in the overall analysis and may represent novel loci that contribute to T2DM. We have identified novel T2DM-susceptibility variants in the African-American population. Notably, T2DM risk was associated with the major allele and implies an interesting genetic architecture in this population. These results suggest that multiple loci underlie T2DM susceptibility in the African-American population and that these loci are distinct from those identified in other ethnic populations.     

Sequence Variants in Three Loci Influence Monocyte Counts and Erythrocyte Volume

Blood cells participate in vital physiological processes, and their numbers are tightly regulated so that homeostasis is maintained. Disruption of key regulatory mechanisms underlies many blood-related Mendelian diseases but also contributes to more common disorders, including atherosclerosis. We searched for quantitative trait loci (QTL) for hematology traits through a whole-genome association study, because these could provide new insights into both hemopoeitic and disease mechanisms. We tested 1.8 million variants for association with 13 hematology traits measured in 6015 individuals from the Australian and Dutch populations. These traits included hemoglobin composition, platelet counts, and red blood cell and white blood cell indices. We identified three regions of strong association that, to our knowledge, have not been previously reported in the literature. The first was located in an intergenic region of chromosome 9q31 near LPAR1, explaining 1.5% of the variation in monocyte counts (best SNP rs7023923, p = 8.9 × 10⁻¹⁴). The second locus was located on chromosome 6p21 and associated with mean cell erythrocyte volume (rs12661667, p = 1.2 × 10⁻⁹, 0.7% variance explained) in a region that spanned five genes, including CCND3, a member of the D-cyclin gene family that is involved in hematopoietic stem cell expansion. The third region was also associated with erythrocyte volume and was located in an intergenic region on chromosome 6q24 (rs592423, p = 5.3 × 10⁻⁹, 0.6% variance explained). All three loci replicated in an independent panel of 1543 individuals (p values = 0.001, 9.9 × 10⁻⁵, and 7 × 10⁻⁵, respectively). The identification of these QTL provides new opportunities for furthering our understanding of the mechanisms regulating hemopoietic cell fate.     

Identification of novel susceptibility Loci for kawasaki disease in a Han chinese population by a genome-wide association study

Kawasaki disease (KD) is an acute systemic vasculitis syndrome that primarily affects infants and young children. Its etiology is unknown; however, epidemiological findings suggest that genetic predisposition underlies disease susceptibility. Taiwan has the third-highest incidence of KD in the world, after Japan and Korea. To investigate novel mechanisms that might predispose individuals to KD, we conducted a genome-wide association study (GWAS) in 250 KD patients and 446 controls in a Han Chinese population residing in Taiwan, and further validated our findings in an independent Han Chinese cohort of 208 cases and 366 controls. The most strongly associated single-nucleotide polymorphisms (SNPs) detected in the joint analysis corresponded to three novel loci. Among these KD-associated SNPs three were close to the COPB2 (coatomer protein complex beta-2 subunit) gene: rs1873668 (p = 9.52×10⁻⁵), rs4243399 (p = 9.93×10⁻⁵), and rs16849083 (p = 9.93×10⁻⁵). We also identified a SNP in the intronic region of the ERAP1 (endoplasmic reticulum amino peptidase 1) gene (rs149481, p_best = 4.61×10⁻⁵). Six SNPs (rs17113284, rs8005468, rs10129255, rs2007467, rs10150241, and rs12590667) clustered in an area containing immunoglobulin heavy chain variable regions genes, with p_best-values between 2.08×10⁻⁵ and 8.93×10⁻⁶, were also identified. This is the first KD GWAS performed in a Han Chinese population. The novel KD candidates we identified have been implicated in T cell receptor signaling, regulation of proinflammatory cytokines, as well as antibody-mediated immune responses. These findings may lead to a better understanding of the underlying molecular pathogenesis of KD.     

The PNPLA3 rs738409 G-allele associates with reduced fasting serum triglyceride and serum cholesterol in Danes with impaired glucose regulation

Background and Aim

Non-alcoholic fatty liver disease (NAFLD) is a common condition, associated with hepatic insulin resistance and the metabolic syndrome including hyperglycaemia and dyslipidemia. We aimed at studying the potential impact of the NAFLD-associated PNPLA3 rs738409 G-allele on NAFLD-related metabolic traits in hyperglycaemic individuals.

Methods

The rs738409 variant was genotyped in the population-based Inter99 cohort examined by an oral glucose-tolerance test, and a combined study-sample consisting of 192 twins (96 twin pairs) and a sub-set of the Inter99 population (n = 63) examined by a hyperinsulinemic euglycemic clamp (n _total = 255). In Inter99, we analyzed associations of rs738409 with components of the WHO-defined metabolic syndrome (n = 5,847) and traits related to metabolic disease (n = 5,663). In the combined study sample we elucidated whether the rs738409 G-allele altered hepatic or peripheral insulin sensitivity. Study populations were divided into individuals with normal glucose-tolerance (NGT) and with impaired glucose regulation (IGR).

Results

The case-control study showed no associations with components of the metabolic syndrome or the metabolic syndrome. Among 1,357 IGR individuals, the rs738409 G-allele associated with decreased fasting serum triglyceride levels (per allele effect(β) = −9.9% [−14.4%;−4.0% (95% CI)], p = 5.1×10⁻⁵) and fasting total cholesterol (β = −0.2 mmol/l [−0.3;−0.01 mmol/l(95% CI)], p = 1.5×10⁻⁴). Meta-analyses showed no impact on hepatic or peripheral insulin resistance in carriers of the rs738409 G-allele.

Conclusion

Our findings suggest that the G-allele of PNPLA3 rs738409 associates with reduced fasting levels of cholesterol and triglyceride in individuals with IGR.     

Functional Analysis of Two PLA2G2A Variants Associated with Secretory Phospholipase A2-IIA Levels

Background

Secretory phospholipase A2 group IIA (sPLA2-IIA) has been identified as a biomarker of atherosclerosis in observational and animal studies. The protein is encoded by the PLA2G2A gene and the aim of this study was to test the functionality of two PLA2G2A non-coding SNPs, rs11573156 C>G and rs3767221 T>G where the rare alleles have been previously associated with higher and lower sPLA2-IIA levels respectively.

Methodology/Principal Findings

Luciferase assays, electrophoretic mobility shift assays (EMSA), and RNA expression by RT-PCR were used to examine allelic differences. For rs3767221 the G allele showed ∼55% lower luciferase activity compared to the T allele (T = 62.1 (95% CI 59.1 to 65.1) G = 27.8 (95% CI 25.0 to 30.6), p = 1.22×10⁻³⁵, and stronger EMSA binding of a nuclear protein compared to the T-allele. For rs11573156 C >G there were no luciferase or EMSA allelic differences seen. In lymphocyte cell RNA, from individuals of known rs11573156 genotype, there was no allelic RNA expression difference for exons 5 and 6, but G allele carriers (n = 7) showed a trend to lower exon 1–2 expression compared to CC individuals. To take this further, in the ASAP study (n = 223), an rs11573156 proxy (r² = 0.91) showed ∼25% higher liver expression of PLA2G2A (1.67×10⁻¹⁷) associated with the G allele. However, considering exon specific expression, the association was greatly reduced for exon 2 (4.5×10⁻⁵) compared to exons 3–6 (10⁻¹⁰ to 10⁻²⁰), suggesting rs11573156 G allele-specific exon 2 skipping.

Conclusion

Both SNPs are functional and provide useful tools for Mendelian Randomisation to determine whether the relationship between sPLA2-IIA and coronary heart disease is causal.     

Genetic Modulation of Lipid Profiles following Lifestyle Modification or Metformin Treatment: The Diabetes Prevention Program

Weight-loss interventions generally improve lipid profiles and reduce cardiovascular disease risk, but effects are variable and may depend on genetic factors. We performed a genetic association analysis of data from 2,993 participants in the Diabetes Prevention Program to test the hypotheses that a genetic risk score (GRS) based on deleterious alleles at 32 lipid-associated single-nucleotide polymorphisms modifies the effects of lifestyle and/or metformin interventions on lipid levels and nuclear magnetic resonance (NMR) lipoprotein subfraction size and number. Twenty-three loci previously associated with fasting LDL-C, HDL-C, or triglycerides replicated (P = 0.04–1×10⁻¹⁷). Except for total HDL particles (r = −0.03, P = 0.26), all components of the lipid profile correlated with the GRS (partial |r| = 0.07–0.17, P = 5×10⁻⁵–1×10⁻¹⁹). The GRS was associated with higher baseline-adjusted 1-year LDL cholesterol levels (β = +0.87, SEE±0.22 mg/dl/allele, P = 8×10⁻⁵, P _interaction = 0.02) in the lifestyle intervention group, but not in the placebo (β = +0.20, SEE±0.22 mg/dl/allele, P = 0.35) or metformin (β = −0.03, SEE±0.22 mg/dl/allele, P = 0.90; P _interaction = 0.64) groups. Similarly, a higher GRS predicted a greater number of baseline-adjusted small LDL particles at 1 year in the lifestyle intervention arm (β = +0.30, SEE±0.012 ln nmol/L/allele, P = 0.01, P _interaction = 0.01) but not in the placebo (β = −0.002, SEE±0.008 ln nmol/L/allele, P = 0.74) or metformin (β = +0.013, SEE±0.008 nmol/L/allele, P = 0.12; P _interaction = 0.24) groups. Our findings suggest that a high genetic burden confers an adverse lipid profile and predicts attenuated response in LDL-C levels and small LDL particle number to dietary and physical activity interventions aimed at weight loss.     

A genome-wide association study of nephrolithiasis in the Japanese population identifies novel susceptible Loci at 5q35.3, 7p14.3, and 13q14.1

Nephrolithiasis is a common nephrologic disorder with complex etiology. To identify the genetic factor(s) for nephrolithiasis, we conducted a three-stage genome-wide association study (GWAS) using a total of 5,892 nephrolithiasis cases and 17,809 controls of Japanese origin. Here we found three novel loci for nephrolithiasis: RGS14-SLC34A1-PFN3-F12 on 5q35.3 (rs11746443; P = 8.51×10⁻¹², odds ratio (OR) = 1.19), INMT-FAM188B-AQP1 on 7p14.3 (rs1000597; P = 2.16×10⁻¹⁴, OR = 1.22), and DGKH on 13q14.1 (rs4142110; P = 4.62×10⁻⁹, OR = 1.14). Subsequent analyses in 21,842 Japanese subjects revealed the association of SNP rs11746443 with the reduction of estimated glomerular filtration rate (eGFR) (P = 6.54×10⁻⁸), suggesting a crucial role for this variation in renal function. Our findings elucidated the significance of genetic variations for the pathogenesis of nephrolithiasis.     

Genome-wide screen identifies rs646776 near sortilin as a regulator of progranulin levels in human plasma

Recent studies suggest progranulin (GRN) is a neurotrophic factor. Loss-of-function mutations in the progranulin gene (GRN) cause frontotemporal lobar degeneration (FTLD), a progressive neurodegenerative disease affecting ∼10% of early-onset dementia patients. Using an enzyme-linked immunosorbent assay, we previously showed that GRN is detectable in human plasma and can be used to predict GRN mutation status. This study also showed a wide range in plasma GRN levels in non-GRN mutation carriers, including controls. We have now performed a genome-wide association study of 313,504 single-nucleotide polymorphisms (SNPs) in 533 control samples and identified on chromosome 1p13.3 two SNPs with genome-wide significant association with plasma GRN levels (top SNP rs646776; p = 1.7 × 10⁻³⁰). The association of rs646776 with plasma GRN levels was replicated in two independent series of 508 controls (p = 1.9 × 10⁻¹⁹) and 197 FTLD patients (p = 6.4 × 10⁻¹²). Overall, each copy of the minor C allele decreased GRN levels by ∼15%. SNP rs646776 is located near sortilin (SORT1), and the minor C allele of rs646776 was previously associated with increased SORT1 mRNA levels. Supporting these findings, overexpression of SORT1 in cultured HeLa cells dramatically reduced GRN levels in the conditioned media, whereas knockdown of SORT1 increased extracellular GRN levels. In summary, we identified significant association of a locus on chromosome 1p13.3 with plasma GRN levels through an unbiased genome-wide screening approach and implicated SORT1 as an important regulator of GRN levels. This finding opens avenues for future research into GRN biology and the pathophysiology of neurodegenerative diseases.     

The role of osteopontin (OPN/SPP1) haplotypes in the susceptibility to Crohn's disease

Background

Osteopontin represents a multifunctional molecule playing a pivotal role in chronic inflammatory and autoimmune diseases. Its expression is increased in inflammatory bowel disease (IBD). The aim of our study was to analyze the association of osteopontin (OPN/SPP1) gene variants in a large cohort of IBD patients.

Methodology/Principal Findings

Genomic DNA from 2819 Caucasian individuals (n = 841 patients with Crohn''s disease (CD), n = 473 patients with ulcerative colitis (UC), and n = 1505 healthy unrelated controls) was analyzed for nine OPN SNPs (rs2728127, rs2853744, rs11730582, rs11739060, rs28357094, rs4754 = p.Asp80Asp, rs1126616 = p.Ala236Ala, rs1126772 and rs9138). Considering the important role of osteopontin in Th17-mediated diseases, we performed analysis for epistasis with IBD-associated IL23R variants and analyzed serum levels of the Th17 cytokine IL-22. For four OPN SNPs (rs4754, rs1126616, rs1126772 and rs9138), we observed significantly different distributions between male and female CD patients. rs4754 was protective in male CD patients (p = 0.0004, OR = 0.69). None of the other investigated OPN SNPs was associated with CD or UC susceptibility. However, several OPN haplotypes showed significant associations with CD susceptibility. The strongest association was found for a haplotype consisting of the 8 OPN SNPs rs2728127-rs2853744-rs11730582-rs11439060-rs28357094-rs112661-rs1126772-rs9138 (omnibus p-value = 2.07×10⁻⁸). Overall, the mean IL-22 secretion in the combined group of OPN minor allele carriers with CD was significantly lower than that of CD patients with OPN wildtype alleles (p = 3.66×10⁻⁵). There was evidence for weak epistasis between the OPN SNP rs28357094 with the IL23R SNP rs10489629 (p = 4.18×10⁻²) and between OPN SNP rs1126616 and IL23R SNP rs2201841 (p = 4.18×10⁻²) but none of these associations remained significant after Bonferroni correction.

Conclusions/Significance

Our study identified OPN haplotypes as modifiers of CD susceptibility, while the combined effects of certain OPN variants may modulate IL-22 secretion.     

A Weighted False Discovery Rate Control Procedure Reveals Alleles at FOXA2 that Influence Fasting Glucose Levels

Association signals in GWAS are usually prioritized solely by p values. Here, we attempt to improve the power of GWAS by using a weighted false discovery rate control procedure to detect associations of low-frequency variants with effect sizes similar to or even larger than those of common variants. We used the Affymetrix Genome-Wide Human SNP Array 6.0 to test for association with fasting glucose levels in the Atherosclerosis Risk in Communities Study (ARIC) population. In addition to finding several previously identified sequence variations, we identified a low-frequency variant (rs1209523; minor allele frequency = 0.043) near FOXA2 that was associated with fasting glucose levels in European Americans (EAs) (n = 7428, p value = 1.3 × 10⁻⁵). The association between rs1209523 and glucose levels was also significant in African Americans (AAs) (n = 2029, p value = 6.7 × 10⁻³) of the ARIC and was confirmed by replication in both EAs and AAs of the Dallas Heart Study (n = 963 and 1571, respectively; p values = 5.3 × 10⁻³ and 5.8 × 10⁻⁴, respectively) and in EAs of the Cooper Center Longitudinal Study (n = 2862; p value = 1.6 × 10⁻²). A meta-analysis of these five populations yielded an estimated effect size of −1.31 mg/dl per minor allele (p value = 2.2 × 10⁻¹¹). This study reveals that there is a cache of less-frequent variants in GWAS arrays that can be identified via analytical approaches accounting for allele frequencies.