Genome-wide association studies and fine mapping of candidate regions have rapidly advanced our understanding of the genetic basis of systemic lupus erythematosus (SLE). More than 20 robust ...associations have now been identified and confirmed, providing insights at the molecular level that refine our understanding of the involvement of host immune response processes. In addition, genes with unknown roles in SLE pathophysiology have been identified. These findings may provide new routes towards improved clinical management of this complex disease.
Since most variants that impact polygenic disease phenotypes localize to non-coding genomic regions, understanding the consequences of regulatory element variants will advance understanding of human ...disease mechanisms. Here, we report that the systemic lupus erythematosus (SLE) risk variant rs2431697 as likely causal for SLE through disruption of a regulatory element, modulating miR-146a expression. Using epigenomic analysis, genome-editing and 3D chromatin structure analysis, we show that rs2431697 tags a cell-type dependent distal enhancer specific for miR-146a that physically interacts with the miR-146a promoter. NF-kB binds the disease protective allele in a sequence-specific manner, increasing expression of this immunoregulatory microRNA. Finally, CRISPR activation-based modulation of this enhancer in the PBMCs of SLE patients attenuates type I interferon pathway activation by increasing miR-146a expression. Our work provides a strategy to define non-coding RNA functional regulatory elements using disease-associated variants and provides mechanistic links between autoimmune disease risk genetic variation and disease etiology.
The relationship of immune dysregulation and autoantibody production that may contribute to systemic lupus erythematosus (SLE) pathogenesis is unknown. This study evaluates the individual and ...combined contributions of autoantibodies, type I interferon (IFN-α) activity, and IFN-associated soluble mediators to disease development leading to SLE.
Serial serum specimens from 55 individuals collected prior to SLE classification (average timespan=4.3 years) and unaffected healthy controls matched by age (±5 years), gender, race and time of sample procurement were obtained from the Department of Defense Serum Repository. Levels of serum IFN-α activity, IFN-associated mediators and autoantibodies were evaluated and temporal relationships assessed by growth curve modelling, path analysis, analysis of covariance and random forest models.
In cases, but not matched controls, autoantibody specificities and IFN-associated mediators accumulated over a period of years, plateauing near the time of disease classification (p<0.001). Autoantibody positivity coincided with or followed type II IFN dysregulation, preceding IFN-α activity in growth curve models, with elevated IFN-α activity and B-lymphocyte stimulator levels occurring shortly before SLE classification (p≤0.005). Cases were distinguished by multivariate random forest models incorporating IFN-γ, macrophage chemoattractant protein (MCP)-3, anti-chromatin and anti-spliceosome antibodies (accuracy 93% >4 years pre-classification; 97% within 2 years of SLE classification).
Years before SLE classification, enhancement of the type II IFN pathway allows for accumulation of autoantibodies and subsequent elevations in IFN-α activity immediately preceding SLE classification. Perturbations in select immunological processes may help identify at-risk individuals for further clinical evaluation or participation in prospective intervention trials.
Ninety percent of systemic lupus erythematosus (SLE) patients are women. X chromosome-dosage increases susceptibility to SLE and primary Sjögren's syndrome (pSS). Chromosome X open reading frame 21
...escapes X-inactivation and is an SLE risk gene of previously unknown function. We undertook the present study to delineate the function of CXorf21 in the immune system as well as investigate a potential role in the sex bias of SLE and pSS.
Western blot protein analysis, qPCR, BioPlex cytokine immunoassay, pHrodo™ assays, as well as
CRISPR-Cas9 knockdown experiments were employed to delineate the role of CXorf21 in relevant immunocytes.
Expressed in monocytes and B cells,
basal Mrna, and protein expression levels are elevated in female primary monocytes, B cells, and EBV-transformed B cells compared to male cells. We also found
mRNA and protein expression is higher in both male and female cells from SLE patients compared to control subjects. TLR7 ligation increased CXorf21 protein expression and
knockdown abrogated TLR7-driven increased
mRNA expression, and reduced secretion of both TNF-alpha and IL-6 in healthy female monocytes. Similarly, we found increased pH in the lysosomes of
-deficient female monocytes.
CXorf21 is more highly expressed in female compared to male cells and is involved in a sexually dimorphic response to TLR7 activation. In addition, CXorf21 expression regulates lysosomal pH in a sexually dimorphic manner. Thus, sexually dimorphic expression of CXorf21 skews cellular immune responses in manner consistent with expected properties of a mediator of the X chromosome dose risk in SLE and pSS.
Explaining the genetics of many diseases is challenging because most associations localize to incompletely characterized regulatory regions. Using new computational methods, we show that ...transcription factors (TFs) occupy multiple loci associated with individual complex genetic disorders. Application to 213 phenotypes and 1,544 TF binding datasets identified 2,264 relationships between hundreds of TFs and 94 phenotypes, including androgen receptor in prostate cancer and GATA3 in breast cancer. Strikingly, nearly half of systemic lupus erythematosus risk loci are occupied by the Epstein-Barr virus EBNA2 protein and many coclustering human TFs, showing gene-environment interaction. Similar EBNA2-anchored associations exist in multiple sclerosis, rheumatoid arthritis, inflammatory bowel disease, type 1 diabetes, juvenile idiopathic arthritis and celiac disease. Instances of allele-dependent DNA binding and downstream effects on gene expression at plausibly causal variants support genetic mechanisms dependent on EBNA2. Our results nominate mechanisms that operate across risk loci within disease phenotypes, suggesting new models for disease origins.
Objective
Macrophage activation syndrome (MAS), a life‐threatening complication of systemic juvenile idiopathic arthritis (JIA), resembles familial hemophagocytic lymphohistiocytosis (HLH), a ...constellation of autosomal‐recessive immune disorders resulting from deficiency in cytolytic pathway proteins. We undertook this study to test our hypothesis that MAS predisposition in systemic JIA could be attributed to rare gene sequence variants affecting the cytotolytic pathway.
Methods
Whole‐exome sequencing was used in 14 patients with systemic JIA and MAS and in their parents to identify protein‐altering single‐nucleotide polymorphisms/indels in known HLH‐associated genes. To discover new candidate genes, the entire whole‐exome sequencing data were filtered to identify protein‐altering, rare recessive homozygous, compound heterozygous, and de novo variants with the potential to affect the cytolytic pathway.
Results
Heterozygous protein‐altering rare variants in the known genes (LYST,MUNC13‐4, and STXBP2) were found in 5 of 14 patients with systemic JIA and MAS (35.7%). This was in contrast to only 4 variants in 4 of 29 patients with systemic JIA without MAS (13.8%). Homozygosity and compound heterozygosity analysis applied to the entire whole‐exome sequencing data in systemic JIA/MAS revealed 3 recessive pairs in 3 genes and compound heterozygotes in 73 genes. We also identified 20 heterozygous rare protein‐altering variants that occurred in at least 2 patients. Many of the identified genes encoded proteins with a role in actin and microtubule reorganization and vesicle‐mediated transport. “Cellular assembly and organization” was the top cellular function category based on Ingenuity Pathways Analysis (P < 3.10 × 10−5).
Conclusion
Whole‐exome sequencing performed in patients with systemic JIA and MAS identified rare protein‐altering variants in known HLH‐associated genes as well as in new candidate genes.
Although much is known about the natural history of systemic lupus erythematosus (SLE), the development of SLE autoantibodies before the diagnosis of the disease has not been extensively explored. We ...investigated the onset and progression of autoantibody development before the clinical diagnosis.
The Department of Defense Serum Repository contains approximately 30 million specimens prospectively collected from more than 5 million U.S. Armed Forces personnel. We evaluated serum samples obtained from 130 persons before they received a diagnosis of SLE, along with samples from matched controls.
In 115 of the 130 patients with SLE (88 percent), at least one SLE autoantibody tested was present before the diagnosis (up to 9.4 years earlier; mean, 3.3 years). Antinuclear antibodies were present in 78 percent (at a dilution of 1:120 or more), anti-double-stranded DNA antibodies in 55 percent, anti-Ro antibodies in 47 percent, anti-La antibodies in 34 percent, anti-Sm antibodies in 32 percent, anti-nuclear ribonucleoprotein antibodies in 26 percent, and antiphospholipid antibodies in 18 percent. Antinuclear, antiphospholipid antibodies, anti-Ro, and anti-La antibodies were present earlier than anti-Sm and anti-nuclear ribonucleoprotein antibodies (a mean of 3.4 years before the diagnosis vs. 1.2 years, P=0.005). Anti-double-stranded DNA antibodies, with a mean onset 2.2 years before the diagnosis, were found later than antinuclear antibodies (P=0.06) and earlier than anti-nuclear ribonucleoprotein antibodies (P=0.005). For many patients, the earliest available serum sample was positive; therefore, these measures of the average time from the first positive antibody test to the diagnosis are underestimates of the time from the development of antibodies to the diagnosis. Of the 130 initial matched controls, 3.8 percent were positive for one or more autoantibodies.
Autoantibodies are typically present many years before the diagnosis of SLE. Furthermore, the appearance of autoantibodies in patients with SLE tends to follow a predictable course, with a progressive accumulation of specific autoantibodies before the onset of SLE, while patients are still asymptomatic.
To explore the potential influence of the polymorphic 8p23.1 inversion on known autoimmune susceptibility risk at or near BLK locus, we validated a new bioinformatics method that utilizes SNP data to ...enable accurate, high-throughput genotyping of the 8p23.1 inversion in a Caucasian population.
Principal components analysis (PCA) was performed using markers inside the inversion territory followed by k-means cluster analyses on 7416 European derived and 267 HapMaP CEU and TSI samples. A logistic regression conditional analysis was performed.
Three subgroups have been identified; inversion homozygous, heterozygous and non-inversion homozygous. The status of inversion was further validated using HapMap samples that had previously undergone Fluorescence in situ hybridization (FISH) assays with a concordance rate of above 98%. Conditional analyses based on the status of inversion were performed. We found that overall association signals in the BLK region remain significant after controlling for inversion status. The proportion of lupus cases and controls (cases/controls) in each subgroup was determined to be 0.97 for the inverted homozygous group (1067 cases and 1095 controls), 1.12 for the inverted heterozygous group (1935 cases 1717 controls) and 1.36 for non-inverted subgroups (924 cases and 678 controls). After calculating the linkage disequilibrium between inversion status and lupus risk haplotype we found that the lupus risk haplotype tends to reside on non-inversion background. As a result, a new association effect between non-inversion status and lupus phenotype has been identified ((p = 8.18×10(-7), OR = 1.18, 95%CI = 1.10-1.26).
Our results demonstrate that both known lupus risk haplotype and inversion status act additively in the pathogenesis of lupus. Since inversion regulates expression of many genes in its territory, altered expression of other genes might also be involved in the development of lupus.
Genetic variants can confer risk to complex genetic diseases by modulating gene expression through changes to the epigenome. To assess the degree to which genetic variants influence epigenome ...activity, we integrate epigenetic and genotypic data from lupus patient lymphoblastoid cell lines to identify variants that induce allelic imbalance in the magnitude of histone post-translational modifications, referred to herein as histone quantitative trait loci (hQTLs). We demonstrate that enhancer hQTLs are enriched on autoimmune disease risk haplotypes and disproportionately influence gene expression variability compared with non-hQTL variants in strong linkage disequilibrium. We show that the epigenome regulates HLA class II genes differently in individuals who carry HLA-DR3 or HLA-DR15 haplotypes, resulting in differential 3D chromatin conformation and gene expression. Finally, we identify significant expression QTL (eQTL) x hQTL interactions that reveal substructure within eQTL gene expression, suggesting potential implications for functional genomic studies that leverage eQTL data for subject selection and stratification.
Genome-wide association studies of Systemic Lupus Erythematosus (SLE) nominate 3073 genetic variants at 91 risk loci. To systematically screen these variants for allelic transcriptional enhancer ...activity, we construct a massively parallel reporter assay (MPRA) library comprising 12,396 DNA oligonucleotides containing the genomic context around every allele of each SLE variant. Transfection into the Epstein-Barr virus-transformed B cell line GM12878 reveals 482 variants with enhancer activity, with 51 variants showing genotype-dependent (allelic) enhancer activity at 27 risk loci. Comparison of MPRA results in GM12878 and Jurkat T cell lines highlights shared and unique allelic transcriptional regulatory mechanisms at SLE risk loci. In-depth analysis of allelic transcription factor (TF) binding at and around allelic variants identifies one class of TFs whose DNA-binding motif tends to be directly altered by the risk variant and a second class of TFs that bind allelically without direct alteration of their motif by the variant. Collectively, our approach provides a blueprint for the discovery of allelic gene regulation at risk loci for any disease and offers insight into the transcriptional regulatory mechanisms underlying SLE.