Growth Hormone-Regulated mRNAs and miRNAs in Chicken Hepatocytes

Growth hormone (GH) is a key regulatory factor in animal growth, development and metabolism. Based on the expression level of the GH receptor, the chicken liver is a major target organ of GH, but the biological effects of GH on the chicken liver are not fully understood. In this work we identified mRNAs and miRNAs that are regulated by GH in primary hepatocytes from female chickens through RNA-seq, and analyzed the functional relevance of these mRNAs and miRNAs through GO enrichment analysis and miRNA target prediction. A total of 164 mRNAs were found to be differentially expressed between GH-treated and control chicken hepatocytes, of which 112 were up-regulated and 52 were down-regulated by GH. A total of 225 chicken miRNAs were identified by the RNA-Seq analysis. Among these miRNAs 16 were up-regulated and 1 miRNA was down-regulated by GH. The GH-regulated mRNAs were mainly involved in growth and metabolism. Most of the GH-upregulated or GH-downregulated miRNAs were predicted to target the GH-downregulated or GH-upregulated mRNAs, respectively, involved in lipid metabolism. This study reveals that GH regulates the expression of many mRNAs involved in metabolism in female chicken hepatocytes, which suggests that GH plays an important role in regulating liver metabolism in female chickens. The results of this study also support the hypothesis that GH regulates lipid metabolism in chicken liver in part by regulating the expression of miRNAs that target the mRNAs involved in lipid metabolism.


Introduction
Growth hormone (GH) is a peptide hormone from the anterior pituitary gland [1,2]. It has many biological effects at both the whole body and tissue levels [3]. GH regulates animal growth, development and metabolism [3,4,5,6]. GH regulates the metabolism of not only protein but also that of lipid and carbohydrates [7]. GH initiates its function by binding to the GH receptor (GHR) [8,9,10]. Binding of GH to the GHR activates the receptor-associated tyrosine kinase JAK2 [11], and JAK2 then activates multiple proteins, including STAT1, STAT3, STAT5, MAPK, and PI3K [12]. These proteins in turn mediate GHcaused changes in gene expression or protein modification.
Liver is a key metabolic organ, and in chickens this is where most of the de novo synthesis of fatty acids occurs [13,14]. microRNAs (miRNAs) are a class of small non-coding RNAs about 22 nucleotides in length, and regulate gene expression by interacting with the 39 untranslated regions (UTRs) of target mRNAs [15]. miRNAs have been shown to play important roles in many biological processes including liver metabolism. For example, miR-122, abundantly expressed in liver, modulates protein metabolism in liver by targeting cationic amino acid transporter 1 [16]; it regulates the synthesis of fatty acids and cholesterol by repressing the expression of aldolase-A, 3-hydroxy-3-methylglutaryl-coenzyme A reductase, and AMP-activated protein kinase [17,18]. miR-33 is another miRNA involved in liver metabolism: it regulates cholesterol efflux and high-density lipoprotein metabolism by targeting ATP-binding cassette, subfamily A (ABC1), member 1 and ATP-binding cassette, sub-family G (WHITE), member 1 [19], and it reduces fatty acid degradation by targeting multiple genes involved in fatty acid b-oxidation [20].
RNA sequencing (RNA-seq) is a novel gene expression profiling technology based on high-throughput DNA sequencing. The benefit of RNA-Seq over other large-scale gene expression profiling methods is its ability to measure mRNA expression in a single assay and, at the same time, reveal new genes and transcripts [21,22,23]. Similarly, RNA-seq can be also used to identify novel miRNAs and detect differentially expressed miRNAs between samples [24].
As in mammals [3,7,25,26], GH has metabolic effects in chickens: GH regulates lipid metabolism in chicken adipose tissue [27] and chicken liver [28]. However, it is unclear whether GH has the same growth-stimulating effect in chickens as in mammals because exogenous GH treatment to chickens causes no growth responses and because plasma GH concentrations in chickens are in general not correlated with growth rates [29].
The mechanism by which GH regulates lipid metabolism in chicken liver is not clear. In this study, we determined the effects of GH on the expression levels of all mRNAs and miRNAs in primary hepatocytes from female chickens by RNA-seq. We analyzed the differentially expressed mRNAs or genes (DEG) and differentially expressed miRNAs (DEM) with multiple bioinformatics tools to correlate the DEG and the DEM to the physiological functions of GH. The main hypothesis to be tested in this study was that GH regulates liver metabolism in the chicken in part by regulating the expression of miRNAs that target mRNAs directly related to liver metabolism.

Culture of primary chicken hepatocytes
All procedures involving animals were approved by Changshu Institute of Technology Institutional Animal Care and Use Committee and conformed to the Guide for the Care and Use of Laboratory Animals of Jiangsu Province. All efforts were made to minimize suffering. Primary chicken hepatocytes were isolated from 4 female, 4-week-old Arbor Acres commercial chickens which were fasted 12 hours (h) before being anaesthetized by intraperitoneal injection of sodium thiopenthal (50 mg/kg) and anticoagulated by intraperitoneal injection of heparin (1750 U/ kg). Livers were isolated and the chickens were sacrificed by removal of the hearts. Hepatocytes were isolated from the livers as previously described [30]. Hepatocytes from individual chickens were cultured separately. Hepatocytes were plated in 24-well plates or 10-cm dishes at a density of 1.3610 6 cells/ml in Willam's E medium (Gibco, Grand Island, NY) supplemented with 5% chicken serum, 100 U/ml penicillin-streptomycin, 10 mg/ml insulin and 30 mmol/L NaCl in a humidified incubator at 37uC with 5% CO 2 . Twenty-eight h later, hepatocytes were serum starved for 8 h, followed by 12 h of treatment with 500 ng/ml chicken GH (chGH) (Prospec, Ness-Ziona, Israel) or an equal volume of PBS. Cells were lysed and total RNA was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA) following the manufacturer's directions.

Real time RT-PCR
Concentrations and quality of RNA samples were determined by NanoDrop ND2000 spectrophotometry (Thermo Scientific, Wilmington, DE) and formaldehyde-agarose gel electrophoresis. Five hundred ng of total RNA was reverse-transcribed to cDNA in a total volume of 10 ml using Takara PrimeScript RT reagent kit (Takara, Dalian, China). The relative expression levels of genes were quantified using SYBR Premix Ex Taq (Takara) on an Applied Biosystems 7500 Thermocycler (Applied Biosystems) according to the manufacturer's directions. In this analysis, GAPDH was used as an internal control. All reactions were run in duplicate. Data are means 6 SEM of 4 independent cell culture experiments (i.e., 4 chickens) and analyzed by student's t-test. The primers for real time RT-PCR were presented in Table S1.
Gene Ontology (GO) enrichment analysis was performed using GOslim (http://www.geneontology.org/page/go-slim-and-subsetguide). The data were presented as 2log (P-value). P-value = 12 where N is the total gene number in genome, n is the number of all DEG, M is the number of genes in genome that are involved in certain GO, and m is the number of DEG that are involved in certain GO. The signaling pathways were analyzed using the KEGG database.

Small RNA sequencing
The RNA samples used for small RNA sequencing were the same as for mRNA sequencing. Small RNA libraries for small RNA sequencing were constructed according to the Miseq small RNA sequencing protocol. Briefly, total RNA was ligated to the 39 adapter. Then it was ligated to the 59 adapter. The ligated RNA was reverse-transcribed to cDNA. The cDNA which had a 39 adapter and a 59 adapter was amplified by PCR. The PCR products with appropriate length were extracted from polyacrylamide gels to construct the small RNA libraries. The libraries of four samples were pooled in equal volumes after first being normalized to 10 nM. Sequencing of the small RNA libraries was performed by Personalbio using the Illumina Miseq system.

miRNA analysis
The raw reads from small RNA sequencing were treated by trimming the adapters and removing the low-quality sequences to obtain the clean reads. The clean reads with 15-30 nt were analyzed by counting and grouping identical sequences as unique reads. The unique reads were mapped to the chicken genome using Bowtie and BLAST searched against the ncRNA database Rfam (10.1) assess the quality of sequences and obtain ncRNA annotation. The sequence reads were first searched against the chicken miRNA database in miRBase (20.0) to identify known chicken miRNAs, then against the miRNA databases of other species in miRBase to identify chicken miRNAs homologous to known miRNAs in other species, and finally against the chicken genome to obtain their ,80 nt genomic sequences flanking the 59 or 39 end and analyzed using the program mireap (http:// sourceforge.net/projects/mireap/) to predict potentially novel chicken miRNAs and their precursors according to the miRNA biogenesis principle [31].
The reads of miRNAs from control (PBS treated) samples and GH treated samples were normalized using reads per million (RPM). RPM = Actual miRNA count/Total count of clean reads*1000000. When the RPM of a certain miRNA in one type of samples was zero, it was revised to 0.01, and if the RPM of a certain miRNA in both types of samples was less than 1, the miRNA would not be used in further DEM analysis. Fold change and P-value were calculated from the RPM. Fold change = RPM (chGH)/RPM (PBS). P-value formula: p xjy . The x and y represent normalized expression level, and the N1 and N2 represent total count of clean reads of a given miRNA in small RNA libraries of chGH and PBS groups, respectively. The miRNAs with fold change (RPM (chGH/PBS)) .1.5 or ,2/3 and P-value ,0.05 were identified as DEM [32]. Target genes of identified known chicken DEM were predicted in the DEG using miRanda with the TargetScan principle [30]. In brief, the target mRNAs of GH up-regulated miRNAs were predicted in GH down-regulated mRNAs; the target mRNAs of GH down-regulated miRNAs were predicted in the GH upregulated mRNAs. Then GO enrichment analysis of GHregulated mRNAs targeted by GH-regulated miRNAs was performed using the GOslim program.

mRNA expression profiling by RNA-seq
The RNA-seq generated 16,474,842 and 10,662,250 paired end reads from chGH-treated group and PBS-treated group, respec- Table 4. List of growth-related mRNAs whose expression in chicken hepatocytes was regulated by chGH.  Table 5. List of metabolism-related mRNAs whose expression in chicken hepatocytes was regulated by chGH.  (Table S2). Table 2 shows the top 10 up-regulated DEG and the top 10 down-regulated DEG.
To confirm the result of RNA-seq, we classified the genes into four groups: 1) fold change (RPKM (chGH/PBS)) .1.5 and P, 0.05, 2) fold change (RPKM (chGH/PBS)) ,2/3 and P,0.05, 3) fold change (RPKM (chGH/PBS)).1.5 and P.0.05, 4) fold change (RPKM (chGH/PBS)) ,2/3 and P.0.05. And real-time RT-PCR was performed on a total of 39 genes that were randomly selected from four groups. Among them, 21 genes that were not differentially expressed between chGH and PBS based on RNAseq were also found to be not differentially expressed by real-time RT-PCR, and 11 of 18 genes that were found to be differentially expressed by RNA-seq were also confirmed to be differentially expressed by real time RT-PCR (Table 3). Thus, our RNA-seq results were in general confirmed by real-time RT-PCR.

Functional analysis of differentially expressed mRNAs
To correlate the differentially expressed mRNAs with biological functions, we analyzed the functional bias of the DEG according to Gene Ontology enrichment. The analysis showed that the DEG included many genes involved in lipid binding, ion binding, transferase activity, oxidoreductase activity, immune system Table 6. List of signal transduction-related mRNAs whose expression in chicken hepatocytes was regulated by chGH. process, lyase activity, and lipid metabolic process. Other enriched GO terms included growth, cellular nitrogen compound metabolic process, carbohydrate metabolic process, and biosynthetic process ( Figure 1). Most of the KEGG Orthology (KO) terms were related to lipid metabolism (Table S3). Among the GH up-regulated genes, three were involved in the pathways controlling lipid metabolism: AKR1D1 was involved in primary bile acid biosynthesis and steroid hormone biosynthesis pathway, FABP1 was involved in PPAR and fat digestion and absorption pathway, and LPIN1 was involved in glycerolipid metabolism pathway. Among the GH down-regulated genes, only one was involved in the pathway related to lipid metabolism: PPAP2B was involved in glycerolipid metabolism and fat digestion and absorption pathway.

Differentially expressed mRNAs involved in
growth. We found that there were 6 GH-regulated genes involved in animal growth. Of them, CISH, ULK2, IRF8, and FOXP2 were up-regulated by chGH, and Bcl6 and ROBO1 were down-regulated by chGH (Table 4). Both Bcl6 and ROBO1 are negative regulators of cell growth [33,34]. This indicates that chGH may regulate chicken growth by both increasing the expression of genes that positively regulate growth and reducing the expression of genes that negatively regulate growth.

Differentially expressed mRNAs related to
metabolism. A total of 46 GH-regulated DEG were found to participate in protein, carbohydrate and lipid metabolism. Among them, 26 were related to protein metabolism as shown in Table 5, of which 20 were up-regulated and 6 were down-regulated by chGH, suggesting that GH plays an important role in protein metabolism in chicken hepatocytes. Further analysis showed that most of the 26 genes either promote protein synthesis or inhibit protein degradation, consistent with a previous study [7]. In addition to protein metabolism, 5 genes were related to carbohydrate metabolism, and all of them were up-regulated by chGH ( Table 5). The specific functions of these genes include carbohydrate binding, transferase activity, and gluconeogenesis. As shown in Table 5, 15 DEG were related to lipid metabolism. Among these 15 genes, 8 showed increased expression after chGH treatment and all have oxidoreductase activity, suggesting that GH may promote lipid oxidation in chicken hepatocytes. Most of the 7 down-regulated genes were related to unsaturated fatty acid and long chain fatty acid biosynthesis, suggesting that GH may inhibit lipid synthesis in chicken hepatocytes.

Differentially expressed mRNAs related to signaling
pathways. We also identified which GH-regulated genes might be involved in intracellular signaling. The analysis indicated that 7 GH-regulated genes have DNA binding activity, of which 5 were up-regulated and 2 were down-regulated by chGH, and 24 GHregulated genes were components of signal transduction pathways, of which 11 were up-regulated and 13 were down-regulated by chGH (Table 6). These pathways include the STAT and MAPK pathways, which are known to be activated by GH [35,36] and pathways that are not known to be activated by GH.

miRNA profiling by RNA-seq
Small RNA libraries of chGH-treated and PBS-treated hepatocytes were deep sequenced, generating 8,248,771 raw reads from PBS-treated hepatocytes and 5,232,471 raw reads from chGH-treated cells. Of these sequences, 5,240,376 and 2,144,128 were clean reads. A total of 3,124,826 clean reads and a total of 1,487,885 clean reads from PBS-treated cells and chGH-treated cells, respectively, were found to be between 15 and 30 nt in length, and these reads were considered as potential miRNAs.
The majority of the 15-30 nt long reads ranged from 18 to 24 nt, and reads of 22 nt long were more abundant than reads of other lengths, indicating that the distribution of the small RNA sequences was consistent with the length range of miRNAs ( Figure 2). In other words, most of the 15-30 nt long reads were likely miRNAs. A further analysis showed that most of the small RNAs could be mapped to the chicken genome. After the 15-30 nt reads were counted and grouped to generate unique reads, total reads or unique reads were compared to the sequences in Rfam database and classified. About 75% of total reads matched the sequences in Rfam. As shown in Figure 3, among the matched reads, ,80% of total reads represented miRNAs, and ,20% of unique reads represented miRNAs, indicating that the sequenced small RNA reads were enriched with miRNAs. The unique reads that matched miRNAs in Rfam were BLAST searched against chicken miRNAs in miRBase (20). Of 996 known chicken miRNAs in miRBase, 225 were identified in the present study. Among them, 219 known chicken miRNAs were identified from PBS-treated and 206 were identified from chGH-treated hepatocytes, which corresponded to 491 and 422 known chicken miRNA precursors, respectively (Tables 7 and S4). The remaining reads that did not match known chicken miRNAs were BLAST searched against all the miRNAs in other species, and 264 homologous miRNAs were identified. Of them, 259 were from PBS-treated cells and 234 were from chGH-treated cells (Table  S5). The small RNA reads that had no matches in Rfam were analyzed to predict novel miRNAs and their precursors. A total of 95 sequences were predicted to be novel miRNAs, of which 93 were sequenced from PBS-treated cells and 73 from chGH-treated cells (Table S6). We estimated the potential of the precursor sequences for these miRNAs to form stable stem-loop hairpin structures ( Figure S1). Interestingly, some of these sequences are related to known chicken miRNAs. For example, the predicted novel chicken miRNA gga-m0015-3p is complementary to the known chicken miRNA gga-miR-126-5p, the 39 portion of gga-m0016-5p overlaps with the 59 portion of gga-miR-219b, and the 59 portion of gga-m0028-3p overlaps with the 39 portion of gga-miR-3525 ( Figure S2).

GH-regulated miRNAs
We estimated the expression levels of miRNAs based on their read numbers. The 10 most abundantly expressed miRNAs in chicken hepatocytes were listed in Table S7. Among the identified known chicken miRNAs, 17 were identified as DEM between PBS-treated and chGH-treated hepatocytes, of which 16 were upregulated and 1 was down-regulated by chGH. Among the identified chicken miRNAs homologous to miRNAs in other species, 21 were DEM between PBS-treated and chGH-treated hepatocytes, of which 15 were up-regulated and 6 were downregulated in chGH-treated cells. Among the predicted novel chicken miRNAs, 7 were DEM between PBS-treated and chGHtreated cells, of which 5 were up-regulated and 2 were downregulated by chGH (Table 8).

Functional analysis of GH-regulated miRNAs
The major function of miRNAs is to down-regulate the expression of target mRNAs. To investigate the function of GHregulated miRNAs, we determined which of the GH-regulated mRNAs could be targeted by the GH-regulated miRNAs. Among the GH down-regulated mRNAs, 32 were predicted as targets of GH up-regulated miRNAs. Among the GH up-regulated mRNAs, 12 were predicted as target genes of GH down-regulated miRNAs ( Table 9). A GO enrichment analysis showed that these miRNA target genes were enriched in lipid metabolism, lipid binding, cell motility, and small molecule metabolic process (Figure 4). The most significant GO term was the lipid metabolic process. As shown in Table 10, the GH up-regulated miRNAs targeted 7 genes related to lipid metabolism and the GH down-regulated miRNAs were predicted to target one gene related to lipid metabolism. Among these GH-regulated miRNAs, miR-15b had more predicted target genes related to lipid metabolism than had any other miRs. Several GH-regulated miRNAs were predicted to target the same genes (Table 10).

Discussion
In this research through RNA-seq, we found a total of 164 DEG between chGH-treated and untreated chicken hepatocytes. Functional analyses showed most of the GH-regulated genes are involved in liver metabolism, indicating that GH regulates liver metabolism. Lipid metabolism is different between birds and mammals. In mammals lipogenesis occurs in liver, adipose tissue, and mammary gland, whereas in birds this occurs mostly in liver [37]. Lipid metabolic process was identified as a significant GO term among chGH-regulated genes, indicating that GH might play an important regulatory role in lipid metabolism in chicken hepatocytes. In the study of human adipose tissue, Zhao et al. (2011) found the GH-regulated DEG include those that stimulate triglyceride (TG)/ free fatty acid (FFA) cycle, and they also found a new TG hydrolase gene called patatin-like phospholipase domain containing 3 (PNPLA3), which could promote TG hydrolysis [38]. Comparing our DEG with those in the above report, we found there was little overlap, indicating the mechanism of GH regulation of lipid metabolism is different between birds and mammals, or between adipose tissue and liver. The chicken PNPLA4 identified in our study is homologous to PNPLA3 in mammals, suggesting this gene may play a role in regulating TG hydrolysis in the chicken liver.
In addition to genes that are involved in animal growth, many of the 164 DEG are involved in biosynthetic process (Table S8). This suggests that GH stimulates body growth in chickens not only by stimulating cell proliferation and hypertrophy but also by stimulating nutrient accumulation.
Recent research indicated that multiple signaling pathways are activated by GH, including the JAK2-STAT pathway, the MAPK-ERK1/2 pathway and the PI3K-AKT pathway [12]. In this study, 31 DEG are transcription factors or signaling molecules. Some of them belong to the pathways known to be activated by GH [35,36] while others belong to pathways that are not known to be activated by GH. Further study of the DEG showed that some of them, for example CAMK1D, had both DNA binding activity and tyrosine kinase activity. Another interesting GH-regulated transcription factor is BCL6. Studies showed that the expression of BCL6 was controlled by GH through STAT5, and as a transcription factor BCL6 could mediate downstream gene expression [39,40]. We found that chicken BCL6 gene has 11 potential STAT5 binding sites (data not shown), implying it could be regulated by chGH through STAT5.
The KEGG pathway analysis of GH-regulated DEG showed the KO terms were mainly related to metabolism, in particular, lipid metabolism. Two DEG were involved in fat digestion and absorption. One of them is FABP1, which is an important lipidrelated gene [41] and has been shown to be up-regulated by GH in young chickens [42]. The other gene PPAP2B is involved in the glycerolipid metabolism pathway.
Among the identified miRNAs in this study, gga-miR-148a is the most abundant miRNA in chicken hepatocytes. This is consistent with the previous report, in which miR-148a was identified as the most abundant miRNA in porcine livers of different breeds [43]. This result indicated miR-148a might play important roles in liver. From many previous reports, miR-122  was known as the most abundant miRNA in liver. But in our study gga-miR-122-5p (previously named gga-miR-122) was identified as the fifth abundant miRNA in primary chicken hepatocytes, similar to the report by Li et al. which found miR-122 was the fourth abundant miRNA in porcine liver [43].
The differentially expressed miRNAs identified in this study might play various important roles in chGH regulation of growth, development and metabolism. miR-223 is expressed higher in the liver of Large White pig (lean type) than in the liver of Tongcheng pig (fatty type) [44]. In this study, we found gga-miR-223 was upregulated by chGH, indicating it might play a role in lipid metabolism in chicken liver. miR-193a and -190 are expressed at higher levels in the liver of Erhualian pig than in the liver of Large White pig (The two breeds have different rate of lipid metabolism) [43]. In this study gga-miR-193a and 190 were also up-regulated by chGH, indicating they may play roles in lipid metabolism in chicken liver too. miR-15b was another up-regulated miRNA by the chGH treatment. This miRNA in pigs is expressed at a lower level in the liver of Erhualian pig than in the liver of Large White pig, opposite to miR-193a and -190 [43], implying it may play different roles between chickens and pigs. In male rat livers, miR-451 and miR-29b are down-regulated by GH [45]. gga-miR-451 was too down-regulated by GH in chicken hepatocytes in our study, indicating gga-miR-451 might mediate GH regulation of chicken liver metabolism but gga-miR-29b played different roles in GH regulation of liver metabolism comparing to miR-29b in rat.
The expression of miR-122 is not changed by GH in male rat liver [45]. We found the same for gga-miR-122 in this study, indicating this miRNA might not be involved in GH regulation of liver metabolism. In this study we also found that some of GHregulated DEG were potential targets of GH-regulated DEM. GO enrichment analysis showed that these DEG were enriched with genes involved in lipid metabolism, indicating that GH may regulate lipid metabolism through miRNAs. Among the GHregulated DEM, miR-15b was predicted to target more GHregulated DEG related to lipid metabolism than any other DEM, including ACSL3, LIPG, PLCD1, PPAP2B and STAR. miR-15b expression was increased in NAFLD models and it lead to inducing the storage of intracellular triglyceride [46]. When ACSL3 increases, triglycerides and free fatty acids reduce [47]. Knock down of ACSL3 decrease hepatic lipogenesis [48]. LIPG has substantial phospholipase activity, and it can reduce plasma concentrations of HDL cholesterol [49]. PLCD1 can catalyze the hydrolysis of membrane lipid phosphatidylinositol 4,5-bisphosphate into second messengers inositol 1,4,5-trisphosphate and diacylglycerol [50]. PPAP2B bioactive lysophospholipids, including lysophosphatidic acid and sphingosine-1-phosphate, and thereby terminates their signaling effects [51]. STAR plays a critical role in the rapid translocation of cholesterol across the outer and inner mitochondrial membranes [52]. It also plays a key role in steroidogenesis by enhancing the metabolism of cholesterol into pregnenolone [53]. This finding indicates miR-15b might be a  Table 9. GH-regulated chicken genes predicted to be targeted by GH-regulated chicken miRNAs. very important miRNA that mediates the effect of GH on lipid metabolism in chicken liver. In mammals, the secretion of GH is characterized by an ultradian rhythm and sexual dimorphism. In the male rat, the major bursts of GH secretion occur at regular 3.3 h intervals separated by intervening trough periods with very low or undetectable basal GH levels [54,55]. GH is also secreted in a pulsatile manner in chickens [56,57]. In this study, hepatocytes from chickens were treated with one concentration of GH in culture. Therefore, a limitation of the study is that the observed GH-regulated mRNA and miRNA expression profiles may not represent those that occur in the chicken liver in vivo.

Conclusions
This study reveals that GH regulates the expression of many mRNAs involved in metabolism in female chicken hepatocytes, which suggests that GH plays an important role in regulating hepatic metabolism besides somatic growth. This study suggests that GH regulates lipid metabolism in female chicken liver in part . GO enrichment analysis of chGH-regulated mRNAs predicted to be targeted by GH-regulated miRNAs. The target mRNAs for GH up-regulated miRNAs were predicted in the GH down-regulated mRNAs; the target mRNAs for GH down-regulated miRNAs were predicted in the GH up-regulated mRNAs. doi:10.1371/journal.pone.0112896.g004 by regulating the expression of miRNAs that target the mRNAs involved in lipid metabolism. Figure S1 Structures of predicted novel miRNAs. For each novel miRNA, the precursor name, position, strand, length and energy are shown in the first line. The sequence of precursor and total reads are described in the second line. The brackets in the third line denote the secondary structure. The sequence of major unique reads, novel miRNA name and the number of reads are shown below with '*'. Then the sequence of each unique reads and the number of reads are shown below with '2'. (TXT) Figure S2 Relationship between predicted novel miR-NAs and known chicken miRNAs. For each novel miRNA, the name is listed in the first line. The sequence of known precursor and its name are described in the second line. The sequence of known miRNA and its name are shown below with '*'. Then the sequence of novel miRNA and its name are shown below with '2'. Note: The sequence of gga-m0015-3p is complementary to gga-miR-126-5p. (TXT)