Distinct transcriptome architectures underlying lupus establishment and exacerbation.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease involving multiple immune cells. To elucidate SLE pathogenesis, it is essential to understand the dysregulated gene expression pattern linked to various clinical statuses with a high cellular resolution. Here, we conducted a large-scale transcriptome study with 6,386 RNA sequencing data covering 27 immune cell types from 136 SLE and 89 healthy donors. We profiled two distinct cell-type-specific transcriptomic signatures: disease-state and disease-activity signatures, reflecting disease establishment and exacerbation, respectively. We then identified candidate biological processes unique to each signature. This study suggested the clinical value of disease-activity signatures, which were associated with organ involvement and therapeutic responses. However, disease-activity signatures were less enriched around SLE risk variants than disease-state signatures, suggesting that current genetic studies may not well capture clinically vital biology. Together, we identified comprehensive gene signatures of SLE, which will provide essential foundations for future genomic and genetic studies.

Novel susceptibility loci for steroid-associated osteonecrosis of the femoral head in systemic lupus erythematosus.

Osteonecrosis of the femoral head (ONFH) involves necrosis of bone and bone marrow of the femoral head caused by ischemia with unknown etiology. Previous genetic studies on ONFH failed to produce consistent results, presumably because ONFH has various causes with different genetic backgrounds and the underlying diseases confounded the associations. Steroid-associated ONFH (S-ONFH) accounts for one-half of all ONFH, and systemic lupus erythematosus (SLE) is a representative disease underlying S-ONFH. We performed a genome-wide association study (GWAS) to identify genetic risk factors for S-ONFH in patients with SLE. We conducted a two-staged GWAS on 636 SLE patients with S-ONFH and 95 588 non-SLE controls. Among the novel loci identified, we determined S-ONFH-specific loci by comparing allele frequencies between SLE patients without S-ONFH and non-SLE controls. We also used Korean datasets comprising 148 S-ONFH cases and 37 015 controls to assess overall significance. We evaluated the functional annotations of significant variants by in silico analyses. The Japanese GWAS identified 4 significant loci together with 12 known SLE susceptibility loci. The four significant variants showed comparable effect sizes on S-ONFH compared with SLE controls and non-SLE controls. Three of the four loci, MIR4293/MIR1265 [odds ratio (OR) = 1.99, P-value = 1.1 × 10-9)], TRIM49/NAALAD2 (OR = 1.65, P-value = 4.8 × 10-8) and MYO16 (OR = 3.91, P-value = 4.9 × 10-10), showed significant associations in the meta-analysis with Korean datasets. Bioinformatics analyses identified MIR4293, NAALAD2 and MYO16 as candidate causal genes. MIR4293 regulates a PPARG-related adipogenesis pathway relevant to S-ONFH. We identified three novel susceptibility loci for S-ONFH in SLE.

Transcriptomic network analysis reveals key drivers of response to anti-TNF biologics in patients with rheumatoid arthritis.

OBJECTIVE: Anti-TNF biologics have been widely used to ameliorate disease activity in patients with rheumatoid arthritis (RA). However, a large fraction of patients show a poor response to these agents. Moreover, no clinically applicable predictive biomarkers have been established. This study aimed to identify response-associated biomarkers using longitudinal transcriptomic data in two independent RA cohorts. METHODS: RNA sequencing data from peripheral blood cell samples of Korean and Caucasian RA cohorts before and after initial treatment with anti-TNF biologics were analyzed to assess treatment-induced expression changes that differed between highly reliable excellent and null responders. Weighted correlation network, immune cell composition, and key driver analyses were performed to understand response-associated transcriptomic networks and cell types and their correlation with disease activity indices. RESULTS: In total, 305 response-associated genes showed significantly different treatment-induced expression changes between excellent and null responders. Co-expression network construction and subsequent key driver analysis revealed that 41 response-associated genes played a crucial role as key drivers of transcriptomic alteration in four response-associated networks involved in various immune pathways: type I interferon signalling, myeloid leucocyte activation, B cell activation, and NK cell/lymphocyte-mediated cytotoxicity. Transcriptomic response scores that we developed to estimate the individual-level degree of expression changes in the response-associated key driver genes were significantly correlated with the changes in clinical indices in independent patients with moderate or ambiguous response outcomes. CONCLUSIONS: This study provides response-specific treatment-induced transcriptomic signatures by comparing the transcriptomic landscape between patients with excellent and null responses to anti-TNF drugs at both gene and network levels.

Biological insights into systemic lupus erythematosus through an immune cell-specific transcriptome-wide association study.

OBJECTIVE: Genome-wide association studies (GWAS) have identified >100 risk loci for systemic lupus erythematosus (SLE), but the disease genes at most loci remain unclear, hampering translation of these genetic discoveries. We aimed to prioritise genes underlying the 110 SLE loci that were identified in the latest East Asian GWAS meta-analysis. METHODS: We built gene expression predictive models in blood B cells, CD4+ and CD8+ T cells, monocytes, natural killer cells and peripheral blood cells of 105 Japanese individuals. We performed a transcriptome-wide association study (TWAS) using data from the latest genome-wide association meta-analysis of 208 370 East Asians and searched for candidate genes using TWAS and three data-driven computational approaches. RESULTS: TWAS identified 171 genes for SLE (p<1.0×10-5); 114 (66.7%) showed significance only in a single cell type; 127 (74.3%) were in SLE GWAS loci. TWAS identified a strong association between CD83 and SLE (p<7.7×10-8). Meta-analysis of genetic associations in the existing 208 370 East Asian and additional 1498 cases and 3330 controls found a novel single-variant association at rs72836542 (OR=1.11, p=4.5×10-9) around CD83. For the 110 SLE loci, we identified 276 gene candidates, including 104 genes at recently-identified SLE novel loci. We demonstrated in vitro that putative causal variant rs61759532 exhibited an allele-specific regulatory effect on ACAP1, and that presence of the SLE risk allele decreased ACAP1 expression. CONCLUSIONS: Cell-level TWAS in six types of immune cells complemented SLE gene discovery and guided the identification of novel genetic associations. The gene findings shed biological insights into SLE genetic associations.

Amino acid signatures of HLA Class-I and II molecules are strongly associated with SLE susceptibility and autoantibody production in Eastern Asians.

Human leukocyte antigen (HLA) is a key genetic factor conferring risk of systemic lupus erythematosus (SLE), but precise independent localization of HLA effects is extremely challenging. As a result, the contribution of specific HLA alleles and amino-acid residues to the overall risk of SLE and to risk of specific autoantibodies are far from completely understood. Here, we dissected (a) overall SLE association signals across HLA, (b) HLA-peptide interaction, and (c) residue-autoantibody association. Classical alleles, SNPs, and amino-acid residues of eight HLA genes were imputed across 4,915 SLE cases and 13,513 controls from Eastern Asia. We performed association followed by conditional analysis across HLA, assessing both overall SLE risk and risk of autoantibody production. DR15 alleles HLA-DRB1*15:01 (P = 1.4x10-27, odds ratio (OR) = 1.57) and HLA-DQB1*06:02 (P = 7.4x10-23, OR = 1.55) formed the most significant haplotype (OR = 2.33). Conditioned protein-residue signals were stronger than allele signals and mapped predominantly to HLA-DRB1 residue 13 (P = 2.2x10-75) and its proxy position 11 (P = 1.1x10-67), followed by HLA-DRB1-37 (P = 4.5x10-24). After conditioning on HLA-DRB1, novel associations at HLA-A-70 (P = 1.4x10-8), HLA-DPB1-35 (P = 9.0x10-16), HLA-DQB1-37 (P = 2.7x10-14), and HLA-B-9 (P = 6.5x10-15) emerged. Together, these seven residues increased the proportion of explained heritability due to HLA to 2.6%. Risk residues for both overall disease and hallmark autoantibodies (i.e., nRNP: DRB1-11, P = 2.0x10-14; DRB1-13, P = 2.9x10-13; DRB1-30, P = 3.9x10-14) localized to the peptide-binding groove of HLA-DRB1. Enrichment for specific amino-acid characteristics in the peptide-binding groove correlated with overall SLE risk and with autoantibody presence. Risk residues were in primarily negatively charged side-chains, in contrast with rheumatoid arthritis. We identified novel SLE signals in HLA Class I loci (HLA-A, HLA-B), and localized primary Class II signals to five residues in HLA-DRB1, HLA-DPB1, and HLA-DQB1. These findings provide insights about the mechanisms by which the risk residues interact with each other to produce autoantibodies and are involved in SLE pathophysiology.

Large-scale meta-analysis across East Asian and European populations updated genetic architecture and variant-driven biology of rheumatoid arthritis, identifying 11 novel susceptibility loci.

OBJECTIVES: Nearly 110 susceptibility loci for rheumatoid arthritis (RA) with modest effect sizes have been identified by population-based genetic association studies, suggesting a large number of undiscovered variants behind a highly polygenic genetic architecture of RA. Here, we performed the largest-ever trans-ancestral meta-analysis with the aim to identify new RA loci and to better understand RA biology underlying genetic associations. METHODS: Genome-wide RA association summary statistics in three large case-control collections consisting of 311 292 individuals of Korean, Japanese and European populations were used in an inverse-variance-weighted fixed-effects meta-analysis. Several computational analyses using public omics resources were conducted to prioritise causal variants and genes, RA variant-implicating features (tissues, pathways and transcription factors) and potentially repurposable drugs for RA treatment. RESULTS: We identified 11 new RA susceptibility loci that explained 6.9% and 1.8% of the single-nucleotide polymorphism-based heritability in East Asians and Europeans, respectively, and confirmed 71 known non-human leukocyte antigens (HLA) susceptibility loci, identifying 90 independent association signals. The RA variants were preferentially located in binding sites of various transcription factors and in cell type-specific transcription-activation histone marks that simultaneously highlighted the importance of CD4+ T-cell activation and the potential role of non-immune organs in RA pathogenesis. A total of 615 plausible effector genes, based on gene-based associations, expression-associated variants and chromatin interaction, included targets of drugs approved for RA treatments and potentially repurposable drugs approved for other indications. CONCLUSION: Our findings provide useful insights regarding RA genetic aetiology and variant-driven RA pathogenesis.

Genetic variants shape rheumatoid arthritis-specific transcriptomic features in CD4+ T cells through differential DNA methylation, explaining a substantial proportion of heritability.

OBJECTIVE: CD4+ T cells have been suggested as the most disease-relevant cell type in rheumatoid arthritis (RA) in which RA-risk non-coding variants exhibit allele-specific effects on regulation of RA-driving genes. This study aimed to understand RA-specific signatures in CD4+ T cells using multi-omics data, interpreting inter-omics relationships in shaping the RA transcriptomic landscape. METHODS: We profiled genome-wide variants, gene expression and DNA methylation in CD4+ T cells from 82 patients with RA and 40 healthy controls using high-throughput technologies. We investigated differentially expressed genes (DEGs) and differential methylated regions (DMRs) in RA and localised quantitative trait loci (QTLs) for expression and methylation. We then integrated these based on individual-level correlations to inspect DEG-regulating sources and investigated the potential regulatory roles of RA-risk variants by a partitioned-heritability enrichment analysis with RA genome-wide association summary statistics. RESULTS: A large number of RA-specific DEGs were identified (n=2575), highlighting T cell differentiation and activation pathways. RA-specific DMRs, preferentially located in T cell regulatory regions, were correlated with the expression levels of 548 DEGs mostly in the same topologically associating domains. In addition, expressional variances in 771 and 83 DEGs were partially explained by expression QTLs for DEGs and methylation QTLs (meQTLs) for DEG-correlated DMRs, respectively. A large number of RA variants were moderately to strongly correlated with meQTLs. DEG-correlated DMRs, enriched with meQTLs, had strongly enriched heritability of RA. CONCLUSION: Our findings revealed that the methylomic changes, driven by RA heritability-explaining variants, shape the differential expression of a substantial fraction of DEGs in CD4+ T cells in patients with RA, reinforcing the importance of a multidimensional approach in disease-relevant tissues.

Meta-analysis of 208370 East Asians identifies 113 susceptibility loci for systemic lupus erythematosus.

OBJECTIVE: Systemic lupus erythematosus (SLE), an autoimmune disorder, has been associated with nearly 100 susceptibility loci. Nevertheless, these loci only partially explain SLE heritability and their putative causal variants are rarely prioritised, which make challenging to elucidate disease biology. To detect new SLE loci and causal variants, we performed the largest genome-wide meta-analysis for SLE in East Asian populations. METHODS: We newly genotyped 10 029 SLE cases and 180 167 controls and subsequently meta-analysed them jointly with 3348 SLE cases and 14 826 controls from published studies in East Asians. We further applied a Bayesian statistical approach to localise the putative causal variants for SLE associations. RESULTS: We identified 113 genetic regions including 46 novel loci at genome-wide significance (p<5×10-8). Conditional analysis detected 233 association signals within these loci, which suggest widespread allelic heterogeneity. We detected genome-wide associations at six new missense variants. Bayesian statistical fine-mapping analysis prioritised the putative causal variants to a small set of variants (95% credible set size ≤10) for 28 association signals. We identified 110 putative causal variants with posterior probabilities ≥0.1 for 57 SLE loci, among which we prioritised 10 most likely putative causal variants (posterior probability ≥0.8). Linkage disequilibrium score regression detected genetic correlations for SLE with albumin/globulin ratio (rg=-0.242) and non-albumin protein (rg=0.238). CONCLUSION: This study reiterates the power of large-scale genome-wide meta-analysis for novel genetic discovery. These findings shed light on genetic and biological understandings of SLE.

Genome-wide association study in a Korean population identifies six novel susceptibility loci for rheumatoid arthritis.

OBJECTIVE: Genome-wide association studies (GWAS) in rheumatoid arthritis (RA) have discovered over 100 RA loci, explaining patient-relevant RA pathogenesis but showing a large fraction of missing heritability. As a continuous effort, we conducted GWAS in a large Korean RA case-control population. METHODS: We newly generated genome-wide variant data in two independent Korean cohorts comprising 4068 RA cases and 36 487 controls, followed by a whole-genome imputation and a meta-analysis of the disease association results in the two cohorts. By integrating publicly available omics data with the GWAS results, a series of bioinformatic analyses were conducted to prioritise the RA-risk genes in RA loci and to dissect biological mechanisms underlying disease associations. RESULTS: We identified six new RA-risk loci (SLAMF6, CXCL13, SWAP70, NFKBIA, ZFP36L1 and LINC00158) with pmeta<5×10-8 and consistent disease effect sizes in the two cohorts. A total of 122 genes were prioritised from the 6 novel and 13 replicated RA loci based on physical distance, regulatory variants and chromatin interaction. Bioinformatics analyses highlighted potentially RA-relevant tissues (including immune tissues, lung and small intestine) with tissue-specific expression of RA-associated genes and suggested the immune-related gene sets (such as CD40 pathway, IL-21-mediated pathway and citrullination) and the risk-allele sharing with other diseases. CONCLUSION: This study identified six new RA-associated loci that contributed to better understanding of the genetic aetiology and biology in RA.

Confirmation of five novel susceptibility loci for systemic lupus erythematosus (SLE) and integrated network analysis of 82 SLE susceptibility loci.

We recently identified ten novel SLE susceptibility loci in Asians and uncovered several additional suggestive loci requiring further validation. This study aimed to replicate five of these suggestive loci in a Han Chinese cohort from Hong Kong, followed by meta-analysis (11,656 cases and 23,968 controls) on previously reported Asian and European populations, and to perform bioinformatic analyses on all 82 reported SLE loci to identify shared regulatory signatures. We performed a battery of analyses for these five loci, as well as joint analyses on all 82 SLE loci. All five loci passed genome-wide significance: MYNN (rs10936599, Pmeta = 1.92 × 10-13, OR = 1.14), ATG16L2 (rs11235604, Pmeta = 8.87 × 10 -12, OR = 0.78), CCL22 (rs223881, Pmeta = 5.87 × 10-16, OR = 0.87), ANKS1A (rs2762340, Pmeta = 4.93 × 10-15, OR = 0.87) and RNASEH2C (rs1308020, Pmeta = 2.96 × 10-19, OR = 0.84) and co-located with annotated gene regulatory elements. The novel loci share genetic signatures with other reported SLE loci, including effects on gene expression, transcription factor binding, and epigenetic characteristics. Most (56%) of the correlated (r2 > 0.8) SNPs from the 82 SLE loci were implicated in differential expression (9.81 × 10-198 < P < 5 × 10-3) of cis-genes. Transcription factor binding sites for p53, MEF2A and E2F1 were significantly (P < 0.05) over-represented in SLE loci, consistent with apoptosis playing a critical role in SLE. Enrichment analysis revealed common pathways, gene ontology, protein domains, and cell type-specific expression. In summary, we provide evidence of five novel SLE susceptibility loci. Integrated bioinformatics using all 82 loci revealed that SLE susceptibility loci share many gene regulatory features, suggestive of conserved mechanisms of SLE etiopathogenesis.

Risk for ACPA-positive rheumatoid arthritis is driven by shared HLA amino acid polymorphisms in Asian and European populations.

Previous studies have emphasized ethnically heterogeneous human leukocyte antigen (HLA) classical allele associations to rheumatoid arthritis (RA) risk. We fine-mapped RA risk alleles within the major histocompatibility complex (MHC) in 2782 seropositive RA cases and 4315 controls of Asian descent. We applied imputation to determine genotypes for eight class I and II HLA genes to Asian populations for the first time using a newly constructed pan-Asian reference panel. First, we empirically measured high imputation accuracy in Asian samples. Then we observed the most significant association in HLA-DRß1 at amino acid position 13, located outside the classical shared epitope (Pomnibus = 6.9 × 10(-135)). The individual residues at position 13 have relative effects that are consistent with published effects in European populations (His > Phe > Arg > Tyr ≅ Gly > Ser)--but the observed effects in Asians are generally smaller. Applying stepwise conditional analysis, we identified additional independent associations at positions 57 (conditional Pomnibus = 2.2 × 10(-33)) and 74 (conditional Pomnibus = 1.1 × 10(-8)). Outside of HLA-DRß1, we observed independent effects for amino acid polymorphisms within HLA-B (Asp9, conditional P = 3.8 × 10(-6)) and HLA-DPß1 (Phe9, conditional P = 3.0 × 10(-5)) concordant with European populations. Our trans-ethnic HLA fine-mapping study reveals that (i) a common set of amino acid residues confer shared effects in European and Asian populations and (ii) these same effects can explain ethnically heterogeneous classical allelic associations (e.g. HLA-DRB1*09:01) due to allele frequency differences between populations. Our study illustrates the value of high-resolution imputation for fine-mapping causal variants in the MHC.

Association of Single Nucleotide Polymorphisms of PADI4 and HLA-DRB1 Alleles with Susceptibility to Rheumatoid Arthritis-Related Lung Diseases.

OBJECTIVE: Lung diseases (LD) are common extra-articular manifestations in rheumatoid arthritis (RA). However, little is known about factors associated with susceptibility to rheumatoid arthritis-related lung diseases (RA-LD). The aim of the present study was to investigate whether the single nucleotide polymorphisms (SNPs) of PADI4 and HLA-DRB1 alleles were associated with RA-LD. METHODS: Blood samples and clinical data were collected from 116 consecutive RA patients who satisfied the 1987 American College of Rheumatology classification criteria. RA-LD was diagnosed using high-resolution computed tomography of the chest. All patients were genotyped for SNPs of PADI4 and HLA-DRB1 alleles and analyzed for full amino acid sequence of the HLA protein corresponding to a 4-digit HLA typing. Data were analyzed by independent t test (or Mann-Whitney test) for continuous variables, Chi-square test (or Fisher's exact test) and trend test for categorical variables, and logistic regression analysis. RESULTS: Ninety-four (81.0 %) RA patients had LD, of which eight (6.9 %) had interstitial lung disease (ILD) and 92 (79.3 %) had airway abnormalities in which 64 (55.2 %) showed bronchiectasis and 47 (40.5 %) revealed bronchial wall thickening. The recessive genotype of padi4_92 was susceptible to airway abnormalities (OR = 2.22, 95 % CI = 1.05-4.49, p = 0.034). Tryptophan at position 9 of HLA-DRB1 sequence was associated with the susceptibility to RA-ILD (OR = 22.89, 95 % CI = 1.20-432.56, p = 0.037). CONCLUSION: PADI4 polymorphisms and HLA-DRB1 alleles could attribute differently to the development of airway abnormalities and ILD, respectively, in RA.

Arhgef16, a novel Elmo1 binding partner, promotes clearance of apoptotic cells via RhoG-dependent Rac1 activation.

Elmo is an evolutionarily conserved mammalian ortholog of Caenorhabditis elegans CED-12 with proposed roles during the removal of apoptotic cells, cell migration, neurite outgrowth, and myoblast fusion (Katoh and Negishi (2003) [1], Park and Tosello (2007) [2], Grimsley et al. (2004) [3], Hamoud et al. (2014) [4]). Elmo mediates these cellular processes by interacting with various proteins located in the plasma membrane, cytoplasm and nucleus, and by modulating their activities although it has no intrinsic catalytic activity (Park and Tosello (2007) [2], Hamoud et al. (2014) [4], Li et al. (2013) [5], Margaron, Fradet and Cote (2013) [6], and Mauldin et al. (2013)[7]). Because there are a limited number of proteins known to interact with Elmo, we performed a yeast two-hybrid screen using Elmo1 as bait to identify Elmo1-interacting proteins and to evaluate their mode of regulation. Arhgef16 was one of the proteins identified through the screen and subsequent analyses revealed that Arhgef16 interacted with Elmo1 in mammalian cells as well. Expression of Arhgef16 in phagocytes promoted engulfment of apoptotic cells, and engulfment mediated by Arhgef16 increased synergistically in the presence of Elmo1 but was abrogated in the absence of Elmo1. In addition, Arhgef16-mediated removal of apoptotic cells was dependent on RhoG, but independent of Dock1. Taken together, this study suggests that the newly identified Elmo1-interacting protein, Arhgef16, functions synergistically with Elmo1 to promote clearance of apoptotic cells in a RhoG-dependent and Dock1-independent manner.

High-density genotyping of immune loci in Koreans and Europeans identifies eight new rheumatoid arthritis risk loci.

OBJECTIVE: A highly polygenic aetiology and high degree of allele-sharing between ancestries have been well elucidated in genetic studies of rheumatoid arthritis. Recently, the high-density genotyping array Immunochip for immune disease loci identified 14 new rheumatoid arthritis risk loci among individuals of European ancestry. Here, we aimed to identify new rheumatoid arthritis risk loci using Korean-specific Immunochip data. METHODS: We analysed Korean rheumatoid arthritis case-control samples using the Immunochip and genome-wide association studies (GWAS) array to search for new risk alleles of rheumatoid arthritis with anticitrullinated peptide antibodies. To increase power, we performed a meta-analysis of Korean data with previously published European Immunochip and GWAS data for a total sample size of 9299 Korean and 45,790 European case-control samples. RESULTS: We identified eight new rheumatoid arthritis susceptibility loci (TNFSF4, LBH, EOMES, ETS1-FLI1, COG6, RAD51B, UBASH3A and SYNGR1) that passed a genome-wide significance threshold (p<5×10(-8)), with evidence for three independent risk alleles at 1q25/TNFSF4. The risk alleles from the seven new loci except for the TNFSF4 locus (monomorphic in Koreans), together with risk alleles from previously established RA risk loci, exhibited a high correlation of effect sizes between ancestries. Further, we refined the number of single nucleotide polymorphisms (SNPs) that represent potentially causal variants through a trans-ethnic comparison of densely genotyped SNPs. CONCLUSIONS: This study demonstrates the advantage of dense-mapping and trans-ancestral analysis for identification of potentially causal SNPs. In addition, our findings support the importance of T cells in the pathogenesis and the fact of frequent overlap of risk loci among diverse autoimmune diseases.

Ethnic specificity of lupus-associated loci identified in a genome-wide association study in Korean women.

OBJECTIVES: To identify novel genetic candidates for systemic lupus erythematosus (SLE) in the Korean population, and to validate the risk loci for SLE identified in previous genome-wide association studies (GWAS). METHODS: We performed a GWAS in 400 Korean female SLE patients and 445 controls. Selected single-nucleotide polymorphisms (SNP) were then replicated in an independent cohort of 385 SLE patients and 583 controls (replication cohort 1), and in a further 811 SLE patients and 1502 controls (replication cohort 2). RESULTS: In the GWAS phase, rs9275428 located near HLA-DQB1 showed the strongest association with SLE (OR 0.50, false discovery rate (FDR) p=3.07×10(-6)). Although no loci reached genome-wide significance outside major histocompatibility complex (MHC), C8orf13-BLK, STAT4, CSMD1, DIAPH3, GLDC and TNFSF4 showed FDR p < 0.05. Our results suggest that STAT4, BLK, IRF5, PTTG1-miR-146a, UBE2L3 and TNFAIP3 are shared susceptibility loci among Caucasians and Asians, while ETS1, IKZF1, SLC15A4 are likely to be Asian-specific loci. In a combined analysis of 1596 SLE patients and 2540 controls for selected 22 candidate SNP, STAT4 and BLK as positive controls showed a strong association with SLE (FDR p=9.85×10(-13) and 2.28×10(-8), respectively). Of these, 16 candidates (PEX5L, TRAJ50, MYO18B, SOS1, ARHGAP26, SMURF1, CADPS, HAND1, FAM78B, DIAPH3, TBL1XR1, CSMD1, ZBTB20, C3orf21, HIPK1 and AP001042.1) showed only nominal significance (7.05×10(-4)≤FDR p≤4.38×10(-2)). CONCLUSIONS: There are similarities and differences in genetic susceptibility for SLE between Caucasian and Asian ethnic groups. Although 16 putative novel loci for SLE have been suggested in the Korean population, further research on a larger sample is required to discriminate truth from error.

Regulatory Variants on the Leukocyte Immunoglobulin-Like Receptor Gene Cluster are Associated with Crohn's Disease and Interact with Regulatory Variants for TAP2.

BACKGROUND AND AIMS: Crohn's disease [CD] has a complex polygenic aetiology with high heritability. There is ongoing effort to identify novel variants associated with susceptibility to CD through a genome-wide association study [GWAS] in large Korean populations. METHODS: Genome-wide variant data from 902 Korean patients with CD and 72 179 controls were used to assess the genetic associations in a meta-analysis with previous Korean GWAS results from 1621 patients with CD and 4419 controls. Epistatic interactions between CD-risk variants of interest were tested using a multivariate logistic regression model with an interaction term. RESULTS: We identified two novel genetic associations with the risk of CD near ZBTB38 and within the leukocyte immunoglobulin-like receptor [LILR] gene cluster [p < 5 × 10-8], with highly consistent effect sizes between the two independent Korean cohorts. CD-risk variants in the LILR locus are known quantitative trait loci [QTL] for multiple LILR genes, of which LILRB2 directly interacts with various ligands including MHC class I molecules. The LILR lead variant exhibited a significant epistatic interaction with CD-associated regulatory variants for TAP2 involved in the antigen presentation of MHC class I molecules [p = 4.11 × 10-4], showing higher CD-risk effects of the TAP2 variant in individuals carrying more risk alleles of the LILR lead variant (odds ratio [OR] = 0.941, p = 0.686 in non-carriers; OR = 1.45, p = 2.51 × 10-4 in single-copy carriers; OR = 2.38, p = 2.76 × 10-6 in two-copy carriers). CONCLUSIONS: This study demonstrated that genetic variants at two novel susceptibility loci and the epistatic interaction between variants in LILR and TAP2 loci confer a risk of CD.

Evaluation of TRAF6 in a large multiancestral lupus cohort.

OBJECTIVE: Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease with significant immune system aberrations resulting from complex heritable genetics as well as environmental factors. We undertook to study the role of TRAF6 as a candidate gene for SLE, since it plays a major role in several signaling pathways that are important for immunity and organ development. METHODS: Fifteen single-nucleotide polymorphisms (SNPs) across TRAF6 were evaluated in 7,490 SLE patients and 6,780 control subjects from different ancestries. Population-based case-control association analyses and meta-analyses were performed. P values, false discovery rate q values, and odds ratios (ORs) with 95% confidence intervals (95% CIs) were calculated. RESULTS: Evidence of associations was detected in multiple SNPs. The best overall P values were obtained for SNPs rs5030437 and rs4755453 (P = 7.85 × 10(-5) and P = 4.73 × 10(-5) , respectively) without significant heterogeneity among populations (P = 0.67 and P = 0.50, respectively, in Q statistic). In addition, SNP rs540386, which was previously reported to be associated with rheumatoid arthritis (RA), was found to be in linkage disequilibrium with these 2 SNPs (r(2) = 0.95) and demonstrated evidence of association with SLE in the same direction (meta-analysis P = 9.15 × 10(-4) , OR 0.89 [95% CI 0.83-0.95]). The presence of thrombocytopenia improved the overall results in different populations (meta-analysis P = 1.99 × 10(-6) , OR 0.57 [95% CI 0.45-0.72], for rs5030470). Finally, evidence of family-based association in 34 African American pedigrees with the presence of thrombocytopenia was detected in 1 available SNP (rs5030437) with a Z score magnitude of 2.28 (P = 0.02) under a dominant model. CONCLUSION: Our data indicate the presence of association of TRAF6 with SLE, consistent with the previous report of association with RA. These data provide further support for the involvement of TRAF6 in the pathogenesis of autoimmunity.

Higher Genetic Risk Loads Confer More Diverse Manifestations and Higher Risk of Lupus Nephritis in Systemic Lupus Erythematosus.

OBJECTIVE: Systemic lupus erythematosus (SLE) is a highly heritable complex disorder with heterogeneous clinical manifestations. In this study, we aimed to identify the genetic risk load using clinical and serological manifestations in SLE patients. METHODS: We genotyped a total of 1,655 Korean patients with SLE (n = 1,243 as a discovery set and n = 412 as a replication set) using a customized genome-wide single-nucleotide polymorphism (SNP) array, KoreanChip. A weighted genetic risk score (wGRS) for an individual was calculated from 112 well-validated non-HLA SNPs and HLA haplotypes of SLE-risk loci. We analyzed associations between individual wGRS and clinical SLE subphenotypes and autoantibodies using multivariable linear or logistic regression adjusted by onset age, sex, and disease duration. RESULTS: Childhood-onset SLE (<16 years) conferred the highest genetic risk compared with adult-onset (16-50 years) or late-onset (>50 years) SLE (P = 6.8 × 10-6 ). High wGRS significantly increased associations with SLE manifestations, regardless of onset age, sex, and disease duration. Individual wGRS significantly correlated positively with more clinical American College of Rheumatology criteria (ß = 0.143, P = 1.8 × 10-6 ). Subphenotype analysis revealed significant associations between the highest and lowest wGRS quartile with risk of renal disorder (hazard ratio [HR] 1.74, P = 2.2 × 10-8 ) and anti-Sm antibody production (HR 1.85, P = 2.8 × 10-5 ). Higher wGRS markedly modulated the pathogenesis of proliferative and membranous lupus nephritis class III or IV (HR 1.98, P = 1.6 × 10-5 ) and class V (HR 2.79, P = 1.0 × 10-3 ), but especially lupus nephritis class V in anti-Sm-positive SLE (area under the curve 0.68, P = 1.8 × 10-4 ). CONCLUSION: Patients with SLE and high wGRS tended to have earlier age of SLE onset, higher anti-Sm antibody positivity, and more diverse clinical phenotypes. Genetic profiling may predict high risk for lupus nephritis and a diverse clinical course in SLE patients.

Variation in the ICAM1-ICAM4-ICAM5 locus is associated with systemic lupus erythematosus susceptibility in multiple ancestries.

OBJECTIVE: Systemic lupus erythematosus (SLE; OMIM 152700) is a chronic autoimmune disease for which the aetiology includes genetic and environmental factors. ITGAM, integrin α(M) (complement component 3 receptor 3 subunit) encoding a ligand for intracellular adhesion molecule (ICAM) proteins, is an established SLE susceptibility locus. This study aimed to evaluate the independent and joint effects of genetic variations in the genes that encode ITGAM and ICAM. METHODS: The authors examined several markers in the ICAM1-ICAM4-ICAM5 locus on chromosome 19p13 and the single ITGAM polymorphism (rs1143679) using a large-scale case-control study of 17 481 unrelated participants from four ancestry populations. The single-marker association and gene-gene interaction were analysed for each ancestry, and a meta-analysis across the four ancestries was performed. RESULTS: The A-allele of ICAM1-ICAM4-ICAM5 rs3093030, associated with elevated plasma levels of soluble ICAM1, and the A-allele of ITGAM rs1143679 showed the strongest association with increased SLE susceptibility in each of the ancestry populations and the trans-ancestry meta-analysis (OR(meta)=1.16, 95% CI 1.11 to 1.22; p=4.88×10(-10) and OR(meta)=1.67, 95% CI 1.55 to 1.79; p=3.32×10(-46), respectively). The effect of the ICAM single-nucleotide polymorphisms (SNPs) was independent of the effect of the ITGAM SNP rs1143679, and carriers of both ICAM rs3093030-AA and ITGAM rs1143679-AA had an OR of 4.08 compared with those with no risk allele in either SNP (95% CI 2.09 to 7.98; p=3.91×10(-5)). CONCLUSION: These findings are the first to suggest that an ICAM-integrin-mediated pathway contributes to susceptibility to SLE.