Identification of novel susceptibility loci for non‐syndromic cleft lip with or without cleft palate

Although several genome‐wide association studies (GWAS) of non‐syndromic cleft lip with or without cleft palate (NSCL/P) have been reported, more novel association signals are remained to be exploited. Here, we performed an in‐depth analysis of our previously published Chinese GWAS cohort study with replication in an extra dbGaP case‐parent trios and another in‐house Nanjing cohort, and finally identified five novel significant association signals (rs11119445: 3’ of SERTAD4, P = 6.44 × 10⁻¹⁴; rs227227 and rs12561877: intron of SYT14, P = 5.02 × 10⁻¹³ and 2.80 × 10⁻¹¹, respectively; rs643118: intron of TRAF3IP3, P = 4.45 × 10⁻⁶; rs2095293: intron of NR6A1, P = 2.98 × 10⁻⁵). The mean (standard deviation) of the weighted genetic risk score (wGRS) from these SNPs was 1.83 (0.65) for NSCL/P cases and 1.58 (0.68) for controls, respectively (P = 2.67 × 10⁻¹⁶). Rs643118 was identified as a shared susceptible factor of NSCL/P among Asians and Europeans, while rs227227 may contribute to the risk of NSCL/P as well as NSCPO. In addition, sertad4 knockdown zebrafish models resulted in down‐regulation of sox2 and caused oedema around the heart and mandibular deficiency, compared with control embryos. Taken together, this study has improved our understanding of the genetic susceptibility to NSCL/P and provided further clues to its aetiology in the Chinese population.     

A Large-Scale Genetic Analysis Reveals a Strong Contribution of the HLA Class II Region to Giant Cell Arteritis Susceptibility

We conducted a large-scale genetic analysis on giant cell arteritis (GCA), a polygenic immune-mediated vasculitis. A case-control cohort, comprising 1,651 case subjects with GCA and 15,306 unrelated control subjects from six different countries of European ancestry, was genotyped by the Immunochip array. We also imputed HLA data with a previously validated imputation method to perform a more comprehensive analysis of this genomic region. The strongest association signals were observed in the HLA region, with rs477515 representing the highest peak (p = 4.05 × 10⁻⁴⁰, OR = 1.73). A multivariate model including class II amino acids of HLA-DRβ1 and HLA-DQα1 and one class I amino acid of HLA-B explained most of the HLA association with GCA, consistent with previously reported associations of classical HLA alleles like HLA-DRB1^∗04. An omnibus test on polymorphic amino acid positions highlighted DRβ1 13 (p = 4.08 × 10⁻⁴³) and HLA-DQα1 47 (p = 4.02 × 10⁻⁴⁶), 56, and 76 (both p = 1.84 × 10⁻⁴⁵) as relevant positions for disease susceptibility. Outside the HLA region, the most significant loci included PTPN22 (rs2476601, p = 1.73 × 10⁻⁶, OR = 1.38), LRRC32 (rs10160518, p = 4.39 × 10⁻⁶, OR = 1.20), and REL (rs115674477, p = 1.10 × 10⁻⁵, OR = 1.63). Our study provides evidence of a strong contribution of HLA class I and II molecules to susceptibility to GCA. In the non-HLA region, we confirmed a key role for the functional PTPN22 rs2476601 variant and proposed other putative risk loci for GCA involved in Th1, Th17, and Treg cell function.     

Genome-wide Association Study Identifies TNFSF15 and POU2AF1 as Susceptibility Loci for Primary Biliary Cirrhosis in the Japanese Population

For the identification of susceptibility loci for primary biliary cirrhosis (PBC), a genome-wide association study (GWAS) was performed in 963 Japanese individuals (487 PBC cases and 476 healthy controls) and in a subsequent replication study that included 1,402 other Japanese individuals (787 cases and 615 controls). In addition to the most significant susceptibility region, human leukocyte antigen (HLA), we identified two significant susceptibility loci, TNFSF15 (rs4979462) and POU2AF1 (rs4938534) (combined odds ratio [OR] = 1.56, p = 2.84 × 10⁻¹⁴ for rs4979462, and combined OR = 1.39, p = 2.38 × 10⁻⁸ for rs4938534). Among 21 non-HLA susceptibility loci for PBC identified in GWASs of individuals of European descent, three loci (IL7R, IKZF3, and CD80) showed significant associations (combined p = 3.66 × 10⁻⁸, 3.66 × 10⁻⁹, and 3.04 × 10⁻⁹, respectively) and STAT4 and NFKB1 loci showed suggestive association with PBC (combined p = 1.11 × 10⁻⁶ and 1.42 × 10⁻⁷, respectively) in the Japanese population. These observations indicated the existence of ethnic differences in genetic susceptibility loci to PBC and the importance of TNF signaling and B cell differentiation for the development of PBC in individuals of European descent and Japanese individuals.     

Gender-specific genetic associations of polymorphisms in ACE,AKR1C2, FTO and MMP2 with weight gain over a 10-year period

Weight gain, when it leads to overweight or obesity, is nowadays one of the major health problems. ACE, FTO, AKR1C2, TIMP4 and MMP2 genes have been implicated in previous studies on weight regulation. This study investigated the contribution of polymorphisms in these five candidate genes to the risk of weight gain over a 10-year time period. Two groups were selected from participants of the Doetinchem cohort study who were followed over a 10-year period: A stable weight group (±2 kg/10 year; n = 259) and a weight gainers group who increased their body weight by roughly 10 % (>8 kg/10 year; n = 237). Starting BMI was between 20 and 35 kg/m² and baseline age between 20 and 45 years. Selected SNPs and insert/deletion in candidate genes were measured in each group. In men, the allelic distribution of FTO rs9939609 (χ²p = 0.005), ACE rs4340 (χ²p = 0.006) and AKR1C2 rs12249281 (χ²p = 0.019) differed between the weight stable and weight gainers group. Interaction between FTO rs9939609 and ACE rs4340 was observed. In women, the allelic distribution of MMP2 rs1132896 differed between the weight stable and weight gainers group (χ²p = 0.00001). The A-allele of FTO was associated with a 1.99× higher risk of gaining weight in men (OR 1.99, p = 0.020), while in women, the C-allele of MMP2 was associated with a 2.50× higher risk of weight gain (OR 2.50, p = 0.001) over the 10-year period. We found that FTO in men and MMP2 in women are associated with weight gain over a 10-year follow-up period.

Electronic supplementary material

The online version of this article (doi:10.1007/s12263-014-0434-2) contains supplementary material, which is available to authorized users.     

Genome-wide Association Study Identifies Five Susceptibility Loci for Follicular Lymphoma outside the HLA Region

Genome-wide association studies (GWASs) of follicular lymphoma (FL) have previously identified human leukocyte antigen (HLA) gene variants. To identify additional FL susceptibility loci, we conducted a large-scale two-stage GWAS in 4,523 case subjects and 13,344 control subjects of European ancestry. Five non-HLA loci were associated with FL risk: 11q23.3 (rs4938573, p = 5.79 × 10⁻²⁰) near CXCR5; 11q24.3 (rs4937362, p = 6.76 × 10⁻¹¹) near ETS1; 3q28 (rs6444305, p = 1.10 × 10⁻¹⁰) in LPP; 18q21.33 (rs17749561, p = 8.28 × 10⁻¹⁰) near BCL2; and 8q24.21 (rs13254990, p = 1.06 × 10⁻⁸) near PVT1. In an analysis of the HLA region, we identified four linked HLA-DRβ1 multiallelic amino acids at positions 11, 13, 28, and 30 that were associated with FL risk (p_omnibus = 4.20 × 10⁻⁶⁷ to 2.67 × 10⁻⁷⁰). Additional independent signals included rs17203612 in HLA class II (odds ratio [OR_per-allele] = 1.44; p = 4.59 × 10⁻¹⁶) and rs3130437 in HLA class I (OR_per-allele = 1.23; p = 8.23 × 10⁻⁹). Our findings further expand the number of loci associated with FL and provide evidence that multiple common variants outside the HLA region make a significant contribution to FL risk.     

Identification of Multiple Genetic Susceptibility Loci in Takayasu Arteritis

Takayasu arteritis is a rare inflammatory disease of large arteries. The etiology of Takayasu arteritis remains poorly understood, but genetic contribution to the disease pathogenesis is supported by the genetic association with HLA-B^∗52. We genotyped ∼200,000 genetic variants in two ethnically divergent Takayasu arteritis cohorts from Turkey and North America by using a custom-designed genotyping platform (Immunochip). Additional genetic variants and the classical HLA alleles were imputed and analyzed. We identified and confirmed two independent susceptibility loci within the HLA region (r² < 0.2): HLA-B/MICA (rs12524487, OR = 3.29, p = 5.57 × 10⁻¹⁶) and HLA-DQB1/HLA-DRB1 (rs113452171, OR = 2.34, p = 3.74 × 10⁻⁹; and rs189754752, OR = 2.47, p = 4.22 × 10⁻⁹). In addition, we identified and confirmed a genetic association between Takayasu arteritis and the FCGR2A/FCGR3A locus on chromosome 1 (rs10919543, OR = 1.81, p = 5.89 × 10⁻¹²). The risk allele in this locus results in increased mRNA expression of FCGR2A. We also established the genetic association between IL12B and Takayasu arteritis (rs56167332, OR = 1.54, p = 2.18 × 10⁻⁸).     

Common Polymorphisms Influencing Serum Uric Acid Levels Contribute to Susceptibility to Gout,but Not to Coronary Artery Disease

Background

Recently, a large meta-analysis including over 28,000 participants identified nine different loci with association to serum uric acid (UA) levels. Since elevated serum UA levels potentially cause gout and are a possible risk factor for coronary artery disease (CAD) and myocardial infarction (MI), we performed two large case-control association analyses with participants from the German MI Family Study. In the first study, we assessed the association of the qualitative trait gout and ten single nucleotide polymorphisms (SNP) markers that showed association to UA serum levels. In the second study, the same genetic polymorphisms were analyzed for association with CAD.

Methods and Findings

A total of 683 patients suffering from gout and 1,563 healthy controls from the German MI Family Study were genotyped. Nine SNPs were identified from a recently performed genome-wide meta-analysis on serum UA levels (rs12129861, rs780094, rs734553, rs2231142, rs742132, rs1183201, rs12356193, rs17300741 and rs505802). Additionally, the marker rs6855911 was included which has been associated with gout in our cohort in a previous study. SNPs rs734553 and rs6855911, located in SLC2A9, and SNP rs2231142, known to be a missense polymorphism in ABCG2, were associated with gout (p = 5.6*10⁻⁷, p = 1.1*10⁻⁷, and p = 1.3*10⁻³, respectively). Other SNPs in the genes PDZK1, GCKR, LRRC16A, SLC17A1-SLC17A3, SLC16A9, SLC22A11 and SLC22A12 failed the significance level. None of the ten markers were associated with risk to CAD in our study sample of 1,473 CAD cases and 1,241 CAD-free controls.

Conclusion

SNP markers in SLC2A9 and ABCG2 genes were found to be strongly associated with the phenotype gout. However, not all SNP markers influencing serum UA levels were also directly associated with the clinical manifestation of gout in our study sample. In addition, none of these SNPs showed association with the risk to CAD in the German MI Family Study.     

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.     

Genetic Variants Associated with Lipid Profiles in Chinese Patients with Type 2 Diabetes

Dyslipidemia is a strong risk factor for cardiovascular disease among patients with type 2 diabetes (T2D). The aim of this study was to identify lipid-related genetic variants in T2D patients of Han Chinese ancestry. Among 4,908 Chinese T2D patients who were not taking lipid-lowering medications, single nucleotide polymorphisms (SNPs) in seven genes previously found to be associated with lipid traits in genome-wide association studies conducted in populations of European ancestry (ABCA1, GCKR, BAZ1B, TOMM40, DOCK7, HNF1A, and HNF4A) were genotyped. After adjusting for multiple covariates, SNPs in ABCA1, GCKR, BAZ1B, TOMM40, and HNF1A were identified as significantly associated with triglyceride levels in T2D patients (P < 0.05). The associations between the SNPs in ABCA1 (rs3890182), GCKR (rs780094), and BAZ1B (rs2240466) remained significant even after correction for multiple testing (P = 8.85×10⁻³, 7.88×10⁻⁷, and 2.03×10⁻⁶, respectively). BAZ1B (rs2240466) also was associated with the total cholesterol level (P = 4.75×10⁻²). In addition, SNP rs157580 in TOMM40 was associated with the low-density lipoprotein cholesterol level (P = 6.94×10⁻³). Our findings confirm that lipid-related genetic loci are associated with lipid profiles in Chinese patients with type 2 diabetes.     

Susceptibility Loci Associated with Specific and Shared Subtypes of Lymphoid Malignancies

The genetics of lymphoma susceptibility reflect the marked heterogeneity of diseases that comprise this broad phenotype. However, multiple subtypes of lymphoma are observed in some families, suggesting shared pathways of genetic predisposition to these pathologically distinct entities. Using a two-stage GWAS, we tested 530,583 SNPs in 944 cases of lymphoma, including 282 familial cases, and 4,044 public shared controls, followed by genotyping of 50 SNPs in 1,245 cases and 2,596 controls. A novel region on 11q12.1 showed association with combined lymphoma (LYM) subtypes. SNPs in this region included rs12289961 near LPXN, (P_LYM = 3.89×10⁻⁸, OR = 1.29) and rs948562 (P_LYM = 5.85×10⁻⁷, OR = 1.29). A SNP in a novel non-HLA region on 6p23 (rs707824, P_NHL = 5.72×10⁻⁷) was suggestive of an association conferring susceptibility to lymphoma. Four SNPs, all in a previously reported HLA region, 6p21.32, showed genome-wide significant associations with follicular lymphoma. The most significant association with follicular lymphoma was for rs4530903 (P_FL = 2.69×10⁻¹², OR = 1.93). Three novel SNPs near the HLA locus, rs9268853, rs2647046, and rs2621416, demonstrated additional variation contributing toward genetic susceptibility to FL associated with this region. Genes implicated by GWAS were also found to be cis-eQTLs in lymphoblastoid cell lines; candidate genes in these regions have been implicated in hematopoiesis and immune function. These results, showing novel susceptibility regions and allelic heterogeneity, point to the existence of pathways of susceptibility to both shared as well as specific subtypes of lymphoid malignancy.     

A Genome-Wide Association Study Suggests Novel Loci Associated with a Schizophrenia-Related Brain-Based Phenotype

Patients with schizophrenia and their siblings typically show subtle changes of brain structures, such as a reduction of hippocampal volume. Hippocampal volume is heritable, may explain a variety of cognitive symptoms of schizophrenia and is thus considered an intermediate phenotype for this mental illness. The aim of our analyses was to identify single-nucleotide polymorphisms (SNP) related to hippocampal volume without making prior assumptions about possible candidate genes. In this study, we combined genetics, imaging and neuropsychological data obtained from the Mind Clinical Imaging Consortium study of schizophrenia (n = 328). A total of 743,591 SNPs were tested for association with hippocampal volume in a genome-wide association study. Gene expression profiles of human hippocampal tissue were investigated for gene regions of significantly associated SNPs. None of the genetic markers reached genome-wide significance. However, six highly correlated SNPs (rs4808611, rs35686037, rs12982178, rs1042178, rs10406920, rs8170) on chromosome 19p13.11, located within or in close proximity to the genes NR2F6, USHBP1, and BABAM1, as well as four SNPs in three other genomic regions (chromosome 1, 2 and 10) had p-values between 6.75×10⁻⁶ and 8.3×10⁻⁷. Using existing data of a very recently published GWAS of hippocampal volume and additional data of a multicentre study in a large cohort of adolescents of European ancestry, we found supporting evidence for our results. Furthermore, allelic differences in rs4808611 and rs8170 were highly associated with differential mRNA expression in the cis-acting region. Associations with memory functioning indicate a possible functional importance of the identified risk variants. Our findings provide new insights into the genetic architecture of a brain structure closely linked to schizophrenia. In silico replication, mRNA expression and cognitive data provide additional support for the relevance of our findings. Identification of causal variants and their functional effects may unveil yet unknown players in the neurodevelopment and the pathogenesis of neuropsychiatric disorders.     

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.     

Two Susceptibility Loci to Takayasu Arteritis Reveal a Synergistic Role of the IL12B and HLA-B Regions in a Japanese Population

Takayasu arteritis (TAK) is an autoimmune systemic vasculitis of unknown etiology. Although previous studies have revealed that HLA-B^∗52:01 has an effect on TAK susceptibility, no other genetic determinants have been established so far. Here, we performed genome scanning of 167 TAK cases and 663 healthy controls via Illumina Infinium Human Exome BeadChip arrays, followed by a replication study consisting of 212 TAK cases and 1,322 controls. As a result, we found that the IL12B region on chromosome 5 (rs6871626, overall p = 1.7 × 10⁻¹³, OR = 1.75, 95% CI 1.42–2.16) and the MLX region on chromosome 17 (rs665268, overall p = 5.2 × 10⁻⁷, OR = 1.50, 95% CI 1.28–1.76) as well as the HLA-B region (rs9263739, a proxy of HLA-B^∗52:01, overall p = 2.8 × 10⁻²¹, OR = 2.44, 95% CI 2.03–2.93) exhibited significant associations. A significant synergistic effect of rs6871626 and rs9263739 was found with a relative excess risk of 3.45, attributable proportion of 0.58, and synergy index of 3.24 (p ≤ 0.00028) in addition to a suggestive synergistic effect between rs665268 and rs926379 (p ≤ 0.027). We also found that rs6871626 showed a significant association with clinical manifestations of TAK, including increased risk and severity of aortic regurgitation, a representative severe complication of TAK. Detection of these susceptibility loci will provide new insights to the basic mechanisms of TAK pathogenesis. Our findings indicate that IL12B plays a fundamental role on the pathophysiology of TAK in combination with HLA-B^∗52:01 and that common autoimmune mechanisms underlie the pathology of TAK and other autoimmune disorders such as psoriasis and inflammatory bowel diseases in which IL12B is involved as a genetic predisposing factor.     

Inherited genetic variants associated with glucocorticoid sensitivity in leukaemia cells

Identification of genetic variants associated with glucocorticoids (GC) sensitivity of leukaemia cells may provide insight into potential drug targets and tailored therapy. In the present study, within 72 leukaemic cell lines derived from Japanese patients with B‐cell precursor acute lymphoblastic leukaemia (ALL), we conducted genome‐wide genotyping of single nucleotide polymorphisms (SNP) and attempted to identify genetic variants associated with GC sensitivity and NR3C1 (GC receptor) gene expression. IC50 measures for prednisolone (Pred) and dexamethasone (Dex) were available using an alamarBlue cell viability assay. IC50 values of Pred showed the strongest association with rs904419 (P = 4.34 × 10⁻⁸), located between the FRMD4B and MITF genes. The median IC50 values of prednisolone for cell lines with rs904419 AA (n = 13), AG (n = 31) and GG (n = 28) genotypes were 0.089, 0.139 and 297 µmol/L, respectively. For dexamethasone sensitivity, suggestive association was observed for SNP rs2306888 (P = 1.43 × 10⁻⁶), a synonymous SNP of the TGFBR3 gene. For NR3C1 gene expression, suggestive association was observed for SNP rs11982167 (P = 6.44 × 10⁻⁸), located in the PLEKHA8 gene. These genetic variants may affect GC sensitivity of ALL cells and may give rise to opportunities in personalized medicine for effective and safe chemotherapy in ALL patients.     

Trans-ancestral fine-mapping of MHC reveals key amino acids associated with spontaneous clearance of hepatitis C in HLA-DQβ1

Spontaneous clearance of acute hepatitis C virus (HCV) infection is associated with single nucleotide polymorphisms (SNPs) on the MHC class II. We fine-mapped the MHC region in European (n = 1,600; 594 HCV clearance/1,006 HCV persistence) and African (n = 1,869; 340 HCV clearance/1,529 HCV persistence) ancestry individuals and evaluated HCV peptide binding affinity of classical alleles. In both populations, HLA-DQβ1Leu26 (p value_Meta = 1.24 × 10⁻¹⁴) located in pocket 4 was negatively associated with HCV spontaneous clearance and HLA-DQβ1Pro55 (p value_Meta = 8.23 × 10⁻¹¹) located in the peptide binding region was positively associated, independently of HLA-DQβ1Leu26. These two amino acids are not in linkage disequilibrium (r² < 0.1) and explain the SNPs and classical allele associations represented by rs2647011, rs9274711, HLA-DQB1^∗03:01, and HLA-DRB1^∗01:01. Additionally, HCV persistence classical alleles tagged by HLA-DQβ1Leu26 had fewer HCV binding epitopes and lower predicted binding affinities compared to clearance alleles (geometric mean of combined IC₅₀ nM of persistence versus clearance; 2,321 nM versus 761.7 nM, p value = 1.35 × 10⁻³⁸). In summary, MHC class II fine-mapping revealed key amino acids in HLA-DQβ1 explaining allelic and SNP associations with HCV outcomes. This mechanistic advance in understanding of natural recovery and immunogenetics of HCV might set the stage for much needed enhancement and design of vaccine to promote spontaneous clearance of HCV infection.     

Association of Parkinson Disease with Structural and Regulatory Variants in the HLA Region

Historically, association of disease with the major histocompatibility complex (HLA) genes has been tested with HLA alleles that encode antigen-binding affinity. The association with Parkinson disease (PD), however, was discovered with noncoding SNPs in a genome-wide association study (GWAS). We show here that several HLA-region SNPs that have since been associated with PD form two blocks tagged by rs3129882 (p = 9 × 10⁻¹¹) and by rs9268515 and/or rs2395163 (p = 3 × 10⁻¹¹). We investigated whether these SNP-associations were driven by HLA-alleles at adjacent loci. We imputed class I and class II HLA-alleles for 2000 PD cases and 1986 controls from the NeuroGenetics Research Consortium GWAS and sequenced a subset of 194 cases and 204 controls. We were therefore able to assess accuracy of two imputation algorithms against next-generation-sequencing while taking advantage of the larger imputed data sets for disease study. Additionally, we imputed HLA alleles for 843 cases and 856 controls from another GWAS for replication. PD risk was positively associated with the B^∗07:02_C^∗07:02_DRB5^∗01_DRB1^∗15:01_DQA1^∗01:02_DQB1^∗06:02 haplotype and negatively associated with the C^∗03:04, DRB1^∗04:04 and DQA1^∗03:01 alleles. The risk haplotype and DQA1^∗03:01 lost significance when conditioned on the SNPs, but C^∗03:04 (OR = 0.72, p = 8 × 10⁻⁶) and DRB1^∗04:04 (OR = 0.65, p = 4 × 10⁻⁵) remained significant. Similarly, rs3129882 and the closely linked rs9268515 and rs2395163 remained significant irrespective of HLA alleles. rs3129882 and rs2395163 are expression quantitative trait loci (eQTLs) for HLA-DR and HLA-DQ (9 × 10⁻⁵ ≥ P_eQTL ≥ 2 × 10⁻⁷⁹), suggesting that HLA gene expression might influence PD. Our data suggest that PD is associated with both structural and regulatory elements in HLA. Furthermore, our study demonstrates that noncoding SNPs in the HLA region can be associated with disease irrespective of HLA alleles, and that observed associations with HLA alleles can sometimes be secondary to a noncoding variant.     

A genome-wide association study follow-up suggests a possible role for PPARG in systemic sclerosis susceptibility

Introduction

A recent genome-wide association study (GWAS) comprising a French cohort of systemic sclerosis (SSc) reported several non-HLA single-nucleotide polymorphisms (SNPs) showing a nominal association in the discovery phase. We aimed to identify previously overlooked susceptibility variants by using a follow-up strategy.

Methods

Sixty-six non-HLA SNPs showing a P value <10^-4 in the discovery phase of the French SSc GWAS were analyzed in the first step of this study, performing a meta-analysis that combined data from the two published SSc GWASs. A total of 2,921 SSc patients and 6,963 healthy controls were included in this first phase. Two SNPs, PPARG rs310746 and CHRNA9 rs6832151, were selected for genotyping in the replication cohort (1,068 SSc patients and 6,762 healthy controls) based on the results of the first step. Genotyping was performed by using TaqMan SNP genotyping assays.

Results

We observed nominal associations for both PPARG rs310746 (P_MH = 1.90 × 10^-6, OR, 1.28) and CHRNA9 rs6832151 (P_MH = 4.30 × 10^-6, OR, 1.17) genetic variants with SSc in the first step of our study. In the replication phase, we observed a trend of association for PPARG rs310746 (P value = 0.066; OR, 1.17). The combined overall Mantel-Haenszel meta-analysis of all the cohorts included in the present study revealed that PPARG rs310746 remained associated with SSc with a nominal non-genome-wide significant P value (P_MH = 5.00 × 10^-7; OR, 1.25). No evidence of association was observed for CHRNA9 rs6832151 either in the replication phase or in the overall pooled analysis.

Conclusion

Our results suggest a role of PPARG gene in the development of SSc.     

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.