Impact Factor 4.137 | CiteScore 4.28
More on impact ›

Original Research ARTICLE

Front. Oncol., 28 February 2020 | https://doi.org/10.3389/fonc.2020.00251

Effect of Non-alcoholic Fatty Liver Disease on the Risk of Synchronous Liver Metastasis: Analysis of 451 Consecutive Patients of Newly Diagnosed Colorectal Cancer

  • Department of Oncology, The Affiliated Zhongda Hospital of Southeast University, Medical School of Southeast University, Nanjing, China

Background: The purpose of this study was to investigate the effect of non-alcoholic fatty liver disease (NAFLD) on the risk of synchronous colorectal liver metastasis (synCRLM).

Methods: A retrospective analysis was performed on 451 consecutive patients with newly diagnosed colorectal cancer (CRC) from January 2014 to January 2019. According to the presence of NAFLD, the CRC patients were divided into two groups, NAFLD group (60 cases) and the control group (391 cases). The clinicopathological features and the prevalence of synCRLM between the two groups were compared. Logistic regression analysis was used to analyze the risk factors of synCRLM. Different non-invasive liver fibrosis scoring models were used to evaluate the effect of advanced fibrosis and cirrhosis stage in NAFLD on the prevalence of synCRLM.

Results: The prevalence of synCRLM was significantly higher in patients with NAFLD than that in patients without NAFLD (18.33 vs. 7.42%; χ2 = 7.669, P = 0.006). A logistic regression analysis indicated that NAFLD, CEA, CA19-9, and lymph node status were risk factors for synCRLM, and NAFLD showed the highest hazard ratio (3.930 [95% confidence interval: 1.616 ~ 9.560]). In NAFLD patients, both fibrosis-4 index (FIB-4) and NAFLD fibrosis score (NFS) were significantly lower in those with synCRLM compared to those without synCRLM [FIB-4: 1.246 (0.833 ~ 1.276) vs. 1.436 (1.016 ~ 2.699), Z = −2.130, P = 0.033; NFS: −1.282 (−2.407 ~ −0.262) vs. −0.255 (−1.582 ~ 0.755), Z = −2.302, P = 0.021; Mann-Whitney test].

Conclusion: NAFLD may be associated with increased liver metastasis, and for NAFLD-related advanced liver fibrosis and cirrhosis may be associated with reduced synchronous liver metastasis in CRC patients. However, the correlation between simple steatosis and steatohepatitis remains to be further determined. Certain factors such as NAFLD, lymph node metastasis, elevated levels of preoperative CEA and CA19-9 are suggesting a high risk of synCRLM.

Introduction

Colorectal cancer (CRC) is one of the most well-established malignancies of the digestive system, which is also the third most common malignant tumor in the world and the fourth most frequent cause of cancer-related deaths. It is estimated that by 2030, this disease burden will increase by nearly 60% (1). There will be 2.2 million new CRC patients and 1.1 million deaths from CRC (2). Distant metastasis is the most important independent risk factors for poor prognosis in CRC patients. The liver is the most common site for distant metastasis of CRC and is often the only organ involved (3). Previous studies have reported that more than 50% of CRC patients will develop liver metastasis in the process of the disease, about 14–25% of CRC patients are found to have synchronous liver metastasis at the time of diagnosis, and about 10–25% of CRC patients develop metachronous liver metastasis during follow-up after the initiation of treatment. Approximately 80–90% of CRC patients with liver metastasis could not initially get radical resection of liver metastases, and these people who unable to get radical resection have a very low 5-year survival rate (48). Therefore, whether liver metastasis occurs or not is the crucial point affecting the prognosis of CRC, and understanding the risk factors of liver metastasis may be one of the effective strategies to reduce the circumstance of liver metastasis.

At present, with the change of lifestyle and dietary structure, the prevalence of insulin resistance and obesity is gradually increasing, and the incidence of adult non-alcoholic fatty liver disease (NAFLD) is moderately increasing, which has become one of the most recurrent causes of chronic liver diseases across the globe. The worldwide prevalence of NAFLD is ~25.24% (9). NAFLD is a chronic and persistent pathological process from hepatic steatosis to obvious hepatic injury. It includes two pathological diagnoses with different prognosis: non-alcoholic fatty liver (NAFL) and non-alcoholic steatohepatitis (NASH) (10). Current studies have shown that NAFL is a disease with or without progression of histopathological changes very slowly over time. NAFL is a relatively benign stage for NAFLD which is easy to be reversed, roughly 10–20% of NAFL can be converted into NASH (11, 12). NASH has a potentially invasive progression process that can lead the way to liver fibrosis, cirrhosis, hepatocellular carcinoma and even induce distant organs carcinogenesis, such as CRC (13, 14).

Most studies have reported that NAFLD is an independent risk factor for the development of CRC (15, 16), while the effect of NAFLD on colorectal liver metastasis (CRLM) remains unclear. Some studies have outlined that NAFLD reduced the risk of CRLM (1720), while others suggested the opposite conclusion (2123). Updated studies that further investigate the mechanism underlying the effects of NAFLD on CRLM are certainly needed. Referring to the well-known “seed-soil” hypothesis in the mechanism of tumor metastasis (24, 25). Herein, CRC cells as “seed” and the liver with NAFLD as “soil” have been taken into consideration to explore whether “soil,” the liver with NAFLD is responsible for metastasis. Our strict definition of “synchronous” ruled out the effect of variations in the follow-up treatment on the likelihood of liver metastasis. Meanwhile, we further analyzed the effects of different stages of NAFLD on synchronous colorectal liver metastasis (synCRLM), which is the first known study to evaluate the impact of the dynamic progression of NAFLD on synCRLM. In this retrospective clinical study, we obtained unique and interesting results and proposed new insights into the association between NAFLD and synCRLM.

Materials and Methods

Study Population

We conducted a retrospective study of all newly diagnosed CRC patients who were continuously admitted to the Affiliated Zhongda Hospital of Southeast University between January 2014 and January 2019.

The inclusion criteria were as follows: (1) all patients confirmed as CRC histopathologically, (2) did not receive any surgery, interventional treatment, chemotherapy or radiotherapy in other hospitals before admission, (3) for the purpose of detecting the presence of NAFLD and distant organ metastasis, all patients screened by preoperative chest, abdominal and pelvic regions imaging, such as ultrasound, computed tomography (CT) and magnetic resonance imaging (MRI), (4) complete medical records.

The exclusion criteria were as follows: (1) CRC patients with synchronous distant metastases except for liver, (2) CRC patients accompanied by malignant tumors in other parts of the body, (3) CRC patients with liver lesions caused by alcohol, viruses, drugs, genetics and other pathogenic factors. In total, 451 patients were included in this study.

The study was approved from an ethics approval by the Ethics Committee of Zhongda Hospital, Southeast University. The study protocol protected the private information of enrolled patients in accordance with the provisions of the Helsinki Declaration.

Data Characteristics

Data collection included demographic information [gender, age, height, weight, body mass index (BMI)] and clinical data, including the following: NAFLD, hypertension, diabetes or impaired fasting glucose (IFG), hepatitis B surface antigen (HBsAg), aspartate aminotransferase (AST), alanine aminotransferase (ALT), glutamyl transferase (GGT), alkaline phosphatase (ALP), total bilirubin (TBIL), direct bilirubin (DBIL), indirect bilirubin (IBIL), triglyceride (TG), platelet (PLT), albumin (ALB) and the following tumor markers: carcinoembryonic antigen (CEA), carbohydrate antigen 19-9 (CA19-9) and alpha-fetoprotein (AFP). Tumor characteristic which based on the guideline of the American Joint Committee on Cancer staging, with detailed descriptions on each patient's pathology report, including primary tumor site, tumor size, tumor type (protuberant, ulcerative, infiltrative), differentiation (well, moderate, poor), tumor (T) status, lymph node (N) status, TNM stage, vascular invasion, nerve invasion, and KRAS, NRAS, BRAF mutation status. The original data for all patients in this study were shown in Supplementary Data Sheet 1.

Imaging Evaluation of the Liver, Metastasis Assessment

Compared with liver biopsy as the gold standard, liver imaging examination is more commonly used in clinical work. Ultrasound, CT and MRI have high sensitivity and specificity in the diagnosis of liver steatosis (26, 27). According to the Chinese guidelines for the prevention and treatment of NAFLD (2018 updated edition) (28), the diagnosis of NAFLD should meet the following requirements: (1) liver biopsy indicates significant hepatic steatosis and/or imaging examination, such as ultrasound, CT and MRI results consistent with the manifestations of fatty liver, which liver fat content is more than 5%, (2) no history of alcohol consumption more than 30 g for men and 20 g for women per day, and (3) rule out other factors or diseases that contribute to hepatic steatosis.

Based on the imaging examination results, patients were divided into the NAFLD group and control (without NAFLD) group. The baseline clinicopathological parameters and prevalence of synCRLM were compared between the two groups. SynCRLM was defined as the synchronous or prior diagnosis of liver metastasis together with primary CRC. All diagnoses of liver metastasis were confirmed by pathological biopsy or independently confirmed by CT or MRI by two senior radiologists.

Calculation of Non-invasive Liver Fibrosis Scoring Models

NAFLD fibrosis score (NFS), Fibrosis-4 Index (FIB-4), AST-to-PLT ratio index (APRI), and BRAD score were calculated and used as the indicators of advanced liver fibrosis or cirrhosis. These scores are calculated as follows:

NFS=1.675+0.037×age (years)+0.094×BMI (Kgm2)                + 1.13×IFG or diabetes(yes=1,no=0)                + 0.99×(AST÷ALT)0.013×PLT (109L)                 0.66×ALB (gdl)
FIB4=age (years)×AST(IUL)÷[PLT (109L)                      ×ALT(IUL)]
APRI=AST(IUL)÷upper limit of the liver biopsy of AST                   ×100÷PLT(109L)

and the upper limit of normal of AST in our hospital is 37IU/L. BRAD score combines the three variables into one weight sum (BMI ≥28 = 1 point, AST/ALT ratio ≥0.8 = 2 points, diabetes = 1 point), which can be used to predict advanced fibrosis quickly and easily.

As reported previously, those patients with NFS higher than 0.676, FIB-4 higher than 1.30, APRI higher than 0.50 and the BRAD score of 2–4 were diagnosed as having an advanced fibrotic liver (2933), patients in NAFLD group were grouped according to high-level (>0.676) or low-level (≤0.676) NFS, high-level (>1.30) or low-level (≤1.30) FIB-4, high-level (>0.50), or low-level (≤0.50) APRI and high-score (2–4 score) or low-score (0–1 score) BRAD. NFS, FIB-4, APRI, and BARD score were compared between the NAFLD patients with or without synCRLM.

Statistical Methods

The correlations between demographics, clinical and pathological data of the two groups were analyzed. Categorical variables were compared by χ2 test, correction χ2 test or Fisher's exact test. Numerical variables conforming to normal distribution were compared using student t test, comparisons for numerical variables with skewed distribution were performed using the Mann-Whitney U test. Significant risk factors for synCRLM were analyzed first by univariate logistic regression analysis and then by multivariate logistic regression analysis. All the statistical tests considered two-sided P value <0.05 as statistically significant. Statistical analysis was performed using SPSS version 25.0 software (IBM Corporation, Armonk, NY, USA).

Results

Baseline Parameters of CRC Patients

A total of 451 patients were confirmed for the analysis during the study period. Among them, 60 (13.30%) patients were diagnosed with NAFLD, and 391 (86.70%) patients were regarded as the control group. The baseline clinicopathological parameters of the two groups are presented in Table 1. The weight and BMI of the NAFLD patients were significantly higher than that of the control patients (weight: P = 0.022; BMI: P < 0.001). NAFLD was found at a higher incidence in patients with diabetes or IFG (41.67 vs. 19.69%, P < 0.001). There were no significant differences in the sex, age, height, hypertension, HBsAg, primary tumor site, tumor size, tumor type, tumor differentiation, T status, LN status, vascular invasion, nerve invasion, and KRAS, NRAS, BRAF mutation status. The prevalence of synCRLM was 18.33% (11/60) in the NAFLD group, which was significantly higher than the prevalence of 7.42% (29/391) in the control group (χ2 = 7.669, P = 0.006). The overall primary disease stage (TNM) was different between the two groups (χ2 = 7.939, P = 0.047), but there was no significant difference between the two groups during stage I to III (χ2 = 0.267, P = 0.862), while there was a significant difference between stage I~III and IV (χ2 = 7.669, P = 0.006), which was attributed to the difference in distant metastasis (M) status between the two groups. Figures 1, 2 showed enhanced CT and enhanced MRI images of liver metastasis, NAFLD and normal liver in CRC patients in this study, respectively. Figure 3 showed the histopathological manifestation of resection of liver metastasis in a CRC patient in this study.

TABLE 1
www.frontiersin.org

Table 1. Clinicopathological parameters of primary colorectal cancer in the NAFLD group and control group.

FIGURE 1
www.frontiersin.org

Figure 1. Enhanced CT images of liver metastasis, NAFLD and normal liver in the CRC patients. (A) NAFLD group, synCRLM+: The shape of the liver was regular, and the proportion of each lobe was within the normal range. The density of the liver parenchyma was lower than that of the spleen in each stage of enhancement, and multiple round-like enhancement shadows could be seen in the liver parenchyma with unclear boundary. (B) Control group (without NAFLD), synCRLM+: A patchy abnormal enhancement foci with unclear boundary was seen in the left inner lobe of the liver, showing slight enhancement, while no abnormal density or enhancement was observed in the remaining liver parenchyma. (C) NAFLD group, synCRLM: The density of liver parenchyma decreased, which was lower than that of the spleen in the same layer. (D) Control group (without NAFLD), synCRLM: No abnormal density and enhancement foci were found in the liver parenchyma.

FIGURE 2
www.frontiersin.org

Figure 2. Enhanced MRI images of liver metastasis, NAFLD and normal liver in the CRC patients. (A) NAFLD group, synCRLM+: The density of liver parenchyma decreased significantly in the reverse position. Multiple patchy abnormal signal shadows were seen in the liver, especially in the hepatic hilum. Necrosis was observed in the center of some lesions, with low enhancement in the center of the enhanced lesions and line-like enhancement at the edges. (B) Control group (without NAFLD), synCRLM+: Enhanced scan showed a round-like long T1 and long T2 signal shadow near the right hepatic vein in the right posterior lobe of the liver. In the arterial phase, the lesion was significantly enhanced, and the surrounding areas could be seen with enhanced flocculent perfusion. (C) NAFLD group, synCRLM: The density of liver parenchyma decreased significantly in the reverse position. (D) Control group (without NAFLD), synCRLM: The liver was in regular shape, and the proportion of each lobe was harmonious. No obvious abnormal signal shadow was observed in the liver parenchyma.

FIGURE 3
www.frontiersin.org

Figure 3. Hepatic histopathology of a case of CRC patient. Accumulation of fat vacuoles (steatosis) and abnormal tumor cells were observed in hepatic lobules. Immunohistochemical result: CK7 (focal +), CK20 (–), CDX2 (+), villin (+), Ki67 (about 60%+), SATB2 (+), CDH17 (+). The liver metastasis was consistent with the origin of colorectal adenocarcinoma.

Risk Factors for synCRLM

Among all patients, 40 CRC patients were diagnosed with synchronous liver metastasis, and the total prevalence of synCRLM was 8.87%. According to whether the occurrence of synCRLM, they were divided into synCRLM+ and synCRLM group. Chi-square test, independent sample t test and non-parametric test were performed on the clinicopathological parameters of the two groups, and the results showed that there were no significant differences in HBsAg, AFP, ALT, AST, ALP, GGT, TBIL, DBIL, IBIL, ALB, TG, PLT, tumor size, tumor type, nerve invasion, and NRAS, BRAF mutation status between the two groups (P > 0.05) (Supplementary Table 1). Following univariate logistic regression analysis, CEA, CA19-9, primary tumor site, differentiation, T status, LN status, vascular invasion and NAFLD were chosen for the subsequent multivariate logistic regression analysis. As the loss of KRAS mutation status was more than 50%, it was not included in the multivariate logistic regression analysis, which finally indicated that NAFLD was considered an independent risk factor for liver metastasis and NAFLD had the highest hazard ratio (3.930[95% confidence interval: 1.616 ~ 9.560]). Moreover, advanced LN status (N1 and N2), elevated levels of preoperative CEA and CA19-9 were linked with increased risks of synCRLM (Table 2). The sensitivity of CEA and CA19-9 in predicting synCRLM were 47.5 and 40.0%, respectively; the specificities were 95.6 and 94.6%, respectively. ROC curve analysis displayed that the optimal critical value of CEA was 34.86 ng/mL and CA19-9 was 53.82 U/mL, which suggested that preoperative CEA levels above 34.86 ng/mL and preoperative CA19-9 levels above 53.82 U/mL were associated with greater risks of synCRLM. The sensitivity and specificity of the combined diagnosis of CEA and CA19-9 for synCRLM were 57.5 and 92.2%, respectively. The combination of CEA and CA19-9 (Area under curve (AUC) [95%CI], 0.739 [0.642~0.835]) exhibited a better predictive value than CEA (AUC [95%CI], 0.705 [0.603~0.807]) or CA19-9 (AUC [95%CI], 0.649 [0.541~0.835]) alone (Figure 4).

TABLE 2
www.frontiersin.org

Table 2. Univariate and multivariate logistic regression analysis of the significant predictors for synchronous colorectal liver metastasis.

FIGURE 4
www.frontiersin.org

Figure 4. Receiver operating characteristic (ROC) analysis for the prediction of synchronous liver metastasis in colorectal cancer patients.

The Influence of Advanced Liver Fibrosis and Cirrhosis on synCRLM

Sixty patients with NAFLD were divided into two groups according to whether synchronous liver metastasis occurred. The median (25–75th percentile) NFS and FIB-4 were −1.282 (−2.407 ~ 0.262) and 1.246 (0.833 ~ 1.276) in those with synCRLM, both were significantly lower compared to the median of −0.255 (−1.582 ~ 0.755) and 1.436 (1.016 ~ 2.699) in those without synCRLM (NFS: Z = −2.302, P = 0.021, Mann-Whitney test; FIB-4: Z = −2.130, P = 0.033, Mann-Whitney test). There was no significant difference in APRI between the two groups (Z = −1.939, P = 0.052, Mann-Whitney test). FIB-4, APRI, NFS, and BRAD score of NAFLD patients were grouped according to the recommended threshold values, respectively, to compare the prevalence of synCRLM in each group at different levels (Table 3). These results demonstrated that the risk of synCRLM in NAFLD patients with advanced fibrosis/cirrhosis was significantly lower than patients without advanced fibrosis/cirrhosis.

TABLE 3
www.frontiersin.org

Table 3. Comparison of synCRLM between different levels of FIB-4, APRI, NFS, and BRAD score.

Discussion

In our study, we found that NAFLD may be associated with increased synCRLM. However, NAFLD is a general term for multiple disease stages, including NAFL, NASH, NAFLD-related fibrosis, cirrhosis and even hepatocellular carcinoma. We need to analyze the different effects that may occur at different stages in the progression of NAFLD. We observed that when NAFLD progressed to advanced fibrosis or cirrhosis were associated with decreased synCRLM. Meanwhile, we used multivariate analysis to suggest that NAFLD, CEA, CA19-9, and LN status were associated with an increased incidence of synchronous liver metastasis in CRC patients.

The “seed-soil” hypothesis suggests that metastatic tumor cells will migrate to an area where the local microenvironment is favorable. Most prior studies have focused on determining how “seeds” (cancer cells) promote metastasis, while ignoring the role of “soil” (target organs). The preferential growth of metastasis in an organ is due to its unique biological characteristics and its endogenous microenvironment containing special cells and molecular components that promote the formation of metastasis. Interventions that target blocking these signals may be effective in inhibiting metastasis. Our study aimed at the “soil” factor, re-recognized the impact of NAFLD on synCRLM and highlighted the role of liver microenvironment changes in the formation of liver metastases.

It is notable that the phenomenon of liver metastasis is one of the main aspects determining the prognosis of CRC patients. We need to pay more attention to the moment when liver metastasis occurs. Synchronous liver metastasis may represent worse biological behavior of the primary tumor and even worse prognosis. Although with the continuous deepening of the medical research and the rapid development of medical technology, the diagnosis and treatments of CRC are becoming more and more standardized, but often because of a variety of reasons, including differences of medical services levels in different regions, limitations of patients' family conditions, differences of patients' individual physical conditions and compliance, etc., there is still a significant imbalance in the treatment of CRC. At the same time, the standard of diagnosis and treatments of CRC are constantly upgrading. Even if the guidelines are strictly followed, the standard treatments received by patients in different periods are not alike. CRC patients who receive curative resection will have a longer expected survival time and a longer follow-up time. While, the disparity in the treatment among CRC patients, there is no reasonable explanation for the difference in the incidence of metachronous liver metastasis, which is not conducive to truly understanding the impact of NAFLD on CRLM. Based on the above situation, it is bound to affect the final outcomes if the synchronous and metachronous liver metastasis are not strictly distinguished. In our study, the prevalence of synCRLM was 8.87%. However, previous clinical studies have reported that the prevalence of synCRLM is about 14–25% (4–8). The decreased prevalence of synCRLM in our study may be related to inconsistencies in the definition of “synchronous.” Previous studies have defined synchronous as the development of a liver metastasis within 6 months and even up to 12 months after the primary CRC has been diagnosed (34). In addition, in recent years, the general population has gradually deepened their understanding of CRC, as well as the progress of early screening and diagnosis technology of CRC. Therefore, the decrease in the incidence of liver metastasis is within the expected range. However, our study defined as liver metastasis detected concurrently or before the primary CRC with a strict definition of synchronous can exclude the effect of subsequent treatments on liver metastasis and reflect the differences of cancer biology more accurately, providing guidance for individualized treatment.

In this study, the prevalence of NAFLD was 13.30%, and the median prevalence of NAFLD in China was 15.3% (11.3–24.6%), which was in line with the prevalence trend of NAFLD in China (35). Based on previous studies, there was an independent interrelation between NAFLD and both diabetes or IFG, and this relationship was closely related to BMI changes (36, 37), which was consistent with this observed trend in our study. Compared with the control group, CRC patients in NAFLD group had a significantly higher rate of synchronous liver metastasis. Further multivariate logistic regression analysis proposed that NALFD was associated with an increased incidence of synchronous liver metastasis.

A meta-analysis included a total of 10,349 CRC patients from 10 studies, the results demonstrated that patients with chronic liver diseases including hepatitis, cirrhosis, and fatty liver had significantly lower incidences of CRLM than those with the normal livers conditions (OR = 0.32, 95% CI 0.26–0.38, P < 0.001). However, in this meta-analysis, synchronous and metachronous liver metastasis were included. Only 3 studies explored the CRLM correlation with fatty liver, and the correlation with NAFL or NASH was not further studied in detail (17). A study explored the relationship between hepatic steatosis (HS) with the incidence of CRLM by objective quantifying steatosis using a CT scan. Patients with liver-spleen attenuation ratio lower than 1.1 were objectively defined as HS, and the results manifested that the incidence of synchronous metastasis in liver observed in HS patients was lower than that in normal liver patients, but there was no statistically significant difference (3.2 vs. 9.6%, P = 0.06) (18). Some retrospective studies suggested that liver metastasis of CRC was less frequently detected in patients with fatty liver, both of these studies did not distinguish between synchronous and metachronous metastasis and were reported a long time ago, whereas fatty liver was diagnosed only by two-dimensional ultrasound (19, 20). These prior studies underestimated the risk of NAFLD for liver metastasis. However, recent studies reported contradictory results that NAFLD patients have a higher risk of liver metastasis and postoperative recurrence of liver metastasis when compared with non-NAFLD patients, although these studies were looking at metachronous liver metastasis (2123).

Some animal research also presented a different conclusion. The murine models revealed that hepatic steatosis induced by a high-fat diet increased the metastatic tumor burden of colon cancer tumors metastatic to the liver, and the changes in the liver microenvironment caused by fatty liver established a more sensitive microenvironment for metastasis when compared with normal liver (3840). However, the opposite conclusion was discovered in a rat model (41). This difference might be attributed to the application of various experimental conditions, such as different tumor cell lines, atypical types of mice, unlike injection sites and methods, etc. And it was still necessary to pay attention to the different effects of the etiology and degree of hepatic steatosis on the results.

To farther explore the disparities in the prevalence of synCRLM in a different stage of NAFLD, we compared FIB-4, APRI, NFS, and BRAD score in NALFD patients with and without synCRLM. Compared with liver biopsy, FIB-4, APRI, NFS, and BRAD score may have moderate diagnostic efficacy for advanced fibrosis or cirrhosis, but these scoring models are non-invasive, inexpensive, widely available, and the dynamic change process of liver fibrosis can be understood (2933). NFS and FIB-4 scoring were most accurately determined severe liver diseases compared with APRI and BRAD scoring (42). Our results demonstrated that both FIB-4 and NFS scoring in NAFLD patients without synCRLM were significantly higher than that in NAFLD patients with synCRLM. Meanwhile, the prevalence of synCRLM was lower in the high-level FIB-4 and the high-level NFS group than in the low-level FIB-4 and the low-level NFS group. APRI and BARD score did not reach the statistical significance. Therefore, the possible explanation for the above results is that advanced fibrosis or cirrhosis may decrease the risk of synCRLM.

Indeed, earlier studies had recommended that cirrhosis reduces the risk of liver metastasis (43, 44). To be specific, our study observed that when NAFLD progressed to advanced fibrosis or cirrhosis, it may be associated with a decrease in synchronous liver metastasis. We consider several hypotheses that might explain this phenomenon: (1) It may be related to morphological changes in the liver. Gradually accumulation of triglycerides in the liver leads to the compression and deformation of intrahepatic blood vessels. With the aggravation of the severity of the disease, when it develops to advanced fibrosis or cirrhosis, leading to the distortion and reconstruction of intrahepatic blood vessels. Especially in cirrhosis, scar tissue and regenerative nodules are formed, which directly compress the blood vessels and shunt the hepatic artery and portal vein into the central vein, thereby affecting the blood supply in the liver, making it difficult for tumor cells to enter the liver and not conducive to the growth of metastatic tumors. (2) Cirrhosis is associated with increased intrahepatic portal vein resistance. With the increasing severity of the degree of cirrhosis, peak venous velocity decreases significantly, and the contraflow of a portal venous system can be observed in patients with cirrhosis. The disruption of venous shunting on portal blood flow may prevent tumor cells from reaching the liver (45). (3) Kupffer cells in cirrhosis are activated to release pro-inflammatory factors, which further promote fibrogenesis and up-regulate the expression of Fas receptor, thus increasing the sensitivity of metastatic tumor cells to FasR-mediated apoptosis, and then being killed and cleared by tumor-infiltrating cells (46). (4) With the aggravation of fibrosis, the expression of metalloproteinase tissue inhibitors and other anti-proteinases in hepatocytes and hepatic stellate cells increased, making the overall activity of metalloproteinase reduced, providing resistance to metastasis (47).

Currently, there is no effective non-invasive scoring model to predict the occurrence of NAFL and NASH. In the progression of NAFLD, NAFL plays a dominant role. About 10–20% of NAFL can develop into NASH, and about 20% of NASH can further progress into hepatic fibrosis. Progressive accumulation of fibrosis can lead to the formation of advanced fibrosis or cirrhosis (1114). Consequently, the proportion of advanced fibrosis and cirrhosis in NAFLD is very low. Based on our findings, we hypothesize that simple steatosis and steatohepatitis might be associated with an increased risk of synchronous liver metastasis. Inflammation renders a key role in tumor progression and contributes to the formation of metastases (48). Obesity is closely related to the development of NAFLD, which is a state of chronic low-grade inflammation. Many cytokines and growth factors synthesized and released by adipocytes have direct carcinogenic effects in the gastrointestinal tract, elevated levels of free fatty acids can induce inflammatory pathways (49). The increase of insulin-like growth factor-I (IGF-I) level can enhance the proliferation and anti-apoptosis effect of tumor cells. IGF-I deficient mice developed smaller primary tumors and reduced the burden of liver metastasis (50). Adiponectin and leptin generally affect cellular behavior in opposite ways. Hypoadiponectinemia is a characteristic feature of NASH and is not linked with insulin resistance (51). Overexpression of adiponectin leads to reduced metastasis, adiponectin can inhibit tumor angiogenesis and down-regulate invasive signaling pathways. Leptin expression elevated significantly with the increase in obesity. Leptin receptors are overexpressed in CRC, and leptin has been shown to reduce the apoptosis of CRC cells and promote tumor growth and metastasis (52). Tumor-related macrophages (TAMs) are a major part of the tumor microenvironment that promotes tumor growth, and NOD-like receptor C4 (NLRC4) is a component of the inflammatory body complex. A recent study further confirmed that a high-fat diet-induced NAFLD promotes the growth of CRLM. More importantly, this study explored the role of NLRC4 and interleukin (IL)-1β in the growth of metastatic liver tumors under NAFLD-related liver microenvironment, emphasizing the effect of inflammation changes in NAFLD on the metastatic microenvironment. The study concluded that in NAFLD, NLRC4 contributes to M2 TAMs polarization, IL-1β and vascular endothelial growth factor (VEGF) production in TAMs, which promotes the growth of metastatic liver tumors in CRC (40).

To be obvious, the present study has several limitations. First, this was a single-center study, and there may be some inevitable selection biases. At the same time, this study was retrospective, so only the association between NAFLD and liver metastasis could be determined, rather than the cause-effect relationship. In addition, the sample size of this study was limited and needs to be verified by prospective studies with a large sample size in multiple centers. Second, all non-invasive liver fibrosis scoring models involved in this study were based on the assessment of common clinical parameters and functional alternations in the liver, these alternations may not accurately reflect the replacement of extracellular matrix or changes in fibroblastic cells and cannot fully assess change in the liver microenvironment. Pathological findings can provide strong evidence of immunology and histology. Third, the effect of simple steatosis and steatohepatitis on the prevalence of synchronous liver metastasis was calculated by the “sum of local effects” method, without any actual data or model support.

Despite these limitations, this study keeps the center of attention on the synchronous liver metastasis which can more accurately reflect the role of the “soil” factor in the mechanism and formation of metastasis, also re-recognize the impact of NAFLD on the microenvironment of metastasis. As a matter of fact, this is the first study to explore the specific effects of NAFLD-related fibrosis/cirrhosis on the development of synCRLM. The results of this clinical research provide evidence for the effect of NAFLD on synchronous liver metastasis in CRC patients and potential explanations for different stages in NAFLD. Future studies need to further verify the effects of simple steatosis and steatohepatitis on synchronous liver metastasis, and more prospective clinical trials and animal model experiments are needed to verify and explore the potential mechanism of the interrelation with liver metastasis. As the incidence of NAFLD is increasing rapidly worldwide, follow-up clinical trials of CRC patients should pay more attention to the prevalence of NAFLD in different populations. The interpretation of whether patients with NAFLD and the extent of liver lesions can help to determine the risk of synCRLM and provide guidance for subsequent treatment.

Conclusion

NAFLD may be associated with increased synchronous liver metastasis in CRC patients, while NAFLD progresses to advanced fibrosis or cirrhosis may be associated with decreased synCRLM. This study provides a sufficient and referable clinical basis for exploring the role of NAFLD-related liver microenvironment changes in tumor metastasis.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Author Contributions

YL completed data collection, data analysis, and manuscript drafting. HZ directed the arrangement of the study and supervised the whole writing of the manuscript and manuscript revision. Both the authors conceived the review and approved of the final analysis and results.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.00251/full#supplementary-material

References

1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2018) 68:394–424. doi: 10.3322/caac.21492

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Arnold M, Sierra MS, Laversanne M, Soerjomayaram I, Jemal A, Bray F, et al. Global patterns and trends in colorectal cancer incidence and mortality. Gut. (2017) 66:683–91. doi: 10.1136/gutjnl-2015-310912

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Valderrama-Trevino AI, Barrera-Mera B, Ceballos-Villalva JC, Montalvo-Jave EE. Hepatic metastasis from colorectal cancer. Euroasian J Hepatogastroenterol. (2017) 7:166–75. doi: 10.5005/jp-journals-10018-1241

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Eichler K, Dufas T, Hammerstingl R, Gruber-Rouh T, Vogl TJ, Zangos S. Hepatic arterial infusion with irinotecan in patients with liver metastases of colorectal cancer: results of an extended phase I study. Chemotherapy. (2013) 59:66–73. doi: 10.1159/000348579

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Siriwardena AK, Mason JM, Mullamitha S, Hancock HC, Jegatheeswaran S. Management of colorectal cancer presenting with synchronous liver metastases. Nat Rev Clin Oncol. (2014) 11:446–59. doi: 10.1038/nrclinonc.2014.90

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Waisberg J, Ivankovics IG. Liver-first approach of colorectal cancer with synchronous hepatic metastases:a reverse strategy. World J Hepatol. (2015) 7:1444–9. doi: 10.4254/wjh.v7.i11.1444

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Manfredi S, Lepage C, Hatem C, Coatmeur O, Faivre J, Bouvier AM. Epidemiology and management of liver metastases from colorectal cancer. Ann Surg. (2006) 244:254–9. doi: 10.1097/01.sla.0000217629.94941.cf

PubMed Abstract | CrossRef Full Text | Google Scholar

8. van der Geest LG, Lam-Boer J, Koopman M, Verhoef C, Elferink MA, de Wilt JH. Nationwide trends in incidence, treatment and survival of colorectal cancer patients with synchronous metastases. Clin Exp Metastasis. (2015) 32:457–65. doi: 10.1007/s10585-015-9719-0

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Younossi ZM, Koenig AB, Abdelatif D, Fazel Y, Henry L, Wymer M. Global epidemiology of nonalcoholic fatty liver disease-meta-analytic assessment of prevalence, incidence, and outcomes. Hepatology. (2016) 64:73–84. doi: 10.1002/hep.28431

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Iqbal U, Perumpail BJ, Akhtar D, Kim D, Ahmed A. The epidemiology, risk profiling and diagnostic challenges of nonalcoholic fatty liver disease. Medicines. (2019) 6:41. doi: 10.3390/medicines6010041

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Ekstedt M, Nasr P, Kechagias S. Natural history of NAFLD/NASH. Curr Hepatol Rep. (2017) 16:391–7. doi: 10.1007/s11901-017-0378-2

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Tanaka N, Kimura T, Fujimori N, Nagaya T, Komatsu M, Tanaka E. Current status, problems, and perspectives of non-alcoholic fatty liver disease research. World J Gastroenterol. (2019) 25:163–77. doi: 10.3748/wjg.v25.i2.163

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Byrne CD, Targher G. NAFLD: a multisystem disease. J Hepatol. (2015) 62:S47–64. doi: 10.1016/j.jhep.2014.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Mikolasevic I, Milic S, Orlic L, Stimac D, Franjic N, Targher G. Factors associated with significant liver steatosis and fibrosis as assessed by transient elastography in patients with one or more components of the metabolic syndrome. J Diabetes Complicat. (2016) 30:1347–53. doi: 10.1016/j.jdiacomp.2016.05.014

CrossRef Full Text | Google Scholar

15. Mantovani A, Dauriz M, Byrne CD, Lonardo A, Zoppini G, Bonora E, et al. Association between nonalcoholic fatty liver disease and colorectal tumours in asymptomatic adults undergoing screening colonoscopy: a systematic review and meta-analysis. Metabolism. (2018) 87:1–12. doi: 10.1016/j.metabol.2018.06.004

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Chen J, Bian DX, Zang SF, Yang ZX, Tian GY, Luo Y, et al. The association between nonalcoholic fatty liver disease and risk of colorectal adenoma and cancer incident and recurrence: a meta-analysis of observational studies. Expert Rev Gastroenterol Hepatol. (2019) 13:385–95. doi: 10.1080/17474124.2019.1580143

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Cai B, Liao K, Song XQ, Wei WY, Zhuang Y, Zhang S. Patients with chronically diseased livers have lower incidence of colorectal liver metastases: a meta-analysis. PLoS ONE. (2014) 9:e108618. doi: 10.1371/journal.pone.0108618

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Murono K, Kitayama J, Tsuno NH, Nozawa H, Kawai K, Sunami E, et al. Hepatic steatosis is associated with lower incidence of liver metastasis from colorectal cancer. Int J Colorectal Dis. (2013) 28:1065–72. doi: 10.1007/s00384-013-1656-2

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Hayashi S, Masuda H, Shigematsu M. Liver metastasis rare in colorectal cancer patients with fatty liver. Hepatogastroenterology. (1997) 44:1069–75.

PubMed Abstract | Google Scholar

20. Xi ZF, Fan ZP, Qiu DK, Zeng MD. The relationship between fatty liver disease and liver metastases from colorectal cancer. Chinese J Digestion. (2009) 29:157–60. doi: 10.3760/cma:j.issn:0254-1432.2009.03.005

CrossRef Full Text | Google Scholar

21. Kondo T, Okabayashi K, Hasegawa H, Tsuruta M, Shigeta K, Kitagawa Y. The impact of hepatic fibrosis on the incidence of liver metastasis from colorectal cancer. Br J Cancer. (2016) 115:34–9. doi: 10.1038/bjc.2016.155

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Hanady ZZR, Rees M, Welsh FK, Toogood GJ, Prasad KR, John TK, et al. Fatty liver disease as a predictor of local recurrence following resection of colorectal liver. Br J Surg. (2013) 100:820–26. doi: 10.1002/bjs.9057

CrossRef Full Text | Google Scholar

23. Brouquet A, Nordlinger B. Metastatic colorectal cancer outcome and fatty liver disease. Nat Rev Gastroenterol Hepatol. (2013) 10:266–7. doi: 10.1038/nrgastro.2013.53

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Paget S. The distribution of secondary growths in cancer of the breast 1889. Cancer Metastasis Rev. (1889) 18:571–3. doi: 10.1016/S0140-6736(00)49915-0

CrossRef Full Text | Google Scholar

25. Fidler IJ, Poste G. The “seed and soil” hypothesis revisited. Lancet Oncol. (2008) 9:808. doi: 10.1016/S1470-2045(08)70201-8

CrossRef Full Text | Google Scholar

26. Kinner S, Reeder SB, Yokoo T. Quantitative imaging biomarkers of NAFLD. Dig Dis Sci. (2016) 61:1337–47. doi: 10.1007/s10620-016-4037-1

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Bohte AE, van Werven JR, Bipat S, Stoker J. The diagnostic accuracy of US, CT, MRI and -1H-MRS for the evaluation of hepatic steatosis compared with liver biopsy: a meta-analysis. Eur Radiol. (2011) 21:87–97. doi: 10.1007/s00330-010-1905-5

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Fan JG, Wei L, Zhuang H. Guidelines of prevention and treatment for nonalcoholic fatty liver disease: a 2018 update. Chin J Hepatol. (2018) 26:195–203. doi: 10.3760/cma.j.issn.1007-3418.2018.03.008

CrossRef Full Text | Google Scholar

29. Angulo P, Hui JM, Marchesini G, Bugianesi E, George J, Farrell GC, et al. The NAFLD fibrosis score: a noninvasive system that identifies liver fibrosis in patients with NAFLD. Hepatology. (2007) 45:846–54. doi: 10.1002/hep.21496

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Sterling RK, Lissen E, Clumeck N, Sola R, Correa MC, Montaner J, et al. Development of a simple noninvasive index to predict significant fibrosis in patients with HIV/HCV coinfection. Hepatology. (2006) 43:1317–25. doi: 10.1002/hep.21178

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Nallagangula KS, Nagaraj SK, Venkataswamy L, Chandrappa M. Liver fibrosis: a compilation on the biomarkers status and their significance during disease progression. Future Sci OA. (2018) 4:FSO250. doi: 10.4155/fsoa-2017-0083

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Harrison SA, Oliver D, Arnold HL, Gogia S, Neuschwander-Tetri BA. Development and validation of a simple NAFLD clinical scoring system for identifying patients without advanced disease. Gut. (2008) 57:1441–7. doi: 10.1136/gut.2007.146019

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Sun W, Cui H, Li N, Wei Y, Lai S, Yang Y, et al. Comparison of FIB-4 index, NAFLD fibrosis score and BARD score for prediction of advanced fibrosis in adult patients with non-alcoholic fatty liver disease: a meta-analysis study. Hepatol Res. (2016) 46:862–70. doi: 10.1111/hepr.12647

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Chan AK, Siriwardena AK. Improving definition of the term “synchronous liver metastases” from colorectal cancer. Hepatobiliary Pancreat Dis Int. (2016) 15:458–60. doi: 10.1016/s1499-3872(16)60125-x

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Fan JG. Epidemiology of alcoholic and nonalcoholic fatty liver disease in China. J Gastroenterol Hepatol. (2013) 28:11–17. doi: 10.1111/jgh.12036

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Liu M, Wang JH, Zeng J, Cao XT, He Y. Association of NAFLD with diabetes and the impact of BMI changes: a 5-year cohort study based on 18,507 elderly. J Clin Endocrinol Metab. (2017) 102:1309–16. doi: 10.1210/jc.2016-3440

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Noureddin M, Rinella ME. Nonalcoholic fatty liver disease, diabetes, obesity, and hepatocellular carcinoma. Clin Liver Dis. (2015) 19:361–79. doi: 10.1016/j.cld.2015.01.012

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Earl TM, Nicoud IB, Pierce JM, Wright JP, Majoras NE, Rubin JE, et al. Silencing of TLR4 decreases liver tumor burden in a murine model of colorectal metastasis and hepatic steatosis. Ann Surg Oncol. (2009) 16:1043–50. doi: 10.1245/s10434-009-0325-8

PubMed Abstract | CrossRef Full Text | Google Scholar

39. VanSaun MN, Lee IK, Washington MK, Matrisian L, Gorden DL. High fat diet induced hepatic steatosis establishes a permissive microenvironment for colorectal metastases and promotes primary dysplasia in a murine model. Am J Pathol. (2009) 175:355–64. doi: 10.2353/ajpath.2009.080703

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Ohashi K, Wang Z, Yang YM, Billet S, Tu W, Pimienta M, et al. NOD-like receptor C4 inflammasome regulates the growth of colon cancer liver metastasis in NAFLD. Hepatology. (2019) 70:1582–99. doi: 10.1002/hep.30693

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Karube H, Masuda H, Hayashi S, Ishii Y, Nemoto N. Fatty liver suppressed the angiogenesis in liver metastatic lesions. Hepatogastroenterology. (2000) 47:1541–5.

PubMed Abstract | Google Scholar

42. Hagstrom H, Nasr P, Ekstedt M, Stal P, Hultcrantz R, Kechagias S. Accuracy of noninvasive scoring systems in assessing risk of death and liver-related endpoints in patients with nonalcoholic fatty liver disease. Clin Gastroenterol Hepatol. (2019) 17:1148–56.e4. doi: 10.1016/j.cgh.2018.11.030

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Dahl E, Rumessen J, Gluud LL. Systematic review with meta-analyses of studies on the association between cirrhosis and liver metastases. Hepatol Res. (2011) 41:618–25. doi: 10.1111/j.1872-034X.2011.00813.x

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Gervaz P, Pak-art R, Nivatvongs S, Wolff B, Larson D, Ringel S. Colorectal adenocarcinoma in cirrhotic patients. J Am Coll Surg. (2003) 196:874–9. doi: 10.1016/S1072-7515(03)00117-0

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Mittal P, Gupta R, Mittal G, Kalia V. Association between portal vein color doppler findings and the severity of disease in cirrhotic patients with portal hypertension. Iran J Radiol. (2011) 8:211–7. doi: 10.5812/iranjradiol.4489

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Song E, Chen J, Ouyang N, Wang M, Exton MS, Heemann U. Kupffer cells of cirrhotic rat livers sensitize colon cancer cells to Fas-mediated apoptosis. Br J Cancer. (2001) 84:1265–71. doi: 10.1054/bjoc.2000.1737

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Consolo M, Amoroso A, Spandidos DA, Mazzarino MC. Matrix metalloproteinases and their inhibitors as markers of inflammation and fibrosis in chronic liver disease (Review). Int J Mol Med. (2009) 24:143–52. doi: 10.3892/ijmm_00000217

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Mantovani A, Allavena P, Sica A, Balkwill F. Cancer-related inflammation. Nature. (2008) 454:436–44. doi: 10.1038/nature07205

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Pathak S, Pandanaboyana S, Daniels I, Smart N, Prasad KR. Obesity and colorectal liver metastases: mechanisms and management. Surg Oncol. (2016) 25:246–51. doi: 10.1016/j.suronc.2016.05.021

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Wu YP, Yakar S, Zhao L, Hennighausen L, LeRoith D. Circulating insulin-like growth factor-I levels regulate colon cancer growth and metastasis. Cancer Res. (2002) 62:1030–35.

PubMed Abstract | Google Scholar

51. Hui JM, Hodge A, Farrell GC, Kench JG, Kriketos A, George J, et al. Beyond insulin resistance in NASH: TNF-alpha or adiponectin? Hepatology. (2004) 40:46–54. doi: 10.1002/hep.20280

CrossRef Full Text | Google Scholar

52. Vansaun MN. Molecular pathways: adiponectin and leptin signaling in cancer. Clin Cancer Res. (2013) 19:1926–32. doi: 10.1158/1078-0432.CCR-12-0930

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: non-alcoholic fatty liver disease, colorectal cancer, synchronous liver metastasis, cirrhosis, fibrosis-4 index, non-alcoholic fatty liver disease fibrosis score

Citation: Lv Y and Zhang H-j (2020) Effect of Non-alcoholic Fatty Liver Disease on the Risk of Synchronous Liver Metastasis: Analysis of 451 Consecutive Patients of Newly Diagnosed Colorectal Cancer. Front. Oncol. 10:251. doi: 10.3389/fonc.2020.00251

Received: 27 June 2019; Accepted: 13 February 2020;
Published: 28 February 2020.

Edited by:

Christina Wu, Emory University, United States

Reviewed by:

Ali Coskun, Izmir Bozyaka Eğitim ve Araştırma Hastanesi, Turkey
Ibrahim Halil Sahin, Emory University, United States

Copyright © 2020 Lv and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Hai-jun Zhang, haijunzhang@seu.edu.cn