Lack of replication of higher genetic risk load in men than in women with systemic lupus erythematosus.

INTRODUCTION: We aimed to replicate a recent study which showed higher genetic risk load at 15 loci in men than in women with systemic lupus erythematosus (SLE). This difference was very significant, and it was interpreted as indicating that men require more genetic susceptibility than women to develop SLE. METHODS: Nineteen SLE-associated loci (thirteen of which are shared with the previous study) were analyzed in 1,457 SLE patients and 1,728 healthy controls of European ancestry. Genetic risk load was calculated as sex-specific sum genetic risk scores (GRS(s)). RESULTS: Our results did not replicate those of the previous study at either the level of individual loci or the global level of GRS(s). GRS(s) were larger in women than in men (4.20 ± 1.07 in women vs. 3.27 ± 0.98 in men). This very significant difference (P < 10(-16)) was more dependent on the six new loci not included in the previous study (59% of the difference) than on the thirteen loci that are shared (the remaining 41%). However, the 13 shared loci also showed a higher genetic risk load in women than in men in our study (P = 6.6 × 10(-7)), suggesting that heterogeneity of participants, in addition to different loci, contributed to the opposite results. CONCLUSION: Our results show the lack of a clear trend toward higher genetic risk in one of the sexes for the analyzed SLE loci. They also highlight several limitations of assessments of genetic risk load, including the possibility of ascertainment bias with loci discovered in studies that have included mainly women.

Further evidence of subphenotype association with systemic lupus erythematosus susceptibility loci: a European cases only study.

INTRODUCTION: Systemic Lupus Erythematosus (SLE) shows a spectrum of clinical manifestations that complicate its diagnosis, treatment and research. This variability is likely related with environmental exposures and genetic factors among which known SLE susceptibility loci are prime candidates. The first published analyses seem to indicate that this is the case for some of them, but results are still inconclusive and we aimed to further explore this question. METHODS: European SLE patients, 1444, recruited at 17 centres from 10 countries were analyzed. Genotypes for 26 SLE associated SNPs were compared between patients with and without each of 11 clinical features: ten of the American College of Rheumatology (ACR) classification criteria (except ANAs) and age of disease onset. These analyses were adjusted for centre of recruitment, top ancestry informative markers, gender and time of follow-up. Overlap of samples with previous studies was excluded for assessing replication. RESULTS: THERE WERE THREE NEW ASSOCIATIONS: the SNPs in XKR6 and in FAM167A-BLK were associated with lupus nephritis (OR=0.76 and 1.30, P(corr) =0.007 and 0.03, respectively) and the SNP of MECP2, which is in chromosome X, with earlier age of disease onset in men. The previously reported association of STAT4 with early age of disease onset was replicated. Some other results were suggestive of the presence of additional associations. Together, the association signals provided support to some previous findings and to the characterization of lupus nephritis, autoantibodies and age of disease onset as the clinical features more associated with SLE loci. CONCLUSION: Some of the SLE loci shape the disease phenotype in addition to increase susceptibility to SLE. This influence is more prominent for some clinical features than for others. However, results are only partially consistent between studies and subphenotype specific GWAS are needed to unravel their genetic component.

Bias in effect size of systemic lupus erythematosus susceptibility loci across Europe: a case-control study.

INTRODUCTION: We aimed to investigate whether the effect size of the systemic lupus erythematosus (SLE) risk alleles varies across European subpopulations. METHODS: European SLE patients (n = 1,742) and ethnically matched healthy controls (n = 2,101) were recruited at 17 centres from 10 different countries. Only individuals with self-reported ancestry from the country of origin were included. In addition, participants were genotyped for top ancestry informative markers and for 25 SLE associated SNPs. The results were used to compare effect sizes between the Central Eureopan and Southern European subgroups. RESULTS: Twenty of the 25 SNPs showed independent association with SLE, These SNPs showed a significant bias to larger effect sizes in the Southern subgroup, with 15/20 showing this trend (P = 0.019) and a larger mean odds ratio of the 20 SNPs (1.46 vs. 1.34, P = 0.02) as well as a larger difference in the number of risk alleles (2.06 vs. 1.63, P = 0.027) between SLE patients and controls than for Central Europeans. This bias was reflected in a very significant difference in the cumulative genetic risk score (4.31 vs. 3.48, P = 1.8 × 10-32). Effect size bias was accompanied by a lower number of SLE risk alleles in the Southern subjects, both patients and controls, the difference being more marked between the controls (P = 1.1 × 10-8) than between the Southern and Central European patients (P = 0.016). Seven of these SNPs showed significant allele frequency clines. CONCLUSION: Our findings showed a bias to larger effect sizes of SLE loci in the Southern Europeans relative to the Central Europeans together with clines of SLE risk allele frequencies. These results indicate the need to study risk allele clines and the implications of the polygenic model of inheritance in SLE.

Association of systemic lupus erythematosus clinical features with European population genetic substructure.

Systemic Lupus Erythematosus (SLE) is an autoimmune disease with a very varied spectrum of clinical manifestations that could be partly determined by genetic factors. We aimed to determine the relationship between prevalence of 11 clinical features and age of disease onset with European population genetic substructure. Data from 1413 patients of European ancestry recruited in nine countries was tested for association with genotypes of top ancestry informative markers. This analysis was done with logistic regression between phenotypes and genotypes or principal components extracted from them. We used a genetic additive model and adjusted for gender and disease duration. Three clinical features showed association with ancestry informative markers: autoantibody production defined as immunologic disorder (P = 6.8×10(-4)), oral ulcers (P = 6.9×10(-4)) and photosensitivity (P = 0.002). Immunologic disorder was associated with genotypes more common in Southern European ancestries, whereas the opposite trend was observed for photosensitivity. Oral ulcers were specifically more common in patients of Spanish and Portuguese self-reported ancestry. These results should be taken into account in future research and suggest new hypotheses and possible underlying mechanisms to be investigated. A first hypothesis linking photosensitivity with variation in skin pigmentation is suggested.

Cis-regulation of IRF5 expression is unable to fully account for systemic lupus erythematosus association: analysis of multiple experiments with lymphoblastoid cell lines.

INTRODUCTION: Interferon regulatory factor 5 gene (IRF5) polymorphisms are strongly associated with several diseases, including systemic lupus erythematosus (SLE). The association includes risk and protective components. They could be due to combinations of functional polymorphisms and related to cis-regulation of IRF5 expression, but their mechanisms are still uncertain. We hypothesised that thorough testing of the relationships between IRF5 polymorphisms, expression data from multiple experiments and SLE-associated haplotypes might provide useful new information. METHODS: Expression data from four published microarray hybridisation experiments with lymphoblastoid cell lines (57 to 181 cell lines) were retrieved. Genotypes of 109 IRF5 polymorphisms, including four known functional polymorphisms, were considered. The best linear regression models accounting for the IRF5 expression data were selected by using a forward entry procedure. SLE-associated IRF5 haplotypes were correlated with the expression data and with the best cis-regulatory models. RESULTS: A large fraction of variability in IRF5 expression was accounted for by linear regression models with IRF5 polymorphisms, but at a different level in each expression data set. Also, the best models from each expression data set were different, although there was overlap between them. The SNP introducing an early polyadenylation signal, rs10954213, was included in the best models for two of the expression data sets and in good models for the other two data sets. The SLE risk haplotype was associated with high IRF5 expression in the four expression data sets. However, there was also a trend towards high IRF5 expression with some protective and neutral haplotypes, and the protective haplotypes were not associated with IRF5 expression. As a consequence, correlation between the cis-regulatory best models and SLE-associated haplotypes, regarding either the risk or protective component, was poor. CONCLUSIONS: Our analysis indicates that although the SLE risk haplotype of IRF5 is associated with high expression of the gene, cis-regulation of IRF5 expression is not enough to fully account for IRF5 association with SLE susceptibility, which indicates the need to identify additional functional changes in this gene.

Lack of replication of genetic predictors for the rheumatoid arthritis response to anti-TNF treatments: a prospective case-only study.

INTRODUCTION: We aimed to replicate the strong associations that a recent genome wide association study (GWAS) has found between 16 single nucleotide polymorphisms (SNPs) and response to anti-tumour necrosis factor (TNF) treatment in 89 patients with rheumatoid arthritis (RA). This study is very important because, according to published simulations, associations as strong as the reported ones will mean that these SNPs could be used as predictors of response at the individual level. METHODS: Disease activity score (DAS28) was evaluated in 151 anti-TNF treated patients with RA of Spanish ancestry at baseline and every 3 months thereafter. Genotypes of the 16 putative predictor SNPs were obtained by single-base extension. Association between the relative change in DAS28 and SNP genotypes was tested by linear regression. In addition, logistic regression was applied to compare genotypes in non-responders (n = 34) versus good-responders (n = 61) following the EULAR response criteria. RESULTS: None of the analyses showed any significant association between the 16 SNPs and response to anti-TNF treatments at 3 or 6 months. Results were also negative when only patients treated with infliximab (66.9% of the total) were separately analyzed. These negative results were obtained in spite of a very good statistical power to replicate the reported strong associations. CONCLUSIONS: We still do not have any sound evidence of genetic variants associated with RA response to anti-TNF treatments. In addition, the possibility we had envisaged of using the results of a recent GWAS for prediction in individual patients should be dismissed.

Rheumatoid arthritis does not share most of the newly identified systemic lupus erythematosus genetic factors.

OBJECTIVE: Rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) share some genetic factors such as HLA, PTPN22, STAT4, and 6q23. The aim of this study was to determine whether 9 other SLE genetic factors are also implicated in RA susceptibility. METHODS: A characteristic single-nucleotide polymorphism (SNP) in each of 9 genetic factors, ITGAM (rs1143679), C8orf13-BLK (rs13277113), TYK2 (rs2304256), 1q25.1 (rs10798269), PXK (rs6445975), KIAA1542 (rs4963128), MECP2 (rs17435), BANK1 (rs17266594), and LY9 (rs509749), was studied in 1,635 patients with RA and 1,906 control subjects from Spain. The rs7574865 SNP in STAT4 was also included. Analyses were conducted globally and after stratification by sex and clinical features (anti-cyclic citrullinated peptide and rheumatoid factor, shared epitope, rheumatoid nodules, radiographic changes, sicca syndrome, and pneumonitis). RESULTS: No association was observed between RA and any of the 9 newly identified SLE genetic factors. A meta-analysis using previous data was consistent with these results. In addition, there were no significant differences between individuals with and those without each of the clinical features analyzed, except the frequency of the minor allele in the C8orf13-BLK locus that was decreased in patients with sicca syndrome (14.6% versus 22.4% in controls; P = 0.003). CONCLUSION: None of the 9 recently identified SLE risk factors showed association with RA. Therefore, common genetic factors affecting the pathogenesis of these 2 disorders seem to be limited, revealing that the genetic component contributes to the different expression of these diseases.

Replication of recently identified systemic lupus erythematosus genetic associations: a case-control study.

INTRODUCTION: We aimed to replicate association of newly identified systemic lupus erythematosus (SLE) loci. METHODS: We selected the most associated SNP in 10 SLE loci. These 10 SNPs were analysed in 1,579 patients with SLE and 1,726 controls of European origin by single-base extension. Comparison of allele frequencies between cases and controls was done with the Mantel-Haenszel approach to account for heterogeneity between sample collections. RESULTS: A previously controversial association with a SNP in the TYK2 gene was replicated (odds ratio (OR) = 0.79, P = 2.5 x 10-5), as well as association with the X chromosome MECP2 gene (OR = 1.26, P = 0.00085 in women), which had only been reported in a single study, and association with four other loci, 1q25.1 (OR = 0.81, P = 0.0001), PXK (OR = 1.19, P = 0.0038), BANK1 (OR = 0.83, P = 0.006) and KIAA1542 (OR = 0.84, P = 0.001), which have been identified in a genome-wide association study, but not found in any other study. All these replications showed the same disease-associated allele as originally reported. No association was found with the LY9 SNP, which had been reported in a single study. CONCLUSIONS: Our results confirm nine SLE loci. For six of them, TYK2, MECP2, 1q25.1, PXK, BANK1 and KIAA1542, this replication is important. The other three loci, ITGAM, STAT4 and C8orf13-BLK, were already clearly confirmed. Our results also suggest that MECP2 association has no influence in the sex bias of SLE, contrary to what has been proposed. In addition, none of the other associations seems important in this respect.