Stratification of Systemic Lupus Erythematosus Patients Into Three Groups of Disease Activity Progression According to Longitudinal Gene Expression

OBJECTIVE

The highly heterogeneous clinical presentation of systemic lupus erythematosus (SLE) is characterized by the unpredictable occurrence of disease flares and organ damage. Attempts to stratify lupus patients have been limited to classification based on clinical information, leading to unsuccessful clinical trials and controversial research results. This study was undertaken to develop and validate a robust method to stratify patients with lupus according to longitudinal disease activity and whole-genome gene expression data in order to establish subgroups of patients who share disease progression mechanisms.

METHODS

We used a cluster-based approach to stratify SLE patients based on the correlation between disease activity scores and longitudinal gene expression information. Clustering robustness was evaluated by the bootstrap method, and the clusters were characterized in terms of clinical and functional features.

RESULTS

We observed a clear partition into 3 different disease clusters in 2 independent sets of patients, one pediatric and one adult, which was not influenced by treatment, race, or other source of bias. Two of the clusters differentiated into a group showing a correlation between the percentage of neutrophils and disease activity progression and a group showing a correlation between the percentage of lymphocytes and disease activity progression. The third cluster, in which the percentage of neutrophils correlated to a lesser degree with disease activity, was functionally more heterogeneous. Patients in the neutrophil-driven clusters had an increased risk of developing proliferative nephritis.

CONCLUSION

Our findings indicate that SLE patients can be stratified into 3 subgroups of patients who show different mechanisms of disease progression and are clinically differentiated. Our results have important implications for treatment options, the design of clinical trials, our understanding of the etiology of the disease, and the prediction of severe glomerulonephritis.

Trans-Ethnic Mapping of BANK1 Identifies Two Independent SLE-Risk Linkage Groups Enriched for Co-Transcriptional Splicing Marks

BANK1 is a susceptibility gene for several systemic autoimmune diseases in several populations. Using the genome-wide association study (GWAS) data from Europeans (EUR) and African Americans (AA), we performed an extensive fine mapping of ankyrin repeats 1 (BANK1). To increase the SNP density, we used imputation followed by univariate and conditional analysis, combined with a haplotypic and expression quantitative trait locus (eQTL) analysis. The data from Europeans showed that the associated region was restricted to a minimal and dependent set of SNPs covering introns two and three, and exon two. In AA, the signal found in the Europeans was split into two independent effects. All of the major risk associated SNPs were eQTLs, and the risks were associated with an increased BANK1 gene expression. Functional annotation analysis revealed the enrichment of repressive B cell epigenomic marks (EZH2 and H3K27me3) and a strong enrichment of splice junctions. Furthermore, one eQTL located in intron two, rs13106926, was found within the binding site for RUNX3, a transcriptional activator. These results connect the local genome topography, chromatin structure, and the regulatory landscape of BANK1 with co-transcriptional splicing of exon two. Our data defines a minimal set of risk associated eQTLs predicted to be involved in the expression of BANK1 modulated through epigenetic regulation and splicing. These findings allow us to suggest that the increased expression of BANK1 will have an impact on B-cell mediated disease pathways.

Metagene projection characterizes GEN2.2 and CAL-1 as relevant human plasmacytoid dendritic cell models

Motivation

Plasmacytoid dendritic cells (pDC) play a major role in the regulation of adaptive and innate immunity. Human pDC are difficult to isolate from peripheral blood and do not survive in culture making the study of their biology challenging. Recently, two leukemic counterparts of pDC, CAL-1 and GEN2.2, have been proposed as representative models of human pDC. Nevertheless, their relationship with pDC has been established only by means of particular functional and phenotypic similarities. With the aim of characterizing GEN2.2 and CAL-1 in the context of the main circulating immune cell populations we have performed microarray gene expression profiling of GEN2.2 and carried out an integrated analysis using publicly available gene expression datasets of CAL-1 and the main circulating primary leukocyte lineages.

Results

Our results show that GEN2.2 and CAL-1 share common gene expression programs with primary pDC, clustering apart from the rest of circulating hematopoietic lineages. We have also identified common differentially expressed genes that can be relevant in pDC biology. In addition, we have revealed the common and differential pathways activated in primary pDC and cell lines upon CpG stimulatio.

Availability and implementation

R code and data are available in the supplementary material.

Contact

pedro.carmona@genyo.es or concepcion.maranon@genyo.es.

Supplementary information

Supplementary data are available at Bioinformatics online.

A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.

Whole Exome Sequencing of Patients from Multicase Families with Systemic Lupus Erythematosus Identifies Multiple Rare Variants

In an effort to identify rare alleles associated with SLE, we have performed whole exome sequencing of the most distantly related affected individuals from two large Icelandic multicase SLE families followed by Ta targeted genotyping of additional relatives. We identified multiple rare likely pathogenic variants in nineteen genes co-segregating with the disease through multiple generations. Gene co-expression and protein-protein interaction analysis identified a network of highly connected genes comprising several loci previously implicated in autoimmune diseases. These genes were significantly enriched for immune system development, lymphocyte activation, DNA repair, and V(D)J gene recombination GO-categories. Furthermore, we found evidence of aggregate association and enrichment of rare variants at the FAM71E1/EMC10 locus in an independent set of 4,254 European SLE-cases and 4,349 controls. Our study presents evidence supporting that multiple rare likely pathogenic variants, in newly identified genes involved in known disease pathogenic pathways, segregate with SLE at the familial and population level.

Epigenome-wide association studies for systemic autoimmune diseases: The road behind and the road ahead

Epigenetics is known to be an important mechanism in the pathogenesis of autoimmune diseases. Epigenetic variations can act as integrators of environmental and genetic exposures and propagate activated states in immune cells. Studying epigenetic alterations by means of genome-wide approaches promises to unravel novel molecular mechanisms related to disease etiology, disease progression, clinical manifestations and treatment responses. This paper reviews what we have learned in the last five years from epigenome-wide studies for three systemic autoimmune diseases, namely systemic lupus erythematosus, primary Sjögren's syndrome, and rheumatoid arthritis. We examine the degree of epigenetic sharing between different diseases and the possible mediating role of epigenetic associations in genetic and environmental risks. Finally, we also shed light into the use of epigenetic markers towards a better precision medicine regarding disease prediction, prevention and personalized treatment in systemic autoimmunity.

Effects of Amerindian Genetic Ancestry on Clinical Variables and Therapy in Patients with Rheumatoid Arthritis

OBJECTIVE

To define whether Amerindian genetic ancestry correlates with clinical and therapeutic variables in admixed individuals with rheumatoid arthritis (RA) from Latin America.

METHODS

Patients with RA (n = 1347) and healthy controls (n = 1012) from Argentina, Mexico, Chile, and Peru were included. Samples were genotyped for the Immunochip v1 using the Illumina platform. Clinical data were obtained through interviews or the clinical history.

RESULTS

Percentage of Amerindian ancestry was comparable between cases and controls. Morning stiffness (p < 0.0001, OR 0.05), rheumatoid factor (RF; p < 0.0001, OR 0.22), radiographic changes (p < 0.0001, OR 0.05), and higher number of criteria were associated with lower Amerindian ancestry after Bonferroni correction. Higher Amerindian ancestry correlated only with weight loss (p_Bonferroni < 0.0001, OR 2.85). Increased Amerindian ancestry correlated with higher doses of azathioprine (p < 0.0001, OR 163.6) and sulfasalazine (p < 0.0001, OR 48.6), and inversely with methotrexate (p = 0.001, OR 0.35), leflunomide (p = 0.001, OR 0.16), and nonsteroidal antiinflammatory drugs (p_Bonferroni = 0.001, OR 0.37). Only the presence of RF and weight loss were modified after confounders adjustment.

CONCLUSION

Amerindian ancestry protects against most major clinical criteria of RA, but regarding the association of RF with increased European ancestry, age, sex, and smoking are modifiers. Ancestry also correlates with the therapeutic profiles.

SLE redefined on the basis of molecular pathways

The implementation of precision medicine requires the recruiting of patients in statistically enough numbers, the possibility of obtaining enough materials, and the integration of data from various platforms, which are all real limitations. These types of studies have been performed extensively in cancer but barely on systemic lupus erythematosus (SLE) or other rheumatic diseases. To consider the practical use of the information obtained from such studies, we have to take into account the best biological fluid to use, the ease to perform the analysis in clinical practice, and its relevance to clinical practice. Here we review the most relevant studies that have performed analyses that attempt to classify or stratify SLE. We focus on two types of studies: those that stratify individuals diagnosed with SLE and those that compare SLE with other autoimmune diseases, defining differences and similarities that may be clinically relevant in the future.

Omics studies: their use in diagnosis and reclassification of SLE and other systemic autoimmune diseases

Omics studies of systemic autoimmune diseases (SADs) in general, and SLE in particular, have delivered isolated information from transcriptome, epigenome, genome, cytokine and metabolome analyses. Such analyses have resulted in the identification of disease susceptibility genes and the description of IFN expression signatures, allowing extensive insight into the mechanisms of disease and the development of new therapies. Access to such technologies allows the recognition of patterns of disease at a pathway level, thereby, to reclassify SLE and other SADs and to develop new therapeutics from a personalized perspective. The use of omic information allows the discovery of correlative patterns involving drugs not currently suspected to be of value in SADs. In this review, we summarize the omics findings for SLE and propose ways of using the data for the identification of new biomarkers, finding new drugs and reclassifying patients not only with SLE, but also with other SADs.

Genetic architecture distinguishes systemic juvenile idiopathic arthritis from other forms of juvenile idiopathic arthritis: clinical and therapeutic implications

OBJECTIVES

Juvenile idiopathic arthritis (JIA) is a heterogeneous group of conditions unified by the presence of chronic childhood arthritis without an identifiable cause. Systemic JIA (sJIA) is a rare form of JIA characterised by systemic inflammation. sJIA is distinguished from other forms of JIA by unique clinical features and treatment responses that are similar to autoinflammatory diseases. However, approximately half of children with sJIA develop destructive, long-standing arthritis that appears similar to other forms of JIA. Using genomic approaches, we sought to gain novel insights into the pathophysiology of sJIA and its relationship with other forms of JIA.

METHODS

We performed a genome-wide association study of 770 children with sJIA collected in nine countries by the International Childhood Arthritis Genetics Consortium. Single nucleotide polymorphisms were tested for association with sJIA. Weighted genetic risk scores were used to compare the genetic architecture of sJIA with other JIA subtypes.

RESULTS

The major histocompatibility complex locus and a locus on chromosome 1 each showed association with sJIA exceeding the threshold for genome-wide significance, while 23 other novel loci were suggestive of association with sJIA. Using a combination of genetic and statistical approaches, we found no evidence of shared genetic architecture between sJIA and other common JIA subtypes.

CONCLUSIONS

The lack of shared genetic risk factors between sJIA and other JIA subtypes supports the hypothesis that sJIA is a unique disease process and argues for a different classification framework. Research to improve sJIA therapy should target its unique genetics and specific pathophysiological pathways.

Support for phosphoinositol 3 kinase and mTOR inhibitors as treatment for lupus using in-silico drug-repurposing analysis

Background

Systemic lupus erythematosus (SLE) is an autoimmune disease with few treatment options. Current therapies are not fully effective and show highly variable responses. In this regard, large efforts have focused on developing more effective therapeutic strategies. Drug repurposing based on the comparison of gene expression signatures is an effective technique for the identification of new therapeutic approaches. Here we present a drug-repurposing exploratory analysis using gene expression signatures from SLE patients to discover potential new drug candidates and target genes.

Methods

We collected a compendium of gene expression signatures comprising peripheral blood cells and different separate blood cell types from SLE patients. The Lincscloud database was mined to link SLE signatures with drugs, gene knock-down, and knock-in expression signatures. The derived dataset was analyzed in order to identify compounds, genes, and pathways that were significantly correlated with SLE gene expression signatures.

Results

We obtained a list of drugs that showed an inverse correlation with SLE gene expression signatures as well as a set of potential target genes and their associated biological pathways. The list includes drugs never or little studied in the context of SLE treatment, as well as recently studied compounds.

Conclusion

Our exploratory analysis provides evidence that phosphoinositol 3 kinase and mammalian target of rapamycin (mTOR) inhibitors could be potential therapeutic options in SLE worth further future testing.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-017-1263-7) contains supplementary material, which is available to authorized users.

CymeR: cytometry analysis using KNIME, docker and R

Summary

Here we present open-source software for the analysis of high-dimensional cytometry data using state of the art algorithms. Importantly, use of the software requires no programming ability, and output files can either be interrogated directly in CymeR or they can be used downstream with any other cytometric data analysis platform. Also, because we use Docker to integrate the multitude of components that form the basis of CymeR, we have additionally developed a proof-of-concept of how future open-source bioinformatic programs with graphical user interfaces could be developed.

Availability and Implementation

CymeR is open-source software that ties several components into a single program that is perhaps best thought of as a self-contained data analysis operating system. Please see https://github.com/bmuchmore/CymeR/wiki for detailed installation instructions.

Contact

brian.muchmore@genyo.es or marta.alarcon@genyo.es.

Tartrate-Resistant Acid Phosphatase Deficiency in the Predisposition to Systemic Lupus Erythematosus

OBJECTIVE

Mutations in the ACP5 gene, which encodes tartrate-resistant acid phosphatase (TRAP), cause the immuno-osseous disorder spondyloenchondrodysplasia, which includes as disease features systemic lupus erythematosus (SLE) and a type I interferon (IFN) signature. Our aims were to identify TRAP substrates, determine the consequences of TRAP deficiency in immune cells, and assess whether ACP5 mutations are enriched in sporadic cases of SLE.

METHODS

Interaction between TRAP and its binding partners was tested by a yeast 2-hybrid screening, confocal microscopy, and immunoprecipitation/Western blotting. TRAP knockdown was performed using small interfering RNA. Phosphorylation of osteopontin (OPN) was analyzed by mass spectrometry. Nucleotide sequence analysis of ACP5 was performed by Sanger sequencing or next-generation sequencing.

RESULTS

TRAP and OPN colocalized and interacted in human macrophages and plasmacytoid dendritic cells (PDCs). TRAP dephosphorylated 3 serine residues on specific OPN peptides. TRAP knockdown resulted in increased OPN phosphorylation and increased nuclear translocation of IRF7 and P65, with resultant heightened expression of IFN-stimulated genes and IL6 and TNF following Toll-like receptor 9 stimulation. An excess of heterozygous ACP5 missense variants was observed in SLE compared to controls (P = 0.04), and transfection experiments revealed a significant reduction in TRAP activity in a number of variants.

CONCLUSION

Our findings indicate that TRAP and OPN colocalize and that OPN is a substrate for TRAP in human immune cells. TRAP deficiency in PDCs leads to increased IFNα production, providing at least a partial explanation for how ACP5 mutations cause lupus in the context of spondyloenchondrodysplasia. Detection of ACP5 missense variants in a lupus cohort suggests that impaired TRAP functioning may increase susceptibility to sporadic lupus.

Genome-Wide Association Study in an Amerindian Ancestry Population Reveals Novel Systemic Lupus Erythematosus Risk Loci and the Role of European Admixture

OBJECTIVE

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease with a strong genetic component. We undertook the present work to perform the first genome-wide association study on individuals from the Americas who are enriched for Native American heritage.

METHODS

We analyzed 3,710 individuals from the US and 4 countries of Latin America who were diagnosed as having SLE, and healthy controls. Samples were genotyped with HumanOmni1 BeadChip. Data on out-of-study controls genotyped with HumanOmni2.5 were also included. Statistical analyses were performed using SNPtest and SNPGWA. Data were adjusted for genomic control and false discovery rate. Imputation was performed using Impute2 and, for classic HLA alleles, HiBag. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated.

RESULTS

The IRF5-TNPO3 region showed the strongest association and largest OR for SLE (rs10488631: genomic control-adjusted P [Pgcadj ] = 2.61 × 10(-29), OR 2.12 [95% CI 1.88-2.39]), followed by HLA class II on the DQA2-DQB1 loci (rs9275572: Pgcadj  = 1.11 × 10(-16), OR 1.62 [95% CI 1.46-1.80] and rs9271366: Pgcadj  = 6.46 × 10(-12), OR 2.06 [95% CI 1.71-2.50]). Other known SLE loci found to be associated in this population were ITGAM, STAT4, TNIP1, NCF2, and IRAK1. We identified a novel locus on 10q24.33 (rs4917385: Pgcadj  = 1.39 × 10(-8)) with an expression quantitative trait locus (eQTL) effect (Peqtl  = 8.0 × 10(-37) at USMG5/miR1307), and several new suggestive loci. SLE risk loci previously identified in Europeans and Asians were corroborated. Local ancestry estimation showed that the HLA allele risk contribution is of European ancestral origin. Imputation of HLA alleles suggested that autochthonous Native American haplotypes provide protection against development of SLE.

CONCLUSION

Our results demonstrate that studying admixed populations provides new insights in the delineation of the genetic architecture that underlies autoimmune and complex diseases.

The genetic basis of systemic lupus erythematosus: What are the risk factors and what have we learned

The genome-wide association study is a free-hypothesis approach based on screening of thousands or even millions of genetic variants distributed throughout the whole human genome in relation to a phenotype. The relevant role of the genome-wide association studies in the last decade is undisputed because it has permitted to elucidate multiple risk genetic factors associated with the susceptibility to several human complex diseases. Regarding systemic lupus erythematosus (SLE) this approach has allowed to identify more than 60 risk loci for SLE susceptibility across populations to date, increasing our understanding on the pathogenesis of this disease. We present the latest findings in the genetic of SLE across populations using genome-wide approaches. These studies revealed that most of the genetic risk is shared across borders and ethnicities. Finally, we focus on describing the most important risk loci for SLE attempting to cover the genetic findings in relation to functional polymorphisms, such as missense single nucleotide polymorphisms (SNPs) or regulatory variants involved in the development of the disease. The functional studies try to identify the causality of some GWAS-associated variants, many of which fall in non-coding regions of the genome, suggesting a regulatory role. Many loci show an environmental interaction, another aspect revealed by the studies of epigenetic modifications and those associated with genetic variants. Finally, new-generation sequencing technologies can open other paths in the research on SLE genetics, the role of rare variants and the detailed identification of causal regulatory variation. The clinical relevance of the genetic factors will be shown when we are able to use them or in combination with other molecular measurements to re-classify a heterogeneous disease such as SLE.

BANK1 Regulates IgG Production in a Lupus Model by Controlling TLR7-Dependent STAT1 Activation

The purpose of our study was to investigate the effects of the adaptor Bank1 in TLR7 signaling using the B6.Sle1.yaa mouse, a lupus model that develops disease through exacerbated TLR7 expression. Crosses of B6.Sle1.yaa with Bank1^-/- mice maintained several B and myeloid cell phenotypes close to normal wild-type levels. Most striking was the reduction in total serum IgG antibodies, but not of IgM, and reduced serum levels of autoantibodies, IL-6, and BAFF. Bank1 deficiency did modify numbers of MZ B cells and total B cell numbers, as well as expression of CXCR4 by follicular helper T cells. Other T cell changes were not observed. Bank1 deficiency did not modify numbers of germinal center B cells or plasma cells or clinical disease outcomes. Purified B cells from Bank1 deficient mice had strongly reduced Ifnb, Ifna4, Irf7, Aicda and Stat1 gene expression following TLR7 agonist stimulation. Interestingly, phosphorylation of Tyr701, but not of Ser727 of STAT1, was impaired in splenic B cells from B6.Sle1.yaa.Bank1-/- mice, as was the nuclear translocation of IRF7 in response to TLR7 agonist stimulation. Further, Bank1 deficiency in B6.Sle1.yaa mice reduced the production of IgG2c after in vitro TLR7 agonist stimulation. Our results demonstrate that Bank1 controls TLR7-mediated type I interferon production. Combined with the control of the nuclear translocation of IRF7, the modulation of STAT1 transcription and phosphorylation, Bank1 contributes to IgG production during development of autoimmune disease.

A combined large-scale meta-analysis identifies COG6 as a novel shared risk locus for rheumatoid arthritis and systemic lupus erythematosus

OBJECTIVES

During the last years, genome-wide association studies (GWASs) have identified a number of common genetic risk factors for rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE). However, the genetic overlap between these two immune-mediated diseases has not been thoroughly examined so far. The aim of the present study was to identify additional risk loci shared between RA and SLE.

METHODS

We performed a large-scale meta-analysis of GWAS data from RA (3911 cases and 4083 controls) and SLE (2237 cases and 6315 controls). The top-associated polymorphisms in the discovery phase were selected for replication in additional datasets comprising 13 641 RA cases and 31 921 controls and 1957 patients with SLE and 4588 controls.

RESULTS

The rs9603612 genetic variant, located nearby the COG6 gene, an established susceptibility locus for RA, reached genome-wide significance in the combined analysis including both discovery and replication sets (p value=2.95E-13). In silico expression quantitative trait locus analysis revealed that the associated polymorphism acts as a regulatory variant influencing COG6 expression. Moreover, protein-protein interaction and gene ontology enrichment analyses suggested the existence of overlap with specific biological processes, specially the type I interferon signalling pathway. Finally, genetic correlation and polygenic risk score analyses showed cross-phenotype associations between RA and SLE.

CONCLUSIONS

In conclusion, we have identified a new risk locus shared between RA and SLE through a meta-analysis including GWAS datasets of both diseases. This study represents the first comprehensive large-scale analysis on the genetic overlap between these two complex disorders.

Identification of a Systemic Lupus Erythematosus Risk Locus Spanning ATG16L2, FCHSD2, and P2RY2 in Koreans

OBJECTIVE

Systemic lupus erythematosus (SLE) is a chronic autoimmune disorder whose etiology is incompletely understood, but likely involves environmental triggers in genetically susceptible individuals. Using an unbiased genome-wide association (GWA) scan and replication analysis, we sought to identify the genetic loci associated with SLE in a Korean population.

METHODS

A total of 1,174 SLE cases and 4,246 population controls from Korea were genotyped and analyzed with a GWA scan to identify single-nucleotide polymorphisms (SNPs) significantly associated with SLE, after strict quality control measures were applied. For select variants, replication of SLE risk loci was tested in an independent data set of 1,416 SLE cases and 1,145 population controls from Korea and China.

RESULTS

Eleven regions outside the HLA exceeded the genome-wide significance level (P = 5 × 10(-8) ). A novel SNP-SLE association was identified between FCHSD2 and P2RY2, peaking at rs11235667 (P = 1.03 × 10(-8) , odds ratio [OR] 0.59) on a 33-kb haplotype upstream of ATG16L2. In the independent replication data set, the SNP rs11235667 continued to show a significant association with SLE (replication meta-analysis P = 0.001, overall meta-analysis P = 6.67 × 10(-11) ; OR 0.63). Within the HLA region, the SNP-SLE association peaked in the class II region at rs116727542, with multiple independent effects observed in this region. Classic HLA allele imputation analysis identified HLA-DRB1*1501 and HLA-DQB1*0602, each highly correlated with one another, as most strongly associated with SLE. Ten previously established SLE risk loci were replicated: STAT1-STAT4, TNFSF4, TNFAIP3, IKZF1, HIP1, IRF5, BLK, WDFY4, ETS1, and IRAK1-MECP2. Of these loci, previously unreported, independent second risk effects of SNPs in TNFAIP3 and TNFSF4, as well as differences in the association with a putative causal variant in the WDFY4 region, were identified.

CONCLUSION

Further studies are needed to identify true SLE risk effects in other loci suggestive of a significant association, and to identify the causal variants in the regions of ATG16L2, FCHSD2, and P2RY2.

IFI44L promoter methylation as a blood biomarker for systemic lupus erythematosus

OBJECTIVE

Systemic lupus erythematosus (SLE) is a clinically heterogeneous disease with limited reliable diagnostic biomarkers. We investigated whether gene methylation could meet sensitivity and specificity criteria for a robust biomarker.

METHODS

IFI44L promoter methylation was examined using DNA samples from a discovery set including 377 patients with SLE, 358 healthy controls (HCs) and 353 patients with rheumatoid arthritis (RA). Two independent sets including 1144 patients with SLE, 1350 HCs, 429 patients with RA and 199 patients with primary Sjögren's syndrome (pSS) were used for validation.

RESULTS

Significant hypomethylation of two CpG sites within IFI44L promoter, Site1 (Chr1: 79 085 222) and Site2 (Chr1: 79 085 250; cg06872964), was identified in patients with SLE compared with HCs, patients with RA and patients with pSS. In a comparison between patients with SLE and HCs included in the first validation cohort, Site1 methylation had a sensitivity of 93.6% and a specificity of 96.8% at a cut-off methylation level of 75.5% and Site2 methylation had a sensitivity of 94.1% and a specificity of 98.2% at a cut-off methylation level of 25.5%. The IFI44L promoter methylation marker was also validated in an European-derived cohort. In addition, the methylation levels of Site1 and Site2 within IFI44L promoter were significantly lower in patients with SLE with renal damage than those without renal damage. Patients with SLE showed significantly increased methylation levels of Site1 and Site2 during remission compared with active stage.

CONCLUSIONS

The methylation level of IFI44L promoter can distinguish patients with SLE from healthy persons and other autoimmune diseases, and is a highly sensitive and specific diagnostic marker for SLE.

Decreased SMG7 expression associates with lupus-risk variants and elevated antinuclear antibody production

OBJECTIVES

Following up the systemic lupus erythematosus (SLE) genome-wide association studies (GWAS) identification of NMNAT2 at rs2022013, we fine-mapped its 150 kb flanking regions containing NMNAT2 and SMG7 in a 15 292 case-control multi-ancestry population and tested functions of identified variants.

METHODS

We performed genotyping using custom array, imputation by IMPUTE 2.1.2 and allele specific functions using quantitative real-time PCR and luciferase reporter transfections. SLE peripheral blood mononuclear cells (PBMCs) were cultured with small interfering RNAs to measure antinuclear antibody (ANA) and cyto/chemokine levels in supernatants using ELISA.

RESULTS

We confirmed association at NMNAT2 in European American (EA) and Amerindian/Hispanic ancestries, and identified independent signal at SMG7 tagged by rs2702178 in EA only (p=2.4×10^-8, OR=1.23 (95% CI 1.14 to 1.32)). In complete linkage disequilibrium with rs2702178, rs2275675 in the promoter region robustly associated with SMG7 mRNA levels in multiple expression quantitative trait locus (eQTL) datasets. Its risk allele was dose-dependently associated with decreased SMG7 mRNA levels in PBMCs of 86 patients with SLE and 119 controls (p=1.1×10^-3 and 6.8×10^-8, respectively) and conferred reduced transcription activity in transfected HEK-293 (human embryonic kidney cell line) and Raji cells (p=0.0035 and 0.0037, respectively). As a critical component in the nonsense-mediated mRNA decay pathway, SMG7 could regulate autoantigens including ribonucleoprotein (RNP) and Smith (Sm). We showed SMG7 mRNA levels in PBMCs correlated inversely with ANA titres of patients with SLE (r=-0.31, p=0.01), and SMG7 knockdown increased levels of ANA IgG and chemokine (C-C motif) ligand 19 in SLE PBMCs (p=2.0×10^-5 and 2.0×10^-4, respectively).

CONCLUSION

We confirmed NMNAT2 and identified independent SMG7 association with SLE. The inverse relationship between levels of the risk allele-associated SMG7 mRNAs and ANA suggested the novel contribution of mRNA surveillance pathway to SLE pathogenesis.

Genomic Insights into the Ancestry and Demographic History of South America

South America has a complex demographic history shaped by multiple migration and admixture events in pre- and post-colonial times. Settled over 14,000 years ago by Native Americans, South America has experienced migrations of European and African individuals, similar to other regions in the Americas. However, the timing and magnitude of these events resulted in markedly different patterns of admixture throughout Latin America. We use genome-wide SNP data for 437 admixed individuals from 5 countries (Colombia, Ecuador, Peru, Chile, and Argentina) to explore the population structure and demographic history of South American Latinos. We combined these data with population reference panels from Africa, Asia, Europe and the Americas to perform global ancestry analysis and infer the subcontinental origin of the European and Native American ancestry components of the admixed individuals. By applying ancestry-specific PCA analyses we find that most of the European ancestry in South American Latinos is from the Iberian Peninsula; however, many individuals trace their ancestry back to Italy, especially within Argentina. We find a strong gradient in the Native American ancestry component of South American Latinos associated with country of origin and the geography of local indigenous populations. For example, Native American genomic segments in Peruvians show greater affinities with Andean indigenous peoples like Quechua and Aymara, whereas Native American haplotypes from Colombians tend to cluster with Amazonian and coastal tribes from northern South America. Using ancestry tract length analysis we modeled post-colonial South American migration history as the youngest in Latin America during European colonization (9–14 generations ago), with an additional strong pulse of European migration occurring between 3 and 9 generations ago. These genetic footprints can impact our understanding of population-level differences in biomedical traits and, thus, inform future medical genetic studies in the region.

South America is home to over 400 million people who share a rich demographic history, including settlement by Native Americans, European colonization, and the African slave trade. We use genomic data to infer which populations from Europe and the Americas contributed to these admixture events. We provide evidence for multiple origins of the Native American ancestry of admixed South American Latinos. The Native American ancestral component correlates strongly with geography, indicating that admixture occurred between European colonists and local Native American populations throughout South America. We also show that the European ancestry of South American Latinos comes mainly from the Iberian peninsula, however, a significant number of Argentinians have European ancestry from other Southern European regions. The genetic signal of European admixture in South American populations is younger than the signal observed in Mexico and the Caribbean. We find evidence for a second pulse of European migration to many regions of South America subsequent to the original colonization. These results demonstrate the heterogeneous nature of the Latino population in South America and help elucidate the complex genetic and admixture events that shaped the population structure of the region.

Concordance of increased B1 cell subset and lupus phenotypes in mice and humans is dependent on BLK expression levels

Polymorphisms in the B lymphoid tyrosine kinase (BLK) gene have been associated with autoimmune diseases, including systemic lupus erythematosus, with risk correlating with reduced expression of BLK. How reduced expression of BLK causes autoimmunity is unknown. Using Blk(+/+) , Blk(+/-) , and Blk(-/-) mice, we show that aged female Blk(+/-) and Blk(-/-) mice produced higher anti-dsDNA IgG Abs and developed immune complex-mediated glomerulonephritis, compared with Blk(+/+) mice. Starting at young age, Blk(+/-) and Blk(-/-) mice accumulated increased numbers of splenic B1a cells, which differentiated into class-switched CD138(+) IgG-secreting B1a cells. Increased infiltration of B1a-like cells into the kidneys was also observed in aged Blk(+/-) and Blk(-/-) mice. In humans, we found that healthy individuals had BLK genotype-dependent levels of anti-dsDNA IgG Abs as well as increased numbers of a B1-like cell population, CD19(+)CD3(-)CD20(+)CD43(+)CD27(+), in peripheral blood. Furthermore, we describe the presence of B1-like cells in the tubulointerstitial space of human lupus kidney biopsies. Taken together, our study reveals a previously unappreciated role of reduced BLK expression on extraperitoneal accumulation of B1a cells in mice, as well as the presence of IgG autoantibodies and B1-like cells in humans.

Lupus in Latin-American patients: lessons from the GLADEL cohort

The need for comprehensive published epidemiologic and clinical data from Latin American systemic lupus erythematosus (SLE) patients motivated the late Dr Alarcón-Segovia and other Latin American professionals taking care of these patients to spearhead the creation of the G: rupo L: atino A: mericano D: e E: studio del L: upus (GLADEL) cohort in 1997. This inception cohort recruited a total of 1480 multiethnic (Mestizo, African-Latin American (ALA), Caucasian and other) SLE patients diagnosed within two years from the time of enrollment from 34 Latin American centers with expertise in the diagnosis and management of this disease. In addition to the initial 2004 description of the cohort, GLADEL has contributed to improving our knowledge about the course and outcome of lupus in patients from this part of the Americas. The major findings from this cohort are highlighted in this review. They have had important clinical implications for the adequate care of SLE patients both in Latin America and worldwide where these patients may have emigrated.

Lupus risk variants in the PXK locus alter B-cell receptor internalization

Genome wide association studies have identified variants in PXK that confer risk for humoral autoimmune diseases, including systemic lupus erythematosus (SLE or lupus), rheumatoid arthritis and more recently systemic sclerosis. While PXK is involved in trafficking of epidermal growth factor Receptor (EGFR) in COS-7 cells, mechanisms linking PXK to lupus pathophysiology have remained undefined. In an effort to uncover the mechanism at this locus that increases lupus-risk, we undertook a fine-mapping analysis in a large multi-ancestral study of lupus patients and controls. We define a large (257kb) common haplotype marking a single causal variant that confers lupus risk detected only in European ancestral populations and spans the promoter through the 3′ UTR of PXK. The strongest association was found at rs6445972 with P < 4.62 × 10⁻¹⁰, OR 0.81 (0.75–0.86). Using stepwise logistic regression analysis, we demonstrate that one signal drives the genetic association in the region. Bayesian analysis confirms our results, identifying a 95% credible set consisting of 172 variants spanning 202 kb. Functionally, we found that PXK operates on the B-cell antigen receptor (BCR); we confirmed that PXK influenced the rate of BCR internalization. Furthermore, we demonstrate that individuals carrying the risk haplotype exhibited a decreased rate of BCR internalization, a process known to impact B cell survival and cell fate. Taken together, these data define a new candidate mechanism for the genetic association of variants around PXK with lupus risk and highlight the regulation of intracellular trafficking as a genetically regulated pathway mediating human autoimmunity.