Epigenetic Mechanisms and Posttranslational Modifications in Systemic Lupus Erythematosus

Research advances in protein lysine 2-hydroxyisobutyrylation: From mechanistic regulation to disease relevance

Histone lysine 2-hydroxyisobutyrylation (Khib) was identified as a novel posttranslational modification in 2014. Significant progress has been made in understanding its roles in reproduction, development, and disease. Although 2-hydroxyisobutyrylation shares some overlapping modification sites and regulatory factors with other lysine residue modifications, its unique structure suggests distinct functions. This review summarizes the latest advancements in Khib, including its regulatory mechanisms, roles in mammalian physiological processes, and its relationship with diseases. This provides direction for further research on Khib and offers new perspectives for developing treatment strategies for related diseases.

Proteomic analysis of lysine 2-hydroxyisobutyryl in SLE reveals protein modification alteration in complement and coagulation cascades and platelet activation Pathways

BACKGROUND

Post-translational modifications (PTMs) are considered to be an important factor in the pathogenesis of Systemic lupus erythematosus (SLE). Lysine 2-hydroxyisobutyryl (Khib), as an emerging post-translational modification of proteins, is involved in some important biological metabolic activities. However, there are poor studies on its correlation with diseases, especially SLE.

OBJECTIVE

We performed quantitative, comparative, and bioinformatic analysis of Khib proteins in Peripheral blood mononuclear cells (PBMCs) of SLE patients and PBMCs of healthy controls. Searching for pathways related to SLE disease progression and exploring the role of Khib in SLE.

METHODS

Khib levels in SLE patients and healthy controls were compared based on liquid chromatography tandem mass spectrometry, then proteomic analysis was conducted.

RESULTS

Compared with healthy controls, Khib in SLE patients was up-regulated at 865 sites of 416 proteins and down-regulated at 630 sites of 349 proteins. The site abundance, distribution and function of Khib protein were investigated further. Bioinformatics analysis showed that Complement and coagulation cascades and Platelet activation in immune-related pathways were significantly enriched, suggesting that differentially modified proteins among them may affect SLE.

CONCLUSION

Khib in PBMCs of SLE patients was significantly up- or down-regulated compared with healthy controls. Khib modification of key proteins in the Complement and coagulation cascades and Platelet activation pathways affects platelet activation and aggregation, coagulation functions in SLE patients. This result provides a new direction for the possible significance of Khib in the pathogenesis of SLE patients.

Advances in the potential roles of Cullin-RING ligases in regulating autoimmune diseases

Cullin-RING ligases (CRLs) are the largest class of E3 ubiquitin ligases regulating the stability and subsequent activity of a large number of important proteins responsible for the development and progression of various diseases, including autoimmune diseases (AIDs). However, the detailed mechanisms of the pathogenesis of AIDs are complicated and involve multiple signaling pathways. An in-depth understanding of the underlying regulatory mechanisms of the initiation and progression of AIDs will aid in the development of effective therapeutic strategies. CRLs play critical roles in regulating AIDs, partially by affecting the key inflammation-associated pathways such as NF-κB, JAK/STAT, and TGF-β. In this review, we summarize and discuss the potential roles of CRLs in the inflammatory signaling pathways and pathogenesis of AIDs. Furthermore, advances in the development of novel therapeutic strategies for AIDs through targeting CRLs are also highlighted.

The circRNA-miRNA-mRNA regulatory network in plasma and peripheral blood mononuclear cells and the potential associations with the pathogenesis of systemic lupus erythematosus

OBJECTIVES

This study aimed to explore the possible role of plasma and peripheral blood mononuclear cells (PBMCs) circular RNA (circRNA) in systemic lupus erythematosus (SLE).

METHOD

Total RNA was extracted from blood plasma samples obtained from 10 patients with SLE and 10 healthy controls and subjected to microarray analysis to define the profile of circRNA expression. The quantitative reverse transcription-polymerase chain reaction (qRT-PCR) amplification was conducted. The overlapped circRNA between PBMCs and plasma was performed, the interactions with microRNAs were predicted, the miRNA target mRNA was predicted, and the GEO database was used. The Gene ontology and pathway analysis was performed.

RESULTS

One hundred thirty-one upregulated and 314 significantly downregulated circRNAs were identified in the plasma of patients with SLE by the Fold change criteria (≥ 2.0) and P < 0.05. The qRT-PCR results showed that the expression of has-circRNA-102531, has-circRNA-103984, and has-circRNA-104262 was increased in plasma of SLE, and the expression of has-circRNA-102972, has-circRNA-102006, has-circRNA-104313 was decreased in plasma of SLE. Twenty-eight upregulated circRNAs and 119 downregulated circRNAs were overlapped from PBMCs and plasma, and ubiquitination was enriched. Furthermore, the circRNA-miRNA-mRNA network was constructed in SLE after analyzing dataset GSE61635 from GEO. The circRNA-miRNA-mRNA network comprises 54 circRNAs, 41 miRNAs, and 580 mRNAs. In addition, the TNF signaling pathway and the MAPK pathway were enriched from the mRNA of the miRNA target.

CONCLUSION

We first revealed the differentially expressed circRNAs in plasma and PBMCs, and then the circRNA-miRNA-mRNA network was constructed. The network's circRNAs could be a potential diagnostic biomarker and potentially play an important role in the pathogenesis and development of SLE. Key Points • This study analyzed the circRNAs expression profiles combined with the plasma and PBMCs, which provided a comprehensive overview of circRNAs expression patterns in SLE. • The network of the circRNA-miRNA-mRNA in SLE was constructed, which contributes to a better understanding of the pathogenesis and development of SLE.

Aberrant expression of JMJD3 in SLE promotes B-cell differentiation

Systemic lupus erythematosus (SLE) is a typical autoimmune disease distinguished by multiple organ dysfunction, which is related to a variety of causative factors. B-cell overactivation is a key factor in SLE. However, the pathogenesis underlying anomalous B cells has not been well elucidated. B-cell fate is regulated in diverse epigenetic ways apart from traditional ways. As one of the mechanisms of epigenetics, histone modification mainly affects transcription and translation by changing the chemical groups on histones by histone modification enzymes. JMJD3, a histone demethylase, can promote T-cell proliferation in SLE patients, which exacerbates SLE. However, the mechanism of JMJD3 in B cells in SLE has not been studied. Here, we found that the mean fluorescence intensity (MFI) of JMJD3 in classical memory B cells (CMBs) was higher than that in naïve B cells (NBs) from human tonsil tissue; JMJD3 was overexpressed in B cells from the peripheral blood of SLE patients compared with healthy controls (HCs). In vitro, our experiment showed that JMJD3 could regulate B-cell differentiation by promoting naïve B-cell differentiation into CD27⁺ B cells, and Blimp-1 and Bcl-6 also decreased after inhibitor treatment. These findings provide a new direction for the pathogenesis of SLE and may supply a new idea for subsequent drug development.

Classification of Systemic Lupus Erythematosus Using Raman Spectroscopy of Blood and Automated Computational Detection Methods: A Novel Tool for Future Diagnostic Testing

The aim of this study was to explore the proof of concept for using Raman spectroscopy as a diagnostic platform in the setting of systemic lupus erythematosus (SLE). We sought to identify unique Raman signatures in serum blood samples to successfully segregate SLE patients from healthy controls (HC). In addition, a retrospective audit was undertaken to assess the clinical utility of current testing platforms used to detect anti-double stranded DNA (dsDNA) antibodies (n = 600). We examined 234 Raman spectra to investigate key variances between SLE patients (n = 8) and HC (n = 4). Multi-variant analysis and classification model construction was achieved using principal component analysis (PCA), PCA-linear discriminant analysis and partial least squares-discriminant analysis (PLS-DA). We achieved the successful segregation of Raman spectra from SLE patients and healthy controls (p-value < 0.0001). Classification models built using PLS-DA demonstrated outstanding performance characteristics with 99% accuracy, 100% sensitivity and 99% specificity. Twelve statistically significant (p-value < 0.001) wavenumbers were identified as potential diagnostic spectral markers. Molecular assignments related to proteins and DNA demonstrated significant Raman intensity changes between SLE and HC groups. These wavenumbers may serve as future biomarkers and offer further insight into the pathogenesis of SLE. Our audit confirmed previously reported inconsistencies between two key methodologies used to detect anti-dsDNA, highlighting the need for improved laboratory testing for SLE. Raman spectroscopy has demonstrated powerful performance characteristics in this proof-of-concept study, setting the foundations for future translation into the clinical setting.

EZH2: Its regulation and roles in immune disturbance of SLE

The pathogenesis of systemic lupus erythematosus (SLE) is related to immune homeostasis imbalance. Epigenetic mechanisms have played a significant role in breaking immune tolerance. Enhancer of zeste homolog 2 (EZH2), the specific methylation transferase of lysine at position 27 of histone 3, is currently found to participate in the pathogenesis of SLE through affecting multiple components of the immune system. This review mainly expounds the mechanisms underlying EZH2-mediated disruption of immune homeostasis in SLE patients, hoping to provide new ideas in the pathogenesis of SLE and new targets for future treatment.

Autoantibodies against specific post-translationally modified proteins are present in patients with lupus and associate with major neuropsychiatric manifestations

Background

Although autoantibodies are an important hallmark of systemic lupus erythematosus (SLE), most are not specific for SLE or any of its clinical manifestations. Autoantibodies against post-translationally modified (PTM) proteins have been studied extensively in rheumatoid arthritis and associate with disease progression. While PTMs have also been detected in patients with SLE, studies on anti-PTM antibodies remain scarce. We studied the presence of anti-PTM antibodies in SLE and neuropsychiatric SLE (NPSLE), a manifestation that lacks serological markers.

Methods

IgG antibody responses against six PTMs (malondialdehyde–acetaldehyde adducts (MAA), advanced glycation end-products (AGE), carbamylation (CarP), citrullination, acetylation and nitration) were tested using ELISA in sera of 349 patients with SLE (mean age 44±13 years; 87% female) and compared with 108 healthy controls. Levels and positivity were correlated with clinical features and SLE manifestations.

Results

Anti-MAA, anti-AGE and anti-CarP antibodies were more prevalent in SLE compared with controls (MAA: 29% vs 3%, AGE: 18% vs 4%, CarP: 14% vs 5%, all p≤0.0001). Anti-MAA and anti-AGE antibodies correlated with clinical manifestations and serological inflammatory markers. Patients with major NPSLE showed higher positivity of anti-MAA (39% vs 24%, p=0.01) and anti-CarP antibodies (20% vs 11%, p=0.04) than patients without major NPSLE. In addition, anti-PTM antibody levels correlated with brain volumes, an objective measure of nervous system involvement.

Conclusions

In our NPSLE cohort, a subset of patients with SLE have anti-PTM antibodies against MAA, AGE and CarP modified proteins. Interestingly, anti-MAA and anti-CarP were more prevalent in NPSLE, a manifestation for which no biomarkers exist.

Docking, Binding Free Energy Calculations and In Vitro Characterization of Pyrazine Linked 2-Aminobenzamides as Novel Class I Histone Deacetylase (HDAC) Inhibitors

Class I histone deacetylases, HDAC1, HDAC2, and HDAC3, represent potential targets for cancer treatment. However, the development of isoform-selective drugs for these enzymes remains challenging due to their high sequence and structural similarity. In the current study, we applied a computational approach to predict the selectivity profile of developed inhibitors. Molecular docking followed by MD simulation and calculation of binding free energy was performed for a dataset of 2-aminobenzamides comprising 30 previously developed inhibitors. For each HDAC isoform, a significant correlation was found between the binding free energy values and in vitro inhibitory activities. The predictive accuracy and reliability of the best preforming models were assessed on an external test set of newly designed and synthesized inhibitors. The developed binding free-energy models are cost-effective methods and help to reduce the time required to prioritize compounds for further studies.

Belimumab Decreases Autophagy and Citrullination in Peripheral Blood Mononuclear Cells from Patients with Systemic Lupus Erythematosus

Belimumab (BLM) is a B lymphocyte stimulator (BLyS) inhibitor approved for the treatment of systemic lupus erythematosus (SLE). Autophagy is a cell survival mechanism involved in the pathogenesis of SLE. Citrullination is a post-translational modification catalyzed by peptidylarginine deiminase (PAD) enzymes. Autophagy and citrullination may generate neoepitopes, evoking an autoimmune response. No previous studies have investigated the connection of these processes, and how BLM could affect them, in SLE. Ex vivo autophagy and protein citrullination were analyzed by western blot in lysates from 26 SLE patients' PBMCs at baseline and after 2, 4, and 12 weeks of BLM administration, and from 16 healthy donors' PBMCs. Autophagic PBMCs were identified by the immunofluorescent detection of the autophagy-associated proteins LC3B (LC3 puncta) and LAMP-1. Autophagosome accumulation was evaluated in CD14- (PBLs) and CD14+ (monocytes) SLE cells. The presence of the BLyS receptors BAFF-R, BCMA, and TACI on SLE CD4+, CD8+ T cells and monocytes, as well as serum IL-18 levels, was also assessed. Following BLM administration, we observed a decrease in autophagy and citrullination, with a lowering of LC3-II, citrullinated vimentin, and PAD4 expression levels in PBMCs from SLE patients. LC3-II levels showed a correlation with the SLE Disease Activity Index 2000 (SLEDAI-2K) after 12 weeks of therapy. The LC3B/LAMP-1 analysis confirmed the reduction in autophagy. A lesser autophagosome accumulation occurred in PBLs and monocytes which, in turn, seemed to be the main cellular populations contributing to autophagy. A reduction in patients' serum IL-18 concentrations occurred. CD4+ and CD8+ cells weakly expressed BAFF receptors; monocytes expressed only BAFF-R. BLM could impact on autophagy and citrullination, offering an opportunity for a deeper understanding of these mechanisms in SLE, and a possible tool for the clinical management of SLE.

Association of macronutrients and dietary patterns with risk of systemic lupus erythematosus in the Black Women's Health Study

BACKGROUND

Systemic lupus erythematosus (SLE) affects African-American (AA) women disproportionately. The few prospective studies assessing dietary intake in relation to risk of SLE have been conducted in predominantly white populations and have been null.

OBJECTIVES

The present study assessed associations of macronutrients and dietary patterns with risk of SLE in AA women.

METHODS

Data from the Black Women's Health Study was collected prospectively via biennial questionnaires starting in 1995. Participants completed a self-administered 68-item FFQ in 1995. Self-reported SLE was verified through medical record review. We used multivariable (MV) Cox regression models to estimate HRs and 95% CIs for macronutrients, carbohydrates, proteins, total fats, PUFAs, ω-3 fatty acids, ω-6 fatty acids, MUFAs, saturated fats, trans fatty acids, Alternative Healthy Eating Index score, vegetable/fruit and meat/fried food dietary patterns, and a reduced rank regression (RRR)-derived dietary pattern in relation to SLE risk.

RESULTS

We confirmed a total of 114 incident cases of SLE among 51,934 women during 1995-2015. MVHRs and 95% CIs for the highest quintile of intake versus the lowest were HR: 1.96, 95% CI: 1.02, 3.67 for carbohydrates; HR: 0.66, 95% CI: 0.37, 1.18 for protein; and HR: 0.54, 95% CI: 0.28, 1.01 for total fats. MUFAs, saturated fatty acids, and trans fatty acids were significantly associated with a lower risk of SLE. An RRR-derived factor, rich in fruits and sugar-sweetened drinks and low in margarines and butter, red and processed meats, fried chicken, poultry, and eggs, which explained 53.4% of the total variation of macronutrients, was the only food pattern associated with increased SLE risk (HR: 1.88, 95% CI: 1.06, 3.35).

CONCLUSION

These analyses suggest that a diet high in carbohydrates and low in fats is associated with increased SLE risk in AA women.

Current paradigms in epigenetic anticancer therapeutics and future challenges

Any alteration at the genetic or epigenetic level, may result in multiplex of diseases including tumorigenesis which ultimately results in the cancer development. Restoration of the normal epigenome by reversing the epigenetic alterations have been reported in tumors paving the way for development of an effective epigenetic treatment in cancer. However, delineating various epigenetic events has been a challenging task so far despite substantial progress in understanding DNA methylation and histone modifications during transcription of genes. Many inhibitors in the form of epigenetic drugs mostly targeting chromatin and histone modifying enzymes including DNA methyltransferase (DNMT) enzyme inhibitors and a histone deacetylases (HDACs) inhibitor, have been in use subsequent to the approval by FDA for cancer treatment. Similarly, other inhibitory drugs, such as FK228, suberoylanilide hydroxamic acid (SAHA) and MS-275, have been successfully tested in clinical studies. Despite all these advancements, still we see a hazy view as far as a promising epigenetic anticancer therapy is concerned. The challenges are to have more specific and effective inhibitors with negligible side effects. Moreover, the alterations seen in tumors are not well understood for which one has to gain deeper insight into the tumor pathology as well. Current review focusses on such epigenetic alterations occurring in cancer and the effective strategies to utilize such alterations for potential therapeutic use and treatment in cancer.

Metabolomics bridging proteomics along metabolites/oncometabolites and protein modifications: Paving the way toward integrative multiomics

Systems biology adopted functional and integrative multiomics approaches enable to discover the whole set of interacting regulatory components such as genes, transcripts, proteins, metabolites, and metabolite dependent protein modifications. This interactome build up the midpoint of protein-protein/PTM, protein-DNA/RNA, and protein-metabolite network in a cell. As the key drivers in cellular metabolism, metabolites are precursors and regulators of protein post-translational modifications [PTMs] that affect protein diversity and functionality. The precisely orchestrated core pattern of metabolic networks refer to paradigm 'metabolites regulate PTMs, PTMs regulate enzymes, and enzymes modulate metabolites' through a multitude of feedback and feed-forward pathway loops. The concept represents a flawless PTM-metabolite-enzyme(protein) regulomics underlined in reprogramming cancer metabolism. Immense interconnectivity of those biomolecules in their spectacular network of intertwined metabolic pathways makes integrated proteomics and metabolomics an excellent opportunity, and the central component of integrative multiomics framework. It will therefore be of significant interest to integrate global proteome and PTM-based proteomics with metabolomics to achieve disease related altered levels of those molecules. Thereby, present update aims to highlight role and analysis of interacting metabolites/oncometabolites, and metabolite-regulated PTMs loop which may function as translational monitoring biomarkers along the reprogramming continuum of oncometabolism.

MicroRNAs overexpressed in Crohn's disease and their interactions with mechanisms of epigenetic regulation explain novel aspects of Crohn's disease pathogenesis

Background

In this review, we were interested to identify the wide universe of enzymes associated with epigenetic modifications, whose gene expression is regulated by miRNAs with a high relative abundance in Crohn's disease (CD) affected tissues, with the aim to determine their impact in the pathogenesis and evolution of the disease.

Methods

We used HMDD and Bibliometrix R-package in order to identify the miRNAs overexpressed in CD. The identified enzymes associated with epigenetic mechanisms and post-translational modifications, regulated by miRNAs upregulated in CD, were analyzed using String v11 database.

Results

We found 190 miRNAs with great abundance in patients with CD, of which 26 miRNAs regulate the gene expression of enzymes known to catalyze epigenetic modifications involved in essentials pathophysiological processes, such as chromatin architecture reorganization, immune response regulation including CD4+ T cells polarization, integrity of gut mucosa, gut microbiota composition and tumorigenesis.

Conclusion

The integrated analysis of miRNAs with a high relative abundance in patients with CD showed a combined and superimposed gene expression regulation of enzymes associated with relevant epigenetic mechanisms and that could explain, in part, the pathogenesis of CD.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01022-8.

Epigenetics in Necrotizing Enterocolitis

Epigenetic alterations in our genetic material can lead to heritable changes in the risk, clinical manifestations, course, and outcomes of many diseases. Understanding these epigenetic mechanisms can help in identifying potential therapeutic targets. This is especially important in necrotizing enterocolitis (NEC), where prenatal as well as postnatal factors impact susceptibility to this devastating condition, but our therapeutic options are limited. Developmental factors affecting intestinal structure and function, our immune system, gut microbiome, and postnatal enteral nutrition are all thought to play a prominent role in this disease. In this manuscript, we have reviewed the epigenetic mechanisms involved in NEC. These include key developmental changes in DNA methylation in the immature intestine, the role of long non-coding RNA (lncRNA) in maintaining intestinal barrier function, epigenetic influences of prenatal inflammation on immunological pathways in NEC pathogenesis such as Toll-Like Receptor 4 (TLR4) and epigenetic changes associated with enteral feeding causing upregulation of pro-inflammatory genes. We have assimilated research findings from our own laboratory with an extensive review of the literature utilizing key terms in multiple databases, including PubMed, EMBASE, and Science Direct.

Role of β2-adrenergic receptors in chronic obstructive pulmonary disease

Beta-2 adrenergic receptors (β2-ARs) have important roles in the pathogenesis and treatment of chronic obstructive pulmonary disease (COPD). In recent years, progress has been made in the study of β2-ARs. Here, we introduce the basic concepts of β2-ARs, related pathways, as well as application of blockers/agonists of β2-ARs, and β2-AR autoantibodies in COPD. Drugs targeting the β2-AR are being developed rapidly, and we expect them to improve the symptoms and prognosis of COPD patients in the future.

Aberrant Non-Coding RNA Expression in Patients with Systemic Lupus Erythematosus: Consequences for Immune Dysfunctions and Tissue Damage

Systemic lupus erythematosus (SLE) is a complex systemic autoimmune disease with heterogeneous clinical manifestations. A diverse innate and adaptive immune dysregulation is involved in the immunopathogenesis of SLE. The dysregulation of immune-related cells may derive from the intricate interactions among genetic, epigenetic, environmental, and immunological factors. Of these contributing factors, non-coding RNAs (ncRNAs), including microRNAs (miRNAs, miRs), and long non-coding RNAs (lncRNAs) play critical roles in the post-transcriptional mRNA expression of cytokines, chemokines, and growth factors, which are essential for immune modulation. In the present review, we emphasize the roles of ncRNA expression in the immune-related cells and cell-free plasma, urine, and tissues contributing to the immunopathogenesis and tissue damage in SLE. In addition, the circular RNAs (circRNA) and their post-translational regulation of protein synthesis in SLE are also briefly described. We wish these critical reviews would be useful in the search for biomarkers/biosignatures and novel therapeutic strategies for SLE patients in the future.

Systemic lupus erythematosus: an expert insight into emerging therapy agents in preclinical and early clinical development

INTRODUCTION

Systemic lupus erythematosus (SLE) is a chronic disease that is potentially fatal. There is no cure for SLE and the medications used are associated with toxic side effects. In the era of revolutionary emerging novel biologic agents, the design and investigation of targeted therapy for these patients is necessary. Novel therapies under investigation in phase II-III clinical trials showed promising results. Therapies can target various pathways involved in SLE including cytokines, signal transduction inhibitors, B-cell depletion and interference with co-stimulation. Of interest is the proof of concept of sequential therapy.

AREAS COVERED

We performed an extensive literature search via PubMed, Medline, Elsevier Science and Springer Link databases between the years 2014-2020 using the following terms: SLE, novel treatments. We have reviewed 232 articles and selected those articles that (i) focus on phase II-III emerging therapies and (ii) offer new findings from existing therapies, which reveal breakthrough concepts in SLE treatment.

EXPERT OPINION

It is still difficult to crack the puzzle of a successful SLE treatment approach. New strategies with potential may encompass the targeting of more than one protein. Another way forward is to identify each SLE patient and personalize therapy by clinical manifestations, disease activity, serology and activated protein.

An integrative Bayesian network approach to highlight key drivers in systemic lupus erythematosus

BACKGROUND

A comprehensive intuition of the systemic lupus erythematosus (SLE), as a complex and multifactorial disease, is a biological challenge. Dealing with this challenge needs employing sophisticated bioinformatics algorithms to discover the unknown aspects. This study aimed to underscore key molecular characteristics of SLE pathogenesis, which may serve as effective targets for therapeutic intervention.

METHODS

In the present study, the human peripheral blood mononuclear cell (PBMC) microarray datasets (n = 6), generated by three platforms, which included SLE patients (n = 220) and healthy control samples (n = 135) were collected. Across each platform, we integrated the datasets by cross-platform normalization (CPN). Subsequently, through BNrich method, the structures of Bayesian networks (BNs) were extracted from KEGG-indexed SLE, TCR, and BCR signaling pathways; the values of the node (gene) and edge (intergenic relationships) parameters were estimated within each integrated datasets. Parameters with the FDR < 0.05 were considered significant. Finally, a mixture model was performed to decipher the signaling pathway alterations in the SLE patients compared to healthy controls.

RESULTS

In the SLE signaling pathway, we identified the dysregulation of several nodes involved in the (1) clearance mechanism (SSB, MACROH2A2, TRIM21, H2AX, and C1Q gene family), (2) autoantigen presentation by MHCII (HLA gene family, CD80, IL10, TNF, and CD86), and (3) end-organ damage (FCGR1A, ELANE, and FCGR2A). As a remarkable finding, we demonstrated significant perturbation in CD80 and CD86 to CD28, CD40LG to CD40, C1QA and C1R to C2, and C1S to C4A edges. Moreover, we not only replicated previous studies regarding alterations of subnetworks involved in TCR and BCR signaling pathways (PI3K/AKT, MAPK, VAV gene family, AP-1 transcription factor) but also distinguished several significant edges between genes (PPP3 to NFATC gene families). Our findings unprecedentedly showed that different parameter values assign to the same node based on the pathway topology (the PIK3CB parameter values were 1.7 in TCR vs - 0.5 in BCR signaling pathway).

CONCLUSIONS

Applying the BNrich as a hybridized network construction method, we highlight under-appreciated systemic alterations of SLE, TCR, and BCR signaling pathways in SLE. Consequently, having such a systems biology approach opens new insights into the context of multifactorial disorders.