Protective Role of the Interleukin 33 rs3939286 Gene Polymorphism in the Development of Subclinical Atherosclerosis in Rheumatoid Arthritis Patients

Objectives

To determine whether the interleukin-33 (IL-33)-interleukin-1 receptor like 1 (IL-1RL1) signaling pathway is implicated in the risk of subclinical atherosclerosis in patients with rheumatoid arthritis (RA).

Methods

A total of 576 Spanish RA patients from Northern Spain were genotyped for 6 well-known IL33-IL1RL1 polymorphisms (IL33 rs3939286, IL33 rs7025417, IL33 rs7044343, IL1RL1 rs2058660, IL1RL1 rs2310173 and IL1RL1 rs13015714) by TaqMan genotyping assay. The presence of subclinical atherosclerosis was determined by the assessment of carotid intima-media thickness (cIMT) by carotid ultrasound (US).

Results

RA patients carrying the TT genotype of the IL33 rs3939286 polymorphism had lower cIMT values than those homozygous for the CC genotype (mean ± standard deviation (SD): 0.71 ± 0.14 mm versus 0.76 ± 0.16 mm, respectively) while patients carrying the CT genotype had intermediate cIMT values (mean ± SD: 0.73 ± 0.17 mm). Moreover, RA patients carrying the mutant allele T of the IL33 rs3939286 polymorphism exhibited significantly lower cIMT values than those carrying the wild allele C (mean ± SD: 0.72 ± 0.16 mm versus 0.75 ± 0.18 mm respectively; p = 0.04). The association of both genotype and allele frequencies of IL33 rs3939286 and cIMT levels remained statistically significant after adjustment for sex, age at the time of US study, follow-up and center (p = 0.006 and p = 0.0023, respectively), evidencing that the potential effect conferred by IL33 rs3939286 may be independent of confounder factors. No association with other IL33-IL1RL1 genetic variants was observed.

Conclusions

In conclusion, our results may suggest a potential protective effect of the IL33 rs3939286 allele T in the risk of subclinical atherosclerosis in patients with RA.     

Identification of IRF8, TMEM39A, and IKZF3-ZPBP2 as susceptibility loci for systemic lupus erythematosus in a large-scale multiracial replication study

Systemic lupus erythematosus (SLE) is a chronic heterogeneous autoimmune disorder characterized by the loss of tolerance to self-antigens and dysregulated interferon responses. The etiology of SLE is complex, involving both heritable and environmental factors. Candidate-gene studies and genome-wide association (GWA) scans have been successful in identifying new loci that contribute to disease susceptibility; however, much of the heritable risk has yet to be identified. In this study, we sought to replicate 1,580 variants showing suggestive association with SLE in a previously published GWA scan of European Americans; we tested a multiethnic population consisting of 7,998 SLE cases and 7,492 controls of European, African American, Asian, Hispanic, Gullah, and Amerindian ancestry to find association with the disease. Several genes relevant to immunological pathways showed association with SLE. Three loci exceeded the genome-wide significance threshold: interferon regulatory factor 8 (IRF8; rs11644034; p_meta-Euro = 2.08 × 10⁻¹⁰), transmembrane protein 39A (TMEM39A; rs1132200; p_meta-all = 8.62 × 10⁻⁹), and 17q21 (rs1453560; p_meta-all = 3.48 × 10⁻¹⁰) between IKAROS family of zinc finger 3 (AIOLOS; IKZF3) and zona pellucida binding protein 2 (ZPBP2). Fine mapping, resequencing, imputation, and haplotype analysis of IRF8 indicated that three independent effects tagged by rs8046526, rs450443, and rs4843869, respectively, were required for risk in individuals of European ancestry. Eleven additional replicated effects (5 × 10⁻⁸ < p_meta-Euro < 9.99 × 10⁻⁵) were observed with CFHR1, CADM2, LOC730109/IL12A, LPP, LOC63920, SLU7, ADAMTSL1, C10orf64, OR8D4, FAM19A2, and STXBP6. The results of this study increase the number of confirmed SLE risk loci and identify others warranting further investigation.     

A Candidate Gene Approach Identifies an IL33 Genetic Variant as a Novel Genetic Risk Factor for GCA

Introduction

Increased expression of IL-33 and its receptor ST2, encoded by the IL1RL1 gene, has been detected in the inflamed arteries of giant cell arteritis (GCA) patients. The aim of the present study was to investigate for the first time the potential influence of the IL33 and IL1RL1 loci on GCA predisposition.

Methods

A total of 1,363 biopsy-proven GCA patients and 3,908 healthy controls from four European cohorts (Spain, Italy, Germany and Norway) were combined in a meta-analysis. Six genetic variants: rs3939286, rs7025417 and rs7044343, within the IL33 gene, and rs2058660, rs2310173 and rs13015714, within the IL1RL1 gene, previously associated with immune-related diseases, were genotyped using predesigned TaqMan assays.

Results

A consistent association between the rs7025417 polymorphism and GCA was evident in the overall meta-analysis, under both allele (P_MH = 0.041, OR = 0.88, CI 95% 0.78–0.99) and recessive (P_MH = 3.40E-03, OR = 0.53, CI 95% 0.35–0.80) models. No statistically significant differences between allele or genotype frequencies for the other IL33 and IL1RL1 genetic variants were detected in this pooled analysis.

Conclusions

Our results clearly evidenced the implication of the IL33 rs7025417 polymorphism in the genetic network underlying GCA.     

Association of genetic variants in enamel-formation genes with dental caries: A meta- and gene-cluster analysis

Previous studies have reported the association between multiple genetic variants in the enamel-formation genes and the risk of dental caries with inconsistent results. We performed a systematic literature search of the PubMed, Cochrane Library, HuGE and Google Scholar databases for studies published before March 21, 2020 and conducted meta-, gene-based and gene-cluster analysis on the association between genetic variants in the enamel-formation genes and the risk of dental caries. We identified 21 relevant publications including a total of 24 studies for analysis. The genetic variant rs17878486 in AMELX was significantly associated with dental caries risk (OR = 1.40, 95% CI: 1.02–1.93, P = 0.037). We found no significant association between the risk of dental caries with rs12640848 in ENAM (OR = 1.15, 95% CI: 0.88–1.52, P = 0.310), rs1784418 in MMP20 (OR = 1.07, 95% CI: 0.76–1.49, P = 0.702) and rs3796704 in ENAM (OR = 1.06, 95% CI: 0.96–1.17, P = 0.228). Gene-based analysis indicated that multiple genetic variants in AMELX showed joint association with the risk of dental caries (6 variants; P < 10⁻⁵), so did genetic variants in MMP13 (3 variants; P = 0.004), MMP2 (3 variants; P < 10⁻⁵), MMP20 (2 variants; P < 10⁻⁵) and MMP3 (2 variants; P < 10⁻⁵). The gene-cluster analysis indicated a significant association between the genetic variants in this enamel-formation gene cluster and the risk of dental caries (P < 10⁻⁵). The present meta-analysis revealed that genetic variant rs17878486 in AMELX was associated with dental caries, and multiple genetic variants in the enamel-formation genes jointly contributed to the risk of dental caries, supporting the role of genetic variants in the enamel-formation genes in the etiology of dental caries.     

A genome-wide association study reveals that variants within the HLA region are associated with risk for nonobstructive azoospermia

A genome-wide association study of Han Chinese subjects was conducted to identify genetic susceptibility loci for nonobstructive azoospermia (NOA). In the discovery stage, 802 azoospermia cases and 1,863 controls were screened for genetic variants in the genome. Promising SNPs were subsequently confirmed in two independent sets of subjects: 818 azoospermia cases and 1,755 controls from northern China, and 606 azoospermia cases and 958 controls from central and southern China. We detected variants at human leukocyte antigen (HLA) regions that were independently associated with NOA (HLA-DRA, rs3129878, p_combine = 3.70 × 10⁻¹⁶, odds ratio [OR] = 1.37; C6orf10 and BTNL2, rs498422, p_combine = 2.43 × 10⁻¹², OR = 1.42). These findings provide additional insight into the pathogenesis of NOA.     

Variants near FOXE1 are associated with hypothyroidism and other thyroid conditions: using electronic medical records for genome- and phenome-wide studies

We repurposed existing genotypes in DNA biobanks across the Electronic Medical Records and Genomics network to perform a genome-wide association study for primary hypothyroidism, the most common thyroid disease. Electronic selection algorithms incorporating billing codes, laboratory values, text queries, and medication records identified 1317 cases and 5053 controls of European ancestry within five electronic medical records (EMRs); the algorithms'' positive predictive values were 92.4% and 98.5% for cases and controls, respectively. Four single-nucleotide polymorphisms (SNPs) in linkage disequilibrium at 9q22 near FOXE1 were associated with hypothyroidism at genome-wide significance, the strongest being rs7850258 (odds ratio [OR] 0.74, p = 3.96 × 10⁻⁹). This association was replicated in a set of 263 cases and 1616 controls (OR = 0.60, p = 5.7 × 10⁻⁶). A phenome-wide association study (PheWAS) that was performed on this locus with 13,617 individuals and more than 200,000 patient-years of billing data identified associations with additional phenotypes: thyroiditis (OR = 0.58, p = 1.4 × 10⁻⁵), nodular (OR = 0.76, p = 3.1 × 10⁻⁵) and multinodular (OR = 0.69, p = 3.9 × 10⁻⁵) goiters, and thyrotoxicosis (OR = 0.76, p = 1.5 × 10⁻³), but not Graves disease (OR = 1.03, p = 0.82). Thyroid cancer, previously associated with this locus, was not significantly associated in the PheWAS (OR = 1.29, p = 0.09). The strongest association in the PheWAS was hypothyroidism (OR = 0.76, p = 2.7 × 10⁻¹³), which had an odds ratio that was nearly identical to that of the curated case-control population in the primary analysis, providing further validation of the PheWAS method. Our findings indicate that EMR-linked genomic data could allow discovery of genes associated with many diseases without additional genotyping cost.     

Combined analysis of genome-wide association studies for Crohn disease and psoriasis identifies seven shared susceptibility loci

Psoriasis (PS) and Crohn disease (CD) have been shown to be epidemiologically, pathologically, and therapeutically connected, but little is known about their shared genetic causes. We performed meta-analyses of five published genome-wide association studies on PS (2,529 cases and 4,955 controls) and CD (2,142 cases and 5,505 controls), followed up 20 loci that showed strongest evidence for shared disease association and, furthermore, tested cross-disease associations for previously reported PS and CD risk alleles in additional 6,115 PS cases, 4,073 CD cases, and 10,100 controls. We identified seven susceptibility loci outside the human leukocyte antigen region (9p24 near JAK2, 10q22 at ZMIZ1, 11q13 near PRDX5, 16p13 near SOCS1, 17q21 at STAT3, 19p13 near FUT2, and 22q11 at YDJC) shared between PS and CD with genome-wide significance (p < 5 × 10⁻⁸) and confirmed four already established PS and CD risk loci (IL23R, IL12B, REL, and TYK2). Three of the shared loci are also genome-wide significantly associated with PS alone (10q22 at ZMIZ1, p_rs1250544 = 3.53 × 10⁻⁸, 11q13 near PRDX5, p_rs694739 = 3.71 × 10⁻⁰⁹, 22q11 at YDJC, p_rs181359 = 8.02 × 10⁻¹⁰). In addition, we identified one susceptibility locus for CD (16p13 near SOCS1, p_rs4780355 = 4.99 × 10⁻⁸). Refinement of association signals identified shared genome-wide significant associations for exonic SNPs at 10q22 (ZMIZ1) and in silico expression quantitative trait locus analyses revealed that the associations at ZMIZ1 and near SOCS1 have a potential functional effect on gene expression. Our results show the usefulness of joint analyses of clinically distinct immune-mediated diseases and enlarge the map of shared genetic risk loci.     

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.     

Association of JAG1 with Bone Mineral Density and Osteoporotic Fractures: A Genome-wide Association Study and Follow-up Replication Studies

Bone mineral density (BMD), a diagnostic parameter for osteoporosis and a clinical predictor of fracture, is a polygenic trait with high heritability. To identify genetic variants that influence BMD in different ethnic groups, we performed a genome-wide association study (GWAS) on 800 unrelated Southern Chinese women with extreme BMD and carried out follow-up replication studies in six independent study populations of European descent and Asian populations including 18,098 subjects. In the meta-analysis, rs2273061 of the Jagged1 (JAG1) gene was associated with high BMD (p = 5.27 × 10⁻⁸ for lumbar spine [LS] and p = 4.15 × 10⁻⁵ for femoral neck [FN], n = 18,898). This SNP was further found to be associated with the low risk of osteoporotic fracture (p = 0.009, OR = 0.7, 95% CI 0.57–0.93, n = 1881). Region-wide and haplotype analysis showed that the strongest association evidence was from the linkage disequilibrium block 5, which included rs2273061 of the JAG1 gene (p = 8.52 × 10⁻⁹ for LS and 3.47 × 10⁻⁵ at FN). To assess the function of identified variants, an electrophoretic mobility shift assay demonstrated the binding of c-Myc to the “G” but not “A” allele of rs2273061. A mRNA expression study in both human bone-derived cells and peripheral blood mononuclear cells confirmed association of the high BMD-related allele G of rs2273061 with higher JAG1 expression. Our results identify the JAG1 gene as a candidate for BMD regulation in different ethnic groups, and it is a potential key factor for fracture pathogenesis.     

Genetic Variation in the 3'-Untranslated Region of NBN Gene Is Associated with Gastric Cancer Risk in a Chinese Population

NBN plays a crucial role in carcinogenesis as a core component for both homologous recombination (HR) and non-homologous end-joining (NHEJ) DNA double-strand breaks (DSBs) repair pathways. Genetic variants in the NBN gene have been associated with multiple cancers risk, suggesting pleiotropic effect on cancer. We hypothesized that genetic variants in the NBN gene may modify the risk of gastric cancer. To test this hypothesis, we evaluated the association between four potentially functional single nucleotide polymorphisms in NBN and gastric cancer risk in a case–control study of 1,140 gastric cancer cases and 1,547 controls in a Chinese population. We found that the A allele of rs10464867 (G>A) was significantly associated with a decreased risk of gastric cancer (odds ratio [OR] = 0.81, 95% confidence interval [95% CI] = 0.71–0.94; P = 4.71×10⁻³). Furthermore, the association between A allele of rs10464867 and decreased risk of gastric cancer was more significantly in elder individuals (per-allele OR = 0.72[0.59–0.88], P = 1.07×10⁻³), and male individuals (per-allele OR = 0.73[0.62–0.87], P = 3.68×10⁻⁴). We further conducted a haplotype analysis and identified that the NBN A_rs10464867G_rs14448G_rs1063053 haplotype conferred stronger protective effect on gastric cancer (OR = 0.76[0.65–0.89], P = 6.39×10⁻⁴). In summary, these findings indicate that genetic variants at NBN gene may contribute to gastric cancer susceptibility and may further advance our understanding of NBN gene in cancer development.     

Abdominal aortic aneurysm is associated with a variant in low-density lipoprotein receptor-related protein 1

Abdominal aortic aneurysm (AAA) is a common cause of morbidity and mortality and has a significant heritability. We carried out a genome-wide association discovery study of 1866 patients with AAA and 5435 controls and replication of promising signals (lead SNP with a p value < 1 × 10⁻⁵) in 2871 additional cases and 32,687 controls and performed further follow-up in 1491 AAA and 11,060 controls. In the discovery study, nine loci demonstrated association with AAA (p < 1 × 10⁻⁵). In the replication sample, the lead SNP at one of these loci, rs1466535, located within intron 1 of low-density-lipoprotein receptor-related protein 1 (LRP1) demonstrated significant association (p = 0.0042). We confirmed the association of rs1466535 and AAA in our follow-up study (p = 0.035). In a combined analysis (6228 AAA and 49182 controls), rs1466535 had a consistent effect size and direction in all sample sets (combined p = 4.52 × 10⁻¹⁰, odds ratio 1.15 [1.10–1.21]). No associations were seen for either rs1466535 or the 12q13.3 locus in independent association studies of coronary artery disease, blood pressure, diabetes, or hyperlipidaemia, suggesting that this locus is specific to AAA. Gene-expression studies demonstrated a trend toward increased LRP1 expression for the rs1466535 CC genotype in arterial tissues; there was a significant (p = 0.029) 1.19-fold (1.04–1.36) increase in LRP1 expression in CC homozygotes compared to TT homozygotes in aortic adventitia. Functional studies demonstrated that rs1466535 might alter a SREBP-1 binding site and influence enhancer activity at the locus. In conclusion, this study has identified a biologically plausible genetic variant associated specifically with AAA, and we suggest that this variant has a possible functional role in LRP1 expression.     

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.     

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⁻⁸).     

Genetic variants at 13q12.12 are associated with high myopia in the Han Chinese population

High myopia, which is extremely prevalent in the Chinese population, is one of the leading causes of blindness in the world. Genetic factors play a critical role in the development of the condition. To identify the genetic variants associated with high myopia in the Han Chinese, we conducted a genome-wide association study (GWAS) of 493,947 SNPs in 1088 individuals (419 cases and 669 controls) from a Han Chinese cohort and followed up on signals that were associated with p < 1.0 × 10⁻⁴ in three independent cohorts (combined, 2803 cases and 5642 controls). We identified a significant association between high myopia and a variant at 13q12.12 (rs9318086, combined p = 1.91 × 10⁻¹⁶, heterozygous odds ratio = 1.32, and homozygous odds ratio = 1.64). Furthermore, five additional SNPs (rs9510902, rs3794338, rs1886970, rs7325450, and rs7331047) in the same linkage disequilibrium (LD) block with rs9318086 also proved to be significantly associated with high myopia in the Han Chinese population; p values ranged from 5.46 × 10⁻¹¹ to 6.16 × 10⁻¹⁶. This associated locus contains three genes—MIPEP, C1QTNF9B-AS1, and C1QTNF9B. MIPEP and C1QTNF9B were found to be expressed in the retina and retinal pigment epithelium (RPE) and are more likely than C1QTNF9B-AS1 to be associated with high myopia given the evidence of retinal signaling that controls eye growth. Our results suggest that the variants at 13q12.12 are associated with high myopia.     

Genome-wide scan identifies TNIP1, PSORS1C1, and RHOB as novel risk loci for systemic sclerosis

Systemic sclerosis (SSc) is an orphan, complex, inflammatory disease affecting the immune system and connective tissue. SSc stands out as a severely incapacitating and life-threatening inflammatory rheumatic disease, with a largely unknown pathogenesis. We have designed a two-stage genome-wide association study of SSc using case-control samples from France, Italy, Germany, and Northern Europe. The initial genome-wide scan was conducted in a French post quality-control sample of 564 cases and 1,776 controls, using almost 500 K SNPs. Two SNPs from the MHC region, together with the 6 loci outside MHC having at least one SNP with a P<10⁻⁵ were selected for follow-up analysis. These markers were genotyped in a post-QC replication sample of 1,682 SSc cases and 3,926 controls. The three top SNPs are in strong linkage disequilibrium and located on 6p21, in the HLA-DQB1 gene: rs9275224, P = 9.18×10⁻⁸, OR = 0.69, 95% CI [0.60–0.79]; rs6457617, P = 1.14×10⁻⁷ and rs9275245, P = 1.39×10⁻⁷. Within the MHC region, the next most associated SNP (rs3130573, P = 1.86×10⁻⁵, OR = 1.36 [1.18–1.56]) is located in the PSORS1C1 gene. Outside the MHC region, our GWAS analysis revealed 7 top SNPs (P<10⁻⁵) that spanned 6 independent genomic regions. Follow-up of the 17 top SNPs in an independent sample of 1,682 SSc and 3,926 controls showed associations at PSORS1C1 (overall P = 5.70×10⁻¹⁰, OR:1.25), TNIP1 (P = 4.68×10⁻⁹, OR:1.31), and RHOB loci (P = 3.17×10⁻⁶, OR:1.21). Because of its biological relevance, and previous reports of genetic association at this locus with connective tissue disorders, we investigated TNIP1 expression. A markedly reduced expression of the TNIP1 gene and also its protein product were observed both in lesional skin tissue and in cultured dermal fibroblasts from SSc patients. Furthermore, TNIP1 showed in vitro inhibitory effects on inflammatory cytokine-induced collagen production. The genetic signal of association with TNIP1 variants, together with tissular and cellular investigations, suggests that this pathway has a critical role in regulating autoimmunity and SSc pathogenesis.     

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

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

Caucasian and Asian specific rheumatoid arthritis risk loci reveal limited replication and apparent allelic heterogeneity in north Indians

Genome-wide association studies and meta-analysis indicate that several genes/loci are consistently associated with rheumatoid arthritis (RA) in European and Asian populations. To evaluate the transferability status of these findings to an ethnically diverse north Indian population, we performed a replication analysis. We investigated the association of 47 single-nucleotide polymorphisms (SNPs) at 43 of these genes/loci with RA in a north Indian cohort comprising 983 RA cases and 1007 age and gender matched controls. Genotyping was done using Infinium human 660w-quad. Association analysis by chi-square test implemented in plink was carried out in two steps. Firstly, association of the index or surrogate SNP (r2>0.8, calculated from reference GIH Hap-Map population) was tested. In the second step, evidence for allelic/locus heterogeneity at aforementioned genes/loci was assessed for by testing additional flanking SNPs in linkage equilibrium with index/surrogate marker.Of the 44 European specific index SNPs, neither index nor surrogate SNPs were present for nine SNPs in the genotyping array. Of the remaining 35, associations were replicated at seven genes namely PTPN22 (rs1217407, p = 3×10⁻³); IL2–21 (rs13119723, p = 0.008); HLA-DRB1 (rs660895, p = 2.56×10⁻⁵; rs6457617, p = 1.6×10⁻⁰⁹; rs13192471, p = 6.7×10⁻¹⁶); TNFA1P3 (rs9321637, p = 0.03); CCL21 (rs13293020, p = 0.01); IL2RA (rs2104286, p = 1.9×10⁻⁴) and ZEB1 (rs2793108, p = 0.006). Of the three Asian specific loci tested, rs2977227 in PADI4 showed modest association (p<0.02). Further, of the 140 SNPs (in LE with index/surrogate variant) tested, association was observed at 11 additional genes: PTPRC, AFF3, CD28, CTLA4, PXK, ANKRD55, TAGAP, CCR6, BLK, CD40 and IL2RB. This study indicates limited replication of European and Asian index SNPs and apparent allelic heterogeneity in RA etiology among north Indians warranting independent GWAS in this population. However, replicated associations of HLA-DRB1, PTPN22 (which confer ∼50% of the heritable risk to RA) and IL2RA suggest that cross-ethnicity fine mapping of such loci is apposite for identification of causal variants.