Skip to main content
Download PDF

Causal relationship between systemic lupus erythematosus and primary liver cirrhosis based on two-sample bidirectional Mendelian randomization and transcriptome overlap analysis

Arthritis Research & Therapy volume 26, Article number: 10 (2024) Cite this article

Abstract

Background

Overlapping cases of systemic lupus erythematosus (SLE) and primary biliary cirrhosis (PBC) are rare and have not yet been fully proven to be accidental or have a common genetic basis.

Methods

Two-sample bidirectional Mendelian randomization (MR) analysis was applied to explore the potential causal relationship between SLE and PBC. The heterogeneity and reliability of MR analysis were evaluated through Cochran’s Q-test and sensitivity test, respectively. Next, transcriptome overlap analysis of SLE and PBC was performed using the Gene Expression Omnibus database to identify the potential mechanism of hub genes. Finally, based on MR analysis, the potential causal relationship between hub genes and SLE or PBC was validated again.

Results

The MR analysis results indicated that SLE and PBC were both high-risk factors for the occurrence and development of the other party. On the one hand, MR analysis had heterogeneity, and on the other hand, it also had robustness. Nine hub genes were identified through transcriptome overlap analysis, and machine learning algorithms were used to verify their high recognition efficiency for SLE patients. Finally, based on MR analysis, it was verified that there was no potential causal relationship between the central gene SOCS3 and SLE, but it was a high-risk factor for the potential risk of PBC.

Conclusion

The two-sample bidirectional MR analysis revealed that SLE and PBC were high-risk factors for each other, indicating that they had similar genetic bases, which could to some extent overcome the limitation of insufficient overlap in case samples of SLE and PBC. The analysis of transcriptome overlapping hub genes provided a theoretical basis for the potential mechanisms and therapeutic targets of SLE with PBC overlapping cases.

Background

Autoimmune diseases are all caused by the loss of immune tolerance, and genetic susceptibility is related to the loss of tolerance and autoimmunity. Autoimmune diseases may have similar genetic susceptibility, and multiple autoimmune diseases can occur in the same patient [1]. Systemic lupus erythematosus (SLE) and primary biliary (PBC) cirrhosis are both chronic autoimmune diseases.

SLE is an immune disease with multiple system involvement, recurrent attacks, diverse clinical manifestations, and high incidence rate in women. SLE is also often accompanied by other immune diseases. For example, the prevalence of SLE with Sjogren’s syndrome exceeded 20% based on two different assessment criteria [2]. The liver is an important target organ involved in multiple systems of SLE. Liver involvement in SLE can also lead to liver dysfunction, which can further lead to cirrhosis and liver cancer. However, the diagnosis and treatment of liver involvement in SLE may not be timely, which masks the progression of liver involvement. In addition, SLE inhibition therapy may also mask the progression of primary liver disease. There are still difficulties in accurately distinguishing SLE from primary liver diseases that are difficult to overcome.

PBC is a primary autoimmune disease of the liver that is more common in women, characterized by cholestasis, bile duct injury, cholangitis, and liver fibrosis, which can also lead to liver dysfunction, cirrhosis, or liver cancer [3]. PBC is also accompanied by extrahepatic autoimmune diseases. For example, a systematic review found that the prevalence of PBC associated with Sjogren’s syndrome ranged from 3.5 to 73%. However, the low incidence of PBC has rendered many existing treatments ineffective for PBC. If combined with extrahepatic autoimmune diseases, this will pose greater challenges for the precise treatment of PBC [4].

The overview of SLE and PBC mentioned above indicates that there are many similarities between the two, so the two may be overlapping diseases. However, there are relatively few overlapping cases of SLE and PBC, with a reported prevalence ranging from 0.5 to 3.7% [5].

The independent occurrence of SLE or PBC was found to be related to genetics. For example, based on genome-wide association studies, 38 new SLE-related loci and incompletely shared genetic structures were identified in a large sample population in China [6]; 50 significant loci were identified across the entire genome of 10,516 PBC patients and 20,772 healthy individuals in 6 countries, and potential gene targeting pathways for new loci were identified [7]. However, the occurrence of overlapping cases of SLE and PBC has not yet been fully demonstrated as an accidental event or having a common genetic basis [8].

On this basis, from the perspective of clinical characteristics, a total of 34 overlapping cases of SLE and PBC publicly published before 2015 were analyzed. It was found that they were still more common in women, mostly in the middle-aged and elderly. PBC often occurred before SLE, and some cases were combined with multiple immune diseases (such as Sjogren’s syndrome and hepatocellular carcinoma) [8]. Unfortunately, there is still a lack of research based on genetics to verify the causal relationship between overlapping cases of SLE and PBC.

Mendelian randomization study (MR) is defined as using single nucleotide polymorphisms (SNPs) as instrumental variables to analyze the causal relationship between exposure factors and disease outcomes [9]. Based on a large sample genome-wide association study (GWAS) database, suitable SNPs were screened to effectively reduce the bias of confounding factors on disease outcomes, which were widely used in exploring the causal relationship between phenotypes and diseases. For example, two-sample MR studies evaluated the relationship between SLE and celiac disease and found that SLE could significantly increase the risk of celiac disease [10]; In the European population, inflammatory bowel disease predicted based on 99 SNPs increases the risk of PBC [11]. These and other studies validate the causal relationship between immune diseases at the genetic level, unlike clinical observational studies, which reduce the impact of many confounding factors. Transcriptomics analysis of overlapping diseases, further analyzing the hub genes most likely involved in the development of SLE and PBC overlap, and using hub genes to identify and design drugs that can be used to treat SLE and PBC. For example, the overlapping genes of SLE and venous thromboembolism were identified and found to be closely related to the immune microenvironment, which providing clues for further exploration of the common mechanisms and interactions between the two diseases [12].

Therefore, in this study, firstly, the systematic evaluation of overlapping cases of SLE and PBC would be re-investigated to evaluate the prevalence of overlapping cases. Subsequently, based on genetic evidence from bidirectional dual samples, MR analysis revealed a causal relationship between SLE and PBC. Then, in the analysis of transcriptome data from two samples, potential overlapping genes between SLE and PBC were identified, and optimal models for identifying SLE were constructed based on 11 machine learning algorithms. Finally, based on MR analysis of the causal relationship between overlapping hub genes and SLE/PBC, the robustness of hub genes was verified again. These research contents will provide a deep understanding of the potential common pathophysiological processes of SLE with PBC, and provide a theoretical basis for the implementation of precise diagnosis and treatment in the future.

Methods

Literature retrieval of SLE with PBC

This study was searched in the PubMed database using the following search term combinations: [systemic lupus erythematosus or SLE] and [primary biliary cirrhosis or primary cholangitis or PBC] to extract information in the abstract or full text: author, year, and cases.

Data retrieval

IEU OPEN GWAS database (https://gwas.mrcieu.ac.uk/) was used to retrieve exposure and outcome data. SLE data were obtained from the GWAS dataset (GWAS ID: ebi-a-GCST003156), which dataset was analyzed and constructed by Bentham J et al., consisting of 14,267 Europeans (5201 cases and 9066 controls) with 7,071,163 SNPs [13]. PBC data also were obtained from the GWAS dataset (GWAS ID: ebi-a-GCST005581), which dataset was analyzed and constructed by Liu JZ et al., consisting of 11,375 Europeans (2861 cases and 8514 controls) with 119,756 SNPs [14].

Genetic instrumental analysis

In order to ensure complete random independence of instrumental variable estimation and exclude the impact of linkage disequilibrium (LD) on outcome factors, a threshold (P < 5e − 8) was set to screen SNPs closely related to exposure factors in a genome-wide sense. Next, when the LD parameter was set (r2 = 0.001 and kb = 10,000), the above SNP data would be filtered for preprocessing to maintain consistency between the exposure effect allele and the outcome effect magnitude. In addition, the F-statistic was calculated to estimate overlapping effects and weak instrumental bias, and F < 10 was considered a suspicious bias and excluded.

MR analysis

Based on the TwoSampleMR software package, the inverse variance weighted (IVW) method was used as the main analysis to evaluate the causal relationship between SLE and PBC. In addition, MR Egger regression, weighted median method, simple estimation based on mode, and weighted estimation based on mode were used as auxiliary analysis methods. Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO), left one sensitivity analysis, and Cochran’s Q-test were used to extensively evaluate MR analysis results. The MR-PRESSO test was used to examine the potential skewness of horizontal pleiotropic by identifying and excluding pleiotropic SNPs with P < 0.05 to reassess causal relationships. The leave-one-out method was used to determine the robustness of MR results by excluding SNPs. Cochran’s Q-test determined the heterogeneity leading to MR results based on IVW and MR Egger estimates [15].

Transcriptome in causal relationship validation between SLE and PBC

The microarray expression and clinical data of GSE65391 (924 SLEs and 72 controls) and GSE93170 (6 PBCs and 6 controls) were downloaded from the Gene Expression Omnibus (GEO) database. Based on the limma software package, overlapping differential genes (DEGs) of SLE and PBC were identified, with a threshold of log2 FC > 0 and P < 0.05. The biological mechanism of overlapping DEGs based on kyoto encyclopedia of genes and genomes (KEGG) annotations. Based on the STRING database (https://string-db.org), fifteen key hub genes were identified and a PPI network was constructed. On the basis of the 15 bub genes mentioned above, models for identifying SLE were jointly developed based on the least absolute shrinkage and selection operator (LASSO) dimensionality reduction using 11 machine learning algorithms [16]. Further validate the expression of 15 hub genes in SLE patients based on single-cell data downloaded from GEO (GSE163121, GSE162577). The integration and removal of batch effects of diverse materials were processed using the “Harmony” R package. The single-cell analysis process adopts the Seurat method, and the marker genes referred to previous study. The difference analysis between cell clusters was performed using the Kruskal–Wallis H method.

Statistical analysis

R software (version 4.2.2) for data analysis and plotting. The R software package mainly includes “VariantAnnotation”, “gwasglue”, and “TwoSampleMR”. The estimated value of the causal correlation effect was represented by the odds ratio (OR) and 95% confidence interval (CI) of the risk change of the output factor responding to each increase in standard deviation (SD) of the exposure factor, and the test level was α = 0.05. The classification performance of the models was evaluated based on sensitivity, specificity, accuracy, receiver operating characteristic curve (ROC), and area under ROC curve (AUC).

Results

Literature retrieval of SLE with PBC

This study included a total of 9 studies and extracted 125 SLE with PBC (Table 1), including 3 from USA (2.4%), 28 from Turkey (22.4%), 37 from Japan (29.6%), 8 from Italy (6.4%), and 49 from China (39.2%) [5, 8, 17,18,19,20,21,22,23].

Table 1 SLE with PBC cases derived from the literature
Full size table

Genetic liability to PBC with SLE

After genetic instrumental analysis, 44 SNPs of SLE were used for MR analysis. Two pleiotropic SNPs (rs35000415 and rs389884) were identified through MR-PRESSO test to evaluate horizontal pleiotropy. After removing outliers from the SNPs, only 16 SNPs were found to match the SNPs for PBC. A significant causal relationship was observed between SLE and PBC in the IVW main analysis (OR 1.308, 95% CI, [1.131, 1.513], P < 0.001). In addition, the weighted median (OR 1.390, 95% CI, [1.265, 1.523], P < 0.001), simple estimation based on mode (OR 1.435, 95% CI [1.231, 1.673], P < 0.001), and weighted estimation (OR 1.408, 95% CI [1.218, 1.626], P < 0.001) also showed the same results as IVW (Table 2, Fig. 1A). Cochran’s Q-test showed significant statistical significance (IVW, P < 0.001, MR Egger, P < 0.001), indicating heterogeneity in MR analysis results (Fig. 1C). The leave-one-out method was conducted to verify the robustness of the results found (P = 0.124) (Fig. 1E).

Fig. 1
figure 1

MR estimation of genetic reliability between SLE and PBC. (A-B) Scatter plots. The estimation effect of MR methods with slope. Bilateral MR analysis showed that both IVW methods were high-risk factors for SLE or PBC. (C-D) Funnel plots. Cochran Q-test showed heterogeneity in the results of bidirectional MR analysis. (E-F) Leave-one-out plots. The black dots were used by the IVW method to evaluate causal effects and excluded individual analyses one by one. The red dot indicates the use of IVW for all SNPs. Bidirectional sensitivity testing verified the robustness of MR analysis

Full size image
Table 2 MR estimation of genetic reliability between SLE and PBC
Full size table

Genetic liability to SLE with PBC

After genetic instrumental analysis, 22 SNPs of PBC were used for MR analysis. Seven pleiotropic SNPs (rs11117431, rs2523882, rs3135024, rs34725611, rs35188261, rs7775055 and rs79513546) were identified through MR-PRESSO test to evaluate horizontal pleiotropy. After removing outliers from the SNPs, 22 SNPs were found to match the SNPs for SLE. A significant causal relationship was observed between SLE and PBC in the IVW main analysis (OR 1.407, 95% CI [1.199, 1.652], P < 0.001). In addition, the weighted median (OR 1.271, 95% CI [1.161, 1.392], P < 0.001), simple estimation based on mode (OR 1.363, 95% CI [1.174, 1.582], P = 0.001), and weighted estimation (OR 1.247, 95% CI [1.101, 1.412], P = 0.002) also showed the same results as IVW (Table 2, Fig. 1B). Cochran’s Q-test showed significant statistical significance (IVW, P < 0.001, MR Egger, P < 0.001), indicating heterogeneity in MR analysis results (Fig. 1D). The leave-one-out method was conducted to verify the robustness of the results found (P = 0.169) (Fig. 1F).

Transcriptome in causal relationship validation between SLE and PBC

In the differential analysis of SLE (GSE65391) and PBC (GSE93170) data, 418 overlapping DEGs were identified. The potential functions of overlapping DEGs had been explored by KEGG, mainly involving the NOD-like receiver signaling pathway, FoxO signaling pathway, cell cycle, TNF signaling pathway, and JAK-STAT signaling pathway. On the basis of the 418 DEGs mentioned above, the PPI network was constructed and imported into the Cytoscape’s Cytohubba plugin to calculate the top 15 hub genes (DDX60, EIF2AK2, IRF9, IFITM3, STAT1, IFITM1, OAS3, HERC5, IFI6, IFI44, IFIT5, IFNGR1, ISG15, STAT3, SOCS3, Fig. 2A). The GSE65391 dataset was divided into a training cohort (n = 498) and a validation cohort (n = 498) in a ratio of 5 to 5. After LASSO further selection of hub genes, 9 genes (DDX60, EIF2AK2, IRF9, IFITM3, STAT1, IFI6, IFNGR1, ISG15, SOCS3, Fig. 2B–C) were used to construct models for identifying SLE through 11 machine learning algorithms. The results of the training cohort and validation cohort showed high discriminative performance (AUC > 0.90, accuracy > 0.87, sensitivity > 0.87, specificity > 0.87, Table 3). Through comprehensive analysis, it was found that the SVM machine learning algorithm was superior in terms of fitting degree and evaluation reliability. The AUC, accuracy, sensitivity, and specificity of the training and validation queues were 0.96 and 0.94, 0.96 and 0.97, 0.96 and 0.97, 0.93 and 0.87, respectively (Fig. 2D). Surprisingly, through single-cell data analysis, among the 15 hub genes, EIF2AK2, IFITM3, STAT1, IFITM1, IFI6, IFI44, IFNGR1, ISG15, and STAT3 were highly expressed in various cells (Endothecal cells, Erythroid cells, Monocytes, NK cells, Platelets). Unfortunately, hub genes were not significantly expressed in B cells, with some being highly expressed in T cells (Fig. 3).

Fig. 2
figure 2

Transcriptome in causal relationship validation between SLE and PBC. A Protein–protein interaction network of the top 15 hub genes. B–C LASSO regression was used for machine learning gene set selection. Coefficient profile diagram (B) and cross-validation diagram (C). D ROC curve of SVM models

Full size image
Table 3 Models for identifying SLE patients using 9 hub genes based on 11 machine learning algorithms
Full size table
Fig. 3
figure 3

Expression analysis of 15 hub genes in single-cell data. Through single-cell data analysis, among the 15 hub genes, EIF2AK2, IFITM3, STAT1, IFITM1, IFI6, IFI44, IFNGR1, ISG15, STAT3 were highly expressed in various cells (endothecal cells, erythroid cells, monocytes, NK cells, platelets)

Full size image

SOCS3 genetic reliability for SLE and PBC

Due to the low significance of SNPs (P < 5e − 8) in the SOCS3 (GWAS ID: prot-a-2797) dataset, the new threshold was set to P < 5e − 5. In the SOCS3’s genetic reliability for SLE, Cochran’s Q-test showed no significant statistical significance (IVW, P = 0.513, MR Egger, P = 0.536), indicating no heterogeneity in MR analysis results. In the IVW main analysis, there was no significant causal relationship between SOCS3 and SLE (OR 0.992, 95% CI [0.936, 1.051], P = 0.787), while other methods showed the same results. The leave-one-out method verified the robustness of the results (P = 0.207). In the SOCS3’s genetic reliability for PBC, Cochran’s Q-test did not show significant statistical significance (IVW, P = 0.916, MR Egger, P = 0.906), indicating no heterogeneity in MR analysis results. In the IVW main analysis, there was a significant causal relationship between SLE and PBC (OR 1.387, 95% CI [1.049, 1.835], P = 0.022, Table 4), while other methods showed the same results without statistical significance. The leave-one-out method verified the robustness of the results (P = 0.632).

Table 4 MR estimation of genetic reliability between SOCS3 and SLE/PBC
Full size table

Discussion

At present, the clinical characteristics and pathophysiology of SLE with PBC overlap cases are not clear. Further analysis is needed to determine whether the overlap between the two will affect the risk of SLE or PBC occurrence and development, as well as whether it interferes with treatment effectiveness. This study investigated the epidemiological characteristics of SLE with PBC, which was similar to literature reports and had a low incidence rate. In addition, the potential causal relationship between SLE with PBC was elucidated from a genetic perspective. The MR analysis results showed that both SLE and PBC were high-risk factors for the occurrence and development of the other party. In addition, hub genes were identified through transcriptome overlap analysis, and their roles in SLE were validated through machine learning algorithms. Finally, based on MR analysis, the causal relationship between hub gene SOCS3 and SLE/PBC was verified. This study provides a new perspective for the precise diagnosis and treatment of SLE with PBC overlapping cases and also provides a theoretical basis for future research on common potential mechanisms and potential therapeutic targets.

Numerous well-known, the prevalence of SLE with PBC is low. According to statistics, a total of 125 cases had been reported, most of which were case reports, with one-third of the cases coming from China, and 75 cases had been reported in the past three years. The clinical characteristics were still the same as previous findings, and were more common in middle-aged and elderly women, with diverse clinical manifestations. However, observational studies had shown that SLE with PBC may have potential independent clinical features compared to PBC, such as anti-mitochondrial antibodies, hemoglobin, and albumin. In addition, SLE with PBC was a high-risk factor for adverse events related to liver cirrhosis and liver cancer. For example, univariate analysis showed that SLE with PBC was associated with the development of HCC (hazard ratio = 5.02) [23]. Unfortunately, no significant difference in survival was found between SLE with PBC and PBC alone [5]. However, the application of observational studies on confounding factors of diversity and reverse causal effects in causal inference. To reduce confounding factors, observational studies can typically use propensity scoring methods to match baseline characteristics of exposed and unexposed patients, resulting in effect estimates similar to those of randomized controlled trials [24]. Unfortunately, matching queues during propensity score analysis results in a large amount of data loss, which may weaken the reliability of the conclusion [25]. Therefore, more new technologies are still needed to explore the causal relationship between SLE and PBC.

MR is a genetic variable analysis that follows mendelian laws of inheritance, allowing SNPs to infer the causal relationship between exposure factors and clinical outcomes. Due to the random segregation of alleles, MR can reduce bias caused by confounding factors. In addition, since genetic variation occurs before the disease, the order of the two cannot be reversed, MR can also avoid interference from reverse causal relationships [26]. GWAS is a database of MR analysis sources that can be used to indicate that SNPs present at multiple locations in the genome can be associated with specific phenotypes [27]. In this study, the genetic causal relationship between SLE and PBC was elucidated in a dual sample bidirectional manner based on large sample whole genome association analysis. In this study, there were 15674 significant loci in SLE and 1135 significant loci in PBC. However, after linkage disequilibrium analysis, there were 44 remaining SNPs in SLE and 22 remaining SNPs in PBC, and these SNPs did not overlap. The use of IVW as the main analysis method not only indicated a causal relationship between the increased risk of PBC and SLE (OR 1.308, 95% CI [1.131, 1.513], P < 0.001), but also indicated a causal relationship between the increased risk of SLE and PBC (OR 1.407, 95% CI [1.199, 1.652], P < 0.001). In addition, the three methods of weighted media, simple estimation based on mode, and weighted estimation also supported the analysis results of IVW. However, the results of bidirectional MR analysis showed significant heterogeneity, which may be related to the number of SNPs detected. Although the sample size of SLE and PBC was similar and both were European populations, there was a significant difference in the number of SNPs detected between the two (SLE, 7071163; PBC, 119756), which led to a large number of SNPs in SLE whose genetic value had not been elucidated. Excitedly, sensitivity testing supported the stability and accuracy of causal results between SLE and PBC.

Transcriptome overlap analysis is a genetic validation method for analyzing the relationship between two diseases. For example, based on three datasets of polycystic ovary syndrome and recurrent implantation failure, 12 disease overlapping genes were identified using weighted gene co expression networks, functional enrichment analysis, and three machine learning algorithms, which potential pathogenic associations and mechanisms between the two were explored [28]; The inflammatory immune pathway of calcified aortic valve disease (CAVD) related to chronic kidney disease (CKD) had been revealed, providing new insights for future serum based diagnosis and treatment interventions for CKD combined with CAVD [29]. In this study, 418 upregulated overlapping genes were identified using the GEO dataset of SLE and PBC, and their biological functions were found to be related to TNF signaling pathway and JAK-STAT signaling pathway, which was similar to previous reports [7]. Subsequently, based on the PPI network, the top 15 most valuable hub genes were identified. Through single-cell data analysis, multiple hub genes had also been validated to exhibit high expression in various cells. After further screening by LASSO, 9 genes (DDX60, EIF2AK2, IRF9, IFITM3, STAT1, IFI6, IFNGR1, ISG15, SOCS3) were selected and presented to 11 machine learning algorithms to validate their value for SLE. The results showed that the machine learning models constructed by these 9 genes were effective in distinguishing SLE patients from non-SLE patients, with high efficiency in AUC, accuracy, sensitivity, and specificity. In addition, among the 9 genes mentioned above, 6 genes (DDX60, EIF2AK2, IRF9, IFI6, IFNGR1, ISG15) have been elucidated to be associated with SLE, but no reports have been found to be related to PBC; Both STAT1 and SOCS3 were found to be associated with SLE and PBC; IFITM3 has not been reported to be associated with SLE or PBC. SOCS3 could be lost by the IL10R1 gene allele, inducing genetic susceptibility to SLE disease in Caucasian populations [30]; The expression of SOCS3 mRNA was more inhibited in the PBC liver, but it was unlikely to explain its loss of Stat3 DNA binding in PBC [31]. The causal relationship between the SOCS3 gene and SLE and PBC will be further confirmed based on MR analysis. The results showed that elevated levels of the SOCS3 gene were a high-risk factor for increased risk of PBC (OR 1.387, 95% CI [1.049, 1.835], P = 0.022), and there was no causal relationship with SLE (OR 0.992, 95% CI [0.936, 1.051], P = 0.787).

This study had some advantages and limitations. The advantages of this study: Firstly, MR analysis could evaluate the bidirectional causal relationship between SLE and PBC in both directions and reduce confounding factors. Then, a novel method based on MR analysis was used to explore the clinical prognostic characteristics of SLE with PBC overlapping cases, which could to some extent overcome the limitation of insufficient case samples. Finally, unlike other MR studies, this study used transcriptome overlap analysis to re-validate the potential common mechanism between the two. Limitations of this study: Firstly, both population groups were European, and although the sample size was similar, there was a significant difference in the number of SNPs tested. Subsequently, the representativeness of SNPs’ significance was insufficient, and cross-validation of SNPs’ significance was still needed in more samples. Then, MR analysis exhibits heterogeneity, which may be related to significant differences in SNPs between the two samples and insufficient utilization of SNPs information. Thirdly, similar to other studies, based on the summary data analysis used by GWAS, individual data could not be obtained and other confounding factors could not be excluded. Fourthly, the microarray expression sequencing platforms for SLE and PBC downloaded from GEO were different, and there was a significant difference in sample size, which was also an important limitation. Finally, sensitivity analysis could not completely eliminate potential level pleiotropy [32, 33].

Conclusions

This study is the first attempt to use MR analysis to explore the bidirectional causal relationship between SLE and PBC, which can to some extent overcome the limitation of insufficient overlapping case samples in SLE with PBC. The results showed that both SLE and PBC were associated with an increased risk of development in the other party. The analysis of transcriptome overlap hub genes provides a new perspective for the precise diagnosis and treatment of SLE with PBC overlap cases and also provides a theoretical basis for future research on common potential mechanisms and potential therapeutic targets.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

SLE:

Systemic lupus erythematosus

PBC:

Primary biliary cirrhosis

MR:

Mendelian randomization

SNPs:

Single nucleotide polymorphisms

GWAS:

Genome-wide association study

IVW:

Inverse variance weighted

MR-PRESSO:

Mendelian randomization pleiotropy residual sum and outlier

GEO:

Gene expression omnibus

DEG:

Overlapping differential gene

GO:

Gene ontology

KEGG:

Kyoto encyclopedia of genes and genomes

LASSO:

Least absolute shrinkage and selection operator

OR:

Odds ratio

SD:

Standard deviation

ROC:

Receiver operating characteristic curve

AUC:

Area under curve

CI:

Confidence interval

CVAD:

Calcified aortic valve disease

CKD:

Chronic kidney disease

LD:

Linkage disequilibrium

References

  1. Zhang L, Yu C, Ren C. Systemic lupus erythematosus combined with Evans syndrome and primary biliary cirrhosis: a rare case report. Asian J Surg. 2023;S1015–9584(23):00850–3.

    Google Scholar 

  2. Gianordoli APE, Laguardia RVRB, Santos MCFS, et al. Prevalence of Sjögren’s syndrome according to 2016 ACR-EULAR classification criteria in patients with systemic lupus erythematosus. Adv Rheumatol. 2023;63(1):11.

    Article  PubMed  Google Scholar 

  3. Yang Y, He X, Rojas M, et al. Mechanism-based target therapy in primary biliary cholangitis: opportunities before liver cirrhosis? Front Immunol. 2023;14:1184252.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Colapietro F, Bertazzoni A, Lleo A. Contemporary Epidemiology of Primary Biliary Cholangitis. Clin Liver Dis. 2022;26(4):555–70.

    Article  PubMed  Google Scholar 

  5. Fan X, Men R, Ni P, et al. Concomitant systemic lupus erythematosus might have a negative impact on the biochemical responses to treatment in patients with primary biliary cholangitis. Clin Rheumatol. 2020;39(3):795–803.

    Article  PubMed  Google Scholar 

  6. Wang YF, Zhang Y, Lin Z, et al. Identification of 38 novel loci for systemic lupus erythematosus and genetic heterogeneity between ancestral groups. Nat Commun. 2021;12(1):772.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Cordell HJ, Fryett JJ, Ueno K, et al. An international genome-wide meta-analysis of primary biliary cholangitis: Novel risk loci and candidate drugs. J Hepatol. 2021;75(3):572–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Shizuma T. Clinical characteristics of concomitant systemic lupus erythematosus and primary biliary cirrhosis: a literature review. J Immunol Res. 2015;2015:713728.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Birney E. Mendelian randomization. Cold Spring Harb Perspect Med. 2022;12(4):a041302.

    PubMed  Google Scholar 

  10. Wang Q, Jia S, Lu Q, et al. Systemic lupus erythematosus is associated with the risk of coeliac disease: a Mendelian randomisation study. Autoimmunity. 2023;56(1):2250103.

    Article  PubMed  Google Scholar 

  11. Zhang H, Chen L, Fan Z, et al. The causal effects of inflammatory bowel disease on primary biliary cholangitis: A bidirectional two-sample Mendelian randomization study. Liver Int. 2023;43(8):1741–8.

    Article  CAS  PubMed  Google Scholar 

  12. Yu J, Yang J, He Q, et al. Comprehensive bioinformatics analysis reveals the crosstalk genes and immune relationship between the systemic lupus erythematosus and venous thromboembolism. Front Immunol. 2023;14:1196064.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Bentham J, Morris DL, Graham DSC, et al. Genetic association analyses implicate aberrant regulation of innate and adaptive immunity genes in the pathogenesis of systemic lupus erythematosus. Nat Genet. 2015;47(12):1457–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Liu JZ, Almarri MA, Gaffney DJ, et al. Dense fine-mapping study identifies new susceptibility loci for primary biliary cirrhosis. Nat Genet. 2012;44(10):1137–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Wei Y, Li Z, Lai H, et al. Instant coffee is negatively associated with telomere length: finding from observational and Mendelian randomization analyses of UK biobank. Nutrients. 2023;15(6):1354.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Bi Q, Wang Y, Deng Y, et al. Different multiparametric MRI-based radiomics models for differentiating stage IA endometrial cancer from benign endometrial lesions: A multicenter study. Front Oncol. 2022;12:939930.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Chowdhary VR, Crowson CS, Poterucha JJ, et al. Liver involvement in systemic lupus erythematosus: case review of 40 patients. J Rheumatol. 2008;35(11):2159–64.

    Article  CAS  PubMed  Google Scholar 

  18. Efe C, Purnak T, Ozaslan E, et al. Autoimmune liver disease in patients with systemic lupus erythematosus: a retrospective analysis of 147 cases. Scand J Gastroenterol. 2011;46(6):732–7.

    Article  CAS  PubMed  Google Scholar 

  19. Takahashi A, Abe K, Saito R, et al. Liver dysfunction in patients with systemic lupus erythematosus. Intern Med. 2013;52(13):1461–5.

    Article  PubMed  Google Scholar 

  20. Floreani A, Franceschet I, Cazzagon N, et al. Extrahepatic autoimmune conditions associated with primary biliary cirrhosis. Clin Rev Allergy Immunol. 2015;48(2–3):192–7.

    Article  CAS  PubMed  Google Scholar 

  21. Cheng C, Wang Z, Wang L, et al. Clinical characteristics and prognosis of concomitant systemic lupus erythematosus and primary biliary cholangitis. Clin Rheumatol. 2021;40(5):1819–26.

    Article  PubMed  Google Scholar 

  22. Liu Y, Han K, Liu C, et al. Clinical characteristics and prognosis of concomitant primary biliary cholangitis and autoimmune diseases: a retrospective study. Can J Gastroenterol Hepatol. 2021;2021:5557814.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Efe C, Torgutalp M, Henriksson I, et al. Extrahepatic autoimmune diseases in primary biliary cholangitis: Prevalence and significance for clinical presentation and disease outcome. J Gastroenterol Hepatol. 2021;36(4):936–42.

    Article  CAS  PubMed  Google Scholar 

  24. Liang J, Hu Z, Zhan C, et al. Using propensity score matching to balance the baseline characteristics. J Thorac Oncol. 2021;16(6):e45–6.

    Article  PubMed  Google Scholar 

  25. Qian B, Hu L, Zhang S, et al. Comparison of clinicopathological features and prognosis of papillary thyroid carcinoma and microcarcinoma: A population-based propensity score matching analysis. Front Endocrinol (Lausanne). 2022;13:944758.

    Article  PubMed  Google Scholar 

  26. Xu J, Zhang S, Tian Y, et al. Genetic causal association between iron status and osteoarthritis: a two-sample mendelian randomization. Nutrients. 2022;14(18):3683.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Bastarache L, Denny JC, Roden DM. Phenome-wide association studies. JAMA. 2022;327(1):75–6.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Chen W, Yang Q, Hu L, et al. Shared diagnostic genes and potential mechanism between PCOS and recurrent implantation failure revealed by integrated transcriptomic analysis and machine learning. Front Immunol. 2023;14:1175384.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Zhu E, Shu X, Xu Z, et al. Screening of immune-related secretory proteins linking chronic kidney disease with calcific aortic valve disease based on comprehensive bioinformatics analysis and machine learning. J Transl Med. 2023;21(1):359.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Hermann J, Gruber S, Neufeld JB, et al. IL10R1 loss-of-function alleles in rheumatoid arthritis and systemic lupus erythematosus. Clin Exp Rheumatol. 2009;27(4):603–8.

    CAS  PubMed  Google Scholar 

  31. Stärkel P, Bishop K, Horsmans Y, et al. Expression and DNA-binding activity of signal transducer and activator of transcription 3 in alcoholic cirrhosis compared to normal liver and primary biliary cirrhosis in humans. Am J Pathol. 2003;162(2):587–96.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Li W, Wang R, Wang W. Exploring the causality and pathogenesis of systemic lupus erythematosus in breast cancer based on Mendelian randomization and transcriptome data analyses. Front Immunol. 2023;13:1029884.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Li P, Wang H, Guo L, et al. Association between gut microbiota and preeclampsia-eclampsia: a two-sample Mendelian randomization study. BMC Med. 2022;20(1):443.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Our thanks go out to the original GWAS participants and investigators for sharing and managing the summary statistics.

Funding

This study received no funding.

Author information

Author notes

  1. Linyong Wu was the first author.

Authors and Affiliations

  1. Department of Medical Ultrasound, Maoming People’s Hospital, Maoming, Guangdong, 525000, People’s Republic of China

    Linyong Wu, Songhua Li, Chaojun Wu, Shaofeng Wu, Yan Lin & Dayou Wei

Authors
  1. Linyong Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Songhua Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Chaojun Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Shaofeng Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yan Lin
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Dayou Wei
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

LYW and DYW conceived and designed the study. LYW, SHL, CJW, SFW and YL collected the data. LYW and DYW analyzed the data and results. LYW and DYW completed the writing of the manuscript. LYW, SHL, CJW, SFW, YL and DYW reviewed and edited the manuscript. All authors contributed to the article and approved the submitted version.

Corresponding authors

Correspondence to Songhua Li or Dayou Wei.

Ethics declarations

Ethics approval and consent to participate

Data were collected from public databases. And there weas no ethical approval necessary.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, L., Li, S., Wu, C. et al. Causal relationship between systemic lupus erythematosus and primary liver cirrhosis based on two-sample bidirectional Mendelian randomization and transcriptome overlap analysis. Arthritis Res Ther 26, 10 (2024). https://doi.org/10.1186/s13075-023-03235-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13075-023-03235-z

Keywords

Arthritis Research & Therapy

ISSN: 1478-6362

Contact us