Systematic Evaluation of Genetic Variants for Polycystic Ovary Syndrome in a Chinese Population

To date, eleven genome-wide significant (GWS) loci (P < 5×10⁻⁸) for polycystic ovary syndrome (PCOS) have been identified through genome-wide association studies (GWAS). Some of the risk loci have been selected for replications and validated in multiple ethnicities, however, few previous studies investigated all loci. Scanning all the GWAS variants would demonstrate a more informative profile of variance they explained. Thus, we analyzed all the 17 single nucleotide polymorphisms (SNPs) mapping to the 11 GWAS loci in an independent sample set of 800 Chinese subjects with PCOS and 1110 healthy controls systematically. Variants of rs3802457 in C9orf3 locus (P = 5.99×10⁻⁴) and rs13405728 in LHCGR locus (P = 3.73×10⁻⁴) were significantly associated with PCOS after the strict Bonferroni correction in our data set. The further haplotype analysis indicated that in the block of C9orf3 gene (rs4385527 and rs3802457), GA haplotype played a protective role in PCOS (8.7 vs 5.0, P = 9.85×10⁻⁶, OR = 0.548, 95%CI = 0.418–0.717), while GG haplotype was found suffering from an extraordinarily increased risk of PCOS (73.6% vs79.2%, P = 3.41×10⁻⁵, OR = 1.394, 95%CI = 1.191–1.632). Moreover, the directions of effects for all SNPs were consistent with previous GWAS reports (P = 1.53×10⁻⁵). Polygenic score analysis demonstrated that these 17 SNPs have a significant capacity on predicting case-control status in our samples (P = 7.17×10⁻⁹), meanwhile all these gathered 17 SNPs explained about 2.40% of variance. Our findings supported that C9orf3 and LHCGR loci variants were vital susceptibility of PCOS.     

Genetic variation at the FADS1-FADS2 gene locus influences delta-5 desaturase activity and LC-PUFA proportions after fish oil supplement

Delta-5 and delta-6 desaturases (D5D and D6D) are key enzymes in endogenous synthesis of long-chain PUFAs. In this sample of healthy subjects (n = 310), genotypes of single nucleotide polymorphisms (SNPs) rs174537, rs174561, and rs3834458 in the FADS1-FADS2 gene cluster were strongly associated with proportions of LC-PUFAs and desaturase activities estimated in plasma and erythrocytes. Minor allele carriage associated with decreased activities of D5D (FADS1) (5.84 × 10⁻¹⁹ ≤ P ≤ 4.5 × 10⁻¹⁸) and D6D (FADS2) (6.05 × 10⁻⁸ ≤ P ≤ 4.20 × 10⁻⁷) was accompanied by increased substrate and decreased product proportions (0.05 ≤ P ≤ 2.49 × 10⁻¹⁶). The significance of haplotype association with D5D activity (P = 2.19 × 10⁻¹⁷) was comparable to that of single SNPs, but haplotype association with D6D activity (P = 3.39 × 10⁻²⁸) was much stronger. In a randomized controlled dietary intervention, increasing eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3) intake significantly increased D5D (P = 4.0 × 10⁻⁹) and decreased D6D activity (P = 9.16 × 10⁻⁶) after doses of 0.45, 0.9, and 1.8 g/day for six months. Interaction of rs174537 genotype with treatment was a determinant of D5D activity estimated in plasma (P = 0.05). In conclusion, different sites at the FADS1-FADS2 locus appear to influence D5D and D6D activity, and rs174537 genotype interacts with dietary EPA+DHA to modulate D5D.     

Trans-Ancestral Studies Fine Map the SLE-Susceptibility Locus TNFSF4

We previously established an 80 kb haplotype upstream of TNFSF4 as a susceptibility locus in the autoimmune disease SLE. SLE-associated alleles at this locus are associated with inflammatory disorders, including atherosclerosis and ischaemic stroke. In Europeans, the TNFSF4 causal variants have remained elusive due to strong linkage disequilibrium exhibited by alleles spanning the region. Using a trans-ancestral approach to fine-map the locus, utilising 17,900 SLE and control subjects including Amerindian/Hispanics (1348 cases, 717 controls), African-Americans (AA) (1529, 2048) and better powered cohorts of Europeans and East Asians, we find strong association of risk alleles in all ethnicities; the AA association replicates in African-American Gullah (152,122). The best evidence of association comes from two adjacent markers: rs2205960-T (P = 1.71×10⁻³⁴, OR = 1.43[1.26–1.60]) and rs1234317-T (P = 1.16×10⁻²⁸, OR = 1.38[1.24–1.54]). Inference of fine-scale recombination rates for all populations tested finds the 80 kb risk and non-risk haplotypes in all except African-Americans. In this population the decay of recombination equates to an 11 kb risk haplotype, anchored in the 5′ region proximal to TNFSF4 and tagged by rs2205960-T after 1000 Genomes phase 1 (v3) imputation. Conditional regression analyses delineate the 5′ risk signal to rs2205960-T and the independent non-risk signal to rs1234314-C. Our case-only and SLE-control cohorts demonstrate robust association of rs2205960-T with autoantibody production. The rs2205960-T is predicted to form part of a decameric motif which binds NF-κBp65 with increased affinity compared to rs2205960-G. ChIP-seq data also indicate NF-κB interaction with the DNA sequence at this position in LCL cells. Our research suggests association of rs2205960-T with SLE across multiple groups and an independent non-risk signal at rs1234314-C. rs2205960-T is associated with autoantibody production and lymphopenia. Our data confirm a global signal at TNFSF4 and a role for the expressed product at multiple stages of lymphocyte dysregulation during SLE pathogenesis. We confirm the validity of trans-ancestral mapping in a complex trait.     

Targeting 160 Candidate Genes for Blood Pressure Regulation with a Genome-Wide Genotyping Array

The outcome of Genome-Wide Association Studies (GWAS) has challenged the field of blood pressure (BP) genetics as previous candidate genes have not been among the top loci in these scans. We used Affymetrix 500K genotyping data of KORA S3 cohort (n = 1,644; Southern-Germany) to address (i) SNP coverage in 160 BP candidate genes; (ii) the evidence for associations with BP traits in genome-wide and replication data, and haplotype analysis. In total, 160 gene regions (genic region±10 kb) covered 2,411 SNPs across 11.4 Mb. Marker densities in genes varied from 0 (n = 11) to 0.6 SNPs/kb. On average 52.5% of the HAPMAP SNPs per gene were captured. No evidence for association with BP was obtained for 1,449 tested SNPs. Considerable associations (P<10⁻³) were detected for the genes, where >50% of HAPMAP SNPs were tagged. In general, genes with higher marker density (>0.2 SNPs/kb) revealed a better chance to reach close to significance associations. Although, none of the detected P-values remained significant after Bonferroni correction (P<0.05/2319, P<2.15×10⁻⁵), the strength of some detected associations was close to this level: rs10889553 (LEPR) and systolic BP (SBP) (P = 4.5×10⁻⁵) as well as rs10954174 (LEP) and diastolic BP (DBP) (P = 5.20×10⁻⁵). In total, 12 markers in 7 genes (ADRA2A, LEP, LEPR, PTGER3, SLC2A1, SLC4A2, SLC8A1) revealed considerable association (P<10⁻³) either with SBP, DBP, and/or hypertension (HYP). None of these were confirmed in replication samples (KORA S4, HYPEST, BRIGHT). However, supportive evidence for the association of rs10889553 (LEPR) and rs11195419 (ADRA2A) with BP was obtained in meta-analysis across samples stratified either by body mass index, smoking or alcohol consumption. Haplotype analysis highlighted LEPR and PTGER3. In conclusion, the lack of associations in BP candidate genes may be attributed to inadequate marker coverage on the genome-wide arrays, small phenotypic effects of the loci and/or complex interaction with life-style and metabolic parameters.     

Fine Mapping on Chromosome 10q22-q23 Implicates Neuregulin 3 in Schizophrenia

Linkage studies have implicated 10q22-q23 as a schizophrenia (SZ) susceptibility locus in Ashkenazi Jewish (AJ) and Han Chinese from Taiwan populations. To further explore our previous linkage signal in the AJ population (NPL score: 4.27, empirical p = 2 × 10⁻⁵), we performed a peakwide association fine mapping study by using 1414 SNPs across ~12.5 Mb in 10q22-q23. We genotyped 1515 AJ individuals, including 285 parent-child trios, 173 unrelated cases, and 487 unrelated controls. We analyzed the binary diagnostic phenotype of SZ and 9 heritable quantitative traits derived from a principal components factor analysis of 73 items from our consensus diagnostic ratings and direct assessment interviews. Although no marker withstood multiple test correction for association with the binary SZ phenotype, we found strong evidence of association by using the “delusion” factor as the quantitative trait at three SNPs (rs10883866, rs10748842, and rs6584400) located in a 13 kb interval in intron 1 of Neuregulin 3 (NRG3). Our best p value from family-based association analysis was 7.26 × 10⁻⁷. We replicated this association in the collection of 173 unrelated AJ cases (p = 1.55 × 10⁻²), with a combined p value of 2.30 × 10⁻⁷. After performing 10,000 permutations of each of the phenotypes, we estimated the empirical study-wide significance across all 9 factors (90,000 permutations) to be p = 2.7 × 10⁻³. NRG3 is primarily expressed in the central nervous system and is one of three paralogs of NRG1, a gene strongly implicated in SZ. These biological properties together with our linkage and association results strongly support NRG3 as a gene involved in SZ.     

Genetic Analysis Identifies DDR2 as a Novel Gene Affecting Bone Mineral Density and Osteoporotic Fractures in Chinese Population

DDR2 gene, playing an essential role in regulating osteoblast differentiation and chondrocyte maturation, may influence bone mineral density (BMD) and osteoporosis, but the genetic variations actually leading to the association remain to be elucidated. Therefore, the aim of this study was to investigate whether the genetic variants in DDR2 are associated with BMD and fracture risk. This study was performed in three samples from two ethnicities, including 1,300 Chinese Han subjects, 700 Chinese Han subjects (350 with osteoporotic hip fractures and 350 healthy controls) and 2,286 US white subjects. Twenty-eight SNPs in DDR2 were genotyped and tested for associations with hip BMD and fractures. We identified 3 SNPs in DDR2 significantly associated with hip BMD in the Chinese population after multiple testing adjustments, which were rs7521233 (P = 1.06×10⁻⁴, β: −0.018 for allele C), rs7553831 (P = 1.30×10⁻⁴, β: −0.018 for allele T), and rs6697469 (P = 1.59×10⁻³, β: −0.015 for allele C), separately. These three SNPs were in high linkage disequilibrium. Haplotype analyses detected two significantly associated haplotypes, including one haplotype in block 2 (P = 9.54×10⁻⁴, β: −0.016) where these three SNPs located. SNP rs6697469 was also associated with hip fractures (P = 0.043, OR: 1.42) in the Chinese population. The effect on fracture risk was consistent with its association with lower BMD. However, in the white population, we didn’t observe significant associations with hip BMD. eQTL analyses revealed that SNPs associated with BMD also affected DDR2 mRNA expression levels in Chinese. Our findings, together with the prior biological evidence, suggest that DDR2 could be a new candidate for osteoporosis in Chinese population. Our results also reveal an ethnic difference, which highlights the need for further genetic studies in each ethnic group.     

Novel association of ABO histo-blood group antigen with soluble ICAM-1: results of a genome-wide association study of 6,578 women

While circulating levels of soluble Intercellular Adhesion Molecule 1 (sICAM-1) have been associated with diverse conditions including myocardial infarction, stroke, malaria, and diabetes, comprehensive analysis of the common genetic determinants of sICAM-1 is not available. In a genome-wide association study conducted among 6,578 participants in the Women's Genome Health Study, we find that three SNPs at the ICAM1 (19p13.2) locus (rs1799969, rs5498 and rs281437) are non-redundantly associated with plasma sICAM-1 concentrations at a genome-wide significance level (P<5×10⁻⁸), thus extending prior results from linkage and candidate gene studies. We also find that a single SNP (rs507666, P=5.1×10⁻²⁹) at the ABO (9q34.2) locus is highly correlated with sICAM-1 concentrations. The novel association at the ABO locus provides evidence for a previously unknown regulatory role of histo-blood group antigens in inflammatory adhesion processes.     

PRKCA: A Positional Candidate Gene for Body Mass Index and Asthma

Asthma incidence and prevalence are higher in obese individuals. A potential mechanistic basis for this relationship is pleiotropy. We hypothesized that significant linkage and candidate-gene association would be found for body mass index (BMI) in a population ascertained on asthma affection status. Linkage analysis for BMI was performed on 657 subjects in eight Costa Rican families enrolled in a study of asthma. Family-based association studies were conducted for BMI with SNPs within a positional candidate gene, PRKCA. SNPs within PRKCA were also tested for association with asthma. Association studies were conducted in 415 Costa Rican parent-child trios and 493 trios participating in the Childhood Asthma Management Program (CAMP). Although only modest evidence of linkage for BMI was obtained for the whole cohort, significant linkage was noted for BMI in females on chromosome 17q (peak LOD = 3.39). Four SNPs in a candidate gene in this region (PRKCA) had unadjusted association p values < 0.05 for BMI in both cohorts, with the joint p value for two SNPs remaining significant after adjustment for multiple comparisons (rs228883 and rs1005651, joint p values = 9.5 × 10⁻⁵ and 5.6 × 10⁻⁵). Similarly, eight SNPs had unadjusted association p values < 0.05 for asthma in both populations, with one SNP remaining significant after adjustment for multiple comparisons (rs11079657, joint p value = 2.6 × 10⁻⁵). PRKCA is a pleiotropic locus that is associated with both BMI and asthma and that has been identified via linkage analysis of BMI in a population ascertained on asthma.     

Genome-Wide Scan Identifies Variant in TNFSF13 Associated with Serum IgM in a Healthy Chinese Male Population

IgM provides a first line of defense during microbial infections. Serum IgM levels are detected routinely in clinical practice. And IgM is a genetically complex trait. We conducted a two-stage genome-wide association study (GWAS) to identify genetic variants affecting serum IgM levels in a Chinese population of 3495, including 1999 unrelated subjects in the first stage and 1496 independent individuals in the second stage. Our data show that a common single nucleotide polymorphism (SNP), rs11552708 located in the TNFSF13 gene was significantly associated with IgM levels (p = 5.00×10⁻⁷ in first stage, p = 1.34×10⁻³ in second stage, and p = 4.22×10⁻⁹ when combined). Besides, smoking was identified to be associated with IgM levels in both stages (P<0.05), but there was no significant interaction between smoking and the identified SNP (P>0.05). It is suggested that TNFSF13 may be a susceptibility gene affecting serum IgM levels in Chinese male population.     

Genome Wide Meta-analysis Highlights the Role of Genetic Variation in RARRES2 in the Regulation of Circulating Serum Chemerin

Chemerin is an adipokine proposed to link obesity and chronic inflammation of adipose tissue. Genetic factors determining chemerin release from adipose tissue are yet unknown. We conducted a meta-analysis of genome-wide association studies (GWAS) for serum chemerin in three independent cohorts from Europe: Sorbs and KORA from Germany and PPP-Botnia from Finland (total N = 2,791). In addition, we measured mRNA expression of genes within the associated loci in peripheral mononuclear cells by micro-arrays, and within adipose tissue by quantitative RT-PCR and performed mRNA expression quantitative trait and expression-chemerin association studies to functionally substantiate our loci. Heritability estimate of circulating chemerin levels was 16.2% in the Sorbs cohort. Thirty single nucleotide polymorphisms (SNPs) at chromosome 7 within the retinoic acid receptor responder 2 (RARRES2)/Leucine Rich Repeat Containing (LRRC61) locus reached genome-wide significance (p<5.0×10⁻⁸) in the meta-analysis (the strongest evidence for association at rs7806429 with p = 7.8×10⁻¹⁴, beta = −0.067, explained variance 2.0%). All other SNPs within the cluster were in linkage disequilibrium with rs7806429 (minimum r² = 0.43 in the Sorbs cohort). The results of the subgroup analyses of males and females were consistent with the results found in the total cohort. No significant SNP-sex interaction was observed. rs7806429 was associated with mRNA expression of RARRES2 in visceral adipose tissue in women (p<0.05 after adjusting for age and body mass index). In conclusion, the present meta-GWAS combined with mRNA expression studies highlights the role of genetic variation in the RARRES2 locus in the regulation of circulating chemerin concentrations.     

Association of Genetic Variants at 3q22 with Nephropathy in Patients with Type 1 Diabetes Mellitus

Diabetic nephropathy (DN) is the primary cause of morbidity and mortality in patients with type 1 diabetes mellitus (T1DM) and affects about 30% of these patients. We have previously localized a DN locus on chromosome 3q with suggestive linkage in Finnish individuals. Linkage to this region has also been reported earlier by several other groups. To fine map this locus, we conducted a multistage case-control association study in T1DM patients, comprising 1822 cases with nephropathy and 1874 T1DM patients free of nephropathy, from Finland, Iceland, and the British Isles. At the screening stage, we genotyped 3072 tag SNPs, spanning a 28 Mb region, in 234 patients and 215 controls from Finland. SNPs that met the significance threshold of p < 0.01 at this stage were followed up by a series of sample sets. A genetic variant, rs1866813, in the noncoding region at 3q22 was associated with increased risk of DN (overall p = 7.07 × 10⁻⁶, combined odds ratio [OR] of the allele = 1.33). The estimated genotypic ORs of this variant in all Finnish samples suggested a codominant effect, resulting in significant association, with a p value of 4.7 × 10⁻⁵ (OR = 1.38; 95% confidence interval = 1.18–1.62). Additionally, an 11 kb segment flanked by rs62408925 and rs1866813, two strongly correlated variants (r² = 0.95), contains three elements highly conserved across multiple species. Independent replication will clarify the role of the associated variants at 3q22 in influencing the risk of DN.     

Fine Mapping and Functional Analysis of the Multiple Sclerosis Risk Gene CD6

CD6 has recently been identified and validated as risk gene for multiple sclerosis (MS), based on the association of a single nucleotide polymorphism (SNP), rs17824933, located in intron 1. CD6 is a cell surface scavenger receptor involved in T-cell activation and proliferation, as well as in thymocyte differentiation. In this study, we performed a haptag SNP screen of the CD6 gene locus using a total of thirteen tagging SNPs, of which three were non-synonymous SNPs, and replicated the recently reported GWAS SNP rs650258 in a Spanish-Basque collection of 814 controls and 823 cases. Validation of the six most strongly associated SNPs was performed in an independent collection of 2265 MS patients and 2600 healthy controls. We identified association of haplotypes composed of two non-synonymous SNPs [rs11230563 (R225W) and rs2074225 (A257V)] in the 2^nd SRCR domain with susceptibility to MS (P _{max(T) permutation} = 1×10⁻⁴). The effect of these haplotypes on CD6 surface expression and cytokine secretion was also tested. The analysis showed significantly different CD6 expression patterns in the distinct cell subsets, i.e. – CD4⁺ naïve cells, P = 0.0001; CD8⁺ naïve cells, P<0.0001; CD4⁺ and CD8⁺ central memory cells, P = 0.01 and 0.05, respectively; and natural killer T (NKT) cells, P = 0.02; with the protective haplotype (RA) showing higher expression of CD6. However, no significant changes were observed in natural killer (NK) cells, effector memory and terminally differentiated effector memory T cells. Our findings reveal that this new MS-associated CD6 risk haplotype significantly modifies expression of CD6 on CD4⁺ and CD8⁺ T cells.     

Variants in TF and HFE Explain ∼40% of Genetic Variation in Serum-Transferrin Levels

Only a small proportion of genetic variation in complex traits has been explained by SNPs from genome-wide association studies (GWASs). We report the results from two GWASs for serum markers of iron status (serum iron, serum transferrin, transferrin saturation with iron, and serum ferritin), which are important in iron overload (e.g., hemochromatosis) and deficiency (e.g., anemia) conditions. We performed two GWASs on samples of Australians of European descent. In the first GWAS, 411 adolescent twins and their siblings were genotyped with 100K SNPs. rs1830084, 10.8 kb 3′ of TF, was significantly associated with serum transferrin (p total association test = 1.0 × 10⁻⁹; p within-family test = 2.2 × 10⁻⁵). In the second GWAS on an independent sample of 459 female monozygotic (MZ) twin pairs genotyped with 300K SNPs, we found rs3811647 (within intron 11 of TF, HapMap CEU r² with rs1830084 = 0.86) was significantly associated with serum transferrin (p = 3.0 × 10⁻¹⁵). In the second GWAS, we found two additional and independent SNPs on TF (rs1799852 and rs2280673) and confirmed the known C282Y mutation in HFE to be independently associated with serum transferrin. The three variants in TF (rs3811647, rs1799852 and rs2280673) plus the HFE C282Y mutation explained ~40% of genetic variation in serum transferrin (p = 7.8 × 10⁻²⁵). These findings are potentially important for our understanding of iron metabolism and of regulation of hepatic protein secretion, and also strongly support the hypothesis that the genetic architecture of some endophenotypes may be simpler than that of disease.     

Genome-Wide Association Study of Late-Onset Myasthenia Gravis: Confirmation of TNFRSF11A and Identification of ZBTB10 and Three Distinct HLA Associations

To investigate the genetics of late-onset myasthenia gravis (LOMG), we conducted a genome-wide association study imputation of >6 million single nucleotide polymorphisms (SNPs) in 532 LOMG cases (anti–acetylcholine receptor [AChR] antibody positive; onset age ≥50 years) and 2,128 controls matched for sex and population substructure. The data confirm reported TNFRSF11A associations (rs4574025, P = 3.9 × 10⁻⁷, odds ratio [OR] 1.42) and identify a novel candidate gene, ZBTB10, achieving genome-wide significance (rs6998967, P = 8.9 × 10⁻¹⁰, OR 0.53). Several other SNPs showed suggestive significance including rs2476601 (P = 6.5 × 10⁻⁶, OR 1.62) encoding the PTPN22 R620W variant noted in early-onset myasthenia gravis (EOMG) and other autoimmune diseases. In contrast, EOMG-associated SNPs in TNIP1 showed no association in LOMG, nor did other loci suggested for EOMG. Many SNPs within the major histocompatibility complex (MHC) region showed strong associations in LOMG, but with smaller effect sizes than in EOMG (highest OR ~2 versus ~6 in EOMG). Moreover, the strongest associations were in opposite directions from EOMG, including an OR of 0.54 for DQA1*05:01 in LOMG (P = 5.9 × 10⁻¹²) versus 2.82 in EOMG (P = 3.86 × 10⁻⁴⁵). Association and conditioning studies for the MHC region showed three distinct and largely independent association peaks for LOMG corresponding to (a) MHC class II (highest attenuation when conditioning on DQA1), (b) HLA-A and (c) MHC class III SNPs. Conditioning studies of human leukocyte antigen (HLA) amino acid residues also suggest potential functional correlates. Together, these findings emphasize the value of subgrouping myasthenia gravis patients for clinical and basic investigations and imply distinct predisposing mechanisms in LOMG.     

A at Single Nucleotide Polymorphism-358 Is Required for G at -420 to Confer the Highest Plasma Resistin in the General Japanese Population

Insulin resistance is a feature of type 2 diabetes. Resistin, secreted from adipocytes, causes insulin resistance in mice. We previously reported that the G/G genotype of single nucleotide polymorphism (SNP) at −420 (rs1862513) in the human resistin gene (RETN) increased susceptibility to type 2 diabetes by enhancing its promoter activity. Plasma resistin was highest in Japanese subjects with G/G genotype, followed by C/G, and C/C. In this study, we cross-sectionally analyzed plasma resistin and SNPs in the RETN region in 2,019 community-dwelling Japanese subjects. Plasma resistin was associated with SNP-638 (rs34861192), SNP-537 (rs34124816), SNP-420, SNP-358 (rs3219175), SNP+299 (rs3745367), and SNP+1263 (rs3745369) (P<10⁻¹³ in all cases). SNP-638, SNP -420, SNP-358, and SNP+157 were in the same linkage disequilibrium (LD) block. SNP-358 and SNP-638 were nearly in complete LD (r² = 0.98), and were tightly correlated with SNP-420 (r² = 0.50, and 0.51, respectively). The correlation between either SNP-358 (or SNP-638) or SNP-420 and plasma resistin appeared to be strong (risk alleles for high plasma resistin; A at SNP-358, r² = 0.5224, P = 4.94×10⁻³²⁴; G at SNP-420, r² = 0.2616, P = 1.71×10⁻¹³³). In haplotypes determined by SNP-420 and SNP-358, the estimated frequencies for C-G, G-A, and G-G were 0.6700, 0.2005, and 0.1284, respectively, and C-A was rare (0.0011), suggesting that subjects with A at −358, generally had G at −420. This G-A haplotype conferred the highest plasma resistin (8.24 ng/ml difference/allele compared to C-G, P<0.0001). In THP-1 cells, the RETN promoter with the G-A haplotype showed the highest activity. Nuclear proteins specifically recognized one base difference at SNP-358, but not at SNP-638. Therefore, A at -358 is required for G at −420 to confer the highest plasma resistin in the general Japanese population. In Caucasians, the association between SNP-420 and plasma resistin is not strong, and A at −358 may not exist, suggesting that SNP-358 could explain this ethnic difference.     

Influence of Genetic Variation on Plasma Protein Levels in Older Adults Using a Multi-Analyte Panel

Proteins, widely studied as potential biomarkers, play important roles in numerous physiological functions and diseases. Genetic variation may modulate corresponding protein levels and point to the role of these variants in disease pathophysiology. Effects of individual single nucleotide polymorphisms (SNPs) within a gene were analyzed for corresponding plasma protein levels using genome-wide association study (GWAS) genotype data and proteomic panel data with 132 quality-controlled analytes from 521 Caucasian participants in the Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort. Linear regression analysis detected 112 significant (Bonferroni threshold p = 2.44×10⁻⁵) associations between 27 analytes and 112 SNPs. 107 out of these 112 associations were tested in the Indiana Memory and Aging Study (IMAS) cohort for replication and 50 associations were replicated at uncorrected p<0.05 in the same direction of effect as those in the ADNI. We identified multiple novel associations including the association of rs7517126 with plasma complement factor H-related protein 1 (CFHR1) level at p<1.46×10⁻⁶⁰, accounting for 40 percent of total variation of the protein level. We serendipitously found the association of rs6677604 with the same protein at p<9.29×10⁻¹¹². Although these two SNPs were not in the strong linkage disequilibrium, 61 percent of total variation of CFHR1 was accounted for by rs6677604 without additional variation by rs7517126 when both SNPs were tested together. 78 other SNP-protein associations in the ADNI sample exceeded genome-wide significance (5×10⁻⁸). Our results confirmed previously identified gene-protein associations for interleukin-6 receptor, chemokine CC-4, angiotensin-converting enzyme, and angiotensinogen, although the direction of effect was reversed in some cases. This study is among the first analyses of gene-protein product relationships integrating multiplex-panel proteomics and targeted genes extracted from a GWAS array. With intensive searches taking place for proteomic biomarkers for many diseases, the role of genetic variation takes on new importance and should be considered in interpretation of proteomic results.     

A Genome-Wide Association Study of Circulating Galectin-3

Galectin-3 is a lectin involved in fibrosis, inflammation and proliferation. Increased circulating levels of galectin-3 have been associated with various diseases, including cancer, immunological disorders, and cardiovascular disease. To enhance our knowledge on galectin-3 biology we performed the first genome-wide association study (GWAS) using the Illumina HumanCytoSNP-12 array imputed with the HapMap 2 CEU panel on plasma galectin-3 levels in 3,776 subjects and follow-up genotyping in an additional 3,516 subjects. We identified 2 genome wide significant loci associated with plasma galectin-3 levels. One locus harbours the LGALS3 gene (rs2274273; P = 2.35×10⁻¹⁸⁸) and the other locus the ABO gene (rs644234; P = 3.65×10⁻⁴⁷). The variance explained by the LGALS3 locus was 25.6% and by the ABO locus 3.8% and jointly they explained 29.2%. Rs2274273 lies in high linkage disequilibrium with two non-synonymous SNPs (rs4644; r² = 1.0, and rs4652; r² = 0.91) and wet lab follow-up genotyping revealed that both are strongly associated with galectin-3 levels (rs4644; P = 4.97×10⁻⁴⁶⁵ and rs4652 P = 1.50×10⁻⁴²¹) and were also associated with LGALS3 gene-expression. The origins of our associations should be further validated by means of functional experiments.     

Fine Mapping of a Region of Chromosome 11q23.3 Reveals Independent Locus Associated with Risk of Glioma

Background

A single nucleotide polymorphism (SNP) at locus 11q23.3 (rs498872) in the near 5′-UTR of the PHLDB1 gene was recently implicated as a risk factor for gliomas in a genome-wide association study, and this involvement was confirmed in three additional studies.

Methodology/Principal Findings

To identify possible causal variants in the region, the authors genotyped 15 tagging SNPs in the 200 kb genomic region at 11q23.3 locus in a Chinese Han population-based case-control study with 983 cases and 1024 controls. We found evidence for an association between two independent loci (both the PHLDB1 and the ACRN1 genes) and a predisposition for gliomas. Among the multiple significant SNPs in the PHLDB1 gene region, the rs17749 SNP was the most significant [P = 1.31×10⁻⁶ in a recessive genetic model]. Additionally, two novel SNPs (rs2236661 and rs494560) that were independent of rs17749 were significantly associated with glioma risk in a recessive genetic model [P = 1.31×10⁻⁵ and P = 3.32×10⁻⁵, respectively]. The second novel locus was within the ARCN1 gene, and it was associated with a significantly reduced risk for glioma.

Conclusions/Significance

Our data strongly support PHLDB1 as a susceptibility gene for glioma, also shedding light on a new potentially candidate gene, ARCN1.