Li Y, Quan X, Tai Y, Wu YT, Wei B, Wu H. Causal association between 731 immunocyte phenotypes and liver cirrhosis: A bidirectional two-sample mendelian randomization analysis. World J Hepatol 2024; 16(8): 1156-1166 [DOI: 10.4254/wjh.v16.i8.1156]
Corresponding Author of This Article
Hao Wu, PhD, Doctor, Professor, Department of Gastroenterology, West China Hospital, Sichuan University, No. 37 Guo Xue Xiang, Chengdu 610041, Sichuan Province, China. hxxhwh@163.com
Research Domain of This Article
Gastroenterology & Hepatology
Article-Type of This Article
Clinical and Translational Research
Open-Access Policy of This Article
This article is an open-access article which was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Ying Li, Xin Quan, Yang Tai, Bo Wei, Hao Wu, Department of Gastroenterology, West China Hospital, Sichuan University, Chengdu 610041, Sichuan Province, China
Yu-Tong Wu, Department of Clinical Medicine, Chongqing Medical University, Chongqing 400016, China
Author contributions: Li Y and Quan X designed the study, reviewed the literature, and wrote the manuscript; Tai Y and Wu YT performed the data analyses; Wei B downloaded the data from the genome-wide association studies database, and Wu H participated in the drafting and editing of the manuscript; All the authors have read and approved the final manuscript for submission.
Supported bythe National Natural Science Foundation of China, No. 82270649.
Conflict-of-interest statement: The authors declare that they have no conflict of interest.
Open-Access: This article is an open-access article that was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution NonCommercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: https://creativecommons.org/Licenses/by-nc/4.0/
Corresponding author: Hao Wu, PhD, Doctor, Professor, Department of Gastroenterology, West China Hospital, Sichuan University, No. 37 Guo Xue Xiang, Chengdu 610041, Sichuan Province, China. hxxhwh@163.com
Received: June 2, 2024 Revised: July 24, 2024 Accepted: August 2, 2024 Published online: August 27, 2024 Processing time: 81 Days and 15.8 Hours
Abstract
BACKGROUND
Liver cirrhosis is a progressive hepatic disease whose immunological basis has attracted increasing attention. However, it remains unclear whether a concrete causal association exists between immunocyte phenotypes and liver cirrhosis.
AIM
To explore the concrete causal relationships between immunocyte phenotypes and liver cirrhosis through a mendelian randomization (MR) study.
METHODS
Data on 731 immunocyte phenotypes were obtained from genome-wide association studies. Liver cirrhosis data were derived from the Finn Gen dataset, which included 214403 individuals of European ancestry. We used inverse variable weighting as the primary analysis method to assess the causal relationship. Sensitivity analyses were conducted to evaluate heterogeneity and horizontal pleiotropy.
RESULTS
The MR analysis demonstrated that 11 immune cell phenotypes have a positive association with liver cirrhosis [P < 0.05, odds ratio (OR) > 1] and that 9 immunocyte phenotypes were negatively correlated with liver cirrhosis (P < 0.05, OR < 1). Liver cirrhosis was positively linked to 9 immune cell phenotypes (P < 0.05, OR > 1) and negatively linked to 10 immune cell phenotypes (P < 0.05; OR < 1). None of these associations showed heterogeneity or horizontally pleiotropy (P > 0.05).
CONCLUSION
This bidirectional two-sample MR study demonstrated a concrete causal association between immunocyte phenotypes and liver cirrhosis. These findings offer new directions for the treatment of liver cirrhosis.
Core Tip: The causal relationship between immunocyte phenotypes and liver cirrhosis has not been fully elucidated. This bidirectional two-sample mendelian randomization study identified a significant causal association between 731 immunocyte phenotypes and liver cirrhosis. We found that 20 immunocyte phenotypes were associated with liver cirrhosis (P < 0.05), whereas liver cirrhosis was associated with 19 immunocyte phenotypes (P < 0.05). None of these associations showed heterogeneity or horizontal pleiotropy (P > 0.05). These findings provide novel directions for the treatment of liver cirrhosis.
Citation: Li Y, Quan X, Tai Y, Wu YT, Wei B, Wu H. Causal association between 731 immunocyte phenotypes and liver cirrhosis: A bidirectional two-sample mendelian randomization analysis. World J Hepatol 2024; 16(8): 1156-1166
Liver cirrhosis represents a chronic and advancing stage of liver disease, characterized by widespread fibrosis and the appearance of pseudolobules as a result of the distortion of normal liver structure by excessive connective tissue. This disease has various aetiologies, and it ultimately leads to multiorgan/system dysfunction[1,2]. It is estimated to be the 11th highest cause of death and the 15th highest contributor to disability-associated life years[3]. The progression from the asymptomatic compensated phase to the symptomatic decompensated phase which is characterized by liver function impairment and portal hypertension results in hospitalization, decreased quality of life, and increased mortality[4,5]. There are currently no treatments that normalize the number and function of liver cells; thus, suppressing the aetiological factor(s) that cause liver inflammation and cirrhosis development is the main method of managing decompensated liver cirrhosis[4]. However, in addition to treatments involving aetiologic factor(s), treatments that can prevent or delay disease progression, complications and multiorgan dysfunction based on the pathophysiologic mechanisms should be identified and implemented[6].
The liver possesses a specialized vascular structure and creates a distinctive immune environment for immune cells[7]. In addition to the sinusoids being composed of macrophages, the liver harbors a diverse population of immune cells, comprising dendritic cells (DCs), innate lymphoid cells, B lymphocytes, and T lymphocytes[8]. Liver inflammation and immune microenvironment changes have been identified as crucial factors in the pathogenesis of cirrhosis. Liver sinusoidal endothelial cells perform a vital function by capturing soluble antigens in a living organism and presenting them to cluster of differentiation (CD) 8 + T cells, resulting in the establishment of tolerance specific to these antigens[9,10]. Hepatic stellate cells (HSCs) have been suggested to present antigens and lead to the activation of CD1d-restricted natural killer T cells (NKTs) and CD4 + and CD8 + T cells[11]. Increasing research has demonstrated that immune cells are pivotal in regulating both the advancement and resolution of liver fibrosis[12,13]. Ramachandran et al[14] profiled the transcriptomes of over 100000 single human cells from both healthy and diseased livers, revealing an expansion of a scar-associated macrophage subpopulation expressing TREM2 + CD9 + during liver fibrosis. Nakamoto et al[15] reported that durable CD8 + T-cell-dependent liver disease of moderate severity plays an essential role in the development of liver cirrhosis and hepatocellular carcinoma in a mouse model. Furthermore, Rueschenbaum et al[16] demonstrated that adaptive cellular immunity plays a role at a relatively early stage in the development of liver cirrhosis and the capacity of CD4 + and CD8 + T cells to produce proinflammatory cytokines is lower in a prospective cohort study of patients with different stages of liver cirrhosis or acute-on-chronic liver failure than in healthy subjects. All of the abovementioned findings indicated that there is a multifaceted relationship between immune cells and liver cirrhosis. A deeper insight into these relationships might foster the innovation of improved therapeutic modalities and assist in pinpointing patients most likely to respond favorably to immunotherapy. Mendelian randomization (MR) is an important analytical method for inferring causal relationships in epidemiology based on Mendelian genetic principles, which are screened for instrumental variables (IVs) strongly associated with exposure factors by strict criteria to evaluate the relationships between exposure factors and outcomes[17,18]. Previous observational studies have suggested that there are complex interactions between immune cell traits and liver cirrhosis, thus providing evidence for their putative association[16]. This study aimed to elucidate the causal link between 731 distinct immunocyte phenotypes and liver cirrhosis through a comprehensive application of bidirectional two-sample MR analysis.
MATERIALS AND METHODS
Study design
We examined the causal associations between 731 immune cell characteristics and liver cirrhosis via a bidirectional two-sample MR analysis. All genetic variations used as IVs must satisfy three fundamental assumptions: (1) IVs are directly associated with the exposure; (2) IVs are unrelated to confounders of the exposure-outcome connection; and (3) IVs do not affect outcome through pathways other than exposure[19]. Figure 1 provides an overview of the study design. The data used in this study were obtained from the Finngen and OPEN genome-wide association studies (GWAS) public databases in this study. This database is publicly accessible, thus eliminating the need for additional ethical approval.
Figure 1 Diagram of this bidirectional mendelian randomization study design.
The arrows indicate bidirectional causal relationships between liver cirrhosis and immunocyte phenotypes, and the causal pathway is blocked if a “fork” is placed in the arrowed line. IVW: Inverse variance weighting; MR: Mendelian randomization; SNPs: Single-nucleotide polymorphisms; IVs: Instrumental variables.
Data sources
Source of immune cell data: The GWAS catalogue (GCST90001391 to GCST90002121) provides summary statistics for all immunologic characteristics[20]. GWASs of immune-related traits involve extensive data derived from 3757 non-overlapping European participants. Through the use of a high-density array generated from Sardinian sequence data, we ascertained approximately 22 million single-nucleotide polymorphisms (SNPs) and examined their associations, taking into account covariates like sex, age, and age squared. In total, 731 immunophenotypic parameters were investigated, comprising 118 absolute cell counts, 389 median fluorescence intensities indicative of surface antigen expression, 32 morphological attributes, and 192 relative cell counts[21].
Source of liver cirrhosis data: We acquired a compendium of GWAS summary statistics pertaining to liver cirrhosis from FinnGen database (https://www.finngen.fi/en). This dataset encompassed a total of 214403 samples (case number = 213592, control number = 811), and the GWAS integrated over 21306350 phenotypic datapoints linked to liver cirrhosis, thereby uncovering in excess of 23 million independent SNPs.
Selection of IVs
These standards must be fulfilled by the candidate IVs. First, SNPs with a P value < 1e-05 were needed. Second, we excluded any SNPs with significant linkage disequilibrium (r² < 0.001 and kb < 10000) to ensure that the analysis results would not be compromised by SNP interdependencies. Moreover, the PhenoScanner online tool was used to exclude SNPs related to immune cell and liver cirrhosis confounders. In addition, to avoid bias from weak IVs and evaluate the statistical power of the relationship between each SNP and the exposures, we utilized the F statistic [calculated as F = R2 (N-2)/(1-R2), with R² being the variance in exposure explained by the genetic variant, approximated as 2 × EAF × (1 - EAF) × β; here, EAF is the effect allele’s frequency, β is the genetic effect estimate on exposure, and N is the sample size from the exposure GWAS][22,23]. IVs with an F statistic > 10 were considered strong instruments. Finally, harmonization of the SNP exposures and outcomes was conducted to ensure that the effect alleles for each SNP on the exposure matched their corresponding effect alleles on the outcome. We excluded SNPs that had mismatched alleles from our analysis. Additionally, palindromic SNPs and those with ambiguous genotyping were identified and excluded during the harmonization process.
Statistical analysis
The MR analyses in this study were performed using R 4.4.0 software (https://www.r-project.org/) using the two sample MR package. First, significant SNPs were identified according to the above criteria. Moreover, the exposure and outcome data were harmonized, and the two-sample MR effect was calculated. Five distinct MR analysis techniques were employed: Inverse variance weighted (IVW), weighted median, simple mode, weighted mode, and MR-Egger regression. The IVW serves as the primary analytical approach, which involves computing a weighted average of the ratio estimates for each genetic variant. Next, sensitivity analyses were conducted to test potential pleiotropy, such as heterogeneity and horizontal multiple validity. We employed Cochran’s Q test to measure the heterogeneity among the IVs and utilized MR-Egger regression for weighted linear regression with intercepts to detect potential horizontal pleiotropy within the IVs. Heterogeneity and pleiotropy were considered not present if the P value was larger than 0.05. In addition, A leave-one-out analysis was conducted by iteratively excluding each genetic variant from the analysis and recalculating the causal effect to determine if any single SNP exerted a substantial influence on the overall causal effect. The outcomes are expressed as odds ratios (OR) accompanied by 95% confidence intervals. Statistical significance was defined as a P value less than 0.05.
RESULTS
Causal effects of immunocytes on liver cirrhosis
From the GWAS data encompassing 731 immunocyte phenotypes, we identified candidate IVs. Our analysis confirmed that the F-statistic for these IVs was above 10, ruling out the likelihood of weak instrument bias. All SNPs with positive results are shown in Table 1. IVW analysis revealed 20 immune cell phenotypes with potential causal relationships with liver cirrhosis. The following 11 phenotypes of immune cells were positively associated with the development of liver cirrhosis (OR > 1, P < 0.05). Treg panel: CD28 + CD45 receptor alpha (RA) + CD8 bright% T cell; B-Cell panel: CD25 on memory B cell, CD25 on sw mem, CD38 on CD20-, and B-cell-activating factor (BAFF)-R on CD20-; Myeloid cell panel: CD33 on CD66b ++ myeloid cell, CD33 on Im myeloid-derived suppressor cells (MDSC), CD33 on CD33dim human leukocyte antigen (HLA) dweller region (DR)-, CD33 on basophil, CD33 bright HLA DR + CD14dim %CD33 bright HLA DR +, CD33 on monocytic (Mo) MDSC. The remaining 9 phenotypes were negatively associated with the progression of liver cirrhosis (OR < 1, P < 0.05). Conventional DC (cDC) panel: C-C motif chemokine receptor 2 (CCR2) on CD62 L + myeloid DC; Treg panel: CD39 on CD39 + secreting Treg, CD28 on CD39 + resting Treg; B panel: BAFF-R on IgD + CD38- naive, IgD- CD27- % lymphocyte, CD20 on IgD +; Myeloid cell: HLA DR on CD33dim HLA DR + CD11b +; T cells, B cells, and natural killer cells (TBNK) panel: CD8 bright NKT absolute count (AC), side scatter (SSC)-A on NK. The results of IVW analysis are provided in Figure 2, and the weighted median, simple mode, weighted mode, and MR-Egger regression methods used for MR analysis are depicted in Supplementary Figure 1.
Figure 2 Forest plots showed the causal associations between liver cirrhosis and immune cell traits.
CI: Confidence interval; OR: Odds ratio; TBNK: T cells, B cells, and natural killer cells; cDC: Conventional dendritic cells; CD: Cluster of differentiation; BAFF: B-cell-activating factor; HLA: Human leukocyte antigen; DR: Dweller region; AC: Absolute count; MDSC: Myeloid-derived suppressor cells; NKT: Natural killer T cells; NK: Natural killer; CCR2: C-C motif chemokine receptor 2.
Table 1 Number of single-nucleotide polymorphisms screened in each step.
ID exposure
Panel
Immune traits
Number of SNPs after LD
Number of SNPs after F > 10
Number of final IVs
Ebi-a-GCST90001433
B cell
IgD- CD27-% lymphocyte
18
18
18
Ebi-a-GCST90001706
B cell
BAFF-R on IgD + CD38- naive
29
29
28
Ebi-a-GCST90001762
B cell
CD20 on IgD +
23
23
19
Ebi-a-GCST90001790
B cell
CD25 on memory B cell
27
27
24
Ebi-a-GCST90001830
B cell
BAFF-R on CD20-
14
14
14
Ebi-a-GCST90001793
B cell
CD25 on sw mem
21
21
20
Ebi-a-GCST90001809
B cell
CD38 on CD20-
20
20
19
Ebi-a-GCST90002110
Myeloid cell
HLA DR on CD33dim HLA DR + CD11b +
24
24
23
Ebi-a-GCST90001954
Myeloid cell
CD33 on basophil
23
23
22
Ebi-a-GCST90001952
Myeloid cell
CD33 on Mo MDSC
21
21
21
Ebi-a-GCST90001951
Myeloid cell
CD33 on CD66b + + myeloid cell
19
19
18
Ebi-a-GCST90001521
Myeloid cell
CD33 bright HLA DR + CD14 dim% CD33 bright HLA DR +
26
26
25
Ebi-a-GCST90001953
Myeloid cell
CD33 on CD33 dim HLA DR-
21
21
19
Ebi-a-GCST90001955
Myeloid cell
CD33 on Im MDSC
25
25
22
Ebi-a-GCST90002031
Treg
CD39 on CD39 + secreting Treg
24
24
23
Ebi-a-GCST90001901
Treg
CD28 on CD39 + resting Treg
23
23
22
Ebi-a-GCST90001688
Treg
CD28 + CD45 RA + CD8 bright% T cell
112
91
88
Ebi-a-GCST90001630
TBNK
CD8 bright NKT AC
27
27
26
Ebi-a-GCST90002076
TBNK
SSC-A on NK
24
24
23
Ebi-a-GCST90002014
cDC
CCR2 on CD62 L + myeloid DC
15
15
14
Forward sensitivity analyses
Table 2 presents the results from our sensitivity analysis, which confirmed that there was no significant heterogeneity present. (P > 0.05 for the Q test) or directional pleiotropy (P > 0.05 for the MR-Egger intercept method) for 20 immune cell immunocytes in the MR analysis of liver cirrhosis patients. Moreover, the leave-one-out method and funnel plots indicated that the data were robust, as shown in Supplementary Figure 2 and Supplementary Figure 3.
The GWAS data for liver cirrhosis patients were screened for IVs, and all selected IVs exhibited F values exceeding 10, indicating a lower possibility of weak instrument bias. The SNPs that showed positive associations are summarized in Table 3.
Table 3 Number of Single-nucleotide polymorphisms screened in each step.
ID exposure
Number of SNPs after LD
Number of SNPs after F > 10
Number of final IVs
Finn-b-CIRRHOSIS_BROAD
31
31
31
The results of IVW for the causal effects of liver cirrhosis on immunocytes are provided in Figure 3, and other analytical methods are provided in Supplementary Figure 4. The immune cell immunocytes with a positive association were as follows (OR > 1, P < 0.05): B-cell panel: CD20 on CD24 + CD27 + and CD20 on sw mem; TBNK panel: CD4 + %leukocyte, forward scatter (FSC)-A on B cell, and SSC-A on B cell; Myeloid cell panel: Im MDSC% CD33 dim HLA DR- CD66b-, Im MDSC AC; Monocyte panel: CX3CR1 on monocyte, CX3CR1 on CD14- CD16 + monocyte. The remaining 10 immunocytes of immune cells with a negative correlation were as follows (OR < 1, P < 0.05): B Cell Panel: IgD + CD38 bright AC, IgD + CD38 dim AC, IgD + AC, naive-mature B cell AC, IgD + CD24- AC; cDC panel: CD62 L-plasmacytoid DC AC; TBNK panel: CD3- lymphocyte% leukocyte; Myeloid cell: CD33- HLA DR + AC; Treg panel: CD8 on CD28- CD8 bright, CD28 + CD45RA- CD8 dim AC.
Figure 3 Forest plots showed the causal associations between immune cell traits and liver cirrhosis.
CI: Confidence interval; OR: Odds ratio; TBNK: T cells, B cells, and natural killer cells; cDC: Conventional dendritic cells; CD: Cluster of differentiation; AC: Absolute count; DR: Dweller region; FSC-A: Forward scatter-A; SSC-A: Side scatter-A; MDSC: Myeloid-derived suppressor cells.
Reverse sensitivity analysis
The results of the sensitivity analysis (Table 4) showed that there was no heterogeneity (P > 0.05 for the Q test) or directional pleiotropy (P > 0.05 for MR-Egger’s intercept method) in the 19 immunocytes from the MR analysis of liver cirrhosis patients. Moreover, the leave-one-out method and funnel plots indicated that the results were robust, as shown in Supplementary Figure 5 and Supplementary Figure 6.
We first performed a bidirectional MR investigation, which examined the causal associations between 731 immune cell phenotypes and liver cirrhosis based on a large amount of publicly available genetic data. We found that 20 immunophenotypes had significant causal effects on liver cirrhosis, and liver cirrhosis had a causal effect on 19 immunophenotypes. Moreover, no single immune cell phenotype was observed to have a bidirectional causal association.
B cells play diverse roles in liver cirrhosis by mediating antigen presentation and cytokine release, which trigger the activation of various immune cell phenotypes[24]. In our study, CD25 on memory B cell, CD25 on sw mem cell, CD38 on CD20- cell, and BAFF-R on CD20- cell in the B-cell panel were shown to be associated with increased liver cirrhosis risk. Moreover, we also found that the risk of liver cirrhosis decreased with increasing proportions of BAFF-R on IgD + CD38- naive cell, IgD- CD27-% lymphocyte cell, and CD20 on IgD + cell. BAFF has been identified as a crucial regulator of peripheral B-cell survival, homeostasis, and antibody-mediated responses[25,26]. Yang et al[27] reported that BAFF levels are elevated in individuals with chronic hepatitis B virus (HBV) infection and that BAFF upregulation in patients with liver cirrhosis is more prominent. The serum levels of BAFF are also elevated in patients with non-alcoholic steatohepatitis (NASH) and are correlated with the degree of fibrosis[28]. Furthermore, there is evidence that CD20 + B cells are also related to the occurrence and progression of cirrhosis and that the depletion of CD20 + B cells plays a crucial role in improving cirrhosis[29,30].
Tregs constitute a population of CD4 + T helper cells that mainly express the transcription factors forkhead box P3 and signal transducer and activator of transcription 5. The number of Tregs is increased in patients with liver disorders[12]. In the liver, CXCR6 + CD8 T cells were found to exhibit reduced activity of the FOXO1 transcription factor. However, they are abundant in patients with NASH and are related to self-directed immune responses[31]. Our study revealed that CD28 + CD45RA + CD8 bright% T cell increased the risk of liver cirrhosis, and CD39 on CD39 + secreting Treg and CD28 on CD39 + resting Treg decreased the risk of liver cirrhosis. A study based on a bile duct ligation rat model suggested that Tregs protect the liver from cholestasis and fibrosis[32]. Additional research indicates that Tregs are essential for protecting the liver from damage in the chronic phase of HBV infection, aiding in the prevention of the disease’s progression to cirrhosis and hepatocellular carcinoma[33]. However, Langhans et al[34] demonstrated that Tregs activated HSCs and promoted liver fibrogenesis by producing interleukin (IL)-8 in chronic hepatitis C virus (HCV). We postulate that the various roles of Tregs in the pathogenesis of liver cirrhosis might be due to alterations in protein synthesis on their cell surface in response to distinct aetiological factors and immune microenvironments, resulting in the secretion of a broad spectrum of cytokines.
MDSCs encompass a heterogeneous mixture of myeloid progenitor cells and immature myeloid cells, and can maintain their suppressive function through the induction of Tregs[35]. MDSCs from inflamed portal tracts and patients with advanced inflammation and hepatic fibrosis are significantly elevated[36]. We found that six MDSC phenotypes, including CD33 on CD66b + myeloid cells, CD33 on Mo MDSCs, CD33 on basophils, CD33 bright HLA DR + CD14 dim% CD33 bright HLA DR +, CD33 on CD33 dim HLA DR-, and CD33 on Im MDSCs have positive causal associations with liver cirrhosis, and HLA DR on CD33dim HLA DR + CD11b + has a negative association with liver cirrhosis. One study revealed an increase in monocytic MDSCs with a phenotype of HLA-DR-/Low CD33 + CD11b + CD14 + CD15- in patients with liver inflammatory diseases compared to healthy controls[36]. In addition, previous studies have shown that the accumulation of CD33 + MDSCs is detected in patients with HCV infection and results in reactive oxygen species-mediated suppression of T-cell responsiveness, which may lead to liver cirrhosis[37]. These observations align with our MR outcomes.
NK cells demonstrate antifibrotic effects by targeting and eliminating HSCs, acting as a critical first line of defense against viral hepatitis. They augment the antiviral immune response by directly killing virus-infected cells or by promoting antigen-specific T cell responses through the release of interferon-γ and tumor necrosis factor-α[38,39]. Our MR analysis findings also suggested that CD8 bright NKT AC and SSC-A on NK cells play a role in improving liver cirrhosis.
DCs, a subset of mononuclear phagocytes that display major histocompatibility complex class II molecules, can be classified into myeloid or classical DCs and plasmacytoid DCs in the blood[40,41]. Our study revealed that CCR2 on CD62 L + myeloid DC was negatively associated with liver cirrhosis. CD62 L, an adhesive molecule from the selection family, facilitates interactions with endothelial cells and migration into tissues. Compared with controls, cirrhotic patients presented a marked decrease in cDCs and greater percentages of monocytes expressing CD62 L[40,42].
Additionally, we found that the causally associated immunocyte phenotypes were not inversely causally associated. Notably, the presence of liver cirrhosis was associated with increased CD20 on CD24+, CD27+; CD20 on sw mem; CD4 +% leukocyte; FSC-A on B cell; SSC-A on B cell; Im MDSC %CD33dim HLA DR- CD66b-; Im MDSC AC; CX3CR1 on monocyte; and CX3CR1 on CD14- CD16 + monocyte levels. Doi et al[43] reported that CD27 + memory B cells were markedly less common in cirrhotic patients and that CD70 upregulation, tumour necrosis factor beta secretion, and IgG production were observed. The proportion of total B cells, characterized by a mature phenotype, was found to be elevated, whereas the percentage of memory B cells was significantly reduced in individuals with decompensated cirrhosis compared with healthy individuals. Furthermore, the serum concentrations of IL-10, IL-21, and IL-4 are substantially decreased in patients with decompensated cirrhosis[44].
The role of various immune cell phenotypes in the pathogenesis of liver cirrhosis is complex and multifaceted. Depending on the underlying aetiology and stage of liver cirrhosis, distinct immune cell phenotypes are generated under the guidance of related genes, leading to the secretion of either pro- or anti-liver fibrosis cytokines, which in turn perform diverse functions in the development of liver cirrhosis. Our study offers novel insights into the immune cells participating in the immune response culminating in liver cirrhosis, potentially paving the way for enhanced treatment strategies and aiding in the identification of patients who would derive the greatest benefit from immunotherapy.
We conducted bidirectional two-sample MR analyses utilizing extensive GWAS datasets with large sample sizes, which significantly enhanced the statistical efficiency of our study. By employing genetic variants as IVs and applying various MR techniques, we were able to make robust causal inferences. These methods mitigate the potential impact of horizontal pleiotropy and confounding factors, reinforcing the reliability of our results. There are several limitations in our study. Firstly, we used a screening threshold of P < 1 × 10-5 to identify IVs, which was not strong enough, although a comprehensive assessment of the association between the immune cell phenotype and liver cirrhosis was conducted. Second, the study was based on a European database and genetically differentiated from other populations, so the conclusions cannot be generalized to other ethnic groups. Additionally, owing to the unavailability of individual data, we were unable to perform further stratified analysis of the population. Moreover, confounders could not be eliminated, although we conducted a sensitivity analysis to exclude SNPs associated with potential confounders as much as possible. Finally, to make clinical inferences, these findings need to be further validated in extensive clinical trials.
CONCLUSION
In summary, our research identified causal links between immunocyte phenotypes and liver cirrhosis through comprehensive bidirectional MR analysis, thus highlighting the intricate interplay between the immune system and liver cirrhosis and providing insight into early interventions and therapeutic strategies for liver cirrhosis.
ACKNOWLEDGEMENTS
We thank FinnGen and the database for sharing the data.
Footnotes
Provenance and peer review: Unsolicited article; Externally peer reviewed.
Peer-review model: Single blind
Specialty type: Gastroenterology and hepatology
Country of origin: China
Peer-review report’s classification
Scientific Quality: Grade B
Novelty: Grade A
Creativity or Innovation: Grade A
Scientific Significance: Grade A
P-Reviewer: Shah SZA S-Editor: Fan M L-Editor: A P-Editor: Chen YX
Kawashima K, Andreata F, Beccaria CG, Iannacone M. Priming and Maintenance of Adaptive Immunity in the Liver.Annu Rev Immunol. 2024;42:375-399.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 2][Reference Citation Analysis (0)]
Limmer A, Ohl J, Kurts C, Ljunggren HG, Reiss Y, Groettrup M, Momburg F, Arnold B, Knolle PA. Efficient presentation of exogenous antigen by liver endothelial cells to CD8+ T cells results in antigen-specific T-cell tolerance.Nat Med. 2000;6:1348-1354.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 540][Cited by in F6Publishing: 534][Article Influence: 22.3][Reference Citation Analysis (0)]
Rueschenbaum S, Ciesek S, Queck A, Widera M, Schwarzkopf K, Brüne B, Welsch C, Wedemeyer H, Zeuzem S, Weigert A, Lange CM. Dysregulated Adaptive Immunity Is an Early Event in Liver Cirrhosis Preceding Acute-on-Chronic Liver Failure.Front Immunol. 2020;11:534731.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 8][Cited by in F6Publishing: 20][Article Influence: 6.7][Reference Citation Analysis (0)]
Burgess S, Mason AM, Grant AJ, Slob EAW, Gkatzionis A, Zuber V, Patel A, Tian H, Liu C, Haynes WG, Hovingh GK, Knudsen LB, Whittaker JC, Gill D. Using genetic association data to guide drug discovery and development: Review of methods and applications.Am J Hum Genet. 2023;110:195-214.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 37][Cited by in F6Publishing: 35][Article Influence: 35.0][Reference Citation Analysis (0)]
Orrù V, Steri M, Sidore C, Marongiu M, Serra V, Olla S, Sole G, Lai S, Dei M, Mulas A, Virdis F, Piras MG, Lobina M, Marongiu M, Pitzalis M, Deidda F, Loizedda A, Onano S, Zoledziewska M, Sawcer S, Devoto M, Gorospe M, Abecasis GR, Floris M, Pala M, Schlessinger D, Fiorillo E, Cucca F. Complex genetic signatures in immune cells underlie autoimmunity and inform therapy.Nat Genet. 2020;52:1036-1045.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 28][Cited by in F6Publishing: 237][Article Influence: 59.3][Reference Citation Analysis (0)]
Wang YX, Zhou CP, Wang DT, Ma J, Sun XH, Wang Y, Zhang YM. Unraveling the causal role of immune cells in gastrointestinal tract cancers: insights from a Mendelian randomization study.Front Immunol. 2024;15:1343512.
[PubMed] [DOI][Cited in This Article: ][Reference Citation Analysis (0)]
Miyake T, Abe M, Tokumoto Y, Hirooka M, Furukawa S, Kumagi T, Hamada M, Kawasaki K, Tada F, Ueda T, Hiasa Y, Matsuura B, Onji M. B cell-activating factor is associated with the histological severity of nonalcoholic fatty liver disease.Hepatol Int. 2013;7:539-547.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 30][Cited by in F6Publishing: 39][Article Influence: 3.3][Reference Citation Analysis (0)]
Faggioli F, Palagano E, Di Tommaso L, Donadon M, Marrella V, Recordati C, Mantero S, Villa A, Vezzoni P, Cassani B. B lymphocytes limit senescence-driven fibrosis resolution and favor hepatocarcinogenesis in mouse liver injury.Hepatology. 2018;67:1970-1985.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 52][Cited by in F6Publishing: 56][Article Influence: 9.3][Reference Citation Analysis (0)]
Petrarca A, Rigacci L, Caini P, Colagrande S, Romagnoli P, Vizzutti F, Arena U, Giannini C, Monti M, Montalto P, Matucci-Cerinic M, Bosi A, Laffi G, Zignego AL. Safety and efficacy of rituximab in patients with hepatitis C virus-related mixed cryoglobulinemia and severe liver disease.Blood. 2010;116:335-342.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 97][Cited by in F6Publishing: 99][Article Influence: 7.1][Reference Citation Analysis (0)]
Katz SC, Ryan K, Ahmed N, Plitas G, Chaudhry UI, Kingham TP, Naheed S, Nguyen C, Somasundar P, Espat NJ, Junghans RP, Dematteo RP. Obstructive jaundice expands intrahepatic regulatory T cells, which impair liver T lymphocyte function but modulate liver cholestasis and fibrosis.J Immunol. 2011;187:1150-1156.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 46][Cited by in F6Publishing: 52][Article Influence: 4.0][Reference Citation Analysis (0)]
Langhans B, Krämer B, Louis M, Nischalke HD, Hüneburg R, Staratschek-Jox A, Odenthal M, Manekeller S, Schepke M, Kalff J, Fischer HP, Schultze JL, Spengler U. Intrahepatic IL-8 producing Foxp3⁺CD4⁺ regulatory T cells and fibrogenesis in chronic hepatitis C.J Hepatol. 2013;59:229-235.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 58][Cited by in F6Publishing: 56][Article Influence: 5.1][Reference Citation Analysis (0)]
Zhang H, Lian M, Zhang J, Bian Z, Tang R, Miao Q, Peng Y, Fang J, You Z, Invernizzi P, Wang Q, Gershwin ME, Ma X. A functional characteristic of cysteine-rich protein 61: Modulation of myeloid-derived suppressor cells in liver inflammation.Hepatology. 2018;67:232-246.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 28][Cited by in F6Publishing: 33][Article Influence: 5.5][Reference Citation Analysis (0)]
Cardoso CC, Matiollo C, Pereira CHJ, Fonseca JS, Alves HEL, da Silva OM, de Souza Menegassi V, Dos Santos CR, de Moraes ACR, de Lucca Schiavon L, Santos-Silva MC. Patterns of dendritic cell and monocyte subsets are associated with disease severity and mortality in liver cirrhosis patients.Sci Rep. 2021;11:5923.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 9][Cited by in F6Publishing: 15][Article Influence: 5.0][Reference Citation Analysis (0)]
Doi H, Iyer TK, Carpenter E, Li H, Chang KM, Vonderheide RH, Kaplan DE. Dysfunctional B-cell activation in cirrhosis resulting from hepatitis C infection associated with disappearance of CD27-positive B-cell population.Hepatology. 2012;55:709-719.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 60][Cited by in F6Publishing: 64][Article Influence: 5.3][Reference Citation Analysis (0)]