Microarray Analysis For Expression Profiles of lncRNAs and circRNAs in Rat Liver after Brain-Dead Donor Liver Transplantation

Chen, Sanyang; Fang, Hongbo; Li, Jie; Shi, Jihua; Zhang, Jiakai; Wen, Peihao; Wang, Zhihui; Yang, Han; Cao, Shengli; Zhang, Huapeng; Tang, Hongwei; Guo, Wenzhi; Zhang, Shuijun

doi:https://doi.org/10.1155/2019/5604843

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 5604843 | https://doi.org/10.1155/2019/5604843

Microarray Analysis For Expression Profiles of lncRNAs and circRNAs in Rat Liver after Brain-Dead Donor Liver Transplantation

Sanyang Chen,^1,2,3,4Hongbo Fang,¹Jie Li,¹Jihua Shi,^1,2,3,4Jiakai Zhang,¹Peihao Wen,¹Zhihui Wang,¹Han Yang,¹Shengli Cao,¹Huapeng Zhang,^1,2,3,4Hongwei Tang,^1,2,3,4and Wenzhi Guo^1,2,3,4 et al.

Academic Editor: Shinichi Aishima

Received05 Jan 2019

Revised01 Jun 2019

Accepted18 Jul 2019

Published07 Nov 2019

Abstract

The mechanisms underlying severe liver injury after brain-dead (BD) donor liver transplantation (BDDLT) remain unclear. In this study, we aimed to explore the roles of lncRNAs and circRNAs in liver injury after BDDLT. Rat liver injury was detected in the sham, BD, control, and BDDLT groups. We examined the expression profiles of lncRNAs and circRNAs in the livers of the BDDLT and control group using microarray analysis. The main functions of the differentially expressed genes were analyzed by gene ontology (GO) and KEGG pathway enrichment analysis. In addition, we used bioinformatic analyses to construct related expression networks. Liver injury was aggravated in the BD and BDDLT groups. We found various mRNAs, lncRNAs, and circRNAs that were differentially expressed in the BDDLT group compared with those in the control group. Coding-noncoding gene co-expression (CNC) network analysis showed that expression of the lncRNA LOC102553657 was associated with that of the apoptosis-related genes including HMOX1 and ATF3. Furthermore, competing endogenous RNAs (ceRNAs) network analysis revealed that the lncRNA LOC103692832 and rno_circRNA_007609 were ceRNAs of rno-miR-135a-5p targeting Atf3, Per2, and Mras. These results suggest that lncRNAs and circRNAs play important roles in the pathogenesis and development of liver injury during BDDLT.

1. Introduction

In current practice, most liver grafts are obtained from brain-dead (BD) donors [1]. Brain death refers to the irreversible and permanent loss of all function throughout the entire brain, including the brainstem [2]. Previous clinical and experimental studies have shown that a series of complex pathophysiological changes occur during brain death, including hemodynamic instability, increased leukocyte infiltration, hormonal imbalance, induction of cytokines and chemokines, and increased expression of heat shock proteins. These changes reduce the tolerance of liver grafts to the further cold preservation/reperfusion injury and aggravate liver injury after liver transplantation (LT), thus reducing the graft survival rate [3–6]. Despite an increasing number of studies on brain death, to date, there have been few in-depth studies on its mechanism, especially at the genetic level. Furthermore, no precise biomarker has been identified which can predict liver injury following BD donor LT (BDDLT).

Noncoding RNAs are RNA molecules that normally do not encode proteins, but instead functionally regulate protein expression [7, 8]. Noncoding RNAs mainly include long noncoding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs) [9]. lncRNAs are transcripts longer than 200 nucleotides which, because of their stability, are ideal biomarkers for disease research. New evidence suggests that lncRNAs play a key role in regulating many cellular biological processes, including cell proliferation, differentiation, apoptosis, migration and invasion, and the cell cycle [10–12]. In recent years, their potential roles in liver injury have been recognized. The lncRNA Gm2199 can promote hepatocyte regeneration and provides protection against liver injury, while the lncRNAs HOTAIR, HULC, and uc003wbd are highly expressed in patients with liver cancer or hepatitis B [13–16]. Recent studies have shown that circRNAs play a key role in the development of human diseases. CircRNAs are a type of RNA formed from a single-stranded, covalently closed loop, which serve mainly as a “sponge” for miRNAs [17, 18]. Additionally, circRNAs can compete with pre-mRNA splicing or affect gene transcription by their association with phosphorylated Pol II [19, 20]. Recently, an increasing number of studies have focused on the role of circRNAs in liver injury [21, 22].

Although Catania et al. investigated changes in mRNA expression during brain death [23]. However, to date, no study has explored the mechanisms by which noncoding RNAs participate in liver injury after BDDLT. Therefore, to investigate the role of noncoding RNAs in liver injury after BDDLT, we measured the differential expression of lncRNAs, circRNAs, and mRNAs in liver tissues after BDDLT, using microarray analysis. Significant differences in expression levels of lncRNAs and mRNAs were further verified by quantitative real-time polymerase chain reaction (qRT-RCR) analysis. In addition, we analyzed related pathways and gene ontology (GO), and furthermore performed network analyses of coding-noncoding gene CNC (CNC), competing endogenous RNAs (ceRNAs), lncRNAs’ adjacent coding genes, and TF-lncRNAs enriched with mRNAs using bioinformatic methods. Our study aimed to elucidate the mechanisms of liver injury after BDDLT to help identify biomarkers of liver injury and new therapeutic targets.

2. Materials and Methods

2.1. Animals and Establishment of BD, Orthotopic LT Model

Male Lewis rats weighing 250–280 g were purchased from Beijing Vital River Laboratory Animal Technology Co. (Beijing, China). All of the experimental procedures were approved by the Animal Ethics Committee of Zhengzhou University.

Brain death was induced by increasing intracranial pressure in a modified, slow, and intermittent fashion [24]. A 3F Fogarty catheter was connected to a pressure pump and then inserted into the epidural space through a hole in the skull. An additional catheter was inserted into the saphenous artery to record the arterial blood pressure, and a blood pressure above 80 mmHg was considered normal. The liver was procured 6 h after brain death induction according to the clinical practice [25, 26].

Livers from BD and nonBD donors were perfused with the University of Wisconsin (UW) preservation solution. The rat models of LT were established using the two-cuff technique without liver artery reconstruction [27]. The anhepatic phase lasted 17–20 min.

Experimental animals were randomized into four groups: sham-operated rats underwent only drilling of the skull and laparotomy, without any other procedure performed; in the BD group, brain death was maintained for 6 h after induction; in the control group, the sham-operated rats were used as donors for orthotopic LT; in the BDDLT group, grafts for orthotopic LT were obtained 6 h after brain death. Serum and tissue samples were obtained before liver harvesting and 12 h after reperfusion.

2.2. Biochemistry and Enzyme-Linked Immunosorbent Assay (ELISA)

Serum alanine aminotransferase (ALT) and aspartate aminotransferase (AST) were measured using a fully automated multi-functional serum biochemical analyzer (7600-020, Hitachi, Japan). Malondialdehyde (MDA) and superoxide dismutase (SOD) activity in the liver tissue were detected using an ELISA kit (Jiancheng, Nanjing, China). Liver tissue myeloperoxidase (MPO), liver tissue 3-nitrotyrosine, bile lactate dehydrogenase (LDH), and gamma-glutamyl transferase (GGT) were detected using a different ELISA kit (CUSABIO, Wuhan, China).

2.3. TUNEL Fluorescence Assay

From the TUNEL kit, appropriate volumes of reagent 1 (TdT) and reagent 2 (dUTP) were mixed at a ratio of 1 : 9 to cover the tissue. The mixture was added to the tissue with DAPI and then incubated at room temperature in the dark. The results were observed under a fluorescence microscope [24].

2.4. RNA Extraction and Quantitative Real-Time PCR Validation

RNA was extracted from tissue using TRIZOL (Invitrogen, America). Synthesis of cDNA was performed with a reverse transcription kit (TaKaRa, Japan) using the Gene Amp PCR System 9700 (Applied Biosystems, America), real-time PCR was performed using SYBR Green PCR Master Mix (TaKaRa) on the Gene Amp PCR System 9700, and the primers used were as described below. The primers of the randomly selected lncRNAs, mRNAs, and internal parameters of GAPDH are shown in Table 1.

2.5. Microarray Assay

Microarray analysis was performed by KangChen Bio-tech, Shanghai, PR China. Briefly, Three pairs of liver tissues were harvested 12 h after LT in the BDDLT group and the control group. Differential expression levels of lncRNAs, circRNAs, and mRNAs were measured using microarray analysis techniques [28]. Arraystar Rat LncRNA Microarray V2.0 was used. Sample labeling and array hybridization were performed according to the Agilent One-Color Microarray-Based Gene Expression Analysis protocol (Agilent Technology) with minor modifications. Briefly, mRNA was purified from total RNA after removal of rRNA (mRNA-ONLY™ Eukaryotic mRNA Isolation Kit, Epicentre). Then, each sample was amplified and transcribed into fluorescent cRNA along the entire length of the transcripts without 3′ bias utilizing a random priming method (Arraystar Flash RNA Labeling Kit, Arraystar). The hybridized arrays were washed, fixed, and scanned using the Agilent DNA Microarray Scanner (part number G2505C). Meanwhile, the Arraystar Rat circRNA Array V1 analysis of the liver samples was conducted. Total RNAs were digested with Rnase R (Epicentre, Inc.) to remove linear RNAs and enrich circular RNAs. Then, the enriched circular RNAs were amplified and transcribed into fluorescent cRNA utilizing a random priming method (Arraystar Super RNA Labeling Kit; Arraystar). The labeled cRNAs were hybridized onto the Arraystar Rat circRNA Array V1 (8 ×∉15 K, Arraystar). After having washed the slides, the arrays were scanned by the Agilent Scanner G2505C. Agilent Feature Extraction software (version 11.0.1.1) was used to analyze acquired array images. Statistically significant differences in the expression of lncRNAs, circRNAs, and mRNAs were evaluated in the BDDLT group and the control group based on -values and FDR filtering. A fold change ≥ 2.0, , and FDR < 0.05 were considered statistically significant. We deposited the microarray data and additional information to the GEO database with accession number GSE14520 and GSE119748, respectively.

2.6. Gene Ontology (GO) and KEGG Pathway Analysis

Gene ontology (GO) analysis provides a controlled vocabulary to describe gene and gene product attributes in any organism. We use topGO to analyze differential genes in order to infer their molecular functions (http://www.geneontology.org). The 10 items with the smallest -values were selected for graphical display; in the case that fewer than 10 items met the statistical criteria, all the items were displayed. We performed KEGG pathway analysis to harvest pathway clusters covering our knowledge on the molecular interaction and reaction networks in differentially regulated gene profiling. Again, the 10 items with the smallest -values were selected for graphical display or all the items if fewer than 10 met the statistical criteria. We also used differentially expressed genes for pathway analysis to infer the pathways they participate in. Pathway analysis is a functional analysis mapping genes to KEGG pathways. The -value (EASE-score, Fisher--value or Hypergeometric--value) denotes the significance of the Pathway correlated to the conditions. Lower the -value, more significant is the Pathway. (The recommend -value cut-off is 0.05.) Significant pathways contain some genes that differently, pathway-gene network was constructed by the relationship between pathway and gene.

2.7. CNC and CeRNA Network Analysis

CNC is a method of correlation analysis between noncoding RNA and mRNA co-expression data, which is based on the Pearson correlation coefficient (PCC) between coding genes and noncoding transcripts according to their expression levels. The criteria for significant effects were set to PCC ≥ 0.90, , and FDR <0.01. The normalized signal intensities of specific mRNA and lncRNA expression levels were used to construct coexpression network. According to the expression data of coding and noncoding RNA, we calculated the co-expression relationship (Pearson correlation coefficient, PCC) between coding and noncoding RNA to construct the network. By analyzing the functions of these known RNA molecules, we can associate unknown RNA molecules with specific signaling pathways or diseases to predict their functions and mechanisms. The network was visualized using Cytoscape software (version 2.8.3) [29]. By analyzing the functions of these known RNA molecules, we can associate unknown RNA molecules with specific signaling pathways or diseases to predict their functions and mechanisms.

As is known to us, the ceRNA regulatory mechanism is very important between the mRNA and ncRNAs, including miRNA, lncRNA, and circRNA. Possible miRNA response elements were searched for based on the lncRNA, circRNA, and mRNA sequences. We predicted target mRNAs and lncRNAs (circRNAs) of these miRNAs based on TargetScan 7.2 and compared these target genes to the differentially expressed genes that were identified in our array results. Then the overlapping of the same miRNA seed sequence binding site both on the lncRNAs/circRNAs, and the mRNA predicted lncRNA/circRNA-miRNA-mRNA interaction.

2.8. Cis- and Trans-Regulation Prediction

lncRNAs whose expression showed significant correlations with the cis- and trans-regulation of mRNAs were processed. For cis-regulation, previous study defined that a cis-regulator is the one that exerts its function on a neighbouring gene located at the same chromosome, so we search for the adjacent genes less than 300 kb upstream or downstream of lncRNAs may serve as potential targets for cis-regulation of lncRNA [28]. For trans-regulation, we were more interested in the way lncRNAs exert their effects via transcription factors (TFs). The TFs were predicted by the TRANSFAC8.3 database based on promoter regions of lncRNA. LncRNA binding transcription factors were predicted by scoring and E value <0.01. Therefore, we enriched those mRNAs co-expressed with lncRNAs that significantly overlapped with the target genes of a given TF and constructed an lncRNA-TF-mRNA network. The enrichment and connectivity was based on Position Frequency Matrix (PFM) performed as described previously [30, 31].

2.9. Statistical Analysis

All statistical analyses were performed in SPSS 21.0 and GraphPad Prism software. Data are expressed as mean ± standard deviation (means ± SD). Statistical significance was defined as value less than 0.05 using two-tailed unpaired Student’s t-test.

3. Results

3.1. Brain Death Promotes Early Liver Graft Injury in Rats

Brain death is characterized by a rapid increase in blood pressure (autonomic storm), followed by a decrease, usually to baseline level or below, and a flat-lined EEG (Supplementary Figure 1). We first measured the serum levels of transaminases (AST and ALT). Compared with those in the control group, AST and ALT levels were significantly increased in the BD group and the BDDLT group. Based on the ELISA results, the BD and BDDLT groups showed greater aggravated oxidative stress and inflammatory injury than the control group, characterized by significantly increased levels of MDA, MPO, and 3-notrotyrosine and decreased SOD levels in the BDDLT group (Figure 1(a)). In addition, liver tissue cell apoptosis was detected by TUNEL staining; compared with the control group, the BDDLT group had increased cell apoptosis (Figure 1(b).

(a)

(b)

Figure 1

Biochemical and ELISA measurements of liver damage markers; TUNEL staining of liver tissue. (a) ALT, AST, MDA, 3-nitrotyrosine, and MPO levels were significantly increased, and the SOD level was significantly decreased in the BD and BDDLT groups. All data are shown as mean ± SD. The levels of statistical significance are indicated as ^∗, ^∗^∗ statistically significant. For statistical analysis, a two-tailed Student’s t test was used. (b) TUNEL staining of liver tissue showed more apoptotic cells in the BDDLT group (scale bar: 50 µm).

3.2. Microarray Analysis Reveals Differentially Expressed lncRNA and circRNA Profiles

Microarray probes successfully detected thousands of transcripts in the liver tissue in the BDDLT group and in the control group. Compared with those in the control group, 228 differentially expressed lncRNAs and 330 differentially expressed mRNAs were detected in three liver tissues in the BDDLT group. Among these differentially expressed transcripts, there were 181 up-regulated lncRNAs, among which LOC102553824 showed the largest difference (fold change: 7.86); there were 47 down-regulated lncRNAs, among which AABR06018038.2 showed the largest difference (fold change: 4.32). There were 164 and 166 up- and down-regulated mRNA transcripts, respectively (fold change ≥ 2.0, , and FDR < 0.05). In addition, a total of 5,501 up-regulated and 5,329 down-regulated circRNAs were detected in the liver tissue chips from both the BDDLT and control groups. Coding gene expression profiling revealed 19 mRNAs with a more than 5-fold difference, among which 12 were up-regulated, with the largest fold change being that for Acpp (fold change: 12.22), and 7 were down-regulated, with the largest fold change being that for Serpina7 (fold change: 7.71). Clustering analysis showed that the expression patterns of lncRNAs, mRNAs, and circRNAs significantly differed between the two groups. Thus, expression levels of the identified lncRNAs, cirRNAs, and mRNAs in the liver tissue after BDDLT were significantly different from those in the control group (Figure 2(a)).

(a)

(b)

(c)

(d)

Figure 2

Heat maps showing expression profiles of lncRNAs, circRNAs, mRNAs, and identification of differentially expressed lncRNAs in the liver tissues of the control and BDDLT group. (a) Heat maps show normalized expression values of significantly changed lncRNAs, circRNAs, and mRNAs with a fold change ≥ 2.0, , and FDR < 0.05. C: control group tissues. B: BDDLT group tissues. (b) Circos images show the distribution of lncRNA in rat chromosomes, the outermost layer of Circos is the distribution of rat chromosomes, the black and white bands are the cytobands of chromosomes, the middle two circles are the up-down distribution of lncRNA, the intermediate outer ring represents the difference of all detected lncRNAs in the chip, and the intermediate inner ring represents the difference of different lncRNAs (fold change is greater than 2.0, ). (c) Types and counts of differentially regulated lncRNAs detected by microarray (fold change ≥ 2.0, , and FDR < 0.05). The lncRNAs are classified into 6 types according to the relationship and genomic loci with their associated coding genes. (d) The Venn diagram represents the overlapping of relationships and the numbers indicate the lncRNA counts.

These lncRNAs are widely distributed on all chromosomes including the X chromosome and the Y chromosome. Circos images show the distribution of lncRNA in rat chromosomes. The outermost layer of Circos shows the distribution of rat chromosomes. The black and white bands represent the cytobands of the chromosomes. The two circles in the middle represent the up-down distribution of lncRNAs. The more intense the red color, the higher the up-down multiple, and the more intense the blue color, the higher the down-down multiple. The intermediate outer ring represents the difference in all detected lncRNAs in the chip, and the intermediate inner ring represents the difference in different lncRNAs (fold change is greater than 2.0, ) (Figure 2(b)). The abnormally expressed lncRNAs were divided into six categories according to their relationships with protein-encoding genes: exon-overlapping (10.9%), intron-overlapping (8.7%), natural antisense (3.8%), and intron antisense (9.3%), and there were a number of crossover events among these four categories, with 60.7% being intergenic and 6.6% being bidirectional. Furthermore, intergenic lncRNAs were the largest category among all of the differentially expressed lncRNAs (60.7%), with both the up-regulated intergenic lncRNAs (45.9%) and down-regulated intergenic lncRNAs (14.8%) accounting for large proportions (Figure 2(c)). The Venn diagram presents the overlapping of relationships and the numbers indicate the lncRNA counts (Figure 2(d)).

3.3. GO and KEGG Pathway Analysis

GO enrichment analysis of significantly differentially expressed mRNAs can reveal the role of differentially regulated lncRNAs. Up-regulated mRNAs have been found to be primarily involved in single-organism metabolic processes and activities of the germ cell nucleus and hydrolase. Furthermore, down-regulated transcripts mainly participate in drug metabolism, extracellular space, and immunoglobulin receptor binding (Figure 3(a)). Our data revealed ten pathways associated with the up- and down-regulation of mRNAs. KEGG enrichment analysis of the up-regulated protein-coding genes showed that the signaling pathway responsible for fatty acid elongation was the most significantly enriched. In contrast, KEGG analysis of the down-regulated transcripts revealed that the PI3K-Akt signaling pathway was the most significantly enriched pathway (Figure 3(b)). According to the relationship between enrichment pathway and mRNA, the corresponding pathway-gene network map was constructed (Figure 3(c)).

(a)

(b)

(c)

3.4. Verification of Differentially Expressed lncRNAs and mRNAs

In our results, mRNAs reported to be involved in hepatic ischemia-reperfusion injury, and their corresponding LncRNA by bioinformatics speculation were selected for qRT-PCR analysis to verify the microarray results, using 10 pairs of liver tissues obtained after LT in the BDDLT group and the control group. The analysis revealed that expression of LOC102555648, LOC103692056, LOC103692832, Ablim3-XR_596973, Ablim3-XR_596974, and LOC102554711 decreased in the liver tissue of the BDDLT group, consistent with the microarray results (Figure 4(a)). Similarly, consistent with the microarray results, eight mRNAs showed altered expression compared with those in the control group, with Per2, Inhba, and Rb1 up-regulated, but Angptl4, Atf3, Hmox1, Ndrg1, and Trib3 down-regulated (Figure 4(b)).

(a)

(b)

3.5. CNC and Functional Prediction of lncRNAs/mRNAs

We selected seven coding genes with significant differential expression to construct the CNC network based on the correlation degree. These mRNAs are involved in many biological processes such as the cell cycle, apoptosis, epithelial-mesenchymal transition (EMT), and angiogenesis. The network revealed that the up-regulated lncRNA AABR06081886.1 was negatively correlated with Atf3, Hmox1, Angptl4, and Trib3, whereas the down-regulated lncRNA LOC100911923 was positively correlated with these four protein-coding genes, which are mainly involved in cell apoptosis (Figure 5).

3.6. Prediction of Cis- and Trans-Regulatory Functions of lncRNAs (Adjacent Coding Genes of lncRNAs and lncRNA-TF-mRNA)

Based on the coexpression, we further explored how the abnormally regulated lncRNAs could achieve cis- or trans-regulation of mRNAs. Ten differentially expressed lncRNAs were selected to conduct searches for their adjacent coding genes. The coexpressed protein-coding genes were defined as cis-regulated genes with one differentially expressed lncRNA within 300 kb on the same chromosome. The number of adjacent coding genes differs among lncRNAs. For example, LOC103692721 has up to three adjacent coding genes, whereas other lncRNAs have only one target gene. The adjacent gene pattern showed that the down-regulated IGH-6 and ABR06046430.3 were cis-regulated by lncRNA LOC103692721 and involved in liver injury after LT from BD donors. The up-regulated lncRNAs-TFs network consisted of 15 TFs and 95 correlated lncRNAs, and they were connected by 172 edges. The down-regulated lncRNAs-TFs network consisted of 13 TFs and 26 correlated lncRNAs, and they were connected by 43 edges (Figures 6(a) and 6(b)).

(a)

(b)

Figure 6

Cis regulation of lncRNAs to nearby coding genes and the network of enrichment transcription factors and query lncRNAs. (a) lncRNAs and their potential cis regulated nearby genes are shown in the network and the distances between lncRNAs and their cis regulated genes are presented. (b) According to the prediction of transcription factor regulation in the promoter region of lncRNA, the transcription factors of rats in the TRANSFAC8.3 database were predicted. The up-regulated lncRNAs-TFs network consists of 15 TFs and 95 correlated lncRNAs, and they are connected by 172 edges. The down-regulated lncRNAs-TFs network consists of 13 TFs and 26 correlated lncRNAs, and they are connected by 43 edges.

Assuming that lncRNAs have transregulatory functions, we analyzed the mRNAs coexpressed with these lncRNAs and mRNAs regulated by TFs. We found that the up-regulated-lncRNA-TF-gene network consisted of 83 lncRNAs, 15 TFs, and 6 correlated genes. The down-regulated-lncRNA-TF-gene network consisted of 22 lncRNAs, 15 TFs, and 6 correlated genes. The lncRNA-TF-mRNA diagram showed that the lncRNA LOC102553657 trans-regulates HMOX1 via the HNF-3beta TF. Next, on the basis of TF-lncRNA binary analysis, we additionally introduced mRNAs to construct the TF-lncRNA-mRNA ternary network. Most lncRNAs, including AP-2, AR, HNF-1, Crx, Sp3, CBF (2), HNF-3alpha, Nkx2-1, and POU1F1a, were found to be involved in TF regulation (Figures 7(a) and 7(b)).

(a)

(b)

3.7. Construction of the ceRNA Network

The ceRNA network was constructed based on microarray results from liver tissues after LT in both the BDDLT group and the control group. Six lncRNAs and 11 circRNAs with significant differential expression were selected. Some of these shared miRNA response elements (MREs); for example, lncRNAs LOC103692832 and rno_circRNA_007609 are ceRNAs of rno-miR-135a-5p targeting Atf3, Per2, and Mras, whereas lncRNAs Ablim3 and rno_circRNA_013693 are ceRNAs of rno-miR-322-5p targeting Per2, Mapkapk3, and Trim59 (Figure 8).

4. Discussion

Our study corroborates previous research demonstrating that the complex pathophysiological changes occurring during brain death have a significant impact on graft function and survival [32, 33]. The “autonomous storm” during brain death is considered the main cause of hemodynamic instability, which can lead to insufficient organ perfusion, activate proinflammatory factors, and thus promote liver injury soon after transplantation [34, 35]. We found that the levels of transaminases were elevated after BDDLT, along with activation of inflammatory factors and an increased number of apoptotic cells. However, the mechanisms underlying liver injury after BDDLT remain unclear.

The traditional view of gene regulation focuses on protein-coding genes; however, this view has changed since the discovery of many noncoding genes such as lncRNAs, circRNAs, and microRNAs [9]. The roles that lncRNAs can play in the diagnosis and treatment of liver diseases have been recognized [36]. Abnormal lncRNA expression levels have been found in different types of liver injury, suggesting that lncRNAs are involved in the pathogenesis of liver injury [37, 38]. MEG3 and GAS5 can inhibit liver fibrosis, while TUG1 can protect liver grafts during cold preservation [39–41]. A potential relationship between liver injury and circRNA has been previously established, suggesting that circRNA inhibits hepatic steatosis [42]. A large number of differentially expressed circRNAs were also found in hepatic ischemia-reperfusion injury models [21]. In this study, we analyzed for the first time the changes in expression of lncRNAs, circRNAs, and mRNAs in liver tissues after BDDLT.

We obtained the expression profiles of lncRNAs and circRNAs by microarray analysis. A large number of differentially expressed lncRNAs and circRNAs were found in liver tissues after BDDLT, with about 60% of the down-regulated lncRNAs being intergenic. Six down-regulated lncRNAs and eight differentially expressed mRNAs were further verified by qRT-PCR. Overall, the qRT-PCR verification results were consistent with the microarray results. The differentially expressed genes were subsequently classified according to the data using thermograms and cluster maps. We also found that the differentially expressed lncRNAs were distributed across all the chromosomes, suggesting that every chromosome in the liver tissue, including the X and Y chromosome, was abnormal after BDDLT. We found that Hmox1, Ndrg1, and Atf3 were significantly down-regulated in liver tissues after BDDLT; thus, these genes may be involved in inducing apoptosis after BDDLT to cause liver injury.

KEGG pathway analysis of differentially expressed mRNAs revealed that multiple pathways (primarily the fatty acid elongation pathway and the PI3K-Akt signaling pathway) were involved in liver injury after BDDLT. GO analysis revealed that the main mechanisms of liver injury after BDDLT included single-organism metabolism, drug metabolism, and immunoglobulin receptor binding. At present, the functions of most lncRNAs are unknown, and constructing a CNC network is one method for predicting lncRNA functions [43, 44]. We found that the expression levels of many lncRNAs were significantly associated with protein-encoding genes through correlation of expression levels. Therefore, we established a CNC network to predict the relationship between lncRNAs and mRNAs. The network revealed that the up-regulated lncRNA AABR06081886.1 was negatively correlated with Atf3, Hmox1, Angptl4, and Trib3, whereas the down-regulated lncRNA LOC100911923 was positively correlated with these four protein-coding genes, which are mainly involved in cell apoptosis.

It has been reported that lncRNA can affect the expression of adjacent or distal genes via cis-regulation and trans-regulation [45–47]. We studied the relationships between ten significantly differentially expressed lncRNAs and their adjacent genes, in an attempt to provide a new theoretical mechanism for liver injury after BDDLT. We found that the number of adjacent coding genes differed among lncRNAs; for example, LOC103692721 had up to three adjacent coding genes, whereas other lncRNAs had only one target gene. In addition, the adjacent gene pattern showed that the down-regulated gene IGH-6 and ABR06046430.3 were cis-regulated by lncRNA LOC103692721 and involved in liver injury after BDDLT.

Furthermore, based on our results, we constructed the lncRNA-TF network. Transcriptional regulation by lncRNAs can be predicted from the regulation of TFs, which they express [48–50]. The combination of TFs and cis-acting elements can also regulate the expression levels of target genes in the promoter region [51, 52]. It has also been reported that TFs participate in the occurrence and development of liver injury [53, 54]. Therefore, a comprehensive analysis of TFs and differentially coexpressed genes can improve our understanding of liver injury after BDDLT. A number of lncRNAs may be involved in specific pathways for TF regulation. Therefore, assuming that lncRNAs regulate transcription, we analyzed mRNAs coexpressed with lncRNA and TF-regulated mRNAs. Using and FDR <0.01 as the statistical criteria, we found that each lncRNA was associated with between one to over a dozen TFs, and that each lncRNA-TF pair resulted from the enrichment of several genes, providing key data for subsequent studies. In the liver tissues from the BDDLT group, we found that the up-regulated lncRNAs-TFs network consists of 15 TFs and 95 correlated lncRNAs, and they are connected by 172 edges. The down-regulated lncRNAs-TFs network consists of 13 TFs and 26 correlated lncRNAs, and they are connected by 43 edges. Furthermore, the lncRNA-TF-mRNA diagram showed that lncRNA LOC102553657 trans-regulated HMOX1 via HNF-3beta TF. Thus, we speculate that these TFs participate in liver injury after BDDLT by regulating the transcription of lncRNAs and mRNAs. However, the relationships between lncRNAs and TFs needs further investigation.

Pandolfi et al. introduced a hypothesis called the ceRNA mechanism, which proposed that transcripts such as mRNAs, pseudogenes, and lncRNAs can serve as natural miRNA sponges by competitively binding to miRNA response elements (MREs) to suppress their expression and function [55]. Indeed, ceRNAs have important influences on gene expression at the post-transcriptional level, and several studies have reported that lncRNA is involved in the pathogenesis of liver injury through ceRNA [56–59]. However, the roles of lncRNA- and circRNA-related ceRNAs in liver injury have not been studied. In the current study, we constructed for the first time the ceRNA network with lncRNA-miRNA-circRNA-mRNA associations based on microarray results. Our results reveal a specific ceRNA network underlying liver injury after BDDLT. We found that the lncRNAs LOC103692832 and rno_circRNA_007609 were ceRNAs of rno-miR-135a-5p targeting Atf3, Per2, and Mras, whereas the lncRNAs Ablim3 and rno_circRNA_013693 were ceRNAs of rno-miR-322-5p targeting Per2, Mapkapk3, and Trim59. Thus, lncRNAs and circRNAs harbor MREs and thus play key roles in liver injury after BDDLT.

5. Conclusions

We found that brain death promotes liver injury after LT. Furthermore, we identified a number of specific abnormally regulated lncRNAs and circRNAs, which could be promising markers of the pathogenesis of liver injury after BDDLT and provide new directions for predicting and diagnosing this condition. Our research paves the way for further studies regarding the roles of lncRNAs and circRNAs after BDDLT.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conficts of interest.

Authors’ Contributions

Sanyang Chen and Hongbo Fang have contributed equally to this work.

Acknowledgments

The authors gratefully acknowledge financial support from the National Natural Science Foundation of China (No. 81571947, 81671958, U1604282). Henan basic and frontier science and technology research project (No. 142300410039).

Supplementary Materials

Supplementary Figure 1: EEG (A) and MAP (B) 6 hours after brain death. (Supplementary Materials)

References

M. Jiménez-Castro, E. Negrete-Sánchez, A. Casillas-Ramírez et al., “The effect of cortisol in rat steatotic and non-steatotic liver transplantation from brain-dead donors,” Clinical Science, vol. 131, no. 8, pp. 733–746, 2017.
View at: Publisher Site | Google Scholar
D. Magnus, B. Wilfond, and A. Caplan, “Accepting brain death,” New England Journal of Medicine, vol. 370, no. 10, pp. 891–894, 2014.
View at: Publisher Site | Google Scholar
R. F. Vieira, A. C. Breithaupt-Faloppa, B. C. Matsubara et al., “17β-Estradiol protects against lung injuries after brain death in male rats,” The Journal of Heart and Lung Transplantation, vol. 37, no. 11, pp. 1381–1387, 2018.
View at: Publisher Site | Google Scholar
R. Carlessi, N. Lemos, A. Dias et al., “Exendin-4 attenuates brain death-induced liver damage in the rat,” Liver Transplantation, vol. 21, no. 11, pp. 1410–1418, 2015.
View at: Publisher Site | Google Scholar
R. Spindler, P. Schnuelle, L. Nickels et al., “N-octanoyl dopamine for donor treatment in a brain-death model of kidney and heart transplantation,” Transplantation, vol. 99, no. 5, pp. 935–941, 2015.
View at: Publisher Site | Google Scholar
M. Jiménez-Castro, N. Meroño, M. Mendes-Braz et al., “The effect of brain death in rat steatotic and non-steatotic liver transplantation with previous ischemic preconditioning,” Journal of Hepatology, vol. 62, no. 1, pp. 83–91, 2015.
View at: Publisher Site | Google Scholar
K. Panir, J. Schjenken, S. Robertson, and M. Hull, “Non-coding RNAs in endometriosis: a narrative review,” Human Reproduction Update, vol. 24, no. 4, pp. 497–515, 2018.
View at: Publisher Site | Google Scholar
C.-M. Wong, F. H.-C. Tsang, and I. O.-L. Ng, “Non-coding RNAs in hepatocellular carcinoma: molecular functions and pathological implications,” Nature Reviews Gastroenterology & Hepatology, vol. 15, no. 3, pp. 137–151, 2018.
View at: Publisher Site | Google Scholar
J. Beermann, M. Piccoli, J. Viereck, and T. Thum, “Non-coding RNAs in development and disease: background, mechanisms, and therapeutic approaches,” Physiological Reviews, vol. 96, no. 4, pp. 1297–1325, 2016.
View at: Publisher Site | Google Scholar
O. Wapinski and H. Chang, “Long noncoding RNAs and human disease,” Trends in Cell Biology, vol. 21, no. 6, pp. 354–361, 2011.
View at: Publisher Site | Google Scholar
F. P. Marchese, I. Raimondi, and M. Huarte, “The multidimensional mechanisms of long noncoding RNA function,” Genome Biology, vol. 18, no. 1, 206 pages, 2017.
View at: Publisher Site | Google Scholar
J. Evans, F. Feng, and A. Chinnaiyan, “The bright side of dark matter: lncRNAs in cancer,” Journal of Clinical Investigation, vol. 126, no. 8, pp. 2775–2782, 2016.
View at: Publisher Site | Google Scholar
Q. Gao, Y. Gu, Y. 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, 602 pages, 2018.
View at: Publisher Site | Google Scholar
Y. Du, G. Kong, X. You et al., “Elevation of highly up-regulated in liver cancer (HULC) by hepatitis B virus X protein promotes hepatoma cell proliferation via down-regulating p18,” Journal of Biological Chemistry, vol. 287, no. 31, pp. 26302–26311, 2012.
View at: Publisher Site | Google Scholar
F. Yang, L. Zhang, X. Huo et al., “Long noncoding RNA high expression in hepatocellular carcinoma facilitates tumor growth through enhancer of zeste homolog 2 in humans,” Hepatology, vol. 54, no. 5, pp. 1679–1689, 2011.
View at: Publisher Site | Google Scholar
J. Lu, F. Xie, L. Geng, W. Shen, C. Sui, and J. Yang, “Investigation of serum lncRNA-uc003wbd and lncRNA-AF085935 expression profile in patients with hepatocellular carcinoma and HBV,” Tumour Biology, vol. 36, no. 5, pp. 3231–3236, 2015.
View at: Publisher Site | Google Scholar
B. Chen, W. Wei, X. Huang et al., “circEPSTI1 as a prognostic marker and mediator of triple-negative breast cancer progression,” Theranostics, vol. 8, no. 14, pp. 4003–4015, 2018.
View at: Publisher Site | Google Scholar
D. Han, J. Li, H. Wang et al., “Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression,” Hepatology, vol. 66, no. 4, pp. 1151–1164, 2017.
View at: Publisher Site | Google Scholar
Y. Zhang, X.-O. Zhang, T. Chen et al., “Circular intronic long noncoding RNAs,” Molecular Cell, vol. 51, no. 6, pp. 792–806, 2013.
View at: Publisher Site | Google Scholar
R. Ashwal-Fluss, M. Meyer, N. Pamudurti et al., “circRNA biogenesis competes with pre-mRNA splicing,” Molecular Cell, vol. 56, no. 1, pp. 55–66, 2014.
View at: Publisher Site | Google Scholar
Z. Ye, Q. Kong, J. Han, J. Deng, M. Wu, and H. Deng, “Circular RNAs are differentially expressed in liver ischemia/reperfusion injury model,” Journal of Cellular Biochemistry, vol. 119, no. 9, pp. 7397–7405, 2018.
View at: Publisher Site | Google Scholar
J. Guo, Y. Zhou, Y. Cheng et al., “Metformin-induced changes of the coding transcriptome and non-coding rnas in the livers of non-alcoholic fatty liver disease mice,” Cellular Physiology Biochemistry, vol. 45, no. 4, pp. 1487–1505, 2018.
View at: Publisher Site | Google Scholar
G. Colombo, S. Gatti, F. Turcatti et al., “Alteration in the transcriptional profile of livers from brain-dead organ donors,” Transplantation, vol. 82, no. 1, pp. 69–79, 2006.
View at: Publisher Site | Google Scholar
W. Guo, S. Cao, B. Yan et al., “Myocardial protective effects of a c-Jun N-terminal kinase inhibitor in rats with brain death,” Journal of Cellular and Molecular Medicine, vol. 20, no. 7, pp. 1214–1218, 2016.
View at: Publisher Site | Google Scholar
J. Van der Hoeven, S. Lindell, R. van Schilfgaarde et al., “Donor brain death reduces survival after transplantation in rat livers preserved for 20 hr,” Transplantation, vol. 72, no. 10, pp. 1632–1636, 2001.
View at: Publisher Site | Google Scholar
J. A. B. van der Hoeven, G. J. Ter Horst, G. Molema et al., “Effects of brain death and hemodynamic status on function and immunologic activation of the potential donor liver in the rat,” Annals of Surgery, vol. 232, no. 6, pp. 804–813, 2000.
View at: Publisher Site | Google Scholar
N. Kamada and R. Y. Calne, “Orthotopic liver transplantation in the rat. technique using cuff for portal vein anastomosis and biliary drainage,” Transplantation, vol. 28, no. 1, pp. 47–50, 1979.
View at: Publisher Site | Google Scholar
M. Huang, Z. Zhong, M. Lv, J. Shu, Q. Tian, and J. Chen, “Comprehensive analysis of differentially expressed profiles of lncRNAs and circRNAs with associated co-expression and ceRNA networks in bladder carcinoma,” Oncotarget, vol. 7, no. 30, pp. 47186–47200, 2016.
View at: Publisher Site | Google Scholar
X. Li, J. Ao, and J. Wu, “Systematic identification and comparison of expressed profiles of lncRNAs and circRNAs with associated co-expression and ceRNA networks in mouse germline stem cells,” Oncotarget, vol. 8, no. 16, pp. 26573–26590, 2017.
View at: Publisher Site | Google Scholar
M. Fu, G. Huang, Z. Zhang et al., “Expression profile of long noncoding RNAs in cartilage from knee osteoarthritis patients,” Osteoarthritis and Cartilage, vol. 23, no. 3, pp. 423–432, 2015.
View at: Publisher Site | Google Scholar
M. Guttman and J. Rinn, “Modular regulatory principles of large non-coding RNAs,” Nature, vol. 482, no. 7385, pp. 339–346, 2012.
View at: Publisher Site | Google Scholar
H. Fang, S. Zhang, W. Guo et al., “Cobalt protoporphyrin protects the liver against apoptosis in rats of brain death,” Clinics and Research in Hepatology and Gastroenterology, vol. 39, no. 4, pp. 475–481, 2015.
View at: Publisher Site | Google Scholar
S. Cao, T. Wang, B. Yan, Y. Lu, W. Guo, and S. Zhang, “Protective effects of SP600125 in brain death-induced liver injury,” Clinics and Research in Hepatology and Gastroenterology, vol. 38, no. 5, pp. 577–582, 2014.
View at: Publisher Site | Google Scholar
J. Pratschke, M. Wilhelm, M. Kusaka et al., “Brain death and its influence on donor organ quality and outcome after transplantation,” Transplantation, vol. 67, no. 3, pp. 343–348, 1999.
View at: Publisher Site | Google Scholar
M. Cornide-Petronio, E. Negrete-Sánchez, M. Mendes-Braz et al., “The effect of high-mobility group box 1 in rat steatotic and nonsteatotic liver transplantation from donors after brain death,” American Journal of Transplantation, vol. 16, no. 4, pp. 1148–1159, 2016.
View at: Publisher Site | Google Scholar
S. Roy, C. Trautwein, T. Luedde, and C. Roderburg, “A general overview on non-coding RNA-based diagnostic and therapeutic approaches for liver diseases,” Frontiers in Pharmacology, vol. 9, Article ID 805, 2018.
View at: Publisher Site | Google Scholar
X. Li, R. Liu, Z. Huang et al., “Cholangiocyte-derived exosomal long noncoding RNA H19 promotes cholestatic liver injury in mouse and humans,” Hepatology, vol. 68, no. 2, pp. 599–615, 2018.
View at: Publisher Site | Google Scholar
W. Yang, J. Chen, Y. Meng, Z. Chen, and J. Yang, “Novel targets for treating ischemia-reperfusion injury in the liver,” International Journal of Molecular Sciences, vol. 19, no. 5, p. 1302, 2018.
View at: Publisher Site | Google Scholar
F. Yu, J. Zheng, Y. Mao et al., “Long non-coding RNA growth arrest-specific transcript 5 (gas5) inhibits liver fibrogenesis through a mechanism of competing endogenous RNA,” Journal of Biological Chemistry, vol. 290, no. 47, pp. 28286–28298, 2015.
View at: Publisher Site | Google Scholar
L. Zhang, Z. Yang, J. Trottier, O. Barbier, and L. Wang, “Long noncoding RNA MEG3 induces cholestatic liver injury by interaction with PTBP1 to facilitate shp mRNA decay,” Hepatology, vol. 65, no. 2, pp. 604–615, 2017.
View at: Publisher Site | Google Scholar
S. Su, J. Liu, K. He et al., “Overexpression of the long noncoding RNA TUG1 protects against cold-induced injury of mouse livers by inhibiting apoptosis and inflammation,” The Febs Journal, vol. 283, no. 7, pp. 1261–1274, 2016.
View at: Publisher Site | Google Scholar
X.-Y. Guo, F. Sun, J.-N. Chen, Y.-Q. Wang, Q. Pan, and J.-G. Fan, “circRNA_0046366 inhibits hepatocellular steatosis by normalization of PPAR signaling,” World Journal of Gastroenterology, vol. 24, no. 3, pp. 323–337, 2018.
View at: Publisher Site | Google Scholar
Y. Li, C. Zeng, J. Hu et al., “Long non-coding RNA-SNHG7 acts as a target of miR-34a to increase GALNT7 level and regulate PI3K/Akt/mTOR pathway in colorectal cancer progression,” Journal of Hematology & Oncology, vol. 11, no. 1, 89 pages, 2018.
View at: Publisher Site | Google Scholar
Q. Liao, C. Liu, X. Yuan et al., “Large-scale prediction of long non-coding RNA functions in a coding-non-coding gene co-expression network,” Nucleic Acids Research, vol. 39, no. 9, pp. 3864–3878, 2011.
View at: Publisher Site | Google Scholar
S. Lykke-Andersen, B. K. Ardal, A. K. Hollensen, C. K. Damgaard, and T. H. Jensen, “Box C/D snoRNP autoregulation by a cis-acting snoRNA in the NOP56 pre-mRNA,” Molecular Cell, vol. 72, no. 1, pp. 99–111.e5, 2018.
View at: Publisher Site | Google Scholar
W. Xiong, Y. Jiang, Y. Ai et al., “Microarray analysis of long non-coding RNA expression profile associated with 5-fluorouracil-based chemoradiation resistance in colorectal cancer cells,” Asian Pacific Journal of Cancer Prevention, vol. 16, no. 8, pp. 3395–3402, 2015.
View at: Publisher Site | Google Scholar
I. Cajigas, A. Chakraborty, K. R. Swyter et al., “The Evf2 ultraconserved enhancer lncRNA functionally and spatially organizes megabase distant genes in the developing forebrain,” Molecular Cell, vol. 71, no. 6, pp. 956–972.e9, 2018.
View at: Publisher Site | Google Scholar
Y. Zhang, M. Manjunath, S. Zhang, D. Chasman, S. Roy, and J. Song, “Integrative genomic analysis predicts causative -regulatory mechanisms of the breast cancer-associated genetic variant rs4415084,” Cancer Research, vol. 78, no. 7, pp. 1579–1591, 2018.
View at: Publisher Site | Google Scholar
Y. Li, J. Zhang, C. Huo et al., “Dynamic organization of lncRNA and circular RNA regulators collectively controlled cardiac differentiation in humans,” EBioMedicine, vol. 24, pp. 137–146, 2017.
View at: Publisher Site | Google Scholar
M. Herriges, D. Tischfield, Z. Cui et al., “The NANCI-Nkx2.1 gene duplex buffers Nkx2.1 expression to maintain lung development and homeostasis,” Genes & Development, vol. 31, no. 9, pp. 889–903, 2017.
View at: Publisher Site | Google Scholar
L. Tian, J. Cao, X. Deng et al., “Unveiling transcription factor regulation and differential co-expression genes in duchenne muscular dystrophy,” Diagnostic Pathology, vol. 9, no. 1, Article ID 210, 2014.
View at: Publisher Site | Google Scholar
Z. H. Luo, Y. Li, X. F. Liu et al., “Systems biology of myasthenia gravis, integration of aberrant lncRNA and mRNA expression changes,” BMC Medical Genomics, vol. 8, no. 1, 13 pages, 2015.
View at: Publisher Site | Google Scholar
P.-E. Rautou, K. Tatsumi, S. Antoniak et al., “Hepatocyte tissue factor contributes to the hypercoagulable state in a mouse model of chronic liver injury,” Journal of Hepatology, vol. 64, no. 1, pp. 53–59, 2016.
View at: Publisher Site | Google Scholar
R. T. Stravitz, R. Bowling, R. L. Bradford et al., “Role of procoagulant microparticles in mediating complications and outcome of acute liver injury/acute liver failure,” Hepatology, vol. 58, no. 1, pp. 304–313, 2013.
View at: Publisher Site | Google Scholar
L. Salmena, L. Poliseno, Y. Tay, L. Kats, and P. Pandolfi, “A ceRNA hypothesis: the rosetta stone of a hidden RNA language?” Cell, vol. 146, no. 3, pp. 353–358, 2011.
View at: Publisher Site | Google Scholar
X. Han, Y. Hong, and K. Zhang, “TUG1 is involved in liver fibrosis and activation of HSCs by regulating miR-29b,” Biochemistry and Biophysical Research Communications, vol. 503, no. 3, pp. 1394–1400, 2018.
View at: Publisher Site | Google Scholar
X. Tu, Y. Zhang, X. Zheng et al., “TGF-β-induced hepatocyte lincRNA-p21 contributes to liver fibrosis in mice,” Scientific Reports, vol. 7, no. 1, 2957 pages, 2017.
View at: Publisher Site | Google Scholar
M. Gao, C. Li, M. Xu, Y. Liu, M. Cong, and S. Liu, “LncRNA MT1DP aggravates cadmium-induced oxidative stress by repressing the function of Nrf2 and is dependent on interaction with miR-365,” Advanced Science, vol. 5, no. 7, p. 1800087, 2018.
View at: Publisher Site | Google Scholar
L. Shen, Q. Wang, R. Liu et al., “LncRNA lnc-RI regulates homologous recombination repair of DNA double-strand breaks by stabilizing RAD51 mRNA as a competitive endogenous RNA,” Nucleic Acids Research, vol. 46, no. 2, pp. 717–729, 2018.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Sanyang Chen 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.

PDF Download Citation

Download other formats

Order printed copies

Views

989

Downloads

1358

Citations