- Research article
- Open Access
Identification and characterization of long noncoding RNAs and mRNAs expression profiles related to postnatal liver maturation of breeder roosters using Ribo-zero RNA sequencing
BMC Genomics volume 19, Article number: 498 (2018)
The liver is mainly hematopoietic in the embryo, and converts into a major metabolic organ in the adult. Therefore, it is intensively remodeled after birth to adapt and perform adult functions. Long non-coding RNAs (lncRNAs) are involved in organ development and cell differentiation, likely they have potential roles in regulating postnatal liver development. Herein, in order to understand the roles of lncRNAs in postnatal liver maturation, we analyzed the lncRNAs and mRNAs expression profiles in immature and mature livers from one-day-old and adult (40 weeks of age) breeder roosters by Ribo-Zero RNA-Sequencing.
Around 21,939 protein-coding genes and 2220 predicted lncRNAs were expressed in livers of breeder roosters. Compared to protein-coding genes, the identified chicken lncRNAs shared fewer exons, shorter transcript length, and significantly lower expression levels. Notably, in comparison between the livers of newborn and adult breeder roosters, a total of 1570 mRNAs and 214 lncRNAs were differentially expressed with the criteria of log2fold change > 1 or < − 1 and P values < 0.05, which were validated by qPCR using randomly selected five mRNAs and five lncRNAs. Further GO and KEGG analyses have revealed that the differentially expressed mRNAs were involved in the hepatic metabolic and immune functional changes, as well as some biological processes and pathways including cell proliferation, apoptotic and cell cycle that are implicated in the development of liver. We also investigated the cis- and trans- regulatory effects of differentially expressed lncRNAs on its target genes. GO and KEGG analyses indicated that these lncRNAs had their neighbor protein coding genes and trans-regulated genes associated with adapting of adult hepatic functions, as well as some pathways involved in liver development, such as cell cycle pathway, Notch signaling pathway, Hedgehog signaling pathway, and Wnt signaling pathway.
This study provides a catalog of mRNAs and lncRNAs related to postnatal liver maturation of chicken, and will contribute to a fuller understanding of biological processes or signaling pathways involved in significant functional transition during postnatal liver development that differentially expressed genes and lncRNAs could take part in.
Liver is the most important metabolic organ exhibiting both endocrine and exocrine properties . However, most of the liver functions are not mature at birth and many changes are extensively remodeled during postnatal liver development to rapidly adapt and perform adult functions . Functional adaption during postnatal liver development is relied on finely programmed alteration of gene expression . Cui et al.  and Peng et al.  reported that the Cytochrome P450 gene isoforms and their alternative transcripts were closely related to postnatal liver development of mouse. Some other important genes, such as yes-associated protein, aryl, estrogen, and hydroxysteroid sulfotransferases, also could regulate the process of postnatal liver development [6, 7]. Using recently developed RNA sequencing techniques, several genes have been found to involve in postnatal liver development of mouse . However, the differentially expressed genes in chicken liver, which involved in postnatal liver development process and the adaption of mature hepatic metabolic and immune functions, have not been reported before.
It is well-known that protein-coding genes account for only approximately 1.5% of the genome, meaning that most of transcripts have little translation potential [9,10,11], and these non-coding RNAs could play crucial roles in regulating target gene expressions [11, 12]. As the most important non-protein coding transcripts longer than 200 base pairs, long noncoding RNAs (LncRNAs) are extensively expressed in several species of animals [8, 13,14,15]. Studies on lncRNAs have been shown to control several levels of the gene expression program, including DNA methylation , mRNA expression and degradation , and the effective concentration of miRNAs , by acting as signal, decoys, guides, scaffolds , or competing endogenous RNA . These results have highlighted the regulatory roles of lncRNAs in regulating epigenetic modification and gene expression. However, compared with the extensive characterization of DNA methylation and alternative pre-mRNA splicing related to fetal-to-adult liver maturation [19, 20], the roles of lncRNAs in regulating postnatal liver maturation need to pay more attentions.
RNA-seq technology has rapidly developed to enable discovery and analysis of non-coding RNA, and differential methods have been developed to identify novel lncRNAs using RNA-seq data. With oligo (dT) selection of poly (A)+ mRNA, recent studies in mice have focused on the changes of lncRNAs or mRNAs related to the regulation of postnatal live development and functional changes [8, 21, 22]. However, not all of the lncRNAs contains 3’polyadenylation . Thus it could ignore some lncRNAs information by using poly (A) selection methods to prepare cDNA library. Compared with poly (A) selection Sequencing, Ribo-Zero RNA Sequencing, which provided equivalent rRNA removal efficiency and coverage uniformity but exhibited a highly technical reproducibility, can help study poly (A)− mRNA, immature transcripts, and the lncRNAs [24, 25]. In contrast to the small RNAs, which are highly conserved and involved in transcriptional and posttranscriptional gene silencing through specific base pairing with their targets, lncRNAs are poorly conserved and regulate gene expression by diverse mechanisms, suggesting primarily lineage-specific functions .
Hence, we hypothesized that postnatal liver maturation of chicken were related to their differential lncRNA expression profiles. Herein, we investigated the expression profiles of lncRNAs and mRNAs related to postnatal liver maturation of chicken by Ribo-Zero RNA-Seq  by using 6 livers transcriptome libraries from arbor acres roosters, which is one of the most used broiler breeder roosters. Taken together, these expression profiles could especially clarify the changes of mRNAs related to postnatal liver maturation and the roles of lncRNAs in postnatal liver development, metabolism, and other liver functions.
RNA sequencing and identification of mRNA and lncRNAs in chicken liver
To systematically identify mRNAs and lncRNAs expressed in the mature and immature livers of breeder roosters, we generated six RNA expression profiles of liver tissues with an average of 87 million 150 bp paired end raw reads. After initial processing, the average of 58 million valid reads were obtained from each sequencing library (Table 1). By using TopHat , nearly 86% of the reads were mapped to the Gallus gallus reference genome (Galgal5). The mapped sequences in each library were assembled and annotated using the StringTie . In the present study, a total of 20,158, 19,726, 19,751, 19,351, 19,543, and 19,687 unique genes from the six libraries were respectively identified (Table 1 and Additional file 1). According to length and coding potentials (see methods), we identified 1918, 1927, 1932, 2019, 2042, and 2022 unique lncRNAs from those six libraries (Table 1 and Additional file 2). These 2220 unique lncRNAs were distributed across the chromosomes in Gallus gallus related to the length of the chromosomes (Fig. 1a). According to the locations of lncRNAs in the genome, 80 antisense, 49 sense, 1922 intergenic, 32 intronic, and 137 intron lncRNAs were identified among our identified lncRNAs (Fig. 1b). Moreover, the sequence information of all identified lncRNAs were listed in Additional file 3.
In the present study, the average length of lncRNAs was 1718 bp compared to more than 2936 bp for protein-coding genes, which indicated that lncRNAs were shorter than protein coding transcripts (Fig. 1c). Genes of lncRNAs tend to contain fewer exon: lncRNAs identified in our study had only 1.35 exons per transcript on average while protein-coding genes had averaged 10.15 exons (Fig. 1d). Furthermore, the lncRNAs in chicken livers tended to be shorter in length of open reading frame than protein coding genes (Fig. 1e and f). Overall, lncRNAs identified in our study were characterized with fewer exons, shorter transcript length, and significantly lower expression levels, compared to protein-coding genes, which is consistent with previous studies in other species as well as in chickens [30,31,32].
Developmental changes of protein-coding genes and lncRNAs in postnatal liver development of chicken
We gathered 1570 differentially expressed mRNAs between immature and mature livers meeting the criteria of P < 0.05 and log2foldchange > 1 or < − 1. Of these, 381 genes were up-regulated and 1189 mRNAs were down-regulated in the mature livers relative to the immature livers in breeder roosters (Fig. 2a and Additional file 4). In addition, 214 differentially expressed lncRNAs between immature and mature livers meeting the criteria of P < 0.05 and log2foldchange > 1 or < − 1 were obtained. Compared with the immature livers, 34 up-regulated lncRNAs and 180 down-regulated lncRNAs were obtained in mature livers (Fig. 2b and Additional file 5).
The expression levels of ten randomly selected lncRNAs and mRNAs were determined by quantitative real-time PCR. The results confirmed that these lncRNAs and mRNAs were expressed at both mature and immature livers (Fig. 3) and showed differential expression at different stages. In addition, the qRT-PCR confirmed that the expression patterns of these lncRNAs and mRNAs were consistent with their expression levels calculated from the RNA-seq data (Fig. 3).
Enrichment analysis of differentially expressed mRNAs
GO enrichments  of the differentially expressed mRNAs were categorized into 519 functional annotations that met the criteria of P < 0.05 (Additional file 6). Results of the biological process analysis revealed that these differentially expressed mRNAs could mainly regulate the biological process relative to metabolism, immune response, oxidation-reduction process, as well as some process related to the growth and functional maturation of liver, including the cell division, proliferation, apoptotic, and cell cycle. Moreover, analyses based on cellular components and molecular functions showed that differentially expressed mRNAs were also involved in the aforementioned process relative to the growth and functional maturation of liver. Further KEGG pathway  analyses revealed that these mRNAs were mostly involved in several pathways affecting the metabolism of livers, including the amino acid metabolism, glycometabolism, and lipid metabolism, as being listed in Table 2. Furthermore, pathways relative to the growth of liver, such as cell cycle pathway, were also found to take part in the postnatal liver maturation in the present study.
Cis-regulatory roles and trans-regulatory roles of differentially expressed lncRNAs in postnatal liver maturation of chicken
To investigate the possible functions of the lncRNAs, we predicted the potential targets of lncRNAs in cis-regulatory and trans- regulatory relationships. In the present study, 4 and 214 differentially expressed lncRNAs could respectively cis- and trans-regulate 4 and 1394 differentially expressed mRNAs in the present study (Additional file 7). GO analysis based on these cis- and trans- regulated targets was performed and 174 significant enriched (P < 0.05) GO terms were obtained (Additional file 8). As the main function of mature liver after birth, several biological processes and molecular function related to metabolic and immune process have also been detected. Results also revealed that these differentially expressed lncRNAs could cis- or trans- regulate targets involved in liver development, including liver development process, cell adhesion, cell proliferation, and apoptotic processes. KEGG pathway analyses showed that a total of 153 KEGG pathways were annotated based on the cis- and trans- regulated differentially expressed targets of differentially expressed lncRNAs. Of these, 25 KEGG pathways were annotated with P < 0.05 (Table 3), most of which were related to the regulation of mature hepatic metabolic and immune functions. Moreover, some pathways, which have been proved to regulate liver development, have also been detected in the KEGG analyses of differentially expressed lncRNAs, such as the cell cycle, Notch signaling pathway, Wnt signaling pathway, and Hedgehog signaling pathway. These analyses indicated that the differentially expressed lncRNAs and their annotated pathways should play cis- or trans- regulated roles in postnatal liver maturation. Moreover, the interaction between lncRNAs and their cis- or trans- regulated target genes which involved in different GO terms and KEGG pathways were also respectively listed in Additional file 7.
Integrated analyses of differentially expressed lncRNAs and mRNAs involved in postnatal liver development
Several differentially expressed protein-coding genes involved in liver development, such as SNX, CYP7A1, CYP39A1, HNF4α, and IGF2BP1, could be cis-regulated or trans-regulated (Table 4) by several differentially expressed lncRNAs. Furthermore, we can find that several key genes of some important pathways involved in liver development, including the cell cycle, Notch signaling pathway, Wnt signaling pathway, and Hedgehog signaling pathway, could also be regulated by the lncRNAs in cis-role or trans-role (Additional file 9). Overall, we suspected that these lncRNAs most probably participated in the postnatal liver development, although its underlying mechanisms require additional investigations.
Breeder roosters and their offspring broilers (Gallus gallus) are famous for its high feed conversion efficiency, which means a high efficient metabolic process in chicken liver . There are several serious metabolic diseases occurred during the feeding of breeder roosters, such as fatty liver and ascites syndrome , which are induced by incorrect metabolic regulation and could further influence the usability of breeder roosters. It’s important to clarify the postnatal liver maturation process and the hepatic function changes after birth, so that the hepatic metabolic or immune condition of breeder roosters could be better regulated and the breeder roosters could be used maximally. However, compared with abundant researches on liver development of human or mice [8, 21, 22, 37,38,39,40], limited researches have focused on the postnatal liver maturation process in chicken. In the present study, we identified the differentially expressed mRNAs and lncRNAs in immature and mature livers of breeder roosters by using RIBO-zero RNA-Sequencing, which is the first report to systematically identify mRNAs and lncRNAs expression profiles during chicken postnatal liver development.
Postnatal liver development from the newborn to the adult stage consists of a series of exquisitely regulated and orchestrated changes in the expression of many genes. In the previous studies, we can easily find that several differentially expressed mRNAs have been found to be involved in liver maturation in other species of animals, such as sorting nexin (SNX) in zebra fish , insulin like growth factor 2 mRNA binding protein 1 (IGF2BP1) in human , as well as hepatocyte nuclear factor 4-alpha (HNF4α)  and cytochrome P450 gene isoforms [4, 5] in mice. These genes could be detected in our identified differentially expressed mRNAs, which indicted that these differentially expressed mRNAs have important roles in regulating livers development of chickens, as being shown in other animal species. Further GO and KEGG analyses indicated that differentially expressed mRNAs regulated many metabolic and immune processes in liver, which are related to normal mature liver functions [1, 44, 45] and were consist with the results of previous studies in mice [8, 21, 22]. Except for the metabolic and immune changes with the maturation of liver, we found that there were several GO terms involved in the oxidation-reduction process and antioxidant ability, which was also the main function of mature liver . Moreover, the majority of genes in the liver cell proliferation process, such as the genes involved in cell cycle, DNA replication, cytokine binding, cell proliferation, and cell division(including cytokinesis, chromosome segregation, and mitotic nuclear division), displayed remarkable changes in expression during chicken liver development, which highlighted the importance of gene products in postnatal liver growth process [47, 48]. Specifically, the cell cycle pathway have been implied in the liver regeneration process [49, 50]. The liver regeneration is the compensatory growth of the liver, which indicated that the growth of liver could also be influenced by the cell cycle pathway. Sadler et al.  have reported that the uhrf1 gene, a cell cycle regulator, is required for physiologic liver growth in both embryos and adults in zebrafish. The uhrf1 was also been detect as the differentially expressed genes in the present study. These results again proved that the cell cycle pathway were involved in the postnatal liver growth process. To sum up, GO and KEGG analyses based on differentially expressed mRNAs indicated that livers underwent hypertrophic growth and maturation via large-scale changes in metabolic and immune functions after birth.
LncRNAs are a group of endogenous RNAs involved in developmental and physiological processes [52,53,54]. We obtained 300 up-regulated lncRNAs and 322 down-regulated lncRNAs in mature livers. These lncRNAs may have specific biological roles in postnatal liver development in chickens. Several recent studies have proved that lncRNAs could play crucial roles in liver development by using the RNA sequencing [8, 55, 56]. Compared with these previous studies, lncRNAs with or without poly (A) tails were obtained using RIBO-zero RNA-Sequencing in the present study. More types of lncRNAs, including those of sense, antisense, intronic and intergenic lncRNAs, were identified, while some previous studies only obtained the information of long intergenic non-coding RNAs. However, limit research proved the roles of single lncRNA in liver development process. Only the lncRNA-LALR1 were proved to enhance hepatocyte proliferation by promoting progression of the cell cycle, and further accelerate mouse hepatocyte growth and cell cycle progression during liver growth . In addition, lncRNA-LALR1 facilitated cyclin D1 expression through activation of Wnt/β-catenin signaling by way of suppression of Axin1 . Another study revealed that a set of lncRNAs highly correlated with expression of cytochrome P450, which alter their expressions during liver development and have critical functions in liver to metabolize xeno−/endo-biotics . .Therefore, the differentially expressed lncRNAs reported in the present study can also be considered as important novel regulators of chicken postnatal liver development process.
On the one hand, most evidence suggests that the expression of lncRNAs can regulate and have high correlations with expression of neighboring mRNAs in animals, microorganism, and plants [42, 59, 60]. On the other hand, many lncRNAs can also function in trans mode to target gene loci distant from where the lncRNAs are transcribed [11, 61]. In this study, the cis- and trans-regulated targets of the differentially expressed lncRNAs were obtained. There were several differentially expressed mRNAs related to liver development, such as SNX , IGF2BP1 , CYP7A1, CYP39A1 [4, 5], and HNF4α , could be regulated by differentially expressed lncRNAs in either cis or trans roles. These findings indicated that lncRNAs could cis- and trans-regulated the protein-coding genes associated with postnatal liver maturation. Further GO and KEGG analyses for trans- and cis-regulatory roles of differentially expressed lncRNAs were performed and found these lncRNAs could regulate the cell proliferation, cell cycle, as well as several liver function, such as metabolism, immunity, and antioxidant, which was consist with the results of functional annotation of differentially mRNA. These results indicated that the hepatic mature function could be regulated by the differentially expressed lncRNAs found in our study. We further found that the differentially expressed lncRNAs could regulate several pathways that have proven to be involved in liver development, such as Notch pathway, Hedgehog pathway, adherens junction pathway, and Wnt signaling pathway. In the previous studies, Zong et al.  and Kodama et al.  have proven that the Notch pathway were involved in the liver development by regulating biliary differentiation. Tanimizu and Miyajima  suggested that Notch signaling could control hepatoblast differentiation by altering the expression of liver-enriched transcription factors. As for hedgehog pathway, it is well known for its mitogenic and morphogenic functions during development . Reactivation of Hedgehog, a signaling pathway that controls hepatic progenitor cell fate and tissue construction, have also been linked to the regulation of adult liver tissue homeostasis, repair, and development [66, 67]. Further research have also proven that the Notch and Hedgehog pathways could interact to control the fate of hepatic key cell types involved in adult liver repair . Moreover, adherens junction pathway and Wnt signaling pathway were been linked to liver development as well [69,70,71]. Therefore, the differentially expressed lncRNAs involved in these above pathways (Additional file 9) could play important roles in postnatal liver development and were worthy of further research to illuminate their roles in liver development.
In conclusion, we firstly obtained high-quality mRNA and lncRNA expression profiles in chicken liver based on a RIBO-zero RNA-seq approach and have identified differentially expressed mRNAs and lncRNAs related to postnatal liver maturation of chicken. Moreover, we firstly reported the functional annotation of differentially expressed mRNAs and found that these some genes potentially play an important role in the development of chicken liver and the regulation of the functions of mature livers. Moreover, bioinformatics analysis suggests that some lncRNAs are involved in important biological processes and pathways associated with liver development such as cell cycle pathway, Notch signaling pathway, Hedgehog signaling pathway, and Wnt signaling pathway, and also could play an important role in regulating the gene expression of mature liver functions. Our results not only reveal new information regarding the development of chicken liver but also provide a broad and novel vision for future research at the molecular level in chicken.
Animals and sample collections
All experimental protocols and animals’ managements in the study were approved by the Institutional Animal Care and Use Committee (IACUC) of the Northwest A&F University (Yangling, Shaanxi, China). Five adult healthy Arbor Acres breeder roosters (40 week of age) and Twenty one-day-old healthy Arbor Acres breeder roosters were collected from Experimental Teaching Center of Animal Science of the Northwest A&F University (Yangling, Shaanxi, China). Here, the Arbor Acres breeder roosters in 40 week age were both somatic and sexual matured, which ensured their livers were assuredly matured. All 40-week-old breeder roosters were kept in an environmentally controlled henhouse with double-floor metabolism cages and exposed to a 16 h photoperiod. Water was available ad libitum and food was available according to the feeding standard of Arbor Acres breeder roosters (Aviagen, Alabama, USA).
These randomly selected roosters were fed deprived for 12 h, then euthanized by exsanguination after intravenous (IV) administration of 3% sodium pentobarbital (25 mg/kg; Sigma, USA) and immediately dissected. All efforts were made to minimize animals’ suffering. The whole left side livers were collected into Eppendorf tubes, and frozen immediately in liquid nitrogen. For each liver sample we gathered, the liver was grinded and homogenized using liquid nitrogen, and then all the homogenized liver samples were stored at − 80 °C until be analyzed. The livers of three adult roosters and twelve one-day-old roosters were selected randomly for RNA isolation and next generation sequencing analyses; furthermore, all livers samples from five adult roosters and twenty one-day-old roosters were used to extract total RNA and perform quantitative RT-PCR validation.
RNA isolation and sequencing
Total RNA from 3 adult roosters’ and 12 one-day-old rooster’ livers for RNA sequencing were extracted using Trizol reagent (Invitrogen, CA, USA) according to the manufacturer’s procedure. Specifically, the DNaseI was used during the RNA isolation process to avoid contamination with genomic DNA. The quantity and purity of total RNA were analyzed by a NanoDrop® ND-1000 spectrophotometer (Thermo Scientific, MA, USA), and integrity of RNA was accessed with Bioanalyzer 2100 and RNA Nano6000 LabChip Kit (Agilent, CA, USA). Only samples that had the OD260/280 > 1.8, OD260/230 > 2.0, and the RNA Integrity Number (RIN) > 7.0 were used for further sequencing (Additional file 10). Four RNA samples of one-day-old chicks’ liver were mixed equally together as a pooled RNA sample according to the purity of total RNA, and the RNA samples of adult chickens’ livers were directly used for library construction. In total, we gathered three RNA samples from three adult roosters and three pooled RNA samples from 12 one-day-old roosters for further library construction.
Approximately 3 μg of total RNA was used to prepare an LncRNA library. According to protocol of Epicentre Ribo-zero™ Gold Kit (Illumine, San Diego, USA), ribosomal RNA was removed and the rRNA-depleted RNA (Poly A+ and Poly A− RNA) were collected [24, 25, 42]. Subsequently, high strand-specificity libraries were generated using the rRNA-depleted RNA and a NEBNext Ultra Directional RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA) following the manufacturer’s recommendations. Briefly, the rRNA-depleted RNA was fragmented using divalent cations under elevated temperature in NEBNext. First-strand cDNA was synthesized using random hexamer primers and M-MuLV reverse transcriptase (RNase H−). Subsequently, second-strand cDNA synthesis was performed using second-strand synthesis reaction buffer, DNA polymerase I, and RNase H. Remaining overhangs were converted into blunt ends by exonuclease/polymerase activity. After adenylation of the 3′ ends of the DNA fragments, NEBNext adaptors with hairpin loop structures were ligated to the fragments to prepare them for hybridization. To select cDNA fragments that are 150–200 bp in length, the fragments in each of the library were purified with an AMPure XP system (Beckman Coulter, Brea, CA, USA). Then 3 μl USER Enzyme (NEB, Ipswich, MA, USA) was used with size-selected, adaptor-ligated cDNA at 37 °C for 15 min followed by 5 min at 95 °C before PCR. The qPCRs were performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers, and Index (X) Primer. The PCR products were purified (AMPure XP system) and library quality was assessed on an Agilent Bioanalyzer 2100 system. Clustering of the index-coded samples was performed on a cBot Cluster Generation System using a TruSeq PE Cluster Kit v3-cBot-HS (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. After cluster generation, the paired-end sequencing (2*150 bp) were performed on an Illumina Hiseq2500 at the LC-BIO (Hangzhou, China).
Reads mapping and transcriptome assembling
The 150 bp paired-end raw reads were firstly processed through FastQC to obtain the clean data, by removing the reads that contain sequencing adapter contaminations or poly-N and the low quality reads whose Q value were less than 20. At the same time, Q20, Q30 and GC content of the clean data were calculated. The clean reads from six cDNA libraries were merged and mapped to the Gallus gallus 5 (http://www.ensembl.org/Gallus_gallus/Info/Index) using TopHat and Bowtie v2.0.6 [28, 72]. The mapped reads of each sample were assembled using StringTie .Then, all transcriptomes were merged to reconstruct a comprehensive transcriptome using perl scripts. After the final transcriptome was generated, StringTie  and Ballgown  were used to estimate the expression levels of all transcripts. Specifically, ribosomal RNA in sequencing data was removed. We aligned all reads to ribosomal RNA of chicken download from ensemble by bowtie2, then those aligned reads were removed from fastq files.
Coding potential and identification of lncRNAs
The known protein-coding transcripts and the transcripts whose length were smaller than 200 bp were firstly removed, and the remaining unknown transcripts were used to screen for putative lncRNAs. Then, the coding potential for the remaining transcripts was calculated by CNCI  and CPC . A transcript was deemed to be lncRNA if the coding potentials were scored to be less than − 1 by CPC software and the coding potentials were scored to be less than 0 by CNCI software, which suggest that this transcript has no capacity of coding for proteins. Briefly, those candidate transcripts whose length were longer than 200 nt and the intersection between both coding prediction tools CPC and CNCI were deemed to be lncRNA.
Analysis of differential expression patterns
Expression levels of all transcripts, including putative lncRNAs and mRNAs, were quantified as FPKM using the StringTie . Based on Negative binomial distribution, differential gene expression was determined using DESeq with a P value < 0.05 and log2foldchange > 1 or < − 1 .
GO and KEGG enrichment analysis of differentially expressed mRNA
Functional annotation of differentially expressed mRNA were performed based on Gene Ontology (GO) database and enriched pathways were analyzed using Kyoto Encyclopedia of Genes and Genomes database (KEGG).GO enrichment analysis for the screened differentially expressed mRNAs was carried out using GOseq platform . The KEGG pathway enrichment analysis for the differentially expressed mRNAs was performed by using KOBAS software . In these two analyses, P < 0.05 were defined as significantly enriched GO terms or KEGG pathways.
Target gene prediction and functional enrichment analysis
To explore the function of lncRNAs, we predicted the target genes of lncRNAs in cis and trans. Cis-acting lncRNAs target neighboring genes . We searched for coding genes 100-kb upstream and downstream of all the identified differentially expressed lncRNAs by using python script, according to the previously described method . Some lncRNAs play trans-roles in regulating target genes through complementary base pairing (part region, default max trace back 50 nt). In the present study, based on the sequences of our identified differentially expressed mRNAs and lncRNAs, trans-regulated targets of the lncRNAs were identified according to the free energy (< − 11 kcal/mol) that needed to form the secondary structure using RIsearch . Moreover, those lncRNAs and genes should house in different chromosomes. Briefly, the parameters of RIsearch were set as: “ext_len=50, ext_penalty=30, RIsearch_energy= -11”; and the linux shell command was ~/Softare/RIsearch1–1/RIsearch -q lncRNA.fasta -t mRNA.fasta -d 30 -l 50 > RIsearch_run_out.txt.
Then, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of the target genes for lncRNAs . GO terms were enriched when P value was less than 0.05 using GOseq platform  and The KEGG pathways with a P < 0.05 were defined as significantly enriched pathways using KOBAS software .
Real-time quantitative PCR
We selected 5 lncRNAs and 5 mRNAs represent different expression levels for further qRT-PCR analysis. Total RNA from 5 adult roosters’ and 20 one-day-old rooster’ livers were extracted using Trizol reagent (TaKaRa, Dalian, China). The RNA was quantified using a NanoDrop® ND-1000 spectrophotometer (Thermo Scientific, MA, USA) with the OD value set at 260 nm; the purity was assessed by determining the OD260/OD280 ratio and the quality (RNA degradation and contamination) was further assessed using formaldehyde-agarose gel electrophoresis. RNA samples from the 5 adult chickens and 20 one-day-old chicks were analyzed by qPCR.
About 1 μg of total RNA was reverse transcribed using the PrimeScript™ RT reagent Kit with gDNA eraser (TaKaRa, Dalian, China). qPCR were performed using SYBR® Green PCR Master Mix (TaKaRa, Dalian, China).A 25 μL PCR mixture was quickly prepared from 12.5 μL of SYBR® Premix ExTaq II (2×), 1 μL of forward primer (10 μM/L), 1 μL of reverse primer (10 μM/L), 1 μL of cDNA, and 9.5 μL of double-distilled water. Primers for β-actin (internal control genes) as well as the differentially expressed mRNAs and lncRNAs (Additional file 11) were designed using Primer-BLAST (https://www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi?LINK_LOC=BlastHome). The PCR with amplifications was conducted in an iCycler iQ5 multicolor real-time PCR detection system (Bio-Rad Laboratories) and programmed as follows: 95 °C for 10 min; 40 cycles of 95 °C for 10 s; 60 °C for 30 s; 72 °C for 30 s; and 72 °C for 5 min. All samples were examined in triplicate.
All data were analyzed using the 2−ΔΔCt method . The statistical evaluation of experimental results was analyzed by Student’s T test using SPSS 20.0 statistical software. All data were expressed as means with standard error (SE). Differences were considered to be statistically significant at P < 0.05.
Si-Tayeb K, Lemaigre FP, Duncan SA. Organogenesis and development of the liver. Dev Cell. 2010;18:175–89.
Bhate A, Parker DJ, Bebee TW, Ahn J, Arif W, Rashan EH, et al. ESRP2 controls an adult splicing programme in hepatocytes to support postnatal liver maturation. Nat Commun. 2015;6:8768. https://doi.org/10.1038/ncomms9768.
Costa RH, Kalinichenko VV, Holterman AXL, Wang X. Transcription factors in liver development, differentiation, and regeneration. Hepatology. 2003;38:1331–47.
Cui JY, Renaud HJ, Klaassen CD. Ontogeny of novel cytochrome P450 gene isoforms during postnatal liver maturation in mice. Drug Metab Dispos. 2012;40:1226–37.
Peng L, Yoo B, Gunewardena SS, Lu H, Klaassen CD, Zhong XB. RNA sequencing reveals dynamic changes of mRNA abundance of cytochromes P450 and their alternative transcripts during mouse liver development. Drug Metab Dispos. 2012;40:1198–209.
Septer S, Edwards G, Gunewardena S, Wolfe A, Li H, Daniel J, et al. Yes-associated protein is involved in proliferation and differentiation during postnatal liver development. Am J Physiol Gastrointest Liver Physiol. 2012;302:G493–503.
Duanmu Z, Weckle A, Koukouritaki SB, Hines RN, Falany JL, Falany CN, et al. Developmental expression of aryl, estrogen, and hydroxysteroid sulfotransferases in pre-and postnatal human liver. J Pharmacol Exp Ther. 2006;316:1310–7.
Peng L, Paulson A, Li H, Piekos S, He X, Li L, et al. Developmental programming of long non-coding RNAs during postnatal liver maturation in mice. PLoS One. 2014;9:e114917. https://doi.org/10.1124/jpet.105.093633.
Ponting CP, Hardison RC. What fraction of the human genome is functional? Genome Res. 2011;21:1769–76.
Lindblad-Toh K, Garber M, Zuk O, Lin MF, Parker BJ, Washietl S, et al. A high-resolution map of human evolutionary constraint using 29 mammals. Nature. 2011;478:476–82.
Wang KC, Chang HY. Molecular mechanisms of long noncoding RNAs. Mol Cell. 2011;43:904–14.
He L, Hannon GJ. MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004;5:522–31.
Spurlock CF III, Tossberg JT, Guo Y, Collier SP, Crooke PS III, Aune TM. Expression and functions of long noncoding RNAs during human T helper cell differentiation. Nat Commun. 2015;6:6932. https://doi.org/10.1038/ncomms7932.
Shen Y, Mao H, Huang M, Chen L, Chen J, Cai Z, et al. Long noncoding RNA and mRNA expression profiles in the thyroid gland of two phenotypically extreme pig breeds using Ribo-zero RNA sequencing. Genes (Basel). 2016;7:34. https://doi.org/10.3390/genes7070034.
Li T, Wang S, Wu R, Zhou X, Zhu D, Zhang Y. Identification of long non-protein coding RNAs in chicken skeletal muscle using next generation sequencing. Genomics. 2012;99:292–8.
Bao X, Wu H, Zh X, Guo X, Hutchins AP, Luo Z, et al. The p53-induced lincRNA-p21 derails somatic cell reprogramming by sustaining H3K9me3 and CpG methylation at pluripotency gene promoters. Cell Res. 2015;25:80–92.
Mohammad F, Mondal T, Kanduri C. Epigenetics of imprinted long noncoding RNAs. Epigenetics. 2009;4:277–86.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PPA. ceRNA hypothesis: the Rosetta stone of a hidden RNA language? Cell. 2011;146:353–8.
Waterland RA, Kellermayer R, Rached MT, Tatevian N, Gomes MV, Zhang J, et al. Epigenomic profiling indicates a role for DNA methylation in early postnatal liver development. Hum Mol Genet. 2009;18:3026–38.
Sen S, Jumaa H, Webster NJ. Splicing factor SRSF3 is crucial for hepatocyte differentiation and metabolic function. Nat Commun. 2013;4:1336. https://doi.org/10.1038/ncomms2342.
Li T, Wan B, Huang J, Zhang X. Comparison of gene expression in hepatocellular carcinoma, liver development, and liver regeneration. Mol Gen Genomics. 2010;283:485–92.
Li T, Huang J, Jiang Y, Zeng Y, He F, Zhang MQ, et al. Multi-stage analysis of gene expression and transcription regulation in C57/B6 mouse liver development. Genomics. 2009;93:235–42.
Di C, Yuan J, Wu Y, Li J, Lin H, Hu L, et al. Characterization of stress-responsive lncRNAs in Arabidopsis thaliana by integrating expression, epigenetic and structural features. Plant J. 2014;80:848–61.
Adiconis X, Borges-Rivera D, Satija R, DeLuca DS, Busby MA, Berlin AM, et al. Comparative analysis of RNA sequencing methods for degraded or low-input samples. Nat Methods. 2013;10:623–9.
Zhao W, He X, Hoadley KA, Parker JS, Hayes DN, Perou CM. Comparison of RNA-Seq by poly (a) capture, ribosomal RNA depletion, and DNA microarray for expression profiling. BMC Genomics. 2014;15:419. https://doi.org/10.1186/1471-2164-15-419.
Bracken AP, Helin K. Polycomb group proteins: navigators of lineage pathways led astray in cancer. Nat Rev Cancer. 2009;9:773–84.
Sun X, Li M, Sun Y, Cai H, Li R, Wei X, et al. The developmental transcriptome landscape of bovine skeletal muscle defined by Ribo-zero ribonucleic acid sequencing. J Anim Sci. 2015;93:5648–58.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.
Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015;33:290–5.
Wang Y, Xue S, Liu X, Liu H, Hu T, Qiu X, et al. Analyses of long non-coding RNA and mRNA profiling using RNA sequencing during the pre-implantation phases in pig endometrium. Sci Rep. 2016;6:20238. https://doi.org/10.1038/srep20238.
Pauli A, Valen E, Lin MF, Garber M, Vastenhouw NL, Levin JZ, et al. Systematic identification of long noncoding RNAs expressed during zebrafish embryogenesis. Genome Res. 2012;22:577–91.
He Y, Ding Y, Zhan F, Zhang H, Han B, Hu G, et al. The conservation and signatures of lincRNAs in Marek’s disease of chicken. Sci Rep. 2016;6:19422. https://doi.org/10.1038/srep19422.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11:R14. https://doi.org/10.1186/gb-2010-11-2-r14.
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39:W316–22. https://doi.org/10.1093/nar/gkr483.
Konarzewski M, Gavin A, McDevitt R, Wallis IR. Metabolic and organ mass responses to selection for high growth rates in the domestic chicken (Gallus domesticus). Physiol Biochem Zool. 2000;73:237–48.
Kalmar ID, Vanrompay D, Janssens GP. Broiler ascites syndrome: collateral damage from efficient feed to meat conversion. Vet J. 2013;197:169–74. https://doi.org/10.1016/j.tvjl.2013.03.011.
Szunyogova E, Zhou H, Maxwell GK, Powis RA, Francesco M, Gillingwater TH, et al. Survival motor neuron (SMN) protein is required for normal mouse liver development. Sci Rep. 2016;6:35898. https://doi.org/10.1038/srep35898.
Chaudhari P, Tian L, Deshmukh A, Jang YY. Expression kinetics of hepatic progenitor markers in cellular models of human liver development recapitulating hepatocyte and biliary cell fate commitment. Exp Biol Med (Maywood). 2016;241:1653–62. https://doi.org/10.1177/1535370216657901.
Hsu SH, Delgado ER, Otero PA, Teng KY, Kutay H, Meehan KM, et al. MicroRNA-122 regulates polyploidization in the murine liver. Hepatology. 2016;64:599–615.
Hand NJ, Master ZR, Eauclaire SF, Weinblatt DE, Matthews RP, Friedman JR. The microRNA-30 family is required for vertebrate hepatobiliary development. Gastroenterology. 2009;136:1081–90.
Xu L, Yin W, Xia J, Peng M, Li S, Lin S, et al. An antiapoptotic role of sorting nexin 7 is required for liver development in zebrafish. Hepatology. 2012;55:1985–93.
Hämmerle M, Gutschner T, Uckelmann H, Ozgur S, Fiskin E, Gross M, et al. Posttranscriptional destabilization of the liver-specific long noncoding RNA HULC by the IGF2 mRNA-binding protein 1 (IGF2BP1). Hepatology. 2013;58:1703–12.
Qu X, Lam E, Doughman YQ, Chen Y, Chou YT, Lam M, et al. Cited2, a coactivator of HNF4α, is essential for liver development. EMBO J. 2007;26:4445–56.
Girard JR, Cuendet GS, Marliss EB, Kervran A, Rieutort M, Assan R. Fuels, hormones, and liver metabolism at term and during the early postnatal period in the rat. J Clin Invest. 1973;52:3190–200.
Marques PE, Oliveira AG, Chang L, Paula-Neto HA, Menezes GB. Understanding liver immunology using intravital microscopy. J Hepatol. 2015;63:733–42.
Rickett GM, Kelly FJ. Developmental expression of antioxidant enzymes in Guinea pig lung and liver. Development. 1990;108:331–6.
Ryu KY, Maehr R, Gilchrist CA, Long MA, Bouley DM, Mueller B, et al. The mouse polyubiquitin gene UbC is essential for fetal liver development, cell-cycle progression and stress tolerance. EMBO J. 2007;26:2693–706.
Miyaoka Y, Ebato K, Kato H, Arakawa S, Shimizu S, Miyajima A. Hypertrophy and unconventional cell division of hepatocytes underlie liver regeneration. Curr Biol. 2012;22:1166–75.
Fausto N, Campbell JS, Riehle KJ. Liver regeneration. Hepatology. 2006;43:S45–53. https://doi.org/10.1002/hep.20969.
Köhler UA, Kurinna S, Schwitter D, Marti A, Schäfer M, Hellerbrand C, et al. Activated Nrf2 impairs liver regeneration in mice by activation of genes involved in cell-cycle control and apoptosis. Hepatology. 2014;60:670–8. https://doi.org/10.1002/hep.26964.
Sadler KC, Krahn KN, Gaur NA, Ukomadu C. Liver growth in the embryo and during liver regeneration in zebrafish requires the cell cycle regulator, uhrf1. Proc Natl Acad Sci U S A. 2007;104:1570–5. https://doi.org/10.1073/pnas.0610774104.
Grote P, Wittler L, Hendrix D, Koch F, Währisch S, Beisaw A, et al. The tissue-specific lncRNA Fendrr is an essential regulator of heart and body wall development in the mouse. Dev Cell. 2013;24:206–14.
Fatica A, Bozzoni I. Long non-coding RNAs: new players in cell differentiation and development. Nat Rev Genet. 2014;15:7–21.
Zhang B, Arun G, Mao YS, Lazar Z, Hung G, Bhattacharjee G, et al. The lncRNA Malat1 is dispensable for mouse development but its transcription plays a cis-regulatory role in the adult. Cell Rep. 2012;2:111–23.
Li F, Xiao Y, Huang F, Deng W, Zhao H, Shi X, et al. Spatiotemporal-specific lncRNAs in the brain, colon, liver and lung of macaque during development. Mol BioSyst. 2015;11:3253–63.
Lv J, Huang Z, Liu H, Liu H, Cui W, Li B, et al. Identification and characterization of long intergenic non-coding RNAs related to mouse liver development. Mol Gen Genomics. 2014;289:1225–35.
Xu D, Yang F, Yuan JH, Zhang L, Bi HS, Zhou CC, et al. Long noncoding RNAs associated with liver regeneration 1 accelerates hepatocyte proliferation during liver regeneration by activating Wnt/β-catenin signaling. Hepatology. 2013;58:739–51.
Peng L, Paulson A, Li H, He X, Lu H, Klaassen CD, et al. Long noncoding RNAs and transcription of cytochrome P450s in mouse liver during maturation. FASEB J. 2013;27:1102–7.
Li H, Wang Y, Chen M, Xiao P, Hu C, Zeng Z, et al. Genome-wide long non-coding RNA screening, identification and characterization in a model microorganism Chlamydomonas reinhardtii. Sci Rep. 2016;6:34109. https://doi.org/10.1038/srep34109.
Wang TZ, Liu M, Zhao MG, Chen R, Zhang WH. Identification and characterization of long non-coding RNAs involved in osmotic and salt stress in using genome-wide high-throughput sequencing. BMC Plant Biol. 2015;15:131. https://doi.org/10.1186/s12870-015-0530-5.
Wenzel A, Akbaşli E, Gorodkin J. RIsearch: fast RNA–RNA interaction search using a simplified nearest-neighbor energy model. Bioinformatics. 2012;28:2738–46.
Zong Y, Panikkar A, Xu J, Antoniou A, Raynaud P, Lemaigre F, et al. Notch signaling controls liver development by regulating biliary differentiation. Development. 2009;136:1727–39.
Kodama Y, Hijikata M, Kageyama R, Shimotohno K, Chiba T. The role of notch signaling in the development of intrahepatic bile ducts. Gastroenterology. 2004;127:1775–86. https://doi.org/10.1093/bioinformatics/bts519.
Tanimizu N, Miyajima A. Notch signaling controls hepatoblast differentiation by altering the expression of liver-enriched transcription factors. J Cell Sci. 2004;117(Pt 15):3165–74.
Omenetti A, Choi S, Michelotti G, Diehl AM. Hedgehog signaling in the liver. J Hepatol. 2011;54:366–73.
Ochoa B, Syn WK, Delgado I, Karaca GF, Jung Y, Wang J, et al. Hedgehog signaling is critical for normal liver regeneration after partial hepatectomy in mice. Hepatology. 2010;51:1712–23.
Petrova R, Joyner AL. Roles for hedgehog signaling in adult organ homeostasis and repair. Development. 2014;141:3445–57.
Xie G, Karaca G, Swiderska-Syn M, Michelotti GA, Krüger L, Chen Y, et al. Cross-talk between notch and hedgehog regulates hepatic stellate cell fate in mice. Hepatology. 2013;58:1801–13.
McLin VA, Rankin SA, Zorn AM. Repression of Wnt/β-catenin signaling in the anterior endoderm is essential for liver and pancreas development. Development. 2007;134:2207–17.
Matsui T, Kinoshita T, Morikawa Y, Tohya K, Katsuki M, Ito Y, et al. K-Ras mediates cytokine-induced formation of E-cadherin-based adherens junctions during liver development. EMBO J. 2002;21:1021–30.
Birchmeier W. Orchestrating Wnt signalling for metabolic liver zonation. Nat Cell Biol. 2016;18:463–5.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–9.
Frazee AC, Pertea G, Jaffe AE, Langmead B, Salzberg SL, Leek JT. Ballgown bridges the gap between transcriptome assembly and expression analysis. Nat Biotechnol. 2015;33:243–6.
Sun L, Luo H, Bu D, Zhao G, Yu K, Zhang C, et al. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013;41:e166. https://doi.org/10.1093/nar/gkt646.
Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L, et al. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007;35:W345–9. https://doi.org/10.1093/nar/gkm391.
Zhang L, An XP, Liu XR, Fu MZ, Han P, Peng JY, et al. Characterization of the transcriptional complexity of the receptive and pre-receptive endometria of dairy goats. Sci Rep. 2015;5:14244. https://doi.org/10.1038/srep14244.
Ponjavic J, Oliver PL, Lunter G, Ponting CP. Genomic and transcriptional co-localization of protein-coding and long non-coding RNA pairs in the developing brain. PLoS Genet. 2009;5:e1000617.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25:402–8.
We thank Mr. Shen Lize, Miss. Lin Shuhong, and their colleagues at LC-BIO (Hangzhou, China) for assistance in data processing.
The work was supported by the National Key Research and Development Projects (2017YFD0502200, 20170500500), and the Program for Shaanxi Science & Technology (2017TSCXL-NY-04-04, 2015KTCQ02–19).
Availability of data and materials
The data set supporting the results of this article is available in the NCBI Gene Expression Omnibus repository, GSE101394 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE101394). All data generated or analyzed during this study are included in this published article and its Additional files.
Ethics approval and consent to participate
All experimental protocols and animals’ managements in the study were approved by the Institutional Animal Care and Use Committee (IACUC) of the Northwest A&F University (Yangling, Shaanxi, China).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Summary of protein-coding genes identified in the chicken liver libraries. (XLSX 4278 kb)
Table S2. Summary of lncRNAs identified in the chicken liver libraries. 1. Class code: “x” represent the antisense lncRNAs, “o” represent the sense lncRNAs, “j” represent the intronic lncRNAs, “i” represent the intron lncRNAs, and “u” represent the intergenic lncRNAs; 2. Known/Novel: We used reference genome download on Ensembl (galgal 5). However, there weas few annotated lncRNA in gga5 gtf files. As a result, all novel assembled transcripts (> 200 bp) by cufflinks with non protein-coding potential were labelled as novel lncRNAs. (XLSX 290 kb)
Table S3. Sequence information of all expressed lncRNAs found in the present study. (FA 69436 kb)
Table S4. Summary of differentially expressed protein-coding genes between mature livers and immature livers in chickens. (XLSX 498 kb)
Table S5. Summary of differentially expressed lncRNAs between mature livers and immature livers in chickens. (XLSX 32 kb)
Table S6. GO enrichment analysis (P < 0.05) of the differentially expressed mRNAs. S gene number: the number of significant differentially expressed mRNAs which match to a GO term; TS gene number: the number of significant differentially expressed mRNAs which have GO annotations; B gene number: the number of detected mRNAs which match to a GO term; TB gene number: the number of all detected mRNAs which have GO annotations. (XLSX 39 kb)
Table S7. The differentially expressed target protein-coding genes of the differentially expressed lncRNA in the present study. (XLSX 3850 kb)
Table S8. GO enrichment analysis (P < 0.05) of differentially expressed protein-coding genes targeted by differentially expressed lncRNAs in either trans- or cis- regulatory roles. S gene number: the number of significant differentially expressed lncRNAs whose cis- or trans-regulated mRNAs match to a GO term; TS gene number: the number of significant differentially expressed mRNAs lncRNAs whose cis- or trans-regulated mRNAs have GO annotations; B gene number: the number of detected lncRNAs whose cis- or trans-regulated mRNAs match to a GO term; TB gene number: the number of all detected lncRNAs whose cis- or trans-regulated mRNAs have GO annotations. (XLSX 19 kb)
Table S9. LncRNAs and its potential target genes that are involved in some key pathways related to the postnatal liver development. (XLSX 897 kb)
Table S10. Quality parameters of the RNA samples used in Ribo-zero RNA Sequencing. (XLSX 9 kb)
Table S11. Primers used in the qRT-PCR analysis. (XLSX 10 kb)
About this article
Cite this article
Wu, S., Liu, Y., Guo, W. et al. Identification and characterization of long noncoding RNAs and mRNAs expression profiles related to postnatal liver maturation of breeder roosters using Ribo-zero RNA sequencing. BMC Genomics 19, 498 (2018). https://doi.org/10.1186/s12864-018-4891-7
- Postnatal liver maturation
- Differential expression