- 1Laboratory of Gastroenterology and Hepatology, State Key Laboratory of Biotherapy, West China Hospital, Sichuan University, Chengdu, China
- 2Department of Gastroenterology, West China Hospital, Sichuan University, Chengdu, China
Liver fibrosis results from the imbalance between extracellular matrix (ECM) production and degradation, which is a common pathological consequence of various chronic liver diseases. Although many miRNAs have been reported in liver fibrosis progression, miRNA-mRNA interactions in its reversal process remain to be elucidated. In the current study, we performed an integrated analysis of miRNA and mRNA expression profiles in the mouse model with the spontaneous reversal potency of liver fibrosis. A total of 102 miRNA and 2,845 mRNAs showed significant differential expression in reversal mice compared to fibrotic mice. Moreover, 3,769 putative negatively correlated miRNA-mRNA pairs were revealed to be potentially implicated in the biological function regulation of small molecule metabolism and ECM organization. By integrating miRNA-mRNA regulatory networks, mmu-miR-1843a-5p, mmu-miR-193a-5p, mmu-miR-194-2-3p, and mmu-miR-30c-2-3p were identified as lysyl oxidases-specific miRNAs that were correlated with fibrosis reversal. Our results provided potential candidate targets for the treatment of liver fibrosis.
Introduction
Liver fibrosis, characterized by the excessive deposition of extracellular matrix (ECM), is a common pathological consequence of various chronic liver diseases, including viral hepatitis, alcoholic and metabolic associated liver diseases (Tsochatzis et al., 2014). Essentially, liver fibrosis is an abnormal wound-healing process response to chronic liver injuries due to various etiological factors (Kisseleva and Brenner, 2021). It was ever believed that liver fibrosis was not reversible. Nevertheless, several studies with sequential liver biopsies have already shown that hepatic fibrosis could be reversed in patients with hepatitis B, hepatitis C, non-alcoholic steatohepatitis (NASH), autoimmune hepatitis and secondary biliary fibrosis after the underlying causes of liver injury are removed (Hammel et al., 2001; Kweon et al., 2001; Arthur, 2002; Czaja and Carpenter, 2004; Dixon et al., 2004). However, the biological molecular mechanism regarding the spontaneous reversal of liver fibrosis is unclear.
MiRNAs are a class of evolutionarily conserved small (∼23 nucleotides in length), single-stranded non-coding RNA generated from the endogenous transcripts following a series of cleavages and modifications (Bartel, 2009). Mature miRNAs act as post-transcriptional regulators via the degradation or translational repression of the target mRNAs by binding to the complementary sequences. It has been estimated that over 30% of human genes are regulated by miRNAs, and one single miRNA may modulate hundreds of target mRNAs (Lewis et al., 2005). Recently, a competitive endogenous RNA hypothesis has been proposed that long non-coding RNAs and circular RNAs may serve as “miRNA sponges” by preferentially occupying the miRNA response elements to regulate mRNAs expression (Thomson and Dinger, 2016), highlighting the central position of miRNAs in cellular function. Accumulated evidence has demonstrated that the dysregulated miRNAs play a crucial role in the development of liver fibrosis by interfering with genes expression and downstream signaling pathways (Wang et al., 2021). Serous miR-122 may serve as a biomarker of inflammation and fibrotic severity in patients with alcoholic liver disease (Satishchandran et al., 2018) and hepatitis C (Bihrer et al., 2011; Trebicka et al., 2013). miR-155 and miR-378 exert pro-inflammatory effects by positively regulating the NF-κB signaling pathway in alcoholic and NASH induced fibrosis, respectively (Bala et al., 2016; Zhang et al., 2019). Thus far, little is known about miRNAs expression profile and their potential functions in the spontaneous reversal of liver fibrosis.
Ongoing liver fibrosis associated with angiogenesis may progress into irreversible cirrhosis with life-threatening complications even inactivation of primary diseases (Tsuchida et al., 2018). Therefore, it is very important to understand the biological molecular mechanism especial in expression of miRNAs and their target mRNAs related to fibrosis regression during the spontaneous reversal of liver fibrosis, by which it would be benefit for prevention of liver cirrhosis on the basis of etiological treatment.
In the present study, we analyzed the expression profiles of miRNAs and mRNAs in the mouse model with the spontaneous reversal potency of liver fibrosis on cessation of liver injuries after given sufficient time. Furthermore, the candidate miRNA-mRNA interactive networks were also defined preliminarily.
Materials and Methods
Animal Model and Sample Collection
C57BL/6 mice (male, 8-week-old), weighing 20–22 g, were obtained from the Experimental Animal Center of Sichuan University (Chengdu, China). All animals were kept under the 12-h light/dark cycles with free access to food and water. The procedure for animal experiments was approved by the Animal Use and Care Committee of West China Hospital, Sichuan University (No. 2017005A).
Liver fibrosis of mice was induced by intraperitoneal injection of thioacetamide (TAA, 200 mg/kg, Sigma-Aldrich, United States) every 3 days for 8 weeks. Control mice received the same volume of normal saline. Control and TAA-fibrotic mice were sacrificed under anesthesia 72 h after the final TAA injection. The liver fibrosis reversal model of mice was firstly established with the same 8-week TAA injection as above. Then, the animals underwent spontaneous recovery for additional 8 weeks. Liver tissues were fixed with 4% paraformaldehyde (PFA) for histological and immunohistochemical examinations, or immediately snap-frozen in liquid nitrogen and stored at −80°C for further gene and protein analysis.
Histological Examination
Liver samples fixed with 4% PFA were sectioned (5 μm thick) after embedded in paraffin and stained with hematoxylin and eosin (H&E) and Sirius Red. Five images per section (at × 100 magnification) from each mouse were selected and assessed by two experienced pathologists independently blinded to the experimental groups. The severity of liver fibrosis was graded according to the Ishak scoring system (Ishak et al., 1995).
Immunohistochemistry Staining
Paraffin-embedded liver sections (3 μm thick) were deparaffinized in xylene and rehydrated with graded ethanol dilutions. Antigen retrieval was performed at high temperature under high pressure in sodium citrate buffer (10 mM, pH = 6.0) for 20 min. After blocking with H2O2 and 10% goat serum, the sections were incubated with anti-collagen I (1:200, Abcam, United Kingdom) or anti-α-smooth muscle actin (α-SMA, 1:200, Abcam) overnight at 4°C followed by incubation with biotin-streptavidin-horseradish peroxidase (HRP) detection system (ZSGB-BIO, Beijing, China) at room temperature. Finally, the sections were stained with a solution of 3, 3’-diaminobenzidine (DAB, ZSGB-BIO) and counterstained with hematoxylin.
Hepatic Hydroxyproline Measurement
Hepatic hydroxyproline content, a major component of collagen, was measured in 400–500 mg liver samples using the Hydroxyproline Assay Kit (Nanjing Jiancheng Bioengineering Institute, Nanjing, China) according to the manufacturer’s instructions and expressed as microgram per gram tissue (μg/g liver).
Western Blot Analysis
Whole proteins from liver tissues were extracted using a protein extraction kit (Nanjing KeyGen Biotech Co., Ltd., Nanjing, China). 50 μg (7.2 μg/μL, 7 μL) of proteins for each sample were resolved by 10% SDS-PAGE, transferred to PVDF membrane (Merck Millipore, United States), and blocked with 5% non-fat dry milk before incubated with anti-collagen I (1:1,000, Abcam), anti-α-SMA (1:4,000, Abcam) or anti-HSC70 (1:1,000, Santa Cruz Biotechnology, United States) overnight at 4°C. Then, the blots were washed and incubated with HRP-conjugated secondary antibodies (1:20,000, ZSGB-BIO) at room temperature for 2 h. Protein bands were visualized by chemiluminescence using Western Blotting Luminol Reagent (Merck Millipore), and densitometric analyses were made using Quantity One software (v4.6.2). Protein levels were normalized against HSC70 and were shown as fold changes relative to the control group.
RNA Extraction and Quality Control
To obtain miRNA and mRNA expression profiles, total RNA was extracted from liver samples using TRIzol reagent (Thermo Fisher Scientific, United States). Next, the RNA samples were qualified and quantified as follows: RNA degradation and contamination were monitored on 1% agarose gels; RNA purity and concentration were then examined using the NanoPhotometer spectrophotometer (Implen, Germany); RNA integrity and concentration were finally assessed and quantified using the RNA Nano 6,000 Assay Kit of the Bioanalyzer 2,100 system (Agilent Technologies, United States). The data on the quality of RNA were provided in Supplementary Figure 1.
Small RNA Sequencing (Small RNA-Seq) and Data Analysis
Small RNA libraries were constructed using the NEBNext Multiplex Small RNA Library Prep Kit for Illumina (NEB, United States), following the manufacturer’s instructions. Then, the purified libraries were quantified by Qubit 2.0 Fluorometer (Thermo Fisher Scientific), and validated by Agilent 2,100 bioanalyzer (Agilent Technologies). After cluster generation, the final library preparations were submitted for small RNA-seq on the Illumina-HiSeq 2500 platform (Illumina, United States) and 50 bp single-end reads were generated. Raw data were processed through custom perl and python scripts, and were then filtered by Q20, Q30, and GC contents. The acquired clean data were mapped to the mouse reference genome (mm10) and aligned against miRbase (Release 22.1) by Bowtie (v1.3.0) (Langmead et al., 2009) to identify mature miRNAs. For each sample, the miRNA expression profiles were normalized by transcript per million (TPM), as previously reported (Zhou et al., 2010). Differentially expressed miRNAs (DEMs) were analyzed using DESeq2 (v1.16.1), with a significance threshold of false discovery rate (FDR) < 0.05 and |fold change| ≥ 2.
mRNA Sequencing (RNA-Seq) and Data Analysis
For RNA-seq, 1 μg of total RNA per sample was used for library preparation. The NEBNext Ultra RNA Library Prep Kit for Illumina (NEB) was used for library construction, according to the manufacturer’s protocols. After purification and quantification, the libraries were subjected to 150bp paired-end sequencing on the Illumina NovaSeq platform (Illumina). As described above, quality-filtered reads were then mapped to the mouse reference genome (mm10) using HISAT2 (v2.2.1) (Kim et al., 2019). Read counts for each gene were obtained with featureCounts (v1.5.0-p3) and were standardized as fragments per kilobase million (FPKM). DESeq2 (v1.16.1) was used to analyze differentially expressed genes (DEGs) with a filter criterion of FDR < 0.05 and |fold change| ≥ 2. Heatmaps and volcano plots were generated to visualize gene expression patterns using TB (Toolbox for Biologists) tools (v1.082) (Chen C. et al., 2020) and GraphPad Prism 8.4.0 (GraphPad, United States), respectively. Small RNA-seq and RNA-seq were performed by Novogene (Beijing, China).
miRNA-mRNA Integrated Analysis
The target genes of DEMs were predicted by an online tool miRWalk1 (Sticht et al., 2018), which collates data from multiple prediction programs (TargetScan, miRDB and miRTarBase). The intersections between the predicted target genes and DEGs were screened out for further Pearson correlation analysis. Since miRNA mainly inhibits the translational activity of its target mRNA, negatively paired miRNA-mRNA correlations were selected, according to the screening criteria of a p-value < 0.05 and a Pearson correlation coefficient (PCC) ≥ 0.8. DEGs involved in the negative miRNA-mRNA pairs were subjected to subsequent analysis. Then, the miRNA-mRNA interaction networks were constructed and visualized using Cytoscape software (v3.8.2).
Functional Enrichment Analysis
To explore the functional roles of these significant DEGs, GO enrichment and KEGG pathway analysis were performed using Metascape2 (Zhou et al., 2019). The GO terms and KEGG pathways with FDR less than 0.05 were considered significantly enriched, and visualized by online bioinformatics platforms.3 Additionally, the enriched Reactome pathway interaction networks were constructed using Cytoscape plug-in ClueGO (v2.5.7).
Protein-Protein Interaction (PPI) Network
The differentially expressed target genes were mapped to the STRING database4 (Szklarczyk et al., 2019) to screen PPI, and results were visualized by Cytoscape software as a network structure.
Quantitative Real-Time PCR (qRT-PCR)
Total RNA was extracted from the frozen liver tissues using an RNA Isolation Kit (Foregene, Chengdu, China), according to the manufacturer’s instructions. Then, equal amounts of 1 μg RNA were reverse-transcribed into cDNA using RevertAidTM First-Strand cDNA Synthesis Kit (Thermo Fisher Scientific) or miRcute Plus miRNA First-Strand cDNA Synthesis Kit (Tiangen, Beijing, China), respectively. qRT-PCR was performed in triplicate using SYBR Green qPCR Master Mix (Bimake, United States) on the CFX96 Real-Time PCR Detection System (Bio-Rad). Expression of miRNA and mRNA were normalized to U6 snRNA or Gapdh, and were shown as fold changes relative to the control group. Primer sequences were listed in Supplementary Table 1.
Statistical Analysis
All data were expressed as mean ± standard deviation and analyzed by SPSS 19.0 software (SPSS, United States). The Student’s t-test was performed for comparisons between two groups, and one-way ANOVA (Students-Newman-Keuls test) was used for comparisons of multiple groups. A p-value < 0.05 was considered statistically significant.
Results
Model Establishment of the Spontaneous Reversed Liver Fibrosis
Following the animal treatment procedure shown in Figure 1A, we successfully established a mouse model of spontaneous reversed liver fibrosis. At the end of 8-week modeling, liver tissues in the TAA group showed a typical advanced fibrosis appearance with disordered hepatic lobular structure and continuous fibrotic septa (Figure 1B). However, the pathological characteristics of liver fibrogenesis and liver injury were remarkably regressed after TAA withdrawal for additional 8 weeks (Figure 1B). In support of this, Ishak score and hepatic hydroxyproline content in the recovery group were 43.6% and 25.2% lower than those in the TAA group, respectively, p < 0.05 (Figures 1C,D). Additionally, the protein levels of collagen I and α-SMA were also significantly decreased in the recovery group compared with the TAA group (Figures 1B,E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/706341/fgene-12-706341-HTML-r1/image_m/fgene-12-706341-g001.jpg)
Figure 1. Model establishment of the spontaneous reversed liver fibrosis. (A) Flow diagram of the animal treatment procedure for establishing the spontaneous reversed liver fibrosis model. (B–D) Representative liver histology of H&E and Sirius Red staining, immunohistochemistry of collagen I and α-SMA (B), and quantification of hepatic fibrosis as measured by the Ishak score (C) and hepatic hydroxyproline content (D). (E) Expression of collagen I and α-SMA were determined by Western blot. Data were presented as mean ± SD. n = 6/group, *p < 0.05. α-SMA, α-smooth muscle actin; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; H&E, hematoxylin and eosin; TAA, thioacetamide.
Identification of DEMs and DEGs
In the miRNA and mRNA expression profiling by deep sequencing, 9 mice were selected from the control (n = 3), TAA (n = 3) and recovery group (n = 3). DEMs and DEGs were visualized in heatmaps (Figures 2A, 3A) and volcano plots (Figures 2B, 3B). Information on the top 20 DEMs and top 20 DEGs in detail was listed in Tables 1, 2, respectively. According to the screening criteria of FDR < 0.05 and |fold change| ≥ 2, 206 miRNAs were differentially expressed in the TAA group compared with the control group, among which 102 DEMs were reversed in the recovery group (Supplementary Table 2). Furthermore, 53 miRNAs were up-regulated in TAA-fibrotic mice while their expression was suppressed in mice with reversed liver fibrosis. There were 49 miRNAs down-regulated in TAA-fibrotic mice while rescued in mice with reversed liver fibrosis (Figure 2C). Among these DEMs, mmu-miR-592-5p was overexpressed with the highest fold change of 10.67. In contrast, mmu-miR-344d-3p was downregulated with the highest fold change of 18.96.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/706341/fgene-12-706341-HTML-r1/image_m/fgene-12-706341-g002.jpg)
Figure 2. Identification of differentially expressed miRNAs (DEMs). (A) Heatmap of miRNAs via hierarchical cluster analysis. Colors from green to red represented the miRNA expression abundance from poor to rich. (B) Volcano plots of miRNAs between TAA and Control, and Recovery and TAA, respectively. Red dots represented up-regulated miRNAs, green dots represented down-regulated miRNAs, and gray dots represented miRNAs with no significant difference. (C) Venn diagram of DEMs between TAA-fibrotic mice and reversal mice. DEMs, differentially expressed miRNAs; FDR, false discovery rate; TAA, thioacetamide.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/706341/fgene-12-706341-HTML-r1/image_m/fgene-12-706341-g003.jpg)
Figure 3. Identification of differentially expressed genes (DEGs). (A) Heatmap of mRNA via hierarchical cluster analysis. Colors from green to red represented the mRNA expression abundance from poor to rich. (B) Volcano plots of mRNAs between TAA and Control, and Recovery and TAA, respectively. Red dots represented up-regulated mRNAs, green dots represented down-regulated mRNAs, and gray dots represented mRNAs with no significant difference. (C) Venn diagram of DEGs between TAA-fibrotic mice and reversal mice. DEGs, differentially expressed genes; FDR, false discovery rate; TAA, thioacetamide.
Using the same criteria as for miRNAs, a total of 4,118 DEGs were identified in liver fibrosis. Notably, 2,100 out of 3,025 up-regulated mRNA transcripts upon TAA stimulation were suppressed on cessation of TAA exposure, while 745 out of 1,093 down-regulated mRNA transcripts were rescued simultaneously followed by the withdrawal of repeated liver damages (Figure 3C and Supplementary Table 3). Within the identified 2,845 DEGs, Cyp3a16, Rdh16f1, Sprr1a, Acta1, Myh2, etc. were decreased, while Mup19, Mup15, Mup8, Mup1, Meig1, etc. were increased. Overall, miRNA and mRNA with similar expression patterns in control and reversal mice and opposite to fibrotic mice may potentially relate to the progression and regression of liver fibrosis.
Target Gene Identification and Functional Enrichment
One hundred and two DEMs and 2,845 DEGs identified in the previous steps were considered for the miRNA-mRNA integrated analysis. A total of 16,706 genes were consistently predicted as potential targets of 102 DEMs using an online bioinformatics database miRWalk. Then a Venn diagram showed that 2,492 overlapped genes were obtained between the predicted target genes and the identified DEGs (Figure 4A), which were selected for further correlation analysis. According to the standard of absolute PCC ≥ 0.8, we finally identified 11,434 differentially expressed miRNA-mRNA pairs, including 3,769 negative correlating pairs and 7,665 positive correlating pairs (Figure 4B). Since miRNAs generally suppress the expression of target mRNAs, those pairs where the miRNA and mRNA expression levels changed in the same direction were filtered out. In the case of 3,769 negative miRNA-mRNA interaction pairs, 2,168 pairs that composed 23 miRNA and 1,056 mRNA were up-regulated miRNAs vs. down-regulated mRNAs, while 1,601 pairs composed 53 miRNA and 434 mRNA were down-regulated miRNAs vs. up-regulated mRNAs (Supplementary Table 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/706341/fgene-12-706341-HTML-r1/image_m/fgene-12-706341-g004.jpg)
Figure 4. Target gene identification and functional enrichment. (A) Venn diagram between identified DEGs and predicted target genes of DEMs. (B) Pearson correlation analysis to identify the negative correlating miRNA-mRNA pairs. (C) Top 10 enriched GO terms of the corresponding biological process, cellular component and molecular function for the up- and down-regulated DEGs, respectively. (D) Top 20 enriched KEGG pathways for the up- and down-regulated DEGs, respectively. DEGs, differentially expressed genes; DEMs, differentially expressed miRNAs; FDR, false discovery rate; PCC, Pearson correlation coefficient.
To understand the biological function of these dysregulated target genes, GO enrichment and KEGG pathway analysis were performed for the 1,490 DEGs. GO analysis revealed that the up-regulated DEGs were enriched in the metabolic process of lipid and amino acids, peroxisome and membrane components, and protein binding and enzyme activity, while the down-regulated DEGs were enriched in terms related to ECM organization, such as cell-cell adhesion, response to growth factor, and polymeric cytoskeletal fiber (Supplementary Table 5). Ten most significantly enriched GO terms, involving biological processes, cellular component and molecular function, were shown in Figure 4C. The correlating transcripts were also analyzed for enrichment in KEGG pathway categories (Supplementary Table 6), and the top 20 KEGG pathways were identified for the up- or down-regulated DEGs (Figure 4D). Consistent with GO analysis, lipid and amino acid metabolism were significantly enriched by KEGG pathway analysis for the up-regulated genes (Figure 4D, left panel). In contrast, the top 20 down-regulated genes were mainly enriched in regulation of actin cytoskeleton, cell cycle, ECM-receptor interaction, protein digestion and absorption, and several signaling pathways, including the Hedgehog signaling pathway, Rap1 signaling pathway, and PI3K-Akt signaling pathway (Figure 4D, right panel).
Regulatory Networks of Putative miRNA-mRNA Interaction
To further reveal the regulatory role of DEMs and DEGs in the spontaneous reversal process of liver fibrosis, we constructed miRNA-mRNA interaction networks based on the Reactome pathway database. Indeed, most up-regulated DEGs were enriched in the metabolic pathways, mainly focusing on biological oxidations, lipid metabolism, and amino acid metabolism (Figure 5A). By integrating these biological pathways for both miRNA and mRNA categories (Figure 5B), we found that multiple miRNAs and their target mRNAs were involved in respective metabolic regulation, totally covering more than 90% of the identified DEMs (50/53). Remarkably, distinct pathways might interact with each other through the hub genes. Lipid metabolism and biological oxidations might interact with each other via Cytochrome P450 genes (Cyp2e1, Cyp8b1, Cyp3a13, etc.), while biological oxidations might interact with amino acid metabolism through Gm4737, Mat1a, Gstz1, and Ahcy. Unanimously, Acat1 might serve as a linker between lipid metabolism and amino acid metabolism.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/706341/fgene-12-706341-HTML-r1/image_m/fgene-12-706341-g005.jpg)
Figure 5. Regulatory networks of down-regulated miRNAs and up-regulated mRNAs interaction. (A) Enriched Reactome pathway interaction networks for up-regulated DEGs. Each node represented a Reactome pathway, and the colors indicated functional clusters of related biological pathways. The node size correlated with the enrichment significance, while the edge thickness reflected the strength of correlation between individual pathways. (B) Schematic overview of the miRNA-mRNA interaction networks in specific pathways. The triangles represented miRNAs and the circles represented the target mRNAs. The sizes of triangles and circles represented the fold changes in reversal mice compared to TAA mice, while the edge thickness represented Pearson correlation coefficient between miRNA and mRNA. DEGs, differentially expressed genes; TAA, thioacetamide.
Additionally, 1,056 down-regulated DEGs were mapped preferentially to pathways associated with the extracellular matrix, cell proliferation, biosynthesis, and kinase signaling pathway (Figure 6A). Since liver fibrosis reversal is a dynamic process involving coordinated changes in ECM synthesis and degradation, pathways related to ECM organization (collagen formation, elastic fiber formation and ECM degradation) were selected to validate the miRNA-mRNA regulatory networks (Figure 6B). 14 out of 23 DEMs were directly involved in ECM degradation, of which 12 DEMs regulated collagen formation, and 8 DEMs regulated elastic fiber formation. Additionally, the genes encoding collagen (Col1a1, Col7a1, Col8a1, Col12a1, etc.) and elastic fiber (Eln) were involved in the ECM degradation by regulation of collagen and elastic fiber formation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/706341/fgene-12-706341-HTML-r1/image_m/fgene-12-706341-g006.jpg)
Figure 6. Regulatory networks of up-regulated miRNAs and down-regulated mRNAs interaction. (A) Enriched Reactome pathway interaction networks for down-regulated DEGs. Each node represented a Reactome pathway, and the colors indicated functional clusters of related biological pathways. The node size correlated with the enrichment significance, while the edge thickness reflected the strength of correlation between individual pathways. (B) Schematic overview of the miRNA-mRNA interaction networks in specific pathways. The triangles represented miRNAs and the circles represented the target mRNAs. The sizes of triangles and circles represented the fold changes in reversal mice compared to TAA mice, while the edge thickness represented Pearson correlation coefficient between miRNA and mRNA. DEGs, differentially expressed genes; TAA, thioacetamide.
Key miRNA-mRNA Interaction Networks in Fibrosis Reversal
Collagen and elastin are major components of fibrotic tissues in the liver. Therefore, the DEGs intersections with collagen formation or elastin formation were assessed. 4 hub genes that belong to the lysyl oxidase (LOX) family (Lox, Loxl1, Loxl2, Loxl3, Figure 7A) were identified. Using the STRING database, PPI network analysis was performed for all down-regulated mRNAs, and a sub-network was further constructed by highlighting the first-degree neighbors of these key hub genes obtained from the PPI analysis (Figure 7B). Interestingly, we found that genes encoding ECM proteins and activated hepatic stellate cell (HSC) markers were closely clustered within the PPI network, suggesting their coordinate regulation on liver fibrosis reversal. To further clarify the miRNA-mRNA regulatory interactions, 4 overlapped miRNAs (mmu-miR-1843a-5p, mmu-miR-193a-5p, mmu-miR-194-2-3p, and mmu-miR-30c-2-3p) were identified based on a Venn diagram of DEMs targeting ECM component, HSCs activation and LOXs regulation (Figure 7C). These miRNAs and their target genes constructed a regulatory networks for miRNA-mRNA interactions (Figure 7D). A single miRNA might interact with different genes, while the same gene might be a target of different miRNAs. For example, mmu-miR-1843a-5p modulated Lox, Loxl1, and Loxl3, while Loxl1 was the target of mmu-miR-1843a-5p, mmu-miR-194-2-3p, and mmu-miR-30c-2-3p. These data indicated complex interactions of miRNA-mRNA in the regulation of the spontaneous reversal of liver fibrosis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/706341/fgene-12-706341-HTML-r1/image_m/fgene-12-706341-g007.jpg)
Figure 7. Identification and validation of key miRNAs and mRNAs in fibrosis reversal. (A) Venn diagram of DEGs between collagen formation and elastin formation to identify the hub genes. (B) PPI network of the hub genes (Lox, Loxl1, Loxl2, Loxl3) and their first-degree neighbors. The node size represented node degree in the PPI network, while the edge thickness reflected the combined scores between different nodes. (C) Venn diagram of DEMs targeting ECM component, HSCs activation and LOXs regulation. (D) miRNA-mRNA regulatory network of the overlapped miRNAs (mmu-miR-1843a-5p, mmu-miR-193a-5p, mmu-miR-194-2-3p, and mmu-miR-30c-2-3p) and their target genes. (E) Expression of miRNAs and mRNAs were validated by qRT-PCR in mice livers. Data were presented as mean ± SD. n = 6/group, *p < 0.05. DEGs, differentially expressed genes; DEMs, differentially expressed miRNAs; ECM, extracellular matrix; HSCs, hepatic stellate cells; LOXs, lysyl oxidases; Loxl1, lysyl oxidase-like 1; Loxl2, lysyl oxidase-like 2; Loxl3, lysyl oxidase-like 3; PPI, protein-protein interaction; qRT-PCR, quantitative real-time PCR; TAA, thioacetamide.
Validation of Key miRNA and mRNA
To validate the RNA sequencing results, we selected key miRNAs (mmu-miR-1843a-5p, mmu-miR-193a-5p, mmu-miR-194-2-3p, and mmu-miR-30c-2-3p) and mRNAs (Lox, Loxl1, Loxl2, and Loxl3)for qRT-PCR analysis in livers from the TAA group and the recovery group. As shown in Figure 7E, the RNA expression patterns of real-time PCR data were consistent with that of RNA sequencing data. Compared with TAA-fibrotic mice, miRNAs expression was significantly up-regulated while mRNAs expression was down-regulated in reversal mice, p < 0.05 (Figure 7E).
Discussion
As an important pathological feature and an intermediate development stage for most chronic liver diseases, liver fibrosis develops and progresses due to the imbalance between ECM production and degradation (Karsdal et al., 2015). There is growing clinical evidence that early to moderate hepatic fibrosis can regress and possibly even resolve in a number of liver diseases based on histologic assessment (Ellis and Mann, 2012; Kisseleva and Brenner, 2021). In recent years, high-throughput sequencing technology enables a deep and efficient exploration at the molecular level for the regulatory mechanisms in fibrosis reversal (Stark et al., 2019). This study performed an integrative analysis of miRNA and mRNA expression profiles in the mouse model induced by TAA. A series of putative miRNAs and their target mRNAs were screened out. Their interactions and potential functions in the spontaneous reversal of liver fibrosis were revealed.
It has been known that miRNAs are epigenetic modulators that inhibit the expression of target genes. Although a lot of miRNAs have been identified in liver fibrosis progression, there is less report on miRNA-mRNA interactions and their regulatory mechanisms involved in fibrosis reversal. 102 differentially expressed miRNAs (DEMs), with 49 significantly up-regulated and 53 significantly down-regulated were identified in this study. Among them, miR-200a-5p and miR-122-3p have been reported as potential biomarkers associated with fibrosis progression (Van Keuren-Jensen et al., 2016). Interestingly, the miRNAs that regulate HSCs phenotype were consistently reversed with a concomitant decrease in hepatic ECM deposition. miRNAs that activate HSCs (miR-199a-3p and miR-214-3p) were significantly suppressed (Ma et al., 2018; Yang et al., 2020), while those restrain the proliferation and activation of HSCs (miR-378a-3p, miR-139-5p, miR-455-3p, and miR-193a-3p) were remarkably increased during the fibrosis regression (Hyun et al., 2016; Ju et al., 2019; Wei et al., 2019; He et al., 2021). Among the differentially expressed genes (DEGs), we obtained 2,845 differential genes between reversal mice and TAA-fibrotic mice. According to the integrative analysis, 3,769 putative negatively correlated miRNA-mRNA pairs containing 76 DEMs and 1,490 DEGs were revealed to be potentially implicated in the pathogenesis of this spontaneous reversal process. Nevertheless, 26 out of 102 DEMs failed to regulate the target genes, and only 52.4% of the DEGs were regulated by the corresponding miRNAs. These results suggested that a single miRNA could target multiple mRNAs, while multiple miRNAs could target a single mRNA. However, not all target genes had contrasting expression patterns compared to miRNA, which might be due to the synergistic or antagonistic regulatory effects of individual miRNAs with concordant or opposed function on the same intronic targets (Fabian et al., 2010).
GO enrichment and pathway analysis indicated that small-molecule metabolic pathways related to biological oxidation were strongly up-regulated during the reversal process of fibrosis. Substantial evidence has shown that lipotoxicity, originated from the imbalance between lipid uptake and utilization, is a key contributor to the development of fatty liver diseases by inducing endoplasmic reticulum stress, mitochondrial dysfunction and oxidative stress (Marra and vegliati-Baroni, 2018). As a canonical pathway involved in lipid metabolism, PPAR activation exerts potent anti-inflammatory, anti-steatosis, and anti-fibrosis effects due to its extended protection for hepatocytes, liver sinusoidal endothelial cells (LSECs), and HSCs via multi-pathways (Han et al., 2020). Several clinical trials of NASH indicated that PPAR agonists may improve steatosis, inflammation and fibrosis (Ratziu et al., 2016; Gawrieh et al., 2021). Our results provided valuable evidence that the PPAR pathway took part in the fibrosis regression regulation, and suggested that stimulation of PPAR pathway might serve as a potential therapeutic strategy for human liver fibrosis. Besides, the down-regulated DEGs were mainly enriched in signal transduction pathways related to ECM organization. Among these, Hedgehog and PI3K-Akt are two important stress-responsive pathways in the liver. During liver wound healing, both pathways are activated, ultimately leading to the transdifferentiation of HSCs into myofibroblasts, LSECs capillarization, and M2-polarization of macrophages (Yu et al., 2015; Machado and Diehl, 2018; Wang et al., 2019). Our results validated that deregulation of Hedgehog and PI3K-Akt signaling pathways was consistent with the ameliorative pathological characteristics of liver fibrogenesis in reversal mice. Although the small molecule inhibitors of Hedgehog and PI3K-Akt have been developed for use as anticancer agents, it is worthwhile to investigate their anti-fibrotic efficacy in future clinical trials.
Through the integration of the deregulated signaling pathways, we identified a series of genes in LOX family that affect collagen and elastin formation. Interestingly, all of these genes were significantly associated with those encoding ECM structures and activated HSCs markers in the PPI network. Supportively, LOXs-mediated cross-linking of both collagen and elastin throughout the initiation, progression, and regression of liver fibrosis (Chen W. et al., 2020). Over-expression of LOXs leads to an abnormal ECM synthesis and subsequent cross-linking, forming a vicious circle. The resultant ECM stabilization contributes to fibrosis progression and retards fibrolysis, for example, through matrix metalloproteinases-induced proteolytic degradation (Ikenaga et al., 2017; Zhao et al., 2018). As LOXs-specific miRNAs, mmu-miR-1843a-5p, mmu-miR-193a-5p, mmu-miR-194-2-3p, and mmu-miR-30c-2-3p identified in this study might be the potential pharmacological targets for the reversal of liver fibrosis.
In conclusion, the expression profiles of both miRNAs and mRNAs in the mouse model with the spontaneous reversal potency of liver fibrosis have been identified. Some differentially expressed miRNAs and their target mRNAs were correlated with fibrosis regression. A core functional miRNA-mRNA regulatory network which was meaningful to understand the molecular mechanisms for liver fibrosis reversal was constructed. The potential candidate targets for the treatment of liver fibrosis provided in this study need further elucidation.
Data Availability Statement
The raw miRNA and mRNA sequencing data have been deposited in the NCBI Gene Expression Omnibus (GSE173961 and GSE173962).
Ethics Statement
The animal study was reviewed and approved by the Animal Use and Care Committee of West China Hospital, Sichuan University.
Author Contributions
CT and JG conceived and supervised the study. YT, CZ, TL, LZ, and YX performed the experiments. YT, HT, and RL analyzed the data. YT, JG, and CT wrote the manuscript with input from all the authors. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (81670551, U1702281, 81873584, 82000613, and 82000574), National Key R&D Program of China (2017YFA0205404), Sichuan Science and Technology Program (2020YJ0084 and 2021YFS0147), and Post-Doctor Research Project, West China Hospital, Sichuan University (2019HXBH074).
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/fgene.2021.706341/full#supplementary-material
Supplementary Figure 1 | Quality control and quantification of RNA samples.
Supplementary Table 1 | Primer sequences for qRT-PCR.
Supplementary Table 2 | Differentially expressed miRNAs (DEMs) in fibrosis reversal.
Supplementary Table 3 | Differentially expressed genes (DEGs) in fibrosis reversal.
Supplementary Table 4 | Putative negative miRNA-mRNA pairs by integrated analysis.
Supplementary Table 5 | GO enrichment analysis of differentially expressed genes (DEGs) in fibrosis reversal.
Supplementary Table 6 | KEGG pathway analysis of differentially expressed genes (DEGs) in fibrosis reversal.
Footnotes
- ^ http://mirwalk.umm.uni-heidelberg.de/
- ^ https://metascape.org/
- ^ https://cloud.oebiotech.cn and http://www.bioinformatics.com.cn
- ^ https://string-db.org/
References
Arthur, M. J. (2002). Reversibility of liver fibrosis and cirrhosis following treatment for hepatitis C. Gastroenterology 122, 1525–1528. doi: 10.1053/gast.2002.33367
Bala, S., Csak, T., Saha, B., Zatsiorsky, J., Kodys, K., Catalano, D., et al. (2016). The pro-inflammatory effects of miR-155 promote liver fibrosis and alcohol-induced steatohepatitis. J. Hepatol. 64, 1378–1387. doi: 10.1016/j.jhep.2016.01.035
Bartel, D. P. (2009). MicroRNAs: target recognition and regulatory functions. Cell 136, 215–233. doi: 10.1016/j.cell.2009.01.002
Bihrer, V., Friedrich-Rust, M., Kronenberger, B., Forestier, N., Haupenthal, J., Shi, Y., et al. (2011). Serum miR-122 as a biomarker of necroinflammation in patients with chronic hepatitis C virus infection. Am. J. Gastroenterol. 106, 1663–1669. doi: 10.1038/ajg.2011.161
Chen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y., et al. (2020). TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 13, 1194–1202. doi: 10.1016/j.molp.2020.06.009
Chen, W., Yang, A., Jia, J., Popov, Y. V., Schuppan, D., and You, H. (2020). Lysyl oxidase (LOX) family members: rationale and their potential as therapeutic targets for liver fibrosis. Hepatology 72, 729–741. doi: 10.1002/hep.31236
Czaja, A. J., and Carpenter, H. A. (2004). Decreased fibrosis during corticosteroid therapy of autoimmune hepatitis. J. Hepatol. 40, 646–652. doi: 10.1016/j.jhep.2004.01.009
Dixon, J. B., Bhathal, P. S., Hughes, N. R., and O’Brien, P. E. (2004). Nonalcoholic fatty liver disease: improvement in liver histological analysis with weight loss. Hepatology 39, 1647–1654. doi: 10.1002/hep.20251
Ellis, E. L., and Mann, D. A. (2012). Clinical evidence for the regression of liver fibrosis. J. Hepatol. 56, 1171–1180. doi: 10.1016/j.jhep.2011.09.024
Fabian, M. R., Sonenberg, N., and Filipowicz, W. (2010). Regulation of mRNA translation and stability by microRNAs. Annu. Rev. Biochem. 79, 351–379. doi: 10.1146/annurev-biochem-060308-103103
Gawrieh, S., Noureddin, M., Loo, N., Mohseni, R., Awasty, V., Cusi, K., et al. (2021). Saroglitazar, a PPAR-alpha/gamma agonist, for treatment of nonalcoholic fatty liver disease: a randomized controlled double-blind phase 2 trial. Hepatology. doi: 10.1002/hep.31843 [Epub ahead of print].
Hammel, P., Couvelard, A., O’Toole, D., Ratouis, A., Sauvanet, A., Flejou, J. F., et al. (2001). Regression of liver fibrosis after biliary drainage in patients with chronic pancreatitis and stenosis of the common bile duct. N. Engl. J. Med. 344, 418–423. doi: 10.1056/NEJM200102083440604
Han, X., Wu, Y., Yang, Q., and Cao, G. (2020). Peroxisome proliferator-activated receptors in the pathogenesis and therapies of liver fibrosis. Pharmacol. Ther. 222:107791. doi: 10.1016/j.pharmthera.2020.107791
He, C., Shu, B., Zhou, Y., Zhang, R., and Yang, X. (2021). The miR-139-5p/peripheral myelin protein 22 axis modulates TGF-beta-induced hepatic stellate cell activation and CCl4-induced hepatic fibrosis in mice. Life Sci. 276:119294. doi: 10.1016/j.lfs.2021.119294
Hyun, J., Wang, S., Kim, J., Rao, K. M., Park, S. Y., Chung, I., et al. (2016). MicroRNA-378 limits activation of hepatic stellate cells and liver fibrosis by suppressing Gli3 expression. Nat. Commun. 7:10993. doi: 10.1038/ncomms10993
Ikenaga, N., Peng, Z. W., Vaid, K. A., Liu, S. B., Yoshida, S., Sverdlov, D. Y., et al. (2017). Selective targeting of lysyl oxidase-like 2 (LOXL2) suppresses hepatic fibrosis progression and accelerates its reversal. Gut 66, 1697–1708. doi: 10.1136/gutjnl-2016-312473
Ishak, K., Baptista, A., Bianchi, L., Callea, F., De Groote, J., Gudat, F., et al. (1995). Histological grading and staging of chronic hepatitis. J. Hepatol. 22, 696–699.
Ju, B., Nie, Y., Yang, X., Wang, X., Li, F., Wang, M., et al. (2019). miR-193a/b-3p relieves hepatic fibrosis and restrains proliferation and activation of hepatic stellate cells. J. Cell. Mol. Med. 23, 3824–3832. doi: 10.1111/jcmm.14210
Karsdal, M. A., Manon-Jensen, T., Genovese, F., Kristensen, J. H., Nielsen, M. J., Sand, J. M., et al. (2015). Novel insights into the function and dynamics of extracellular matrix in liver fibrosis. Am. J. Physiol. Gastrointest. Liver Physiol. 308, G807–G830. doi: 10.1152/ajpgi.00447.2014
Kim, D., Paggi, J. M., Park, C., Bennett, C., and Salzberg, S. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915. doi: 10.1038/s41587-019-0201-4
Kisseleva, T., and Brenner, D. (2021). Molecular and cellular mechanisms of liver fibrosis and its regression. Nat. Rev. Gastroenterol. Hepatol. 18, 151–166. doi: 10.1038/s41575-020-00372-7
Kweon, Y. O., Goodman, Z. D., Dienstag, J. L., Schiff, E. R., Brown, N. A., Burchardt, E., et al. (2001). Decreasing fibrogenesis: an immunohistochemical study of paired liver biopsies following lamivudine therapy for chronic hepatitis B. J. Hepatol. 35, 749–755. doi: 10.1016/s0168-8278(01)00218-5
Langmead, B., Trapnell, C., Pop, M., and Salzberg, S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10:R25. doi: 10.1186/gb-2009-10-3-r25
Lewis, B. P., Burge, C. B., and Bartel, D. P. (2005). Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 120, 15–20. doi: 10.1016/j.cell.2004.12.035
Ma, L., Yang, X., Wei, R., Ye, T., Zhou, J. K., Wen, M., et al. (2018). MicroRNA-214 promotes hepatic stellate cell activation and liver fibrosis by suppressing Sufu expression. Cell Death Dis. 9:718. doi: 10.1038/s41419-018-0752-1
Machado, M. V., and Diehl, A. M. (2018). Hedgehog signalling in liver pathophysiology. J. Hepatol. 68, 550–562. doi: 10.1016/j.jhep.2017.10.017
Marra, F., and vegliati-Baroni, G. S. (2018). Lipotoxicity and the gut-liver axis in NASH pathogenesis. J. Hepatol. 68, 280–295. doi: 10.1016/j.jhep.2017.11.014
Ratziu, V., Harrison, S. A., Francque, S., Bedossa, P., Lehert, P., Serfaty, L., et al. (2016). Elafibranor, an agonist of the peroxisome proliferator-activated receptor-alpha and -delta, induces resolution of nonalcoholic steatohepatitis without fibrosis worsening. Gastroenterology 150, 1147–1159.e5. doi: 10.1053/j.gastro.2016.01.038
Satishchandran, A., Ambade, A., Rao, S., Hsueh, Y. C., Iracheta-Vellve, A., Tornai, D., et al. (2018). MicroRNA 122, regulated by GRLH2, protects livers of mice and patients from ethanol-induced liver disease. Gastroenterology 154, 238–252.e7. doi: 10.1053/j.gastro.2017.09.022
Stark, R., Grzelak, M., and Hadfield, J. (2019). RNA sequencing: the teenage years. Nat. Rev. Genet. 20, 631–656. doi: 10.1038/s41576-019-0150-2
Sticht, C., De La Torre, C., Parveen, A., and Gretz, N. (2018). miRWalk: an online resource for prediction of microRNA binding sites. PLoS One 13:e0206239. doi: 10.1371/journal.pone.0206239
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47, D607–D613. doi: 10.1093/nar/gky1131
Thomson, D. W., and Dinger, M. E. (2016). Endogenous microRNA sponges: evidence and controversy. Nat. Rev. Genet. 17, 272–283. doi: 10.1038/nrg.2016.20
Trebicka, J., Anadol, E., Elfimova, N., Strack, I., Roggendorf, M., Viazov, S., et al. (2013). Hepatic and serum levels of miR-122 after chronic HCV-induced fibrosis. J. Hepatol. 58, 234–239. doi: 10.1016/j.jhep.2012.10.015
Tsochatzis, E. A., Bosch, J., and Burroughs, A. K. (2014). Liver cirrhosis. Lancet 383, 1749–1761. doi: 10.1016/S0140-6736(14)60121-5
Tsuchida, T., Lee, Y. A., Fujiwara, N., Ybanez, M., Allen, B., Martins, S., et al. (2018). A simple diet- and chemical-induced murine NASH model with rapid progression of steatohepatitis, fibrosis and liver cancer. J. Hepatol. 69, 385–395. doi: 10.1016/j.jhep.2018.03.011
Van Keuren-Jensen, K. R., Malenica, I., Courtright, A. L., Ghaffari, L. T., Starr, A. P., Metpally, R. P., et al. (2016). microRNA changes in liver tissue associated with fibrosis progression in patients with hepatitis C. Liver Int. 36, 334–343. doi: 10.1111/liv.12919
Wang, L., Feng, Y., Xie, X., Wu, H., Su, X. N., Qi, J., et al. (2019). Neuropilin-1 aggravates liver cirrhosis by promoting angiogenesis via VEGFR2-dependent PI3K/Akt pathway in hepatic sinusoidal endothelial cells. EBioMedicine 43, 525–536. doi: 10.1016/j.ebiom.2019.04.050
Wang, X., He, Y., Mackowiak, B., and Gao, B. (2021). MicroRNAs as regulators, biomarkers and therapeutic targets in liver diseases. Gut 70, 784–795. doi: 10.1136/gutjnl-2020-322526
Wei, S., Wang, Q., Zhou, H., Qiu, J., Li, C., Shi, C., et al. (2019). miR-455-3p alleviates hepatic stellate cell activation and liver fibrosis by suppressing HSF1 expression. Mol. Ther. Nucleic Acids 16, 758–769. doi: 10.1016/j.omtn.2019.05.001
Yang, X., Ma, L., Wei, R., Ye, T., Zhou, J., Wen, M., et al. (2020). Twist1-induced miR-199a-3p promotes liver fibrosis by suppressing caveolin-2 and activating TGF-beta pathway. Signal Transduct. Target Ther. 5:75. doi: 10.1038/s41392-020-0169-z
Yu, D. K., Zhang, C. X., Zhao, S. S., Zhang, S. H., Zhang, H., Cai, S. Y., et al. (2015). The anti-fibrotic effects of epigallocatechin-3-gallate in bile duct-ligated cholestatic rats and human hepatic stellate LX-2 cells are mediated by the PI3K/Akt/Smad pathway. Acta Pharmacol. Sin. 36, 473–482. doi: 10.1038/aps.2014.155
Zhang, T., Hu, J., Wang, X., Zhao, X., Li, Z., Niu, J., et al. (2019). MicroRNA-378 promotes hepatic inflammation and fibrosis via modulation of the NF-kappaB-TNFalpha pathway. J. Hepatol. 70, 87–96. doi: 10.1016/j.jhep.2018.08.026
Zhao, W., Yang, A., Chen, W., Wang, P., Liu, T., Cong, M., et al. (2018). Inhibition of lysyl oxidase-like 1 (LOXL1) expression arrests liver fibrosis progression in cirrhosis by reducing elastin crosslinking. Biochim. Biophys. Acta Mol. Basis Dis. 1864(Pt A), 1129–1137. doi: 10.1016/j.bbadis.2018.01.019
Zhou, L., Chen, J., Li, Z., Li, X., Hu, X., Huang, Y., et al. (2010). Integrated profiling of microRNAs and mRNAs: microRNAs located on Xq27.3 associate with clear cell renal cell carcinoma. PLoS One 5:e15224. doi: 10.1371/journal.pone.0015224
Keywords: miRNA, mRNA, liver fibrosis, spontaneous reversal, extracellular matrix, lysyl oxidase
Citation: Tai Y, Zhao C, Lan T, Zhang L, Xiao Y, Tong H, Liu R, Tang C and Gao J (2021) Integrated Analysis of Hepatic miRNA and mRNA Expression Profiles in the Spontaneous Reversal Process of Liver Fibrosis. Front. Genet. 12:706341. doi: 10.3389/fgene.2021.706341
Received: 07 May 2021; Accepted: 23 June 2021;
Published: 22 July 2021.
Edited by:
Dominique Belin, Université de Genève, SwitzerlandReviewed by:
Marie-Luce Bochaton-Piallat, Université de Genève, SwitzerlandMalgorzata Kloc, Houston Methodist Research Institute, United States
Copyright © 2021 Tai, Zhao, Lan, Zhang, Xiao, Tong, Liu, Tang and Gao. 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: Chengwei Tang, shcqcdmed@163.com; Jinhang Gao, gjh731@foxmail.com