Two functional lupus-associated BLK promoter variants control cell-type- and developmental-stage-specific transcription

Epstein-Barr virus and genetic risk variants as determinants of T-bet+ B cell-driven autoimmune diseases

B cells expressing the transcription factor T-bet are found to have a protective role in viral infections, but are also considered major players in the onset of different types of autoimmune diseases. Currently, the exact mechanisms driving such 'atypical' memory B cells to contribute to protective immunity or autoimmunity are unclear. In addition to general autoimmune-related factors including sex and age, the ways T-bet+ B cells instigate autoimmune diseases may be determined by the close interplay between genetic risk variants and Epstein-Barr virus (EBV). The impact of EBV on T-bet+ B cells likely relies on the type of risk variants associated with each autoimmune disease, which may affect their differentiation, migratory routes and effector function. In this hypothesis-driven review, we discuss the lines of evidence pointing to such genetic and/or EBV-mediated influence on T-bet+ B cells in a range of autoimmune diseases, including systemic lupus erythematosus (SLE) and multiple sclerosis (MS). We provide examples of how genetic risk variants can be linked to certain signaling pathways and are differentially affected by EBV to shape T-bet+ B-cells. Finally, we propose options to improve current treatment of B cell-related autoimmune diseases by more selective targeting of pathways that are critical for pathogenic T-bet+ B-cell formation.

Haplotype-specific chromatin looping reveals genetic interactions of regulatory regions modulating gene expression in 8p23.1

A major goal of genetics research is to elucidate mechanisms explaining how genetic variation contributes to phenotypic variation. The genetic variants identified in genome-wide association studies (GWASs) generally explain only a small proportion of heritability of phenotypic traits, the so-called missing heritability problem. Recent evidence suggests that additional common variants beyond lead GWAS variants contribute to phenotypic variation; however, their mechanistic underpinnings generally remain unexplored. Herein, we undertake a study of haplotype-specific mechanisms of gene regulation at 8p23.1 in the human genome, a region associated with a number of complex diseases. The FAM167A-BLK locus in this region has been consistently found in the genome-wide association studies (GWASs) of systemic lupus erythematosus (SLE) in all major ancestries. Our haplotype-specific chromatin interaction (Hi-C) experiments, allele-specific enhancer activity measurements, genetic analyses, and epigenome editing experiments revealed that: 1) haplotype-specific long-range chromatin interactions are prevalent in 8p23.1; 2) BLK promoter and cis-regulatory elements cooperatively interact with haplotype-specificity; 3) genetic variants at distal regulatory elements are allele-specific modifiers of the promoter variants at FAM167A-BLK; 4) the BLK promoter interacts with and, as an enhancer-like promoter, regulates FAM167A expression and 5) local allele-specific enhancer activities are influenced by global haplotype structure due to chromatin looping. Although systemic lupus erythematosus causal variants at the FAM167A-BLK locus are thought to reside in the BLK promoter region, our results reveal that genetic variants at distal regulatory elements modulate promoter activity, changing BLK and FAM167A gene expression and disease risk. Our results suggest that global haplotype-specific 3-dimensional chromatin looping architecture has a strong influence on local allelic BLK and FAM167A gene expression, providing mechanistic details for how regional variants controlling the BLK promoter may influence disease risk.

Genome-wide association study identifies Sjögren's risk loci with functional implications in immune and glandular cells

Sjögren's disease is a complex autoimmune disease with twelve established susceptibility loci. This genome-wide association study (GWAS) identifies ten novel genome-wide significant (GWS) regions in Sjögren's cases of European ancestry: CD247, NAB1, PTTG1-MIR146A, PRDM1-ATG5, TNFAIP3, XKR6, MAPT-CRHR1, RPTOR-CHMP6-BAIAP6, TYK2, SYNGR1. Polygenic risk scores yield predictability (AUROC = 0.71) and relative risk of 12.08. Interrogation of bioinformatics databases refine the associations, define local regulatory networks of GWS SNPs from the 95% credible set, and expand the implicated gene list to >40. Many GWS SNPs are eQTLs for genes within topologically associated domains in immune cells and/or eQTLs in the main target tissue, salivary glands.

Combining SNP-to-gene linking strategies to identify disease genes and assess disease omnigenicity

Disease-associated single-nucleotide polymorphisms (SNPs) generally do not implicate target genes, as most disease SNPs are regulatory. Many SNP-to-gene (S2G) linking strategies have been developed to link regulatory SNPs to the genes that they regulate in cis. Here, we developed a heritability-based framework for evaluating and combining different S2G strategies to optimize their informativeness for common disease risk. Our optimal combined S2G strategy (cS2G) included seven constituent S2G strategies and achieved a precision of 0.75 and a recall of 0.33, more than doubling the recall of any individual strategy. We applied cS2G to fine-mapping results for 49 UK Biobank diseases/traits to predict 5,095 causal SNP-gene-disease triplets (with S2G-derived functional interpretation) with high confidence. We further applied cS2G to provide an empirical assessment of disease omnigenicity; we determined that the top 1% of genes explained roughly half of the SNP heritability linked to all genes and that gene-level architectures vary with variant allele frequency.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

Systemic lupus erythematosus as a genetic disease

Systemic lupus erythematosus is the prototypical systemic autoimmune disease, as it is characterized both by protean multi-organ system manifestations and by the uniform presence of pathogenic autoantibodies directed against components of the nucleus. Prior to the modern genetic era, the diverse clinical manifestations of SLE suggested to many that SLE patients were unlikely to share a common genetic risk basis. However, modern genetic studies have revealed that SLE usually arises when an environmental exposure occurs in an individual with a collection of genetic risk variants passing a liability threshold. Here, we summarize the current state of the field aimed at: (1) understanding the genetic architecture of this complex disease, (2) synthesizing how this genetic risk architecture impacts cellular and molecular disease pathophysiology, (3) providing illustrative examples that highlight the rich complexity of the pathobiology of this prototypical autoimmune disease and (4) communicating this complex etiopathogenesis to patients.

STAT3-mediated allelic imbalance of novel genetic variant Rs1047643 and B-cell-specific super-enhancer in association with systemic lupus erythematosus

Mapping of allelic imbalance (AI) at heterozygous loci has the potential to establish links between genetic risk for disease and biological function. Leveraging multi-omics data for AI analysis and functional annotation, we discovered a novel functional risk variant rs1047643 at 8p23 in association with systemic lupus erythematosus (SLE). This variant displays dynamic AI of chromatin accessibility and allelic expression on FDFT1 gene in B cells with SLE. We further found a B-cell restricted super-enhancer (SE) that physically contacts with this SNP-residing locus, an interaction that also appears specifically in B cells. Quantitative analysis of chromatin accessibility and DNA methylation profiles further demonstrated that the SE exhibits aberrant activity in B cell development with SLE. Functional studies identified that STAT3, a master factor associated with autoimmune diseases, directly regulates both the AI of risk variant and the activity of SE in cultured B cells. Our study reveals that STAT3-mediated SE activity and cis-regulatory effects of SNP rs1047643 at 8p23 locus are associated with B cell deregulation in SLE.

The Role of Clinical Features and Serum Biomarkers in Identifying Patients with Incomplete Lupus Erythematosus at Higher Risk of Transitioning to Systemic Lupus Erythematosus: Current Perspectives

Discovery of antinuclear antibodies (ANA) enabled earlier diagnosis of systemic lupus erythematosus (SLE) and other ANA⁺ connective tissue diseases (CTD). Rheumatologists increasingly encounter high referral volume of ANA⁺ patients. It has been estimated that only a small percentage of these patients will eventually transition to either SLE or other specified CTD. Incomplete lupus erythematosus (ILE) has been defined as a subset of patients who have some SLE-specific clinical manifestations but do not meet currently accepted classification criteria for SLE. Several studies have been performed with the goal of identifying clinical features, serum and tissue biomarkers that can distinguish those patients with ILE at risk of transitioning to SLE from those who will not. Increased autoantibody diversity, presence of anti-double-stranded DNA (dsDNA) antibodies, high expression of type I and type II interferon (IFN)-gene products, increased serum levels of B-cell-activating factor of the TNF family (BAFF), and certain serum cytokines and complement products have been identified as markers with positive predictive value, particularly when combined together. Once this patient population is better characterized biochemically, clinical trials should be considered with the primary objective to completely halt or slow down the transition from ILE to SLE. Hydroxychloroquine (HCQ) appears to be a promising agent due to its good tolerability and low toxicity profile and open-label studies in ILE patients have already shown its ability to delay the onset of SLE. Other therapeutics, like those targeting abnormal type I and type II IFN-signatures, B-cell specific signaling pathways, complement activation pathways and high BAFF levels should also be evaluated, but the risk to benefit ratio must be carefully determined before they can be considered.

SLE non-coding genetic risk variant determines the epigenetic dysfunction of an immune cell specific enhancer that controls disease-critical microRNA expression

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.

Src Family Protein Kinase Controls the Fate of B Cells in Autoimmune Diseases

There are more than 80 kinds of autoimmune diseases known at present, including rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), systemic sclerosis (SSc), inflammatory bowel disease (IBD), as well as other disorders. Autoimmune diseases have a characteristic of immune responses directly attacking own tissues, leading to systematic inflammation and subsequent tissue damage. B cells play a vital role in the development of autoimmune diseases and differentiate into plasma cells or memory B cells to secrete high-affinity antibody or provide long-lasting function. Drugs targeting B cells show good therapeutic effects for the treatment of autoimmune diseases, such as rituximab (anti-CD20 antibody). Src family protein kinases (SFKs) are believed to play important roles in a variety of cellular functions such as growth, proliferation, and differentiation of B cell via B cell antigen receptor (BCR). Lck/Yes-related novel protein tyrosine kinase (LYN), BLK (B lymphocyte kinase), and Fyn are three different kinds of SFKs mainly expressed in B cells. LYN has a dual role in the BCR signal. On the one hand, positive signals are beneficial to the development and maturation of B cells. On the other hand, LYN can also inhibit excessively activated B cells. BLK is involved in the proliferation, differentiation, and immune tolerance of B lymphocytes, and further affects the function of B cells, which may lead to autoreactive or regulatory cellular responses, increasing the risk of autoimmune diseases. Fyn may affect the development of autoimmune disorders via the differentiation of B cells in the early stage of B cell development. This article reviews the recent advances of SFKs in B lymphocytes in autoimmune diseases.

Polymorphisms of BLK are associated with renal disorder in patients with systemic lupus erythematosus

Previous studies are inconclusive on the relationships between BLK gene polymorphisms and clinical features of systemic lupus erythematosus (SLE). The present study aimed to estimate association between BLK loci and SLE clinical features in Chinese population. Associations between BLK single-nucleotide polymorphisms (SNPs) and susceptibility to SLE in this study were estimated using data of 1205 health controls previously reported in the same population. And a total of 814 SLE patients recruited from two different sources according with ACR criteria were analyzed for genotype-phenotype associations. A meta-analysis was conducted of the associations between BLK loci and renal disorder in SLE. The expression quantitative trait locus (eQTL) data were also extracted from the public databases for the selected SNPs. Significant associations were observed between these SNPs and susceptibility to SLE. In addition, the data showed that rs2618479 and rs7812879 were associated with renal disorder [OR = 1.51 (95% CI: 1.15, 1.99) and 1.61 (95% CI: 1.21, 2.14), P_corr = 0.033 and 0.011, respectively] and proteinuria [OR = 1.47 (95% CI: 1.12, 1.95) and 1.52 (95% CI: 1.14, 2.03), P_corr = 0.048 and 0.040, respectively]. The consistent associations were observed in two independent centers as well as new cases group. The result of meta-analysis for rs2618479 was also significant [OR = 1.35 (95% CI: 1.12, 1.62)]. In addition, bioinformatics analysis demonstrated that the two SNPs were significantly associated with the expression of BLK in whole blood and several immune cells. Our data support that variant loci of BLK are associated with presence of renal disorder in patients with SLE.

Association of BLK and BANK1 Polymorphisms and Interactions With Rheumatoid Arthritis in a Latin-American Population

Introduction

BLK has been identified as a risk factor to rheumatoid arthritis (RA) primarily in Asian or European-derived populations. However, this finding has not been evaluated in other populations such as Latin-Americans, except for Colombians. On the other hand, BANK1 single nucleotide variants (SNVs) have been scarcely studied in RA patients.

Objective

The aim of this study was to determine whether the BLK rs2736340T/C, rs13277113A/G, and BANK1 rs10516487G/A (R61H) and rs3733197G/A (A383T) polymorphisms are risk factors to RA in a sample of patients from Central Mexico.

Materials and Methods

We studied 957 women; 487 controls and 470 patients with RA by means of a TaqMan® SNP genotyping assay with fluorescent probes for the BLK rs13277113A/G, rs2736340T/C and BANK1 10516487G/A (R61H) and rs3733197G/A (A383T) variants.

Result

The BLK rs2736340T/C and rs13277113A/G variants were associated with risk for RA: C vs T; OR 1.39, p = 0.001, and G vs A; OR 1.37, p = 0.004, respectively. In addition, there was also an association between BANK1 R61H and RA: A vs G; OR 1.49, p = 0.003, but no with BANK1 A383T. We also identified an interaction significant between genotypes of BLK rs2736340T/C-BANK1 rs10516487G/A and RA: OR 1.65, p = 0.0001.

Conclusions

Our data suggest that both BLK and BANK1 confer susceptibility to RA in Mexican patients. The individual association of BANK1 rs1054857G/A with RA had not been previously reported in a particular population (except for pooled patients from several countries), therefore, our study presents the first evidence of association between this BANK1 variant and RA.

Rare variants in non-coding regulatory regions of the genome that affect gene expression in systemic lupus erythematosus

Personalized medicine approaches are increasingly sought for diseases with a heritable component. Systemic lupus erythematosus (SLE) is the prototypic autoimmune disease resulting from loss of immunologic tolerance, but the genetic basis of SLE remains incompletely understood. Genome wide association studies (GWAS) identify regions associated with disease, based on common single nucleotide polymorphisms (SNPs) within them, but these SNPs may simply be markers in linkage disequilibrium with other, causative mutations. Here we use an hierarchical screening approach for prediction and testing of true functional variants within regions identified in GWAS; this involved bioinformatic identification of putative regulatory elements within close proximity to SLE SNPs, screening those regions for potentially causative mutations by high resolution melt analysis, and functional validation using reporter assays. Using this approach, we screened 15 SLE associated loci in 143 SLE patients, identifying 7 new variants including 5 SNPs and 2 insertions. Reporter assays revealed that the 5 SNPs were functional, altering enhancer activity. One novel variant was linked to the relatively well characterized rs9888739 SNP at the ITGAM locus, and may explain some of the SLE heritability at this site. Our study demonstrates that non-coding regulatory elements can contain private sequence variants affecting gene expression, which may explain part of the heritability of SLE.

Associations of secreted phosphoprotein 1 and B lymphocyte kinase gene polymorphisms with autoimmune thyroid disease

BACKGROUND

Dysregulation of the type 1 interferon (IFN)-related signalling pathway predisposes one to autoimmune diseases. Possible associations of single-nucleotide polymorphisms (SNPs) of secreted phosphoprotein 1 (SPP1) and B lymphocyte kinase (BLK) of the type 1 IFN-related signalling pathway with autoimmune thyroid disease (AITD) in an ethnic Chinese (ie Taiwanese) population were tested.

METHODS

Totally, 83 Hashimoto's thyroiditis (HT) patients, 319 Graves' disease (GD) patients and 369 controls were enrolled. Genotypes of the two SNPs (rs1126772 and rs1126616) of SPP1 and two SNPs (rs13277113 and rs2736340) of BLK were determined.

RESULTS

Our results showed reduced percentages of the G allele of rs13277113 of BLK in GD (P = 0.037, odds ratio [OR] = 0.78, 95% confidence interval [CI] = 0.62-0.99) and HT (P = 0.002, OR = 0.54, 95% CI = 0.36-0.81), compared to the controls. At the same time, lower frequencies of the C allele of rs2736340 of BLK in GD (P = 0.025, OR = 0.76, 95% CI = 0.60-0.97) and HT (P = 0.003, OR = 0.53, 95% CI = 0.35-0.81) than the controls were also observed. There were significantly higher AT haplotype frequencies of rs1327713 and rs2736340 in GD and HT patients than in the controls (P = 0.025, OR = 1.31, 95% CI = 1.03-1.67, and P = 0.003, OR = 1.89, 95% CI = 1.24-2.87, respectively). Moreover, the anti-microsomal antibody titre was associated with rs2736340.

CONCLUSIONS

Genetic variants of rs13277113 and rs2736340 of BLK were associated with susceptibility to GD, HT and AITD in an ethnic Chinese population. Our results suggest the BLK may participate in the pathogenesis of GD, HT and AITD.

High-resolution genetic mapping of putative causal interactions between regions of open chromatin

Physical interaction of regulatory elements in three-dimensional space poses a challenge for studies of disease because non-coding risk variants may be great distances from the genes they regulate. Experimental methods to capture these interactions, such as chromosome conformation capture, usually cannot assign causal direction of effect between regulatory elements, an important component of fine-mapping studies. We developed a Bayesian hierarchical approach that uses two-stage least squares and applied it to an ATAC-seq (assay for transposase-accessible chromatin using sequencing) data set from 100 individuals, to identify over 15,000 high-confidence causal interactions. Most (60%) interactions occurred over <20 kb, where chromosome conformation capture-based methods perform poorly. For a fraction of loci, we identified a single variant that alters accessibility across multiple regions, and experimentally validated the BLK locus, which is associated with multiple autoimmune diseases, using CRISPR genome editing. Our study highlights how association genetics of chromatin state is a powerful approach for identifying interactions between regulatory elements.

Genetic variants in systemic lupus erythematosus susceptibility loci, XKR6 and GLT1D1 are associated with childhood-onset SLE in a Korean cohort

Impact of genetic variants on the age of systemic lupus erythematosus (SLE) onset was not fully understood. We investigated a cumulative effect of SLE-risk variants on the age of SLE onset and scanned genome-wide SNPs to search for new risk loci of childhood-onset SLE (cSLE). We analyzed 781 Korean single-center SLE subjects who previously genotyped by both Immunochip and genome-wide SNP arrays. Individual genetic risk scores (GRS) from well-validated SLE susceptibility loci were calculated and tested for their association with cSLE (<16 years at onset). Single-variant association tests were performed using a multivariable logistic regression adjusting for population stratification. GRS from SLE susceptibility loci was significantly higher in cSLE than aSLE (p = 1.23 × 10⁻³). Two SNPs, rs7460469 in XKR6 (p = 1.26 × 10⁻⁸, OR = 5.58) and rs7300146 in GLT1D1 p = 1.49 × 10⁻⁸, OR = 2.85), showed the most significant associations with cSLE. The model consisting of GRS of SLE and two newly identified loci showed an area under curve (AUC) of 0.71 in a receiver operating characteristics (ROC) curve for prediction of cSLE. In conclusion, cSLE is associated with a high cumulative SLE-risk effect and two novel SNPs rs7460469 and rs7300146, providing the first predictive model for cSLE in Koreans.

The Post-GWAS Era: From Association to Function

During the past 12 years, genome-wide association studies (GWASs) have uncovered thousands of genetic variants that influence risk for complex human traits and diseases. Yet functional studies aimed at delineating the causal genetic variants and biological mechanisms underlying the observed statistical associations with disease risk have lagged. In this review, we highlight key advances in the field of functional genomics that may facilitate the derivation of biological meaning post-GWAS. We highlight the evidence suggesting that causal variants underlying disease risk often function through regulatory effects on the expression of target genes and that these expression effects might be modest and cell-type specific. We moreover discuss specific studies as proof-of-principle examples for current statistical, bioinformatic, and empirical bench-based approaches to downstream elucidation of GWAS-identified disease risk loci.

Updates in Lupus Genetics

PURPOSE OF REVIEW

Our understanding on genetic basis of SLE has been advanced through genome-wide association studies. We review recent progress in lupus genetics with a focus on SLE-associated loci that have been functionally characterized, and discuss the potential for clinical translation of genetics data.

RECENT FINDINGS

Over 100 loci have been confirmed to show robust association with SLE and many share with other immune-mediated diseases. Although causative variants captured at these established loci are limited, they guide biological studies of gene targets for functional characterization which highlight the importance of aberrant recognition of self-nucleic acid, type I interferon overproduction, and defective immune cell signaling underlying the pathogenesis of SLE. Increasing examples illustrate a predictive value of genetic findings in susceptibility/prognosis prediction, clinical classification, and pharmacological implication. Genetic findings provide a foundation for better understanding of disease pathogenic mechanisms and opportunities for target selection in lupus drug development.

Profiling RNA-Seq at multiple resolutions markedly increases the number of causal eQTLs in autoimmune disease

Genome-wide association studies have identified hundreds of risk loci for autoimmune disease, yet only a minority (~25%) share genetic effects with changes to gene expression (eQTLs) in immune cells. RNA-Seq based quantification at whole-gene resolution, where abundance is estimated by culminating expression of all transcripts or exons of the same gene, is likely to account for this observed lack of colocalisation as subtle isoform switches and expression variation in independent exons can be concealed. We performed integrative cis-eQTL analysis using association statistics from twenty autoimmune diseases (560 independent loci) and RNA-Seq data from 373 individuals of the Geuvadis cohort profiled at gene-, isoform-, exon-, junction-, and intron-level resolution in lymphoblastoid cell lines. After stringently testing for a shared causal variant using both the Joint Likelihood Mapping and Regulatory Trait Concordance frameworks, we found that gene-level quantification significantly underestimated the number of causal cis-eQTLs. Only 5.0-5.3% of loci were found to share a causal cis-eQTL at gene-level compared to 12.9-18.4% at exon-level and 9.6-10.5% at junction-level. More than a fifth of autoimmune loci shared an underlying causal variant in a single cell type by combining all five quantification types; a marked increase over current estimates of steady-state causal cis-eQTLs. Causal cis-eQTLs detected at different quantification types localised to discrete epigenetic annotations. We applied a linear mixed-effects model to distinguish cis-eQTLs modulating all expression elements of a gene from those where the signal is only evident in a subset of elements. Exon-level analysis detected disease-associated cis-eQTLs that subtly altered transcription globally across the target gene. We dissected in detail the genetic associations of systemic lupus erythematosus and functionally annotated the candidate genes. Many of the known and novel genes were concealed at gene-level (e.g. IKZF2, TYK2, LYST). Our findings are provided as a web resource.

Genetics of Sjögren's syndrome

The pathogenesis of Sjögren's syndrome has not been elucidated. There has been evidence that genetics play an important role in the development of this disease from earlier studies. However, till now only a number of genes have been identified to be associated with SS, and these have only a weak or moderate effect. In this review we summarize the findings of the genetics studies and emphasize the need of large multicenter projects that will increase the sample sizes to provide more meaningful associations, as is the case in other common autoimmune diseases.

B-cell lymphocyte kinase polymorphisms rs13277113, rs2736340, and rs4840568 and risk of autoimmune diseases: A meta-analysis

BACKGROUND

B-cell lymphocyte kinase (BLK) is an inhibitor of B cells that has an important influence on several autoimmune diseases, but there is a lack of comprehensive analysis of its association with autoimmune diseases. Hence, it is meaningful to conduct a comprehensive analysis.

METHODS

A systematic literature search was performed on the PubMed, ScienceDirect, and Web of Science databases up to June 30, 2016. The data were extracted and quality-assessed before conducting the meta-analysis. The odds ratios (ORs) and 95% confidence intervals (95% CIs) were assessed with the STATA version 12.0 software. Subgroup and sensitivity analysis were conducted to explore potential sources of heterogeneity.

RESULTS

Altogether, 33 studies with 68,874 cases and 90,684 controls, 24 studies with 31,095 cases and 39,077 controls for rs13277113, 21 studies with 26,388 cases and 40,635 controls for rs2736340, and 4 studies with 11,391 cases and 10,972 controls for rs4840568 were included in this meta-analysis. The results revealed that the BLK rs13277113 and rs2736340 polymorphisms increased the risk of autoimmune diseases in the total analysis (A vs G: OR = 1.33, 95% CI = 1.27-1.39, P < .01; T vs C: OR = 1.34, 95% CI = 1.27-1.41, P < .01), and rs4840568 was positively associated with systemic lupus erythematosus (SLE) (A vs G: OR = 1.32, 95% CI = 1.22-1.43, P = .01).

CONCLUSION

This meta-analysis shows that the BLK (rs13277113, rs2736340, rs4840568) polymorphisms may be a risk factor for developing autoimmune diseases, especially for Asian populations and SLE.

Mapping eQTLs with RNA-seq reveals novel susceptibility genes, non-coding RNAs and alternative-splicing events in systemic lupus erythematosus

Studies attempting to functionally interpret complex-disease susceptibility loci by GWAS and eQTL integration have predominantly employed microarrays to quantify gene-expression. RNA-Seq has the potential to discover a more comprehensive set of eQTLs and illuminate the underlying molecular consequence. We examine the functional outcome of 39 variants associated with Systemic Lupus Erythematosus (SLE) through the integration of GWAS and eQTL data from the TwinsUK microarray and RNA-Seq cohort in lymphoblastoid cell lines. We use conditional analysis and a Bayesian colocalisation method to provide evidence of a shared causal-variant, then compare the ability of each quantification type to detect disease relevant eQTLs and eGenes. We discovered the greatest frequency of candidate-causal eQTLs using exon-level RNA-Seq, and identified novel SLE susceptibility genes (e.g. NADSYN1 and TCF7) that were concealed using microarrays, including four non-coding RNAs. Many of these eQTLs were found to influence the expression of several genes, supporting the notion that risk haplotypes may harbour multiple functional effects. Novel SLE associated splicing events were identified in the T-reg restricted transcription factor, IKZF2, and other candidate genes (e.g. WDFY4) through asQTL mapping using the Geuvadis cohort. We have significantly increased our understanding of the genetic control of gene-expression in SLE by maximising the leverage of RNA-Seq and performing integrative GWAS-eQTL analysis against gene, exon, and splice-junction quantifications. We conclude that to better understand the true functional consequence of regulatory variants, quantification by RNA-Seq should be performed at the exon-level as a minimum, and run in parallel with gene and splice-junction level quantification.

Multiple signals at the extended 8p23 locus are associated with susceptibility to systemic lupus erythematosus

BACKGROUND

A major systemic lupus erythematosus (SLE) susceptibility locus lies within a common inversion polymorphism region (encompassing 3.8 - 4.5  Mb) located at 8p23. Initially implicated genes included FAM167A-BLK and XKR6, of which BLK received major attention due to its known role in B-cell biology. Recently, additional SLE risk carried in non-inverted background was also reported.

OBJECTIVE AND METHODS

In this case -control study, we further investigated the 'extended' 8p23 locus (~ 4  Mb) where we observed multiple SLE signals and assessed these signals for their relation to the inversion affecting this region. The study involved a North American discovery data set (~ 1200  subjects) and a replication data set (> 10 000  subjects) comprising European-descent individuals.

RESULTS

Meta-analysis of 8p23 SNPs, with p < 0.05 in both data sets, identified 51 genome-wide significant SNPs (p < 5.0 × 10^-8). While most of these SNPs were related to previously implicated signals (XKR6-FAM167A-BLK subregion), our results also revealed two 'new' SLE signals, including SGK223-CLDN23-MFHAS1 (6.06 × 10^-9 ≤ meta p ≤ 4.88 × 10^-8) and CTSB (meta p = 4.87 × 10^-8) subregions that are located > 2 Mb upstream and ~ 0.3  Mb downstream from previously reported signals. Functional assessment of relevant SNPs indicated putative cis-effects on the expression of various genes at 8p23. Additional analyses in discovery sample, where the inversion genotypes were inferred, replicated the association of non-inverted status with SLE risk and suggested that a number of SLE risk alleles are predominantly carried in non-inverted background.

CONCLUSIONS

Our results implicate multiple (known+novel) SLE signals/genes at the extended 8p23 locus, beyond previously reported signals/genes, and suggest that this broad locus contributes to SLE risk through the effects of multiple genes/pathways.

Dependence of Glomerulonephritis Induction on Novel Intraglomerular Alternatively Activated Bone Marrow-Derived Macrophages and Mac-1 and PD-L1 in Lupus-Prone NZM2328 Mice

Glomerular damage mediated by glomerulus-infiltrating myeloid-derived cells is a key pathogenic event in lupus nephritis (LN), but the process is poorly understood. Confocal microscopy of kidney sections and flow cytometry analysis of glomerular cells from magnetic bead-purified glomeruli have identified glomerulus-infiltrating leukocyte populations in NZM2328 (NZM) lupus-prone mice with spontaneous chronic glomerulonephritis (GN) and anti-glomerular basement membrane-induced nephritis. The occurrence of a major glomerulus-infiltrating CD11b⁺F4/80^-I-A^- macrophage population exhibiting the markers programmed death ligand-1 (PD-L1), Mac-2, and macrophage mannose receptor (CD206) and producing Klf4, Il10, Retnla, Tnf, and Il6 mRNA, which are known to be expressed by alternatively activated (M2b) macrophages, correlated with proteinuria status. In NZM mice with spontaneous LN, glomerular macrophage infiltration is predominant. CD11b⁺F4/80^-I-A^- intraglomerular macrophages and polymorphonuclear neutrophils (PMN) are important in inducing GN, as anti-CD11b and -ICAM-1 mAb inhibited both proteinuria and macrophage and PMN infiltration. The predominant and high expression of PD-L1 by CD11b⁺F4/80^-I-A^- glomerular macrophages in kidneys of mice with GN and the inhibition of proteinuria by anti-PD-L1 mAb supported the pathogenic role of these macrophages but not the PD-L1^- PMN in GN development and in inducing podocyte damage. In NZM mice with spontaneous chronic GN and severe proteinuria, few glomerulus-infiltrating PMN were found, leaving macrophages and, to a less extent, dendritic cells as the major infiltrating leukocytes. Taken together, these data support the important pathogenic effect of CD11b⁺F4/80^-I-A^- M2b-like glomerulus-infiltrating macrophages in LN and reinforce macrophages as a promising target for GN treatment.

Functional variant in the promoter region of IL-27 alters gene transcription and confers a risk for ulcerative colitis in northern Chinese Han

Ulcerative colitis (UC) is a chronic inflammatory disorder of unknown etiology and a polygenic disease. IL-27 encodes p28, a subunit of IL-12 family cytokines, and has been implicated in the pathogenesis of UC. The aims of the present study were to evaluate the genetic association of a variant of the IL-27 gene with UC and to further characterize the functional variant in the IL-27 gene that influences the risk for UC. Our data demonstrated that the genetic variant rs153109 in the 5' upstream region of IL-27 is significantly associated with UC in Chinese Han individuals. Analysis of IL-27 transcripts demonstrated that individuals carrying the risk allele of rs153109 display reduced transcription of IL-27 in PBMCs. Luciferase activity assays demonstrated that the risk allele rs153109 results in decreased promoter activity compared to a non-risk allele in a tissue specific manner. Mechanistic characterization of histone modifications in the promoter region revealed that the risk haplotype tagged by the risk allele of rs153109 reduces the levels of H3K3me3 and H3K27ac.

Direct Identification of Hundreds of Expression-Modulating Variants using a Multiplexed Reporter Assay

Although studies have identified hundreds of loci associated with human traits and diseases, pinpointing causal alleles remains difficult, particularly for non-coding variants. To address this challenge, we adapted the massively parallel reporter assay (MPRA) to identify variants that directly modulate gene expression. We applied it to 32,373 variants from 3,642 cis-expression quantitative trait loci and control regions. Detection by MPRA was strongly correlated with measures of regulatory function. We demonstrate MPRA's capabilities for pinpointing causal alleles, using it to identify 842 variants showing differential expression between alleles, including 53 well-annotated variants associated with diseases and traits. We investigated one in detail, a risk allele for ankylosing spondylitis, and provide direct evidence of a non-coding variant that alters expression of the prostaglandin EP4 receptor. These results create a resource of concrete leads and illustrate the promise of this approach for comprehensively interrogating how non-coding polymorphism shapes human biology.

Genetics of Lupus Nephritis: Clinical Implications

Systemic lupus erythematosus is a heterogeneous autoimmune disease marked by the presence of pathogenic autoantibodies, immune dysregulation, and chronic inflammation that may lead to increased morbidity and early mortality from end-organ damage. More than half of all systemic lupus erythematosus patients will develop lupus nephritis. Genetic-association studies have identified more than 50 polymorphisms that contribute to lupus nephritis pathogenesis, including genetic variants associated with altered programmed cell death and defective immune clearance of programmed cell death debris. These variants may support the generation of autoantibody-containing immune complexes that contribute to lupus nephritis. Genetic variants associated with lupus nephritis also affect the initial phase of innate immunity and the amplifying, adaptive phase of the immune response. Finally, genetic variants associated with the kidney-specific effector response may influence end-organ damage and the progression to end-stage renal disease and death. This review discusses genetic insights of key pathogenic processes and pathways that may lead to lupus nephritis, as well as the clinical implications of these findings as they apply to recent advances in biologic therapies.

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.

Genetics in Sjögren Syndrome

The genes associated with Sjögren syndrome (SS) can be assigned to the NF-kB pathway, the IFN signaling pathway, lymphocyte signaling, and antigen presentation. The frequencies of risk variants show they are common with modest genetic effects. The strongest genetic association outside the human leukocyte antigen region is in IRF5, a gene relevant in the IFN signaling pathway and for B cell differentiation. Although no association has been found with the NF-kB gene itself, associations in TNFAIP3 and TNIP1 (both genome-wide significant), VCAM1 and IRAK1BP (both suggestive), point to genetic explanations for dysregulation of the NF-kB pathway in SS.

Association of FAM167A-BLK rs2736340 Polymorphism with Susceptibility to Autoimmune Diseases: A Meta-Analysis

OBJECTIVE

The purpose of this study is to evaluate the correlation between family with sequence similarity 167A-B lymphoid tyrosine kinase (FAM167A-BLK) rs2736340 polymorphism and autoimmune diseases.

METHODS

Databases including PubMed, EMBASE, Chinese National Knowledge Infrastructure (CNKI), Chinese Biomedical Literature database (CBM) and Chinese database, Wan Fang database were used in searching eligible studies from January 1, 1966 to October 2, 2015. The odds ratios (ORs) and their 95% confidence intervals (CIs) were pooled to estimate the strength of the association.

RESULTS

A total of 25 studies with 30,217 patients and 44,754 controls were included in the meta-analysis. The overall results showed FAM167A-BLK rs2736340 T allele was a risk allele for autoimmune diseases (OR 1.36, 95% CI 1.28-1.44, p < 0.001). In the subgroup by ethnicities, the results suggested T allele was an increased risk in North America, Europe, and Asia (OR 1.33, 95% CI 1.10-1.60, p = 0.004; OR 1.26, 95% CI 1.22-1.31, p < 0.001; and OR 1.46, 95% CI 1.40-1563, p < 0.001, respectively), but not in Africa. Subgroup analysis in different genetic models (recessive, dominant, and additive) revealed significant association between rs2736340 and autoimmune diseases in Asia and North America, but not the recessive model in Europe or Africa, or the additive model in Africa. Stratification analysis by diseases suggested FAM167A-BLK rs2736340 had a positive association with rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), systemic sclerosis (SSc) and Kawasaki disease, primary Sjogren's syndrome (pSS), primary antiphosholipid syndrome (APS), and myositis.

CONCLUSION

The current meta-analysis suggested that FAM167A-BLK rs2736340 polymorphism is associated with several autoimmune diseases.

Autoimmune disease-associated haplotypes of BLK exhibit lowered thresholds for B cell activation and expansion of Ig class-switched B cells

OBJECTIVE

B lymphoid kinase (BLK) is associated with rheumatoid arthritis (RA) and several other B cell-associated autoimmune disorders. BLK risk variants are consistently associated with reduced BLK expression, but the mechanisms by which reduced expression alters human B cell function to confer autoimmune disease susceptibility are unknown. This study was undertaken to characterize the BLK risk haplotype and to determine associated B cell functional phenotypes involved in autoimmunity.

METHODS

The BLK risk haplotype association with RA (determined using whole-genome sequencing data) was confirmed in 2,526 RA cases and 2,134 controls. Peripheral blood mononuclear cells (PBMCs) from RA patients, healthy adults, and umbilical cord blood were used to study B cell functional phenotypes associated with the BLK risk genotype. Association of the BLK haplotype with B cell phenotypes was analyzed using cell culture and flow cytometry.

RESULTS

Two insertion/deletions were found on the RA risk haplotype in BLK, and the reduction in BLK expression associated with the risk haplotype was confirmed in primary B lymphocytes. Carriers of the RA-associated haplotype had evidence of lower basal B cell receptor (BCR) signaling activity, yet their B cells were hyperactivatable, with enhanced up-regulation of CD86 after BCR crosslinking and greater T cell stimulatory capacity. The number of isotype-switched memory B cells was also significantly increased in subjects carrying the risk haplotype.

CONCLUSION

A major mechanism underlying the BLK association with autoimmune disease involves lowered thresholds for BCR signaling, enhanced B cell-T cell interactions, and altered patterns of isotype switching.

The SLE variant Ala71Thr of BLK severely decreases protein abundance and binding to BANK1 through impairment of the SH3 domain function

The B-lymphocyte kinase (BLK) gene is associated genetically with several human autoimmune diseases including systemic lupus erythematosus. We recently described that the genetic risk is given by two haplotypes: one covering several strongly linked single-nucleotide polymorphisms within the promoter of the gene that correlated with low transcript levels, and a second haplotype that includes a rare nonsynonymous variant (Ala71Thr). Here we show that this variant, located within the BLK SH3 domain, is a major determinant of protein levels. In vitro analyses show that the 71Thr isoform is hyperphosphorylated and promotes kinase activation. As a consequence, BLK is ubiquitinated, its proteasomal degradation enhanced and the average life of the protein is reduced by half. Altogether, these findings suggest that an intrinsic autoregulatory mechanism previously unappreciated in BLK is disrupted by the 71Thr substitution. Because the SH3 domain is also involved in protein interactions, we sought for differences between the two isoforms in trafficking and binding to protein partners. We found that binding of the 71Thr variant to the adaptor protein BANK1 is severely reduced. Our study provides new insights on the intrinsic regulation of BLK activation and highlights the dominant role of its SH3 domain in BANK1 binding.

High-density genotyping of immune-related loci identifies new SLE risk variants in individuals with Asian ancestry

Systemic lupus erythematosus (SLE) has a strong but incompletely understood genetic architecture. We conducted an association study with replication in 4,478 SLE cases and 12,656 controls from six East Asian cohorts to identify new SLE susceptibility loci and better localize known loci. We identified ten new loci and confirmed 20 known loci with genome-wide significance. Among the new loci, the most significant locus was GTF2IRD1-GTF2I at 7q11.23 (rs73366469, Pmeta = 3.75 × 10(-117), odds ratio (OR) = 2.38), followed by DEF6, IL12B, TCF7, TERT, CD226, PCNXL3, RASGRP1, SYNGR1 and SIGLEC6. We identified the most likely functional variants at each locus by analyzing epigenetic marks and gene expression data. Ten candidate variants are known to alter gene expression in cis or in trans. Enrichment analysis highlights the importance of these loci in B cell and T cell biology. The new loci, together with previously known loci, increase the explained heritability of SLE to 24%. The new loci share functional and ontological characteristics with previously reported loci and are possible drug targets for SLE therapeutics.

Genetic association analyses implicate aberrant regulation of innate and adaptive immunity genes in the pathogenesis of systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a genetically complex autoimmune disease characterized by loss of immune tolerance to nuclear and cell surface antigens. Previous genome-wide association studies (GWAS) had modest sample sizes, reducing their scope and reliability. Our study comprised 7,219 cases and 15,991 controls of European ancestry, constituting a new GWAS, a meta-analysis with a published GWAS and a replication study. We have mapped 43 susceptibility loci, including ten new associations. Assisted by dense genome coverage, imputation provided evidence for missense variants underpinning associations in eight genes. Other likely causal genes were established by examining associated alleles for cis-acting eQTL effects in a range of ex vivo immune cells. We found an over-representation (n = 16) of transcription factors among SLE susceptibility genes. This finding supports the view that aberrantly regulated gene expression networks in multiple cell types in both the innate and adaptive immune response contribute to the risk of developing SLE.

Use of next-generation DNA sequencing to analyze genetic variants in rheumatic disease

Next-generation DNA sequencing has revolutionized the field of genetics and genomics, providing researchers with the tools to efficiently identify novel rare and low frequency risk variants, which was not practical with previously available methodologies. These methods allow for the sequence capture of a specific locus or small genetic region all the way up to the entire six billion base pairs of the diploid human genome.

Rheumatic diseases are a huge burden on the US population, affecting more than 46 million Americans. Those afflicted suffer from one or more of the more than 100 diseases characterized by inflammation and loss of function, mainly of the joints, tendons, ligaments, bones, and muscles. While genetics studies of many of these diseases (for example, systemic lupus erythematosus, rheumatoid arthritis, and inflammatory bowel disease) have had major successes in defining their genetic architecture, causal alleles and rare variants have still been elusive. This review describes the current high-throughput DNA sequencing methodologies commercially available and their application to rheumatic diseases in both case–control as well as family-based studies.

Immunogenetics of systemic lupus erythematosus: A comprehensive review

Our understanding of the genetic basis of systemic lupus erythematosus has progressed rapidly in recent years. While many genetic polymorphisms have been associated with disease susceptibility, the next major step involves integrating these genetic polymorphisms into the molecular mechanisms and cellular immunology of the human disease. In this review, we summarize some recent work in this area, including the genetics of the type I IFN response in SLE, including polygenic and monogenic factors, as well as epigenetic influences. Contributions of both HLA and non-HLA polymorphisms to the complex genetics of SLE are reviewed. We also review recent reports of specific gene deficits leading to monogenic SLE-like syndromes. The molecular functions of common SLE-risk variants are reviewed in depth, including regulatory variations in promoter and enhancer elements and coding-change polymorphisms, and studies which are beginning to define the molecular and cellular functions of these polymorphisms in the immune system. We discuss epigenetic influences on lupus, with an emphasis on micro-RNA expression and binding, as well as epigenetic modifications that regulate the expression levels of various genes involved in SLE pathogenesis and the ways epigenetic marks modify SLE susceptibility genes. The work summarized in this review provides a fascinating window into the biology and molecular mechanisms of human SLE. Understanding the functional mechanisms of causal genetic variants underlying the human disease greatly facilitates our ability to translate genetic associations toward personalized care, and may identify new therapeutic targets relevant to human SLE disease mechanisms.

An in silico Approach Reveals Associations between Genetic and Epigenetic Factors within Regulatory Elements in B Cells from Primary Sjögren's Syndrome Patients

Recent advances in genetics have highlighted several regions and candidate genes associated with primary Sjögren’s syndrome (SS), a systemic autoimmune epithelitis that combines exocrine gland dysfunctions, and focal lymphocytic infiltrations. In addition to genetic factors, it is now clear that epigenetic deregulations are present during SS and restricted to specific cell type subsets, such as lymphocytes and salivary gland epithelial cells. In this study, 72 single nucleotide polymorphisms (SNPs) associated with 43 SS gene risk factors were selected from publicly available and peer reviewed literature for further in silico analysis. SS risk variant location was tested revealing a broad distribution in coding sequences (5.6%), intronic sequences (55.6%), upstream/downstream genic regions (30.5%), and intergenic regions (8.3%). Moreover, a significant enrichment of regulatory motifs (promoter, enhancer, insulator, DNAse peak, and expression quantitative trait loci) characterizes SS risk variants (94.4%). Next, screening SNPs in high linkage disequilibrium (r² ≥ 0.8 in Caucasians) revealed 645 new variants including 5 SNPs with missense mutations, and indicated an enrichment of transcriptionally active motifs according to the cell type (B cells > monocytes > T cells ≫ A549). Finally, we looked at SS risk variants for histone markers in B cells (GM12878), monocytes (CD14⁺) and epithelial cells (A548). Active histone markers were associated with SS risk variants at both promoters and enhancers in B cells, and within enhancers in monocytes. In conclusion and based on the obtained in silico results that need further confirmation, associations were observed between SS genetic risk factors and epigenetic factors and these associations predominate in B cells, such as those observed at the FAM167A–BLK locus.

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.

Advances in lupus genetics and epigenetics

PURPOSE OF REVIEW

Genome-wide association studies have identified more than 50 robust loci associated with systemic lupus erythematosus (SLE) susceptibility, and follow-up studies help reveal candidate causative genetic variants and their biological relevance contributing to the development of SLE. Epigenetic modulation is emerging as an important mechanism for understanding how the implicated genes interact with environmental factors. We review recent progress toward identifying causative variants of SLE-associated loci and epigenetic impact on lupus, especially genetic-epigenetic interactions that modulate expression levels of SLE susceptibility genes.

RECENT FINDINGS

A few SLE-risk loci have been refined to localize likely causative variants responsible for the observed genome-wide association study signals. Few of such variants disrupt coding sequences resulting in gain or loss of function for the encoded protein, whereas most fall in noncoding regions with potential to regulate gene expression through alterations in transcriptional activity, splicing, mRNA stability and epigenetic modifications. Multiple key pathways related to the SLE pathogenesis have been indicated by the identified genetic risk factors, including type I interferon signaling pathway that can also be regulated by epigenetic changes occurred in SLE.

SUMMARY

These findings provide novel insights into the disease pathogenesis and promise better diagnostic accuracy and new therapeutic targets for patient management.

TYK2 protein-coding variants protect against rheumatoid arthritis and autoimmunity, with no evidence of major pleiotropic effects on non-autoimmune complex traits

Despite the success of genome-wide association studies (GWAS) in detecting a large number of loci for complex phenotypes such as rheumatoid arthritis (RA) susceptibility, the lack of information on the causal genes leaves important challenges to interpret GWAS results in the context of the disease biology. Here, we genetically fine-map the RA risk locus at 19p13 to define causal variants, and explore the pleiotropic effects of these same variants in other complex traits. First, we combined Immunochip dense genotyping (n = 23,092 case/control samples), Exomechip genotyping (n = 18,409 case/control samples) and targeted exon-sequencing (n = 2,236 case/controls samples) to demonstrate that three protein-coding variants in TYK2 (tyrosine kinase 2) independently protect against RA: P1104A (rs34536443, OR = 0.66, P = 2.3x10^-21), A928V (rs35018800, OR = 0.53, P = 1.2x10^-9), and I684S (rs12720356, OR = 0.86, P = 4.6x10^-7). Second, we show that the same three TYK2 variants protect against systemic lupus erythematosus (SLE, P_omnibus = 6x10^-18), and provide suggestive evidence that two of the TYK2 variants (P1104A and A928V) may also protect against inflammatory bowel disease (IBD; P_omnibus = 0.005). Finally, in a phenome-wide association study (PheWAS) assessing >500 phenotypes using electronic medical records (EMR) in >29,000 subjects, we found no convincing evidence for association of P1104A and A928V with complex phenotypes other than autoimmune diseases such as RA, SLE and IBD. Together, our results demonstrate the role of TYK2 in the pathogenesis of RA, SLE and IBD, and provide supporting evidence for TYK2 as a promising drug target for the treatment of autoimmune diseases.

Genetics and pathogenesis of systemic lupus erythematosus and lupus nephritis

Systemic lupus erythematosus (SLE) is a multisystem autoimmune disorder that has a broad spectrum of effects on the majority of organs, including the kidneys. Approximately 40-70% of patients with SLE will develop lupus nephritis. Renal assault during SLE is initiated by genes that breach immune tolerance and promote autoantibody production. These genes might act in concert with other genetic factors that augment innate immune signalling and IFN-I production, which in turn can generate an influx of effector leucocytes, inflammatory mediators and autoantibodies into end organs, such as the kidneys. The presence of cognate antigens in the glomerular matrix, together with intrinsic molecular abnormalities in resident renal cells, might further accentuate disease progression. This Review discusses the genetic insights and molecular mechanisms for key pathogenic contributors in SLE and lupus nephritis. We have categorized the genes identified in human studies of SLE into one of four pathogenic events that lead to lupus nephritis. We selected these categories on the basis of the cell types in which these genes are expressed, and the emerging paradigms of SLE pathogenesis arising from murine models. Deciphering the molecular basis of SLE and/or lupus nephritis in each patient will help physicians to tailor specific therapies.

The effect of inversion at 8p23 on BLK association with lupus in Caucasian population

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. Methods: 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. Results: 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⁻⁷, OR = 1.18, 95%CI = 1.10–1.26). Conclusion: 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.

Genetics of systemic lupus erythematosus: immune responses and end organ resistance to damage

Systemic lupus erythematosus (SLE) is a prototypic systemic autoimmune disorder. Considerable progress has been made to delineate the genetic control of this complex disorder. In this review, selected aspects of human and mouse genetics related to SLE are reviewed with emphasis on genes that contribute to both innate and adaptive immunity and to genes that contribute directly to susceptibility to end organ damage. It is concluded that the interactions among these two major pathways will provide further insight into the pathogenesis of SLE. An interactive model of the two major pathways is proposed without emphasis on the importance of breaking tolerance to autoantigens.

Illumina human exome genotyping array clustering and quality control

With the rise of high-throughput sequencing technology, traditional genotyping arrays are gradually being replaced by sequencing technology. Against this trend, Illumina has introduced an exome genotyping array that provides an alternative approach to sequencing, especially suited to large-scale genome-wide association studies (GWASs). The exome genotyping array targets the exome plus rare single-nucleotide polymorphisms (SNPs), a feature that makes it substantially more challenging to process than previous genotyping arrays that targeted common SNPs. Researchers have struggled to generate a reliable protocol for processing exome genotyping array data. The Vanderbilt Epidemiology Center, in cooperation with Vanderbilt Technologies for Advanced Genomics Analysis and Research Design (VANGARD), has developed a thorough exome chip-processing protocol. The protocol was developed during the processing of several large exome genotyping array-based studies, which included over 60,000 participants combined. The protocol described herein contains detailed clustering techniques and robust quality control procedures, and it can benefit future exome genotyping array-based GWASs.