
==== Front
World J Hepatol
WJH
World Journal of Hepatology
1948-5182
Baishideng Publishing Group Inc

39221101
jWJH.v16.i8.pg1156
10.4254/wjh.v16.i8.1156
97523
Clinical and Translational Research
Causal association between 731 immunocyte phenotypes and liver cirrhosis: A bidirectional two-sample mendelian randomization analysis
Li Y et al. Immunocyte phenotypes-liver cirrhosis link
Li Ying Department of Gastroenterology, West China Hospital, Sichuan University, Chengdu 610041, Sichuan Province, China

Quan Xin Department of Gastroenterology, West China Hospital, Sichuan University, Chengdu 610041, Sichuan Province, China

Tai Yang Department of Gastroenterology, West China Hospital, Sichuan University, Chengdu 610041, Sichuan Province, China

Wu Yu-Tong Department of Clinical Medicine, Chongqing Medical University, Chongqing 400016, China

Wei Bo Department of Gastroenterology, West China Hospital, Sichuan University, Chengdu 610041, Sichuan Province, China

Wu Hao Department of Gastroenterology, West China Hospital, Sichuan University, Chengdu 610041, Sichuan Province, China. hxxhwh@163.com

Co-first authors: Ying Li and Xin Quan.

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 by the National Natural Science Foundation of China, No. 82270649 .

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

27 8 2024
27 8 2024
16 8 11561166
2 6 2024
24 7 2024
2 8 2024
©The Author(s) 2024. Published by Baishideng Publishing Group Inc. All rights reserved.
2024
https://creativecommons.org/licenses/by-nc/4.0/ 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.
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.

Liver cirrhosis
Immune cell
Immunocyte phenotype
Mendelian analysis
Causal association
==== Body
pmc 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.

INTRODUCTION

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 CD62L++ 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 CD62L++ myeloid DC	15	15	14	
TBNK: T cells, B cells, and natural killer cells; ID: Identification; cDC: Conventional dendritic cells; CD: Cluster of differentiation; DR: Dweller region; AC: Absolute count; SSC-A: Side scatter-A; LD: Linkage disequilibrium; BAFF: B-cell-activating factor; HLA: Human leukocyte antigen; MDSC: Myeloid-derived suppressor cells; NKT: Natural killer T cells; NK: Natural killer; CCR2: C-C motif chemokine receptor 2; SNPs: Single-nucleotide polymorphisms; IVs: Instrumental variables.

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.

Table 2 Forward mendelian randomization sensitivity analysis

Exposure panel	Exposure traits	Inverse variance weighted	MR-Egger	
Q value	P value	Intercept	P value	
B cell	IgD- CD27-% lymphocyte	11.20	0.846	0.008	0.707	
B cell	BAFF-R on IgD+ CD38- naive	19.98	0.831	0.002	0.880	
B cell	CD20 on IgD+	18.76	0.407	-0.004	0.662	
B cell	CD25 on memory B cell	21.52	0.549	-0.002	0.854	
B cell	BAFF-R on CD20-	9.74	0.715	-0.005	0.795	
B cell	CD25 on sw mem	17.72	0.541	0.004	0.771	
B cell	CD38 on CD20-	14.73	0.680	0.012	0.521	
Myeloid cell	HLA DR on CD33dim HLA DR+ CD11b+	25.66	0.266	-0.011	0.303	
Myeloid cell	CD33 on basophil	29.56	0.101	0.014	0.248	
Myeloid cell	CD33 on Mo MDSC	23.26	0.276	-0.007	0.490	
Myeloid cell	CD33 on CD66b+ myeloid cell	18.20	0.376	-0.029	0.065	
Myeloid cell	CD33 bright HLA DR+ CD14 dim% CD33 bright HLA DR+	29.63	0.197	0.011	0.294	
Myeloid cell	CD33 on CD33 dim HLA DR-	18.45	0.427	0.018	0.184	
Myeloid cell	CD33 on Im MDSC	30.80	0.077	0.013	0.354	
Treg	CD39 on CD39+ secreting Treg	10.76	0.978	-0.005	0.554	
Treg	CD28 on CD39+ resting Treg	21.59	0.423	0.013	0.220	
Treg	CD28+ CD45 RA+ CD8 bright% T cell	80.46	0.676	0.009	0.252	
TBNK	CD8 bright NKT AC	25.19	0.452	-0.014	0.225	
TBNK	SSC-A on NK	27.35	0.198	0.004	0.771	
cDC	CCR2 on CD62L+ myeloid DC	9.06	0.768	0.011	0.547	
TBNK: T cells, B cells, and natural killer; cDC: Conventional dendritic cells; CD: Cluster of Differentiation; AC: Absolute count; BAFF: B-cell-activating factor; HLA: Human leukocyte antigen; MDSC: Myeloid-derived suppressor cells; NKT: Natural killer T cells; NK: Natural killer; CCR2: C-C motif chemokine receptor 2; SNPs: Single-nucleotide polymorphisms; IVs: Instrumental variables; MR: Mendelian randomization.

Causal effects of liver cirrhosis on immunocytes

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	
ID: Identification; LD: Linkage disequilibrium; SNPs: Single-nucleotide polymorphisms; IVs: Instrumental variables.

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: CD62L-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.

Table 4 Reverse mendelian randomization sensitivity analysis

Outcome panel	Outcome traits	Inverse variance weighted	MR-Egger	
Q value	P value	intercept	P value	
B cell	IgD+ CD38 bright AC	35.51	0.224	-0.014	0.401	
B cell	IgD+ CD38 dim AC	30.73	0.429	-0.030	0.068	
B cell	IgD+ AC	33.15	0.316	-0.028	0.094	
B cell	Naive-mature B cell AC	27.47	0.599	-0.022	0.160	
B cell	IgD+ CD24- AC	29.82	0.475	-0.019	0.250	
cDC	CD62L- plasmacytoid DC AC	29.39	0.497	0.029	0.075	
Myeloid cell	Im MDSC AC	29.23	0.506	0.019	0.249	
Myeloid cell	Im MDSC% CD33 dim HLA DR- CD66b-	24.09	0.768	0.009	0.568	
Myeloid cell	CD33- HLA DR+ AC	30.39	0.446	-0.004	0.790	
TBNK	CD4+% leukocyte	16.12	0.982	0.007	0.670	
TBNK	CD3- lymphocyte% leukocyte	17.10	0.971	0.013	0.443	
Treg	CD28+ CD45RA- CD8 dim AC	37.39	0.166	0.000	0.997	
B cell	CD20 on CD24+ CD27+	22.52	0.834	0.026	0.231	
B cell	CD20 on sw mem	32.59	0.341	0.020	0.381	
TBNK	FSC-A on B cell	25.61	0.695	0.017	0.345	
Monocyte	CX3CR1 on monocyte	39.69	0.111	0.002	0.891	
Monocyte	CX3CR1 on CD14- CD16+ monocyte	19.80	0.922	0.008	0.641	
TBNK	SSC-A on B cell	32.70	0.336	-0.003	0.844	
Treg	CD8 on CD28- CD8 bright	37.59	0.106	-0.003	0.877	
TBNK: T cells, B cells, and natural killer cells; DC: Dendritic cells; AC: Absolute count; FSC-A: Forward scatter-A; SSC-A: Side scatter-A; CD: Cluster of differentiation.

DISCUSSION

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 CD62L++ myeloid DC was negatively associated with liver cirrhosis. CD62L, 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 CD62L[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.

Conflict-of-interest statement: The authors declare that they have no conflict of interest.

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
==== Refs
1 Bernardi M Moreau R Angeli P Schnabl B Arroyo V Mechanisms of decompensation and organ failure in cirrhosis: From peripheral arterial vasodilation to systemic inflammation hypothesis J Hepatol 2015 63 1272 1284 26192220
2 Tsochatzis EA Bosch J Burroughs AK Liver cirrhosis Lancet 2014 383 1749 1761 24480518
3 Ginès P Krag A Abraldes JG Solà E Fabrellas N Kamath PS Liver cirrhosis Lancet 2021 398 1359 1376 34543610
4 European Association for the Study of the Liver EASL Clinical Practice Guidelines for the management of patients with decompensated cirrhosis J Hepatol 2018 69 406 460 29653741
5 Garcia-Pagan JC Francoz C Montagnese S Senzolo M Mookerjee RP Management of the major complications of cirrhosis: Beyond guidelines J Hepatol 2021 75 Suppl 1 S135 S146 34039484
6 Alegre F Pelegrin P Feldstein AE Inflammasomes in Liver Fibrosis Semin Liver Dis 2017 37 119 127 28564720
7 Crispe IN The liver as a lymphoid organ Annu Rev Immunol 2009 27 147 163 19302037
8 Kawashima K Andreata F Beccaria CG Iannacone M Priming and Maintenance of Adaptive Immunity in the Liver Annu Rev Immunol 2024 42 375 399 38360545
9 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 11100119
10 Sørensen KK Simon-Santamaria J McCuskey RS Smedsrød B Liver Sinusoidal Endothelial Cells Compr Physiol 2015 5 1751 1774 26426467
11 Winau F Hegasy G Weiskirchen R Weber S Cassan C Sieling PA Modlin RL Liblau RS Gressner AM Kaufmann SH Ito cells are liver-resident antigen-presenting cells for activating T cell responses Immunity 2007 26 117 129 17239632
12 Pellicoro A Ramachandran P Iredale JP Fallowfield JA Liver fibrosis and repair: immune regulation of wound healing in a solid organ Nat Rev Immunol 2014 14 181 194 24566915
13 Liu Y Dong Y Wu X Wang X Niu J Identification of Immune Microenvironment Changes and the Expression of Immune-Related Genes in Liver Cirrhosis Front Immunol 2022 13 918445 35903097
14 Ramachandran P Dobie R Wilson-Kanamori JR Dora EF Henderson BEP Luu NT Portman JR Matchett KP Brice M Marwick JA Taylor RS Efremova M Vento-Tormo R Carragher NO Kendall TJ Fallowfield JA Harrison EM Mole DJ Wigmore SJ Newsome PN Weston CJ Iredale JP Tacke F Pollard JW Ponting CP Marioni JC Teichmann SA Henderson NC Resolving the fibrotic niche of human liver cirrhosis at single-cell level Nature 2019 575 512 518 31597160
15 Nakamoto Y Suda T Momoi T Kaneko S Different procarcinogenic potentials of lymphocyte subsets in a transgenic mouse model of chronic hepatitis B Cancer Res 2004 64 3326 3333 15126377
16 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 33574809
17 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 36736292
18 Spiga F Gibson M Dawson S Tilling K Davey Smith G Munafò MR Higgins JPT Tools for assessing quality and risk of bias in Mendelian randomization studies: a systematic review Int J Epidemiol 2023 52 227 249 35900265
19 Emdin CA Khera AV Kathiresan S Mendelian Randomization JAMA 2017 318 1925 1926 29164242
20 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 32929287
21 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 38533503
22 Bourne RR Stevens GA White RA Smith JL Flaxman SR Price H Jonas JB Keeffe J Leasher J Naidoo K Pesudovs K Resnikoff S Taylor HR Vision Loss Expert Group Causes of vision loss worldwide, 1990-2010: a systematic analysis Lancet Glob Health 2013 1 e339 e349 25104599
23 Burgess S Thompson SG CRP CHD Genetics Collaboration Avoiding bias from weak instruments in Mendelian randomization studies Int J Epidemiol 2011 40 755 764 21414999
24 Novobrantseva TI Majeau GR Amatucci A Kogan S Brenner I Casola S Shlomchik MJ Koteliansky V Hochman PS Ibraghimov A Attenuated liver fibrosis in the absence of B cells J Clin Invest 2005 115 3072 3082 16276416
25 Rahman ZS Manser T B cells expressing Bcl-2 and a signaling-impaired BAFF-specific receptor fail to mature and are deficient in the formation of lymphoid follicles and germinal centers J Immunol 2004 173 6179 6188 15528355
26 Do RK Chen-Kiang S Mechanism of BLyS action in B cell immunity Cytokine Growth Factor Rev 2002 13 19 25 11750877
27 Yang C Li N Wang Y Zhang P Zhu Q Li F Han Q Lv Y Yu L Wei P Liu Z Serum levels of B-cell activating factor in chronic hepatitis B virus infection: association with clinical diseases J Interferon Cytokine Res 2014 34 787 794 24811691
28 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 26201785
29 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 29105104
30 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 20308604
31 Dudek M Pfister D Donakonda S Filpe P Schneider A Laschinger M Hartmann D Hüser N Meiser P Bayerl F Inverso D Wigger J Sebode M Öllinger R Rad R Hegenbarth S Anton M Guillot A Bowman A Heide D Müller F Ramadori P Leone V Garcia-Caceres C Gruber T Seifert G Kabat AM Mallm JP Reider S Effenberger M Roth S Billeter AT Müller-Stich B Pearce EJ Koch-Nolte F Käser R Tilg H Thimme R Boettler T Tacke F Dufour JF Haller D Murray PJ Heeren R Zehn D Böttcher JP Heikenwälder M Knolle PA Auto-aggressive CXCR6(+) CD8 T cells cause liver immune pathology in NASH Nature 2021 592 444 449 33762736
32 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 21697460
33 Shen XH Xu P Yu X Song HF Chen H Zhang XG Wu MY Wang XF Discrepant Clinical Significance of CD28(+)CD8(-) and CD4(+)CD25(high) Regulatory T Cells During the Progression of Hepatitis B Virus Infection Viral Immunol 2018 31 548 558 30117787
34 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 23624000
35 Preston JE Segal MB The steady-state amino acid fluxes across the perfused choroid plexus of the sheep Brain Res 1990 525 275 279 2123729
36 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 28777871
37 Tacke RS Lee HC Goh C Courtney J Polyak SJ Rosen HR Hahn YS Myeloid suppressor cells induced by hepatitis C virus suppress T-cell responses through the production of reactive oxygen species Hepatology 2012 55 343 353 21953144
38 Muhanna N Abu Tair L Doron S Amer J Azzeh M Mahamid M Friedman S Safadi R Amelioration of hepatic fibrosis by NK cell activation Gut 2011 60 90 98 20660699
39 Sajid M Liu L Sun C The Dynamic Role of NK Cells in Liver Cancers: Role in HCC and HBV Associated HCC and Its Therapeutic Implications Front Immunol 2022 13 887186 35669776
40 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 33723292
41 O'Keeffe M Mok WH Radford KJ Human dendritic cell subsets and function in health and disease Cell Mol Life Sci 2015 72 4309 4325 26243730
42 Lundahl J Halldén G Hallgren M Sköld CM Hed J Altered expression of CD11b/CD18 and CD62L on human monocytes after cell preparation procedures J Immunol Methods 1995 180 93 100 7534806
43 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 21932384
44 Jhun JY Kim HY Byun JK Chung BH Bae SH Yoon SK Kim DG Yang CW Cho ML Choi JY B-cell-associated immune profiles in patients with decompensated cirrhosis Scand J Gastroenterol 2015 50 884 891 25861705
