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.

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.

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.

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.

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.

Recent advances in understanding the genetic basis of systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a polygenic chronic autoimmune disease leading to multiple organ damage. A large heritability of up to 66% is estimated in SLE, with roughly 180 reported susceptibility loci that have been identified mostly by genome-wide association studies (GWASs) and account for approximately 30% of genetic heritability. A vast majority of risk variants reside in non-coding regions, which makes it quite challenging to interpret their functional implications in the SLE-affected immune system, suggesting the importance of understanding cell type-specific epigenetic regulation around SLE GWAS variants. The latest genetic studies have been highly fruitful as several dozens of SLE loci were newly discovered in the last few years and many loci have come to be understood in systemic approaches integrating GWAS signals with other biological resources. In this review, we summarize SLE-associated genetic variants in both the major histocompatibility complex (MHC) and non-MHC loci, examining polygenetic risk scores for SLE and their associations with clinical features. Finally, variant-driven pathogenetic functions underlying genetic associations are described, coupled with discussion about challenges and future directions in genetic studies on SLE.

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.

EphA2 Interacts with Tim-4 through Association between Its FN3 Domain and the IgV Domain of Tim-4.

Tim-4 promotes the engulfment of apoptotic cells or exogenous particles by securing them on phagocytes. It is unable to transduce signals by itself but helps other engulfment receptors sense and internalize them. However, the identity of the engulfment receptors collaborating with Tim-4 is still incompletely understood. In this study, we searched for a candidate transmembrane protein with a FN3 domain, important for interaction with Tim-4, in silico and investigated whether it indeed interacts with Tim-4 and is involved in Tim-4-mediated phagocytosis. We found that EphA2 containing a FN3 domain in the extracellular region interacted with Tim-4, which was mediated by the IgV domain of Tim-4 and the FN3 domain of EphA2. Nevertheless, we found that EphA2 expression failed to alter Tim-4-mediated phagocytosis of apoptotic cells or polystyrene beads. Taken together, our findings suggest that EphA2, a new Tim-4 interacting protein, may intervene in a Tim-4-mediated cellular event even if it is not phagocytosis of endogenous or exogenous particles and vice versa.

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.

Clinical and Genetic Risk Factors Associated With the Presence of Lupus Nephritis.

Objective: To elucidate whether clinical features and the weighted genetic risk score (wGRS) were associated with the presence of lupus nephritis (LN). Methods: We retrospectively divided patients with systemic lupus erythematosus (SLE, n=1,078) into biopsy-proven LN (n=507) and non-LN groups (non-LN, n=571) Baseline clinical features, serologic markers, and the wGRS were collected The wGRS was calculated from 112 non-human leukocyte antigen (non-HLA) loci and HLA-DRß1 amino acid haplotypes for SLE Associations among clinical features, wGRS, and the presence of LN were identified. Results: In the multivariate analysis, patients with LN were younger at diagnosis (odds ratio [OR]=0.97, p<0.001), had more pleuritis (OR=2.44, p<0.001) and pericarditis (OR=1.62, p=0.029), had a higher detection rate of anti-double stranded deoxyribonucleic acid (anti-dsDNA antibodies, OR=2.22, p<0.001), anti-Smith antibodies (anti-Sm antibodies, OR=1.70, p=0.002), low level of complement (OR=1.37, p=0.043) and absence of antiphospholipid antibodies (aPL antibodies, OR=1.60, p=0.002), and had higher wGRS (OR=1.16, p=0.012) Mediation analysis suggested that anti-Sm antibodies and low complement could be mediators in the relationship between high wGRS and the presence of LN. Conclusion: Onset age, pleuritis, pericarditis, several serologic markers, and wGRS were associated with the presence of LN Anti-Sm antibodies and low complement appeared to mediate the indirect relationship between wGRS and the presence of LN.

Allele-Specific Quantification of HLA-DRB1 Transcripts Reveals Imbalanced Allelic Expression That Modifies the Amino Acid Effects in HLA-DRß1.

OBJECTIVE: HLA association fine-mapping studies have shown the effects of missense variants in HLA-DRB1 on rheumatoid arthritis (RA) susceptibility, prognosis, and autoantibody production. However, the phenotypic effects of expression changes in HLA-DRB1 remain poorly understood. Therefore, we investigated the allele-specific expression of HLA-DRB1 and its effect on an HLA-DRß1 structure-associated trait in RA. METHODS: We quantified the allele-specific expression of each HLA-DRB1 3-field classic allele in 48 Korean RA patients with anti-citrullinated protein antibodies (ACPAs) and 319 healthy European subjects by using both RNA sequencing and HLA-DRB1 genotype data to calculate the relative expression strength of multiple HLA-DRB1 alleles (n = 14 in Koreans and n = 25 in Europeans) in each population. The known association between ACPA level and alanine at position 74 of HLA-DRß1 in ACPA-positive RA was revisited to understand the phenotypic effect of allele-specific expression of HLA-DRB1 by modeling multivariate logistic regression with the genomic dosage or relative expression dosage of Ala-74 in 2 independent sets of 1,723 Korean RA patients with ACPA. RESULTS: The relative expression strength was highly allele-specific, causing imbalanced allelic expression in HLA-DRB1 heterozygotes. The association between HLA-DRß1 Ala-74 and ACPA level in RA was better explained by relative expression dosage of Ala-74 than by the genomic dosage (change in Akaike's information criterion = -6.98). Moreover, the expression variance of Ala-74 in Ala-74 heterozygotes with no genomic variance of Ala-74 was significantly associated with ACPA level (P = 2.26 × 10-3 ). CONCLUSION: Our findings illustrate the advantage of integrating quantitative and qualitative changes in HLA-DRB1 into a single model for understanding HLA-DRB1 associations.

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.

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.

Host Genetic and Gut Microbial Signatures in Familial Inflammatory Bowel Disease.

INTRODUCTION: The family history of inflammatory bowel disease (IBD) has been strongly associated with risk of developing IBD. This study aimed to identify the host genetic and gut microbial signatures in familial IBD. METHODS: Genetic analyses using genome-wide single nucleotide polymorphism genotyping and whole exome sequencing were performed to calculate weighted genetic risk scores from known IBD-associated common variants and to identify rare deleterious protein-altering variants specific to patients with familial IBD in 8 Korean families that each included more than 2 affected first-degree relatives (FDRs) and their unaffected FDR(s). In parallel, gut microbial community was analyzed by 16S rRNA sequencing of stools from the sample individuals. RESULTS: The risk of familial IBD was not well explained by the genetic burden from common IBD-risk variants, suggesting the presence of family-shared genetic and environmental disease-risk factors. We identified 17 genes (AC113554.1, ACE, AKAP17A, AKAP9, ANK2, ASB16, ASIC3, DNPH1, DUS3L, FAM200A, FZD10, LAMA5, NUTM2F, PKN1, PRR26, WDR66, and ZC3H4) that each contained rare, potentially deleterious variants transmitted to the affected FDRs in multiple families. In addition, metagenomic analyses revealed significantly different diversity of gut microbiota and identified a number of differentially abundant taxa in affected FDRs, highlighting 22 novel familial disease-associated taxa with large abundance changes and the previously reported gut dysbiosis including low alpha diversity in IBD and 16 known IBD-specific taxa. DISCUSSION: This study identified familial IBD-associated rare deleterious variants and gut microbial dysbiosis in familial IBD.

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.

MHC associations of ankylosing spondylitis in East Asians are complex and involve non-HLA-B27 HLA contributions.

BACKGROUND: The association of HLA-B*27 with AS is amongst the strongest of any known association of a common variant with any human disease. Nonetheless, there is strong evidence indicating that other HLA-B alleles are involved in the disease. European ethnicity studies have demonstrated risk associations with HLA-B*40 and multiple other HLA-B, HLA-A, and HLA class II alleles, and demonstrated that in that ethnic group, the amino acid sequence at position 97 in HLA-B is the key determinant of HLA associations with AS. A recent study in Korean AS cases and controls additionally identified association at HLA-C*15:02. In the current study, we examined the MHC associations of AS in an expanded East Asian cohort. METHODS: A total of 1637 Chinese, Taiwanese, and Korean AS cases meeting the modified New York Criteria for AS, and 1589 ethnically matched controls, were genotyped with the Illumina Immunochip, including a dense coverage of the MHC region. HLA genotypes and amino acid composition were imputed using the SNP2HLA programme using the Han-MHC reference panel based on the data of Han Chinese subjects (n = 9689), and association tested using logistic regression controlling for population stratification effects. RESULTS: A strong association was seen with HLA-B*27 (odds ratio (OR) = 205.3, P = 5.76 × 10-244). Controlling for this association, the strongest risk association is seen with HLA-C*15 at genome-wide significant level (OR = 7.62, P = 9.30 × 10-19), and confirmed association is also seen with HLA-B*40 at suggestive level (OR = 1.65, P = 2.54 × 10-4). At amino acid level, the strongest association seen in uncontrolled analysis was with histidine at position 114 in HLA-B (P = 7.24 × 10-241), but conditional analyses suggest that the primary amino acid associations are with lysine at position 70 and asparagine at position 97. Restriction of the ERAP1 association with HLA-B27-positive AS, previously reported in European subjects, was confirmed in East Asians. CONCLUSIONS: This study confirms in East Asians that the HLA associations of AS are multiple, including previously reported associations at HLA-B*27, HLA-B*40, and HLA-C*15, as well as novel association with HLA-DQB1*04. The HLA-B associations are driven by the amino acids at positions 70 and 97, in the B pocket of HLA-B.

Identifying damage clusters in patients with systemic lupus erythematosus.

AIM: Systemic lupus erythematosus (SLE) causes irreversible damage to organ systems. Recently, evidence has been obtained for subphenotypes of SLE. This study aimed to identify damage clusters and compare the associated clinical manifestations, SLE disease activity, mortality, and genetic risk scores (GRS). METHODS: The study was conducted on the Hanyang BAE lupus cohort. Patients with disease duration <5 years were excluded to minimize confounding effects of disease duration. They were grouped into 3 clusters based on the Systemic Lupus International Collaborating Clinics Damage Index using k-means cluster analysis. RESULTS: Among the 1130 analyzed patients, musculoskeletal damage was most prevalent (20.2%), followed by ocular (11.4%), renal (10.5%), and neuropsychiatric damage (10.2%). Three significantly different damage clusters were identified. Patients in cluster 1 (n = 824) showed the least damage. Cluster 2 (n = 195) was characterized by frequent renal (55.4%) and ocular (58.0%) damage, and cluster 3 (n = 111) was dominated by neuropsychiatric (100%) and musculoskeletal damage (35.1%). Cluster 2 had the highest adjusted mean AMS (adjusted mean SLE Disease Activity Index score; mean ± SD: 5.4 ± 2.9), while cluster 3 had the highest mortality (14.4%). Weighted GRS did not differ significantly between the clusters. CONCLUSION: Patients in prevalent renal and ocular damage cluster had the highest AMS scores, while the cluster with frequent neuropsychiatric damage had the highest mortality.

Update on the Genetics of Systemic Lupus Erythematosus: Genome-Wide Association Studies and Beyond.

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease of complex etiology that primarily affects women of childbearing age. The development of SLE is attributed to the breach of immunological tolerance and the interaction between SLE-susceptibility genes and various environmental factors, resulting in the production of pathogenic autoantibodies. Working in concert with the innate and adaptive arms of the immune system, lupus-related autoantibodies mediate immune-complex deposition in various tissues and organs, leading to acute and chronic inflammation and consequent end-organ damage. Over the past two decades or so, the impact of genetic susceptibility on the development of SLE has been well demonstrated in a number of large-scale genetic association studies which have uncovered a large fraction of genetic heritability of SLE by recognizing about a hundred SLE-susceptibility loci. Integration of genetic variant data with various omics data such as transcriptomic and epigenomic data potentially provides a unique opportunity to further understand the roles of SLE risk variants in regulating the molecular phenotypes by various disease-relevant cell types and in shaping the immune systems with high inter-individual variances in disease susceptibility. In this review, the catalogue of SLE susceptibility loci will be updated, and biological signatures implicated by the SLE-risk variants will be critically discussed. It is optimistically hoped that identification of SLE risk variants will enable the prognostic and therapeutic biomarker armamentarium of SLE to be strengthened, a major leap towards precision medicine in the management of the condition.

Endothelium-Dependent Vasorelaxant Effect of Prunus Persica Branch on Isolated Rat Thoracic Aorta.

Peach (Prunus persica (L.) Batsch) is a popular fruit consumed by people worldwide, owing to its pleasant flavor and high mineral nutrient content. A few plants from the genus Prunus, such as Prunus yedoensis, Prunus cerasus, and Prunus serotina have shown vasorelaxant and vasodilatory effects, to date, no study has investigated the vasorelaxation effects of the P. persica branch extract (PPE). The vasorelaxant effect of PPE was endothelium-dependent, and it was related to the NO-sGC-cGMP, vascular prostacyclin, and muscarinic receptor transduction pathway. K+ channels, such as the BKCa, KV, and KATP channels, were partially associated with PPE-induced vasorelaxation. PPE was effective in relaxing serotonin (5-HT)- or angiotensin II-induced contraction; furthermore, PPE attenuated Ca2+-induced vasoconstriction by IP3 receptors in the SR membrane, but its vasorelaxant effect was not associated with the influx of extracellular Ca2+ via receptor-operative Ca2+ channels or voltage-dependent Ca2+ channels. Recognizing the rising use of functional foods for hypertension treatment, our findings imply that PPE may be a natural antihypertensive agent.