Evidence for more than one Parkinson's disease-associated variant within the HLA region

Parkinson''s disease (PD) was recently found to be associated with HLA in a genome-wide association study (GWAS). Follow-up GWAS''s replicated the PD-HLA association but their top hits differ. Do the different hits tag the same locus or is there more than one PD-associated variant within HLA? We show that the top GWAS hits are not correlated with each other (0.00≤r²≤0.15). Using our GWAS (2000 cases, 1986 controls) we conducted step-wise conditional analysis on 107 SNPs with P<10⁻³ for PD-association; 103 dropped-out, four remained significant. Each SNP, when conditioned on the other three, yielded P_SNP1 = 5×10⁻⁴, P_SNP2 = 5×10⁻⁴, P_SNP3 = 4×10⁻³ and P_SNP4 = 0.025. The four SNPs were not correlated (0.01≤r²≤0.20). Haplotype analysis (excluding rare SNP2) revealed increasing PD risk with increasing risk alleles from OR = 1.27, P = 5×10⁻³ for one risk allele to OR = 1.65, P = 4×10⁻⁸ for three. Using additional 843 cases and 856 controls we replicated the independent effects of SNP1 (P_{conditioned-on-SNP4} = 0.04) and SNP4 (P_{conditioned-on-SNP1} = 0.04); SNP2 and SNP3 could not be replicated. In pooled GWAS and replication, SNP1 had OR_{conditioned-on-SNP4} = 1.23, P_{conditioned-on-SNP4} = 6×10⁻⁷; SNP4 had OR_{conditioned-on-SNP1} = 1.18, P_{conditioned-on-SNP1} = 3×10⁻³; and the haplotype with both risk alleles had OR = 1.48, P = 2×10⁻¹². Genotypic OR increased with the number of risk alleles an individual possessed up to OR = 1.94, P = 2×10⁻¹¹ for individuals who were homozygous for the risk allele at both SNP1 and SNP4. SNP1 is a variant in HLA-DRA and is associated with HLA-DRA, DRB5 and DQA2 gene expression. SNP4 is correlated (r² = 0.95) with variants that are associated with HLA-DQA2 expression, and with the top HLA SNP from the IPDGC GWAS (r² = 0.60). Our findings suggest more than one PD-HLA association; either different alleles of the same gene, or separate loci.     

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.     

TRAIL/TRAIL receptor system and susceptibility to multiple sclerosis

The TNF-related apoptosis inducing ligand (TRAIL)/TRAIL receptor system participates in crucial steps in immune cell activation or differentiation. It is able to inhibit proliferation and activation of T cells and to induce apoptosis of neurons and oligodendrocytes, and seems to be implicated in autoimmune diseases. Thus, TRAIL and TRAIL receptor genes are potential candidates for involvement in susceptibility to multiple sclerosis (MS). To test whether single-nucleotide polymorphisms (SNPs) in the human genes encoding TRAIL, TRAILR-1, TRAILR-2, TRAILR-3 and TRAILR-4 are associated with MS susceptibility, we performed a candidate gene case-control study in the Spanish population. 59 SNPs in the TRAIL and TRAIL receptor genes were analysed in 628 MS patients and 660 controls, and validated in an additional cohort of 295 MS patients and 233 controls. Despite none of the SNPs withstood the highly conservative Bonferroni correction, three SNPs showing uncorrected p values<0.05 were successfully replicated: rs4894559 in TRAIL gene, p = 9.8×10⁻⁴, OR = 1.34; rs4872077, in TRAILR-1 gene, p = 0.005, OR = 1.72; and rs1001793 in TRAILR-2 gene, p = 0.012, OR = 0.84. The combination of the alleles G/T/A in these SNPs appears to be associated with a reduced risk of developing MS (p = 2.12×10⁻⁵, OR = 0.59). These results suggest that genes of the TRAIL/TRAIL receptor system exerts a genetic influence on MS.     

Common variants at 9p21 and 8q22 are associated with increased susceptibility to optic nerve degeneration in glaucoma

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63–0.75], p = 1.86×10⁻¹⁸), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21–1.43], p = 3.87×10⁻¹¹). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50–0.67], p = 1.17×10⁻¹²) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53–0.72], p = 8.88×10⁻¹⁰). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41–0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54–1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.     

FTO gene associated fatness in relation to body fat distribution and metabolic traits throughout a broad range of fatness

Background

A common single nucleotide polymorphism (SNP) of FTO (rs9939609, T/A) is associated with total body fatness. We investigated the association of this SNP with abdominal and peripheral fatness and obesity-related metabolic traits in middle-aged men through a broad range of fatness present already in adolescence.

Methodology/Principal Findings

Obese young Danish men (n = 753, BMI≥31.0 kg/m²) and a randomly selected group (n = 879) from the same population were examined in three surveys (mean age 35, 46 and 49 years, respectively). The traits included anthropometrics, body composition, oral glucose tolerance test, blood lipids, blood pressure, fibrinogen and aspartate aminotransferase. Logistic regression analysis was used to assess the age-adjusted association between the phenotypes and the odds ratios for the FTO rs9939609 (TT and TA genotype versus the AA genotype), for anthropometrics and body composition estimated per unit z-score. BMI was strongly associated with the AA genotype in all three surveys: OR = 1.17, p = 1.1*10⁻⁶, OR = 1.20, p = 1.7*10⁻⁷, OR = 1.17, p = 3.4*10⁻³, respectively. Fat body mass index was also associated with the AA genotype (OR = 1.21, p = 4.6*10⁻⁷ and OR = 1.21, p = 1.0*10⁻³). Increased abdominal fatness was associated with the AA genotype when measured as waist circumference (OR = 1.21, p = 2.2*10⁻⁶ and OR = 1.19, p = 5.9*10⁻³), sagittal abdominal diameter (OR = 1.17, p = 1.3*10⁻⁴ and OR = 1.18, p = 0.011) and intra-abdominal adipose tissue (OR = 1.21, p = 0.005). Increased peripheral fatness measured as hip circumference (OR = 1.19, p = 1.3*10⁻⁵ and OR = 1.18, p = 0.004) and lower body fat mass (OR = 1.26, p = 0.002) was associated with the AA genotype. The AA genotype was significantly associated with decreased Stumvoll insulin sensitivity index (OR = 0.93, p = 0.02) and with decreased non-fasting plasma HDL-cholesterol (OR = 0.57, p = 0.037), but not with any other of the metabolic traits. However, all significant results for both body fat distribution and metabolic traits were explained by a mediating effect of total fat mass.

Conclusion

The association of the examined FTO SNP to general fatness throughout the range of fatness was confirmed, and this association explains the relation between the SNP and body fat distribution and decreased insulin sensitivity and HDL-cholesterol. The SNP was not significantly associated with other metabolic traits suggesting that they are not derived from the general accumulation of body fat.     

Genetic association study identifies HSPB7 as a risk gene for idiopathic dilated cardiomyopathy

Dilated cardiomyopathy (DCM) is a structural heart disease with strong genetic background. Monogenic forms of DCM are observed in families with mutations located mostly in genes encoding structural and sarcomeric proteins. However, strong evidence suggests that genetic factors also affect the susceptibility to idiopathic DCM. To identify risk alleles for non-familial forms of DCM, we carried out a case-control association study, genotyping 664 DCM cases and 1,874 population-based healthy controls from Germany using a 50K human cardiovascular disease bead chip covering more than 2,000 genes pre-selected for cardiovascular relevance. After quality control, 30,920 single nucleotide polymorphisms (SNP) were tested for association with the disease by logistic regression adjusted for gender, and results were genomic-control corrected. The analysis revealed a significant association between a SNP in HSPB7 gene (rs1739843, minor allele frequency 39%) and idiopathic DCM (p = 1.06×10⁻⁶, OR = 0.67 [95% CI 0.57–0.79] for the minor allele T). Three more SNPs showed p < 2.21×10⁻⁵. De novo genotyping of these four SNPs was done in three independent case-control studies of idiopathic DCM. Association between SNP rs1739843 and DCM was significant in all replication samples: Germany (n = 564, n = 981 controls, p = 2.07×10⁻³, OR = 0.79 [95% CI 0.67–0.92]), France 1 (n = 433 cases, n = 395 controls, p = 3.73×10⁻³, OR = 0.74 [95% CI 0.60–0.91]), and France 2 (n = 249 cases, n = 380 controls, p = 2.26×10⁻⁴, OR = 0.63 [95% CI 0.50–0.81]). The combined analysis of all four studies including a total of n = 1,910 cases and n = 3,630 controls showed highly significant evidence for association between rs1739843 and idiopathic DCM (p = 5.28×10⁻¹³, OR = 0.72 [95% CI 0.65–0.78]). None of the other three SNPs showed significant results in the replication stage.This finding of the HSPB7 gene from a genetic search for idiopathic DCM using a large SNP panel underscores the influence of common polymorphisms on DCM susceptibility.     

Genetic determinants of cardiovascular events among women with migraine: a genome-wide association study

Background

Migraine is associated with an increased risk for cardiovascular disease (CVD). Both migraine and CVD are highly heritable. However, the genetic liability for CVD among migraineurs is unclear.

Methods

We performed a genome-wide association study for incident CVD events during 12 years of follow-up among 5,122 migraineurs participating in the population-based Women''s Genome Health Study. Migraine was self-reported and CVD events were confirmed after medical records review. We calculated odds ratios (OR) and 95% confidence intervals (CI) and considered a genome-wide p-value <5×10⁻⁸ as significant.

Results

Among the 5,122 women with migraine 164 incident CVD events occurred during follow-up. No SNP was associated with major CVD, ischemic stroke, myocardial infarction, or CVD death at the genome-wide level; however, five SNPs showed association with p<5×10⁻⁶. Among migraineurs with aura rs7698623 in MEPE (OR = 6.37; 95% CI 3.15–12.90; p = 2.7×10⁻⁷) and rs4975709 in IRX4 (OR = 5.06; 95% CI 2.66–9.62; p = 7.7×10⁻⁷) appeared to be associated with ischemic stroke, rs2143678 located close to MDF1 with major CVD (OR = 3.05; 95% CI 1.98–4.69; p = 4.3×10⁻⁷), and the intergenic rs1406961 with CVD death (OR = 12.33; 95% CI 4.62–32.87; p = 5.2×10⁻⁷). Further, rs1047964 in BACE1 appeared to be associated with CVD death among women with any migraine (OR = 4.67; 95% CI 2.53–8.62; p = 8.0×10⁻⁷).

Conclusion

Our results provide some suggestion for an association of five SNPs with CVD events among women with migraine; none of the results was genome-wide significant. Four associations appeared among migraineurs with aura, two of those with ischemic stroke. Although our population is among the largest with migraine and incident CVD information, these results must be treated with caution, given the limited number of CVD events among women with migraine and the low minor allele frequencies for three of the SNPs. Our results await independent replication and should be considered hypothesis generating for future research.     

Genetic variants at chromosomes 2q35, 5p12, 6q25.1, 10q26.13, and 16q12.1 influence the risk of breast cancer in men

Male breast cancer accounts for approximately 1% of all breast cancer. To date, risk factors for male breast cancer are poorly defined, but certain risk factors and genetic features appear common to both male and female breast cancer. Genome-wide association studies (GWAS) have recently identified common single nucleotide polymorphisms (SNPs) that influence female breast cancer risk; 12 of these have been independently replicated. To examine if these variants contribute to male breast cancer risk, we genotyped 433 male breast cancer cases and 1,569 controls. Five SNPs showed a statistically significant association with male breast cancer: rs13387042 (2q35) (odds ratio (OR)  = 1.30, p = 7.98×10⁻⁴), rs10941679 (5p12) (OR = 1.26, p = 0.007), rs9383938 (6q25.1) (OR = 1.39, p = 0.004), rs2981579 (FGFR2) (OR = 1.18, p = 0.03), and rs3803662 (TOX3) (OR = 1.48, p = 4.04×10⁻⁶). Comparing the ORs for male breast cancer with the published ORs for female breast cancer, three SNPs—rs13387042 (2q35), rs3803662 (TOX3), and rs6504950 (COX11)—showed significant differences in ORs (p<0.05) between sexes. Breast cancer is a heterogeneous disease; the relative risks associated with loci identified to date show subtype and, based on these data, gender specificity. Additional studies of well-defined patient subgroups could provide further insight into the biological basis of breast cancer development.     

Association of six single nucleotide polymorphisms with gestational diabetes mellitus in a Chinese population

Background

To investigate whether the candidate genes that confer susceptibility to type 2 diabetes mellitus are also correlated with gestational diabetes mellitus (GDM) in pregnant Chinese women.

Methodology/Principal Findings

In this study, 1764 unrelated pregnant women were recruited, of which 725 women had GDM and 1039 served as controls. Six single nucleotide polymorphisms (rs7754840 in CDKAL1, rs391300 in SRR, rs2383208 in CDKN2A/2B, rs4402960 in IGF2BP2, rs10830963 in MTNR1B, rs4607517 in GCK) were genotyped using TaqMan allelic discrimination assays. The genotype and allele distributions of each SNP between the GDM cases and controls and the combined effects of alleles for the risk of developing GDM were analyzed. We found that the rs4402960, rs2383208 and rs391300 were statistically associated with GDM (OR = 1.207, 95%CI = 1.029–1.417, p = 0.021; OR = 1.242, 95%CI = 1.077–1.432, p = 0.003; OR = 1.202, 95%CI = 1.020–1.416, P = 0.028, respectively). In addition, the effect was greater under a recessive model in rs391300 (OR = 1.820, 95%CI = 1.226–2.701, p = 0.003). Meanwhile, the joint effect of these three loci indicated an additive effect of multiple alleles on the risk of developing GDM with an OR of 1.196 per allele (p = 1.08×10⁻⁴). We also found that the risk alleles of rs2383208 (b = −0.085, p = 0.003), rs4402960 (b = −0.057, p = 0.046) and rs10830963 (b = −0.096, p = 0.001) were associated with HOMA-B, while rs7754840 was associated with decrease in insulin AUC during a 100 g OGTT given at the time of GDM diagnosis (b = −0.080, p = 0.007).

Conclusions/Significance

Several risk alleles of type 2 diabetes were associated with GDM in pregnant Chinese women. The effects of these SNPs on GDM might be through the impairment of beta cell function and these risk loci contributed additively to the disease.     

HLA-DRB1-DQB1 haplotypes confer susceptibility and resistance to multiple sclerosis in Sardinia

Introduction

Genetic predisposition to multiple sclerosis (MS) in Sardinia (Italy) has been associated with five DRB1*-DQB1* haplotypes of the human leukocyte antigen (HLA). Given the complexity of these associations, an in-depth re-analysis was performed with the specific aims of confirming the haplotype associations; establishing the independence of the associated haplotypes; and assessing patients'' genotypic risk of developing MS.

Methods and Results

A transmission disequilibrium test (TDT) of the DRB1*-DQB1* haplotypes in 943 trio families, confirmed a higher than expected transmission rate (over-transmission) of the *13:03-*03:01 (OR = 2.9, P = 7.6×10⁻³), *04:05-*03:01 (OR = 2.4, P = 4.4×10⁻⁶) and *03:01-*02:01 (OR = 2.1, P = 1.0×10⁻¹⁵) haplotype. In contrast, the *16:01-*05:02 (OR = 0.5, P = 5.4×10⁻¹¹) and the *15:02-*06:01 (OR = 0.3, P = 1.5×10⁻³) haplotypes exhibited a lower than expected transmission rate (under-transmission). The independence of the transmission of each positively and negatively associated haplotype was confirmed relative to all positively associated haplotypes, and to the negatively associated *16:01-*05:02 haplotype. In patients, carriage of two predisposing haplotypes, or of protective haplotypes, respectively increased or decreased the patient''s risk of developing MS. The risk of MS followed a multiplicative model of genotypes, which was, in order of decreasing ORs: *04:05-*0301/*03:01-*02:01 (OR = 4.5); *03:01-*02:01/*03:01-*02:01 (OR = 4.1); and the *16:01-*05:02/*16:01-*0502 (OR = 0.2) genotypes. Analysis of DRB1 and DQB1 protein chain residues showed that the Val/Gly residue at position 86 of the DRB1 chain was the only difference between the protective *16:01- *15:02 alleles and the predisposing *15:01 one. Similarly, the Ala/Val residue at position 38 of the DQB1 chain differentiated the positively associated *06:02 allele and the negatively associated *05:02, *06:01 alleles.

Conclusions

These findings show that the association of specific, independent DRB1*-DQB1* haplotypes confers susceptibility or resistance to MS in the MS-prone Sardinian population. The data also supports a functional role for specific residues of the DRB1 and DQB1 proteins in predisposing patients to MS.     

Identification of a novel idiopathic epilepsy locus in Belgian Shepherd dogs

Epilepsy is the most common neurological disorder in dogs, with an incidence ranging from 0.5% to up to 20% in particular breeds. Canine epilepsy can be etiologically defined as idiopathic or symptomatic. Epileptic seizures may be classified as focal with or without secondary generalization, or as primary generalized. Nine genes have been identified for symptomatic (storage diseases) and one for idiopathic epilepsy in different breeds. However, the genetic background of common canine epilepsies remains unknown. We have studied the clinical and genetic background of epilepsy in Belgian Shepherds. We collected 159 cases and 148 controls and confirmed the presence of epilepsy through epilepsy questionnaires and clinical examinations. The MRI was normal while interictal EEG revealed abnormalities and variable foci in the clinically examined affected dogs. A genome-wide association study using Affymetrix 50K SNP arrays in 40 cases and 44 controls mapped the epilepsy locus on CFA37, which was replicated in an independent cohort (81 cases and 88 controls; combined p = 9.70×10⁻¹⁰, OR = 3.3). Fine mapping study defined a ∼1 Mb region including 12 genes of which none are known epilepsy genes or encode ion channels. Exonic sequencing was performed for two candidate genes, KLF7 and ADAM23. No variation was found in KLF7 but a highly-associated non-synonymous variant, G1203A (R387H) was present in the ADAM23 gene (p = 3.7×10⁻⁸, OR = 3.9 for homozygosity). Homozygosity for a two-SNP haplotype within the ADAM23 gene conferred the highest risk for epilepsy (p = 6.28×10⁻¹¹, OR = 7.4). ADAM23 interacts with known epilepsy proteins LGI1 and LGI2. However, our data suggests that the ADAM23 variant is a polymorphism and we have initiated a targeted re-sequencing study across the locus to identify the causative mutation. It would establish the affected breed as a novel therapeutic model, help to develop a DNA test for breeding purposes and introduce a novel candidate gene for human idiopathic epilepsies.     

GWAS identifies novel susceptibility loci on 6p21.32 and 21q21.3 for hepatocellular carcinoma in chronic hepatitis B virus carriers

Genome-wide association studies (GWAS) have recently identified KIF1B as susceptibility locus for hepatitis B virus (HBV)–related hepatocellular carcinoma (HCC). To further identify novel susceptibility loci associated with HBV–related HCC and replicate the previously reported association, we performed a large three-stage GWAS in the Han Chinese population. 523,663 autosomal SNPs in 1,538 HBV–positive HCC patients and 1,465 chronic HBV carriers were genotyped for the discovery stage. Top candidate SNPs were genotyped in the initial validation samples of 2,112 HBV–positive HCC cases and 2,208 HBV carriers and then in the second validation samples of 1,021 cases and 1,491 HBV carriers. We discovered two novel associations at rs9272105 (HLA-DQA1/DRB1) on 6p21.32 (OR = 1.30, P = 1.13×10⁻¹⁹) and rs455804 (GRIK1) on 21q21.3 (OR = 0.84, P = 1.86×10⁻⁸), which were further replicated in the fourth independent sample of 1,298 cases and 1,026 controls (rs9272105: OR = 1.25, P = 1.71×10⁻⁴; rs455804: OR = 0.84, P = 6.92×10⁻³). We also revealed the associations of HLA-DRB1*0405 and 0901*0602, which could partially account for the association at rs9272105. The association at rs455804 implicates GRIK1 as a novel susceptibility gene for HBV–related HCC, suggesting the involvement of glutamate signaling in the development of HBV–related HCC.     

Lysyl oxidase polymorphisms and susceptibility to osteosarcoma

Despite the knowledge of many genetic alterations present in osteosarcoma, the complexity of this disease precludes placing its biology into a simple conceptual framework. Lysyl oxidase (LOX) catalyzes the cross-linking of elastin and collagen, which is essential for the structural integrity and function of bone tissue. In the current study, we performed genomic sequencing on all seven exons -including the intron-exon splice sites, and the putative promoter region of LOX gene - followed by luciferase reporter assay to analyze the function of newly identified polymorphisms. Associations between LOX polymorphisms and osteosarcoma were then evaluated. Our sequencing data revealed three polymorphisms (−22G/C, 225C/G, and 473G/A) in the exons and promoter region of LOX. The −22G/C polymorphism lies in the downstream core promoter element (DPE) region and caused a decrease in promoter activity of LOX. The prevalence of the −22C allele and 473A allele were significantly increased in osteosarcoma patients compared to controls (odds ratio [OR] = 3.88, 95% confidence interval [CI]  = 1.94−7.78, p = 4.18×10⁻⁵, and OR = 1.38, 95%CI = 1.07−1.78, p = 0.013; p 0.0167 was considered significant after Bonferroni correction). Analyzing haplotype showed that the frequency of CCG haplotype (−22, 225, 473) was significantly higher in osteosarcoma cases than in healthy controls after Bonferroni correction (p = 4.46×10⁻⁴). These results indicate that the −22G/C polymorphism may affect the expression of LOX, and that −22G/C and 473G/A polymorphisms may be new risk factors for osteosarcoma. These findings reveal a potential new pathway by which genetic polymorphisms may affect human diseases.     

Identification of novel genetic markers associated with clinical phenotypes of systemic sclerosis through a genome-wide association strategy

The aim of this study was to determine, through a genome-wide association study (GWAS), the genetic components contributing to different clinical sub-phenotypes of systemic sclerosis (SSc). We considered limited (lcSSc) and diffuse (dcSSc) cutaneous involvement, and the relationships with presence of the SSc-specific auto-antibodies, anti-centromere (ACA), and anti-topoisomerase I (ATA). Four GWAS cohorts, comprising 2,296 SSc patients and 5,171 healthy controls, were meta-analyzed looking for associations in the selected subgroups. Eighteen polymorphisms were further tested in nine independent cohorts comprising an additional 3,175 SSc patients and 4,971 controls. Conditional analysis for associated SNPs in the HLA region was performed to explore their independent association in antibody subgroups. Overall analysis showed that non-HLA polymorphism rs11642873 in IRF8 gene to be associated at GWAS level with lcSSc (P = 2.32×10⁻¹², OR = 0.75). Also, rs12540874 in GRB10 gene (P = 1.27 × 10⁻⁶, OR = 1.15) and rs11047102 in SOX5 gene (P = 1.39×10⁻⁷, OR = 1.36) showed a suggestive association with lcSSc and ACA subgroups respectively. In the HLA region, we observed highly associated allelic combinations in the HLA-DQB1 locus with ACA (P = 1.79×10⁻⁶¹, OR = 2.48), in the HLA-DPA1/B1 loci with ATA (P = 4.57×10⁻⁷⁶, OR = 8.84), and in NOTCH4 with ACA P = 8.84×10⁻²¹, OR = 0.55) and ATA (P = 1.14×10⁻⁸, OR = 0.54). We have identified three new non-HLA genes (IRF8, GRB10, and SOX5) associated with SSc clinical and auto-antibody subgroups. Within the HLA region, HLA-DQB1, HLA-DPA1/B1, and NOTCH4 associations with SSc are likely confined to specific auto-antibodies. These data emphasize the differential genetic components of subphenotypes of SSc.     

A functional variant in ERAP1 predisposes to multiple sclerosis

The ERAP1 gene encodes an aminopeptidase involved in antigen processing. A functional polymorphism in the gene (rs30187, Arg528Lys) associates with susceptibility to ankylosying spondylitis (AS), whereas a SNP in the interacting ERAP2 gene increases susceptibility to another inflammatory autoimmune disorder, Crohn''s disease (CD). We analysed rs30187 in 572 Italian patients with CD and in 517 subjects suffering from multiple sclerosis (MS); for each cohort, an independent sex- and age-matched control group was genotyped. The frequency of the 528Arg allele was significantly higher in both disease cohorts compared to the respective control population (for CD, OR = 1.20 95%CI: 1.01–1.43, p = 0.036; for RRMS, OR = 1.26; 95%CI: 1.04–1.51, p = 0.01). Meta-analysis with the Wellcome Trust Cases Control Consortium GWAS data confirmed the association with MS (p_meta = 0.005), but not with CD. In AS, the rs30187 variant has a predisposing effect only in an HLA-B27 allelic background. It remains to be evaluated whether interaction between ERAP1 and distinct HLA class I alleles also affects the predisposition to MS, and explains the failure to provide definitive evidence for a role of rs30187 in CD. Results herein support the emerging concept that a subset of master-regulatory genes underlay the pathogenesis of autoimmunity.     

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.     

Risk alleles for systemic lupus erythematosus in a large case-control collection and associations with clinical subphenotypes

Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. Recent studies have greatly expanded the number of established SLE risk alleles, but the distribution of multiple risk alleles in cases versus controls and their relationship to subphenotypes have not been studied. We studied 22 SLE susceptibility polymorphisms with previous genome-wide evidence of association (p<5×10⁻⁸) in 1919 SLE cases from 9 independent Caucasian SLE case series and 4813 independent controls. The mean number of risk alleles in cases was 15.1 (SD 3.1) while the mean in controls was 13.1 (SD 2.8), with trend p = 4×10⁻¹²⁸. We defined a genetic risk score (GRS) for SLE as the number of risk alleles with each weighted by the SLE risk odds ratio (OR). The OR for high-low GRS tertiles, adjusted for intra-European ancestry, sex, and parent study, was 4.4 (95% CI 3.8–5.1). We studied associations of individual SNPs and the GRS with clinical manifestations for the cases: age at diagnosis, the 11 American College of Rheumatology classification criteria, and double-stranded DNA antibody (anti-dsDNA) production. Six subphenotypes were significantly associated with the GRS, most notably anti-dsDNA (OR_high-low = 2.36, p = 9e−9), the immunologic criterion (OR_high-low = 2.23, p = 3e−7), and age at diagnosis (OR_high-low = 1.45, p = 0.0060). Finally, we developed a subphenotype-specific GRS (sub-GRS) for each phenotype with more power to detect cumulative genetic associations. The sub-GRS was more strongly associated than any single SNP effect for 5 subphenotypes (the above plus hematologic disorder and oral ulcers), while single loci are more significantly associated with renal disease (HLA-DRB1, OR = 1.37, 95% CI 1.14–1.64) and arthritis (ITGAM, OR = 0.72, 95% CI 0.59–0.88). We did not observe significant associations for other subphenotypes, for individual loci or the sub-GRS. Thus our analysis categorizes SLE subphenotypes into three groups: those having cumulative, single, and no known genetic association with respect to the currently established SLE risk loci.     

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 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.