The muscle tissue transcriptome of idiopathic inflammatory myopathy reflects the muscle damage process by monocytes and presence of skin lesions.

OBJECTIVE: To investigate transcriptomic and immunophenotypic features of muscle specimens from patients with idiopathic inflammatory myopathy (IIM). METHODS: Bulk RNA-sequencing was performed on muscle biopsy samples from 16 patients with dermatomyositis (DM) and 9 patients with polymyositis (PM). Seven tested positive for anti-aminoacyl t-RNA synthetase antibodies in the DM patients (ARS-DM). We conducted weighted gene coexpression network analysis (WGCNA), differentially expressed gene (DEG) analysis, and gene set variation analysis (GSVA) to assess contributions of specific pathways. Cell proportions in muscle specimens were estimated using a deconvolution approach. RESULTS: WGCNA revealed significant positive correlations between serum creatine kinase (CK) levels and gene modules involved in cellular respiration, phagocytosis, and oxidative phosphorylation (OXPHOS). Significant positive correlations were also observed between CK levels and proportions of CD16-positive and -negative monocytes and myeloid dendritic cells. Notably, DM patients demonstrated enrichment of complement and interferon-α and -γ pathway genes compared to those with PM. Furthermore, ARS-DM demonstrated a higher proportion of Th1 cells and DEGs related to OXPHOS. Additionally, serum Krebs von den Lungen-6 levels correlated with gene modules associated with extracellular matrix and transforming growth factor-ß signaling pathway. CONCLUSION: Our study highlights a significant involvement of monocytes in muscle damage and delineates pathological differences among IIM subtypes. DM was characterized by complement, interferon-α and -γ signaling, whilst ARS-DM was associated with OXPHOS. Distinctive gene expression variations in muscle specimens suggest that different pathologic mechanisms underlie muscle damage in each IIM phenotype.

Pannexin1 channel-dependent secretome from apoptotic tumor cells shapes immune-escape microenvironment.

Apoptotic cell death is a critical step in organism development and tissue homeostasis. Apoptotic cells affect immune cell activities in normal tissues. It is not clear whether similar cell death machinery causes tumor environments to evade anti-tumor immune responses. Here, using a mouse transplant model, we found a large number of tumor cells undergoing intrinsic apoptosis in tumors derived from the 4T1 breast cancer cell line, where neutrophils significantly accumulated. Interestingly, these apoptotic 4T1 tumor cells directly induced neutrophil extracellular traps (NETs) in a pannexin 1 (Panx1) channel-dependent manner, and knockdown of Panx1 in 4T1 cells led to a reduction in tumor size. Spermidine released through Panx1 from apoptotic 4T1 cells induced NETs in bone marrow-derived neutrophils in vitro. In addition, inhibition of spermidine synthesis suppressed tumor growth in the mouse transplant model. Collectively, our data suggested a new immune-escape mechanism for tumors by Panx1-mediated secretome from intrinsic apoptotic cells, which may provide a new therapeutic target for cancer.

Progesterone stimulates cortisol production in the maturing bovine cumulus-oocyte complex.

In the bovine cumulus oophorus, 11ß-hydroxysteroid dehydrogenase type 1 (HSD11B1)-mediated cortisol production dramatically increases during the periovulatory period. This event is closely associated with increased progesterone (P4) production, implying a functional connection between these C21 steroids. In this study, we investigated the mutual regulation of P4 and cortisol production in the bovine cumulus oophorus. Bovine cumulus-oocyte complexes (COCs) were aspirated from follicles 2-5 mm in diameter and subjected to in vitro maturation (IVM) for 24 h in an M199 supplemented with fetal calf serum (FCS) and follicle-stimulating hormone (FSH). COCs were treated with trilostane (0, 0.1, 1, 10 mM), an inhibitor of P4 synthesis, RU486 (0, 0.1, 1, 10 mM), a receptor antagonist for the progesterone receptor (PR) and glucocorticoid receptor (GR), and various concentrations of a synthetic progestogen nomegestrol acetate (NA; 0, 0.001, 0.01, 0.1, 1, 10 mM) to examine effect of P4. The effects of cortisol (0, 0.1, 1, 10 mM) were also examined in the presence or absence of trilostane. Trilostane and RU486 suppressed cumulus expansion, cortisol production, and HSD11B1 but not hexose-6-phosphate dehydrogenase (H6PDH) expression. Concomitant treatment with NA reversed the effects of trilostane. Unlike NA, cortisol did not alter the antagonistic effects of trilostane on cumulus expansion and HSD11B1 expression. Cortisol did not affect P4 production or steroidogenic acute regulatory protein (STAR), cholesterol side-chain cleavage enzyme (CYP11A1), 3ß-hydroxysteroid dehydrogenase type 1 (HSD3B1), and HSD11B1 expression. Collectively, these results indicate that locally produced P4 is crucial in regulating the local glucocorticoid environment through PRtg in the maturing bovine cumulus oophorus. Cortisol, however, does not appear to regulate P4 or its production.

Role of Lysine-Specific Demethylase 1 in Metabolically Integrating Osteoclast Differentiation and Inflammatory Bone Resorption Through Hypoxia-Inducible Factor 1α and E2F1.

OBJECTIVE: Hypoxia occurs in tumors, infections, and sites of inflammation, such as in the affected joints of patients with rheumatoid arthritis (RA). It alleviates inflammatory responses and increases bone resorption in inflammatory arthritis by enhancing osteoclastogenesis. The mechanism by which the hypoxia response is linked to osteoclastogenesis and inflammatory bone resorption is unclear. This study was undertaken to evaluate whether the protein lysine-specific demethylase 1 (LSD1) metabolically integrates inflammatory osteoclastogenesis and bone resorption in a state of inflammatory arthritis. METHODS: LSD1-specific inhibitors and gene silencing with small interfering RNAs were used to inhibit the expression of LSD1 in human osteoclast precursor cells derived from CD14-positive monocytes, with subsequent assessment by RNA-sequencing analysis. In experimental mouse models of arthritis, inflammatory osteolysis, or osteoporosis, features of accelerated bone loss and inflammatory osteolysis were analyzed. Furthermore, in blood samples from patients with RA, cis-acting expression quantitative trait loci (cis-eQTL) were analyzed for association with the expression of hypoxia-inducible factor 1α (HIF-1α), and associations between HIF-1α allelic variants and extent of bone erosion were evaluated. RESULTS: In human osteoclast precursor cells, RANKL induced the expression of LSD1 in a mechanistic target of rapamycin-dependent manner. Expression of LSD1 was higher in synovium from RA patients than in synovium from osteoarthritis patients. Inhibition of LSD1 in human osteoclast precursors suppressed osteoclast differentiation. Results of transcriptome analysis identified several LSD1-mediated hypoxia and cell-cycle pathways as key genetic pathways involved in human osteoclastogenesis. Furthermore, HIF-1α protein, which is rapidly degraded by the proteasome in a normoxic environment, was found to be expressed in RANKL-stimulated osteoclast precursor cells. Induction of LSD1 by RANKL stabilized the expression of HIF-1α protein, thereby promoting glycolysis, in conjunction with up-regulation of the transcription factor E2F1. Analyses of cis-eQTL revealed that higher HIF-1α expression was associated with increased bone erosion in patients with RA. Inhibition of LSD1 decreased pathologic bone resorption in mice, both in models of accelerated osteoporosis and models of arthritis and inflammatory osteolysis. CONCLUSION: LSD1 metabolically regulates osteoclastogenesis in an energy-demanding inflammatory environment. These findings provide potential new therapeutic strategies targeting osteoclasts in the management of inflammatory arthritis, including in patients with RA.

Parsing multiomics landscape of activated synovial fibroblasts highlights drug targets linked to genetic risk of rheumatoid arthritis.

OBJECTIVES: Synovial fibroblasts (SFs) are one of the major components of the inflamed synovium in rheumatoid arthritis (RA). We aimed to gain insight into the pathogenic mechanisms of SFs through elucidating the genetic contribution to molecular regulatory networks under inflammatory condition. METHODS: SFs from RA and osteoarthritis (OA) patients (n=30 each) were stimulated with eight different cytokines (interferon (IFN)-α, IFN-γ, tumour necrosis factor-α, interleukin (IL)-1ß, IL-6/sIL-6R, IL-17, transforming growth factor-ß1, IL-18) or a combination of all 8 (8-mix). Peripheral blood mononuclear cells were fractioned into five immune cell subsets (CD4+ T cells, CD8+ T cells, B cells, natural killer (NK) cells, monocytes). Integrative analyses including mRNA expression, histone modifications (H3K27ac, H3K4me1, H3K4me3), three-dimensional (3D) genome architecture and genetic variations of single nucleotide polymorphisms (SNPs) were performed. RESULTS: Unstimulated RASFs differed markedly from OASFs in the transcriptome and epigenome. Meanwhile, most of the responses to stimulations were shared between the diseases. Activated SFs expressed pathogenic genes, including CD40 whose induction by IFN-γ was significantly affected by an RA risk SNP (rs6074022). On chromatin remodelling in activated SFs, RA risk loci were enriched in clusters of enhancers (super-enhancers; SEs) induced by synergistic proinflammatory cytokines. An RA risk SNP (rs28411362), located in an SE under synergistically acting cytokines, formed 3D contact with the promoter of metal-regulatory transcription factor-1 (MTF1) gene, whose binding motif showed significant enrichment in stimulation specific-SEs. Consistently, inhibition of MTF1 suppressed cytokine and chemokine production from SFs and ameliorated mice model of arthritis. CONCLUSIONS: Our findings established the dynamic landscape of activated SFs and yielded potential therapeutic targets associated with genetic risk of RA.

Chromosomal alterations among age-related haematopoietic clones in Japan.

The extent to which the biology of oncogenesis and ageing are shaped by factors that distinguish human populations is unknown. Haematopoietic clones with acquired mutations become common with advancing age and can lead to blood cancers1-10. Here we describe shared and population-specific patterns of genomic mutations and clonal selection in haematopoietic cells on the basis of 33,250 autosomal mosaic chromosomal alterations that we detected in 179,417 Japanese participants in the BioBank Japan cohort and compared with analogous data from the UK Biobank. In this long-lived Japanese population, mosaic chromosomal alterations were detected in more than 35.0% (s.e.m., 1.4%) of individuals older than 90 years, which suggests that such clones trend towards inevitability with advancing age. Japanese and European individuals exhibited key differences in the genomic locations of mutations in their respective haematopoietic clones; these differences predicted the relative rates of chronic lymphocytic leukaemia (which is more common among European individuals) and T cell leukaemia (which is more common among Japanese individuals) in these populations. Three different mutational precursors of chronic lymphocytic leukaemia (including trisomy 12, loss of chromosomes 13q and 13q, and copy-neutral loss of heterozygosity) were between two and six times less common among Japanese individuals, which suggests that the Japanese and European populations differ in selective pressures on clones long before the development of clinically apparent chronic lymphocytic leukaemia. Japanese and British populations also exhibited very different rates of clones that arose from B and T cell lineages, which predicted the relative rates of B and T cell cancers in these populations. We identified six previously undescribed loci at which inherited variants predispose to mosaic chromosomal alterations that duplicate or remove the inherited risk alleles, including large-effect rare variants at NBN, MRE11 and CTU2 (odds ratio, 28-91). We suggest that selective pressures on clones are modulated by factors that are specific to human populations. Further genomic characterization of clonal selection and cancer in populations from around the world is therefore warranted.

Large-scale genome-wide association study in a Japanese population identifies novel susceptibility loci across different diseases.

The overwhelming majority of participants in current genetic studies are of European ancestry. To elucidate disease biology in the East Asian population, we conducted a genome-wide association study (GWAS) with 212,453 Japanese individuals across 42 diseases. We detected 320 independent signals in 276 loci for 27 diseases, with 25 novel loci (P < 9.58 × 10-9). East Asian-specific missense variants were identified as candidate causal variants for three novel loci, and we successfully replicated two of them by analyzing independent Japanese cohorts; p.R220W of ATG16L2 (associated with coronary artery disease) and p.V326A of POT1 (associated with lung cancer). We further investigated enrichment of heritability within 2,868 annotations of genome-wide transcription factor occupancy, and identified 378 significant enrichments across nine diseases (false discovery rate < 0.05) (for example, NKX3-1 for prostate cancer). This large-scale GWAS in a Japanese population provides insights into the etiology of complex diseases and highlights the importance of performing GWAS in non-European populations.

Citrullination of RGG Motifs in FET Proteins by PAD4 Regulates Protein Aggregation and ALS Susceptibility.

Recent proteome analyses have provided a comprehensive overview of various posttranslational modifications (PTMs); however, PTMs involving protein citrullination remain unclear. We performed a proteomic analysis of citrullinated proteins, and we identified more than 100 PAD4 (peptidyl arginine deiminase 4) substrates. Approximately one-fifth of the PAD4 substrates contained an RG/RGG motif, and PAD4 competitively inhibited the methylation of the RGG motif in FET proteins (FUS, EWS, and TAF15) and hnRNPA1, which are causative genes for ALS (amyotrophic lateral sclerosis). PAD4-mediated citrullination significantly inhibited the aggregation of FET proteins, a frequently observed feature in neurodegenerative diseases. FUS protein levels in arsenic-induced stress granules were significantly increased in Padi4-/- mouse embryonic fibroblasts (MEFs). Moreover, rs2240335 was associated with low expression of PADI4 in the brain and a high risk of ALS (p = 0.0381 and odds ratio of 1.072). Our findings suggest that PAD4-mediated RGG citrullination plays a key role in protein solubility and ALS pathogenesis.

Macrophage extracellular trap formation promoted by platelet activation is a key mediator of rhabdomyolysis-induced acute kidney injury.

Rhabdomyolysis is a serious syndrome caused by skeletal muscle injury and the subsequent release of breakdown products from damaged muscle cells into systemic circulation. The muscle damage most often results from strenuous exercise, muscle hypoxia, medications, or drug abuse and can lead to life-threatening complications, such as acute kidney injury (AKI). Rhabdomyolysis and the AKI complication can also occur during crush syndrome, an emergency condition that commonly occurs in victims of natural disasters, such as earthquakes, and man-made disasters, such as wars and terrorism. Myoglobin released from damaged muscle is believed to trigger renal dysfunction in this form of AKI. Recently, macrophages were implicated in the disease pathogenesis of rhabdomyolysis-induced AKI, but the precise molecular mechanism remains unclear. In the present study, we show that macrophages released extracellular traps (ETs) comprising DNA fibers and granule proteins in a mouse model of rhabdomyolysis. Heme-activated platelets released from necrotic muscle cells during rhabdomyolysis enhanced the production of macrophage extracellular traps (METs) through increasing intracellular reactive oxygen species generation and histone citrullination. Here we report, for the first time to our knowledge, this unanticipated role for METs and platelets as a sensor of myoglobin-derived heme in rhabdomyolysis-induced AKI. This previously unknown mechanism might be targeted for treatment of the disease. Finally, we found a new therapeutic tool for prevention of AKI after rhabdomyolysis, which might rescue some sufferers of this pathology.

Splicing variant of WDFY4 augments MDA5 signalling and the risk of clinically amyopathic dermatomyositis.

OBJECTIVES: Idiopathic inflammatory myopathies (IIMs) are a heterogeneous group of rare autoimmune diseases in which both genetic and environmental factors play important roles. To identify genetic factors of IIM including polymyositis, dermatomyositis (DM) and clinically amyopathic DM (CADM), we performed the first genome-wide association study for IIM in an Asian population. METHODS: We genotyped and tested 496 819 single nucleotide polymorphism for association using 576 patients with IIM and 6270 control subjects. We also examined the causal mechanism of disease-associated variants by in silico analyses using publicly available data sets as well as by in in vitro analyses using reporter assays and apoptosis assays. RESULTS: We identified a variant in WDFY4 that was significantly associated with CADM (rs7919656; OR=3.87; P=1.5×10-8). This variant had a cis-splicing quantitative trait locus (QTL) effect for a truncated WDFY4isoform (tr-WDFY4), with higher expression in the risk allele. Transexpression QTL analysis of this variant showed a positive correlation with the expression of NF-κB associated genes. Furthermore, we demonstrated that both WDFY4 and tr-WDFY4 interacted with pattern recognition receptors such as TLR3, TLR4, TLR9 and MDA5 and augmented the NF-κB activation by these receptors. WDFY4 isoforms also enhanced MDA5-induced apoptosis to a greater extent in the tr-WDFY4-transfected cells. CONCLUSIONS: As CADM is characterised by the appearance of anti-MDA5 autoantibodies and severe lung inflammation, the WDFY4 variant may play a critical role in the pathogenesis of CADM.

Polygenic burdens on cell-specific pathways underlie the risk of rheumatoid arthritis.

Recent evidence suggests that a substantial portion of complex disease risk alleles modify gene expression in a cell-specific manner. To identify candidate causal genes and biological pathways of immune-related complex diseases, we conducted expression quantitative trait loci (eQTL) analysis on five subsets of immune cells (CD4+ T cells, CD8+ T cells, B cells, natural killer (NK) cells and monocytes) and unfractionated peripheral blood from 105 healthy Japanese volunteers. We developed a three-step analytical pipeline comprising (i) prediction of individual gene expression using our eQTL database and public epigenomic data, (ii) gene-level association analysis and (iii) prediction of cell-specific pathway activity by integrating the direction of eQTL effects. By applying this pipeline to rheumatoid arthritis data sets, we identified candidate causal genes and a cytokine pathway (upregulation of tumor necrosis factor (TNF) in CD4+ T cells). Our approach is an efficient way to characterize the polygenic contributions and potential biological mechanisms of complex diseases.

Decreased severity of experimental autoimmune arthritis in peptidylarginine deiminase type 4 knockout mice.

BACKGROUND: Peptidylarginine deiminase type 4 (PADI4) has been identified as a susceptibility gene for rheumatoid arthritis (RA) by genome-wide association studies. PADI4 is highly expressed in the bone marrow, macrophages, neutrophils, and monocytes. Peptidyl citrulline is an interesting molecule in RA because it is a target antigen for anti-citrullinated peptide antibodies, and only PADs (translated proteins from PADI genes) can provide peptidyl citrulline via the modification of protein substrates. The aim of this study was to evaluate the importance of the PADI4 gene in the progression of RA. METHODS: We generated Padi4 knockout (Padi4(-/-)) DBA1J mice. The Padi4(-/-) DBA1J and wild-type mice were immunized with bovine type II collagen (CII) to develop collagen-induced arthritis (CIA). The expression of various inflammatory cytokines and Padi genes in immune cells was detected by the real-time TaqMan assay. Cytokine concentrations in sera were measured by enzyme-linked immunosorbent assays. Localization of the PAD4 and PAD2 proteins was indicated by immunohistochemistry. RESULTS: We demonstrated that the clinical disease score was significantly decreased in the Padi4(-/-) mice and Padi4 expression was induced by CII immunization. In the Padi4(-/-) mice, serum anti-type II collagen (CII) immunoglobulin M (IgM), IgG, and inflammatory cytokine levels were significantly decreased compared with those in the wild-type mice. Padi2 expression was induced in the immune cells of the Padi4(-/-) mice as a compensation for the defect in Padi4. CONCLUSIONS: Padi4 affected disease severity in the CIA mice and was involved in the enhancement of the collagen-initiated inflammatory responses.

EAF2 mediates germinal centre B-cell apoptosis to suppress excessive immune responses and prevent autoimmunity.

Regulated apoptosis of germinal centre (GC) B cells is critical for normal humoral immune responses. ELL-associated factor 2 (EAF2) regulates transcription elongation and has been shown to be an androgen-responsive potential tumour suppressor in prostate by inducing apoptosis. Here we show that EAF2 is selectively upregulated in GC B cells among various immune cell types and promotes apoptosis of GC B cells both in vitro and in vivo. EAF2 deficiency results in enlarged GCs and elevated antibody production during a T-dependent immune response. After immunization with type II collagen, mice lacking EAF2 produce high levels of collagen-specific autoantibodies and rapidly develop severe arthritis. Moreover, the mutant mice spontaneously produce anti-dsDNA, rheumatoid factor and anti-nuclear antibodies as they age. These results demonstrate that EAF2-mediated apoptosis in GC B cells limits excessive humoral immune responses and is important for maintaining self-tolerance.

Peptidylarginine deiminase type 4 deficiency reduced arthritis severity in a glucose-6-phosphate isomerase-induced arthritis model.

Peptidyl arginine deiminase 4 (PAD4) is an enzyme that is involved in protein citrullination, and is a target for anti-citrullinated peptide antibodies (ACPAs) in rheumatoid arthritis (RA). Genetic polymorphisms in the PADI4 gene encoding PAD4 are associated with RA susceptibility. We herein analyzed the roles of PADI4 in inflammatory arthritis using a glucose-6-phosphate isomerase (GPI)-induced arthritis (GIA) model in Padi4 knockout (KO) mice. Arthritis severity, serum anti-GPI antibody titers, and IL-6 concentrations were significantly reduced in Padi4 KO mice. The frequency of Th17 cells was decreased in GPI-immunized Padi4 KO mice, whereas WT and Padi4-deficient naïve CD4(+) T cells displayed the same efficiencies for Th17 cell differentiation in vitro. In addition, the numbers of myeloid lineage cells were reduced with the increased expression of pro-apoptotic genes in GPI-immunized Padi4 KO mice. Furthermore, the survival of Padi4-deficient neutrophils was impaired in vitro. Our results suggest that PADI4 exacerbates arthritis with diverse immunological modifications.

Genetic basis of rheumatoid arthritis: a current review.

Rheumatoid arthritis (RA) is one of the most common autoimmune diseases. As with other complex traits, genome-wide association studies (GWASs) have tremendously enhanced our understanding of the complex etiology of RA. In this review, we describe the genetic architecture of RA as determined through GWASs and meta-analyses. In addition, we discuss the pathologic mechanism of the disease by examining the combined findings of genetic and functional studies of individual RA-associated genes, including HLA-DRB1, PADI4, PTPN22, TNFAIP3, STAT4, and CCR6. Moreover, we briefly examine the potential use of genetic data in clinical practice in RA treatment, which represents a challenge in medical genetics in the post-GWAS era.

Genetics of rheumatoid arthritis contributes to biology and drug discovery.

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ∼10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2 - 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses--as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes--to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

PAD4 regulates proliferation of multipotent haematopoietic cells by controlling c-myc expression.

Peptidylarginine deiminase 4 (PAD4) functions as a transcriptional coregulator by catalyzing the conversion of histone H3 arginine residues to citrulline residues. Although the high level of PAD4 expression in bone marrow cells suggests its involvement in haematopoiesis, its precise contribution remains unclear. Here we show that PAD4, which is highly expressed in lineage(-) Sca-1(+) c-Kit(+) (LSK) cells of mouse bone marrow compared with other progenitor cells, controls c-myc expression by catalyzing the citrullination of histone H3 on its promoter. Furthermore, PAD4 is associated with lymphoid enhancer-binding factor 1 and histone deacetylase 1 at the upstream region of the c-myc gene. Supporting these findings, LSK cells, especially multipotent progenitors, in PAD4-deficient mice show increased proliferation in a cell-autonomous fashion compared with those in wild-type mice. Together, our results strongly suggest that PAD4 regulates the proliferation of multipotent progenitors in the bone marrow by controlling c-myc expression.

An association analysis of HLA-DRB1 with systemic lupus erythematosus and rheumatoid arthritis in a Japanese population: effects of *09:01 allele on disease phenotypes.

Objective. To re-evaluate the roles of HLA-DRB1 alleles in susceptibility to SLE and RA and their effects on autoantibody status in large-scale Japanese cohorts. Methods. A total of 656 SLE, 2410 RA and 911 control subjects, who were all Japanese, were genotyped for HLA-DRB1 alleles using sequence-specific oligonucleotide probes. The association of alleles with disease susceptibility was tested by logistic regression analysis and by the relative predispositional effect method. The association with autoantibody status was examined by the standard χ(2) test. Results. HLA-DRB1*15:01, *09:01, *08:02 and *04:01 were significantly associated with SLE susceptibility, while shared epitope (SE) alleles and DRB1*09:01 were associated with RA susceptibility. The compound heterozygote of DRB1*09:01/*15:01 conferred an increased risk for SLE compared with the homozygotes for DRB1*09:01 and *15:01 and was associated with earlier onset of disease, whereas the compound effect of DRB1-SE/*09:01 was not clear in RA. DRB1*09:01 was significantly associated with the appearance of anti-Sm antibody in SLE as well as ACPA in RA, while protectively associated with anti-dsDNA antibody in SLE. No significant interaction was observed between DRB1*09:01 and smoking status for the appearance of ACPA, unlike that observed in SE alleles in RA. Conclusion. We identified HLA-DRB1 alleles associated with SLE and RA in a Japanese population and demonstrated a shared susceptibility of DRB1*09:01 between the diseases as well as its effect on autoantibody production.

Functional variants in NFKBIE and RTKN2 involved in activation of the NF-κB pathway are associated with rheumatoid arthritis in Japanese.

Rheumatoid arthritis is an autoimmune disease with a complex etiology, leading to inflammation of synovial tissue and joint destruction. Through a genome-wide association study (GWAS) and two replication studies in the Japanese population (7,907 cases and 35,362 controls), we identified two gene loci associated with rheumatoid arthritis susceptibility (NFKBIE at 6p21.1, rs2233434, odds ratio (OR) = 1.20, P = 1.3 × 10(-15); RTKN2 at 10q21.2, rs3125734, OR = 1.20, P = 4.6 × 10(-9)). In addition to two functional non-synonymous SNPs in NFKBIE, we identified candidate causal SNPs with regulatory potential in NFKBIE and RTKN2 gene regions by integrating in silico analysis using public genome databases and subsequent in vitro analysis. Both of these genes are known to regulate the NF-κB pathway, and the risk alleles of the genes were implicated in the enhancement of NF-κB activity in our analyses. These results suggest that the NF-κB pathway plays a role in pathogenesis and would be a rational target for treatment of rheumatoid arthritis.

Use of a multiethnic approach to identify rheumatoid- arthritis-susceptibility loci, 1p36 and 17q12.

We have previously shown that rheumatoid arthritis (RA) risk alleles overlap between different ethnic groups. Here, we utilize a multiethnic approach to show that we can effectively discover RA risk alleles. Thirteen putatively associated SNPs that had not yet exceeded genome-wide significance (p < 5 × 10(-8)) in our previous RA genome-wide association study (GWAS) were analyzed in independent sample sets consisting of 4,366 cases and 17,765 controls of European, African American, and East Asian ancestry. Additionally, we conducted an overall association test across all 65,833 samples (a GWAS meta-analysis plus the replication samples). Of the 13 SNPs investigated, four were significantly below the study-wide Bonferroni corrected p value threshold (p < 0.0038) in the replication samples. Two SNPs (rs3890745 at the 1p36 locus [p = 2.3 × 10(-12)] and rs2872507 at the 17q12 locus [p = 1.7 × 10(-9)]) surpassed genome-wide significance in all 16,659 RA cases and 49,174 controls combined. We used available GWAS data to fine map these two loci in Europeans and East Asians, and we found that the same allele conferred risk in both ethnic groups. A series of bioinformatic analyses identified TNFRSF14-MMEL1 at the 1p36 locus and IKZF3-ORMDL3-GSDMB at the 17q12 locus as the genes most likely associated with RA. These findings demonstrate empirically that a multiethnic approach is an effective strategy for discovering RA risk loci, and they suggest that combining GWASs across ethnic groups represents an efficient strategy for gaining statistical power.