Original Research ARTICLE
Integrated expression profiles of mRNA and microRNA in the liver of Fructus Meliae Toosendan water extract injured mice
- 1Pharmaceutical Informatics Institute, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, China
- 2Department of Pharmacology and Toxicology, Beijing Institute of Radiation Medicine, Beijing, China
Liver toxicity is a severe problem associated with Traditional Chinese Medicine (TCM). Fructus Meliae Toosendan (FMT) is a known hepatotoxic TCM, however, the toxicological mechanisms of liver injury caused by FMT treatment still remain largely unknown. In this study, we aimed to reveal possible mechanisms of FMT water extract-induced liver injury using a systemic approach. After three consecutive daily dosing of FMT water extract, significant increases of alanine transaminase, aspartate transaminase, and alkaline phosphatase activities, along with elevated total bilirubin and total cholesterol levels and a decrease of triglyceride level, were detected in mice serum. Moreover, hydropic degeneration was observed in hepatocytes, suggesting the presence of FMT-induced liver injury. mRNA and microRNA expression profiles of liver samples from injured mice were analyzed and revealed 8 miRNAs and 1,723 mRNAs were significantly changed after FMT water extract treatment. For the eight differentially expressed miRNAs, their predicted target genes were collected and a final set of 125 genes and 4 miRNAs (miR-139-5p, miR-199a-5p, miR-2861, and miR-3960) was selected to investigate important processes involved in FMT hepatotoxicity. Our results demonstrated several cellular functions were disordered after FMT treatment, such as cellular growth and proliferation, gene expression and cellular development. We hypothesized that liver cell necrosis was the main liver toxicity of FMT water extract, which was possibly caused by oxidative stress responses.
Drug-induced liver injury (DILI) is the main cause of drugs being removed from the market and it is a universal side effect of drugs (Halegoua-De and Navarro, 2008). Over 80% FDA-approved prescription drugs with the black box warning are related to severe DILI (Chen et al., 2011). Since liver plays a very important role in the detoxification of xenobiotics and it is prone to generate chemically reactive metabolites, it makes the liver being the target organ of DILI (Hornby et al., 2014). The most common histological patterns of DILI in clinic are the acute (21%) and chronic hepatitis (14%), acute (9%), and chronic cholestasis (10%), and cholestatic hepatitis (29%) (Kleiner et al., 2014). Among them, biochemical presentations of hepatocellular injury, cholestatic injury, and both are all detected in these histological patterns (Kleiner et al., 2014). It is worth to mention that the reports of Traditional Chinese Medicine (TCM)-induced liver injury are on the rise from many countries worldwide as it has been used more frequently. Liver injury becomes one of the main toxicities caused by TCM. Due to the characteristics of TCM, e.g., multi-components and multi-targets, it is difficult to understand the mechanism of TCM-induced liver injury.
Microarray provides more sensitive and detailed insights into the mechanisms of toxic action at an earlier molecular level and has been employed as a useful method to study the mechanisms of DILI (Lizarraga et al., 2012). Results from mRNA microarray analysis can also correlate to the observations from histopathology and clinical chemistry (Waring et al., 2001), which could offer further insight into the mechanisms of DILI.
In recent years, microRNAs (miRNAs) have drawn considerable attentions because of their critical effects on cellular processes, pathological processes, and physiological processes (Kloosterman and Plasterk, 2006). miRNAs are short noncoding single-stranded RNAs of approximately 22 nucleotides that decrease the expressions of target mRNAs (Giraldez et al., 2006; Guo et al., 2010) or influence protein production (Wightman et al., 1993; Olsen and Ambros, 1999) by binding to target mRNAs. miRNAs have been demonstrated to play an important role in DILI. For example, miR-122, the most abundant miRNA in liver, interacts with many targets involved in metabolism and stress-response pathways (Bala et al., 2012) and could be a potential biomarker of DILI (Starkey et al., 2011). miRNA microarray is a powerful tool to detect global alteration of miRNA levels and their functions during liver injury. Wang et al (Wang et al., 2009) showed that some specific miRNA species in blood, such as miR-122 and miR-192, have detectable dose- and exposure duration-dependent changes in APAP-induced liver injury in mice using miRNA microarray. In addition, these miRNAs were found to change at an even earlier stage than the traditional chemical markers used in clinic, such as alanine transaminase (ALT).
Integration of mRNA microarray and miRNA microarray would provide a more comprehensive understanding on mechanisms of DILI. This approach has been examined by Lizarraga et al. (2012) to evaluate benzo[α]pyrene (BaP)-induced toxicity at genetic level and eight differentially expressed miRNAs seemed to be involved in the pathways of apoptotic signaling, cell cycle arrest, DNA damage response, and DNA damage repair, thus leading to liver injury.
In this study, Fructus Meliae Toosendan (FMT) was used as a model drug, which is a typically hepatotoxic TCM. FMT has been used generally as an insecticidal and medicinal plant for a long time in China and Korea to treat stomach ache, cholelithiasis, cholecystitis, gastritis, mastitis, and ascariasis (Xie et al., 2008). Liver injury caused by FMT has been reported frequently in recent years (Yuen et al., 2006). Most attention has been concentrated on the single component of FMT, namely toosendanin, which has been identified as one of the toxic compounds based on the increase in serum glutamic pyruvic transaminase (SGPT) and changes in histopathology. Several toxicity mechanisms of toosendanin have been proposed, such as mitochondrial dysfunction, caspase induction, and ROS and MAP kinases (Zhang et al., 2008). However, it remains a question if toosendanin is solely responsible for FMT-induced liver injury since the clinical application form of FMT comprises of complex components. Moreover, TCM prescriptions are usually water-based solutions as well as FMT. Thus, it is interesting to investigate toxicological mechanisms of FMT water extract-induced liver injury. Here we applied an integrated approach of miRNA and mRNA expression profiles to systemically elucidate the potential toxicological mechanisms of FMT water extract-induced liver injury. Figure 1 summarized the overall research strategy and experimental design for this study.
Materials and Methods
Extraction of FMT
Dried whole FMT (6 kg, purchased from Zhejiang Chinese Medical University Medicine Plant, Hangzhou, China) was infused overnight with eightfold water (w:v), decocted for two times (90 min per time), and then filtered. The filtrates were concentrated to 6–7 L by rotary vacuum evaporation and precipitated by adding 95% ethanol to a final concentration of 75%. After laying the mixture overnight, the supernatant was concentrated to 500 mL with a rotary evaporator. The final extract was stored in -20°C.
Animals Treatment and Samples Collection
Male BALB/c mice (18–22 g, Silaike Co. Shanghai, China) were housed in controlled environmental conditions (12:12 h light/dark cycle and 25 ± 1°C) with conventional access to food and water. The mice were acclimatized for a period of 3 days prior to the start of the experiment. With the method of random grouping, mice were divided into two groups: the vehicle control group (n = 8, 1% sodium carboxymethyl cellulose) and the FMT water extract group (n = 16, 240 g/kg, an equivalent amount of the crude drug). Treatment was given at a frequency of intragastrical administration once a day for 3 days. In addition, individual body weight was recorded every day before oral gavage. Orbital blood samples were obtained after 6 h of the last oral gavage on the third day and centrifuged at 4°C, 4000 rpm for 15 min. The resulted serum was used for clinical biochemical analysis. 200 mg excised samples from the left lateral liver lobe were directly snap frozen in liquid nitrogen for microarray analysis (n = 5, chose the first five mice with higher ALT levels in FMT water extract treatment group, and randomly selected five mice in control group). The remaining liver lobes were used for histopathological examinations. All the protocols and studies involving the animals were conducted according to the Guiding Principles in the Use of Animals in Toxicology and the Animal Care and Use Committee of Zhejiang University School of Medicine.
Biochemical Assay and Histopathological Examinations
The several parameters in serum were detected using Cobas C8000 system (Roche Diagnostics, Germany) according to the manufacturer’s instructions, including ALT, aspartate transaminase (AST), alkaline phosphatase (ALP), total bilirubin (TBILI), total cholesterol (TCHOL), and triglyceride (TG).
All livers were fixed with 10% formalin, embedded in paraffin, sectioned into 4-μm thick slices, mounted on poly(L-lysine)-treated slides, stained with hematoxylin and eosin, and examined by Olympus optical microscope.
Total RNA, including miRNA, was extracted from frozen liver tissue using mirVanaTM miRNA Isolation Kit (Ambion, USA) following the manufacturer’s instructions and analyzed on an Agilent Bioanalyzer 2100 (Agilent technologies, USA). Only RNA Integrity Number (RIN) was greater than or equal to 7.0 and 28S/18S was greater than or equal to 0.7, the RNA sample was used for microarray analysis.
miRNA Expression Analysis
The Agilent Mouse miRNA (8∗60 K) V19.0 array containing 1,247 miRNAs was used for miRNA analysis. After isolation, miRNA molecular in total RNA was labeled using miRNA Complete Labeling and Hyb Kit (Agilent technologies, USA) and hybridized in hybridization Oven (Agilent technologies, USA), following the manufacturer’s instructions. Raw data were transformed by log2 and then normalized by Quantile algorithm within Gene Spring GX 12.6.1 (Agilent technologies, USA). Data analysis was carried out using Welch t-test. A filter on low gene expression was used to keep only the probes expressed greater than 4 of the log2 transformed mean intensity in at least one sample. The miRNAs with p-value less than 0.05 and an absolute fold change more than 1.5-fold (Chang et al., 2011) between FMT water extract and vehicle control groups were considered to be significantly changed.
mRNA Expression Analysis
mRNA expression profiles were obtained using the Affymetrix mouse 430 2.0 array as described in Lu et al. (2013). Raw data were imported into the ArrayTrack v 3.1.5 for further analysis. Microarray data were normalized by MAS 5.0 algorithm and were further normalized per chip to the same median intensity value of 1000. Welch t-test within ArrayTrack® was used to identify the differentially expressed genes (DEGs) between FMT water extract and vehicle control groups with cutoffs of the p-value < 0.05, absolute fold change > 2, and the mean channel intensity more than 250.
Real-time Quantitative PCR
To validate the results from microarray analysis, real-time quantitative PCR was carried out with the same RNA samples used for microarray analysis.
After reverse transcribed with miScript II RT Kit (Qiagen, Germany), three miRNAs were selected to validate the results from miRNAs microarray using real-time quantitative PCR with miScript SYBR Green PCR Kit (Qiagen, Germany) and individual specific primers (HuaDa, China) for miR-21a-3p (5′-CAACAGCAGTCGATGGGCTGTC-3′), miR-199a-5p (5′-CCCAGTGTTCAGACTACCTGTTC-3′), miR-139-5p (5′-TCTACAGTGCACGTGTCTCCAGT-3′), and U6 snRNA (5′- CTCGCTTCGGCAGCACA-3′). Furthermore, 11 mRNAs were chosen to validate the results from mRNAs microarray. After mRNAs were reverse transcribed with Oligo (dT)15 Primer (Promega, USA) and SuperScript II reverse transcriptase (Invitrogen, USA), real-time quantitative PCR was implemented for these testing mRNAs with QuantiFast SYBR Green PCR Kit (Qiagen, Germany) and individual specific primers from HuaDa (Table 1).
All real-time quantitative PCR was performed on an Eppendorf Mastercycler ep realplex4 according to the manu facturer’s instructions. 2-ΔΔCt method was used to generate the fold changes based on the previous study (Yang et al., 2011). Housekeeping gene U6 was a normalizer for these three miRNAs while GADPH was used to normalize the data for mRNAs.
Ingenuity Pathways Analysis
The functions of differentially expressed miRNAs were further analyzed using Ingenuity Pathway Analysis software (IPA, Ingenuity Systems, USA). First, the target genes of differentially expressed miRNAs were predicted by IPA. Next, the intersection of the predicted target genes of differentially expressed miRNAs and the DEGs (p-value < 0.05, absolute fold change > 2) of FMT water extract treatment group was identified by IPA. The target mRNAs changed in the opposite direction of the corresponding miRNAs were selected for analysis by IPA Systems. Lastly, the biological and molecular functions of these genes were explored in IPA, such as relevant diseases and disorders, top canonical pathways, and top toxic lists. The canonical pathways of the DEGs were also performed by IPA as well as the network of the specific pathway.
Statistical analysis was performed by unpaired and two-tailed Student’s t-test. The level of significance was set at p-value < 0.05. Values are represented as mean ± standard deviation (SD) of triplicate determinations.
Effects of FMT Water Extract on Mouse Body Weight and Liver Histopathology
Compared with control group, there were no significant changes observed in body weight in FMT water extract group at any time point (Figure 2A). Two mice were dead in FMT water extract group at the third day after exposure to FMT water extract. Moreover, hydropic degeneration of hepatocytes was detected in the livers of FMT water extract-treated mice, whereas the livers of vehicle control group exhibited normal histology in histopathological examinations (Figure 2B).
FIGURE 2. (A) The body weights of mice treated with Fructus Meliae Toosendan (FMT) water extract. The results showed that there were no significant changes in body weight between vehicle control (1% sodium carboxymethyl cellulose, n = 8) and FMT water extract groups (240 g/kg, n = 16) by oral gavage for 3 days; (B) Histological analysis of the liver tissues of FMT water extract-treated mice. Hydropic degeneration of hepatocytes was detected in FMT water extract-treated mice. The scale bar is 50 μm.
Effects of FMT Water Extract on Serum Biochemistry
The biochemical parameters, including ALT, AST, ALP, TG, TBIL, and TCHOL, were used as the biochemical markers to monitor FMT water extract-induced liver injury. As shown in Figure 3, intragastric administration of FMT water extract (240 g/kg) caused significant increases in serum TBIL and TCHOL levels as well as ALT, AST, and ALP activities compared with control group, whereas remarkable decrease of TG level in FMT water extract group was detected.
FIGURE 3. Effects of FMT water extract on serum biochemical parameters after intragastric administration of 240 g/kg FMT water extract for 3 days in mice (control group: n = 8; FMT water extract group: n = 14). Values are represented as mean ± SD (∗p < 0.05, ∗∗p < 0.01 compared with vehicle control group).
Effects of FMT Water Extract on miRNA Expression Profiles
miRNA expression profiles were analyzed with Agilent Mouse miRNA (8∗60 K) V19.0 array. Eight miRNAs (approximately 0.64%) of total 1,247 miRNAs investigated were differentially expressed in FMT water extract-treated group (Figure 4A). Among these miRNAs, six miRNAs were upregulated over 1.5-fold, whereas two miRNAs were downregulated over 1.5-fold. Hierarchical cluster analysis (HCA) within ArrayTrack® was used to visualize clusters of the liver samples corresponding to the differentially expressed miRNAs (Figure 4B). The results showed that all samples were grouped into two main branches (FMT water extract-treated mice versus control mice). Three of these eight miRNAs were verified by real-time quantitative PCR. It was demonstrated that the expression of miR-21a-3p was upregulated and miR-199a-5p and miR-139-5p were downregulated in FMT water extract-treated group when compared with control group (Figure 5), which was consistent with the results from microarray analysis.
FIGURE 4. (A) The differentially expressed miRNAs in FMT water extract group. There were eight miRNAs (0.64%) of total 1,247 studied miRNAs differentially expressed in FMT water extract group when compared with control group (n = 5). six miRNAs were upregulated (fold change > 1.5) and 2 miRNAs were downregulated (fold change < -1.5); (B) Hierarchical Cluster Analysis (HCA) within ArrayTrack® was used to visualize clusters of samples corresponding to the differentially expressed miRNAs. There were two main branches detected (FMT water extract-treated mice versus control mice).
FIGURE 5. Validation of microarray results by real-time quantitative PCR. Gray bars (RT) indicate the mean fold change (±SD) derived from triplicate real-time quantitative PCR reactions (n = 5). Black bars (ARRAY) indicate microarray data.
Effects of FMT Water Extract on mRNA Expression Profiles
We observed that approximately 3.8% (1723 mRNAs) of total genes (45,101 mRNAs) were significantly changed in the liver of FMT water extract-treated mice (Figure 6A). Among them, 1,111 mRNAs were upregulated more than twofold (p < 0.05) and 612 mRNAs were downregulated more than twofold (p < 0.05) in FMT water extract group compared with control group. Moreover, HCA results displayed that all samples were grouped into two main branches (FMT water extract-treated mice versus control mice, Figure 6B). To validate the results obtained from mRNA microarray, the expression levels of 11 DEGs were detected by real-time quantitative PCR. As a result, the general trends observed in the microarray analysis of all genes tested (Gadd45a, Cidec, Cyp4a14, Cd14, Tnfrsf12a, Cxcl2, S100a9, Jun, Map3k5, Prkca, and Car3) were verified by real-time quantitative PCR (Figure 5).
FIGURE 6. (A) 1723 mRNAs (3.8%) of total 45,101 genes were significantly changed by FMT water extract administration detected by Affymetrix mouse 430 2.0 array. Compared with control group, 1,111mRNAs were upregulated (fold change > 2) as well as 612 mRNAs downregulated (fold change < -2); (B) HCA within ArrayTrack was used to visualize clusters of samples corresponding to the DEGs. There were two main branches detected (FMT water extract-treated mice versus control mice).
Prediction of Biological Function and Disease Establishment Associated with Intersection of miRNAs Targets and DEGs
Based on Affymetrix probe ID, 1594 of 1723 DEGs in FMT water extract treatment could be mapped in IPA. Moreover, the target genes of four miRNAs could be predicted by target filter analysis in IPA database, including miR-139-5p, miR-199a-5p, miR-2861, and miR-3960. As miRNAs usually decrease the expressions of target mRNAs, the intersection of the 1594 DEGs and the predicted target genes of these four miRNAs was identified, which was based on the opposite expression direction between miRNAs and mRNAs. As the results, 125 mRNAs showed an opposite expression trend compared with corresponding miRNA. Furthermore, the intersected 125 mRNAs were analyzed to predict biological functions using IPA bioinformatics tool. We found that most of the genes were involved in molecular and cellular functions, such as cellular development, growth and proliferation, gene expression, cell morphology, and cellular assembly and organization (Table 2). In addition, the top 5 canonical pathways (according to the p-value for these 125 mRNAs were involved in tight junction signaling, retinoid acid receptor (RAR) activation, NF-E2-related factor 2 (NRF2)-mediated oxidative stress response, hepatocyte growth factor (HGF) signaling, and planar cell polarity (PCP) pathway (Table 3). In top 5 toxic lists (according to the p-value), liver necrosis or cell death was the most relevant toxicity for FMT water extract-induced liver injury (Table 4).
TABLE 2. Top five cellular functions predicted by Ingenuity Pathway Analysis (IPA) that were corresponding to 125 mRNAs of the intersection of the four differentially expressed miRNAs predicted target genes and differentially expressed genes (DEGs) under Fructus Meliae Toosendan (FMT) water extract treatment.
TABLE 3. Top five canonical pathways predicted by IPA that were corresponding to 125 mRNAs of the intersection of the four differentially expressed miRNAs predicted target genes and DEGs under FMT water extract treatment.
TABLE 4. Top five toxic lists predicted by IPA that were corresponding to 125 mRNAs of the intersection of the four differentially expressed miRNAs predicted target genes and DEGs under FMT water extract treatment.
The Network of NRF2-mediated Oxidative Stress Response Associated with DEGs and Differentially Expressed miRNAs
The biological functions of the 1594 DEGs in FMT water extract treatment were analyzed by IPA. As shown in Table 5, the NRF2-mediated oxidative stress response was the most significant canonical pathway enriched in this dataset, which was also predicted by the 125 intersected mRNAs (Table 3). With the combination of the affected miRNAs and mRNAs, a network of NRF2-mediated oxidative stress response was generated by IPA and shown in Figure 7. As the results, lots of signaling molecules in this pathway were interfered by FMT water extract treatment.
TABLE 5. Top five canonical pathways predicted by IPA that were corresponding to the 1594 DEGs of FMT water extract treatment.
FIGURE 7. A network of NRF2-mediated oxidative stress response associated with DEGs and differentially expressed miRNAs. The yellow color means the up-regulation of the mRNAs and miRNAs after exposure to FMT water extract compared with control group, while the green color indicates the down-regulation of these molecules after exposure to FMT water extract. Different shapes outside of the molecules are illustrated at the figure.
Traditional Chinese Medicine has been widely employed around the world in many diseases, however, liver injury is one of the main side effects after exposure to TCM and it has been increased in recent years (Wai, 2006). As the components of TCM are very plentiful, it is complicated to explore the mechanism of TCM-induced liver injury. Shedding some light on this issue, miRNA microarray has already provided insights into the mechanisms and biomarkers of chemical DILI (Fukushima et al., 2007), which could be a new approach to study TCM-induced liver injury. Thus, in this study, we first applied an integrated approach using mRNA and miRNA microarray to study the molecular mechanisms of TCM-induced liver injury and chose FMT as a model drug.
Serum aminotransferases are the standard biomarkers for liver injury. When liver injured, the activities of ALP, ALT, and AST in serum could increase. Studies reported that the serum levels of ALT and AST were correlated with the degree of the liver injury (Myers et al., 2003; Lee et al., 2004). In this study, traditional serum transaminases (AST, ALT, and ALP) from FMT water extract group were higher than the control group (p < 0.01), indicating liver injury induced by FMT water extract treatment in mice. It may be caused by cell necrosis as the increases of cell membrane permeability which allowed the transaminases leaking from the cells. Consistent with that, hydropic degeneration of hepatocytes was detected in FMT water extract-treated mice. When hepatocytes were swelling, it might decrease albumin production and result in increasing bilirubin in the blood at the later stage (Lu et al., 2011). Concordantly, an increased level of serum TBIL was observed by FMT water extract treatment. The decreased level of TG was also detected after FMT water extract treatment, indicating the impact of FMT water extract on lipid metabolism. Thus, changes in serum biochemistry and histopathology indicated that FMT water extract could cause liver injury.
There were eight miRNAs differentially expressed in livers of FMT water extract-treated mice (Figure 4A). Interestingly, miR-122 and miR-192 did not show any differences between these two groups in our experiment though they are the biomarkers for liver injury (Su et al., 2012). Similarly, Bala et al. (2012) reported that hepatic miR-122 did not change while circulating miR-122 increased in DILI. The possible explanation was that the distinct experiment samples could lead to the contrary results. In our experiment, we chose livers as the source of miRNAs while these two biomarkers (miR-122 and miR-192) were proved to be useful in circulating miRNA (Starkey et al., 2011). The circulating miRNAs are transported by exosomes, microparticles, lipoproteins, and other ribonucleoprotein complexes (Kosaka et al., 2010; Arroyo et al., 2011; Boon and Vickers, 2013). It demonstrated that the miRNA profiles of extracellular vesicles and lipoproteins are not representative of their parent cell type but are distinct sets of miRNAs (Valadi et al., 2007; Wang et al., 2010; Zhang et al., 2010), and probably reach many distant tissues to modulate systemic homeostasis or promote the disease (Boon and Vickers, 2013). In this study, the cardiac hypertrophy and acute renal failure panel were also predicted in the top 5 toxic lists by IPA, indicating these distant organs might be also affected by FMT water extract or the secondary responses of extracellular miRNAs in liver toxicity. It partly corresponded with the reports that FMT can induce hepatotoxicity, cardiotoxicity, and nephrotoxicity (Lihua et al., 2006; Shuangyan et al., 2008).
Although only hydropic degeneration of hepatocytes rather than hepatocytes necrosis was detected in the livers of FMT water extract-treated mice by histopathological examinations, the toxic lists in IPA showed that the main toxicity of FMT water extract was liver necrosis or cell death. The possible reason for the phenomenon is that the alterations in gene expression are always earlier than the changes in protein and histopathological levels (Beer et al., 2002). In this study, we focused on liver necrosis or cell death as it was the most significant pathway enriched in toxic lists and it corresponded to the other canonical pathways predicted by IPA after exposure to FMT water extract (Tables 3 and 4), such as HGF signaling, RAR activation, NRF2-mediated oxidative stress response, tight junction signaling, and PCP pathway. More specifically, HGF was a multifunctional cytokine that stimulated cell necrosis (Suzuki et al., 2000). Wang et al. (2007) reported that activation of RAR alpha (retinoic acid receptors) could inhibit Nrf2, which regulated detoxification processes (Chan et al., 2001). NRF2-mediated oxidative stress related to liver diseases (Cederbaum et al., 2009) and participated in cell death, proliferation and differentiation (Richter and Kass, 1991). In addition, tight junction, which was detected in flutamide-induced liver toxicity (Legendre et al., 2014), played an important role in cell apoptosis (Zhang et al., 2011). PCP pathway was a non-canonical Wnt signaling pathway, which directed processes such as cell fate determination, proliferation, and differentiation (Teoa and Kahn, 2010). Thus, it indicated that these pathways related to cell necrosis or death and may play important roles in FMT water extract-induced liver injury.
Corresponding to the canonical pathways revealed by the 125 intersected mRNAs, the NRF2-mediated oxidative stress response was also enriched as the most significant canonical pathway by the 1594 DEGs of FMT water extract group. Thus, we speculated that the oxidative stress might be the primary cause of FMT-induced liver damage. As shown in Figure 7, lots of signaling molecules in the NRF2-mediated oxidative stress response were interfered by FMT water extract treatment. We noticed that ASK1, PKC, and JUN affected by the differentially expressed miRNAs of FMT water extract, i.e., miR-139-5p and miR-199-5p. The formers were the upstream molecules of NRF2 and played important roles in oxidative stress response. In particular, ASK1, known as Map3k5, was the target of miR-139-5p. It mediated signal transduction of various stressors such as oxidative stress (Noguchi et al., 2005), and played a crucial role in cell survival (Zhang et al., 1999) and hepatic proliferation (Noguchi et al., 2014). PRKCA, a member of PRC, was the target of miR-139-5p. It was activated in oxidative stress-related diseases, such as cancer, cerebral ischemia-reperfusion injury, hepatic damage, and involved in positive and negative regulation of cell proliferation, apoptosis, differentiation, and cell injury (Mandil et al., 2001; Nitti et al., 2008; Kamo et al., 2011). JUN was the common molecule of NRF2-mediated oxidative stress response and cell necrosis from toxic lists of the 125 intersected mRNAs. It was regulated by miR-139-5p and miR-199a-5p and activated by the oxidative stress (Lee et al., 2001). Previous studies reported that JUN worked as a regulator of cell fate, such as hepatocyte survival, proliferation, and liver tumorigenesis (Kayahara et al., 2005; Fuest et al., 2012). Taken together, oxidative stress and subsequent liver cell necrosis may be the toxicological mechanisms of FMT-induced liver injury.
A literature-based analysis was performed to increase evidence of the IPA results for the four miRNAs with targets in our experimental data. These miRNAs were involved in some processes such as cellular growth and proliferation, gene expression, cancer, and cellular development. For example, miR-3960 regulated cellular growth and proliferation through a regulatory feedback loop with miR-2861 (Hu et al., 2011; Wong et al., 2011), meanwhile, these two miRNAs could relate with oxidative stress response (Dalle et al., 2012). The deregulation of miR-2861 and miR-21a-3p was also related to cancer diseases (Wang et al., 2013; Gordillo et al., 2014). In addition, miR-139-5p was a tumor-suppressor miRNA and had been found downregulated in hepatoma cells samples when compared with adjacent benign tissues (Jiang et al., 2008; Au et al., 2012; Huang et al., 2012). Similarly, miR-199a-5p, one of the most abundant miRNAs in hepatocytes, prevented hepatocyte apoptosis (Dai et al., 2013). Thus, the downregulation of miR-139-5p and miR-199a-5p in this study may indicate the potential risk of hepatic carcinoma caused by FMT water extract. Moreover, the dysregulation of miR-21a and miR-199a-5p have been also reported by other liver injury models. Specifically, the increase of miR-21a expression was detected in alcoholic liver injury (Francis et al., 2014), and the expression of miR-199a-5p was elevated in both bile acid- and thapsigargin-stimulated cultured hepatocytes, as well as in the liver of bile duct-ligated mice (Dai et al., 2013). In addition, since the microarray results might not come exclusively from hepatocytes as the liver tissue used for this analysis, we further discussed the possible contribution of the other cell types in liver to the pathology and the expression profile of miRNAs. Among these four miRNAs, the miR-199a-5p and miRNA-21 were found to be expressed by other cell types of liver and contributed to liver diseases. For instance, hepatic stellate cells had been found to express miR-199, which influenced stellate cell activity by regulating collagen synthesis and represented a general mechanism contributing to hepatic fibrosis (Lino et al., 2013); the liver sinusoidal endothelial cells also express the miR-199 and this miRNA might serve as a negative regulator of the function of liver sinusoidal endothelial cells in some chemical exposures, such as alcohol (Wang et al., 2012); the increased miR-21 expression was detected in the hepatic stellate cells with the treatment of ethanol and interleukin-6 (Francis et al., 2014).
Fructus Meliae Toosendan water extract can induce the liver injury in a mice model, and integrated microRNA-mRNA approach provided insights into the molecular mechanisms of FMT water extract-induced liver injury, which demonstrated the feasibility of applying a microRNA-mRNA network based approach to study TCM-induced liver injury.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was financially supported by the Zhejiang Provincial Natural Science Foundation of China (No.LY13B070002), the Specialized Research Fund for the Doctoral Program of Higher Education (No. 20120101120137), and the Fundamental Research Funds for the Central Universities (No. 2015QNA7038).
Arroyo, J. D., Chevillet, J. R., Kroh, E. M., Ruf, I. K., Pritchard, C. C., Gibson, D. F., et al. (2011). Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc. Natl. Acad. Sci. U.S.A. 108, 5003–5008. doi: 10.1073/pnas.1019055108
Au, S. L., Wong, C. C., Lee, J. M., Fan, D. N., Tsang, F. H., Ng, I. O., et al. (2012). Enhancer of zeste homolog 2 epigenetically silences multiple tumor suppressor microRNAs to promote liver cancer metastasis. Hepatology 56, 622–631. doi: 10.1002/hep.25679
Bala, S., Petrasek, J., Mundkur, S., Catalano, D., Levin, I., Ward, J., et al. (2012). Circulating microRNAs in exosomes indicate hepatocyte injury and inflammation in alcoholic, drug-induced, and inflammatory liver diseases. Hepatology 56, 1946–1957. doi: 10.1002/hep.25873
Beer, D. G., Kardia, S., Huang, C. C., Giordano, T. J., Levin, A. M., Misek, D. E., et al. (2002). Gene-expression profiles predict survival of patients with lung adenocarcinoma. Nat. Med. 8, 816–824.
Chan, K., Han, X. D., and Kan, Y. W. (2001). An important function of Nrf2 in combating oxidative stress: detoxification of acetaminophen. Proc. Natl. Acad. Sci. U.S.A. 98, 4611–4616. doi: 10.1073/pnas.081082098
Chang, S. J., Weng, S. L., Hsieh, J. Y., Wang, T. Y., Chang, M. D., and Wang, H. W. (2011). MicroRNA-34a modulates genes involved in cellular motility and oxidative phosphorylation in neural precursors derived from human umbilical cord mesenchymal stem cells. BMC Med. Genomics 4:65. doi: 10.1186/1755-8794-4-65
Chen, M., Vijay, V., Shi, Q., Liu, Z., Fang, H., and Tong, W. (2011). FDA-approved drug labeling for the study of drug-induced liver injury. Drug Discov. Today 16, 697–703. doi: 10.1016/j.drudis.2011.05.007
Dai, B. H., Geng, L., Wang, Y., Sui, C. J., Xie, F., Shen, R. X., et al. (2013). microRNA-199a-5p protects hepatocytes from bile acid-induced sustained endoplasmic reticulum stress. Cell Death Dis. 4:e604. doi: 10.1038/cddis.2013.134
Francis, H., McDaniel, K., Han, Y., Liu, X., Kennedy, L., Yang, F., et al. (2014). Regulation of the extrinsic apoptotic pathway by microRNA-21 in alcoholic liver injury. J. Biol. Chem. 289, 27526–27539. doi: 10.1074/jbc.M114.602383
Fuest, M., Willim, K., MacNelly, S., Fellner, N., Resch, G. P., Blum, H. E., et al. (2012). The transcription factor c-Jun protects against sustained hepatic endoplasmic reticulum stress thereby promoting hepatocyte survival. Hepatology 55, 408–418. doi: 10.1002/hep.24699
Fukushima, T., Hamada, Y., Yamada, H., and Horii, I. (2007). Changes of micro-RNA expression in rat liver treated by acetaminophen or carbon tetrachloride–regulating role of micro-RNA for RNA expression. J. Toxicol. Sci. 32, 401–409. doi: 10.2131/jts.32.401
Giraldez, A. J., Mishima, Y., Rihel, J., Grocock, R. J., Van Dongen, S., Inoue, K., et al. (2006). Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science 312, 75–79. doi: 10.1126/science.1122689
Gordillo, G. M., Biswas, A., Khanna, S., Pan, X., Sinha, M., Roy, S., et al. (2014). Dicer knockdown inhibits endothelial cell tumor growth via microRNA 21a-3p targeting of Nox-4. J. Biol. Chem. 289, 9027–9038. doi: 10.1074/jbc.M113.519264
Hornby, R. J., Lewis, P. S., Dear, J., Goldring, C., and Park, B. K. (2014). MicroRNAs as potential circulating biomarkers of drug-induced liver injury: key current and future issues for translation to humans. Exp. Rev. Clin. Pharmacol. 7, 349–362. doi: 10.1586/17512433.2014.904201
Hu, R., Liu, W., Li, H., Yang, L., Chen, C., Xia, Z. Y., et al. (2011). A Runx2/miR-3960/miR-2861 regulatory feedback loop during mouse osteoblast differentiation. J. Biol. Chem. 286, 12328–12339. doi: 10.1074/jbc.M110.176099
Huang, Y., Chen, H. C., Chiang, C. W., Yeh, C. T., Chen, S. J., and Chou, C. K. (2012). Identification of a two-layer regulatory network of proliferation-related microRNAs in hepatoma cells. Nucleic Acids Res. 40, 10478–10493. doi: 10.1093/nar/gks789
Jiang, J., Gusev, Y., Aderca, I., Mettler, T. A., Nagorney, D. M., Brackett, D. J., et al. (2008). Association of MicroRNA expression in hepatocellular carcinomas with hepatitis infection, cirrhosis, and patient survival. Clin. Cancer Res. 14, 419–427. doi: 10.1158/1078-0432.CCR-07-0523
Kamo, N., Shen, X. D., Ke, B., Busuttil, R. W., and Kupiec-Weglinski, J. W. (2011). Sotrastaurin, a protein kinase C inhibitor, ameliorates ischemia and reperfusion injury in rat orthotopic liver transplantation. Am. J. Transplant. 11, 2499–2507. doi: 10.1111/j.1600-6143.2011.03700.x
Kayahara, M., Wang, X., and Tournier, C. (2005). Selective regulation of c-jun gene expression by mitogen-activated protein kinases via the 12-o-tetradecanoylphorbol-13-acetate- responsive element and myocyte enhancer factor 2 binding sites. Mol. Cell. Biol. 25, 3784–3792. doi: 10.1128/MCB.25.9.3784-3792.2005
Kleiner, D. E., Chalasani, N. P., Lee, W. M., Fontana, R. J., Bonkovsky, H. L., Watkins, P. B., et al. (2014). Hepatic histological findings in suspected drug-induced liver injury: systematic evaluation and clinical associations. Hepatology 59, 661–670. doi: 10.1002/hep.26709
Kosaka, N., Iguchi, H., and Ochiya, T. (2010). Circulating microRNA in body fluid: a new potential biomarker for cancer diagnosis and prognosis. Cancer Sci. 101, 2087–2092. doi: 10.1111/j.1349-7006.2010.01650.x
Lee, J. Y., Chae, D. W., Kim, S. M., Nam, E. S., Jang, M. K., Lee, J. H., et al. (2004). Expression of FasL and perforin/granzyme B mRNA in chronic hepatitis B virus infection. J. Viral. Hepat. 11, 130–135. doi: 10.1046/j.1365-2893.2003.00486.x
Lee, S. A., Dritschilo, A., and Jung, M. (2001). Role of ATM in oxidative stress-mediated c-Jun phosphorylation in response to ionizing radiation and CdCl2. J. Biol. Chem. 276, 11783–11790. doi: 10.1074/jbc.M004517200
Legendre, A., Jacques, S., Dumont, F., Cotton, J., Paullier, P., Fleury, M. J., et al. (2014). Investigation of the hepatotoxicity of flutamide: pro-survival/apoptotic and necrotic switch in primary rat hepatocytes characterized by metabolic and transcriptomic profiles in microfluidic liver biochips. Toxicol. In Vitro 28, 1075–1087. doi: 10.1016/j.tiv.2014.04.008
Lihua, X., Yanping, C., Zhi, Z., Yuzeng, S., Zhimin, Y., and Aiping, L. (2006). [Influence of 24 Traditional Chinese Herbs to Viscera Index of Rat in Long Term Toxicity Test]. Chin. J. Basic Med. Tradit. Chin. Med. 12, 35–36.
Lino, C. C., Henaoui, I. S., Courcot, E., Roderburg, C., Cauffiez, C., Aubert, S., et al. (2013). miR-199a-5p Is upregulated during fibrogenic response to tissue injury and mediates TGFbeta-induced lung fibroblast activation by targeting caveolin-1. PLoS Genet. 9:e1003291. doi: 10.1371/journal.pgen.1003291
Lizarraga, D., Gaj, S., Brauers, K. J., Timmermans, L., Kleinjans, J. C., and van Delft, J. H. (2012). Benzo[a]pyrene-induced changes in microRNA-mRNA networks. Chem. Res. Toxicol. 25, 838–849. doi: 10.1021/tx2003799
Lu, X., Jin, T., Jin, Y., Wu, L., Hu, B., Tian, Y., et al. (2013). Toxicogenomic analysis of the particle dose- and size-response relationship of silica particles-induced toxicity in mice. Nanotechnology 24:15106. doi: 10.1088/0957-4484/24/1/015106
Lu, X., Zhao, Q., Tian, Y., Xiao, S., Jin, T., and Fan, X. (2011). A metabonomic characterization of (+)-usnic acid-induced liver injury by gas chromatography-mass spectrometry-based metabolic profiling of the plasma and liver in rat. Int. J. Toxicol. 30, 478–491. doi: 10.1177/1091581811414436
Mandil, R., Ashkenazi, E., Blass, M., Kronfeld, I., Kazimirsky, G., Rosenthal, G., et al. (2001). Protein kinase Calpha and protein kinase Cdelta play opposite roles in the proliferation and apoptosis of glioma cells. Cancer Res. 61, 4612–4619.
Myers, R. P., Tainturier, M. H., Ratziu, V., Piton, A., Thibault, V., Imbert-Bismut, F., et al. (2003). Prediction of liver histological lesions with biochemical markers in patients with chronic hepatitis B. J. Hepatol. 39, 222–230. doi: 10.1016/S0168-8278(03)00171-5
Noguchi, H., Yamada, S., Nabeshima, A., Guo, X., Tanimoto, A., Wang, K. Y., et al. (2014). Depletion of apoptosis signal-regulating kinase 1 prevents bile duct ligation-induced necroinflammation and subsequent peribiliary fibrosis. Am. J. Pathol. 184, 644–661. doi: 10.1016/j.ajpath.2013.11.030
Noguchi, T., Takeda, K., Matsuzawa, A., Saegusa, K., Nakano, H., Gohda, J., et al. (2005). Recruitment of tumor necrosis factor receptor-associated factor family proteins to apoptosis signal-regulating kinase 1 signalosome is essential for oxidative stress-induced cell death. J. Biol. Chem. 280, 37033–37040. doi: 10.1074/jbc.M506771200
Olsen, P. H., and Ambros, V. (1999). The lin-4 regulatory RNA controls developmental timing in Caenorhabditis elegans by blocking LIN-14 protein synthesis after the initiation of translation. Dev. Biol. 216, 671–680. doi: 10.1006/dbio.1999.9523
Richter, C., and Kass, G. E. (1991). Oxidative stress in mitochondria: its relationship to cellular Ca2+ homeostasis, cell death, proliferation, and differentiation. Chem. Biol. Interact. 77, 1–23. doi: 10.1016/0009-2797(91)90002-O
Starkey, L. P., Dear, J., Platt, V., Simpson, K. J., Craig, D. G., Antoine, D. J., et al. (2011). Circulating microRNAs as potential markers of human drug-induced liver injury. Hepatology 54, 1767–1776. doi: 10.1002/hep.24538
Su, Y. W., Chen, X., Jiang, Z. Z., Wang, T., Wang, C., Zhang, Y., et al. (2012). A panel of serum microRNAs as specific biomarkers for diagnosis of compound- and herb-induced liver injury in rats. PLoS ONE 7:e37395. doi: 10.1371/journal.pone.0037395
Suzuki, A., Hayashida, M., Kawano, H., Sugimoto, K., Nakano, T., and Shiraki, K. (2000). Hepatocyte growth factor promotes cell survival from fas-mediated cell death in hepatocellular carcinoma cells via Akt activation and Fas-death-inducing signaling complex suppression. Hepatology 32, 796–802. doi: 10.1053/jhep.2000.17738
Teoa, J., and Kahn, M. (2010). The Wnt signaling pathway in cellular proliferation and differentiation: a tale of two coactivators. Adv. Drug Deliver. Rev. 62, 1149–1155. doi: 10.1016/j.addr.2010.09.012
Valadi, H., Ekstrom, K., Bossios, A., Sjostrand, M., Lee, J. J., and Lotvall, J. O. (2007). Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat. Cell Biol. 9, 654–672. doi: 10.1038/ncb1596
Wang, K., Zhang, S., Marzolf, B., Troisch, P., Brightman, A., Hu, Z., et al. (2009). Circulating microRNAs, potential biomarkers for drug-induced liver injury. Proc. Natl. Acad. Sci. U.S.A. 106, 4402–4407. doi: 10.1073/pnas.0813371106
Wang, X. J., Hayes, J. D., Henderson, C. J., and Wolf, C. R. (2007). Identification of retinoic acid as an inhibitor of transcription factor Nrf2 through activation of retinoic acid receptor alpha. Proc. Natl. Acad. Sci. U.S.A. 104, 19589–19594. doi: 10.1073/pnas.0709483104
Wang, Z., Zhang, H., Zhang, P., Li, J., Shan, Z., and Teng, W. (2013). Upregulation of miR-2861 and miR-451 expression in papillary thyroid carcinoma with lymph node metastasis. Med. Oncol. 30:577. doi: 10.1007/s12032-013-0577-9
Waring, J. F., Jolly, R. A., Ciurlionis, R., Lum, P. Y., Praestgaard, J. T., Morfitt, D. C., et al. (2001). Clustering of hepatotoxins based on mechanism of toxicity using gene expression profiles. Toxicol. Appl. Pharmacol. 175, 28–42. doi: 10.1006/taap.2001.9243
Wightman, B., Ha, I., and Ruvkun, G. (1993). Posttranscriptional regulation of the heterochronic gene lin-14 by lin-4 mediates temporal pattern formation in C. elegans. Cell 75, 855–862. doi: 10.1016/0092-8674(93)90530-4
Wong, C. C., Wong, C. M., Tung, E. K., Au, S. L., Lee, J. M., Poon, R. T., et al. (2011). The microRNA miR-139 suppresses metastasis and progression of hepatocellular carcinoma by down-regulating Rho-kinase 2. Gastroenterology 140, 322–331. doi: 10.1053/j.gastro.2010.10.006
Xie, F., Zhang, M., Zhang, C. F., Wang, Z. T., Yu, B. Y., and Kou, J. P. (2008). Anti-inflammatory and analgesic activities of ethanolic extract and two limonoids from Melia toosendan fruit. J. Ethnopharmacol. 117, 463–466. doi: 10.1016/j.jep.2008.02.025
Yang, B., Guo, H., Zhang, Y., Chen, L., Ying, D., and Dong, S. (2011). MicroRNA-145 regulates chondrogenic differentiation of mesenchymal stem cells by targeting Sox9. PLoS ONE 6:e21679. doi: 10.1371/journal.pone.0021679
Yuen, M. F., Tam, S., Fung, J., Wong, D. K., Wong, B. C., and Lai, C. L. (2006). Traditional Chinese medicine causing hepatotoxicity in patients with chronic hepatitis B infection: a 1-year prospective study. Aliment. Pharmacol. Ther. 24, 1179–1186. doi: 10.1111/j.1365-2036.2006.03111.x
Zhang, G., Ding, L., Renegar, R., Wang, X., Lu, Q., Huo, S., et al. (2011). Hydroxycamptothecin-loaded Fe3O4 nanoparticles induce human lung cancer cell apoptosis through caspase-8 pathway activation and disrupt tight junctions. Cancer Sci. 102, 1216–1222. doi: 10.1111/j.1349-7006.2011.01930.x
Zhang, L., Chen, J., and Fu, H. (1999). Suppression of apoptosis signal-regulating kinase 1-induced cell death by 14-3-3 proteins. Proc. Natl. Acad. Sci. U.S.A. 96, 8511–8515. doi: 10.1073/pnas.96.15.8511
Zhang, Y., Liu, D., Chen, X., Li, J., Li, L., Bian, Z., et al. (2010). Secreted monocytic miR-150 enhances targeted endothelial cell migration. Mol. Cell. 39, 133–144. doi: 10.1016/j.molcel.2010.06.010
Keywords: DILI, Fructus Meliae Toosendan (FMT), miRNA, mRNA
Citation: Zheng J, Ji C, Lu X, Tong W, Fan X and Gao Y (2015) Integrated expression profiles of mRNA and microRNA in the liver of Fructus Meliae Toosendan water extract injured mice. Front. Pharmacol. 6:236. doi: 10.3389/fphar.2015.00236
Received: 02 June 2015; Accepted: 02 October 2015;
Published: 19 October 2015.
Edited by:Adolfo Andrade-Cetto, Universidad Nacional Autónoma de México, Mexico
Reviewed by:Maria Luisa Del Moral, University of Jaen, Spain
Xiao Yu Tian, Houston Methodist Research Institute, USA
Copyright © 2015 Zheng, Ji, Lu, Tong, Fan and Gao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.