1. Introduction
The liver is an important organ of higher animals and contains hepatocytes, bile duct epithelial cells, oval cells, astrocytes, sinus endothelial cells, Kupffer cells, lacuna cells, dendritic cells, and so forth. Of them, the hepatocytes account for about 80% of the total number of liver cells and 75% of the dry weight of the liver [Reference Foster1]. The liver carries out various functions such as metabolism, detoxification, defense, and so forth [Reference Foster1, Reference Forbes and Newsome2]. Additionally, the liver has a strong regeneration ability [Reference Michalopoulos and DeFrances3]. When the liver is injured by surgery, trauma, infection, necrosis, and so forth, the number of hepatocytes decreases sharply, and various feedback signals stimulate the hepatocytes to rapidly change from a static state to a proliferation state, so that the lost or damaged liver tissue is restored, including its structure and function, which is called as liver regeneration [Reference Fausto, Campbell and Riehle4].
The noncoding RNA (ncRNA) is a class of RNA found in etonymeal organisms with important physiological functions [Reference Hombach and Kretz5]. Of these, microRNA (miRNA) is a group of single-stranded RNA with a length of about 22 nucleotides, which can bind to mRNA and inhibit its role [Reference Chipman and Pasquinelli6–Reference Chen, Song and Chen8]. For instance, Chen et al. discovered that miR-21 bound to phosphatidylinositol 3,4,5-trisphosphate 3-phosphatase and dual-specificity protein phosphatase (PTEN) mRNA, inhibited PTEN formation, and promoted hepatocyte proliferation and LR [Reference Chen, Song and Chen8]. Tao et al. reported that miR-612 inhibited the proliferation and migration of hepatocellular carcinoma by acting on AKT2 mRNA [Reference Tao, Wan and Zeng9]. Chen et al. found that miR-1 promoted cell proliferation by targeting histone deacetylase 4, but miR-133 enhanced cell proliferation by inhibiting the serum reactive factor [Reference Chen, Mandel and Thomson10].
Meanwhile, another class of ncRNA was named as circular RNA (circRNA). They are formed by the reverse splicing of the 5′ and 3′ ends of linear RNA and contain multiple miRNA binding sites, which can serve as miRNA sponges [Reference Panda11]. Furthermore, circRNA exhibits important physiological functions [Reference Hansen, Jensen and Clausen12]. For instance, Guo et al. found that circ_03848, circ_08236, circ_13398, and circ_15013 regulated cell proliferation through binding with several miRNAs [Reference Guo, Jin and Chang13]. Li et al. reported that circ137 and circ2270 regulated hepatocyte proliferation by interacting with miR-127 [Reference Li, Guo, Chen, Chang and Xu14]. Shang et al. found that hsa_circ_0005075 and hsa-miR-23b-5p regulated proliferation, invasion, and metastasis of the liver cancer cells through the circRNA-miRNA-mRNA axis [Reference Shang, Li and Liu15].
Neurogenic loci notch homologous protein 3 (NOTCH3), a transcript factor with an N-terminal extracellular domain, an intermediate transmembrane domain, and a C-terminal intracellular domain [Reference Fleming16–Reference Palermo, Checquolo, Bellavia, Talora and Screpanti18], is one of the NOTCH family. One of its functions is to regulate cell proliferation [Reference Tang, Cao and Peng19–Reference Serafin, Persano and Moserle24]. For example, Tang et al. found that the inhibition of NOTCH3 reduced the rate of cell proliferation in osteosarcoma cells [Reference Tang, Cao and Peng19]. Jing et al. found that increased NOTCH3 expression led to goblet cell proliferation, while decreased NOTCH3 expression decreased the proliferation of goblet cells [Reference Jing, Gimenes and Mishra20]. Su et al. found that NOTCH3 methylation reduced cell viability and tumor cell proliferation [Reference Su, Yang and Deng21]. Hassan et al. found that nonsmall-cell lung carcinoma cells displayed high NOTCH3 levels and cell cycle arrest when NOTCH3 expression was reduced [Reference Hassan, Yoshida, Kudoh, Motooka and Ito22, Reference Jung, Stoeck and Guan23]. Serafin et al. found that NOTCH3 overexpression promoted the proliferation of colorectal cancer cells, but the cell proliferation was blocked when NOTCH3 was inhibited [Reference Serafin, Persano and Moserle24].
This study detected the expression changes of competing endogenous RNAs (ceRNAs) by high-throughput biotechnology, analyzed their expression correlation by bioinformatics and systemic biological methods, constructed their interaction networks using Cytoscape 3.2 software, and revealed their role following the aforementioned detection and analysis to investigate their regulatory role to LR. It was found that rno-Rmdn2_0006 and miR-369-3 influenced the expression of the G0 phase- and G1 phase-related genes, which were regulated by NOTCH3 and the status of hepatocytes in the G0 and G1 phases.
2. Materials and Methods
2.1. Preparation of the Rat Liver Regeneration Model Induced by Two-Third Hepatectomy
The two-third hepatectomy (partial hepatectomy, PH) was performed as described by Higgins [Reference Higgins and Anderson25]. The male Sprague–Dawley rats (aged 10 weeks and weighing 250 ± 10 g) were used in the experiments. The liver right lobes of six rats were taken 0, 6, and 24 h after PH and mixed at the corresponding time points. At the same time, the sham operation (SO) was set up as the control. All experimental procedures in this study were carried out according to “The guidelines for the protection and use of experimental animals” published by the Ministry of Science and Technology of People’s Republic of China.
2.2. Isolation and Identification of Rat Hepatocytes
According to Xu and Zhang [Reference Xu and Zhang26], the rats after PH and SO were perfused under aseptic conditions at 9:00–11:00 in the morning. The hepatocytes were digested with collagenase, collected, placed in 60% saline, and centrifuged using Percoll (200 g, 5 min). The cell viability was evaluated using trypan blue staining. The cell purity was evaluated by the fluorescence immunochemistry of Cy3-labeled albumin (ALB) and G6P. The G0-phase hepatocytes were identified by the fluorescence immunochemistry of fluorescein (FITC)-labeled proliferating cell nuclear antigen (PCNA) and Cy3-labeled G6P. The G1-phase hepatocytes were identified by the fluorescence immunochemistry of FITC-labeled cyclin D (CCND1) and Cy3-labeled G6P. S-phase hepatocytes were identified by fluorescent immunochemistry with CCNA2 and Cy3-labeled G6P.
2.3. Large-Scale Quantitative Detection of mRNA
Total RNA was extracted, the ribosomal RNA was removed, the integrity of RNA was evaluated, and the mRNAs were detected. The rat genome 2302.0 chip was used for detection, and the detection was repeated three times. The ratio of gene expression at 0 and 6 h was calculated with the mRNA signal value of regenerating liver cells 24 h after PH recovery as the control [Reference Cimica, Batusic and Haralanova-Ilieva27]. The difference in mRNA expression between the PH and SO groups was calculated using the F test (P value). The genes with P value ≤0.05 were regarded as liver regeneration-related genes. The t test was used to calculate the difference in mRNA expression 0 and 6 h after PH (P value) [Reference Nishimoto28]. When the P value was ≤0.05, the difference was significant; when the P value was ≤0.01, the difference was extremely significant.
2.4. Large-Scale Quantitative Detection of miRNA
The total RNA was extracted, the ribosomal RNA was removed, the integrity of RNA was evaluated, the agarose electrophoresis of samples was performed, the 25-bp RNA was recovered, the library was built, and the single-end sequencing was performed according to TruSeq Stranded Total RNA with Ribo-Zero Gold (Illumina, CA, USA). Then, the quality control of Q20, removal of splice sequences and sequences less than 15 bp and more than 41 bp in length, filtering out of reads containing N base, database analysis, sequence annotation, and quantitative detection of miRNA were conducted. The ratio of miRNA, correlation of liver regeneration (P value), and expression difference (P value) was calculated according to “materials and methods 3.”
2.5. Large-Scale Quantitative Detection of circRNA
Total RNA was extracted; the ribosomal and linear RNAs were removed; the integrity of RNA was evaluated; the library was built; the circRNA was sequenced; the 150/125 bp terminal pairs (reads) were obtained; the sequence was read, matched, quantified, and annotated; the sequence reliability was verified; the chromosome was located; and the circRNA was quantified. Then, the ratio value, correlation of liver regeneration (P value), and expression difference (P value) were determined according to “materials and methods 3.”
2.6. Prediction of G0 Phase- and G1 Phase-Related Genes of Liver Regeneration
The National Center for Biotechnology Information (NCBI) website (https://www.ncbi.nlm.nih.gov/) and Ingenuity Pathway Analysis (IPA) software were used to predict G0 phase- and G1 phase-related genes. Then, the two sets of predicted genes were integrated to get the type and quantity list of G0- and G1-related genes. Subsequently, they were compared with the meaningful expression, differential expression, and liver regeneration-related genes (tables) of “materials and methods 3,” and the genes in both tables were obtained. They were regarded as the G0/G1 phase-related genes of liver regeneration.
2.7. Screening of the Transcription Factors of Liver Regeneration, Which Regulated the Cellular G0 and G1 Phases
The list of transcription factors was obtained by screening the NCBI website. Then, their mRNAs were compared with the mRNA list of “materials and methods 6” of this manuscript to obtain their significant expression; a significant difference was observed in 0 and 6 h during LRI. This study selectively analyzed the NOTCH3 mRNA level and its regulation by miRNA and circRNA.
2.8. Prediction of G0 Phase- and G1 Phase-Related Genes Regulated by NOTCH3
The Cistrome Data Browser website and IPA software were used to predict the G0 phase- and G1 phase-related genes of NOTCH3 regulation. Then, the genes obtained from the aforementioned two tools were integrated to obtain the list of types and numbers of genes related to the G0 and G1 phases regulated by NOTCH3. Subsequently, they were compared with the meaningful expression, differential expression, and liver regeneration-related genes (tables) of “materials and methods 3” in this manuscript, and the genes in both tables were obtained. These genes were regarded as the G0/G1 phase-related genes regulated by NOTCH3 and related to liver regeneration.
2.9. Prediction of miRNA Binding with NOTCH3 mRNA
The IPA software and the miRwalk website (https://mirwalk.umm.uni-heidelberg.de/) were used to predict miRNA binding with NOTCH3 mRNA. The total number miRNAs, named as “theoretical miRNAs,” was listed by integrating the aforementioned two prediction lists. Finally, the “theoretical miRNA” was compared with the detection results of “materials and methods 4” in this manuscript. The miRNA with significant or extremely significant expression difference in 0 and/or 6 h was found, which was called “detection miRNA.” The distribution of the “detection miRNA” in the NOTCH family was obtained by comparing it with miRNA binding to the mRNA of each member of the NOTCH family. The miRNA that only bound to NOTCH3 mRNA was selected as “target miRNA” for analysis.
2.10. Interaction Network Construction of ceRNAs
In this study, miRNA-bound circRNA was predicted using miRanda software. The destination miRNA sequence was input into the miRNA (s) box of miRanda software, followed by clicking on Species, Rat, and Go to get the matching degree and binding energy between the destination miRNA and circRNA. The miRNA/circRNA pairs with matching degree (Max Score) ≥150 and binding energy (Max Energy) ≤–30 were considered as miRNA/circrna pairs. Then, compared with the results of “materials and methods 5,” circRNA with significant or very significant expression difference after 0 and/or 6 h was found, which was called “detection circRNA.” Comparing the distribution of “detection circRNA” in “destination miRNA,” we selected the circRNA that only combined with the latter, called “destination circRNA,” for follow-up analysis.
2.11. Interaction Network Construction of ceRNAs
The ceRNA interaction networks were constructed with Cytoscape 3.2 software. To summarize, the matches of NOTCH3 mRNA and miRNAs and their sponge molecular circRNAs bound with it were listed, and the documents were saved. Then, we clicked File-Import-Network-File and loaded the aforementioned documents into the analysis bar. At the same time, column1 was set as the analysis condition of Source Node and column2 as Target Node, following which we clicked Apply box to obtain their interaction network diagram. Then, we clicked Style-Shape to adjust the shape of the network graph, clicked Fill Color to adjust the color of the network graph, and clicked File and Import to get the network graph of interaction.
2.12. Statistical Analysis
In this study, the ratio values of ceRNAs were calculated in which the gene expression signal values of the control were divided by that of the experimental group. The relative ceRNAs of liver regeneration were determined using the ratio values of ceRNAs and the F test of SPSS 17.0 software. The expression differences of ceRNAs 0 h and 6 h after PH were analyzed using the ratios of ceRNAs and t tests in SPSS 17.0 software [Reference Nishimoto28].
3. Results
3.1. Rat Liver Regeneration and Hepatocytes
In this study, a model of rat two-third hepatectomy (PH) was prepared by the method proposed by Higgins et al. [Reference Higgins and Anderson25]. The liver right lobe of the rat was taken 0, 6, and 24 h after PH (Figure 1(a)). The hepatocytes were separated, and their mRNAs, miRNAs, and circRNAs were detected quantitatively by high-throughput biotechnology (Figure 1(b)). The results showed that the growth index of the regenerated liver was consistent with the reported findings. Also, the activity and purity of the isolated liver cells were ≥95%, and the cytochemical characteristics of the G0-, G1-, and S-phase liver cells were consistent with the report (Figure 1).
3.2. Interactions of NOTCH3 mRNA with miRNAs, circRNAs, and Other mRNAs of Hepatocytes
The combination of ceRNAs was analyzed using miRwalk website and miRanda software. The results showed that 131 miRNAs, which were bound with NOTCH3 mRNA, were significantly expressed 0 and at 6 h after PH. Of these, 55 miRNAs inhibited each other after 0 h and 59 miRNAs after 6 h; 33 miRNAs promoted each other after 0 h and 57 miRNAs after 6 h. On the contrary, 527 circRNAs, which were bound with miRNA, were significantly expressed 0 and 6 h after PH. Of these, 189 circRNAs inhibited each other in 0 h and 195 circRNAs after 6 h; 200 circRNAs promoted each other at 0 h and 288 circRNAs at 6 h. However, 21 genes were related to the G0 or G1 phase and were regulated by NOTCH3. Though all these genes are included in the results earlier, not all of them are discussed in detail in this article. Of these, 8 genes were inhibited by NOTCH3 at 0 h and 10 genes at 6 h, but 4 genes were promoted by NOTCH3 at 0 h and 7 genes at 6 h (Table 1).
* refers to the mRNAs of genes regulated by NOTCH3.
3.3. Interaction of NOTCH3 mRNA with miRNAs, circRNAs, and Other mRNAs of Rat Hepatocytes
The interaction of the aforementioned ceRNAs was analyzed using Cytoscape 3.2 software. The results showed that 59 miRNAs interacted with NOTCH3 mRNA at 0 h (Figure 2(a)) and 38 at 6 h (Figure 2(d)). Furthermore, 59 miRNAs and 144 circRNAs interacted to form 819 interaction pairs at 0 h (Figure 2(b)), and 38 miRNAs and 99 circRNAs interacted to form 303 interaction pairs at 6 h (Figure 2(e)). Continually, NOTCH3 mRNA, 59 miRNAs, and 144 circRNAs interacted to form 819 interaction pairs at 0 h (Figure 2(c)); however, NOTCH3 mRNA, 38 miRNAs, and 99 circRNAs interacted to form 303 interaction pairs at 6 h (Figure 2(f)).
3.4. Expression Correlation of NOTCH3 mRNAs, miRNAs, and circRNAs of Hepatocytes in G0 and G1 Phases
The expression correlation of miRNAs, circRNAs, and NOTCH3 mRNA of hepatocytes in the G0 and G1 phases was analyzed by system biology methods. NOTCH3 mRNA (miR-369-3p, which combined with NOTCH3 mRNA) and rno-Rmdn2_0006 (which combined with miR-369-3p) did not show significant expression change at 0 h, but rno-Rmdn2_0006 was upregulated, and miR-369-3p was downregulated at 6 h (Table 2).
** Extremely significant.
The base sequence of miR-369-3p was analyzed using the miRbase software; its sequence was aauaauacaugguugaucuuu (Table 3. Meanwhile, the mother source gene of rno-Rmdn2_0006 was analyzed by its site on the chromosome and found to be RMDN2.
* Significant; ** extremely significant.
3.5. Expression Correlation of NOTCH3 and NOTCH3-Regulated G0 Phase- and G1 Phase-Related Genes of Hepatocytes
The expression correlation of NOTCH3 and the NOTCH3-regulated G0 phase- and G1 phase-related genes of hepatocytes was analyzed by the system biology methods. NOTCH3 of hepatocytes did not show significant expression change; the expression of G0 phase-inhibited gene CDKN1c, which was regulated by NOTCH3, was upregulated 0 h after PH. However, the expression of NOTCH3 of hepatocytes was upregulated; also, the expression of G1 phase-promoted genes CHUK, DDX24, HES1, NET1, and STAT3, which were regulated by NOTCH3, was upregulated 6 h after PH (Table 2).
3.6. Correlation with the Role of ceRNA and NOTCH3-Regulated G0 Phase- and G1 Phase-Related Genes of Hepatocytes
This study showed that the expression of NOTCH3 mRNA, miR-369-3p, and rno-Rmdn2_0006 did not change 0 h after PH, but the expression of G0 phase-related gene CDKN1a was downregulated, which was inhibited by NOTCH3. Also, the expression of G1 phase-related genes CHUK, DDX24, HES1, NET1, and STAT3 promoted by NOTCH3 was upregulated, suggesting that all of them matched with the physiological status of hepatocytes in the G0 phase (in 0 h, Figure 3). On the contrary, the expression of miR-369-3p, which inhibited NOTCH3 mRNA, was downregulated 6 h after PH. However, the expression of NOTCH3 mRNA and rno-Rmdn2_0006 was upregulated. Also, the expression of the G0 phase-related gene CDKN1c was upregulated, which was promoted by NOTCH3, while the expression of G1 phase-related gene PSEN2 inhibited by NOTCH3 was downregulated, suggesting that all of them matched with the physiological status of hepatocytes in the G1 phase (in 6 h, Figure 3).
4. Discussion
A large number of genes, RNAs, and proteins exist in cells. Also, most physiological activities, including liver regeneration, are controlled by the content changes of the aforementioned elements and their complex interaction and regulation. Although the data obtained using biological high-throughput technology help understand the mechanism of a biological process, the related technology and methods still have limitations. In this study, the aforementioned data were used to reveal the mechanism regulating hepatocytes in the G0 or G1 phase during the rat LRI, providing inspiring results.
The hepatocytes are the main cells of liver structure and function, accounting for about 80% of the total number of hepatic cells and 75% of the total weight of the liver; they can reflect and represent most functions of the liver tissue [Reference Foster1]. In general, most hepatocytes of adult rat liver are in rest, named as the G0 phase. After PH was performed in rats, the hepatocytes of the remnant liver were rapidly activated and switched to the G1 phase synchronously 6 h after PH; 24 h after PH, DNA synthesis in hepatocytes continued synchronously, which was named as the S phase [Reference Fausto, Campbell and Riehle29]. In this study, the two-third hepatectomy (PH), hepatocyte isolation, and biological high-throughput detection of ceRNAs were performed in a sterile room to ensure experimental reliability and avoid contamination interference. The blood in the liver tissue was removed by perfusing with PBS buffer to eliminate the blood interference in the experimental results. The ceRNA expression in hepatocytes in the S phase (24 h after PH) was used as the control of the G0 phase (in 0 h) and G1 phase (in 6 h) to analyze whether the ceRNA expression changes were meaningful in rat liver regeneration. The SO was used as the control to PH to exclude the operation influence on liver regeneration [Reference Xu and Zhang26, Reference Taub30].
Alquda et al. found that highly expressed NOTCH3 promoted cell proliferation through cell cycle proteins [Reference Alqudah, Agarwal, Al-Keilani, Sibenaller, Ryken and Assem31]. Kurakazu et al. found that the downregulation of CDKN1a led to cell cycle arrest [Reference Kurakazu, Akasaki and Hayashida32, Reference Zhang, Qin and Pan33]. Mademtzoglou et al. found that CDKN1c could be used as a cell cycle repressor. The cell proliferation stopped when CDKN1c was highly expressed, but it was enhanced when CDKN1c expression was low [Reference Mademtzoglou, Asakura and Borok34, Reference Chen, Subedi and Chakraborty35]. Liu et al. found that CHUK could promote the G1 phase of the cells [Reference Liu, Wang and Cheng36]. Shi et al. found that DDX24 was overexpressed in proliferating cells and underexpressed in arrested cells [Reference Shi, Dai, Qin and Gu37]. Rani et al. found that HES1 induced cell proliferation. As an important member of the NOTCH signaling pathway, its reduced expression can lead to a decrease in liver regeneration levels [Reference Rani, Greenlaw, Smith and Galustian38–Reference Li, Xu and Wang40]. Ahmad et al. found that the overexpression of NET led to cell proliferation, but underexpression inhibited cell proliferation [Reference Ahmad, Muiwo and Ramachandran41, Reference Zong, Feng and Jiang42]. Bai et al. found that the overexpression of STAT3 promoted cell proliferation, but the proliferation was arrested when STAT3 was inhibited [Reference Bai, Zhou and Xu43, Reference Lu, Zhu and Zhang44]. Janicki et al. found that the overexpression of PSEN2 led to cell arrest, but underexpression promoted cell proliferation [Reference Janicki, Stabler and Monteiro45–Reference Janicki and Monteiro47]. The aforementioned results showed that the NOTCH3 expression and NOTCH3 content influenced the cell phase.
Li et al. found that the low expression of miR-369-3p reduced the proliferation of chondrocytes, but the high expression of miR-369-3p promoted the proliferation of chondrocytes [Reference Li, Dong and Wang48]. Bing et al. found that miR-369-3p acted on CXCR4, and its downregulated expression increased the number of migratory bone marrow mesenchymal stem cells [Reference Bing, Pang and Qu49]. Meanwhile, the prediction using IPA software, miRanda software, and miRwalk website showed that miR-369-3p bound with rno-Rmdn2_0006. However, the role of rno-Rmdn2_0006 was not reported. We noted that RMDN2, as the parent gene of rno-Rmdn2_0006, could regulate microtubule dynamics proteins, and research has shown that RMDN2 is associated with the pathogenesis of lung cancer [Reference Jiang, Chen and Stinnett50]. Based on this, it can be inferred that rno-Rmdn2_0006 may be closely related to cell proliferation.
The aforementioned results suggested that NOTCH3 mRNA was released 6 h after PH, which was inhibited by miR-369-3p. Furthermore, miR-369-3p could be used as an mRNA sponge to inhibit NOTCH3 mRNA translation. Hence, decreased expression of miR-369-3p resulted in more NOTCH3 protein. The latter promoted the expression of G1 phase-related genes and inhibited the expression of G0 phase-related genes; hence, the hepatocytes remained in the G1 phase. On the contrary, the NOTCH3 mRNA combined with miR-369-3p 0 h after PH. Consequently, NOTCH3 was not formed, and the promotion of the expression of G1 phase-related genes and the inhibition of the expression of G0 phase-related genes by NOTCH3 became difficult. As a result, the hepatocytes remained in the G0 phase. The results might help understand the mechanism of circRNA, miRNA, and mRNA regulating the G0 or G1 phase of liver cells.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Disclosure
The authors declare that this manuscript has been reposited in Research Square [Reference Xiayan, Zihui, Yafei and Han51].
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
ZXY and XCS conceived the study. ZXY and WZH designed the study. LYF, GH, GJL, JW, CCF, LJT, ZKC, and XCS performed the study. ZXY and WZH analyzed the data. ZXY, WZH, and LYF wrote the manuscript.
Acknowledgments
This study was supported by the National Natural Science Foundation of China (31572270 and 81771226) and the National Youth Natural Science Foundation (U1804186, U1404312, 31601038, 31201093, and 81200317).