CD4+ and B Lymphocyte Expression Quantitative Traits at Rheumatoid Arthritis Risk Loci in Patients With Untreated Early Arthritis: Implications for Causal Gene Identification

Multi-omics data integration reveals novel genes related to autoimmune hypothyroidism in the brain: A molecular basis for the brain-thyroid axis

BACKGROUND

The mechanisms underlying the complex relationship between autoimmune hypothyroidism and neurological disorders remain unclear. We conducted a comprehensive analysis of associations between alternative splicing, transcriptomics, and proteomics data and autoimmune hypothyroidism.

METHODS

Splicing-wide association studies (SWAS), proteome-wide association studies (PWAS), and transcriptome-wide association studies (TWAS) were used to identify genes and proteins that regulate autoimmune hypothyroidism within the brain axis. We performed TWAS on GTEx V8 thyroid tissue data to identify autoimmune hypothyroidism-associated thyroid axis genes. A FUSION analysis of overlapping genes in the brain and thyroid axes and brain splicing weights was conducted to determine the influence of alternative splicing in the brain on thyroid tissue gene expression.

RESULTS

SWAS identified 223 alternative splicing events, TWAS identified 270 genes, and PWAS revealed five genes (FDPS, PPIL3, PEX6, MMAB, and ALDH2) encoding proteins associated with autoimmune hypothyroidism. Neuroimaging analyses revealed distinct brain-imaging phenotypes associated with these five genes. TWAS of thyroid tissue identified four genes (FDPS, PPIL3, MMAB, and ALDH2) associated with the brain axis related to thyroid tissue. A FUSION analysis indicated that alternative splicing changes in ALDH2 in brain tissue influenced its expression in thyroid tissue.

CONCLUSION

Integrating brain splicing, proteomic, and transcriptomic data supports the association between specific genes and proteins in the brain and autoimmune hypothyroidism. Additionally, ALDH2 alternative splicing in brain tissue influences its thyroid tissue expression. These findings provide new insights into the molecular basis of autoimmune hypothyroidism, facilitating future pathogenesis research.

Expression quantitative trait loci analysis in rheumatoid arthritis identifies tissue specific variants associated with severity and outcome

OBJECTIVE

Genome-wide association studies have successfully identified more than 100 loci associated with susceptibility to rheumatoid arthritis (RA). However, our understanding of the functional effects of genetic variants in causing RA and their effects on disease severity and response to treatment remains limited.

METHODS

In this study, we conducted expression quantitative trait locus (eQTL) analysis to dissect the link between genetic variants and gene expression comparing the disease tissue against blood using RNA-Sequencing of synovial biopsies (n=85) and blood samples (n=51) from treatment-naïve patients with RA from the Pathobiology of Early Arthritis Cohort.

RESULTS

This identified 898 eQTL genes in synovium and genes loci in blood, with 232 genes in common to both synovium and blood, although notably many eQTL were tissue specific. Examining the HLA region, we uncovered a specific eQTL at HLA-DPB2 with the critical triad of single-nucleotide polymorphisms (SNPs) rs3128921 driving synovial HLA-DPB2 expression, and both rs3128921 and HLA-DPB2 gene expression correlating with clinical severity and increasing probability of the lympho-myeloid pathotype.

CONCLUSIONS

This analysis highlights the need to explore functional consequences of genetic associations in disease tissue. HLA-DPB2 SNP rs3128921 could potentially be used to stratify patients to more aggressive treatment immediately at diagnosis.

New insights into the suppression of inflammation and lipid accumulation by JAZF1

Atherosclerosis is one of the leading causes of disease and death worldwide. The identification of new therapeutic targets and agents is critical. JAZF1 is expressed in many tissues and is found at particularly high levels in adipose tissue (AT). JAZF1 suppresses inflammation (including IL-1β, IL-4, IL-6, IL-8, IL-10, TNFα, IFN-γ, IAR-20, COL3A1, laminin, and MCP-1) by reducing NF-κB pathway activation and AT immune cell infiltration. JAZF1 reduces lipid accumulation by regulating the liver X receptor response element (LXRE) of the SREBP-1c promoter, the cAMP-response element (CRE) of HMGCR, and the TR4 axis. LXRE and CRE sites are present in many cytokine and lipid metabolism gene promoters, which suggests that JAZF1 regulates these genes through these sites. NF-κB is the center of the JAZF1-mediated inhibition of the inflammatory response. JAZF1 suppresses NF-κB expression by suppressing TAK1 expression. Interestingly, TAK1 inhibition also decreases lipid accumulation. A dual-targeting strategy of NF-κB and TAK1 could inhibit both inflammation and lipid accumulation. Dual-target compounds (including prodrugs) 1-5 exhibit nanomolar inhibition by targeting NF-κB and TAK1, EGFR, or COX-2. However, the NF-κB suppressing activity of these compounds is relatively low (IC50 > 300 nM). Compounds 6-14 suppress NF-κB expression with IC50 values ranging from 1.8 nM to 38.6 nM. HS-276 is a highly selective, orally bioavailable TAK1 inhibitor. Combined structural modifications of compounds using a prodrug strategy may enhance NF-κB inhibition. This review focused on the role and mechanism of JAZF1 in inflammation and lipid accumulation for the identification of new anti-atherosclerotic targets.

Juvenile idiopathic arthritis-associated genetic loci exhibit spatially constrained gene regulatory effects across multiple tissues and immune cell types

Juvenile idiopathic arthritis (JIA) is an autoimmune, inflammatory joint disease with complex genetic etiology. Previous GWAS have found many genetic loci associated with JIA. However, the biological mechanism behind JIA remains unknown mainly because most risk loci are located in non-coding genetic regions. Interestingly, increasing evidence has found that regulatory elements in the non-coding regions can regulate the expression of distant target genes through spatial (physical) interactions. Here, we used information on the 3D genome organization (Hi-C data) to identify target genes that physically interact with SNPs within JIA risk loci. Subsequent analysis of these SNP-gene pairs using data from tissue and immune cell type-specific expression quantitative trait loci (eQTL) databases allowed the identification of risk loci that regulate the expression of their target genes. In total, we identified 59 JIA-risk loci that regulate the expression of 210 target genes across diverse tissues and immune cell types. Functional annotation of spatial eQTLs within JIA risk loci identified significant overlap with gene regulatory elements (i.e., enhancers and transcription factor binding sites). We found target genes involved in immune-related pathways such as antigen processing and presentation (e.g., ERAP2, HLA class I and II), the release of pro-inflammatory cytokines (e.g., LTBR, TYK2), proliferation and differentiation of specific immune cell types (e.g., AURKA in Th17 cells), and genes involved in physiological mechanisms related to pathological joint inflammation (e.g., LRG1 in arteries). Notably, many of the tissues where JIA-risk loci act as spatial eQTLs are not classically considered central to JIA pathology. Overall, our findings highlight the potential tissue and immune cell type-specific regulatory changes contributing to JIA pathogenesis. Future integration of our data with clinical studies can contribute to the development of improved JIA therapy.

Characterization of age-associated B cells in early drug-naïve rheumatoid arthritis patients

Age-associated B cells (ABCs) are an immune cell subset linked to autoimmunity, infection and ageing, and whose pathophysiological importance was recently highlighted using single cell synovial tissue profiling. To elucidate their pathophysiological relevance, peripheral blood (PB) ABCs from early rheumatoid arthritis (eRA) patients naïve to disease-modifying anti-rheumatic drugs (DMARDs) were compared with their synovial fluid (SF) counterparts, and to PB ABCs from psoriatic arthritis patients and healthy controls. PB and SF B-cell subsets were phenotyped by multi-parameter flow cytometry, sorted and subjected to gene expression profiling (NanoString nCounter® Immunology V2 Panel) and functional characterization (stimulated cytokine measurements by immunoassay). PB ABCs of eRA patients, which are transcriptionally distinct from those of control cohorts, express chemokine receptors and adhesion molecules, such as CXCR3, that favour homing to inflammatory sites over lymphoid tissue. These cells are an activated, class-switched B-cell subset expressing high levels of HLA-DR, co-stimulatory molecules and T-bet. Their secretion profile includes IL-12p70 and IL-23 but low levels of IL-10. High surface expression of FcRL family members, including FcRL3, furthermore suggests a role for these cells in autoimmunity. Finally, and unlike in the periphery where they are rare, ABCs are the predominant B-cell subsets in SF. These observations indicate the predilection of ABCs for inflammatory tissue in RA, where their propensity for antigen presentation and pro-inflammatory phenotype may support autoimmune pathology. Their potential as a therapeutic target therefore warrants further study.

CD83 Regulates the Immune Responses in Inflammatory Disorders

Activating the immune system plays an important role in maintaining physiological homeostasis and defending the body against harmful infections. However, abnormalities in the immune response can lead to various immunopathological responses and severe inflammation. The activation of dendritic cells (DCs) can influence immunological responses by promoting the differentiation of T cells into various functional subtypes crucial for the eradication of pathogens. CD83 is a molecule known to be expressed on mature DCs, activated B cells, and T cells. Two isotypes of CD83, a membrane-bound form and a soluble form, are subjects of extensive scientific research. It has been suggested that CD83 is not only a ubiquitous co-stimulatory molecule but also a crucial player in monitoring and resolving inflammatory reactions. Although CD83 has been involved in immunological responses, its functions in autoimmune diseases and effects on pathogen immune evasion remain unclear. Herein, we outline current immunological findings and the proposed function of CD83 in inflammatory disorders.

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

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

Methylation quantitative trait loci are largely consistent across disease states in Crohn’s disease

Recently, we identified 1,189 CpG sites whose DNA methylation level in blood associated with Crohn’s disease. Here, we examined associations between DNA methylation and genetic variants to identify methylation quantitative trait loci across disease states in (1) 402 blood samples from 164 newly diagnosed pediatric Crohn’s disease patients taken at 2 time points (diagnosis and follow-up), and 74 non-inflammatory bowel disease controls, (2) 780 blood samples from a non-Crohn’s disease adult population, and (3) 40 ileal biopsies (17 Crohn’s disease cases and 23 non-inflammatory bowel disease controls) from group (1). Genome-wide DNAm profiling and genotyping were performed using the Illumina MethylationEPIC and Illumina Multi-Ethnic arrays. SNP-CpG associations were identified via linear models adjusted for age, sex, disease status, disease subtype, estimated cell proportions, and genotype-based principal components. In total, we observed 535,448 SNP-CpG associations between 287,881 SNPs and 12,843 CpG sites (P < 8.21 × 10−14). Associations were highly consistent across different ages, races, disease states, and tissue types, suggesting that the majority of these methylation quantitative trait loci participate in common gene regulation. However, genes near CpGs associated with inflammatory bowel disease SNPs were enriched for 18 KEGG pathways relevant to inflammatory bowel disease-linked immune function and inflammatory responses. We observed suggestive evidence for a small number of tissue-specific associations and disease-specific associations in ileum, though larger studies will be needed to confirm these results. Our study concludes that the vast majority of blood-derived methylation quantitative trait loci are common across individuals, though a subset may be involved in processes related to Crohn’s disease. Independent cohort studies will be required to validate these findings.

Regulation of activated T cell survival in rheumatic autoimmune diseases

Adaptive immune responses rely on the proliferation of T lymphocytes able to recognize and eliminate pathogens. The magnitude and duration of the expansion of activated T cell clones are finely regulated to minimize immunopathology and avoid autoimmunity. In patients with rheumatic autoimmune diseases, such as systemic lupus erythematosus and rheumatoid arthritis, activated lymphocytes survive and exert effector functions for prolonged periods, defying the mechanisms that normally curb their capacities during acute and chronic infections. Here, we review the molecular mechanisms that limit the duration of immune responses in health and discuss the factors that alter such regulation in the setting of systemic lupus erythematosus and rheumatoid arthritis. We highlight defects that could contribute to the development and progression of autoimmune disease and describe how chronic inflammation can alter the regulation of activated lymphocyte survival, promoting its perpetuation. These concepts might contribute to the understanding of the mechanisms that underlie the chronicity of inflammation in the context of autoimmunity.

Genomic Insights into Myasthenia Gravis Identify Distinct Immunological Mechanisms in Early and Late Onset Disease

OBJECTIVE

The purpose of this study was to identify disease relevant genes and explore underlying immunological mechanisms that contribute to early and late onset forms of myasthenia gravis.

METHODS

We used a novel genomic methodology to integrate genomewide association study (GWAS) findings in myasthenia gravis with cell-type specific information, such as gene expression patterns and promotor-enhancer interactions, in order to identify disease-relevant genes. Subsequently, we conducted additional genomic investigations, including an expression quantitative analysis of 313 healthy people to provide mechanistic insights.

RESULTS

We identified several genes that were specifically linked to early onset myasthenia gravis including TNIP1, ORMDL3, GSDMB, and TRAF3. We showed that regulators of toll-like receptor 4 signaling were enriched among these early onset disease genes (fold enrichment = 3.85, p = 6.4 × 10^-3 ). In contrast, T-cell regulators CD28 and CTLA4 were exclusively linked to late onset disease. We identified 2 causal genetic variants (rs231770 and rs231735; posterior probability = 0.98 and 0.91) near the CTLA4 gene. Subsequently, we demonstrated that these causal variants result in low expression of CTLA4 (rho = -0.66, p = 1.28 × 10^-38 and rho = -0.52, p = 7.01 × 10^-22 , for rs231735 and rs231770, respectively).

INTERPRETATION

The disease-relevant genes identified in this study are a unique resource for many disciplines, including clinicians, scientists, and the pharmaceutical industry. The distinct immunological pathways linked to early and late onset myasthenia gravis carry important implications for drug repurposing opportunities and for future studies of drug development. ANN NEUROL 2021;90:455-463.

METTL21B is a prognostic biomarker and potential therapeutic target in low-grade gliomas

A considerable amount of literature has demonstrated that eukaryotic translation elongation factor 1A (eEF1A) is closely related to tumors. As a newly identified lysine specific methyltransferase targeting eEF1A at Lys-165, too little attention has been paid to the function of METTL21B. To determine the potential significance and prognostic value of METTL21B in low grade glioma (LGG), we analyzed the expression, methylation level and copy number variations (CNV) of METTL21B and its effect on prognosis in patients with LGG by 4 public databases in conjunction with experimental examination of LGG patient samples. As a result, we found that high expression, hypomethylation and gain/amplification of CNV of METTL21B were associated with poor prognosis in LGG. The potential functions of METTL21B in LGG may be involved in cell adhesion, angiogenesis and cell proliferation of tumor by enrichment analysis. In addition, METTL21B may facilitate immune evasion of tumor and affect prognosis by mediating macrophage polarization from M1 to M2 and regulating expression of immune checkpoints. Nevertheless, patients with high METTL21B level are likely to have better response to immune checkpoints blockage therapy. Because of its substrate specificity, METTL21B is expected to be a promising target for the treatment of glioma.

Functional genomics atlas of synovial fibroblasts defining rheumatoid arthritis heritability

BACKGROUND

Genome-wide association studies have reported more than 100 risk loci for rheumatoid arthritis (RA). These loci are shown to be enriched in immune cell-specific enhancers, but the analysis so far has excluded stromal cells, such as synovial fibroblasts (FLS), despite their crucial involvement in the pathogenesis of RA. Here we integrate DNA architecture, 3D chromatin interactions, DNA accessibility, and gene expression in FLS, B cells, and T cells with genetic fine mapping of RA loci.

RESULTS

We identify putative causal variants, enhancers, genes, and cell types for 30-60% of RA loci and demonstrate that FLS account for up to 24% of RA heritability. TNF stimulation of FLS alters the organization of topologically associating domains, chromatin state, and the expression of putative causal genes such as TNFAIP3 and IFNAR1. Several putative causal genes constitute RA-relevant functional networks in FLS with roles in cellular proliferation and activation. Finally, we demonstrate that risk variants can have joint-specific effects on target gene expression in RA FLS, which may contribute to the development of the characteristic pattern of joint involvement in RA.

CONCLUSION

Overall, our research provides the first direct evidence for a causal role of FLS in the genetic susceptibility for RA accounting for up to a quarter of RA heritability.

Genome editing to define the function of risk loci and variants in rheumatic disease

Discoveries in human genetic studies have revolutionized our understanding of complex rheumatic and autoimmune diseases, including the identification of hundreds of genetic loci and single nucleotide polymorphisms that potentially predispose individuals to disease. However, in most cases, the exact disease-causing variants and their mechanisms of action remain unresolved. Functional follow-up of these findings is most challenging for genomic variants that are in non-coding genomic regions, where the large majority of common disease-associated variants are located, and/or that probably affect disease progression via cell type-specific gene regulation. To deliver on the therapeutic promise of human genetic studies, defining the mechanisms of action of these alleles is essential. Genome editing technology, such as CRISPR-Cas, has created a vast toolbox for targeted genetic and epigenetic modifications that presents unprecedented opportunities to decipher disease-causing loci, genes and variants in autoimmunity. In this Review, we discuss the past 5-10 years of progress in resolving the mechanisms underlying rheumatic disease-associated alleles, with an emphasis on how genomic editing techniques can enable targeted dissection and mechanistic studies of causal autoimmune risk variants.

Expression Quantitative Trait Locus Analysis in Systemic Sclerosis Identifies New Candidate Genes Associated With Multiple Aspects of Disease Pathology

OBJECTIVE

To identify the genetic variants that affect gene expression (expression quantitative trait loci [eQTLs]) in systemic sclerosis (SSc) and to investigate their role in the pathogenesis of the disease.

METHODS

We performed an eQTL analysis using whole-blood sequencing data from 333 SSc patients and 524 controls and integrated them with SSc genome-wide association study (GWAS) data. We integrated our findings from expression modeling, differential expression analysis, and transcription factor binding site enrichment with key clinical features of SSc.

RESULTS

We detected 49,123 validated cis-eQTLs from 4,539 SSc-associated single-nucleotide polymorphisms (SNPs) (P_GWAS < 10^-5 ). A total of 1,436 genes were within 1 Mb of the 4,539 SSc-associated SNPs. Of those 1,436 genes, 565 were detected as having ≥1 eQTL with an SSc-associated SNP. We developed a strategy to prioritize disease-associated genes based on their expression variance explained by SSc eQTLs (r² > 0.05). As a result, 233 candidates were identified, 134 (58%) of them associated with hallmarks of SSc and 105 (45%) of them differentially expressed in the blood cells, skin, or lung tissue of SSc patients. Transcription factor binding site analysis revealed enriched motifs of 24 transcription factors (5%) among SSc eQTLs, 5 of which were found to be differentially regulated in the blood cells (ELF1 and MGA), skin (KLF4 and ID4), and lungs (TBX4) of SSc patients. Ten candidate genes (4%) can be targeted by approved medications for immune-mediated diseases, of which only 3 have been tested in clinical trials in patients with SSc.

CONCLUSION

The findings of the present study indicate a new layer to the molecular complexity of SSc, contributing to a better understanding of the pathogenesis of the disease.

Joo JW. MARS: leveraging allelic heterogeneity to increase power of association testing

In standard genome-wide association studies (GWAS), the standard association test is underpowered to detect associations between loci with multiple causal variants with small effect sizes. We propose a statistical method, Model-based Association test Reflecting causal Status (MARS), that finds associations between variants in risk loci and a phenotype, considering the causal status of variants, only requiring the existing summary statistics to detect associated risk loci. Utilizing extensive simulated data and real data, we show that MARS increases the power of detecting true associated risk loci compared to previous approaches that consider multiple variants, while controlling the type I error.

Lymphocytes sense antibodies through human FCRL proteins: Emerging roles in mucosal immunity

Members of the Fc receptor-like (FCRL) family modulate B and T cell responses, yet their functional roles remain enigmatic. Nevertheless, FCRL3 promoter polymorphism that alters gene expression has been associated with autoimmune disease risk, indicating physiologic importance. Providing essential functional context, human FCRL3, FCRL4, and FCRL5 have recently been identified as secretory IgA (SIgA), dimeric IgA, and IgG receptors, respectively, revealing novel ways lymphocytes can interact with antibodies. FCRL3 and FCRL4 are able to distinguish the mucosal and systemic origin of IgA-containing immune complexes, respectively, with clear implications in guiding mucosal responses. SIgA can signal mucosal breach through FCRL3, driving the functional plasticity of regulatory T cells toward inflammatory to help control invading pathogens. Conversely, recognition of dimeric IgA by FCRL4 on memory B cells located in mucosa-associated lymphoid tissues could promote tolerance to commensals. Memory B cells that accumulate under conditions of chronic antigen presence frequently express FCRL4 and FCRL5, and antibody ligands could provide functional feedback to the cells. FCRL5 apparently recognizes the age of the IgG molecule, using deamidation as a molecular clock, conceivably playing regulatory roles in chronic antibody responses. A framework of FCRL3, FCRL4, and FCRL5 operating as sensors of antibodies in immune complexes is proposed. Sensing the spatial origin and age of immune complexes can shape lymphocyte functional attributes and inform their participation in mucosal immune responses. The potential contributions of FCRL3 and SIgA to the pathogenesis of autoimmune diseases are discussed.

Using functional genomics to advance the understanding of psoriatic arthritis

Psoriatic arthritis (PsA) is a complex disease where susceptibility is determined by genetic and environmental risk factors. Clinically, PsA involves inflammation of the joints and the skin, and, if left untreated, results in irreversible joint damage. There is currently no cure and the few treatments available to alleviate symptoms do not work in all patients. Over the past decade, genome-wide association studies (GWAS) have uncovered a large number of disease-associated loci but translating these findings into functional mechanisms and novel targets for therapeutic use is not straightforward. Most variants have been predicted to affect primarily long-range regulatory regions such as enhancers. There is now compelling evidence to support the use of chromatin conformation analysis methods to discover novel genes that can be affected by disease-associated variants. Here, we will review the studies published in the field that have given us a novel understanding of gene regulation in the context of functional genomics and how this relates to the study of PsA and its underlying disease mechanism.

ORMDL3 and Asthma: Linking Sphingolipid Regulation to Altered T Cell Function

Orosomucoid like 3 (ORMDL3) encodes an ER-resident transmembrane protein that regulates the activity of serine palmitoyltransferase (SPT), the first and rate-limiting enzyme for sphingolipid biosynthesis in cells. A decade ago, several genome wide association studies revealed single nucleotide polymorphisms associated with increased ORMDL3 protein expression and susceptibility to allergic asthma. Since that time, numerous studies have investigated how altered ORMDL3 expression might predispose to asthma and other autoimmune/inflammatory diseases. In this brief review, we focus on growing evidence suggesting that heightened ORMDL3 expression specifically in CD4⁺ T lymphocytes, the central orchestrators of adaptive immunity, constitutes a major underlying mechanism of asthma pathogenesis by skewing their differentiation and function. Furthermore, we explore how sphingolipid modulation in T cells might be responsible for these effects, and how further studies may interrogate this intriguing hypothesis.

Two genetic variants explain the association of European ancestry with multiple sclerosis risk in African-Americans

Epidemiological studies have suggested differences in the rate of multiple sclerosis (MS) in individuals of European ancestry compared to African ancestry, motivating genetic scans to identify variants that could contribute to such patterns. In a whole-genome scan in 899 African-American cases and 1155 African-American controls, we confirm that African-Americans who inherit segments of the genome of European ancestry at a chromosome 1 locus are at increased risk for MS [logarithm of odds (LOD) = 9.8], although the signal weakens when adding an additional 406 cases, reflecting heterogeneity in the two sets of cases [logarithm of odds (LOD) = 2.7]. The association in the 899 individuals can be fully explained by two variants previously associated with MS in European ancestry individuals. These variants tag a MS susceptibility haplotype associated with decreased CD58 gene expression (odds ratio of 1.37; frequency of 84% in Europeans and 22% in West Africans for the tagging variant) as well as another haplotype near the FCRL3 gene (odds ratio of 1.07; frequency of 49% in Europeans and 8% in West Africans). Controlling for all other genetic and environmental factors, the two variants predict a 1.44-fold higher rate of MS in European-Americans compared to African-Americans.

Emerging insights on the role of gasdermins in infection and inflammatory diseases

The gasdermins, family of pore-forming proteins, are emerging key regulators of infection, autoinflammation and antitumor immunity. Multiple studies have recently characterised their crucial roles in driving pyroptosis, a lytic pro-inflammatory type of cell death. Additionally, gasdermins also act as key effectors of NETosis, secondary necrosis and apoptosis. In this review, we will address current understanding of the mechanisms of gasdermin activation and further describe the protective and detrimental roles of gasdermins in host defence and autoinflammatory diseases. These data suggest that gasdermins play a prominent role in innate immunity and autoinflammatory disorders, thereby providing potential new therapeutic avenues for the treatment of infection and autoimmune disease.

Functional genomics in autoimmune diseases

Associations between genetic loci and increased susceptibility to autoimmune disease have been well characterized, however, translating this knowledge into mechanistic insight and patient benefit remains a challenge. While improvements in the precision, completeness and accuracy of our genetic understanding of autoimmune diseases will undoubtedly be helpful, meeting this challenge will require two interlinked problems to be addressed: first which of the highly correlated variants at an individual locus is responsible for increased disease risk, and second what are the downstream effects of this variant. Given that the majority of loci are thought to affect non-coding regulatory elements, the second question is often reframed as what are the target gene(s) and pathways affected by causal variants. Currently, these questions are being addressed using a wide variety of novel techniques and datasets. In many cases, these approaches are complementary and it is likely that the most accurate picture will be generated by consolidating information relating to transcription, regulatory activity, chromatin accessibility, chromatin conformation and readouts from functional experiments, such as genome editing and reporter assays. It is clear that it will be necessary to gather this information from disease relevant cell types and conditions and that by doing so our understanding of disease etiology will be improved. This review is focused on the field of autoimmune disease functional genomics with a particular focus on the most exciting and significant research to be published within the last couple of years.

The importance of differences; On environment and its interactions with genes and immunity in the causation of rheumatoid arthritis

The current review uses rheumatoid arthritis (RA) as a prominent example for how studies on the interplay between environmental and genetic factors in defined subsets of a disease can be used to formulate aetiological hypotheses that subsequently can be tested for causality using molecular and functional studies. Major discussed findings are that exposures to airways from many different noxious agents including cigarette smoke, silica dust and more interact with major susceptibility genes, mainly HLA-DR genetic variants in triggering antigen-specific immune reactions specific for RA. We also discuss how several other environmental and lifestyle factors, including microbial, neural and metabolic factors, can influence risk for RA in ways that are different in different subsets of RA.The description of these processes in RA provides the best example so far in any immune-mediated disease of how triggering of immunity at one anatomical site in the context of known environmental and genetic factors subsequently can lead to symptoms that precede the classical inflammatory disease symptoms and later contribute also to the classical RA joint inflammation. The findings referred to in the review have led to a change of paradigms for very early therapy and prevention of RA and to efforts towards what we have named 'personalized prevention'. We believe that the progress described here for RA will be of relevance for research and practice also in other immune-mediated diseases.

The CD83 Molecule - An Important Immune Checkpoint

The CD83 molecule has been identified to be expressed on numerous activated immune cells, including B and T lymphocytes, monocytes, dendritic cells, microglia, and neutrophils. Both isoforms of CD83, the membrane-bound as well as its soluble form are topic of intensive research investigations. Several studies revealed that CD83 is not a typical co-stimulatory molecule, but rather plays a critical role in controlling and resolving immune responses. Moreover, CD83 is an essential factor during the differentiation of T and B lymphocytes, and the development and maintenance of tolerance. The identification of its interaction partners as well as signaling pathways have been an enigma for the last decades. Here, we report the latest data on the expression, structure, and the signaling partners of CD83. In addition, we review the regulatory functions of CD83, including its striking modulatory potential to maintain the balance between tolerance versus inflammation during homeostasis or pathologies. These immunomodulatory properties of CD83 emphasize its exceptional therapeutic potential, which has been documented in specific preclinical disease models.

Lymphocyte DNA methylation mediates genetic risk at shared immune-mediated disease loci

Background

Defining regulatory mechanisms through which noncoding risk variants influence the cell-mediated pathogenesis of immune-mediated disease (IMD) has emerged as a priority in the post–genome-wide association study era.

Objectives

With a focus on rheumatoid arthritis, we sought new insight into genetic mechanisms of adaptive immune dysregulation to help prioritize molecular pathways for targeting in this and related immune pathologies.

Methods

Whole-genome methylation and transcriptional data from isolated CD4⁺ T cells and B cells of more than 100 genotyped and phenotyped patients with inflammatory arthritis, all of whom were naive to immunomodulatory treatments, were obtained. Analysis integrated these comprehensive data with genome-wide association study findings across IMDs and other publicly available resources.

Results

We provide strong evidence that disease-associated DNA variants regulate cis-CpG methylation in CD4⁺ T and/or B cells at 37% RA loci. Using paired, cell-specific transcriptomic data and causal inference testing, we identify examples where site-specific DNA methylation in turn mediates gene expression, including FCRL3 in both cell types and ORMDL3/GSDMB, IL6ST/ANKRD55, and JAZF1 in CD4⁺ T cells. A number of genes regulated in this way highlight mechanisms common to RA and other IMDs including multiple sclerosis and asthma, in turn distinguishing them from osteoarthritis, a primarily degenerative disease. Finally, we corroborate the observed effects experimentally.

Conclusions

Our observations highlight important mechanisms of genetic risk in RA and the wider context of immune dysregulation. They confirm the utility of DNA methylation profiling as a tool for causal gene prioritization and, potentially, therapeutic targeting in complex IMD.

Exploring inflammatory signatures in arthritic joint biopsies with Spatial Transcriptomics

Lately it has become possible to analyze transcriptomic profiles in tissue sections with retained cellular context. We aimed to explore synovial biopsies from rheumatoid arthritis (RA) and spondyloarthritis (SpA) patients, using Spatial Transcriptomics (ST) as a proof of principle approach for unbiased mRNA studies at the site of inflammation in these chronic inflammatory diseases. Synovial tissue biopsies from affected joints were studied with ST. The transcriptome data was subjected to differential gene expression analysis (DEA), pathway analysis, immune cell type identification using Xcell analysis and validation with immunohistochemistry (IHC). The ST technology allows selective analyses on areas of interest, thus we analyzed morphologically distinct areas of mononuclear cell infiltrates. The top differentially expressed genes revealed an adaptive immune response profile and T-B cell interactions in RA, while in SpA, the profiles implicate functions associated with tissue repair. With spatially resolved gene expression data, overlaid on high-resolution histological images, we digitally portrayed pre-selected cell types in silico. The RA displayed an overrepresentation of central memory T cells, while in SpA effector memory T cells were most prominent. Consequently, ST allows for deeper understanding of cellular mechanisms and diversity in tissues from chronic inflammatory diseases.

Identification of rare coding variants in TYK2 protective for rheumatoid arthritis in the Japanese population and their effects on cytokine signalling

OBJECTIVE

Although genome-wide association studies (GWAS) have identified approximately 100 loci for rheumatoid arthritis (RA), the disease mechanisms are not completely understood. We evaluated the pathogenesis of RA by focusing on rare coding variants.

METHODS

The coding regions of 98 candidate genes identified by GWAS were sequenced in 2294 patients with RA and 4461 controls in Japan. An association analysis was performed using cases and controls for variants, genes and domains of TYK2. Cytokine responses for two associated variants (R231W, rs201917359; and R703W, rs55882956) in TYK2 as well as a previously reported risk variant (P1004A, rs34536443) for multiple autoimmune diseases were evaluated by reporter assays.

RESULTS

A variant in TYK2 (R703W) showed a suggestive association (p=5.47×10-8, OR=0.48). We observed more accumulation of rare coding variants in controls in TYK2 (p=3.94×10-12, OR=0.56). The four-point-one, ezrin, radixin, moesin (FERM; 2.14×10-3, OR=0.66) and pseudokinase domains (1.63×10-8, OR=0.52) of TYK2 also showed enrichment of variants in controls. R231W in FERM domain especially reduced interleukin (IL)-6 and interferon (IFN)-γ signalling, whereas P1104A in kinase domain reduced IL-12, IL-23 and IFN-α signalling. R703W in pseudokinase domain reduced cytokine signals similarly to P1104A, but the effects were weaker than those of P1104A.

CONCLUSIONS

The FERM and pseudokinase domains in TYK2 were associated with the risk of RA in the Japanese population. Variants in TYK2 had different effects on cytokine signalling, suggesting that the regulation of selective cytokine signalling is a target for RA treatment.

A Combined Transcriptomic and Genomic Analysis Identifies a Gene Signature Associated With the Response to Anti-TNF Therapy in Rheumatoid Arthritis

Background: Rheumatoid arthritis (RA) is the most frequent autoimmune disease involving the joints. Although anti-TNF therapies have proven effective in the management of RA, approximately one third of patients do not show a significant clinical response. The objective of this study was to identify new genetic variation associated with the clinical response to anti-TNF therapy in RA.

Methods: We performed a sequential multi-omic analysis integrating different sources of molecular information. First, we extracted the RNA from synovial biopsies of 11 RA patients starting anti-TNF therapy to identify gene coexpression modules (GCMs) in the RA synovium. Second, we analyzed the transcriptomic association between each GCM and the clinical response to anti-TNF therapy. The clinical response was determined at week 14 using the EULAR criteria. Third, we analyzed the association between the GCMs and anti-TNF response at the genetic level. For this objective, we used genome-wide data from a cohort of 348 anti-TNF treated patients from Spain. The GCMs that were significantly associated with the anti-TNF response were then tested for validation in an independent cohort of 2,706 anti-TNF treated patients. Finally, the functional implication of the validated GCMs was evaluated via pathway and cell type epigenetic enrichment analyses.

Results: A total of 149 GCMs were identified in the RA synovium. From these, 13 GCMs were found to be significantly associated with anti-TNF response (P < 0.05). At the genetic level, we detected two of the 13 GCMs to be significantly associated with the response to adalimumab (P = 0.0015) and infliximab (P = 0.021) in the Spain cohort. Using the independent cohort of RA patients, we replicated the association of the GCM associated with the response to adalimumab (P = 0.0019). The validated module was found to be significantly enriched for genes involved in the nucleotide metabolism (P = 2.41e-5) and epigenetic marks from immune cells, including CD4+ regulatory T cells (P = 0.041).

Conclusions: These findings show the existence of a drug-specific genetic basis for anti-TNF response, thereby supporting treatment stratification in the search for response biomarkers in RA.

Expression of STAT3-regulated genes in circulating CD4+ T cells discriminates rheumatoid arthritis independently of clinical parameters in early arthritis

OBJECTIVES

Dysregulated signal transduction and activator of transcription-3 (STAT3) signalling in CD4+ T cells has been proposed as an early pathophysiological event in RA. We sought further evidence for this observation, and to determine its clinical relevance.

METHODS

Microarray technology was used to measure gene expression in purified peripheral blood CD4+ T cells from treatment-naïve RA patients and disease controls newly recruited from an early arthritis clinic. Analysis focused on 12 previously proposed transcripts, and concurrent STAT3 pathway activation was determined in the same cells by flow cytometry. A pooled analysis of previous and current gene expression findings incorporated detailed clinical parameters and employed multivariate analysis.

RESULTS

In an independent cohort of 161 patients, expression of 11 of 12 proposed signature genes differed significantly between RA patients and controls, robustly validating the earlier findings. Differential regulation was most pronounced for the STAT3 target genes PIM1, BCL3 and SOCS3 (>1.3-fold difference; P < 0.005), each of whose expression correlated strongly with paired intracellular phospho-STAT3. In a meta-analysis of 279 patients the same three genes accounted for the majority of the signature's ability to discriminate RA patients, which was found to be independent of age, joint involvement or acute phase response.

CONCLUSION

The STAT3-mediated dysregulation of BCL3, SOCS3 and PIM1 in circulating CD4+ T cells is a discriminatory feature of early RA that occurs independently of acute phase response. The mechanistic and functional implications of this observation at a cellular level warrant clarification.

Bioinformatics Analysis and Identification of Genes and Molecular Pathways Involved in Synovial Inflammation in Rheumatoid Arthritis

BACKGROUND Rheumatoid arthritis (RA) has a high prevalence in the elderly population. The genes and pathways in the inflamed synovium in patients with RA are poorly understood. This study aimed to identify differentially expressed genes (DEGs) linked to the progression of synovial inflammation in RA using bioinformatics analysis. MATERIAL AND METHODS Gene expression profiles of datasets GSE55235 and GSE55457 were acquired from the Gene Expression Omnibus (GEO) database. DEGs were identified using Morpheus software, and co-expressed DEGs were identified with Venn diagrams. Protein-protein interaction (PPI) networks were assembled with Cytoscape software and separated into subnetworks using the Molecular Complex Detection (MCODE) algorithm. The functions of the top module were assessed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID). The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed. RESULTS DEGs that were upregulated were significantly enhanced in protein binding, the cell cytosol, organization of the extracellular matrix (ECM), regulation of RNA transcription, and cell adhesion. DEGs that were downregulated were associated with control of the immune response, B-cell and T-cell receptor signaling pathway regulation. KEGG pathway analysis showed that upregulated DEGs enhanced pathways associated with the cell adherens junction, osteoclast differentiation, and hereditary cardiomyopathies. Downregulated DEGs were enriched in primary immunodeficiency, cell adhesion molecules (CAMs), cytokine-cytokine receptor interaction, and hematopoietic cell lineages. CONCLUSIONS The findings from this bioinformatics network analysis study identified molecular mechanisms and the key hub genes that may contribute to synovial inflammation in patients with RA.

Identification of rheumatoid arthritis causal genes using functional genomics

Over the past decade, genome-wide association studies have contributed a wealth of knowledge to our understanding of polygenic disorders such as rheumatoid arthritis. As the size of sample cohorts has improved so too have the computational and experimental methods used to robustly define variants associated with disease susceptibility. The challenge now remains to translate these findings into improved understanding of disease aetiology and patient care. Whilst much of the focus of translating the findings of genome-wide association studies has been on global analysis of all variants identified, careful functional study of individual disease susceptibility loci will be required in order to refine our understanding of how individual variants contribute to disease risk. Here, we present the argument behind such an approach and describe some of the novel tools being used to investigate risk loci. This includes the use of chromosomal conformation capture techniques and modifications of the CRISPR-Cas9 system, with several examples of their implementation being described.

Genome-Wide Expression Quantitative Trait Loci Analysis Using Mixed Models

Expression quantitative trait loci (eQTLs) are important for understanding the genetic basis of cellular activities and complex phenotypes. Genome-wide eQTL analyses can be effectively conducted by employing a mixed model. The mixed model includes random polygenic effects with variability, which can be estimated by the covariance structure of pairwise genomic similarity among individuals based on genotype information for nucleotide sequence variants. This increases the accuracy of identifying eQTLs by avoiding population stratification. Its extensive use will accelerate our understanding of the genetics of gene expression and complex phenotypes. An overview of genome-wide eQTL analyses using mixed model methodology is provided, including discussions of both theoretical and practical issues. The advantages of employing mixed models are also discussed in this review.

The rheumatic disease-associated FAM167A-BLK locus encodes DIORA-1, a novel disordered protein expressed highly in bronchial epithelium and alveolar macrophages

Triggering of autoimmunity that leads to rheumatic disease has been suggested to depend upon gene-environment interactions occurring in epithelial barriers and associated immune cells. Genetic studies have identified associations of the FAM167A-BLK locus with rheumatoid arthritis, systemic lupus erythematosus (SLE) and Sjögren's syndrome. While BLK (B lymphocyte kinase) has a well-established role in B cells, family with sequence similarity to 167 member A (FAM167A) and its gene family remain uncharacterized. To begin to understand the role of FAM167A in rheumatic disease pathogenesis, we explored this gene family and cloned and investigated the gene products. Expression of quantitative trait locus analysis was performed in immune cells. FAM167A and FAM167B were cloned from human peripheral blood mononuclear cells (PBMC). Gene conservation and protein properties were analysed by online tools, mRNA expression measured in mouse organs by quantitative polymerase chain reaction (qPCR) and protein expression investigated in human tissues by immunohistochemistry. We found that autoimmune risk genotypes within the FAM167A-BLK locus lead to increased expression of FAM167A. The FAM167 gene family includes two members, FAM167A and FAM167B, which are not homologous to any other annotated gene but are evolutionarily conserved. The encoded proteins, which we denote 'disordered autoimmunity' (DIORA)-1 and DIORA-2, respectively, are characterized by a high content of intrinsic disorder. Notably, DIORA-1 has its highest expression in the lung, detectable in both bronchial epithelium and alveolar macrophages with an endosomal localization pattern. In summary, the FAM167A gene is associated with several rheumatic diseases and encodes a novel disordered protein, DIORA-1, which is expressed highly in the lung, consistent with a potential role in disease pathogenesis.