The Tohoku Journal of Experimental Medicine
Online ISSN : 1349-3329
Print ISSN : 0040-8727
ISSN-L : 0040-8727
Regular Contribution
Receptor-Interacting Protein Kinase 3 (RIPK3) mRNA Levels Are Elevated in Blood Mononuclear Cells of Patients with Poor Prognosis of Acute-on-Chronic Hepatitis B Liver Failure
Liyan HanYue TengYuchen FanShuai GaoFeng LiKai Wang
Author information
JOURNAL FREE ACCESS FULL-TEXT HTML

2019 Volume 247 Issue 4 Pages 237-245

Details
Abstract

Necroptosis refers to a programmed form of necrosis, which involves the receptor-interacting protein kinase 3 (RIPK3) and mixed lineage kinase domain-like protein (MLKL). In this study, to investigate the role of necroptosis in the pathogenesis of acute-on-chronic hepatitis B liver failure (ACHBLF), we retrospectively analyzed 122 patients with ACHBLF, 131 patients with chronic hepatitis B (CHB), and 35 healthy controls (HCs). Using quantitative real-time polymerase chain reaction (RT-qPCR), we measured RIPK3 mRNA levels in peripheral blood mononuclear cells (PBMCs). ELISA was performed to measure the serum levels of MLKL, TNF-α and caspase-8. We found that RIPK3 mRNA levels were significantly higher in patients with ACHBLF than those with CHB or HCs. RIPK3 mRNA levels in patients with ACHBLF were positively correlated with serum levels of TNF-α or MLKL and negatively correlated with caspase-8 levels. Univariate and multivariate analysis revealed that RIPK3 mRNA level was predictive of 3-month mortality of ACHBLF. The area under receiver operating characteristic curve (AUC) of RIPK3 mRNA levels was 0.810 (95% CI: 0.729-0.876), which was higher than that of MELD scores (0.766, 95% CI: 0.681-0.838). The optimal cut-off point of 8.81 was determined for RIPK3 mRNA levels, which showed a sensitivity of 80.7% and a negative predictive value of 80.4%. These results indicate that elevated RIPK3 mRNA levels in PBMCs are associated with poor prognosis of ACHBLF. We thus propose that necroptosis may play an important role in pathogenesis of ACHBLF.

Introduction

Acute-on-chronic liver failure (ACLF) refers to a rapid and acute deterioration of liver function in patients with chronic liver diseases (Sarin et al. 2014). In Asian, hepatitis B virus (HBV) infection is the main cause of ACLF and acute-on-chronic hepatitis B liver failure (ACHBLF) accounts for more than 70% of ACLF (Gao et al. 2015). ACHBLF usually progress rapidly and result in multiple organ failure within 4 weeks (Katoonizadeh et al. 2010; Jalan et al. 2012). A previous study showed that the short- and medium-term mortality of ACHBLF was about 50-90%, which resulted in the death of more than 120,000 patients every year (Polson et al. 2005). Although liver transplantation remains the only effective therapy for patients with ACHBLF (Fan et al. 2015), it cannot be widely applied because of the shortage of liver donors (Arulraj and Neuberger 2011). Consequently, new prognostic markers are urgently needed for mortality prediction and severity discrimination in patients with ACHBLF (Laleman et al. 2011).

Necroptosis refers to a programmed form of necrosis, which results in programmed cell death via apoptosis (Linkermann and Green 2014). In the process of necroptosis, tumor necrosis factor-α (TNF-α) leads to stimulation of its receptor TNFR1. Then, TNFR1 recruits the TNF receptor-associated death domain (TRADD) and signals to receptor-interacting protein kinase 1 (RIPK1). In the absence of caspase-8, RIPK1 recruits RIPK3 and forms the necrosome (Vanden Berghe et al. 2014). The necrosome then activates mixed lineage kinase domain-like protein (MLKL) through phosphorylation, eventually leading to expulsion of cellular contents into the extracellular space (Wang et al. 2014). Necroptosis was shown to play an important role in viral defense and inflammatory diseases (Gunther et al. 2011). Meanwhile, necroptosis was also found to be involved in the pathogenesis and progression of severe liver diseases (Gunther et al. 2016). We therefore hypothesized that RIPK3 might participate in the pathogenesis of ACHBLF and serve as its prognostic marker.

In this study, we investigated the role of necroptosis in ACHBLF. Using quantitative real-time polymerase chain reaction (RT-qPCR), we measured RIPK3 mRNA levels in peripheral blood mononuclear cells (PBMCs). ELISA was performed to measure serum levels of TNF-α, MLKL and caspase-8. We thus investigated the role of necroptosis in ACHBLF and evaluated the predictive value of RIPK3 mRNA levels in PBMCs for the prognosis.

Methods

Study design and population

One hundred and twenty-two patients with ACHBLF, 131 patients with CHB and 35 healthy controls (HCs) were retrospectively enrolled from July 2014 to November 2015 at the Department of Hepatology, Qilu Hospital of Shandong University. CHB was diagnosed by a positive hepatitis B surface antigen (HBsAg) for more than 6 months before the enrollment (Lok and McMahon 2009). ACHBLF was defined according to the guideline of Asian Pacific Association for the Study of the Liver (APASL) (Sarin et al. 2009).

Exclusion criteria included co-infection with hepatitis A virus, hepatitis C virus, hepatitis D virus, hepatitis E virus, human immunodeficiency virus, autoimmune liver disease, alcoholic liver disease, pregnancy, and hepatocellular carcinoma.

Written informed consents were signed by all participants under protocols approved by the local Research and Ethics Committee at Qilu Hospital of Shandong University, in accordance with the guidelines of the 1975 Declaration of Helsinki.

Isolation of peripheral blood mononuclear cells (PBMCs)

Venous peripheral blood of 5 ml was drawn from each participant at the first day of diagnosis, using EDTA as an anticoagulant agent. PBMCs were isolated by gradient centrifugation from the blood via Ficoll-Plaque Plus (Healthcare, Uppsala, Sweden). Then, PBMCs were washed three times with phosphate-buffered saline and stored at –20°C until use.

RNA extraction and RT-qPCR

Total RNA of PBMCs was extracted using TRIzol (Invitrogen, Carlsbad, CA, USA). The RNA concentration was determined by the Eppendorf Biophotometer (Brinkmann Instruments, Westbury, NY). RNAs were then converted into cDNAs via first-strand cDNA synthesis kit (Fermentas, Vilnius, Lithuania). The RIPK3 mRNA levels were determined by RT-qPCR, which was performed on the Lightcycler 480 (Roche Diagnostic, Mannheim, Germany) with SYBR Green (Toyobo, Osaka, Japan). β-Actin was used as the endogenous control. Amplification was performed in a total volume of 20 μL, which contained 0.5 mM of cDNA, 0.5 mM of each primer and 10 × SYBR Green. The primers were described as follows: RIPK3 forward: 5′-CAAGATCGTAAACTCGAAGG-3′ and reverse: 5′-CCGTTCTCCATGAATTTAGT-3′ (Geserick et al. 2015); β-actin forward: 5′-ATGGGTCAGAAGGATTCCTATGTG-3′, and reverse: 5′-CTTCATGAGGTAGTCAGTCAGGTC-3′. The reaction of PCR was performed as follows: the initial step was 95°C for 30 s, followed by 40 cycles of 95°C for 5 s, 55°C for 30 s, and a final step of 72°C for 30 s. The RIPK3 mRNA levels were calculated by the 2(–∆∆Ct) method.

ELISA for serum levels of TNF-α, MLKL and caspase-8

All blood samples were centrifuged at 2,000 rpm for 5 min to obtain the serum. The serum was then stored at –80°C until subsequent experiments. Serum levels of TNF-α (EIAab, Wuhan, China), caspase-8 (eBioscience, San Diego, CA, USA) and MLKL (Lengton Bioscience Co., Shanghai, China) were determined by ELISA. The lower detection limit of the TNF-α, MLKL and caspase-8 was 3.9 pg/mL, 1 µg/mL and 0.1 ng/mL, respectively.

Clinical and pathological data

The serum biochemical markers included alanine aminotransferase (ALT), aspartate aminotransferase (AST), total bilirubin (TBIL), albumin (ALB), prothrombin activity (PTA), international normalized ratio (INR) and creatinine (Cr). Hepatitis B e antigen (HBeAg) was determined using an automatic analyzer (Roche Diagnostics, Basel, Switzerland). The serum HBV DNA level was quantified using a real-time PCR System (Applied Biosystems, Foster city, CA, USA). All markers were measured according to the standard methodologies at the Department of Laboratory Medicine, Qilu Hospital of Shandong University.

Model for end-stage liver diseases (MELD)

MELD score was calculated by the following formula:

MELD score = 9.57 × loge [creatinine (mg/dl)] + 3.78 × loge [bilirubin (mg/dl)] + 11.2 × loge (INR) + 6.43 × (etiology: 0 if cholestatic or alcoholic, 1 otherwise) (Kamath et al. 2001).

Follow-up of patients with ACHBLF

All the patients with ACHBLF were followed up from the first date of diagnosis. They were followed up for at least 3 months for monitoring the outcome (death or survival).

Statistical analysis

Statistical analyses of the data were performed with SPSS version 19.0 (SPSS Inc., Chicago, IL, USA). Quantitative variables were expressed as median (centile 25; centile 75). Categorical variables were expressed as number (%). Mann-Whitney U-test was performed to compare the quantitative variables. Chi-square test was applied to compare the categorical data. Spearman test was applied for correlation analysis. Diagnostic value of RIPK3 mRNA and MELD score in predicting 3-month mortality for ACHBLF patients was assessed by the area under the receiver operating characteristic curve (AUC). From the receiver operating characteristic (ROC) curve coordinates, the optimal cut-off point associated with the maximum of the Youden index was determined. Youden index-based cut-off point was determined. Sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV) was used to assess the diagnostic accuracy. Survival curves of ACHBLF patients were drawn using the Kaplan-Meier method. All statistical analyses were two-sided and P value < 0.05 was considered statistically significant.

Results

Basic characteristic of the study population

From July 2014 to November 2015, 329 participants (139 ACHBLF patients, 155 CHB patients and 35 HCs) were screened at the Department of Hepatology, Qilu Hospital of Shandong University. Seventeen ACHBLF patients were excluded from this study for co-infection with hepatitis C virus (n = 5), hepatocellular carcinoma (n = 8), alcohol abuse (n = 3) and autoimmune liver disease (n = 1). Twenty-four CHB patients were excluded from this study for co-infection with hepatitis C virus (n = 6), hepatocellular carcinoma (n = 11), alcohol abuse (n = 5) and pregnancy (n = 2), which left 288 participants (122 ACHBLF patients, 131 CHB patients and 35 HCs) (Fig. 1).

Table 1 showed the baseline characteristics of the participants. There was no significant difference between ACHBLF and CHB patients in sex (P = 0.100), age (P = 0.200), HBsAg (P = 0.085), HBeAg (P = 0.310) and log10[HBV DNA] (P = 0.080), respectively. However, significant differences were found between ACHBLF and CHB patients in ALT, AST, TBIL, ALB, PTA, INR and Cr (P < 0.001; respectively). Liver cirrhosis (LC), hepatic encephalopathy (HE) and ascites were found in 65, 39 and 63 patients with ACHBLF. In addition, the 3-month mortality of patients with ACHBLF was 46.72%. After 3-month follow-up, patients with ACHBLF were divided into survivors and non-survivors. The baseline characteristics are shown in Table 2.

Fig. 1.

Flow diagram depicting the participants’ selection process.

Table 1.

The baseline characteristics of enrolled participants.

P value refers to the difference between patients with ACHBLF and CHB.

*Significant difference (P < 0.05).

Table 2.

The baseline characteristics of survivors and non-survivors in patients with ACHBLF.

*Significant difference (P < 0.05).

RIPK3 mRNA levels in PBMCs

We measured RIPK3 mRNA levels in PBMCs of patients with ACHBLF, CHB and HCs. As shown in Fig. 2a, RIPK3 mRNA levels in patients with ACHBLF [9.64 (3.99, 22.68)] were significantly higher than those with CHB [1.81 (1.02, 3.78); P < 0.001] and HCs [0.78 (0.67, 1.35); P < 0.001]. Meanwhile, RIPK3 mRNA levels in patients with CHB were significantly higher than those in HCs (P < 0.001).

Fig. 2.

The association of RIPK3 mRNA levels in PBMCs and clinical parameters.

a. The comparison among RIPK3 mRNA levels in patients with ACHBLF (n = 122), CHB (n = 131) and HCs (n = 35).

b. No significant differences (P = 0.522) were found between RIPK3 mRNA levels in HBeAg(+) group (n = 72) and HBeAg (–) group (n = 50).

c. No significant correlation (P = 0.581) was found between RIPK3 mRNA level and log10 [HBV DNA]. The solid line refers to the best-fit line, and the broken lines refer to the 95% confidence bands of the best-fit line.

d. No significant correlation (P = 0.164) was found between RIPK3 mRNA level and HBsAg.

e. Significant correlation (P < 0.001) was found between RIPK3 mRNA level and TBIL.

f. No significant correlation (P = 0.308) was found between RIPK3 mRNA level and ALB.

g. Significant correlation (P < 0.05) was found between RIPK3 mRNA level and Cr.

h. Significant correlation (P < 0.001) was found between RIPK3 mRNA level and INR.

i. Significant correlation (P < 0.001) was found between RIPK3 mRNA level and MELD score.

Correlation between RIPK3 mRNA levels and clinical parameters in patients with ACHBLF

Among patients with ACHBLF, we compared the RIPK3 mRNA levels in PBMCs between HBeAg (+) and HBeAg (–) patients. No significant difference was found between them (P = 0.522) (Fig. 2b). We also analyzed the correlations between RIPK3 mRNA levels in PBMCs and clinical parameters in patients with ACHBLF. There was no significant correlation between RIPK3 mRNA levels and log10 [HBV DNA] (r = 0.051, P = 0.581), HBsAg (r = 0.127, P = 0.164), or ALB (r = 0.093, P = 0.308) (Fig. 2c, d, f). By contrast, RIPK3 mRNA levels were positively correlated with TBIL (r = 0.377, P < 0.001), Cr (r = 0.214, P < 0.05), INR (r = 0.364, P < 0.001), and MELD scores (r = 0.406, P < 0.001), respectively (Fig. 2e, g, h, i).

Furthermore, we compared RIK3 mRNA levels between patients with LC and those without LC, HE or ascites. Patients with LC, HE or ascites had significantly higher RIPK3 levels than those without (P < 0.05, respectively; Fig. 3).

Fig. 3.

The correlations between RIPK3 mRNA levels and complications of ACHBLF.

a. RIPK3 mRNA levels were significantly higher (P < 0.05) in patients with LC (n = 65) than those without (n = 57).

b. RIPK3 mRNA levels were significantly higher (P < 0.05) in patients with hepatic encephalopathy (HE) (n = 39) than those without (n = 83).

c. RIPK3 mRNA levels wre significantly higher (P < 0.05) in patients with ascites (n = 63) than those without (n = 59).

Serum levels of TNF-α, MLKL and caspase-8

We also measured the serum levels of TNF-α, MLKL and caspase-8 (Fig. 4a-c). We found that the TNF-α levels in patients with ACHBLF [21.81 (12.51-36.90)] were significantly higher than those with CHB [9.34 (4.83-15.36), P < 0.001] and HCs [7.38 (5.29-8.83), P < 0.001]. There was no significant difference between patients with CHB and HCs (P = 0.058). The MLKL levels were significantly higher in patients with ACHBLF [45.71 (26.03-70.36)] than those with CHB [17.00 (11.21-29.31), P < 0.001] and HCs [11.01 (7.26-13.85), P < 0.001]. Meanwhile, Patients with CHB showed significantly higher MLKL levels than HCs (P < 0.001). In addition, patients with ACHBLF [0.38 (0.22-0.74)] showed decreased levels of caspase-8 than those with CHB [0.78 (0.58-0.97), P < 0.001] and HCs [0.67 (0.55-0.97), P < 0.001]. However, no significant difference was found between patients with CHB and HCs (P = 0.705).

Importantly, RIPK3 mRNA levels in patients with ACHBLF were positively correlated with TNF-α levels (r = 0.412, P < 0.001) and MLKL levels (r = 0.625, P < 0.001), respectively. By contrast, RIPK3 mRNA levels were negatively correlated with caspase-8 levels (r = –0.279, P < 0.05) (Fig. 4d-f).

Fig. 4.

The relationships between RIPK3 and RIPK3-associated cytokines.

a. The comparison among TNF-α levels in patients with ACHBLF, CHB and HCs.

b. The comparison among MLKL levels in patients with ACHBLF, CHB and HCs.

c. The comparison among caspase-8 levels in patients with ACHBLF, CHB and HCs.

d. Significant correlation (P < 0.001) were found between RIPK3 mRNA levels and TNF-α levels.

e. Significant correlation (P < 0.001) were found between RIPK3 mRNA levels and MLKL levels.

f. Significant correlation (P < 0.05) was found between RIPK3 mRNA levels and caspase-8 levels.

RIPK3 mRNA levels as an independent predictive factor for the outcomes of ACHBLF patients

Univariate and multivariate regression analysis were performed to investigate the association between clinicopathological parameters and the outcomes of ACHBLF patients (Table 3). Univariate analysis identified that the outcomes of ACHBLF patients were significantly associated with TBIL (OR = 1.007, 95% CI: 1.004-1.010, P < 0.05), INR (OR = 4.269, 95% CI: 2.254-8.086, P < 0.05), PTA (OR = 0.865, 95% CI: 0.811-0.922, P < 0.05), ascites (OR = 2.775, 95% CI: 1.329-5.795, P < 0.05), LC (OR = 2.445, 95% CI: 1.175-5.087, P < 0.05), HE (OR = 4.740, 95% CI: 2.065-10.877, P < 0.05), RIPK3 mRNA levels (OR = 1.106, 95% CI: 1.063-1.151, P < 0.05), or MELD score (OR = 1.207, 95% CI: 1.110-1.312, P < 0.05). Then, multivariate regression analysis was performed with those variables. The RIPK3 mRNA level (OR = 1.082, 95% CI: 1.037-1.129, P < 0.05) was identified as an independent predictive factor for the outcome of ACHBLF patients.

Table 3.

Univariate and multivariate logistic analysis of clinical parameters for predicting the outcome of ACHBLF.

*Significant difference (P < 0.05).

Predictive value of RIPK3 mRNA levels for prognosis of ACHBLF

In this study, the 3-month mortality of ACHBLF patients was 46.7% (57/122). After 3-month follow-up, patients with ACHBLF were then divided into survivors and non-survivors. The 3-month mortality was significantly associated with ALT (P = 0.002), TBIL (P < 0.001), PTA (P < 0.001), INR (P < 0.001), LC (P = 0.016), HE (P < 0.001), ascites (P = 0.006) or MELD scores (P < 0.001). No significant association was found between 3-month mortality and gender (P = 0.480), age (P = 0.210), HBsAg (P = 0.633), HBeAg (P = 0.616), Log10 [HBV DNA] (P = 0.460), AST (P = 0.277), ALB (P = 0.226) or Cr (P = 0.125) (Table 2). Meanwhile, The RIPK3 mRNA levels in survivors were significantly higher than those in non-survivors (P < 0.001) (Fig. 5a).

Fig. 5b showed Kaplan-Meier survival analysis for ACHBLF patients with RIPK3 mRNA levels above 8.81 or not. Compared with patients with RIPK3 mRNA levels ≤ 8.81, survival rate in those with RIPK3 mRNA levels > 8.81 were significantly lower (log rank analysis, P < 0.001). The mean survival time for ACHBLF patients with RIPK3 mRNA levels ≤ 8.81 or > 8.81 was 78.750 ± 3.232 days or 38.820 ± 4.395 days. Then, we performed the receiver operating characteristic (ROC) curve to evaluate the predictive value of RIPK3 mRNA levels for predicting 3-month mortality in patients with ACHBLF. As shown in Fig. 5c, the AUC of RIPK3 mRNA levels (0.810, 95% CI: 0.729-0.876) was higher than that of MELD scores (0.766, 95% CI: 0.681-0.838). With a cut-off point of 8.81, RIPK3 mRNA levels had a sensitivity of 80.7%, a specificity of 69.2%, a positive predictive value of 69.7%, and a negative predictive value of 80.4%.

Fig. 5.

ROC curve and the Kaplan–Meier graph of RIPK3 mRNA level in predicting 3-month mortality of patients with ACHBLF.

a. RIPK3 mRNA levels were significantly higher (P < 0.001) in non-survivors (n = 57) than survivors (n = 65).

b. Kaplan-Meier graph showing survival probability at the end of 3-month follow-up in ACHBLF patients with RIPK3 mRNA levels > 8.81 and ≤ 8.81.

c. ROC curves of RIPK3 mRNA level and MELD score in predicting 3-month mortality of ACHBLF patients.

Discussion

ACHBLF is characterized by its high morbidity and mortality. Accurate and effective noninvasive markers to predict the prognosis of ACHBLF are urgently needed (Wlodzimirow et al. 2013; Zheng et al. 2013). In this study, we found that RIPK3 mRNA levels were significantly higher in patients with ACHBLF than those with CHB or HCs. Meanwhile, the RIPK3 expression in ACHBLF patients was positively correlated with serum TBIL, INR and MELD score, which were often used to predict liver injury. A previous study already demonstrated that TBIL and INR were positively associated with mortality in patients with ACLF (Wlodzimirow et al. 2013). Meanwhile, MELD score was widely accepted as a model to evaluate short-term mortality in patients with end-stage liver diseases (Kamath et al. 2001; Northup et al. 2015).

In this study, we found that RIPK3 mRNA levels in PBMCs were associated with poor prognosis of ACHBLF. Therefore, it might serve as a prognostic marker in further clinical work. Meanwhile, patients with LC, HE or ascites had significantly higher RIPK3 mRNA levels than those without. A previous study also indicated that cirrhosis, HE as well as ascites were associated with adverse outcomes in patients with ACHBLF (Seto et al. 2012).

It was reported that necroptosis might participate in the pathogenesis of acute liver failure (ALF) and was associated with its prognosis (Rastogi et al. 2011). TNF-α, MLKL and caspase-8 were essential for the process of necroptosis (Holler et al. 2000). TNF-α is a cytokine involved in systemic inflammation. In fulminant liver failure, the serum level of TNF-α was significantly increased correlated with serious complications, such as infection and HE (de la Mata et al. 1990; Zhao et al. 2015). In the process of necroptosis, MLKL was phosphorylated and insert into the plasma membranes and organelles, which resulted in the expulsion of cellular contents into the extracellular space. Caspase-8 could negatively regulate the RIPK3 dependent necroptosis and participate in the pathogenesis of viral infection and inflammation (Remijsen et al. 2014). In this study, we found that serum TNF-α, MLKL levels were significantly elevated in ACHBLF. Serum caspase-8 level was significantly decreased in ACHBLF. Meanwhile, RIPK3 mRNA levels in patients with ACHBLF were positively correlated with TNF-α, MLKL levels and negatively correlated with caspase-8 levels. We, therefore, demonstrated that necroptosis might play an important role in the pathogenesis of ACHBLF.

Because of the coagulopathy, liver biopsies were hard to be performed in patients with ACHBLF. PBMCs were easy to be obtained from venous peripheral blood and suitable to be used as non-invasive markers in clinical work. Therefore, we evaluated the RIPK3 mRNA levels in PBMCs in this study. We demonstrated the potential predictive value of RIPK3 mRNA levels in PBMCs for prognosis of ACHBLF. However, the RIPK3 mRNA levels in PBMCs might not reflect the expression levels in hepatocytes. The expression of RIPK3 in hepatocytes and its potential role in the necroptosis of ACHBLF still need to be analyzed in further studies.

There were also several limitations in this study. Firstly, the molecular mechanism that how RIPK3 was involved in the development and progression of ACHBLF remained unclear and might be investigated in our further study. Secondly, the enrolled participants were from a single unit and data from large-scale, multi-center cohort might be more valuable. Therefore, these findings still need further validation prior to its clinical usage.

In conclusion, this study revealed that elevated RIPK3 mRNA levels in PBMCs were associated with poor prognosis of ACHBLF. We thus propose that necroptosis may play an important role in the pathogenesis of ACHBLF.

Acknowledgments

This work was supported by the grants from the Key Project of Chinese Ministry of Science and Technology (2017ZX10202202, 2018ZX10302206-003-003, 2018ZX10302206-006-002), National Natural Science Foundation of China (81171579, 81201287, 81371832, 81600494), Natural Science Foundation of Shandong Province (ZR2016HQ42), Science and Technology Development Plan of Shandong Province (2014GSF118068), and the Fundamental Research Funds of Shandong University-Clinical Research Project of Qilu Hospital (2014QLKY11).

Conflict of Interest

The authors declare no conflict of interest.

References
 
© 2019 Tohoku University Medical Press
feedback
Top