BioMed Research International
Volume 2019 (2019), Article ID 8597953, 11 pages
https://doi.org/10.1155/2019/8597953

Research Article

lncRNA Expression Reveals the Potential Regulatory Roles in Hepatocyte Proliferation during Rat Liver Regeneration

Haijing Bai,1,2 Wei Jin,1,2 Jianlin Guo,1,2 Yi Ding,1,2 Cuifang Chang,1,2 Xueqiang Guo,1,2 Yaping Song,1,2 Jingbo Zhang,³ and Cunshuan Xu1,2

¹College of Life Science, Henan Normal University, Xinxiang 453007, Henan Province, China
²State Key Laboratory Cultivation Base for Cell Differentiation Regulation and Henan Engineering Laboratory for Bioengineering and Drug Development, Xinxiang 453007, Henan Province, China
³Institute of Basic Medical Sciences, Chinese Academy of Medical Sciences, Beijing 100005, China

Correspondence should be addressed to Jingbo Zhang; jingbozh@hotmail.com and Cunshuan Xu; cellkeylab@126.com

Received 10 June 2019; Revised 29 August 2019; Accepted 6 September 2019; Published 11 November 2019

Academic Editor: Monica Fedele

Copyright © 2019 Haijing Bai et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Liver regeneration is a tissue growth process after loss or injury of liver tissue, which is a compensatory hyperplasia rather than true regeneration, mainly depending on hepatocyte proliferation. Currently, a large number of studies on hepatocyte proliferation have been conducted. However, studies on the regulation of long noncoding RNA (lncRNA) on hepatocyte proliferation are still limited. To identify specially expressed lncRNA during rat liver regeneration, high-throughput sequencing technology was performed, and a total of 2446 lncRNAs and 4091 mRNAs were identified as significantly differentially expressed. Gene ontology (GO) enrichment analysis was performed to analyze the role of differentially expressed mRNAs, and 695 mRNAs were identified to be related to cell proliferation. Then, an lncRNA-mRNA coexpression network based on the differentially expressed lncRNAs and proliferation-related genes was constructed to analyze the potential function of lncRNAs on hepatocyte proliferation, and ten lncRNAs, NONRATT003557.2, NONRATT005357.2, NONRATT003292.2, NONRATT001466.2, NONRATT003289.2, NONRATT001047.2, NONRATT005180.2, NONRATT004419.2, NONRATT005336.2, and NONRATT005335.2, were selected as key regulatory factors, which may play crucial roles in hepatocyte proliferation during rat liver regeneration. Finally, a protein-protein interaction (PPI) network was established to illuminate the interaction between proliferation-related genes, and ten hub genes (Aurkb, Cdk1, Cdc20, Bub1b, Mad2l1, Kif11, Prc1, Ccna2, Top2a, and Ccnb1) were screened with the MCC method in the PPI network, which may be important biomarkers involved in the hepatocyte proliferation during rat liver regeneration. These results may provide clues for a more comprehensive understanding of the molecular mechanism of hepatocyte proliferation during rat liver regeneration.

1. Introduction

The vast majority of the eukaryotic genomes are transcribed into noncoding RNAs, which can be divided into small noncoding RNAs (<200 bp) and long noncoding RNAs (lncRNAs; ≥200 nt) based on transcript size [1]. lncRNAs can be divided into five categories: sense, antisense, bidirectional, intronic, and intergenic [2]. Initially, lncRNAs were considered to be “dark matter,” as byproduct of transcription of RNA polymerase II, with no biological function. In the past 20 years, genome-wide identification of lncRNAs has become possible with the development of high-throughput technology of RNA-seq, many of which are involved in various biological functions [3]. Increasing lncRNAs have been found to play a critical role in biological processes, like development [4], gene transcriptional regulation [5], chromatin regulation [6], epithelial-to-mesenchymal transition (EMT) [7], and cell proliferation [8].

In rodents and humans, the liver can grow rapidly after partial hepatectomy (PH) or acute chemical injury. This growth process is known as LR, which is a compensatory hyperplasia rather than true regeneration [9]. During LR, quiescent hepatocytes undergo one or two rounds of replication and then return to a nonproliferative state [10]. This process is very complex and regulated by a variety of growth factors, cytokines and noncoding RNAs [11, 12]. Therefore, the study of the molecular mechanism of hepatocyte proliferation is crucially important to understand the process of LR and provide clues for the treatment of liver diseases. Several recent studies have shown that lncRNAs play a critical role in hepatocyte proliferation [12–14]. However, the study of hepatocyte proliferation during LR is still largely unknown.

In the present study, high-throughput sequencing technology was used to identify DE lncRNAs and mRNAs during rat LR. Then, functional enrichment analysis of DE mRNAs was performed to screen proliferation-related genes. Finally, the lncRNA-mRNA coexpression network and PPI network were constructed based on DE lncRNAs and proliferation-related genes to elucidate the molecular mechanism of hepatocyte proliferation during LR. These results lay a foundation for understanding the regulatory function of lncRNAs on hepatocyte proliferation and provide an important clue for the study of the LR process.

2. Materials and Methods

2.1. Preparation of 2/3 Hepatectomy Model

The healthy adult male Sprague Dawley (SD) rats weighing 210∼250 g were provided by the Laboratory Animal Center of Zhengzhou University (Zhengzhou, China). These rats were raised in a controlled temperature room of 19∼23°C with a relative humidity of 50∼70% and an illumination time of 12 h/d (8 : 00 to 20 : 00) and permitted to freely have water and food. A total of 60 rats were taken for the experiment with six rats per group: nine PH groups and one normal group (CG). The rats in PH groups were conducted 2/3 PH according to the method of Xu et al. The rats were anesthetized and condemned to death at 0, 2, 6, 12, 24, 30, 36, 72, 120, and 168 h after operation. The right liver lobes of six rats were mixed at each time point and stored at −80°C. All operations conformed to the Animal Protection Law of China and Animal Ethics.

2.2. RNA Sequencing

RNA sequencing was performed by the Shanghai OE Biotech (Shanghai, China). In brief, the mirVana miRNA Isolation Kit (Ambion) was used to extract the total RNA from liver tissues. The TruSeq Stranded Total RNA with Ribo-Zero Gold (Illumina) was used to construct cDNA libraries. The purified cDNA libraries were sequenced on Illumina HiSeq 2500 following the manufacturer’s instruction. After filtrating the adaptor and low-quality reads, clean reads were obtained for subsequent analysis. The reads were matched to the rat reference genome using the hisat2 (v2.2.1.0) software. The StringTie2 (v1.3.3b) software was used to splice the aligned reads. lncRNA identification included two categories: one is known lncRNA, which completely matches with the known lncRNAs, and the other is candidate lncRNA lacking protein-coding ability, whose length is greater than 200 bp and exon is greater than or equal to 2. The software CPC (v0.9-r2), CNCI (v1.0), Pfam (v30), and PLEK (v1.2) were used to predict the protein-coding ability of transcripts. The expression of transcription was calculated by the fragments per kilobase of exon per million reads mapped (FPKM) method using the bowtie2 (v2.2.9) and eXpress (v1.5.1) software.

2.3. Identification of Differentially Expressed lncRNAs

The counts of lncRNAs in each sample were standardized by the baseMean value using the DESeq (1.18.0) software. Differentially expressed (DE) lncRNAs were identified with fold change ≥2 or ≤0.5 and as the threshold. All DE lncRNAs in nine PH groups underwent hierarchical clustering analysis using the cluster3.0 and treeview software.

2.4. GO Enrichment Analysis

Gene ontology (GO) enrichment analysis of the DE mRNAs was conducted using David Bioinformatics Resources 6.8 (https://david.ncifcrf.gov/). The enrichment analysis consisted of three parts: biological process (BP), molecular function (MF), and cellular component (CC). is considered statistically significant, which was calculated by the EASE score.

2.5. lncRNA-mRNA Coexpression Analysis

To explore the relationship between DE lncRNAs and proliferation-related genes, a coexpression analysis was performed. Pearson’s correlation coefficients (PCCs) were calculated between the DE lncRNAs and the proliferation-related genes, and only lncRNA-mRNA pairs with PCC ≥0.8 and were selected and considered as coexpression. Then, these lncRNA-mRNA pairs were used to construct a coexpression network, which was visualized by the Cytoscape v3.6.1 software. The node degree was determined by the number of directly connected neighbors to the topological property of the network.

2.6. Construction of PPI Network and Screening of the Key Gene

To illustrate interactions between proliferation-related DE mRNAs, the string database (https://string-db.org/) was used to construct a protein-protein interaction (PPI) network. Only the interacting pairs with combined score ≥0.4 were selected and considered to be significant. The PPI network was visualized by using the Cytoscape v3.6.1 software. Then, a Cytoscape plugin cytoHubba was used to identify the key genes adopting the MCC method.

2.7. RT-PCR Validation

RNA-seq results were validated by RT-PCR, and the primers are listed in Table 1. cDNA was synthesized using the cDNA Reverse Transcription Kit (Takara, Tokyo, Japan). The qRT-PCR was performed using Q-SYBR Green Supermix (Bio-Rad). Primers were also designed to amplify β-actin as an endogenous control. The expression of each lncRNA was represented as fold change using 2−ΔΔCt methods.

Table 1: Primer sequence.

Gene name
Sense primer
Antisense primer
Product length

NONRATT003289.2
AATGCCAGGCCATGCTAAGGAC
GCTCTGCCAGGTGACTGCTTC
180

NONRATT001466.2
TCTGCTGTTGACATTGGCGAAGG
CTAGCATGTGAGAGGTGACGTGAC
189

NONRATT004419.2
AGCCTCCTGAGTCCTGGAATTCTG
GTGAGTCGTGAGTGAGCTGAAGTG
133

NONRATT005336.2
AAGCTCAACACTGCCTGAGTCTTC
GCATGAGCCTTGGAGGACATCTG
112

3. Results

3.1. Sequencing and Identification of lncRNAs during Rat LR

To identify the expression of lncRNAs during rat LR, 10 cDNA libraries were constructed from the regeneration rat liver at different time points after surgery. The Illumina HiSeq X Ten platform was used to sequence these cDNA libraries, and a total of 1116M raw reads were produced. After filtering adaptor sequences and low-quality reads, 1084.09M clean reads were obtained. The percentage of clean reads varied from 90.73 to 95.29%, and the percentage of GC content varied from 48.77 to 51.08% in each library. Approximately 95.91–97.83% of clean reads were selected for further research after mapping the clean reads to the rat reference genome (Table 2).

Table 2: Overview of RNA-seq.

Sample
raw_reads
raw_bases
clean_reads
clean_bases
valid_bases (%)
GC (%)
Mapped reads

0 h
119.65M
17.95G
115.99M
17.05G
94.99
50.92
111716504 (96.32%)

2 h
120.14M
18.02G
116.69M
17.17G
95.25
50.57
113657852 (97.40%)

6 h
119.57M
17.94G
116.13M
17.09G
95.29
50.26
113238482 (97.51%)

12 h
119.91M
17.99G
117.68M
16.68G
92.75
49.88
114912884 (97.65%)

24 h
119.27M
17.89G
117.28M
16.23G
90.73
48.77
114629815 (97.74%)

30 h
88.95M
13.34G
85.96M
12.64G
94.74
50.29
83564746 (97.21%)

36 h
120.15M
18.02G
116.44M
17.16G
95.22
50.62
113908105 (97.83%)

72 h
91.48M
13.72G
88.07M
12.9G
93.98
51.08
86032773 (97.68%)

120 h
120.82M
18.12G
116.34M
17.12G
94.47
50.97
112635056 (96.82%)

168 h
96.06M
14.41G
93.51M
13.73G
95.28
50.71
89681397 (95.91%)

3.2. Identification of DE lncRNAs and DE mRNAs

The expression abundance of the lncRNAs was evaluated by FPKM (fragments per kB per million reads) using DESeq. A total of 2446 lncRNAs were determined to be differentially expressed during rat LR, with 1120 upregulated, 731 downregulated, and 595 up/downregulated lncRNAs (Figure 1(a) and Table S1). To explore the similarity of gene expression, hierarchical clustering was adopted to analyze the expression of DE lncRNAs (Figure 2(a)). To further explore the interactions of DE lncRNAs at different stages: initial stage (2–6 h), proliferation stage (12–72 h), and termination stage (120–168 h), a Venn diagram was constructed using these DE lncRNAs (Figure 2(c)). Among them, 272 DE lncRNAs were common to all three stages.

(a)

(b)

Figure 1: Number of DE lncRNAs and DE mRNAs at different time points during rat LR. (a) Number of upregulated and downregulated lncRNAs at nine time points. (b) Number of upregulated and downregulated mRNAs at nine time points.

(a)

(b)

(c)

(d)

Figure 2: Analyses of DE lncRNAs in the RNA-seq libraries. (a) Hierarchical clustering analysis of DE lncRNAs at nine time points of rat LR. (b) Hierarchical clustering analysis of DE mRNAs at nine time points of rat LR. (c) Venn diagram showing the DE lncRNAs at three stages of rat LR. (d) Venn diagram showing the DE mRNAs at three stages of rat LR.

Through high-throughput RNA-seq, the expression profile of 28635 mRNAs was measured. Among them, 4091 mRNAs were found to be differentially expressed, of which 2,256 were upregulated, 1,686 were downregulated, and 149 were up/downregulated (Figure 1(b) and Table S2). Hierarchical clustering was employed to analyze the expression similarity of DE mRNAs (Figure 2(b)). Venn analysis was conducted to explore the differences of DE mRNAs at different stages (Figure 2(d)).

3.3. GO and KEGG Enrichment Analysis

GO enrichment analysis was employed to determine the function of deregulated genes during rat LR, and 689 significant GO terms (451 under BP, 116 under CC, and 122 under MF) were enriched. The top 30 GO terms of the three groups are listed in Figure 3(a) (Table S3). The most enriched BP terms were response to drug, cell division, and chromosome segregation. As for CC, the most enriched terms were cytoplasm, nucleus, and nucleoplasm. The most enriched MF terms were related to binding activity. Of these BP terms, 41 were associated with cell proliferation involving 695 mRNAs, corresponding to 585 genes (Figure 3(b) and Table S3).

(a)

(b)

Figure 3: GO enrichment analysis of the DE mRNAs. (a) The top 30 significant GO terms in GO enrichment analysis at value <0.05. Green represents biological processes, red represents cellular components, and blue represents molecular functions. (b) GO terms associated with cell proliferation at the biological process level.

3.4. Coexpression Network Construction Based on DE lncRNAs and Proliferation-Related DE mRNAs

The coexpression analysis of the screened proliferation-related DE mRNAs and DE lncRNAs found 18343 significant coexpression pairs, of which 17547 pairs were positively correlated (PCC ≥ 0.8) and 796 pairs were negatively correlated (PCC ≤ −0.8) (Table S4). Subsequently, a coexpressed network was constructed using screened lncRNA-mRNA pairs, and it was found that some lncRNAs could interact with multiple mRNAs. According to nodes and connectivity, the top 10 lncRNAs were selected, which may exert important regulatory roles in hepatocyte proliferation, lncRNAs NONRATT003557.2 (degree = 178), NONRATT005357.2 (degree = 178), NONRATT003292.2 (degree = 177), NONRATT001466.2 (degree = 176), NONRATT003289.2 (degree = 176), NONRATT001047.2 (degree = 175), NONRATT005180.2 (degree = 175), NONRATT004419.2 (degree = 173), NONRATT005336.2 (degree = 173), and NONRATT005335.2 (degree = 169) (Figure 4). The expression levels of the top 10 lncRNAs during rat LR are shown in Table 3. Among them, NONRATT005357.2 had the highest difference at 24 h in regenerating the liver, and the fold change was 119.

Figure 4: The top 10 lncRNAs and their interaction mRNAs. V and ellipse represent lncRNAs and mRNAs, respectively. Upregulated genes are labeled in red, downregulated genes are labeled in green, and up/downregulated genes are labeled in yellow.

Table 3: Expression levels of the top 10 lncRNAs during rat LR.

lncRNAs
0 h
2 h
6 h
12 h
24 h
30 h
36 h
72 h
120 h
168 h

NONRATT003557.2
1
1
1
1
14.67
15.41
10.44
1
1
1

NONRATT005357.2
1
1
1
9.79
119.12
59.69
71.42
32.46
1
1

NONRATT003292.2
1
1
1
1
43.68
40.42
28.50
18.78
1
1

NONRATT001466.2
1
1
1
1
20.37
23.28
24.66
13.95
1
1

NONRATT003289.2
1
1
1
1
69.44
53.97
38.69
23.26
6.55
1

NONRATT001047.2
1
1
1
1
1
1
4.86
1
1
1

NONRATT005180.2
1
1
1
1
9.81
9.02
10.33
1
1
1

NONRATT004419.2
1
1
1
1
1
6.87
6.87
7.5
1
1

NONRATT005336.2
1
1
1
1
1
6.83
1
1
1
1

NONRATT005335.2
1
1
1
1
6.89
6.44
5.27
1
1
1

Text in bold denotes the expression level higher than the control.

3.5. Construction of PPI Network Based on Proliferation-Related Genes

To better elucidate the interaction network of proliferation-related genes, a PPI network containing 551 nodes with scores greater than or equal to 0.4 was constructed using the string database (Figure 5(a)). Moreover, a cytoHubba plugin was used to select the top 10 key genes from the PPI network using the MCC method. The 10 key genes were Aurkb, Cdk1, Cdc20, Bub1b, Mad2l1, Kif11, Prc1, Ccna2, Top2a, and Ccnb1 (Figure 5(b)). The expression levels of the 10 key genes during rat LR are shown in Table 4. The median of multiple mRNAs corresponding to one gene is taken as its expression level. Among them, Cdk1 had the highest difference at 24 h in regenerating the liver, and the fold change was 109.

(a)

(b)

Figure 5: PPI network of proliferation-related genes. (a) PPI network of 551 genes. (b) Top 10 hub DE mRNAs of the PPI network. The node color changes gradually from yellow to red in the ascending order according to scores with the MCC method.

Table 4: Expression levels of 10 key genes during rat LR.

Genes
0 h
2 h
6 h
12 h
24 h
30 h
36 h
72 h
120 h
168 h

Aurkb
1
0.02
1
1
41.88
26.57
23.96
14.00
3.72
1

Cdk1
1
1
1
1
109.01
66.43
80.87
42.31
3.63
1

Cdc20
1
1
1
1
36.85
27.36
25.84
13.01
2.83
1

Bub1b
1
1
1
1
38.34
26.78
22.37
12.92
5.08
1

Mad2l1
1
1
1
1
20.26
13.71
13.11
7.02
1
1

Kif11
1
1
1
1
20.48
13.24
11.22
7.31
2.12
1

Prc1
1
1
1
1
46.88
29.75
27.15
16.27
3.59
1

Ccna2
1
0.05
1
1
49.97
31.85
30.16
15.45
3.92
1

Ccnb1
1
1
1
1
60.31
37.31
37.82
18.90
3.45
1

Top2a
1
0.18
1
1
49.62
34.69
23.24
14.75
3.03
1

Text in bold and text in italics denote the expression level higher and lower than the control, respectively.

3.6. Quantitative Real-Time PCR Validation

We validated the high-throughput RNA-seq results by performing qRT-PCR analysis of differentially expressed lncRNAs; the expression patterns revealed similar conclusions (Figure 6). Our RNA-seq results showed lncRNAs NONRATT003289.2, NONRATT001466.2, NONRATT004419.2, and NONRATT005336.2 were upregulated during rat LR.

(a)

(b)

(c)

(d)

Figure 6: Validation of differentially expressed lncRNAs by qRT-PCR.

4. Discussion

In rodents and humans, the liver can grow rapidly after partial hepatectomy or acute chemical injury. This growth process is known as LR, which is a compensatory hyperplasia rather than true regeneration, mainly depending on the proliferation of hepatocytes [9]. To explore the regulatory roles of lncRNAs in hepatocyte proliferation during rat LR, high-throughput RNA-seq was performed to identify lncRNAs and mRNAs. In the present research, 2446 DE lncRNAs and 4091 DE mRNAs were identified. To investigate the function of DE mRNAs during rat LR, GO enrichment analysis was performed. The result indicated that a large number of GO terms were associated with response to stress, cell proliferation, oxidation-reduction, regulation of transcription, metabolism, and apoptosis, which were considered to be important activities in LR [15–17]. Of these GO terms, 41 were associated with cell proliferation, involving 695 mRNAs.

Through coexpression network analysis, a total of 18343 coexpression pairs were obtained based on DE lncRNAs and proliferation-related DE mRNAs. Among these coexpression pairs, 95.7% were positively correlated and 4.3% were negatively correlated. These results indicate that lncRNAs may regulate the expression of genes mainly in a positive way. It was also found that some lncRNAs could be coexpressed with multiple mRNAs, suggesting that lncRNAs could regulate multiple mRNAs. According to nodes and connections, the top 10 lncRNAs (NONRATT003557.2, NONRATT005357.2, NONRATT003292.2, NONRATT001466.2, NONRATT003289.2, NONRATT001047.2, NONRATT005180.2, NONRATT004419.2, NONRATT005336.2, and NONRATT005335.2) were selected. These lncRNAs were connected to 158 common DE mRNAs, respectively, such as Slc25a16, Cdk1, Rfc3, Brca2, Mcm7, Cdca7, Hat1, Chaf1b, Rfc4, and Mcm4. These findings indicated that these lncRNAs and their interaction genes may play an important role in hepatocyte proliferation during rat LR.

cytoHubba, a Cytoscape plugin, could provide 11 topological analysis methods to explore important nodes in biological networks including maximal clique centrality (MCC). MCC, a new method, has a better precision in identifying hub proteins. Through the MCC method, 10 hub mRNAs were selected including Aurkb, Cdk1, Cdc20, Bub1b, Mad2l1, Kif11, Prc1, Ccna2, Top2a, and Ccnb1. Aurkb (aurora kinase B) was involved in regulation of chromatin, glycolysis metabolism, regulation of telomerase activity, regulation of zygote development, and turnover of kinetochore microtubules [18–22]. It was also associated with AKT signaling pathway, MAPK signaling pathway, inflammatory pathway, DNA damage response pathway, Wnt signaling pathway, and PI3K/Akt/NF-κB signaling pathway [21, 23–27]. Cdk1 (cyclin-dependent kinase 1) was involved in cell cycle, mitosis, DNA end resection, homologous recombination, DNA damage checkpoint, anaphase spindle dynamics, cytokinesis, chromosome biorientation, apoptosis, cell proliferation and survival, and DNA synthesis [28–35]. It also played an important role in AMPK signaling pathway, AKT signaling pathway, Golgi checkpoint signaling, PDK1-PI3K/Akt signaling pathway, type I interferon signaling, and p53 signaling pathway [36–40]. Cdc20 (cell division cycle 20), an APC activator protein, could regulate mitosis, cell cycle, and presynaptic differentiation [41–44]. It was also involved in Wnt/β-catenin signaling pathway and p38 MAPK signaling pathway. Bub1b (BUB1 mitotic checkpoint serine/threonine kinase B) and Mad2l1 (mitotic arrest deficient 2 like 1) were involved in the spindle checkpoint during mitosis [45–48]. Kif11 (kinesin family member 11), a microtubule motor, played a vital role in regulating the transport of β-actin mRNA and cell motility through physically interacting with ZBP1, which could govern the direction of migration by responding to directional cues in chemotaxis [49, 50]. Prc1 (protein regulator of cytokinesis 1), a cell cycle protein, played an important role in ensuring proper cell division through directly acting on protein (FIP), and P53 could regulate the transcription of Prc1 [51, 52]. Ccna2 (cyclin A2) was a key factor in cell cycle, and Ccna2 repression regulated by miR-22 could inhibit HCC cell proliferation and tumorigenesis [53, 54]. It was also a prognostic biomarker for several cancers, such as ER+ breast cancer, pancreatic ductal adenocarcinoma, and colorectal cancer [55–57]. Top2a (DNA topoisomerase II alpha), a key enzyme in DNA replication, was recruited to ultrafine anaphase bridges (UFBs) by TopBP1 to ensure faithful separation of sister chromatids, and these proteins played an important role in maintaining genome stability [58]. Top2a was overexpressed in hepatocellular carcinoma and associated with early age onset, shorter patient survival, and chemoresistance [59]. Ccnb1 (cyclin B1) was a key factor in cell cycle and regulated by STAT3 via the E2F modulating G2-M phase checkpoint [60]. High expression of Ccnb1 was associated with poor prognosis in HCC patients, and knockdown of Ccnb1 could significantly inhibit cell proliferation, migration, and invasion in HCC [61]. Therefore, these genes may play a vital role in hepatocyte proliferation during LR.

5. Conclusions

In this study, the comprehensive expression abundance of lncRNAs and mRNAs was identified by RNA-seq analysis during rat LR. The lncRNA-mRNA coexpression network and PPI network based on lncRNAs and proliferation-related genes were constructed, and 10 key lncRNAs and 10 key mRNAs were determined that may play crucial roles in hepatocyte proliferation. Our study provides a new idea to better understand the mechanism of LR.

Data Availability

The data used to support the findings of this study are included within the supplementary information file(s).

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This study was financially supported by the National Natural Science Foundation of China (no. 31572270) and the National Natural Science Foundation of China Youth Science Fund Project (no. 31601038).

Supplementary Materials

Table S1: results of analyzing differentially expressed lncRNAs at 2, 6, 12, 24, 30, 36, 72,120, and 168 h after PH was compared with CG. Table S2: results of analyzing differentially expressed mRNAs at 2, 6, 12, 24, 30, 36, 72,120, and 168 h after PH was compared with CG. Table S3: results of GO enrichment analysis of differentially expressed mRNAs. Table S4: screened coexpression pairs between DE lncRNAs and proliferation-related DE mRNAs based on the Pearson correlation coefficient. (Supplementary Materials)

References

K. Zhang, Z.-M. Shi, Y.-N. Chang, Z.-M. Hu, H.-X. Qi, and W. Hong, “The ways of action of long non-coding RNAs in cytoplasm and nucleus,” Gene, vol. 547, no. 1, pp. 1–9, 2014.

C. P. Ponting, P. L. Oliver, and W. Reik, “Evolution and functions of long noncoding RNAs,” Cell, vol. 136, no. 4, pp. 629–641, 2009.

C. Lee and N. Kikyo, “Strategies to identify long noncoding RNAs involved in gene regulation,” Cell & Bioscience, vol. 2, no. 1, p. 37, 2012.

A. N. Kallen, X.-B. Zhou, J. Xu et al., “The imprinted H19 lncRNA antagonizes let-7 microRNAs,” Molecular Cell, vol. 52, no. 1, pp. 101–112, 2013.

L. Yang, C. Lin, C. Jin et al., “lncRNA-dependent mechanisms of androgen-receptor-regulated gene activation programs,” Nature, vol. 500, no. 7464, pp. 598–602, 2013.

J. M. Engreitz, A. Pandya-Jones, P. McDonel et al., “The Xist lncRNA exploits three-dimensional genome architecture to spread across the X chromosome,” Science, vol. 341, no. 6147, p. 1237973, 2013.

W.-C. Liang, W.-M. Fu, C.-W. Wong et al., “The lncRNA H19 promotes epithelial to mesenchymal transition by functioning as miRNA sponges in colorectal cancer,” Oncotarget, vol. 6, no. 26, pp. 22513–22525, 2015.

Z. H. Bian, L. G. Jin, J. W. Zhang et al., “LncRNA-UCA1 enhances cell proliferation and 5-fluorouracil resistance in colorectal cancer by inhibiting miR-204-5p,” Scientific Reports, vol. 6, no. 1, 2016.

N. Fausto and J. S. Campbell, “The role of hepatocytes and oval cells in liver regeneration and repopulation,” Mechanisms of Development, vol. 120, no. 1, pp. 117–130, 2003.

N. Fausto, A. D. Laird, and E. M. Webber, “Liver regeneration. 2. Role of growth factors and cytokines in hepatic regeneration,” The FASEB Journal, vol. 9, no. 15, pp. 1527–1536, 1995.

C.-G. Huh, V. M. Factor, A. Sanchez, K. Uchida, E. A. Conner, and S. S. Thorgeirsson, “Hepatocyte growth factor/c-met signaling pathway is required for efficient liver regeneration and repair,” Proceedings of the National Academy of Sciences, vol. 101, no. 13, pp. 4477–4482, 2004.

D. Xu, F. Yang, J.-H. Yuan et al., “Long noncoding RNAs associated with liver regeneration 1 accelerates hepatocyte proliferation during liver regeneration by activating Wnt/β-Catenin signaling,” Hepatology, vol. 58, no. 2, pp. 739–751, 2013.

Q. Gao, Y. Y. Gu, Y. N. Jiang et al., “Long non-coding RNA Gm2199 rescues liver injury and promotes hepatocyte proliferation through the upregulation of ERK1/2,” Cell Death & Disease, vol. 9, no. 6, 2018.

L. L. Huang, S. S. Damle, S. Booten et al., “Partial hepatectomy induced long noncoding RNA inhibits hepatocyte proliferation during liver regeneration,” PLoS One, vol. 10, no. 7, 2015.

R. Dayoub, A. Vogel, J. Schuett et al., “Nrf2 activates augmenter of liver regeneration (ALR) via antioxidant response element and links oxidative stress to liver regeneration,” Molecular Medicine, vol. 19, no. 1, pp. 237–244, 2013.

M. A. Fernandez-Rojo, C. Restall, C. Ferguson et al., “Caveolin-1 orchestrates the balance between glucose and lipid-dependent energy metabolism: implications for liver regeneration,” Hepatology, vol. 55, no. 5, pp. 1574–1584, 2012.

Y. Iimuro, T. Nishiura, C. Hellerbrand et al., “NFkappaB prevents apoptosis and liver dysfunction during liver regeneration,” Journal of Clinical Investigation, vol. 101, no. 4, pp. 802–811, 1998.

L. L. Hall, M. Byron, G. Pageau, and J. B. Lawrence, “AURKB-mediated effects on chromatin regulate binding versus release of XIST RNA to the inactive chromosome,” The Journal of Cell Biology, vol. 186, no. 4, pp. 491–507, 2009.

R. Zhou, Y. Zhang, G. Du et al., “Down-regulated let-7b-5p represses glycolysis metabolism by targeting AURKB in asthenozoospermia,” Gene, vol. 663, pp. 83–87, 2018.

J.-P. Mallm and K. Rippe, “Aurora kinase B regulates telomerase activity via a centromeric RNA in stem cells,” Cell Reports, vol. 11, no. 10, pp. 1667–1678, 2015.

L. Xu, T. Liu, F. Han et al., “AURKB and MAPK involvement in the regulation of the early stages of mouse zygote development,” Science China Life Sciences, vol. 55, no. 1, pp. 47–56, 2012.

D. Cimini, X. Wan, C. B. Hirel, and E. D. Salmon, “Aurora kinase promotes turnover of kinetochore microtubules to reduce chromosome segregation errors,” Current Biology, vol. 16, no. 17, pp. 1711–1718, 2006.

J. Yang, T. Ikezoe, C. Nishioka, K. Udaka, and A. Yokoyama, “Bcr-Abl activates AURKA and AURKB in chronic myeloid leukemia cellsviaAKT signaling,” International Journal of Cancer, vol. 134, no. 5, pp. 1183–1194, 2014.

A. Madejón, J. Sheldon, I. Francisco-Recuero et al., “Hepatitis C virus-mediated Aurora B kinase inhibition modulates inflammatory pathway and viral infectivity,” Journal of Hepatology, vol. 63, no. 2, pp. 312–319, 2015.

A. Zekri, S. H. Ghaffari, M. Yaghmaie et al., “Inhibitor of aurora kinase B induces differentially cell death and polyploidy via DNA damage response pathways in neurological malignancy: shedding new light on the challenge of resistance to AZD1152-HQPA,” Molecular Neurobiology, vol. 53, no. 3, pp. 1808–1823, 2016.

K. B. Lee, H. Jin, S. Ye, B. H. Park, and S. M. Kim, “Recombinant human bone morphogenetic protein-2 inhibits gastric cancer cell proliferation by inactivating Wnt signaling pathway via c-Myc with aurora kinases,” Oncotarget, vol. 7, no. 45, pp. 73473–73485, 2016.

L. B. Zhu, J. Jiang, X. P. Zhu et al., “Knockdown of Aurora-B inhibits osteosarcoma cell invasion and migration via modulating PI3K/Akt/NF-kappaB signaling pathway,” International Journal of Clinical and Experimental Pathology, vol. 7, no. 7, pp. 3984–3991, 2014.

D. Santamaría, C. Barrière, A. Cerqueira et al., “Cdk1 is sufficient to drive the mammalian cell cycle,” Nature, vol. 448, no. 7155, pp. 811–815, 2007.

O. Gavet and J. Pines, “Progressive activation of CyclinB1-Cdk1 coordinates entry to mitosis,” Developmental Cell, vol. 18, no. 4, pp. 533–543, 2010.

G. Ira, A. Pellicioli, A. Balijja et al., “DNA end resection, homologous recombination and DNA damage checkpoint activation require CDK1,” Nature, vol. 431, no. 7011, pp. 1011–1017, 2004.

S. P. Wheatley, E. H. Hinchcliffe, M. Glotzer, A. A. Hyman, G. Sluder, and Y.-L. Wang, “CDK1 inactivation regulates anaphase spindle dynamics and cytokinesis in vivo,” The Journal of Cell Biology, vol. 138, no. 2, pp. 385–393, 1997.

T. Tsukahara, Y. Tanno, and Y. Watanabe, “Phosphorylation of the CPC by Cdk1 promotes chromosome bi-orientation,” Nature, vol. 467, no. 7316, pp. 719–723, 2010.

R. M. Golsteyn, “Cdk1 and Cdk2 complexes (cyclin dependent kinases) in apoptosis: a role beyond the cell cycle,” Cancer Letters, vol. 217, no. 2, pp. 129–138, 2005.

P. Liu, T. P. Kao, and H. Huang, “Cdk1 promotes cell proliferation and survival via phosphorylation and inhibition of FOXO1 transcription factor,” Oncogene, vol. 27, no. 34, pp. 4733–4744, 2008.

M. Pagano, “Control of DNA synthesis and mitosis by the Skp2-p27-Cdk1/2 axis,” Molecular Cell, vol. 14, no. 4, pp. 414–416, 2004.

W.-W. Huang, S.-C. Tsai, S.-F. Peng et al., “Kaempferol induces autophagy through AMPK and AKT signaling molecules and causes G2/M arrest via downregulation of CDK1/cyclin B in SK-HEP-1 human hepatic cancer cells,” International Journal of Oncology, vol. 42, no. 6, pp. 2069–2077, 2013.

C. Preisinger, R. Körner, M. Wind, W. D. Lehmann, R. Kopajtich, and F. A. Barr, “Plk1 docking to GRASP65 phosphorylated by Cdk1 suggests a mechanism for Golgi checkpoint signalling,” The EMBO Journal, vol. 24, no. 4, pp. 753–765, 2005.

X. Q. Wang, C. M. Lo, L. Chen, E. S.-W. Ngan, A. Xu, and R. Y. Poon, “Cdk1-PDK1-PI3K/Akt signaling pathway regulates embryonic and induced pluripotency,” Cell Death & Differentiation, vol. 24, no. 1, pp. 38–48, 2017.

L. L. Wu, B. Qu, Y. T. Qin et al., “Identification of cyclin-dependent kinase 1 as a novel regulator for controlling type I Interferon signaling in systemic lupus erythematosus,” Arthritis & Rheumatology, vol. 67, 2015.

W. Gan, H. Zhao, T. G. Li, K. Liu, and J. Huang, “CDK1 interacts with iASPP to regulate colorectal cancer cell proliferation through p53 pathway,” Oncotarget, vol. 8, no. 42, pp. 71618–71629, 2017.

R. Visintin, S. Prinz, and A. Amon, “CDC20 and CDH1: a family of substrate-specific activators of APC-dependent proteolysis,” Science, vol. 278, no. 5337, pp. 460–463, 1997.

E. R. Kramer, N. Scheuringer, A. V. Podtelejnikov, M. Mann, and J.-M. Peters, “Mitotic regulation of the APC activator proteins CDC20 and CDH1,” Molecular Biology of the Cell, vol. 11, no. 5, pp. 1555–1569, 2000.

H. Yu, “Cdc20: a WD40 activator for a cell cycle degradation machine,” Molecular Cell, vol. 27, no. 1, pp. 3–16, 2007.

Y. Yang, A. H. Kim, T. Yamada et al., “A Cdc20-APC ubiquitin signaling pathway regulates presynaptic differentiation,” Science, vol. 326, no. 5952, pp. 575–578, 2009.

M. V. Hadjihannas, D. B. Bernkopf, M. Brückner, and J. Behrens, “Cell cycle control of Wnt/β-catenin signalling by conductin/axin2 through CDC20,” EMBO Reports, vol. 13, no. 4, pp. 347–354, 2012.

A.-H. Yen and J.-L. Yang, “Cdc20 proteolysis requires p38 MAPK signaling and Cdh1-independent APC/C ubiquitination during spindle assembly checkpoint activation by cadmium,” Journal of Cellular Physiology, vol. 223, no. 2, pp. 327–334, 2010.

Y. Ding, C. G. Hubert, J. Herman et al., “Cancer-Specific requirement for BUB1B/BUBR1 in human brain tumor isolates and genetically transformed cells,” Cancer Discovery, vol. 3, no. 2, pp. 198–211, 2013.

Y.-X. Shi, T. Zhu, T. Zou et al., “Prognostic and predictive values of CDK1 and MAD2L1 in lung adenocarcinoma,” Oncotarget, vol. 7, no. 51, pp. 85235–85243, 2016.

T. Song, Y. Zheng, Y. Wang et al., “Specific interaction of KIF11 with ZBP1 regulates the transport of -actin mRNA and cell motility,” Journal of Cell Science, vol. 128, no. 5, pp. 1001–1010, 2015.

F. Wang and S. L. Lin, “Knockdown of kinesin KIF11 abrogates directed migration in response to epidermal growth factor-mediated chemotaxis,” Biochemical and Biophysical Research Communications, vol. 452, no. 3, pp. 642–648, 2014.

Z. T. Swider, R. K. Ng, R. Varadarajan, C. J. Fagerstrom, and N. M. Rusan, “Fascetto (PRC1) interacting protein (FIP) ensures proper cytokinesis and ploidy,” Molecular Biology of the Cell, vol. 30, no. 8, pp. 992–1007, 2019.

C. Li, M. Lin, and J. Liu, “Identification of PRC1 as the p53 target gene uncovers a novel function of p53 in the regulation of cytokinesis,” Oncogene, vol. 23, no. 58, pp. 9336–9347, 2004.

H. W. Chaudhry, N. H. Dashoush, H. Tang et al., “Cyclin A2 mediates cardiomyocyte mitosis in the postmitotic myocardium,” Journal of Biological Chemistry, vol. 279, no. 34, pp. 35858–35866, 2004.

F. Yang, J. Gong, G. Wang, P. Chen, L. Yang, and Z. Wang, “Waltonitone inhibits proliferation of hepatoma cells and tumorigenesis via FXR-miR-22-CCNA2 signaling pathway,” Oncotarget, vol. 7, no. 46, pp. 75165–75175, 2016.

T. Gao, Y. Han, L. Yu, S. Ao, Z. Li, and J. Ji, “CCNA2 is a prognostic biomarker for ER+ breast cancer and tamoxifen resistance,” PLoS One, vol. 9, no. 3, Article ID e91771, 2014.

S. Dong, F. Huang, H. Zhang, and Q. Chen, “Overexpression of BUB1B, CCNA2, CDC20, and CDK1 in tumor tissues predicts poor survival in pancreatic ductal adenocarcinoma,” Bioscience Reports, vol. 39, no. 2, Article ID BSR20182306, 2019.

Y. Gan, Y. Li, T. Li, G. Shu, and G. Yin, “CCNA2 acts as a novel biomarker in regulating the growth and apoptosis of colorectal cancer,” Cancer Management and Research, vol. 10, pp. 5113–5124, 2018.

R. Broderick, J. Nieminuszczy, A. N. Blackford et al., “TOPBP1 recruits TOP2A to ultra-fine anaphase bridges to aid in their resolution,” Nature Communications, vol. 6, no. 1, p. 6572, 2015.

N. Wong, W. Yeo, W.-L. Wong et al., “TOP2A overexpression in hepatocellular carcinoma correlates with early age onset, shorter patients survival and chemoresistance,” International Journal of Cancer, vol. 124, no. 3, pp. 644–652, 2009.

J. Sun, Y. Du, Q. Song et al., “E2F is required for STAT3-mediated upregulation of cyclin B1 and Cdc2 expressions and contributes to G2-M phase transition,” Acta Biochimica et Biophysica Sinica, vol. 51, no. 3, pp. 313–322, 2019.

J. Gu, X. Liu, J. Li et al., “MicroRNA-144 inhibits cell proliferation, migration and invasion in human hepatocellular carcinoma by targeting CCNB1,” Cancer Cell International, vol. 19, no. 1, p. 15, 2019.