Genetic influences on susceptibility to rheumatoid arthritis in African-Americans

Dissolving microneedles: A transdermal drug delivery system for the treatment of rheumatoid arthritis

Rheumatoid arthritis (RA) is an autoimmune disorder that impacts around 1% of the global population. Up to 20% of people become disabled within a year, which has a severely negative impact on their health and quality of life. RA has a complicated pathogenic mechanism, which initially affects small joints and progresses to larger ones over time. It can damage the skin, eyes, heart, kidney, and lung. Oral medications, intra-articular injections, and other treatments are being used; nevertheless, they have drawbacks, including low bioavailability, numerous adverse effects, and poor patient compliance. Dissolving microneedles (DMNs) are a safe and painless method for transdermal drug delivery, achieved through their ability to physically penetrate the epidermal barrier. They enable targeted drug delivery, significantly enhancing the bioavailability of medications and improving patient compliance. DMNs are particularly effective in delivering both lipophilic and high molecular weight biomolecules. The superior bioavailability of DMNs is demonstrated by the fact that low-dose DMN administration can achieve up to 25.8 times higher bioavailability compared to oral administration. This paper provides a comprehensive review of recent advancements in the use of DMNs for RA treatment, encompassing various materials (such as hyaluronic acid, chitosan, etc.), fabrication techniques (such as the two-step casting method, photopolymerization), and performance evaluations (including morphology, mechanical properties, skin penetration capability, solubility, and pharmacodynamics). Additionally, a thorough safety assessment has been conducted, revealing that DMNs cause minimal skin irritation and exhibit low cytotoxicity, ensuring their safety for clinical application. DMNs provide a highly effective and promising alternative to oral and injectable drug delivery systems, offering a novel therapeutic approach for RA patients that significantly improves treatment outcomes and enhances their quality of life.

Prevention of rheumatoid arthritis using a familial predictive medicine approach

Most of the chronic-degenerative diseases deserve a very early recognition of symptoms and signs for the earliest secondary prevention, which could be also very useful in many cases for the most precocious clinical approach. The periodic monitoring of a subject at risk of a specific disease, because of genomic predisposition by predictive medicine approach, may help to earlier detection of onset and/or the progression of the pathology itself, through intra-individual monitoring. This is particularly the case of rheumatoid arthritis (RA) for which an early diagnosis is undoubtedly the first step to ensure the most proper therapy for the patient. Thus, the earlier identification of individuals at high risk of RA could lead to ultra-preventive strategies to start for the best lifestyle performances and/or for any other effective therapeutic interventions to contrast the onset, and/or the evolution of the putative RA. This will also optimize both costs and medical resources, according to the health care policies of many countries.

The Veterans Affairs Rheumatoid Arthritis Registry: A unique population in rheumatoid arthritis research

BACKGROUND

As the largest integrated healthcare system in the U.S., the Veterans Affairs (VA) provides a unique context for the conduct of clinical and clinical-translational research in rheumatoid arthritis (RA).

OBJECTIVES

To review attributes of the VA Rheumatoid Arthritis Registry (RA) and highlight its research contributions.

FINDINGS

With >3,600 participants enrolled from 19 VA medical centers across the U.S., VARA includes longitudinally collected clinical data and a central biorepository that includes serum, plasma, and DNA collected at enrollment. VARA research capacity is enhanced via active linkages with internal data including the VA's Corporate Data Warehouse and elements captured during oncology care. This capacity is further enabled via active linkages with the National Death Index and Centers for Medicare & Medicaid Services (CMS) data.

CONCLUSION

As a highly unique study population with comprehensive data annotation available to researchers, VARA is poised to continue address impactful questions in RA for years to come.

Structural variants in linkage disequilibrium with GWAS-significant SNPs

With the recent expansion of structural variant identification in the human genome, understanding the role of these impactful variants in disease architecture is critically important. Currently, a large proportion of genome-wide-significant genome-wide association study (GWAS) single nucleotide polymorphisms (SNPs) are functionally unresolved, raising the possibility that some of these SNPs are associated with disease through linkage disequilibrium with causal structural variants. Hence, understanding the linkage disequilibrium between newly discovered structural variants and statistically significant SNPs may provide a resource for further investigation into disease-associated regions in the genome. Here we present a resource cataloging structural variant-significant SNP pairs in high linkage disequilibrium. The database is composed of (i) SNPs that have exhibited genome-wide significant association with traits, primarily disease phenotypes, (ii) newly released structural variants (SVs), and (iii) linkage disequilibrium values calculated from unphased data. All data files including those detailing SV and GWAS SNP associations and results of GWAS-SNP-SV pairs are available at the SV-SNP LD Database and can be accessed at 'https://github.com/hliang-SchrodiLab/SV_SNPs. Our analysis results represent a useful fine mapping tool for interrogating SVs in linkage disequilibrium with disease-associated SNPs. We anticipate that this resource may play an important role in subsequent studies which investigate incorporating disease causing SVs into disease risk prediction models.

Massively parallel reporter assay confirms regulatory potential of hQTLs and reveals important variants in lupus and other autoimmune diseases

We designed a massively parallel reporter assay (MPRA) in an Epstein-Barr virus transformed B cell line to directly characterize the potential for histone post-translational modifications, i.e., histone quantitative trait loci (hQTLs), expression QTLs (eQTLs), and variants on systemic lupus erythematosus (SLE) and autoimmune (AI) disease risk haplotypes to modulate regulatory activity in an allele-dependent manner. Our study demonstrates that hQTLs, as a group, are more likely to modulate regulatory activity in an MPRA compared with other variant classes tested, including a set of eQTLs previously shown to interact with hQTLs and tested AI risk variants. In addition, we nominate 17 variants (including 11 previously unreported) as putative causal variants for SLE and another 14 for various other AI diseases, prioritizing these variants for future functional studies in primary and immortalized B cells. Thus, we uncover important insights into the mechanistic relationships among genotype, epigenetics, and gene expression in SLE and AI disease phenotypes.

Eomes-dependent mitochondrial regulation promotes survival of pathogenic CD4+ T cells during inflammation

The mechanisms whereby Eomes controls tissue accumulation of T cells and strengthens inflammation remain ill-defined. Here, we show that Eomes deletion in antigen-specific CD4+ T cells is sufficient to protect against central nervous system (CNS) inflammation. While Eomes is dispensable for the initial priming of CD4+ T cells, it is required for long-term maintenance of CNS-infiltrating CD4+ T cells. We reveal that the impact of Eomes on effector CD4+ T cell longevity is associated with sustained expression of multiple genes involved in mitochondrial organization and functions. Accordingly, epigenetic studies demonstrate that Eomes supports mitochondrial function by direct binding to either metabolism-associated genes or mitochondrial transcriptional modulators. Besides, the significance of these findings was confirmed in CD4+ T cells from healthy donors and multiple sclerosis patients. Together, our data reveal a new mechanism by which Eomes promotes severity and chronicity of inflammation via the enhancement of CD4+ T cell mitochondrial functions and resistance to stress-induced cell death.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

Analysis of the Relationship between Rheumatoid Arthritis and Osteoporosis Based on Mendelian Randomization

BACKGROUND

As the global population ages, the World Health Organization has found a yearly increase in the incidence of rheumatoid arthritis and osteoporosis. This trend poses a challenge to public health and healthcare and calls for the implementation of more preventive and treatment measures to address these health issues.

OBJECTIVE

This study aims to investigate the causal relationship between rheumatoid arthritis (RA) and osteoporosis (OP) using the Mendelian randomization (MR) method.

METHODS

OP diagnosis was based on the gold standard of bone mineral density (BMD). Single nucleotide polymorphisms (SNPs) were identified from the genome-wide association research database formed by RA and BMD, with a parameter setting of P < 5×10-8, chain imbalance r2<0.01, and kb = 10,000. Five complementary MR methods, including inverse variance weighted (IVW), MR-Egger regression, weighted median, simple mode estimation based on mode, and weighted estimation based on mode, were used to evaluate the causal relationship between RA and OP/BMD using odds ratio (OR) values and 95% confidence intervals (CI). Sensitivity analyses were performed using heterogeneity tests, horizontal pleiotropy, and individual rejection tests.

RESULTS

A total of 78 instrumental variables were identified that were closely related to both RA and BMD in mixed populations, while 14 instrumental variables were identified in the European population and 38 instrumental variables were identified in the Asian population. Using IVW as the main analysis method, the MR analysis results of RA and BMD showed the following: mixed population OR = 0.96, 95%CI: 0.93-1.00; European population OR = 0.55, 95%CI: 0.27-1.12; and Asian population OR = 0.95, 95%CI: 0.90-1.01. Sensitivity analyses showed that the MR results were robust.

CONCLUSION

The study found insufficient evidence of a causal relationship between RA and OP/BMD, suggesting that RA may not have a direct effect on OP/BMD.

An update on the roles of transcription factor Ets1 in autoimmune diseases

Transcription factors are crucial to regulate gene expression in immune cells and in other cell types. In lymphocytes, there are a large number of different transcription factors that are known to contribute to cell differentiation and the balance between quiescence and activation. One such transcription factor is E26 oncogene homolog 1 (Ets1). Ets1 expression is high in quiescent B and T lymphocytes and its levels are decreased upon activation. The human ETS1 gene has been identified as a susceptibility locus for many autoimmune and inflammatory diseases. In accord with this, gene knockout of Ets1 in mice leads to development of a lupus-like autoimmune disease, with enhanced activation and differentiation of both B cells and T cells. Prior reviews have summarized functional roles for Ets1 based on studies of Ets1 knockout mice. In recent years, numerous additional studies have been published that further validate ETS1 as a susceptibility locus for human diseases where immune dysregulation plays a causative role. In this update, new information that further links Ets1 to human autoimmune diseases is organized and collated to serve as a resource. This update also describes recent studies that seek to understand molecularly how Ets1 regulates immune cell activation, either using human cells and tissues or mouse models. This resource is expected to be useful to investigators seeking to understand how Ets1 may regulate the human immune response, particularly in terms of its roles in autoimmunity and inflammation. This article is categorized under: Immune System Diseases > Genetics/Genomics/Epigenetics Immune System Diseases > Molecular and Cellular Physiology.

Massively Parallel Reporter Assay Confirms Regulatory Potential of hQTLs and Reveals Important Variants in Lupus and Other Autoimmune Diseases

Objective

To systematically characterize the potential for histone post-translational modifications, i.e., histone quantitative trait loci (hQTLs), expression QTLs (eQTLs), and variants on systemic lupus erythematosus (SLE) and autoimmune (AI) disease risk haplotypes to modulate gene expression in an allele dependent manner.

Methods

We designed a massively parallel reporter assay (MPRA) containing ~32K variants and transfected it into an Epstein-Barr virus transformed B cell line generated from an SLE case.

Results

Our study expands our understanding of hQTLs, illustrating that epigenetic QTLs are more likely to contribute to functional mechanisms than eQTLs and other variant types, and a large proportion of hQTLs overlap transcription start sites (TSS) of noncoding RNAs. In addition, we nominate 17 variants (including 11 novel) as putative causal variants for SLE and another 14 for various other AI diseases, prioritizing these variants for future functional studies primary and immortalized B cells.

Conclusion

We uncover important insights into the mechanistic relationships between genotype, epigenetics, gene expression, and SLE and AI disease phenotypes.

Modulation of the immunity and inflammation by autophagy

Autophagy, a highly conserved cellular self-degradation pathway, has emerged with novel roles in the realms of immunity and inflammation. Genome-wide association studies have unveiled a correlation between genetic variations in autophagy-related genes and heightened susceptibility to autoimmune and inflammatory diseases. Subsequently, substantial progress has been made in unraveling the intricate involvement of autophagy in immunity and inflammation through functional studies. The autophagy pathway plays a crucial role in both innate and adaptive immunity, encompassing various key functions such as pathogen clearance, antigen processing and presentation, cytokine production, and lymphocyte differentiation and survival. Recent research has identified novel approaches in which the autophagy pathway and its associated proteins modulate the immune response, including noncanonical autophagy. This review provides an overview of the latest advancements in understanding the regulation of immunity and inflammation through autophagy. It summarizes the genetic associations between variants in autophagy-related genes and a range of autoimmune and inflammatory diseases, while also examining studies utilizing transgenic animal models to uncover the in vivo functions of autophagy. Furthermore, the review delves into the mechanisms by which autophagy dysregulation contributes to the development of three common autoimmune and inflammatory diseases and highlights the potential for autophagy-targeted therapies.

Single Nucleotide Polymorphisms Associated with Rheumatoid Arthritis in Saudi Patients

Rheumatoid arthritis (RA) is a complex, multifactorial disorder with an autoimmune etiology. RA is highly heritable and is associated with both human leucocyte antigen (HLA) and non-HLA genes. We investigated the associations of 33 single nucleotide polymorphisms (SNPs) with RA in the Saudi population.

METHODS

This study included 105 patients with RA and an equal number of age- and sex-matched controls. The patients with RA attended outpatient clinics at King Khalid University Hospital in Riyadh, Saudi Arabia. Blood samples were collected, and DNA was extracted using Qiagen kits. Primers were designed for the 33 selected SNPs using the MassEXTEND primers program, and samples were genotyped on the Sequenom MassARRAY iPLEX platform. The allele frequencies and genotypes were determined for each SNP, and the results obtained for the patients were compared to those for the controls.

RESULTS

The allele and genotype frequencies of six SNPs were significantly associated with RA: rs1188934, rs10919563, rs3087243, rs1980422, rs10499194, and rs629326. The minor alleles of rs1188934, rs10919563, rs10499194, and rs629326 were protective, with odds ratios of 0.542, 0.597, 0.589, and 0.625, and p-values of 0.002, 0.023, 0.013 and 0.036, respectively. In addition, the heterozygote frequencies of two SNPs (rs6859219 and rs11586238) were significantly higher in the controls than in the patients.

CONCLUSIONS

There is considerable heterogeneity in the genetics of RA in different populations, and the SNPs that are associated with RA in some populations are not in others. We studied 33 SNPs and only eight were associated with RA. The remaining SNPs showed no allelic or genotypic associations with RA.

Multi-ancestry genome-wide meta-analysis of 56,241 individuals identifies LRRC4C, LHX5-AS1 and nominates ancestry-specific loci PTPRK , GRB14 , and KIAA0825 as novel risk loci for Alzheimer's disease: the Alzheimer's Disease Genetics Consortium

Limited ancestral diversity has impaired our ability to detect risk variants more prevalent in non-European ancestry groups in genome-wide association studies (GWAS). We constructed and analyzed a multi-ancestry GWAS dataset in the Alzheimer's Disease (AD) Genetics Consortium (ADGC) to test for novel shared and ancestry-specific AD susceptibility loci and evaluate underlying genetic architecture in 37,382 non-Hispanic White (NHW), 6,728 African American, 8,899 Hispanic (HIS), and 3,232 East Asian individuals, performing within-ancestry fixed-effects meta-analysis followed by a cross-ancestry random-effects meta-analysis. We identified 13 loci with cross-ancestry associations including known loci at/near CR1 , BIN1 , TREM2 , CD2AP , PTK2B , CLU , SHARPIN , MS4A6A , PICALM , ABCA7 , APOE and two novel loci not previously reported at 11p12 ( LRRC4C ) and 12q24.13 ( LHX5-AS1 ). Reflecting the power of diverse ancestry in GWAS, we observed the SHARPIN locus using 7.1% the sample size of the original discovering single-ancestry GWAS (n=788,989). We additionally identified three GWS ancestry-specific loci at/near ( PTPRK ( P =2.4×10 -8 ) and GRB14 ( P =1.7×10 -8 ) in HIS), and KIAA0825 ( P =2.9×10 -8 in NHW). Pathway analysis implicated multiple amyloid regulation pathways (strongest with P adjusted =1.6×10 -4 ) and the classical complement pathway ( P adjusted =1.3×10 -3 ). Genes at/near our novel loci have known roles in neuronal development ( LRRC4C, LHX5-AS1 , and PTPRK ) and insulin receptor activity regulation ( GRB14 ). These findings provide compelling support for using traditionally-underrepresented populations for gene discovery, even with smaller sample sizes.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, while context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell type proportions, we demonstrate that cell type iQTLs could be considered as proxies for cell type-specific QTL effects. The interpretation of age iQTLs, however, warrants caution as the moderation effect of age on the genotype and molecular phenotype association may be mediated by changes in cell type composition. Finally, we show that cell type iQTLs contribute to cell type-specific enrichment of diseases that, in combination with additional functional data, may guide future functional studies. Overall, this study highlights iQTLs to gain insights into the context-specificity of regulatory effects.

Identifying shared genetic loci and common risk genes of rheumatoid arthritis associated with three autoimmune diseases based on large-scale cross-trait genome-wide association studies

Introduction

Substantial links between autoimmune diseases have been shown by an increasing number of studies, and one hypothesis for this comorbidity is that there is a common genetic cause.

Methods

In this paper, a large-scale cross-trait Genome-wide Association Studies (GWAS) was conducted to investigate the genetic overlap among rheumatoid arthritis, multiple sclerosis, inflammatory bowel disease and type 1 diabetes.

Results and discussion

Through the local genetic correlation analysis, 2 regions with locally significant genetic associations between rheumatoid arthritis and multiple sclerosis, and 4 regions with locally significant genetic associations between rheumatoid arthritis and type 1 diabetes were discovered. By cross-trait meta-analysis, 58 independent loci associated with rheumatoid arthritis and multiple sclerosis, 86 independent loci associated with rheumatoid arthritis and inflammatory bowel disease, and 107 independent loci associated with rheumatoid arthritis and type 1 diabetes were identified with genome-wide significance. In addition, 82 common risk genes were found through genetic identification. Based on gene set enrichment analysis, it was found that shared genes are enriched in exposed dermal system, calf, musculoskeletal, subcutaneous fat, thyroid and other tissues, and are also significantly enriched in 35 biological pathways. To verify the association between diseases, Mendelian randomized analysis was performed, which shows possible causal associations between rheumatoid arthritis and multiple sclerosis, and between rheumatoid arthritis and type 1 diabetes. The common genetic structure of rheumatoid arthritis, multiple sclerosis, inflammatory bowel disease and type 1 diabetes was explored by these studies, and it is believed that this important discovery will lead to new ideas for clinical treatment.

Identification of genetic variants that impact gene co-expression relationships using large-scale single-cell data

BACKGROUND

Expression quantitative trait loci (eQTL) studies show how genetic variants affect downstream gene expression. Single-cell data allows reconstruction of personalized co-expression networks and therefore the identification of SNPs altering co-expression patterns (co-expression QTLs, co-eQTLs) and the affected upstream regulatory processes using a limited number of individuals.

RESULTS

We conduct a co-eQTL meta-analysis across four scRNA-seq peripheral blood mononuclear cell datasets using a novel filtering strategy followed by a permutation-based multiple testing approach. Before the analysis, we evaluate the co-expression patterns required for co-eQTL identification using different external resources. We identify a robust set of cell-type-specific co-eQTLs for 72 independent SNPs affecting 946 gene pairs. These co-eQTLs are replicated in a large bulk cohort and provide novel insights into how disease-associated variants alter regulatory networks. One co-eQTL SNP, rs1131017, that is associated with several autoimmune diseases, affects the co-expression of RPS26 with other ribosomal genes. Interestingly, specifically in T cells, the SNP additionally affects co-expression of RPS26 and a group of genes associated with T cell activation and autoimmune disease. Among these genes, we identify enrichment for targets of five T-cell-activation-related transcription factors whose binding sites harbor rs1131017. This reveals a previously overlooked process and pinpoints potential regulators that could explain the association of rs1131017 with autoimmune diseases.

CONCLUSION

Our co-eQTL results highlight the importance of studying context-specific gene regulation to understand the biological implications of genetic variation. With the expected growth of sc-eQTL datasets, our strategy and technical guidelines will facilitate future co-eQTL identification, further elucidating unknown disease mechanisms.

Enhancer variants on chromosome 2p14 regulating SPRED2 and ACTR2 act as a signal amplifier to protect against rheumatoid arthritis

Genome-wide association studies (GWASs) have repeatedly reported multiple non-coding single-nucleotide polymorphisms (SNPs) at 2p14 associated with rheumatoid arthritis (RA), but their functional roles in the pathological mechanisms of RA remain to be explored. In this study, we integrated a series of bioinformatics and functional experiments and identified three intronic RA SNPs (rs1876518, rs268131, and rs2576923) within active enhancers that can regulate the expression of SPRED2 directly. At the same time, SPRED2 and ACTR2 influence each other as a positive feedback signal amplifier to strengthen the protective role in RA by inhibiting the migration and invasion of rheumatoid fibroblast-like synoviocytes (FLSs). In particular, the transcription factor CEBPB preferentially binds to the rs1876518-T allele to increase the expression of SPRED2 in FLSs. Our findings decipher the molecular mechanisms behind the GWAS signals at 2p14 for RA and emphasize SPRED2 as a potential candidate gene for RA, providing a potential target and direction for precise treatment of RA.

The lingering effects of Neanderthal introgression on human complex traits

The genetic variants introduced into the ancestors of modern humans from interbreeding with Neanderthals have been suggested to contribute an unexpected extent to complex human traits. However, testing this hypothesis has been challenging due to the idiosyncratic population genetic properties of introgressed variants. We developed rigorous methods to assess the contribution of introgressed Neanderthal variants to heritable trait variation and applied these methods to analyze 235,592 introgressed Neanderthal variants and 96 distinct phenotypes measured in about 300,000 unrelated white British individuals in the UK Biobank. Introgressed Neanderthal variants make a significant contribution to trait variation (explaining 0.12% of trait variation on average). However, the contribution of introgressed variants tends to be significantly depleted relative to modern human variants matched for allele frequency and linkage disequilibrium (about 59% depletion on average), consistent with purifying selection on introgressed variants. Different from previous studies (McArthur et al., 2021), we find no evidence for elevated heritability across the phenotypes examined. We identified 348 independent significant associations of introgressed Neanderthal variants with 64 phenotypes. Previous work (Skov et al., 2020) has suggested that a majority of such associations are likely driven by statistical association with nearby modern human variants that are the true causal variants. Applying a customized fine-mapping led us to identify 112 regions across 47 phenotypes containing 4303 unique genetic variants where introgressed variants are highly likely to have a phenotypic effect. Examination of these variants reveals their substantial impact on genes that are important for the immune system, development, and metabolism.

Multitrait genome-wide analyses identify new susceptibility loci and candidate drugs to primary sclerosing cholangitis

Primary sclerosing cholangitis (PSC) is a rare autoimmune bile duct disease that is strongly associated with immune-mediated disorders. In this study, we implemented multitrait joint analyses to genome-wide association summary statistics of PSC and numerous clinical and epidemiological traits to estimate the genetic contribution of each trait and genetic correlations between traits and to identify new lead PSC risk-associated loci. We identified seven new loci that have not been previously reported and one new independent lead variant in the previously reported locus. Functional annotation and fine-mapping nominated several potential susceptibility genes such as MANBA and IRF5. Network-based in silico drug efficacy screening provided candidate agents for further study of pharmacological effect in PSC.

Role of Liquid Biopsies in Rheumatoid Arthritis

Rheumatoid arthritis (RA) is a chronic autoimmune disease caused by genetic and environmental factors. Early diagnosis is crucial for effective therapy and prognosis of RA, while biomarkers play important roles in early diagnosis. Traditional laboratory tests include rheumatoid factor, anti-cyclic citrullinated peptide antibody, which are inadequate in the ability of early diagnosis. Liquid biopsy technology is a technique using biomarkers found in the blood, urine, and other biological samples from patients, including DNA, RNA, exosome, etc. Evidence indicates that these biomarkers are involved in pathological and physiological conditions of RA. We reviewed the effects of liquid biopsy technology in the early diagnosis of RA and may provide new ideas for effective and precise treatment.

Biomarkers (mRNAs and non-coding RNAs) for the diagnosis and prognosis of rheumatoid arthritis

In recent years, diagnostic and therapeutic approaches for rheumatoid arthritis (RA) have continued to improve. However, in the advanced stages of the disease, patients are unable to achieve long-term clinical remission and often suffer from systemic multi-organ damage and severe complications. Patients with RA usually have no overt clinical manifestations in the early stages, and by the time a definitive diagnosis is made, the disease is already at an advanced stage. RA is diagnosed clinically and with laboratory tests, including the blood markers C-reactive protein (CRP) and erythrocyte sedimentation rate (ESR) and the autoantibodies rheumatoid factor (RF) and anticitrullinated protein antibodies (ACPA). However, the presence of RF and ACPA autoantibodies is associated with aggravated disease, joint damage, and increased mortality, and these autoantibodies have low specificity and sensitivity. The etiology of RA is unknown, with the pathogenesis involving multiple factors and clinical heterogeneity. The early diagnosis, subtype classification, and prognosis of RA remain challenging, and studies to develop minimally invasive or non-invasive biomarkers in the form of biofluid biopsies are becoming more common. Non-coding RNA (ncRNA) molecules are composed of long non-coding RNAs, small nucleolar RNAs, microRNAs, and circular RNAs, which play an essential role in disease onset and progression and can be used in the early diagnosis and prognosis of RA. In this review of the diagnostic and prognostic approaches to RA disease, we provide an overview of the current knowledge on the subject, focusing on recent advances in mRNA-ncRNA as diagnostic and prognostic biomarkers from the biofluid to the tissue level.

Multi-ancestry genome-wide association analyses identify novel genetic mechanisms in rheumatoid arthritis

Rheumatoid arthritis (RA) is a highly heritable complex disease with unknown etiology. Multi-ancestry genetic research of RA promises to improve power to detect genetic signals, fine-mapping resolution and performances of polygenic risk scores (PRS). Here, we present a large-scale genome-wide association study (GWAS) of RA, which includes 276,020 samples from five ancestral groups. We conducted a multi-ancestry meta-analysis and identified 124 loci (P < 5 × 10-8), of which 34 are novel. Candidate genes at the novel loci suggest essential roles of the immune system (for example, TNIP2 and TNFRSF11A) and joint tissues (for example, WISP1) in RA etiology. Multi-ancestry fine-mapping identified putatively causal variants with biological insights (for example, LEF1). Moreover, PRS based on multi-ancestry GWAS outperformed PRS based on single-ancestry GWAS and had comparable performance between populations of European and East Asian ancestries. Our study provides several insights into the etiology of RA and improves the genetic predictability of RA.

Heavy Chain Constant Region Usage in Antibodies to Peptidylarginine Deiminase 4 as a Marker of Disease Subsets in Rheumatoid Arthritis

OBJECTIVE

The study of autoantibody isotypes in autoimmune diseases is useful for identifying clinically relevant endotypes. This study was undertaken to study the prevalence and clinical significance of different isotypes and IgG subclasses of anti-peptidylarginine deiminase 4 (anti-PAD4) autoantibodies in individuals with rheumatoid arthritis (RA).

METHODS

In 196 RA subjects and 64 healthy controls, anti-PAD4 antibody types were determined using enzyme-linked immunosorbent assay. We investigated associations between anti-PAD4 antibodies and clinical outcomes, and relevant features were confirmed in an independent RA cohort.

RESULTS

Anti-PAD4 IgG1, anti-PAD4 IgG2, anti-PAD4 IgG3, anti-PAD4 IgG4, anti-PAD4 IgA, and anti-PAD4 IgE antibodies were more frequent in RA patients than healthy controls (P < 0.001). Anti-PAD4 IgG1, anti-PAD4 IgG3, and anti-PAD4 IgE were associated with distinct clinical features. Anti-PAD4 IgG1 was predictive of progressive radiographic joint damage (odds ratio [OR] 4.88, P = 0.005), especially in RA patients without baseline joint damage (40% versus 0%, P = 0.003) or in those negative for anti-cyclic citrullinated peptide and/or rheumatoid factor (OR 32; P = 0.009). IgG1 was also associated with higher levels of C-reactive protein (P = 0.006) and interleukin-6 (P = 0.021). RA patients with anti-PAD4 IgG3 had higher baseline joint damage scores (median Sharp/van der Heijde score 13 versus 7, P = 0.046), while those with anti-PAD4 IgE had higher Disease Activity Score in 28 joints (median 4.0 versus 3.5, P = 0.025), more frequent rheumatoid nodules (31% versus 16%, P = 0.025), and more frequent interstitial lung disease (ground-glass opacification) (24% versus 9%, P = 0.014). Anti-PAD4 IgG1 antibody associations with joint damage were corroborated in an independent RA cohort.

CONCLUSION

Anti-PAD4 IgG1, anti-PAD4 IgG3, and anti-PAD4 IgE antibodies identify discrete disease subsets in RA, suggesting that heavy chain usage drives distinct effector mechanisms of anti-PAD4 antibodies in RA.

Evolution of immune genes is associated with the Black Death

Infectious diseases are among the strongest selective pressures driving human evolution1,2. This includes the single greatest mortality event in recorded history, the first outbreak of the second pandemic of plague, commonly called the Black Death, which was caused by the bacterium Yersinia pestis3. This pandemic devastated Afro-Eurasia, killing up to 30-50% of the population4. To identify loci that may have been under selection during the Black Death, we characterized genetic variation around immune-related genes from 206 ancient DNA extracts, stemming from two different European populations before, during and after the Black Death. Immune loci are strongly enriched for highly differentiated sites relative to a set of non-immune loci, suggesting positive selection. We identify 245 variants that are highly differentiated within the London dataset, four of which were replicated in an independent cohort from Denmark, and represent the strongest candidates for positive selection. The selected allele for one of these variants, rs2549794, is associated with the production of a full-length (versus truncated) ERAP2 transcript, variation in cytokine response to Y. pestis and increased ability to control intracellular Y. pestis in macrophages. Finally, we show that protective variants overlap with alleles that are today associated with increased susceptibility to autoimmune diseases, providing empirical evidence for the role played by past pandemics in shaping present-day susceptibility to disease.

rs2013278 in the multiple immunological-trait susceptibility locus CD28 regulates the production of non-functional splicing isoforms

Background

Ligation of CD28 with ligands such as CD80 or CD86 provides a critical second signal alongside antigen presentation by class II major histocompatibility complex expressed on antigen-presenting cells through the T cell antigen receptor for naïve T cell activation. A number of studies suggested that CD28 plays an important role in the pathogenesis of various human diseases. Recent genome-wide association studies (GWASs) identified CD28 as a susceptibility locus for lymphocyte and eosinophil counts, multiple sclerosis, ulcerative colitis, celiac disease, rheumatoid arthritis, asthma, and primary biliary cholangitis. However, the primary functional variant and molecular mechanisms of disease susceptibility in this locus remain to be elucidated. This study aimed to identify the primary functional variant from thousands of genetic variants in the CD28 locus and elucidate its functional effect on the CD28 molecule.

Results

Among the genetic variants exhibiting stronger linkage disequilibrium (LD) with all GWAS-lead variants in the CD28 locus, rs2013278, located in the Rbfox binding motif related to splicing regulation, was identified as a primary functional variant related to multiple immunological traits. Relative endogenous expression levels of CD28 splicing isoforms (CD28i and CD28Δex2) compared with full-length CD28 in allele knock-in cell lines generated using CRISPR/Cas9 were directly regulated by rs2013278 (P < 0.05). Although full-length CD28 protein expressed on Jurkat T cells showed higher binding affinity for CD80/CD86, both CD28i and CD28Δex2 encoded loss-of-function isoforms.

Conclusion

The present study demonstrated for the first time that CD28 has a shared disease-related primary functional variant (i.e., rs2013278) that regulates the CD28 alternative splicing that generates loss-of-function isoforms. They reduce disease risk by inducing anergy of effector T cells that over-react to autoantigens and allergens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-022-00419-7.

Towards a global view of multiple sclerosis genetics

Multiple sclerosis (MS) is a neuroimmunological disorder of the CNS with a strong heritable component. The genetic architecture of MS susceptibility is well understood in populations of European ancestry. However, the extent to which this architecture explains MS susceptibility in populations of non-European ancestry remains unclear. In this Perspective article, we outline the scientific arguments for studying MS genetics in ancestrally diverse populations. We argue that this approach is likely to yield insights that could benefit individuals with MS from all ancestral groups. We explore the logistical and theoretical challenges that have held back this field to date and conclude that, despite these challenges, inclusion of participants of non-European ancestry in MS genetics studies will ultimately be of value to all patients with MS worldwide.

A genome-wide association study for rheumatoid arthritis replicates previous HLA and non-HLA associations in a cohort from South Africa

The complex pathogenesis of rheumatoid arthritis (RA) is not fully understood, with few studies exploring the genomic contribution to RA in patients from Africa. We report a genome-wide association study (GWAS) of South-Eastern Bantu-Speaking South Africans (SEBSSAs) with seropositive RA (n = 531) and population controls (n = 2653). Association testing was performed using PLINK (logistic regression assuming an additive model) with sex, age, smoking and the first three principal components as covariates. The strong association with the Human Leukocyte Antigen (HLA) region, indexed by rs602457 (near HLA-DRB1), was replicated. An additional independent signal in the HLA region represented by the lead SNP rs2523593 (near the HLA-B gene; Conditional P-value = 6.4 × 10-10) was detected. Although none of the non-HLA signals reached genome-wide significance (P < 5 × 10-8), 17 genomic regions showed suggestive association (P < 5 × 10-6). The GWAS replicated two known non-HLA associations with MMEL1 (rs2843401) and ANKRD55 (rs7731626) at a threshold of P < 5 × 10-3 providing, for the first time, evidence for replication of non-HLA signals for RA in sub-Saharan African populations. Meta-analysis with summary statistics from an African-American cohort (CLEAR study) replicated three additional non-HLA signals (rs11571302, rs2558210 and rs2422345 around KRT18P39-NPM1P33, CTLA4-ICOS and AL645568.1, respectively). Analysis based on genomic regions (200 kb windows) further replicated previously reported non-HLA signals around PADI4, CD28 and LIMK1. Although allele frequencies were overall strongly correlated between the SEBSSA and the CLEAR cohort, we observed some differences in effect size estimates for associated loci. The study highlights the need for conducting larger association studies across diverse African populations to inform precision medicine-based approaches for RA in Africa.

A Summary on the Genetics of Systemic Lupus Erythematosus, Rheumatoid Arthritis, Systemic Sclerosis, and Sjögren's Syndrome

Systemic lupus erythematosus, systemic sclerosis, rheumatoid arthritis, and Sjögren's syndrome are four major autoimmune rheumatic diseases characterized by the presence of autoantibodies, caused by a dysregulation of the immune system that leads to a wide variety of clinical manifestations. These conditions present complex etiologies strongly influenced by multiple environmental and genetic factors. The human leukocyte antigen (HLA) region was the first locus identified to be associated and still represents the strongest susceptibility factor for each of these conditions, particularly the HLA class II genes, including DQA1, DQB1, and DRB1, but class I genes have also been associated. Over the last two decades, the genetic component of these disorders has been extensively investigated and hundreds of non-HLA risk genetic variants have been uncovered. Furthermore, it is widely accepted that autoimmune rheumatic diseases share molecular disease pathways, such as the interferon (IFN) type I pathways, which are reflected in a common genetic background. Some examples of well-known pleiotropic loci for autoimmune rheumatic diseases are the HLA region, DNASEL13, TNIP1, and IRF5, among others. The identification of the causal molecular mechanisms behind the genetic associations is still a challenge. However, recent advances have been achieved through mouse models and functional studies of the loci. Here, we provide an updated overview of the genetic architecture underlying these four autoimmune rheumatic diseases, with a special focus on the HLA region.

Single-cell eQTL models reveal dynamic T cell state dependence of disease loci

Non-coding genetic variants may cause disease by modulating gene expression. However, identifying these expression quantitative trait loci (eQTLs) is complicated by differences in gene regulation across fluid functional cell states within cell types. These states-for example, neurotransmitter-driven programs in astrocytes or perivascular fibroblast differentiation-are obscured in eQTL studies that aggregate cells1,2. Here we modelled eQTLs at single-cell resolution in one complex cell type: memory T cells. Using more than 500,000 unstimulated memory T cells from 259 Peruvian individuals, we show that around one-third of 6,511 cis-eQTLs had effects that were mediated by continuous multimodally defined cell states, such as cytotoxicity and regulatory capacity. In some loci, independent eQTL variants had opposing cell-state relationships. Autoimmune variants were enriched in cell-state-dependent eQTLs, including risk variants for rheumatoid arthritis near ORMDL3 and CTLA4; this indicates that cell-state context is crucial to understanding potential eQTL pathogenicity. Moreover, continuous cell states explained more variation in eQTLs than did conventional discrete categories, such as CD4+ versus CD8+, suggesting that modelling eQTLs and cell states at single-cell resolution can expand insight into gene regulation in functionally heterogeneous cell types.

Association Study of Anticitrullinated Peptide Antibody Status with Clinical Manifestations and SNPs in Patients Affected with Rheumatoid Arthritis: A Pilot Study

Introduction. Rheumatoid arthritis (RA) is an autoimmune disease of unknown etiology that leads to disability due to articular and extra-articular damage. RA prevalence is variable. The disease is most common among females with a 3 : 1 ratio. The interaction of environmental and host factors contributes to RA development. Currently, the genome-wide association studies (GWAS) give the opportunity to uncover the RA genetic background. Anticitrullinated peptide antibody (ACPA) is a highly specific RA antibody, associated with poor prognosis and severe course of RA, and regulated by numerous genes. Our study is aimed at investigating whether there are any clinical and genetic aspects correlate with ACPA presence in Kazakhstani patients with RA. Indeed, the available studies on this subject are focused on Caucasian and East Asian populations (mainly Japanese and Chinese), and there are scarce data from Central Asia. Methods. Our study included 70 RA patients. Patients’ blood samples were collected and genotyped for 14 SNPs by real-time polymerase chain reaction (RT-PCR). General examination, anamnestic, and clinical and laboratory data collection were carried out. Statistical analysis was performed using R statistics. Results and Conclusion. Our study revealed a significant association of ACPA positivity with Fc receptor-like 3 (FCRL3) and ACPA negativity with signal transducer and activator of transcription 4 (STAT4) genes, but not with T cell activation Rho GTPase activating protein (TAGAP). In addition, ACPA positivity was associated with radiographic progression, rheumatoid factor (RF), erythrocyte sedimentation rate (ESR), age of RA onset, the patient global assessment, body mass index (BMI), and Gamma globulin. Conclusion. Remained 11 earlier identified significantly associated in Caucasian and Asian population SNPs were not replicated in our cohort. Further studies on larger cohorts are needed to confirm our findings with higher confidence levels and stronger statistical power.

Worse outcomes linked to ethnicity for early inflammatory arthritis in England and Wales: a national cohort study

OBJECTIVE

To assess variability in care quality and treatment outcomes across ethnicities in early inflammatory arthritis (EIA).

METHODS

We conducted an observational cohort study in England and Wales from May 2018 to March 2020, including patients with a suspected/confirmed EIA diagnosis. Care quality was assessed against six metrics defined by national guidelines. Clinical outcomes were measured using DAS28. Outcomes between ethnic groups ('White', 'Black', 'Asian', 'Mixed', 'Other') were compared, and adjusted for confounders.

RESULTS

A total of 35 807 eligible patients were analysed. Of those, 30 643 (85.6%) were White and 5164 (14.6%) were from ethnic minorities: 1035 (2.8%) Black; 2617 (7.3%) Asian; 238 (0.6%) Mixed; 1274 (3.5%) Other. In total, 12 955 patients had confirmed EIA, of whom 11 315 were White and 1640 were from ethnic minorities: 314 (2.4%) Black; 927 (7.1%) Asian; 70 (0.5%) Mixed; 329 (2.5%) Other. A total of 14 803 patients were assessed by rheumatology within three weeks, and 5642 started treatment within six weeks of referral. There were no significant differences by ethnicity. Ethnic minority patients had lower odds of disease remission at three months [adjusted odds ratio 0.79 (95% CI: 0.65, 0.96)] relative to White patients. Ethnic minorities were significantly less likely to receive initial treatment withMTX[0.68 (0.52, 0.90)] or with glucocorticoids [0.63 (0.49, 0.80)].

CONCLUSION

We demonstrate that some ethnic minorities are less likely to achieve disease remission in three months following EIA diagnosis. This is not explained by delays in referral or time to treatment. Our data highlight the need for investigation into the possible drivers of these inequitable outcomes and reappraisal of EIA management pathways.

Early DNA methylation changes in children developing beta cell autoimmunity at a young age

AIMS/HYPOTHESIS

Type 1 diabetes is a chronic autoimmune disease of complex aetiology, including a potential role for epigenetic regulation. Previous epigenomic studies focused mainly on clinically diagnosed individuals. The aim of the study was to assess early DNA methylation changes associated with type 1 diabetes already before the diagnosis or even before the appearance of autoantibodies.

METHODS

Reduced representation bisulphite sequencing (RRBS) was applied to study DNA methylation in purified CD4+ T cell, CD8+ T cell and CD4-CD8- cell fractions of 226 peripheral blood mononuclear cell samples longitudinally collected from seven type 1 diabetes-specific autoantibody-positive individuals and control individuals matched for age, sex, HLA risk and place of birth. We also explored correlations between DNA methylation and gene expression using RNA sequencing data from the same samples. Technical validation of RRBS results was performed using pyrosequencing.

RESULTS

We identified 79, 56 and 45 differentially methylated regions in CD4+ T cells, CD8+ T cells and CD4-CD8- cell fractions, respectively, between type 1 diabetes-specific autoantibody-positive individuals and control participants. The analysis of pre-seroconversion samples identified DNA methylation signatures at the very early stage of disease, including differential methylation at the promoter of IRF5 in CD4+ T cells. Further, we validated RRBS results using pyrosequencing at the following CpG sites: chr19:18118304 in the promoter of ARRDC2; chr21:47307815 in the intron of PCBP3; and chr14:81128398 in the intergenic region near TRAF3 in CD4+ T cells.

CONCLUSIONS/INTERPRETATION

These preliminary results provide novel insights into cell type-specific differential epigenetic regulation of genes, which may contribute to type 1 diabetes pathogenesis at the very early stage of disease development. Should these findings be validated, they may serve as a potential signature useful for disease prediction and management.

Genetic Associations and Differential mRNA Expression Levels of Host Genes Suggest a Viral Trigger for Endemic Pemphigus Foliaceus

The long search for the environmental trigger of the endemic pemphigus foliaceus (EPF, fogo selvagem) has not yet resulted in any tangible findings. Here, we searched for genetic associations and the differential expression of host genes involved in early viral infections and innate antiviral defense. Genetic variants could alter the structure, expression sites, or levels of the gene products, impacting their functions. By analyzing 3063 variants of 166 candidate genes in 227 EPF patients and 194 controls, we found 12 variants within 11 genes associated with differential susceptibility (p < 0.005) to EPF. The products of genes TRIM5, TPCN2, EIF4E, EIF4E3, NUP37, NUP50, NUP88, TPR, USP15, IRF8, and JAK1 are involved in different mechanisms of viral control, for example, the regulation of viral entry into the host cell or recognition of viral nucleic acids and proteins. Only two of nine variants were also associated in an independent German cohort of sporadic PF (75 patients, 150 controls), aligning with our hypothesis that antiviral host genes play a major role in EPF due to a specific virus−human interaction in the endemic region. Moreover, CCL5, P4HB, and APOBEC3G mRNA levels were increased (p < 0.001) in CD4+ T lymphocytes of EPF patients. Because there is limited or no evidence that these genes are involved in autoimmunity, their crucial role in antiviral responses and the associations that we observed support the hypothesis of a viral trigger for EPF, presumably a still unnoticed flavivirus. This work opens new frontiers in searching for the trigger of EPF, with the potential to advance translational research that aims for disease prevention and treatment.

Genome-Wide Association Study of Fluorescent Oxidation Products Accounting for Tobacco Smoking Status in Adults from the French EGEA Study

Oxidative stress (OS) is the main pathophysiological mechanism involved in several chronic diseases, including asthma. Fluorescent oxidation products (FlOPs), a global biomarker of damage due to OS, is of growing interest in epidemiological studies. We conducted a genome-wide association study (GWAS) of the FlOPs level in 1216 adults from the case-control and family-based EGEA study (mean age 43 years old, 51% women, and 23% current smokers) to identify genetic variants associated with FlOPs. The GWAS was first conducted in the whole sample and then stratified according to smoking status, the main exogenous source of reactive oxygen species. Among the top genetic variants identified by the three GWAS, those located in BMP6 (p = 3 × 10⁻⁶), near BMPER (p = 9 × 10⁻⁶), in GABRG3 (p = 4 × 10⁻⁷), and near ATG5 (p = 2 × 10⁻⁹) are the most relevant because of both their link to biological pathways related to OS and their association with several chronic diseases for which the role of OS in their pathophysiology has been pointed out. BMP6 and BMPER are of particular interest due to their involvement in the same biological pathways related to OS and their functional interaction. To conclude, this study, which is the first GWAS of FlOPs, provides new insights into the pathophysiology of chronic OS-related diseases.

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.

Sequence complementarity between human noncoding RNAs and SARS-CoV-2 genes: What are the implications for human health?

Objectives

To investigate in silico the presence of nucleotide sequence complementarity between the RNA genome of Severe Acute Respiratory Syndrome CoronaVirus-2 (SARS-CoV-2) and human non-coding (nc)RNA genes.

Methods

The FASTA sequence (NC_045512.2) of each of the 11 SARS-CoV-2 isolate Wuhan-Hu-1 genes was retrieved from NCBI.nlm.nih.gov/gene and the Ensembl.org library interrogated for any base-pair match with human ncRNA genes. SARS-CoV-2 gene-matched human ncRNAs were screened for functional activity using bioinformatic analysis. Finally, associations between identified ncRNAs and human diseases were searched in GWAS databases.

Results

A total of 252 matches were found between the nucleotide sequence of SARS-CoV-2 genes and human ncRNAs. With the exception of two small nuclear RNAs, all of them were long non-coding (lnc)RNAs expressed mainly in testis and central nervous system under physiological conditions. The percentage of alignment ranged from 91.30% to 100% with a mean nucleotide alignment length of 17.5 ± 2.4. Thirty-three (13.09%) of them contained predicted R-loop forming sequences, but none of these intersected the complementary sequences of SARS-CoV-2. However, in 31 cases matches fell on ncRNA regulatory sites, whose adjacent coding genes are mostly involved in cancer, immunological and neurological pathways. Similarly, several polymorphic variants of detected non-coding genes have been associated with neuropsychiatric and proliferative disorders.

Conclusion

This pivotal in silico study shows that SARS-CoV-2 genes have Watson-Crick nucleotide complementarity to human ncRNA sequences, potentially disrupting ncRNA epigenetic control of target genes. It remains to be elucidated whether this could result in the development of human disease in the long term.

Genomic Multiple Sclerosis Risk Variants Modulate the Expression of the ANKRD55-IL6ST Gene Region in Immature Dendritic Cells

Intronic single-nucleotide polymorphisms (SNPs) in the ANKRD55 gene are associated with the risk for multiple sclerosis (MS) and rheumatoid arthritis by genome-wide association studies (GWAS). The risk alleles have been linked to higher expression levels of ANKRD55 and the neighboring IL6ST (gp130) gene in CD4+ T lymphocytes of healthy controls. The biological function of ANKRD55, its role in the immune system, and cellular sources of expression other than lymphocytes remain uncharacterized. Here, we show that monocytes gain capacity to express ANKRD55 during differentiation in immature monocyte-derived dendritic cells (moDCs) in the presence of interleukin (IL)-4/granulocyte-macrophage colony-stimulating factor (GM-CSF). ANKRD55 expression levels are further enhanced by retinoic acid agonist AM580 but downregulated following maturation with interferon (IFN)-γ and lipopolysaccharide (LPS). ANKRD55 was detected in the nucleus of moDC in nuclear speckles. We also analyzed the adjacent IL6ST, IL31RA, and SLC38A9 genes. Of note, in healthy controls, MS risk SNP genotype influenced ANKRD55 and IL6ST expression in immature moDC in opposite directions to that in CD4+ T cells. This effect was stronger for a partially correlated SNP, rs13186299, that is located, similar to the main MS risk SNPs, in an ANKRD55 intron. Upon analysis in MS patients, the main GWAS MS risk SNP rs7731626 was associated with ANKRD55 expression levels in CD4+ T cells. MoDC-specific ANKRD55 and IL6ST mRNA levels showed significant differences according to the clinical form of the disease, but, in contrast to healthy controls, were not influenced by genotype. We also measured serum sgp130 levels, which were found to be higher in homozygotes of the protective allele of rs7731626. Our study characterizes ANKRD55 expression in moDC and indicates monocyte-to-dendritic cell (Mo-DC) differentiation as a process potentially influenced by MS risk SNPs.

Hao1 Is Not a Pathogenic Factor for Ectopic Ossifications but Functions to Regulate the TCA Cycle In Vivo

Ossification of the posterior longitudinal ligament (OPLL), a disease characterized by the ectopic ossification of a spinal ligament, promotes neurological disorders associated with spinal canal stenosis. While blocking ectopic ossification is mandatory to prevent OPLL development and progression, the mechanisms underlying the condition remain unknown. Here we show that expression of hydroxyacid oxidase 1 (Hao1), a gene identified in a previous genome-wide association study (GWAS) as an OPLL-associated candidate gene, specifically and significantly decreased in fibroblasts during osteoblast differentiation. We then newly established Hao1-deficient mice by generating Hao1-flox mice and crossing them with CAG-Cre mice to yield global Hao1-knockout (CAG-Cre/Hao1flox/flox; Hao1 KO) animals. Hao1 KO mice were born normally and exhibited no obvious phenotypes, including growth retardation. Moreover, Hao1 KO mice did not exhibit ectopic ossification or calcification. However, urinary levels of some metabolites of the tricarboxylic acid (TCA) cycle were significantly lower in Hao1 KO compared to control mice based on comprehensive metabolomic analysis. Our data indicate that Hao1 loss does not promote ectopic ossification, but rather that Hao1 functions to regulate the TCA cycle in vivo.

Molecular and Cellular Heterogeneity in Rheumatoid Arthritis: Mechanisms and Clinical Implications

Rheumatoid arthritis is an autoimmune disease that exhibits significant clinical heterogeneity. There are various treatments for rheumatoid arthritis, including disease-modifying anti-rheumatic drugs (DMARDs), glucocorticoids, non-steroidal anti-inflammatory drugs (NSAIDs), and inflammatory cytokine inhibitors (ICI), typically associated with differentiated clinical effects and characteristics. Personalized responsiveness is observed to the standard treatment due to the pathophysiological heterogeneity in rheumatoid arthritis, resulting in an overall poor prognosis. Understanding the role of individual variation in cellular and molecular mechanisms related to rheumatoid arthritis will considerably improve clinical care and patient outcomes. In this review, we discuss the source of pathophysiological heterogeneity derived from genetic, molecular, and cellular heterogeneity and their possible impact on precision medicine and personalized treatment of rheumatoid arthritis. We provide emphasized description of the heterogeneity derived from mast cells, monocyte cell, macrophage fibroblast-like synoviocytes and, interactions within immune cells and with inflammatory cytokines, as well as the potential as a new therapeutic target to develop a novel treatment approach. Finally, we summarize the latest clinical trials of treatment options for rheumatoid arthritis and provide a suggestive framework for implementing preclinical and clinical experimental results into clinical practice.

MicroRNA Variants and HLA-miRNA Interactions are Novel Rheumatoid Arthritis Susceptibility Factors

Genome-wide association studies have identified >100 genetic risk factors for rheumatoid arthritis. However, the reported genetic variants could only explain less than 40% heritability of rheumatoid arthritis. The majority of the heritability is still missing and needs to be identified with more studies with different approaches and populations. In order to identify novel function SNPs to explain missing heritability and reveal novel mechanism pathogenesis of rheumatoid arthritis, 4 HLA SNPs (HLA-DRB1, HLA-DRB9, HLA-DQB1, and TNFAIP3) and 225 common SNPs located in miRNA, which might influence the miRNA target binding or pre-miRNA stability, were genotyped in 1,607 rheumatoid arthritis and 1,580 matched normal individuals. We identified 2 novel SNPs as significantly associated with rheumatoid arthritis including rs1414273 (miR-548ac, OR = 0.84, p = 8.26 × 10^-4) and rs2620381 (miR-627, OR = 0.77, p = 2.55 × 10^-3). We also identified that rs5997893 (miR-3928) showed significant epistasis effect with rs4947332 (HLA-DRB1, OR = 4.23, p = 0.04) and rs2967897 (miR-5695) with rs7752903 (TNFAIP3, OR = 4.43, p = 0.03). In addition, we found that individuals who carried 8 risk alleles showed 15.38 (95%CI: 4.69-50.49, p < 1.0 × 10^-6) times more risk of being affected by RA. Finally, we demonstrated that the targets of the significant miRNAs showed enrichment in immune related genes (p = 2.0 × 10^-5) and FDA approved drug target genes (p = 0.014). Overall, 6 novel miRNA SNPs including rs1414273 (miR-548ac, p = 8.26 × 10^-4), rs2620381 (miR-627, p = 2.55 × 10^-3), rs4285314 (miR-3135b, p = 1.10 × 10^-13), rs28477407 (miR-4308, p = 3.44 × 10^-5), rs5997893 (miR-3928, p = 5.9 × 10^-3) and rs45596840 (miR-4482, p = 6.6 × 10^-3) were confirmed to be significantly associated with RA in a Chinese population. Our study suggests that miRNAs might be interesting targets to accelerate understanding of the pathogenesis and drug development for rheumatoid arthritis.

The immunological significance of tumor necrosis factor receptor-associated factors (TRAFs)

The tumor necrosis factor receptor (TNFR)-associated factor (TRAF) family of molecules are intracellular signaling adaptors and control diverse signaling pathways mediated not only by the TNFR superfamily and the Toll-like receptor/interleukin-1 receptor superfamily but also by unconventional cytokine receptors such as IL-6 and IL-17 receptors. There are seven family members, TRAF1 to TRAF7, in mammals. Exaggerated immune responses induced through TRAF signaling downstream of these receptors often lead to inflammatory and autoimmune diseases including rheumatoid arthritis, inflammatory bowel disease, psoriasis and autoinflammatory syndromes, and thus those signals are major targets for therapeutic intervention. For this reason, it has been very important to understand signaling mechanisms regulated by TRAFs that greatly impact on life/death decisions and the activation, differentiation and survival of cells of the innate and adaptive immune systems. Accumulating evidence suggests that dysregulated cellular expression and/or signaling of TRAFs causes overproduction of proinflammatory cytokines, which facilitates aberrant activation of immune cells. In this review, I will explain the structural and functional aspects that are responsible for the cellular activity and disease outcomes of TRAFs, and summarize the findings of recent studies on TRAFs in terms of how individual TRAF family molecules regulates biological and disease processes in the body in both positive and negative ways. This review also discusses how TRAF mutations contribute to human disease.

Cell-Free DNA in Rheumatoid Arthritis

Endogenous DNA derived from the nuclei or mitochondria is released into the bloodstream following cell damage or death. Extracellular DNA, called cell-free DNA (cfDNA), is associated with various pathological conditions. Recently, multiple aspects of cfDNA have been assessed, including cfDNA levels, integrity, methylation, and mutations. Rheumatoid arthritis (RA) is the most common form of autoimmune arthritis, and treatment of RA has highly varied outcomes. cfDNA in patients with RA is elevated in peripheral blood and synovial fluid and is associated with disease activity. Profiling of cfDNA in patients with RA may then be utilized in various aspects of clinical practice, such as the prediction of prognosis and treatment responses; monitoring disease state; and as a diagnostic marker. In this review, we discuss cfDNA in patients with RA, particularly the sources of cfDNA and the correlation of cfDNA with RA pathogenesis. We also highlight the potential of analyzing cfDNA profiles to guide individualized treatment approaches for RA.

A novel Mendelian randomization method with binary risk factor and outcome

BACKGROUND

Mendelian randomization (MR) applies instrumental variable (IV) methods to observational data using a genetic variant as an IV. Several Monte-Carlo studies investigate the performance of MR methods with binary outcomes, but few consider them in conjunction with binary risk factors.

OBJECTIVE

To develop a novel MR estimator for scenarios with a binary risk factor and outcome; and compare to existing MR estimators via simulations and real data analysis.

METHODS

A bivariate Bernoulli distribution is adapted to the IV setting. Empirical bias and asymptotic coverage probabilities are estimated via simulations. The proposed method is compared to the Wald method, two-stage predictor substitution (2SPS), two-stage residual inclusion (2SRI), and the generalized method of moments (GMM). An analysis is performed using existing data from the CLEAR study to estimate the potential causal effect of smoking on rheumatoid arthritis risk in African Americans.

RESULTS

Bias was low for the proposed method and comparable to 2SPS. The Wald method was often biased towards the null. Coverage was adequate for the proposed method, 2SPS, and 2SRI. Coverage for the Wald and GMM methods was poor in several scenarios. The causal effect of ever smoking on rheumatoid arthritis risk was not statistically significant using a variety of genetic instruments.

CONCLUSIONS

Simulations suggest the proposed MR method is sound with binary risk factors and outcomes, and comparable to 2SPS and 2SRI in terms of bias. The proposed method also provides more natural framework for hypothesis testing compared to 2SPS or 2SRI, which require ad-hoc variance adjustments.

Integrative genomic analysis in African American children with asthma finds three novel loci associated with lung function

Bronchodilator (BD) drugs are commonly prescribed for treatment and management of obstructive lung function present with diseases such as asthma. Administration of BD medication can partially or fully restore lung function as measured by pulmonary function tests. The genetics of baseline lung function measures taken before BD medication have been extensively studied, and the genetics of the BD response itself have received some attention. However, few studies have focused on the genetics of post-BD lung function. To address this gap, we analyzed lung function phenotypes in 1103 subjects from the Study of African Americans, Asthma, Genes, and Environment, a pediatric asthma case-control cohort, using an integrative genomic analysis approach that combined genotype, locus-specific genetic ancestry, and functional annotation information. We integrated genome-wide association study (GWAS) results with an admixture mapping scan of three pulmonary function tests (forced expiratory volume in 1 s [FEV1 ], forced vital capacity [FVC], and FEV1 /FVC) taken before and after albuterol BD administration on the same subjects, yielding six traits. We identified 18 GWAS loci, and five additional loci from admixture mapping, spanning several known and novel lung function candidate genes. Most loci identified via admixture mapping exhibited wide variation in minor allele frequency across genotyped global populations. Functional fine-mapping revealed an enrichment of epigenetic annotations from peripheral blood mononuclear cells, fetal lung tissue, and lung fibroblasts. Our results point to three novel potential genetic drivers of pre- and post-BD lung function: ADAMTS1, RAD54B, and EGLN3.

The genetics of rheumatoid arthritis

RA is a chronic systemic inflammatory disease that primarily affects the small joints of the hands and feet, and results in a mean reduction in life expectancy of 3-10 years. RA is a multigene disorder with a substantial genetic component and a heritability estimate of 60%. Large-scale Genome-Wide Association Studies (GWAS) and meta-analyses have revealed common disease-associated variants in the population that may contribute cumulatively to RA pathogenesis. This review identifies the most significant genetic variants associated with RA susceptibility to date, with particular focus on the contribution of the HLA class II genes across different ethnic groups. Also discussed are the potential applications of pharmacogenomics to RA management by identifying polymorphisms associated with variation in treatment response or toxicity. The use of genetic variants to guide treatment strategy has the potential to not only reduce National Health Service costs, but also drastically improve patient experience and quality of life.

High Bundled Payment Costs Following Total Hip and Knee Arthroplasty in African Americans

BACKGROUND

Concerns exist that minorities who utilize more resources in an episode-of-care following total hip (THA) and knee arthroplasty (TKA) may face difficulties with access to quality arthroplasty care in bundled payment programs. The purpose of this study is to determine if African American patients undergoing TKA or THA have higher episode-of-care costs compared to Caucasian patients.

METHODS

We queried Medicare claims data for a consecutive series of 7310 primary TKA and THA patients at our institution from 2015 to 2018. We compared patient demographics, comorbidities, readmissions, and 90-day episode-of-care costs between African American and Caucasian patients. A multivariate regression analysis was performed to identify the independent effect of race on episode-of-care costs.

RESULTS

Compared to Caucasians, African Americans were younger, but had higher rates of pulmonary disease and diabetes. African American patients had increased rates of discharge to a rehabilitation facility (20% vs 13%, P < .001), with higher subacute rehabilitation ($1909 vs $1284, P < .001), home health ($819 vs $698, P = .022), post-acute care ($5656 vs $4961, P = .008), and overall 90-day episode-of-care costs ($19,457 vs $18,694, P = .001). When controlling for confounding comorbidities, African American race was associated with higher episode-of-care costs of $440 (P < .001).

CONCLUSION

African American patients have increased episode-of-care costs following THA and TKA when compared to Caucasian patients, mainly due to increased rates of home health and rehabilitation utilization. Further study is needed to identify social variables that can help reduce post-acute care resources and prevent reduction in access to arthroplasty care in bundled payment models.

Extreme enrichment of VNTR-associated polymorphicity in human subtelomeres: genes with most VNTRs are predominantly expressed in the brain

The human genome harbors numerous structural variants (SVs) which, due to their repetitive nature, are currently underexplored in short-read whole-genome sequencing approaches. Using single-molecule, real-time (SMRT) long-read sequencing technology in combination with FALCON-Unzip, we generated a de novo assembly of the diploid genome of a 115-year-old Dutch cognitively healthy woman. We combined this assembly with two previously published haploid assemblies (CHM1 and CHM13) and the GRCh38 reference genome to create a compendium of SVs that occur across five independent human haplotypes using the graph-based multi-genome aligner REVEAL. Across these five haplotypes, we detected 31,680 euchromatic SVs (>50 bp). Of these, ~62% were comprised of repetitive sequences with 'variable number tandem repeats' (VNTRs), ~10% were mobile elements (Alu, L1, and SVA), while the remaining variants were inversions and indels. We observed that VNTRs with GC-content >60% and repeat patterns longer than 15 bp were 21-fold enriched in the subtelomeric regions (within 5 Mb of the ends of chromosome arms). VNTR lengths can expand to exceed a critical length which is associated with impaired gene transcription. The genes that contained most VNTRs, of which PTPRN2 and DLGAP2 are the most prominent examples, were found to be predominantly expressed in the brain and associated with a wide variety of neurological disorders. Repeat-induced variation represents a sizeable fraction of the genetic variation in human genomes and should be included in investigations of genetic factors associated with phenotypic traits, specifically those associated with neurological disorders. We make available the long and short-read sequence data of the supercentenarian genome, and a compendium of SVs as identified across 5 human haplotypes.

Biological Role of Unsaturated Fatty Acid Desaturases in Health and Disease

Polyunsaturated fatty acids (PUFAs) are considered one of the most important components of cells that influence normal development and function of many organisms, both eukaryotes and prokaryotes. Unsaturated fatty acid desaturases play a crucial role in the synthesis of PUFAs, inserting additional unsaturated bonds into the acyl chain. The level of expression and activity of different types of desaturases determines profiles of PUFAs. It is well recognized that qualitative and quantitative changes in the PUFA profile, resulting from alterations in the expression and activity of fatty acid desaturases, are associated with many pathological conditions. Understanding of underlying mechanisms of fatty acid desaturase activity and their functional modification will facilitate the development of novel therapeutic strategies in diseases associated with qualitative and quantitative disorders of PUFA.