- Open Access
Translational control plays a prominent role in the hepatocytic differentiation of HepaRG liver progenitor cells
Genome Biologyvolume 9, Article number: R19 (2008)
We investigated the molecular events associated with the differentiation of liver progenitor cells into functional and polarized hepatocytes, using human HepaRG cells that display potent hepatocytic differentiation-inducible properties and share some features with liver progenitor cells.
Profiling of total and of polysome-bound transcripts isolated from HepaRG cells undergoing hepatocytic differentiation was performed. A group of 3,071 probe sets was reproducibly regulated by at least 2-fold in total or in polysome-bound RNA populations, upon differentiation. The fold changes in the total and the polysome-bound RNA populations for these 3,071 probe sets were poorly correlated (R = 0.38). Moreover, while the majority of the regulated polysome-bound RNA probe sets were up-regulated upon differentiation, the majority of the regulated probe sets selected from the total RNA population was down-regulated. Genes translationally up-regulated were associated with cell cycle inhibition, increased susceptibility to apoptosis and innate immunity. In contrast, genes transcriptionally up-regulated during differentiation corresponded in the majority to liver-enriched transcripts involved in lipid homeostasis and drug metabolism. Finally, several epithelial and hepato-specific transcripts were strongly induced in the total RNA population but were translationally repressed.
Translational regulation is the main genomic event associated with hepatocytic differentiation of liver progenitor cells in vitro and targets genes critical for moderating hepatocellular growth, cell death and susceptibility to pathogens. Transcriptional regulation targets specifically liver-enriched transcripts vital for establishing normal hepatic energy homeostasis, cell morphology and polarization. The hepatocytic differentiation is also accompanied by a reduction of the transcript content complexity.
Liver diseases represent a major public health burden worldwide . Upon acute liver injury, the mature hepatocytes demonstrate a major proliferative capacity. However, in chronic liver diseases such as chronic hepatitis B virus and hepatitis C virus infections and alcohol abuse, their regenerative potential is often impaired and liver progenitor cells, also called oval cells, significantly increase both in number and their capability to proliferate [2, 3]. In recent years, liver progenitor cells have drawn special interest not only because of their regenerative capability and, therefore, therapeutic potential but also because of their possible contribution to liver carcinogenesis [4–6]. Rodent and simian liver progenitor cell lines have been established [7–10] and shown to successfully repopulate diseased livers [11–13].
The HepaRG cell line is a naturally immortalized human liver cell line with progenitor properties and bipotent differentiation-inducible capability that has been established from the non-tumoral region of a resected hepatitis C virus-associated hepatocellular carcinoma (HCC) [14, 15]. These bipotent progenitor cells have been found to repopulate uPA/SCID mouse damaged livers . Throughout differentiation, HepaRG cells evolve from a homogeneous dedifferentiated, depolarized, epithelial phenotype showing no specific organization to a structurally well-defined and polarized monolayer closely resembling those formed in primary human hepatocytes in culture, with canaliculi-like structures . At the hepatocytic differentiated state, hepatocytic polarization markers such as ZO-1 and CD26 and liver-specific proteins such as albumin are expressed at levels similar to those found in normal liver biopsies [14, 15]. Finally, iron storage and metabolism, typical features of mature normal hepatocytes, are intact in HepaRG cells . Although this system bears limitations inherent to its pathological origin, it represents to date the only in vitro human model for hepatocytic differentiation.
We used this powerful system to identify the genomic events associated with the development of a functional and polarized hepatocyte-like cell from a previously dedifferentiated epithelial progenitor. A role for translational control in liver development and for translation regulators such as p70S6 kinase and 4E-BP1 upon liver regeneration has been previously reported [18–21]. Therefore, integrating polysome-bound RNA profiling to total RNA profiling not only provides highly relevant phenotypic information, but also provides insight into the role of translational control on the specific biological process studied.
Results and discussion
Total and polysome-bound RNA changes associated with hepatocytic differentiation of HepaRG cells
HepaRG cells were induced to differentiate into morphologically and functionally mature hepatocyte-like cells. Differentiated HepaRG cells showed features of normal hepatocytes, such as refractile cellular borders, clearly delineated nuclei and tridimensional polarization with the appearance of refringent circular canaliculi vertically (Figure 1). In order to identify the genomic events associated with HepaRG cell differentiation, total RNA and polysome-bound RNA were isolated at the proliferative stage and at the end of the differentiation protocol and analyzed on Affymetrix Human Genome U133A arrays (Figure 1). We separated polysomes from free messenger ribonucleoproteins (mRNPs) using sucrose gradient centrifugation with the assumption that translationally inactive mRNAs are present as free cytoplasmic mRNPs, whereas actively translated mRNAs are contained within polysomes. Total RNA was processed in parallel for each sample.
Out of the 22,283 probe sets spotted on the array, 3,071 (13.8%) were modulated by at least 2-fold upon differentiation and in 3 independent experiments, either in the total RNA or the polysome-bound RNA compartments. Total RNA fold changes were plotted against polysome-bound RNA fold changes for these 3,071 probe sets (Figure 2a). The correlation coefficient for the regression curve calculated from all values was 0.38, demonstrating a poor correlation and, therefore, an uncoupling phenomenon between changes in the polysome-bound fractions and changes in total RNA upon differentiation of HepaRG cells. We then determined the distribution of up- and down-regulated transcripts in each RNA population upon differentiation. In the total RNA compartment, 547 and 1,636 probe sets (a total of 2,183) were up-regulated and down-regulated, respectively. In contrast, in the polysome-bound RNA compartment, 1,325 and 124 probe sets (a total of 1,449) were up-regulated and down-regulated, respectively (Figure 2b). Transcription is, therefore, largely down-regulated during HepaRG differentiation while translation of specific genes is up-regulated. Probe sets that are similarly up-regulated or down-regulated in both RNA populations correspond to genes modulated as a result of transcriptional regulation without any subsequent translational control. These probe sets represented only a small number of genes with 359 up-regulated and 88 down-regulated probe sets. They represented 14.6% of the initially selected 3,071 regulated probe sets (Figure 2b, dark portions of the graph bars). On the other hand, 2,624 probe sets (85.4% of the total number of regulated probe sets) were modulated due to translational control (Figure 2b, gray portions of the bar graphs).
A subset of genes was selected for validation. Validation was performed using real-time PCR on the total RNA and the polysome-bound RNA populations, for ten genes: those encoding apolipoprotein H, solute carrier (SLC)27A3, cytochrome P450 isoforms 3A4 and 7B1, vascular endothelial growth factor (VEGF), E-cadherin, insulin receptor, leptin receptor, transforming growth factor (TGF) beta receptor 2 and membrane metallo-endopeptidase (MME). The PCR results obtained on the three independent experiments confirmed the microarray data for all ten genes (Figure 3a). Validation was also performed using real time PCR on each fraction of the sucrose gradient separating free mRNPs and polysomes, for three genes: those encoding latent transforming growth factor beta binding protein 1 (LTBP1), spectrin repeat-containing nuclear envelope 1 (SYNE-1) and matrix metalloproteinase 3 (MMP3). A shift was observed upon HepaRG differentiation for all three transcripts from the free mRNP fractions to the heavier polysome fractions on the sucrose gradient as shown in Figure 3b for LTBP1. These results demonstrate an increased translation of these transcripts and validate the array data indicating no change or a slight decrease in LTBP1, SYNE-1 and MMP3 transcript levels in the total RNA compartment and a strong up-regulation of all three transcripts in the polysome-bound RNA compartment.
All together, these results suggest that translational control plays a prominent role in the hepatocytic differentiation of liver progenitor cells and that the total RNA content may not be representative of the mature phenotype of hepatocyte-like cells. In addition, transcriptional changes did not overlap with translational changes. The large majority of polysome-bound (that is, translated) genes modified were up-regulated whereas the majority of genes modified at the total RNA level were down-regulated, suggesting that the mature hepatocyte phenotype is acquired by increased translation of pre-existing transcripts. The total RNA population can be considered as a stock of translated and untranslated transcripts that can be utilized by the cell rapidly. The more diverse the total RNA population is, the greater the options the cell has in selecting protein expression patterns. Therefore, the extensive down-regulation of genes in the total RNA compartment can be interpreted as a decrease in cellular RNA diversity, consistent with the commitment of a dedifferentiated epithelial progenitor into a defined, in this case hepatocytic, lineage.
Polysome-bound RNA changes associated with HepaRG cell differentiation: the hepatocytic phenotype
To further characterize the differentiated phenotype of HepaRG cells, we selected all polysome-bound up-regulated probe sets (n= 1641) and all polysome-bound down-regulated probe sets (n= 204), regardless of their fold-change status at the total RNA level. The content of these two lists of genes were separately analyzed using the Ingenuity Systems Pathways Knowledge Base . This database enables one to search for gene products' interactions and annotations coming from curated data from publications and peer-reviewed resources. Networks displaying significant overlap between the selected regulated genes found in our study and the software-preselected members were selected. The Ingenuity pathway analysis identified nine networks (networks A-I) and one network (network J) generated from the up-regulated and down-regulated transcripts, respectively (Table 1 and Additional data file 1). These ten networks can be divided into six groups based on their associated biological top functions: cell cycle, cell death, innate immunity, lipid and drug metabolism, cell morphology, and cell environment and movement.
Network A (Additional data file 1, A) was organized around transcription factors with tumor suppressor activities. These included three members of the SMARC tumor suppressor family (SMARCA2, SMARCB1 and SMARCC2), the transcription factors MEF2C and MEF2D and the NF-KB inhibitor NF-KB1A. Interestingly, several of these transcription factors (SMARC, MEF) remain uncharacterized in the liver.
Network B (Additional data file 1, B) was associated with increased susceptibility to apoptosis and included the initiator caspase 8, insulin growth factor-binding protein (IGFBP)1, inhibitor of hepatocytic proliferation in vivo and in vitro , the interferon-induced gene IFI16, an essential mediator of p53 function  and tuberous sclerosis complex protein 2 (TSC2). The presence of Kininogen 1, a component of the coagulation cascade produced by the mature hepatocyte, confirmed the differentiation status of the cells. Cell death was also a top function of network C (Additional data file 1, C) with the presence of another member of the initiator caspase family, caspase 9, and of FOXO3A, known to trigger caspase 9-induced apoptosis. Other members associated with cell death included two strong inducers of apoptosis in human hepatocytes, TNFSF10/TRAIL  and IRF3  and two members of the BCL2 family, BCL2 and BCL2L11. While BCL2 protects cells against apoptosis, BCL2L11 facilitates this process of cell death by neutralizing BCL2 antiapoptotic activity . Therefore, the concomitant upregulation of BCL2 and BCL2L11, together with the pro-apoptotic genes described above, suggest that upon their differentiation, liver progenitor cells become highly susceptible to apoptosis. It has been reported that normal hepatocytes are highly sensitive to cell death upon, for example, drug-induced liver toxicity and that three-dimensional polarization, as occurs in this system (Figure 1), sensitizes hepatocytes to Fas apoptotic signaling . Noteworthy, both up-regulated caspases identified (caspases 8 and 9) belong to the initiator caspase family, while none of the members of the effector caspase family (caspases 3, 6 and 7)  was affected, supporting the observation that the cells did not undergo apoptosis in culture.
Another function associated with network C (Additional data file 1, C) was innate immunity and responses to viral infections, with the presence of two members of the interferon-regulatory factors, IRF3 and IRF7. IRF3 is a key component of innate immunity in the hepatocyte and has been shown to mediate interferon (IFN)β induction upon hepatitis C virus infection . IRF7 is also mandatory for a proper IFNα-dependent antiviral response against hepatitis C virus . Their up-regulation upon differentiation suggests an association between hepatocytic differentiation and innate immunity maturation. Maturation of the innate immunity upon differentiation was also suggested in network D (Additional data file 1, D) with the up-regulation of STAT1, one of the major components of the type I IFN transduction pathway, playing a key role in antiviral defense, inflammation and injury  and the up-regulation of complement C3 with a role in innate immunity as well as in acute phase response . This network also included the EGFR-like receptor ERRB3 associated with cell survival and CDK5 reported to inhibit FAS/STAT3-dependent apoptosis in hepatoma cell lines in vitro and in vivo .
Lipid metabolism and drug metabolism
Network E (Additional data file 1, E) included the peroxisome proliferative activated receptor alpha (PPARA), regulating the expression of several hepatic genes and lipid homeostasis in the liver , as well as CEBPD and STAT3, key players in the control of the acute-phase response as well as in the protection of the hepatocyte upon acute phase-related injury [32, 33, 36]. As expected, apolipoproteins A1 and C3 as well as fibrinogens A, B, and G, markers of functional differentiation of the hepatocyte in relation to lipid metabolism and acute phase response, were strongly upregulated, downstream of PPARA, CEBPD and STAT3. Network F (Additional data file 1, F) included the liver-enriched transcription factors CAAT/enhancer-binding protein alpha (CEPBA), retinoid X receptor alpha (RXRA), and the peroxisome proliferative activated receptor delta (PPARD). CEBPA regulates two aspects of hepatic terminal differentiation: induction of differentiation-specific genes and repression of mitogenesis [37–39]. RXRA regulates cholesterol, fatty acid, bile acid, steroid, and xenobiotic metabolism and homeostasis in the liver. PPARD also plays a role in lipid metabolism, including cholesterol efflux and fatty acid oxidation [40, 41], activates fat metabolism to prevent obesity , and regulates fatty acid synthesis, glucose metabolism and insulin sensitivity . Network G (Additional data file 1, G) included the sterol regulatory element-binding transcription factor-1 (SREBF1), a major regulator of sterol biosynthesis, hepatic gluconeogenesis and lipogenesis in the liver , the liver-enriched transcription factor retinoid X receptor beta (RXRB) , MLXIPL, a glucose-responsive transcription factor that regulates carbohydrate metabolism in the liver , and angiotensinogen, an endocrine product of the hepatocyte regulating blood pressure . ADH1A, ADH1B and ADH1C, mature hepatocyte-specific inducible genes involved in ethanol metabolism , were also included in this network.
Network H (Additional data file 1, H) contained CDC42, a small GTPase involved in cell polarity. STAT6, also included in this network, is involved in the induction of a TH1 immune response to the hepatocyte and protects the normal parenchyma against liver injury . Jak2 participates in transduction of interleukin (IL)6 signaling in case of acute phase reaction, as well as in the signal transduction of IFNγ . The COP proteins (COPE, COPG, COPZ1, COPA, COPB2) mediate transport between the Golgi and the endoplasmic reticulum . Their up-regulation may be associated with the increased flux of secreted proteins en route to the extracellular compartment through the Golgi complex after synthesis in the mature hepatocyte.
Cell environment and movement
Network I (Additional data file 1, I) included fibronectin (FN1), a co-factor of endogenous anti-angiogenic molecules and enhancer of cell attachment , and EGR1. EGR1 controls FIN1 and TGFβ1 gene expression and acts as a cell cycle blocker in vitro and in vivo through p53 . This network also included MMP3, a secreted metalloprotease implicated in metastasis [52, 53], IGFBP2, an insulin growth factor-binding protein associated with hepatocytic proliferation inhibition in vivo and in vitro  and two members of the serine protease inhibitors, SERPINF2 and SERPINA5. Network J (Additional data file 1, J), the only network associated with down-regulated polysome-bound probe sets, was also associated with cellular movement. Notably, the components of this network included several growth factors and secreted proteins implicated in angiogenesis and metastasis, such as hepatocyte growth factor (HGF), VEGF, platelet-derived growth factor (PDGF)-B, CCL2 and IL8. VEGF and PDGF-B are potent mitogenic and angiogenic factors . HGF is the primary agent promoting the proliferation and apoptosis resistance of mature hepatocytes . CCL2 is a monocyte chemoattractant . IL8 is a proinflammatory cytokine and chemoattractant for neutrophils . Therefore, differentiation of hepatocytic progenitors seems to be associated with a progressive disappearance of an inflammation-like state, as shown by the down-regulation of several chemoattractants and proinflammatory messengers.
Taken together, this analysis identified the regulation of functions specific to a differentiated hepatocytic phenotype. Up-regulation of transcripts belonging to the well known liver-enriched transcription factors, such as CEBPA, RXRA, RXRB, and PPARD, as well as down-regulation of NF-KB expression, are correlated with the differentiation of liver progenitor cells into morphologically and functionally mature hepatocyte-like cells. This study also revealed the involvement of lesser known nuclear proteins in the hepatocytic biology, such as SMARC, MEF and EGR1 proteins, and novel associations, such as the role of several IFN-associated or induced proteins in the acquisition of the hepatocytic phenotype. STAT1 is one of the key elements for the induction of the type I IFN response. Its up-regulation, as well as the up-regulation of several other IFN-related transcripts (OAS1, IRF3, IRF7 and IFI16), suggest that acquisition of key elements to innate immunity is associated with hepatocytic differentiation. It would be interesting, therefore, to investigate if the progenitor cell compartment in regenerative livers of chronically hepatitis B or C virus-infected patients is more prone to viral replication because of an immature innate immunity status.
Contribution from translation
Most of the genes identified in this study and contributing to the differentiation phenotype were modulated by translational control. Translationally regulated transcripts are underlined in Table 1 and indicated in blue in Additional data file 1. To investigate whether translational control specifically affects transcripts involved in defined cellular functions, we calculated the percentage of translationally controlled probe sets in each of the ten networks A-J described above. Paired t-tests were performed between groups of networks sharing the same cellular functions (Figure 4). A significantly greater involvement of translational control was observed in networks related to cell cycle and cell death functions than in networks related to lipid metabolism and drug metabolism (p = 0.005). Likewise, a significantly stronger involvement of translational control was found in innate immunity-related networks compared to cell environment and cell movement-related networks (p = 0.027). The high percentage of translationally controlled probe sets in cell cycle and cell death-related networks is in agreement with the ability of the hepatocyte to massively and rapidly proliferate under acute liver injury, as well as with the hypersensitivity of the hepatocyte to cell death in response, for example, to drug-associated toxicity. Translationally regulated transcripts associated with cell cycle included the nuclear proteins SMARCA2 and SMARCB1, the transcription factors MEF2C, MEF2D and EGR1 and the NF-KB inhibitor NFKBIA. Translationally regulated transcripts associated with cell death included oncostatin M receptor/IL6ST and the initiator caspases 8 and 9. Translationally regulated transcripts associated with innate immunity included several interferon-associated genes, such as those encoding OAS1, IRF3 and IFI16. Finally, numerous transcription factors associated with inflammation were translationally upregulated and included the three liver-enriched transcription factors RARA, RXRA and RXRB and STAT6 (Table 2).
Numerous transcription factors were translationally upregulated while left unchanged or even decreased at the total RNA level. Translational control of these transcription factors provides the cell with a means to modify its phenotype in a timely manner, rapidly expressing genes downstream of these transcription factors. The hepatocyte has to be a highly versatile cell because of at least two of its functions: the ability to generate the acute phase reaction and to maintain blood homeostasy after meals as the first line organ downstream of the portal vein that carries nutrients from the digestive tract.
The importance of translational control during liver progenitor cell differentiation raises the question of the identity of the actors involved. We recently reported a functional down-regulation of the mTOR/4E-BP1/p70S6 kinase pathway during differentiation of HepaRG cells . Moreover, forced expression of an activated mutant of mTOR impairs hepatocytic differentiation in this model . This pathway may therefore contribute at least partially to some of the translational events described here.
Contribution from transcription
Some genes were similarly modified upon differentiation of HepaRG cells, in both the total and the polysome-bound RNA populations, indicative of a transcriptional regulation. These include 435 up-regulated and 142 down-regulated probe sets (Figure 2b), indicated in yellow in Additional data file 1 and not underlined in Table 1. These genes corresponded in the majority to liver-enriched transcripts and to genes involved in lipid and drug metabolism. They included those encoding PPARA, PPARD, CEBPA, the hepatic leukemia factor (HLF) and the alcohol dehydrogenases 1B, 1C and 6. Other transcriptionally regulated genes included those encoding plasma proteins synthesized in the liver: the SERPINs A1, A4, F2, several complement system subunits (C1S, C3, C4A, C5 and C6) and three forms of fibrinogen (A, B and G). Finally, several cytokines, chemokines or hormones and their receptors were transcriptionally regulated as well: TNFSF10/TRAIL, IL6R, BMP2 and PDGFB (Table 2).
As the contribution of transcription appeared restricted to selective genes during HepaRG cell differentiation, we sought to investigate the expression levels and phosphorylation status of the canonic hepatocytic transcription factors HNF1α and HNF4α throughout differentiation. HNF1α is a major player in the acquisition of central hepatocytic functions, including gluconeogenesis, carbohydrate synthesis and storage, lipid metabolism (synthesis of cholesterol and apolipoproteins), detoxification (synthesis of cytochrome P450 monooxygenases), and synthesis of serum proteins (albumin, complements, and coagulation factors) . Interestingly, neither total nor polysome-bound RNA levels of HNF1α were modulated (-1.38 and +1.48-fold, respectively). This observation was confirmed by real time PCR (+1.38 ± 0.08 fold (mean ± standard error of the mean (SEM)) in total RNA and +1.02 ± 0.19 fold (mean ± SEM) in polysome-bound RNA; Figure 5a). In addition, no changes were observed at the protein expression level nor in phosphorylation status for HNF1α (55% of HNF1α is phosphorylated at the proliferative stage versus 38% at the differentiated stage; Figure 5b). HNF4α was slightly increased in both total and polysome-bound RNA (+1.89-fold and +1.35-fold, respectively). These slight increases were confirmed by real time PCR (+2.71 ± 0.13 fold (mean ± SEM) in total RNA and +1.74 ± 0.06 fold (mean ± SEM) in polysome-bound RNA; Figure 5c). However, HNF4α phosphorylation was strongly induced upon differentiation (Figure 5d), suggesting that, in contrast to HNF1α, HNF4α may contribute to HepaRG cell differentiation. Mutations of HNF1α associated with metabolic diseases have been described [60, 61] and, therefore, we cannot exclude that the lack of regulation of HNF1α found in this study results from mutation(s) disrupting its biochemical characteristics. However, the patient that gave rise to HepaRG cells was not known to be affected by any of these diseases.
In conclusion, transcriptional control appears to play a highly selective role in the phenotype of liver progenitor cell maturation and specifically targets liver-enriched transcripts characteristic of the mature hepatocytic phenotype. Novel findings suggest that the complement system is induced during maturation following transcriptional regulation.
Several transcripts were strongly transcriptionally induced upon HepaRG cell differentiation while unchanged or induced to a much weaker level in the polysome-bound RNA population, suggesting a translational repression control. Examples include E-cadherin, involved in hepatocytic polarization, cytochrome P450 3A4, a steroid-inducible cytochrome P450 isoform, cytochrome P450 7B1, a cytochrome P450 isoform involved in cholesterol metabolism, cytochrome P450 2A6 and 2C19, cytochrome P450 isoforms involved in drug metabolism, TGF-β receptor 2 and VEGF, an important regulator of angiogenesis and metastasis (Table 2).
Interestingly, four isoforms of cytochrome P450 were strongly up-regulated at the total RNA level but not at the polysome-bound RNA level. Given that cytochromes are inducible proteins involved in drug and lipid metabolism, high levels of untranslated RNA could serve as a stock that may be rapidly translated and used for the detoxification and acute phase-associated functions of the hepatocyte.
The most prominent result of this study is a strong association between translational control and hepatocytic differentiation of liver progenitor cells, as demonstrated by the fact that the great majority of the regulated genes have been identified in the polysome-bound RNA population and not in the total RNA population. Another interesting feature supporting the involvement of translational control in hepatocytic differentiation of liver progenitor cells is that the large majority of polysome-bound transcripts modified upon differentiation were up-regulated whereas the majority of genes modified in the total RNA population were down-regulated. Altogether, these data suggest that the mature hepatocyte phenotype is acquired by increased translation of pre-existing transcripts and is associated with a reduction in the diversity of transcripts that the differentiated cell can utilize, consistent with the commitment of a dedifferentiated epithelial progenitor into a defined hepatocytic lineage. This study increases our knowledge on gene expression regulation of liver progenitor cells upon differentiation, providing novel paths to successfully use liver progenitor cells to repopulate diseased livers.
Materials and methods
The HepaRG cell line was cultured in William's E medium (Invitrogen, Carlsbad, CA, USA) supplemented with 10% fetal calf serum (Mediatech, Manassas, VA, USA), 100 units/ml penicillin, 100 μg/ml streptomycin (Invitrogen), 5 μg/ml insulin (Sigma-Aldrich, St. Louis, MO, USA), and 5 × 10-5 M hydrocortisone hemisuccinate (Sigma-Aldrich). To induce differentiation, a two-step procedure was used as previously described [14, 15]. Cells were seeded at a density of 4 × 104 cells/cm2 and maintained for 2 weeks in the growth medium. Then, the culture medium was supplemented with 1% DMSO (Sigma-Aldrich) and 20 ng/ml EGF (PeproTech, Rocky Hill, NJ, USA) for 2 additional weeks. Cells were harvested either at 2 days (proliferative stage) or at 28 days (differentiation stage) after seeding. Cell culture pictures were taken using a phase contrast microscope (Nikon). Differentiation was evaluated morphologically by counting bile canaliculi (refringent area) at the intersection of two or three hepatocyte-like cells.
Total RNA extraction and polysome fractionation
Total RNA was extracted, precipitated and resuspended in RNAse-free water using Trizol reagent (Invitrogen) according to the manufacturer's instructions. For polysome fractionation, cycloheximide (100 μg/ml) was added to the medium for 3 minutes prior to harvest. The medium was then removed and the cells were washed with ice-cold phosphate-buffered saline containing 100 μg/ml cycloheximide. The cells were then scraped, centrifuged at 800g for 5 minutes at 4°C and cytoplasmic RNA was obtained by lysis of the cell pellet in 1 ml of polysome buffer containing 10 mM Tris-HCl (pH 8.0), 140 mM NaCl, 1.5 mM MgCl2, 0.5% Nonidet P-40, and a ribonuclease inhibitor, RNasin (500 units/ml; Promega, Madison, WI, USA). After the removal of nuclei, the cytosolic supernatant was supplemented with 100 μg/ml cycloheximide, 665 μg/ml heparin, 20 mM dithiothreitol, and 1 mM phenylmethanesulfonyl fluoride. Mitochondria and membrane debris were removed by centrifugation, and the post-mitochondrial supernatant was overlaid onto a 15-40% sucrose gradient and spun at 38,000 rpm for 2 h at 4°C in a SW41Ti rotor (Beckman Coulter, Fullerton, CA, USA). Fractions (750 μl) were collected from the top of each gradient and deproteinated with 100 μg of proteinase K in the presence of 1% SDS and 10 mM EDTA. After acid phenol extraction, RNA integrity was controlled by electrophoresis analysis on 1.2% agarose gel. Densitometry (GelDoc, Bio-Rad Laboratories, Hercules, CA, USA) was used to identify the fractions in which the 28S/18S ratio equals 2 (that is, fractions corresponding to polysome-bound RNA). These fractions were pooled from each sucrose gradient.
Microarray hybridization and data mining
Total and polysome-bound RNAs were purified using the RNeasy mini-kit clean-up protocol (Qiagen, Valencia, CA, USA), RW1 buffer being used to efficiently remove heparin from the samples. The first-strand cDNA, the double-strand cDNA, and cRNA were synthesized, and cRNA was fragmented using Affymetrix kits and guidelines . All cRNA final products were tested in terms of amount and integrity by Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) prior to microarray hybridization. cRNA samples were processed on Affymetrix HGU133A arrays with strict adherence to the labeling, hybridization and staining protocols provided by Affymetrix. A 'present' (P), 'marginal' (M) or 'absent' (A) call was assigned to each probe set using Affymetrix GeneChip Operation Software (GCOS v1.4). Probe sets with an 'absent' (A) call on all arrays were filtered out. Background correction and normalization steps were carried out using the GC-RMA method and the R Bioconductor software . Microarray data have been deposited in the ArrayExpress repository  under the accession number E-MEXP-1082. Three independent differentiation processes were performed and the correlation coefficients between each duplicate at the proliferative and the differentiated level were calculated using scatter plots. For total RNA, the correlation coefficients were 0.98 and 0.95 for proliferative and differentiated cells, respectively. For polysome-bound RNA, the values were 0.98 and 0.97 for proliferative and differentiated cells, respectively. The Ingenuity pathway analysis  was used to analyze selected probe sets. Each gene identifier was mapped to its corresponding gene object in the Ingenuity Pathways Knowledge Base. The application utilizes a right-sided Fisher's exact test to identify networks that had higher odds ratio of containing significant genes. These genes, called Focus Genes, were then overlaid onto a global molecular network. Networks of these Focus Genes were then algorithmically generated.
Cells were lysed in 50 mM Tris-HCl (pH 8), 150 mM NaCl, 0.1% SDS, 1% NP-40 supplemented with protease inhibitors (Complete, Roche Diagnostics, Indianapolis, IN, USA). Thirty micrograms of proteins were resolved on 10% SDS-polyacrylamide gels and electrotransferred onto nitrocellulose membrane (Amersham Biosciences, Piscataway, NJ, USA). Equal loadings and homogeneous blotting were confirmed using Ponceau red staining. Membranes were blocked with 5% nonfat milk in Tris-buffered saline and incubated with primary antibodies (anti-HNF1α and anti-HNF4α, Santa Cruz Biotechnology, dilution 1/500, Santa Cruz, CA, USA) overnight. Horseradish peroxidase-conjugated immunoglobulins (Dako, dilution 1/1,000, Carpinteria, CA, USA) were used as secondary antibodies and proteins were visualized with enhanced ECL chemiluminescent reagent (Amersham Biosciences). Densitometry was performed using the Total Lab TL100 software.
One microgram of DNAse I-treated (Promega) total RNA or polysome-bound RNA was reverse transcribed using Moloney murine leukemia virus reverse transcriptase and random hexamers (Invitrogen) for 50 minutes at 42°C. cDNA mixtures (1/10) were mixed with an equal volume of 2 × iQ SYBR green supermix (Bio-Rad Laboratories). Amplification was then performed at an annealing temperature of 55°C or 60°C and an elongation time of 30 s or 1 minute on a MyIQ real-time PCR apparatus (Bio-Rad Laboratories). Primer sequences were obtained through the Primer Bank website [65, 66] and are described in Additional data file 2. Differential expression ratios between proliferative and differentiated stages were calculated using the Δ(ΔCt) formula. Specificity of all amplicons was assessed by post-run melting curve analysis and agarose gel electrophoresis. For the analysis on sucrose gradient fraction distribution, 100 μl of each fraction harvested from the sucrose gradients were purified using the RNAeasy kit (Qiagen). Ten microliters of collected RNA were DNAse-digested and reverse transcribed as described above. Since the RNA amount in each fraction was different, in order to avoid efficiency differences of the reverse transcriptase, RNA amounts were equalized by adding an appropriate amount of in vitro transcribed irrelevant RNA to each fraction, giving a final amount of 1 μg of RNA on each sample. One tenth of each cDNA reaction was processed by real-time PCR as described above.
Additional data files
CAAT/enhancer-binding protein alpha
hepatocyte growth factor
hepatic leukemia factor
insulin growth factor-binding protein
latent transforming growth factor beta binding protein 1
matrix metalloproteinase 3
platelet-derived growth factor
peroxisome proliferative activated receptor alpha
peroxisome proliferative activated receptor delta
retinoid X receptor alpha
retinoid X receptor beta
standard error of the mean
serine protease inhibitor
sterol regulatory element-binding transcription factor-1
spectrin repeat-containing nuclear envelope 1
transforming growth factor
tuberous sclerosis complex protein 2
vascular endothelial growth factor.
Williams R: Global challenges in liver disease. Hepatology. 2006, 44: 521-526. 10.1002/hep.21347.
Overturf K, al-Dhalimy M, Ou CN, Finegold M, Grompe M: Serial transplantation reveals the stem-cell-like regenerative potential of adult mouse hepatocytes. Am J Pathol. 1997, 151: 1273-1280.
Sell S, Ilic Z: Liver stem cells. Stem cells. Edited by: Potten CS. 1997, San Diego, CA: Academic Press, 233-282.
Matthews VB, Yeoh GC: Liver stem cells. IUBMB Life. 2005, 57: 549-553.
Roskams TA, Libbrecht L, Desmet VJ: Progenitor cells in diseased human liver. Semin Liver Dis. 2003, 23: 385-396. 10.1055/s-2004-815564.
Lee JS, Heo J, Libbrecht L, Chu IS, Kaposi-Novak P, Calvisi DF, Mikaelyan A, Roberts LR, Demetris AJ, Sun Z, Nevens F, Roskams T, Thorgeirsson SS: A novel prognostic subtype of human hepatocellular carcinoma derived from hepatic progenitor cells. Nat Med. 2006, 12: 410-416. 10.1038/nm1377.
Yin L, Sun M, Ilic Z, Leffert HL, Sell S: Derivation, characterization, and phenotypic variation of hepatic progenitor cell lines isolated from adult rats. Hepatology. 2002, 35: 315-324. 10.1053/jhep.2002.31355.
Strick-Marchand H, Weiss MC: Inducible differentiation and morphogenesis of bipotential liver cell lines from wild-type mouse embryos. Hepatology. 2002, 36: 794-804.
Schwartz RE, Reyes M, Koodie L, Jiang Y, Blackstad M, Lund T, Lenvik T, Johnson S, Hu WS, Verfaillie CM: Multipotent adult progenitor cells from bone marrow differentiate into functional hepatocyte-like cells. J Clin Invest. 2002, 109: 1291-1302. 10.1172/JCI200215182.
Allain JE, Dagher I, Mahieu-Caputo D, Loux N, Andreoletti M, Westerman K, Briand P, Franco D, Leboulch P, Weber A: Immortalization of a primate bipotent epithelial liver stem cell. Proc Natl Acad Sci USA. 2002, 99: 3639-3644. 10.1073/pnas.062038599.
Strick-Marchand H, Morosan S, Charneau P, Kremsdorf D, Weiss MC: Bipotential mouse embryonic liver stem cell lines contribute to liver regeneration and differentiate as bile ducts and hepatocytes. Proc Natl Acad Sci USA. 2004, 101: 8360-8365. 10.1073/pnas.0401092101.
Fougere-Deschatrette C, Imaizumi-Scherrer T, Strick-Marchand H, Morosan S, Charneau P, Kremsdorf D, Faust DM, Weiss MC: Plasticity of hepatic cell differentiation: bipotential adult mouse liver clonal cell lines competent to differentiate in vitro and in vivo. Stem Cells. 2006, 24: 2098-2109. 10.1634/stemcells.2006-0009.
Zender L, Spector MS, Xue W, Flemming P, Cordon-Cardo C, Silke J, Fan ST, Luk JM, Wigler M, Hannon GJ, Mu D, Lucito R, Powers S, Lowe SW: Identification and validation of oncogenes in liver cancer using an integrative oncogenomic approach. Cell. 2006, 125: 1253-1267. 10.1016/j.cell.2006.05.030.
Gripon P, Rumin S, Urban S, Le Seyec J, Glaise D, Cannie I, Guyomard C, Lucas J, Trepo C, Guguen-Guillouzo C: Infection of a human hepatoma cell line by hepatitis B virus. Proc Natl Acad Sci USA. 2002, 99: 15655-15660. 10.1073/pnas.232137699.
Parent R, Marion MJ, Furio L, Trépo C, Petit MA: Origin and characterization of a human bipotent liver progenitor cell line. Gastroenterology. 2004, 126: 1147-1156. 10.1053/j.gastro.2004.01.002.
Cerec V, Glaise D, Garnier D, Morosan S, Turlin B, Drenou B, Gripon P, Kremsdorf D, Guguen-Guillouzo C, Corlu A: Transdifferentiation of hepatocyte-like cells from the human hepatoma HepaRG cell line through bipotent progenitor. Hepatology. 2007, 45: 957-967. 10.1002/hep.21536.
Troadec MB, Glaise D, Lamirault G, Le Cunff M, Guérin E, Le Meur N, Détivaud L, Zindy P, Leroyer P, Guisle I, Duval H, Gripon P, Théret N, Boudjema K, Guguen-Guillouzo C, Brissot P, Léger JJ, Loréal O: Hepatocyte iron loading capacity is associated with differentiation and repression of motility in the HepaRG cell line. Genomics. 2006, 87: 93-103. 10.1016/j.ygeno.2005.08.016.
Haga S, Ogawa W, Inoue H, Terui K, Ogino T, Igarashi R, Takeda K, Akira S, Enosawa S, Furukawa H, Todo S, Ozaki M: Compensatory recovery of liver mass by Akt-mediated hepatocellular hypertrophy in liver-specific STAT3-deficient mice. J Hepatol. 2005, 43: 799-807. 10.1016/j.jhep.2005.03.027.
Svegliati-Baroni G, Ridolfi F, Caradonna Z, Alvaro D, Marzioni M, Saccomanno S, Candelaresi C, Trozzi L, Macarri G, Benedetti A, Folli F: Regulation of ERK/JNK/p70S6K in two rat models of liver injury and fibrosis. J Hepatol. 2003, 39: 528-537. 10.1016/S0168-8278(03)00291-5.
Goggin MM, Nelsen CJ, Kimball SR, Jefferson LS, Morley SJ, Albrecht JH: Rapamycin-sensitive induction of eukaryotic initiation factor 4F in regenerating mouse liver. Hepatology. 2004, 40: 537-544. 10.1002/hep.20338.
Das AT, Salvado J, Boon L, Biharie G, Moorman AF, Lamers WH: Regulation of glutamate dehydrogenase expression in the developing rat liver: control at different levels in the prenatal period. Eur J Biochem. 1996, 235: 677-682. 10.1111/j.1432-1033.1996.00677.x.
Ingenuity Pathways Analysis (Ingenuity® Systems). [http://www.ingenuity.com]
Scharf JG, Dombrowski F, Ramadori G: The IGF axis and hepatocarcinogenesis. Mol Pathol. 2001, 54: 138-144. 10.1136/mp.54.3.138.
Fujiuchi N, Aglipay JA, Ohtsuka T, Maehara N, Sahin F, Su GH, Lee SW, Ouchi T: Requirement of IFI16 for the maximal activation of p53 induced by ionizing radiation. J Biol Chem. 2004, 279: 20339-20344. 10.1074/jbc.M400344200.
Rodrigues CM, Steer CJ: Fooling the Fas ligand in death. Hepatology. 2006, 44: 314-317. 10.1002/hep.21320.
Kim TY, Lee KH, Chang S, Chung C, Lee HW, Yim J, Kim TK: Oncogenic potential of a dominant negative mutant of interferon regulatory factor 3. J Biol Chem. 2003, 278: 15272-15278. 10.1074/jbc.M205792200.
O'Connor L, Strasser A, O'Reilly LA, Hausmann G, Adams JM, Cory S, Huang DC: Bim: a novel member of the Bcl-2 family that promotes apoptosis. EMBO J. 1998, 17: 384-395. 10.1093/emboj/17.2.384.
Haouzi D, Baghdiguian S, Granier G, Travo P, Mangeat P, Hibner U: Three-dimensional polarization sensitizes hepatocytes to Fas/CD95 apoptotic signalling. J Cell Sci. 2005, 118: 2763-2773. 10.1242/jcs.02403.
Danial NN, Korsmeyer SJ: Cell death: critical control points. Cell. 2004, 116: 205-219. 10.1016/S0092-8674(04)00046-7.
Foy E, Li K, Wang C, Sumpter R, Ikeda M, Lemon SM, Gale M: Regulation of interferon regulatory factor-3 by the hepatitis C virus serine protease. Science. 2003, 300: 1145-1148. 10.1126/science.1082604.
Zhang T, Lin RT, Li Y, Douglas SD, Maxcey C, Ho C, Lai JP, Wang YJ, Wan Q, Ho WZ: Hepatitis C virus inhibits intracellular interferon alpha expression in human hepatic cell lines. Hepatology. 2005, 42: 819-827. 10.1002/hep.20854.
Gao B: Cytokines, STATs and liver disease. Cell Mol Immunol. 2005, 2: 92-100.
Juan TS, Wilson DR, Wilde MD, Darlington GJ: Participation of the transcription factor C/EBP delta in the acute-phase regulation of the human gene for complement component C3. Proc Natl Acad Sci USA. 1993, 90: 2584-2588. 10.1073/pnas.90.7.2584.
Selvendiran K, Koga H, Ueno T, Yoshida T, Maeyama M, Torimura T, Yano H, Kojiro M, Sata M: Luteolin promotes degradation in signal transducer and activator of transcription 3 in human hepatoma cells: an implication for the antitumor potential of flavonoids. Cancer Res. 2006, 66: 4826-4834. 10.1158/0008-5472.CAN-05-4062.
Wan YJ, An D, Cai Y, Repa JJ, Hung-Po Chen T, Flores M, Postic C, Magnuson MA, Chen J, Chien KR, French S, Mangelsdorf DJ, Sucov HM: Hepatocyte-specific mutation establishes retinoid X receptor alpha as a heterodimeric integrator of multiple physiological processes in the liver. Mol Cell Biol. 2000, 20: 4436-4444. 10.1128/MCB.20.12.4436-4444.2000.
Cantwell CA, Sterneck E, Johnson PF: Interleukin-6-specific activation of the C/EBPdelta gene in hepatocytes is mediated by Stat3 and Sp1. Mol Cell Biol. 1998, 18: 2108-2117.
Hayashi Y, Wang W, Ninomiya T, Nagano H, Ohta K, Itoh H: Liver enriched transcription factors and differentiation of hepatocellular carcinoma. Mol Pathol. 1999, 52: 19-24.
Birkenmeier EH, Gwynn B, Howard S, Jerry J, Gordon JI, Landschulz WH, McKnight SL: Tissue-specific expression, developmental regulation, and genetic mapping of the gene encoding CCAAT/enhancer binding protein. Genes Dev. 1989, 3: 1146-1156. 10.1101/gad.3.8.1146.
Watkins PJ, Condreay JP, Huber BE, Jacobs SJ, Adams DJ: Impaired proliferation and tumorigenicity induced by CCAAT/enhancer-binding protein. Cancer Res. 1996, 56: 1063-1067.
Oliver WR, Shenk JL, Snaith MR, Russell CS, Plunket KD, Bodkin NL, Lewis MC, Winegar DA, Sznaidman ML, Lambert MH, Xu HE, Sternbach DD, Kliewer SA, Hansen BC, Willson TM: A selective peroxisome proliferator-activated receptor delta agonist promotes reverse cholesterol transport. Proc Natl Acad Sci USA. 2001, 98: 5306-5311. 10.1073/pnas.091021198.
Tanaka T, Yamamoto J, Iwasaki S, Asaba H, Hamura H, Ikeda Y, Watanabe M, Magoori K, Ioka RX, Tachibana K, Watanabe Y, Uchiyama Y, Sumi K, Iguchi H, Ito S, Doi T, Hamakubo T, Naito M, Auwerx J, Yanagisawa M, Kodama T, Sakai J: Activation of peroxisome proliferator-activated receptor delta induces fatty acid beta-oxidation in skeletal muscle and attenuates metabolic syndrome. Proc Natl Acad Sci USA. 2003, 100: 15924-15929. 10.1073/pnas.0306981100.
Wang YX, Lee CH, Tiep S, Yu RT, Ham J, Kang H, Evans RM: Peroxisome-proliferator-activated receptor delta activates fat metabolism to prevent obesity. Cell. 2003, 113: 159-170. 10.1016/S0092-8674(03)00269-1.
Lee CH, Olson P, Hevener A, Mehl I, Chong LW, Olefsky JM, Gonzalez FJ, Ham J, Kang H, Peters JM, Evans RM: PPARdelta regulates glucose metabolism and insulin sensitivity. Proc Natl Acad Sci USA. 2006, 103: 3444-3449. 10.1073/pnas.0511253103.
Pedersen TA, Bereshchenko O, Garcia-Silva S, Ermakova O, Kurz E, Mandrup S, Porse BT, Nerlov C: Distinct C/EBPalpha motifs regulate lipogenic and gluconeogenic gene expression in vivo. EMBO J. 2007, 26: 1081-1093. 10.1038/sj.emboj.7601563.
Yu VC, Delsert C, Andersen B, Holloway JM, Devary OV, Naar AM, Kim SY, Boutin JM, Glass CK, Rosenfeld MG: RXR beta: a coregulator that enhances binding of retinoic acid, thyroid hormone, and vitamin D receptors to their cognate response elements. Cell. 1991, 67: 1251-1266. 10.1016/0092-8674(91)90301-E.
Kawaguchi T, Takenoshita M, Kabashima T, Uyeda K: Glucose and cAMP regulate the L-type pyruvate kinase gene by phosphorylation/dephosphorylation of the carbohydrate response element binding protein. Proc Natl Acad Sci USA. 2001, 98: 13710-13715. 10.1073/pnas.231370798.
Fukamizu A, Takahashi S, Seo MS, Tada M, Tanimoto K, Uehara S, Murakami K: Structure and expression of the human angiotensinogen gene. Identification of a unique and highly active promoter. J Biol Chem. 1990, 265: 7576-7582.
Wilfred de Alwis NM, Day CP: Genetics of alcoholic liver disease and nonalcoholic fatty liver disease. Semin Liver Dis. 2007, 27: 44-54. 10.1055/s-2006-960170.
Gurkan C, Stagg SM, Lapointe P, Balch WE: The COPII cage: unifying principles of vesicle coat assembly. Nat Rev Mol Cell Biol. 2006, 7: 727-738. 10.1038/nrm2025.
Yi M, Sakai T, Fassler R, Ruoslahti E: Antiangiogenic proteins require plasma fibronectin or vitronectin for in vivo activity. Proc Natl Acad Sci USA. 2003, 100: 11435-11438. 10.1073/pnas.1635112100.
Krones-Herzig A, Mittal S, Yule K, Liang H, English C, Urcis R, Soni T, Adamson ED, Mercola D: Early growth response 1 acts as a tumor suppressor in vivo and in vitro via regulation of p53. Cancer Res. 2005, 65: 5133-5143. 10.1158/0008-5472.CAN-04-3742.
Monvoisin A, Bisson C, Si-Tayeb K, Balabaud C, Desmouliere A, Rosenbaum J: Involvement of matrix metalloproteinase type-3 in hepatocyte growth factor-induced invasion of human hepatocellular carcinoma cells. Int J Cancer. 2002, 97: 157-162. 10.1002/ijc.1595.
Yu FL, Liao MH, Lee JW, Shih WL: Induction of hepatoma cells migration by phosphoglucose isomerase/autocrine motility factor through the upregulation of matrix metalloproteinase-3. Biochem Biophys Res Commun. 2004, 314: 76-82. 10.1016/j.bbrc.2003.12.056.
Cao R, Bråkenhielm E, Pawliuk R, Wariaro D, Post MJ, Wahlberg E, Leboulch P, Cao Y: Angiogenic synergism, vascular stability and improvement of hind-limb ischemia by a combination of PDGF-BB and FGF-2. Nat Med. 2003, 9: 604-613. 10.1038/nm848.
Son G, Hirano T, Seki E, Iimuro Y, Nukiwa T, Matsumoto K, Nakamura T, Fujimoto J: Blockage of HGF/c-Met system by gene therapy (adenovirus-mediated NK4 gene) suppresses hepatocellular carcinoma in mice. J Hepatol. 2006, 45: 688-695. 10.1016/j.jhep.2006.04.011.
Salazar-Mather TP, Hokeness KL: Cytokine and chemokine networks: pathways to antiviral defense. Curr Top Microbiol Immunol. 2006, 303: 29-46.
Baggiolini M, Walz A, Kunkel SL: Neutrophil-activating peptide-1/interleukin 8, a novel cytokine that activates neutrophils. J Clin Invest. 1989, 84: 1045-1049.
Parent R, Kolippakkam D, Booth GC, Beretta L: Mammalian target of rapamycin activation impairs hepatocytic differentiation and targets genes moderating lipid homeostasis and hepatocellular growth. Cancer Res. 2007, 67: 4337-4345. 10.1158/0008-5472.CAN-06-3640.
Odom DT, Zizlsperger N, Gordon DB, Bell GW, Rinaldi NJ, Murray HL, Volkert TL, Schreiber J, Rolfe PA, Gifford DK, Fraenkel E, Bell GI, Young RA: Control of pancreas and liver gene expression by HNF transcription factors. Science. 2004, 303: 1378-1381. 10.1126/science.1089769.
Vaxillaire M, Rouard M, Yamagata K, Oda N, Kaisaki PJ, Boriraj VV, Chevre JC, Boccio V, Cox RD, Lathrop GM, Dussoix P, Philippe J, Timsit J, Charpentier G, Velho G, Bell GI, Froguel P: Identification of nine novel mutations in the hepatocyte nuclear factor 1 alpha gene associated with maturity-onset diabetes of the young (MODY3). Hum Mol Genet. 1997, 6: 583-586. 10.1093/hmg/6.4.583.
Yamagata K, Oda N, Kaisaki PJ, Menzel S, Furuta H, Vaxillaire M, Southam L, Cox RD, Lathrop GM, Boriraj VV, Chen X, Cox NJ, Oda Y, Yano H, Le Beau MM, Yamada S, Nishigori H, Takeda J, Fajans SS, Hattersley AT, Iwasaki N, Hansen T, Pedersen O, Polonsky KS, Turner RC, Velho G, Chèvre J-C, Froguel P, Bell GI: Mutations in the hepatocyte nuclear factor-1alpha gene in maturity-onset diabetes of the young (MODY3). Nature. 1996, 384: 455-458. 10.1038/384455a0.
Affymetrix: Technical Notes. [http://www.affymetrix.com/support/technical/technotesmain.affx]
Bioconductor: Open Source Software for Bioinformatics. [http://www.bioconductor.org]
ArrayExpress: European Bioinformatics Institute. [http://www.ebi.ac.uk/arrayexpress]
Primer Bank. [http://pga.mgh.harvard.edu/primerbank/]
Wang X, Seed B: A PCR primer bank for quantitative gene expression analysis. Nucleic Acids Res. 2003, 31: e154-10.1093/nar/gng154.
We thank Drs C Trépo and M-A Petit (INSERM Unit 871, Lyon, France) for the gift of the HepaRG cells. We also thank Deepak Kolippakkam and Neha Lohia for assistance in data analysis and Paul Farley for assistance in the preparation of the manuscript.
RP carried out the study, participated in its design and drafted the manuscript. LB conceived the study and finalized the manuscript. Both authors read and approved the final manuscript.