Abstract

Background. Our study aimed to investigate the relationship between extracellular matrix 1 (ECM1) gene polymorphism and progression of liver fibrosis in the Chinese population. Methods. A total 656 patients with hepatitis B virus (HBV) infection and 298 healthy individuals of the Chinese Han population were recruited for a retrospective case-control study. Of the disease group, 104 cases had chronic hepatitis B (CHB), 266 had LC, and 286 had hepatocellular carcinoma (HCC). Subjects were frequency-matched according to age and gender. Polymorphisms of the ECM1 gene were examined using the MassARRAY SNP genotyping method. Results. There were no associations between genotype and allele frequencies of ECM1 rs3737240 and rs13294 loci with the risk of CHB and CHB-related HCC. After adjustment for age, sex, smoking status, and drinking habits, the GT genotype was dramatically related to a reduced risk of chronic HBV infection in both non-HCC (OR = 0.68, 95% CI: 0.49-0.94) and total chronic HBV infection patients (OR = 0.75, 95% CI: 0.56-1.00). Haplotype analyses revealed twelve protective haplotypes against total chronic HBV infection and four against non-HCC chronic HBV infection. Conclusion. ECM1 gene polymorphism in rs3834087 and rs3754217 loci is associated with a reduced risk of chronic HBV infection but not with liver fibrosis development and the occurrence of HCC.

1. Introduction

Approximately 240 million people, especially those in developing countries, are infected with hepatitis B virus (HBV), and nearly half of these have chronic liver disease [1, 2]. Liver fibrosis is a pathological consequence of chronic liver damage and extracellular matrix (ECM) protein accumulation, which may further result in liver diseases such as liver cirrhosis (LC), liver failure, and hepatocellular carcinoma (HCC) [35]. Liver fibrosis and its associated liver diseases have been considered the most important causes of HBV-induced death [6]. Although the molecular mechanisms of liver fibrosis have been well documented, more effective drugs for liver fibrosis need to be explored [7, 8].

Regardless of whether a toxic or metabolic pathway ensues, a hepatitis virus infection involves the induction of immune cell infiltration and hepatocyte injury, triggering the trans-differentiation of hepatic stellate cells (HSCs) into collagen-producing myofibroblasts [9, 10], and causing deposition of ECM which leads to liver tissue damage [11, 12]. Human ECM1 gene, located on chromosome 1q21 [13, 14], encodes for a glycoprotein which promotes ECM protein binding and regulates the integrity and barrier function of epidermal ECM in a variety of tissues [13, 15, 16]. A recent study on patients and mouse models has shown an inverse correlation between the ECM1 production in chronic hepatitis liver tissues with transforming growth factor-β activation which is responsible for promoting fibrosis and inducing the activation of static HSCs [17].

Several studies identified genetic variants that are associated with inflammatory bowel disease, such as ulcerative colitis [18, 19] and Crohn’s disease [20]. However, the correlation between ECM1 gene polymorphisms and fibrogenic disease formation has not been reported. In this study, we screened the ECM1 gene polymorphisms of both regulatory and coding regions that could modify ECM1 transcription or amino acid sequence and analyzed their association with LC progression in Chinese Han population. Our study might aid in understanding the effects of ECM1 gene polymorphism on the initiation and progression of liver fibrosis and on the occurrence of HCC.

2. Methods

2.1. Subjects

A total of 656 patients with chronic HBV infection, diagnosed as chronic hepatitis B (CHB) (n = 104), LC (n = 266), and HCC (n = 286) based on the Guideline of Prevention and Treatment for CHB [21], were enrolled from the First Hospital of Jilin University [22]. Patients within each group were frequency-matched based on age and sex. CHB was defined as a persistent or intermittent elevation of serum alanine aminotransferase (ALT) (≥2xUTL) level, with increased HBV-DNA level (≥500 IU/ML) for at least 6 months, and the absence of hepatitis C virus (HCV) and human immunodeficiency virus (HIV) infections. CHB/LC was diagnosed as active necroinflammatory liver disease with/without fibrosis by imageological inspection. The diagnosis of HCC was confirmed with histopathology. The exclusion criteria involved the presence of HAV/HEV infection, autoimmune diseases, alcoholic liver disease, silt hemorrhagic liver disease, autoimmune liver disease, and intra- and extrahepatic bile duct stones. Another 298 liver disease-free individuals were selected as healthy controls, and they were negative for HBeAg, HBsAg, anti-HBc, anti-HBe, and anti-HCV. Demographic information, including drinking habits and smoking status, was collected for analysis. Drinkers were classified as those who consumed alcohol more than once per week for more than six months, while smokers were those who smoke daily for over a year. All subjects provided informed consent, and the research was approved by the ethical committee of The First Hospital of Jilin University (Approval number: 2020-Linshen No. 2020-355). All methods were carried out in accordance with relevant guidelines and regulations.

2.2. SNP Selection

SNPs from the promoter, 5′ untranslated regions (UTRs), exons, and 3′ UTR of the ECM1 gene (minor allele frequency [MAF] >5% in the Northern Chinese Han population database) were screened using GeneView of NCBI. SNP function prediction was carried out using the online tool (https://snpinfo.niehs.nih.gov/). Linkage disequilibrium (LD) analysis was performed for the selected SNPs using the “Linkage Disequilibrium Calculator” (https://asia.ensembl.org/Homo_sapiens/Tools/LD?db=core). One of the complete linkage sites were chosen if R2 = 1. Validated and hot SNPs with reported relationship with the disease on Google Scholar were also selected and verified for MAF >5% in Beijing Chinese Han population. Of the SNPs found in the promoter region (rs3834087 and rs3754217), the latter has been predicted to be located in the transcription factor binding site (TFBS) and may have effects on genetic transcription. SNPs found in the exon-missense region (rs3737240 and rs13294) may affect the expression and function of the gene. The ECM1 gene location and the selected SNPs are exhibited in Figure 1.

2.3. Genotype Analysis

Genomic DNA was isolated from whole blood of subjects. SNP genotypes of rs3834087, rs3754217, rs3737240, and rs13294 were analyzed by using the MassARRAY SNP genotyping method (BioMiao Biological Technology, Beijing, CHN). Polymerase chain reaction (PCR) primers and the amplification program are shown in Table 1.

2.4. Statistical Methods

All data were analyzed using the SPSS 18.0 software (SPSS, Chicago, USA). The differences in clinical data and demographic information among groups were calculated using χ2 test or Wilcoxon rank-sum test. Independent segregation of alleles was tested by the Hardy–Weinberg equilibrium (H-WE) analysis. Haplotype analysis of polymorphism was performed by using the UNPHASED 3.1.4 software. -value, odds ratios (ORs), and 95% confidence intervals (CIs) after adjusting for age, gender, and environmental factors were calculated by the bivariate logistic regression analysis. Genotype distributions and allele frequencies were analyzed by the χ2 test or Fisher’s exact test. -value < 0.05 indicated statistical significance.

3. Results

Demographic information of patients including gender, smoking status, age, and alcohol intake is summarized in Tables 2 and 3. No statistical differences were found in gender (χ2 test) and age (Wilcoxon rank-sum test) among the groups () except between non-HCC and HCC groups ( and , respectively). However, significant differences were observed in smoking status (χ2 test) between HCC and healthy controls (), HCC and LC (), as well as HCC and non-HCC (). The same was observed in drinking habits (χ2 test) between non-HCC and healthy controls () and non-HCC and HCC (). Furthermore, the H-WE test determined that all four SNPs (rs3834087, rs3754217, rs3737240, and rs13294) of the controls were in equilibrium (, , , and , respectively).

3.1. Analysis in Chronic HBV Infection Patients and Healthy Controls

The allele and genotype frequencies of the ECM1 SNPs in HBV infected-patients and the control group are displayed in Table 3. There were no correlations between the allele and genotype frequencies of the ECM1 gene polymorphisms at rs3737240 and rs13294 loci of patients with chronic HBV infection. There was an obvious association between the GAG/- genotype of rs3834087 with reduced risk of chronic HBV infection (OR = 0.65, 95% CI: 0.45-0.95). After adjustment for age, smoking status, sex, and drinking habits, bivariate logistic regression analyses showed that the GAG/- and -/- combined genotypes were markedly correlated with a declined risk of chronic HBV infection (OR = 0.66, 95% CI: 0.45-0.96). The GT genotype of rs3754217, as compared to the wild GG genotype, was closely related to a reduced risk of chronic HBV infection (adjusted OR = 0.75, 95% CI: 0.56-1.00).

3.2. Analysis in HCC Patients and Healthy Controls

The genotype and allele frequencies of the ECM1 SNPs in HCC patients and the control group are presented in Table 4. There were no differences in the frequencies of all alleles and genotypes (rs3834087, rs3754217, rs3737240, and rs13294).

3.3. The Genotype and Allele Frequencies of ECM1 SNPs in Non-HCC Patients

The genotype and allele frequencies of ECM1 SNPs in non-HCC patients and the control group are presented in Table 4. No obvious effects were identified between the genotype and allele frequencies of ECM1 SNPs at rs3737240 and rs13294 loci and in the risk of non-HCC after adjusting for smoking status, drinking habits, sex, and age.

However, the DEL allele of rs3834087 was markedly correlated with a decreased risk of CHB (OR = 0.60, 95% CI: 0.39-0.92). Compared to GAG/GAG, the GAG/- and -/- genotypes of rs3834087 were greatly associated with a decreased risk of CHB (OR = 0.54, 95% CI: 0.35-0.85; and OR = 0.55, 95% CI: 0.35-0.86, respectively). Additionally, with adjustment for age, sex, smoking status, and drinking habits, the GT and GT plus TT genotypes of rs3754217, when compared to the wild GG genotype, correlated with a reduced risk of CHB (OR = 0.68, 95% CI: 0.49-0.94; OR = 0.71, 95% CI: 0.52-0.98, respectively).

3.4. Distribution Difference between HBV Infection Subgroups

The allele and genotype frequencies of the ECM1 SNPs in non-HCC (including CHB and LC) and HCC patients are summarized in Table 5. There were no differences in the frequencies of all alleles and genotypes (rs3834087, rs3754217, rs3737240, and rs13294) between the LC and CHB, LC and HCC, and non-HCC and HCC patients.

3.5. Haplotype Analysis

The haplotype distributions of healthy individuals against those of chronic HBV infection, non-HCC, and HCC patients were analyzed, and positive results are shown in Table 6. Between chronic HBV infection patients and controls, the following haplotypes greatly correlated with a decreased risk of chronic HBV infection when compared to their wild type alleles: G-T and T-C alleles of rs3754217-rs3737240 (OR = 0.54, 95% CI: 0.31-0.96; OR = 0.64, 95% CI: 0.44-0.93, respectively) and G-G-T and Del-T-C alleles of rs3834087-rs3754217-rs3737240 (OR = 0.54, 95% CI: 0.31-0.96; OR = 0.65, 95% CI: 0.45-0.96, respectively).

Between non-HCC chronic HBV infection patients and healthy controls, the following haplotypes were related to a decreased risk of non-HCC chronic HBV infection: the T-C alleles of rs3754217-rs3737240 (OR = 0.51; 95% CI: 0.33-0.81) and Del-T alleles of rs3834087-rs3754217 (OR = 0.61; 95% CI: 0.40-0.95) (Table 7).

Liver stiffness measurement (LSM, kPa) scores collected using FibroScan on CHB patients with the different genotypes of rs3834087 and rs3754217 are shown in Supplemental Table 1. Patients with the GT and the GT plus TT genotypes demonstrated lower LSM scores as than that of wild GG genotype patients (Supplemental Table 1). No significant differences were found in the Child-Pugh scores of LC between the different genotypes (Supplemental Table 2).

4. Discussion

Mutation of ECM1 may lead to various genetic diseases such as lipoid proteinosis and autosomal recessive genodermatosis [23, 24]. ECM1 has also been identified to be involved in the differentiation and function of immune cells. Its important role in fibrosis was only recently discovered, but published results have been controversial. In a research by Fan et al. on liver fibrosis in mouse models, liver damage was shown to reduce the levels of ECM1 production during fibrogenesis, and re-expression of ECM1 prevented liver fibrosis progression [17]. In contrast, another study on heart fibrosis showed that ECM1 led to cardiac fibrosis in myocardial infarction by acting as an intermediary between inflammation and fibrosis [25]. We therefore conducted this large case-control study to evaluate the effects of ECM1 on the progression of liver fibrosis in vivo. To our knowledge, our study is the first to have investigated the relationships between ECM1 SNPs and the development of chronic liver fibrosis; other articles on ECM1 SNPs have only involved extrahepatic diseases [18, 20, 26, 27].

We found that the GAG/- and GAG/- plus -/- genotypes of rs3834087 were correlated with a reduced risk of chronic HBV infection both in non-HCC (OR = 0.54, 95% CI: 0.35-0.85; OR = 0.55, 95% CI: 0.35-0.86, respectively) and total chronic HBV infection patients (OR = 0.65, 95% CI: 0.45-0.95; OR = 0.66, 95% CI: 0.45-0.96, respectively). The Del allele of rs3834087 was also related to a lower risk of non-HCC chronic HBV infection (OR = 0.60, 95% CI: 0.39-0.92). After adjustment for age, sex, smoking status, and drinking habits, the GT genotype of rs3754217 was significantly related to a declined risk of chronic HBV infection both in non-HCC (OR = 0.68, 95% CI: 0.49-0.94) and in total chronic HBV infection patients (OR = 0.75, 95% CI: 0.56-1.00). The GT plus TT genotype of rs3754217 was also related to a decreased risk of non-HCC chronic HBV infection (OR = 0.71, 95% CI: 0.52-0.98). Haplotype analyses showed significant association of several haplotypes with a reduced risk of chronic HBV infection, including the T-C allele of rs3754217-rs3737240, DE-T-G allele of rs3834087-rs3754217-rs13294, and the T-C-G allele of rs3754217-rs3737240-rs13294. The Del allele of rs3834087 and the T allele of rs3754217 may be the influencing factors of chronic HBV infection. The mechanisms of these genotype differences on HBV susceptibility remain unclear and warrant further investigations. The rs3834087 and rs3754217 loci in the promoter region of ECM1 and the rs3754217 locus initially predicted as the TFBS were speculated to affect the functions of other associated genes. Previous studies have found that the rs3737240 and rs13294 loci exhibited strong associations with the occurrence of ulcerative colitis [20]. However, we did not find the role of these two SNPs in chronic hepatitis B occurrence, progression of liver fibrosis, and HCC.

In the analysis of the subgroups, no effects on the progression from chronic hepatitis to cirrhosis and HCC were found in the four SNP variants (rs3834087, rs3754217, rs3737240, and rs13294). However, LSM, which has displayed excellent diagnostic accuracy in the identification of HBV-associated fibrosis and cirrhosis [6, 28], revealed significant differences between the GT, GT plus TT, and wild TT genotypes of rs3754217. Mutant genotypes showed low LSM scores, which reflected the lower degree of fibrosis. Due to the limited samples, further studies with expanded sample number are required for further confirmation.

5. Limitations

The primary limitation of this study is that this case-control study was a hospital-based study. Thus, selection bias may have occurred. In addition, due to technical problems in functional studies, we did not find direct evidence on whether these two polymorphisms (rs3834087 and rs3754217) located in promoter regulate or influence ECM1 expression. Therefore, further functional studies are warranted.

6. Conclusion

In conclusion, our study preliminarily demonstrated that the ECM1 locus may mediate the chronicity of HBV infections in vivo. The mechanisms involved in inflammation and fibrosis should be further explored. It should be noted that identification of new therapeutic targets is still required to promote the development of new antifibrotic drugs and fibrotic biomarkers to improve the management of fibrosis.

Abbreviations

SNPs:Single nucleotide polymorphisms
ECM1:Extracellular matrix protein 1
CHB:Chronic hepatitis B
LC:Liver cirrhosis
HCC:Hepatocellular carcinoma
HBV:Hepatitis B virus
ECM:Extracellular matrix
HSCs:Hepatic stellate cells
HCV:Hepatitis C virus
HIV:Human immunodeficiency virus
UTR:Untranslated regions
LD:Linkage disequilibrium
TFBS:Transcription factor binding site
PCR:Polymerase chain reaction
H-WE:Hardy–Weinberg equilibrium
ORs:Odds ratios
CIs:Confidence intervals
LSM, kPa:Liver stiffness measurement.

Data Availability

The data are available from the corresponding author on reasonable request.

Ethical Approval

This study was approved by the Ethics Committee of The First Hospital of Jilin University (Approval number: 2020-Linshen No. 2020-355). Written informed consent for participation in the study was obtained from all participants. All methods were carried out in accordance with relevant guidelines and regulations.

All patients provided written informed consent for the publication of any associated data and accompanying images.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

XL conceived and designed the study; XTH acquired the SNP data and was a major contributor of manuscript writing; TL analyzed and interpreted the patient data; RZ interpreted the SNP data and performed the statistical analysis; and XL revised the manuscript critically for important intellectual content. All authors read and approved the final manuscript.

Acknowledgments

This work was supported by the Finance Special Project on Medical and Health Talents of the Finance Department of Jilin Province [grant no. JLSCZD2019-082]; the Scientific Research Planning Project of the Education Department of Jilin Province [grant no. JJKH20201104KJ]; The Provincial Joint Fund Bethune Medical Project of the Science and Technology Department of Jilin Province [grant no. 20200201497JC]; and WBE Liver Fibrosis Foundation (No. 2020009). The funders had no role in the study design, data collection, data analysis, interpretation, or writing of the report.

Supplementary Materials

Supplemental Table 1. The liver stiffness measurement (LSM, kPa) score’s distribution in FibroScan of rs3834087 and rs3754217 polymorphisms in chronic hepatitis B (CHB) group patients. Supplemental Table 2. The Child-Pugh scores and level’s distribution of rs3834087 and rs3754217 polymorphisms in liver cirrhosis. (Supplementary Materials)