Global Transcriptional Response to Hfe Deficiency and Dietary Iron Overload in Mouse Liver and Duodenum

Abstract

Iron is an essential trace element whose absorption is usually tightly regulated in the duodenum. HFE-related hereditary hemochromatosis (HH) is characterized by abnormally low expression of the iron-regulatory hormone, hepcidin, which results in increased iron absorption. The liver is crucial for iron homeostasis as it is the main production site of hepcidin. The aim of this study was to explore and compare the genome-wide transcriptome response to Hfe deficiency and dietary iron overload in murine liver and duodenum. Illumina™ arrays containing over 47,000 probes were used to study global transcriptional changes. Quantitative RT-PCR (Q-RT-PCR) was used to validate the microarray results. In the liver, the expression of 151 genes was altered in Hfe−/− mice while dietary iron overload changed the expression of 218 genes. There were 173 and 108 differentially expressed genes in the duodenum of Hfe−/− mice and mice with dietary iron overload, respectively. There was 93.5% concordance between the results obtained by microarray analysis and Q-RT-PCR. Overexpression of genes for acute phase reactants in the liver and a strong induction of digestive enzyme genes in the duodenum were characteristic of the Hfe-deficient genotype. In contrast, dietary iron overload caused a more pronounced change of gene expression responsive to oxidative stress. In conclusion, Hfe deficiency caused a previously unrecognized increase in gene expression of hepatic acute phase proteins and duodenal digestive enzymes.

About the Authors:

Alejandra Rodriguez

* E-mail: [email protected]

Affiliation: Institute of Medical Technology, University of Tampere and Tampere University Hospital, Tampere, Finland

Tiina Luukkaala

Affiliation: Science Center, Pirkanmaa Hospital District and Tampere School of Public Health, University of Tampere, Tampere, Finland

Robert E. Fleming

Affiliations Department of Pediatrics, Saint Louis University School of Medicine, St. Louis, Missouri, United States of America, Edward A. Doisy Department of Biochemistry and Molecular Biology, Saint Louis University School of Medicine, St. Louis, Missouri, United States of America, Saint Louis University Liver Center, St. Louis, Missouri, United States of America

Robert S. Britton

Affiliations Saint Louis University Liver Center, St. Louis, Missouri, United States of America, Division of Gastroenterology and Hepatology, Department of Internal Medicine, Saint Louis University School of Medicine, St. Louis, Missouri, United States of America

Bruce R. Bacon

Affiliations Saint Louis University Liver Center, St. Louis, Missouri, United States of America, Division of Gastroenterology and Hepatology, Department of Internal Medicine, Saint Louis University School of Medicine, St. Louis, Missouri, United States of America

Seppo Parkkila

Affiliations Institute of Medical Technology, University of Tampere and Tampere University Hospital, Tampere, Finland, School of Medicine, University of Tampere and Tampere University Hospital, Tampere, Finland

Introduction

Iron plays crucial roles in cellular metabolism but, in excess, it can catalyze the formation of free radicals leading to oxidative stress and cell damage [1]. Iron is absorbed in the duodenum, where it crosses the apical and basolateral membranes of absorptive enterocytes to enter the blood stream [2]. There is no regulated mechanism of iron excretion, and thus the absorption of iron must be tightly regulated to maintain iron balance. HFE-related hereditary hemochromatosis (HH, OMIM-235200) is an autosomal recessive disorder in which absorption of iron is inappropriately high [3], [4]. HH is characterized by high transferrin saturation and low iron content in macrophages. Iron is deposited primarily in the parenchymal cells of various organs, particularly the liver, but also the pancreas, heart, skin, and testes, resulting in tissue damage and organ failure. Clinical complications in untreated HH patients include hepatic fibrosis, cirrhosis, hepatocellular carcinoma, diabetes, cardiomyopathy, hypogonadism, and arthritis [4].

HH is characterized by inappropriately low expression of the iron-regulatory hormone hepcidin [5]. Hepcidin, a small peptide hormone expressed mainly in the liver, is a central player in the maintenance of iron balance [6]. The only known molecule capable of transporting iron out of cells is ferroportin [7]–[9]. This iron exporter is located in the plasma membrane of enterocytes, reticuloendothelial cells, hepatocytes, and placental cells [7]. Hepcidin binds to ferroportin and induces its internalization and degradation, therefore suppressing the transport of iron into the circulation [10]. The expression of hepcidin is induced by increased iron stores and inflammation, and suppressed by hypoxia and anemia [11], [12].

Mice homozygous for a null allele of Hfe (Hfe−/−) provide a genetic animal model of HH [13]. There are several animal models of iron overload based on administration of exogenous iron [14]. According to the route of iron delivery, these can be divided into two main types: enteral (i.e. dietary) and parenteral models. For example, dietary supplementation with carbonyl iron in mice reproduces the HH pattern of hepatic iron loading, with predominantly parenchymal iron deposition [14]. Although both Hfe−/− mice and carbonyl iron-fed mice develop iron overload, there are important differences between these two models. Hfe−/− mice lack Hfe protein and therefore have decreased expression of hepcidin [15], [16], while mice with dietary iron overload express functional Hfe protein and their hepcidin expression is elevated [12].

Current RNA microarray technology allows expression profiling of the whole transcriptome. This methodology has been used to explore the effects of Hfe gene disruption on mRNA expression in the liver and duodenum, two organs with crucial roles in iron metabolism [17]. In the present study, we used this approach to study gene expression in the liver and duodenum of Hfe−/− mice and wild-type mice, with or without dietary iron overload. This allowed the identification of genes whose expression is changed during iron overload and those genes whose expression is differentially influenced by lack of Hfe protein.

Results

We used global microarray analysis to study gene expression in the liver and duodenum of Hfe−/− mice and carbonyl iron loaded mice, and comparing it with that of wild-type mice fed a standard diet. This approach allowed the identification of genes whose expression is changed during iron overload and those genes whose expression is differentially influenced by lack of Hfe protein. All the mice used were males and all had the same genetic background (C57BL/6).

Hepatic transcriptional response to Hfe deficiency and dietary iron overload

Hepatic RNA from 3 Hfe−/− mice and 2 wild-type mice was subjected to microarray analysis. The Pearson correlation coefficient between the knock out mice and between the controls was in both cases 0.989. The results revealed 86 induced genes and 65 repressed genes, using a cutoff value of ±1.4-fold (Table 1 and Dataset S1). This cutoff value has been proposed as an adequate compromise above which there is a high correlation between microarray and Q-RT-PCR data, regardless of other factors such as spot intensity and cycle threshold [18]. The fold-changes ranged from 9.83 to −3.47. Functional annotation of the gene lists highlighted the biological processes that may be modified by Hfe deficiency. This analysis revealed enrichment of heat shock proteins and proteins related to inflammatory responses or antigen processing and presentation, among others (Table 2).

[Figure omitted. See PDF.]

Table 1. Number of genes regulated by Hfe deficiency or dietary iron overload in murine liver and duodenum.

https://doi.org/10.1371/journal.pone.0007212.t001

[Figure omitted. See PDF.]

Table 2. Functional annotation of genes regulated in the liver of Hfe−/− mice.

https://doi.org/10.1371/journal.pone.0007212.t002

Another microarray experiment was performed using hepatic RNA from 3 mice with dietary iron overload and 2 mice fed a standard diet. The similarity between samples from individual mice was measured as the Pearson correlation coefficient, which was 0.989 between iron overloaded mice and 0.991 between control mice. The expression of 123 genes was upregulated and that of 95 genes was downregulated, applying a cutoff value of ±1.4-fold (Table 1 and Dataset S2). The fold-changes ranged between 13.58 and −7.46. The list of regulated genes was functionally annotated (Table 3), showing enrichment of cytochrome P450 proteins as well as others involved in glutathione metabolism, acute-phase response, organic acid biosynthetic process and cellular iron homeostasis, among others.

[Figure omitted. See PDF.]

Table 3. Functional annotation of genes regulated in the liver of iron-fed mice.

https://doi.org/10.1371/journal.pone.0007212.t003

There were 11 upregulated and 7 downregulated genes that were affected by both Hfe deficiency and dietary iron overload in similar fashion, while 27 genes were regulated in opposite directions by these two conditions in the liver (Table 4). In some cases, several genes belonging to the same gene family showed divergent regulation (e.g., Saa1, Saa2, Saa3) with upregulation in Hfe−/− mice and downregulation by dietary iron overload.

[Figure omitted. See PDF.]

Table 4. Comparison of hepatic gene regulation by Hfe deficiency or dietary iron overload.

https://doi.org/10.1371/journal.pone.0007212.t004

Altered expression of iron-related genes in the liver.

The expression of 3 iron-related genes was altered in the liver of Hfe−/− mice. The expression of Hamp1 and Tfrc was decreased and that of Lcn2 was induced. We confirmed these results using Q-RT-PCR, and also tested the expression of Hamp2, which was downregulated (Figure 1). Dietary iron overload changed the expression of 5 iron-related genes in the liver. The expression of Hamp1, Hamp2, Lcn2 and Cp were upregulated using both microarray analysis and Q-RT-PCR, while Tfrc expression was down-regulated by 1.7-fold (Figure 2).

[Figure omitted. See PDF.]

Figure 1. Validation of liver microarray data from Hfe−/− mice by Q-RT-PCR.

The expression of various mRNA species in 5 Hfe−/− mice is compared to those in 4 wild-type controls. Each sample was run in triplicate. (mean±SD). *p<0.05; **p<0.025; ***p<0.01.

https://doi.org/10.1371/journal.pone.0007212.g001

[Figure omitted. See PDF.]

Figure 2. Expression of genes affected by dietary iron overload in the liver, as confirmed by Q-RT-PCR.

Samples from 5 mice fed an iron-supplemented diet and 4 mice fed a control diet were used, and each sample was run in triplicate. (mean±SD). *p<0.05; **p<0.025; ***p<0.01.

https://doi.org/10.1371/journal.pone.0007212.g002

Confirmation of hepatic microarray results by Q-RT-PCR.

Microarray analysis for the expression of several genes was confirmed by performing Q-RT-PCR on hepatic samples from 5 Hfe−/− mice, 4 wild-type control mice, 5 iron-fed mice and 4 mice fed a standard diet. For this purpose, we selected iron-related genes and others whose expression was substantially altered in the experimental groups. A total of 29 results from the hepatic microarray data, corresponding to 24 different genes, were tested by Q-RT-PCR, and 27 (93.1%) of them showed concordant results by these two methods (Figures 1 and 2). Changes in Foxq1 and Dmt1 expression were false-positives in the microarray analysis for Hfe−/− mice and dietary iron overload, respectively. The upregulation of Ltf expression by dietary iron overload observed by microarray analysis could not be confirmed by Q-RT-PCR because the expression levels in samples from all but one of the treated mice and all control mice were below the detection threshold.

Duodenal gene expression response to Hfe deficiency and dietary iron supplementation

Microarray analysis of duodenal RNA from 2 Hfe−/− mice and 2 wild-type mice revealed that the expression of 143 genes was upregulated and that of 30 genes was downregulated when a cutoff value of ±1.4-fold was used (Table 1 and Dataset S3). The fold-changes ranged from 15.67 to −3.14. The Pearson correlation coefficient between knockout mice and between controls was 0.976 and 0.971, respectively. Functional categories overrepresented among the genes regulated by Hfe deficiency included proteins with endopeptidase activity, and others involved in lipid catabolism and antimicrobial activity (Table 5).

[Figure omitted. See PDF.]

Table 5. Functional annotation of genes regulated in the duodenum of Hfe−/− mice.

https://doi.org/10.1371/journal.pone.0007212.t005

Global transcriptional regulation was also studied in the duodenum of mice fed an iron-supplemented diet, using 3 treated mice and 2 controls. The Pearson correlation coefficient was 0.985 between treated mice and 0.983 between controls. The expression of 49 genes was induced and 59 genes were repressed, applying a cutoff value of ±1.4-fold (Table 1 and Dataset S4). The fold-changes ranged between 6.07 and −5.64. Functional annotation of the gene list evidenced enrichment of genes involved in glutathione metabolism, antigen processing and presentation and inflammatory response, among others (Table 6).

[Figure omitted. See PDF.]

Table 6. Functional annotation of genes regulated in the duodenum of mice fed an iron-supplemented diet.

https://doi.org/10.1371/journal.pone.0007212.t006

We identified genes whose expression was affected by both Hfe deficiency and dietary iron supplementation in the duodenum. There were 4 genes whose expression was induced in both conditions, 3 genes whose expression was decreased, and 4 genes with opposite regulation (Table 7).

[Figure omitted. See PDF.]

Table 7. Genes regulated in the duodenum of mice by Hfe deficiency or iron-supplemented diet.

https://doi.org/10.1371/journal.pone.0007212.t007

Altered expression of iron-related genes in the duodenum.

In the duodenum of Hfe−/− mice, Hamp2 expression was increased by 2.7-fold using microarray analysis. However, this could not be confirmed by Q-RT-PCR, because Hamp2 mRNA levels in the samples from wild-type mice and in one Hfe−/− sample were below the detection threshold. In mice fed the iron-supplemented diet, the duodenal expression of Tfrc was downregulated and that for Hmox1 was upregulated: both of these results were validated by Q-RT-PCR (Figure 3).

[Figure omitted. See PDF.]

Figure 3. Expression of genes regulated in the duodenum of dietary iron-loaded mice as verified by Q-RT-PCR.

Samples from 5 mice fed an iron-supplemented diet and 4 mice fed a control diet were used, and each sample was run in triplicate. (mean±SD). *p<0.05; **p<0.025; ***p<0.01.

https://doi.org/10.1371/journal.pone.0007212.g003

Confirmation of duodenal microarray results by Q-RT-PCR.

Q-RT-PCR analyses were done on duodenal RNA samples from 5 Hfe−/− mice, 4 wild-type control mice, 5 iron-fed mice and 4 mice fed a standard diet in order to confirm the microarray results. The mRNA expression of a total of 17 different genes was tested and 16 (94.1%) showed concordant results between microarray analysis and Q-RT-PCR (Figures 3 and 4). The sole discrepant result concerned the expression of Ddb1 that was downregulated according to microarray analysis, while Q-RT-PCR revealed a slight induction (1.25-fold) of expression.

[Figure omitted. See PDF.]

Figure 4. Validation of the duodenal microarray results from Hfe−/− mice by Q-RT-PCR.

The Hfe−/− and control groups contained samples from 4 mice, and each sample was tested in triplicate. (mean±SD). *p<0.05; **p<0.025; ***p<0.01.

https://doi.org/10.1371/journal.pone.0007212.g004

Discussion

The goal of this study was to explore and compare the genome-wide transcriptome response to Hfe deficiency and dietary iron overload in murine liver and duodenum. This approach allowed the identification of genes whose expression is changed during iron overload and those genes whose expression is differentially influenced by lack of Hfe protein. The global transcriptional response to Hfe deficiency has been explored previously in the liver and duodenum of two mouse strains [17]. However, it is notable that only a few analogous changes in gene expression are seen when comparing our data with those of the previous study, even for mice of the same genetic background. Two other reports have explored expression of selected genes by using dedicated arrays in Hfe−/− mice and in mice with secondary iron overload produced by intraperitoneal injection of iron-dextran [19], [20]. In one study, duodenum and liver samples were analyzed using an array of iron-related genes [19]. The results for duodenal gene expression in Hfe−/− mice have no concordance with ours. Regulation of hepatic gene expression, on the other hand, is similar for several genes, such as Hamp1, Tfrc and Mt1. The second report focused on gene expression in the duodenum [20], and again, there is little concordance between their observations and ours. The lack of agreement between these studies is probably due to differences in the animal models (parenteral vs. enteral iron loading; mouse strains) and in the microarray methodology.

The hepatic expression of acute phase proteins (APPs) can be induced by inflammatory mediators such as interleukin-6. Interestingly, the liver of Hfe−/− mice has upregulated expression of APPs such as serum amyloids, lipocalins and orosomucoids. Notably, the expression of serum amyloid genes (Saa1, Saa2, Saa3) was upregulated in the Hfe−/− mice compared to being downregulated in dietary iron overload, suggesting that Hfe deficiency induces this gene expression by an iron-independent mechanism. However, hepatic interleukin-6 mRNA expression was not significantly changed by Hfe deficiency, so the potential involvement of this cytokine in the observed upregulation of APPs remains uncertain.

Lipocalin2 (human Ngal from neutrophil gelatinase-associated lipocalin) is an APP with antimicrobial properties through a mechanism of iron deprivation by siderophore binding [21]. It can donate iron to various types of cells [22], [23] and seems to be capable of intracellular iron chelation and iron excretion [24]. Furthermore, a recent study has shown that lipocalin2 is an adipokine with potential importance in insulin resistance associated with obesity [25]. We observed that Lcn2 expression is increased in the liver of both Hfe−/− mice and those with dietary iron overload, suggesting that this induction is iron-related.

Dietary iron overload of the liver led to increased expression of both hepcidin genes (Hamp1, Hamp2) as previously reported [26], [27], and these results were verified by Q-RT-PCR. In the liver of Hfe−/− mice, Hamp1 expression was downregulated as expected [15], [16], [19]. We also examined the levels of Hamp2 mRNA by Q-RT-PCR and found a -1.92-fold change. The low expression of hepatic Hamp1 in Hfe−/− mice is likely responsible for the increased iron absorption and low microphage iron content in these mice [15], [16], [19].

Inhibitor of DNA-binding/differentiation proteins, also known as Id proteins, comprise a family of proteins that heterodimerize with basic-helix-loop-helix (bHLH) transcription factors to inhibit their binding to DNA. Several studies have reported that Id proteins have important roles in differentiation, cell cycle and angiogenesis in various cell types [28]. Expression of Id1, 2, and 3 is increased during liver disease, with levels that escalate as liver disease progresses from hepatitis to cirrhosis. In hepatocellular carcinoma, high expression is observed in well-differentiated tumors, and it decreases as the tumor cells become undifferentiated [29]. In light of these findings, it has been suggested that Id1, 2, and 3 may play a role in the early stages of hepatocarcinogenesis. Given that, it is notable that we found that the expression of Id1, 2, 3, and 4 was increased in the liver of mice with dietary iron overload, but was unaffected in Hfe−/− mice. Increased hepatic expression of Id1 mRNA has previously been reported in mice fed an iron-supplemented diet [30]. The same study showed upregulation of the gene for bone morphogenetic protein 6 (Bmp6) in the same experimental mice. Recent work demonstrates that Bmp6 is a key player in the signalling pathway that controls hepcidin expression [31]. Unexpectedly, upregulation of hepatic Bmp6 mRNA expression by dietary iron overload was not evident in the current study.

The gene expression of several heat shock proteins was downregulated in the liver and duodenum by both Hfe deficiency and dietary iron overload, with a considerably greater number of these genes downregulated in the liver of Hfe−/− mice. Although these genes are induced under certain stress conditions, such as heat shock and ischemia–reperfusion, their expression is decreased by iron overload [19], [27], [32]. Currently, the physiological implications of this downregulation are unknown.

Our results indicate that disruption of the Hfe gene induces the expression of many genes in the duodenum coding for digestive enzymes, such as elastases, carboxypeptidases, trypsins, chymotrypsins, amylases, and lipases. In contrast, feeding mice with an iron-supplemented diet did not affect the expression of any of these genes. The upregulation of gene expression for digestive enzymes in Hfe−/− mice is surprising because overexpression of these enzymes has not been associated with HH.

A common feature of the duodenal response to both Hfe deficiency and dietary iron overload was the transcriptional repression of genes involved in antimicrobial activities, such as cryptdins. In mice fed an iron-supplemented diet, there was also a decrease in mRNA expression for genes involved in antigen processing and presentation, such as some genes of the MHC class II family.

The solute carrier molecules constitute a large family of proteins involved in membrane transport of diverse molecules. The gene expression of many family members was affected by Hfe deficiency or dietary iron overload. In the duodenum, the expression of the sodium-coupled neutral amino acid transporter Slc38a5 was induced in Hfe−/− mice and repressed in mice fed an iron-supplemented diet. In the liver, the expression of Slc46a3 was upregulated in Hfe−/− mice. This gene belongs to the Slc46 subfamily of heme transporters. It is thus a close relative of Slc46a1 (also known as HCP1), a recently identified, although controversial, heme transporter [33], [34]. The iron transporter Dmt1, encoded by Slc11a2, contains an iron-responsive element (IRE) in the 3′UTR of its mRNA. This permits the regulation of Dmt1 mRNA levels according to the cellular labile iron pool by mediation of the iron regulatory proteins, IRP1 and IRP2. Under iron-replete conditions, IRP activity is reduced rendering the Dmt1 mRNA vulnerable to degradation. The opposite is true under iron-deficient conditions, which is believed to be the situation inside the enterocytes of HH patients and Hfe−/− mice [35], [36]. Accordingly, in some studies, increased expression of Dmt1 has been observed in the duodenum of HH patients [37] as well as in Hfe−/− mice [38]. However, we did not find a significant change in the expression of Dmt1 in the Hfe−/− duodenum. This may be explained by the inability of our microarray probes and PCR primers to discriminate between IRE-positive and IRE–negative transcripts.

The post-transcriptional regulation of Tfrc (transferrin receptor 1) by iron is also mediated through the IRE/IRP system [39]. Tfrc is involved in the uptake of transferrin-bound iron by cells. Analogous to our observations, suppression of Tfrc expression in the duodenum [40] and liver [27] of mice fed an iron-supplemented diet, and in the liver of Hfe−/− mice [19] has been reported previously. Our microarray analysis indicates that the expression of Tfrc was not significantly changed in the duodenum of Hfe−/− mice, a result that agrees with a previous report [19].

Excess free iron increases oxidant production [1]. Subsequently, some antioxidant defense mechanisms are upregulated in order to provide resistance to iron-related toxicity. It is notable from our data that this response is elicited in both liver and duodenum, as seen in the upregulation of glutathione S-transferase genes. Interestingly, dietary iron overload seems to induce a stronger response than Hfe deficiency, especially in the regulation of enzymes involved in glutathione-related detoxification of reactive intermediates.

In conclusion, Hfe deficiency results in increased gene expression of hepatic APPs and duodenal digestive enzymes. In contrast, dietary iron overload causes a more pronounced change of gene expression responsive to oxidative stress.

Materials and Methods

Ethics Statement

The animal protocols were approved by the Animal Care and Use Committees of Saint Louis University and the University of Oulu (permission No 102/05).

Animal care and animal models

Five male C57BL/6 mice homozygous for a disruption of the Hfe gene and 4 male wild-type control mice were fed a standard rodent diet (250 ppm of iron) and sacrificed at approximately 10 weeks of age. The generation of the Hfe−/− mice has been described elsewhere [13]. In addition, 5 male C57BL/6 mice fed an iron-supplemented diet (2% carbonyl iron) and 4 male control mice fed a standard diet (200 ppm of iron) for 6 weeks were used [27]. The mice with dietary iron overload had a hepatic iron concentration that was approximately 2.5 times higher than the Hfe−/− mice. The duodenum and liver samples were immediately collected from anesthetized mice and immersed in RNAlater solution (Ambion, Huntingdon, UK).

RNA isolation

Total RNA extraction and quality control have been described previously [27].

Microarray analysis

All microarray data reported in the present article are described in accordance with MIAME guidelines, have been deposited in NCBI's Gene Expression Omnibus public repository [41], and are accessible through GEO Series accession number GSE17969 [42]. Microarray experiments were performed in the Finnish DNA Microarray Centre at Turku Centre for Biotechnology using Illumina's Sentrix Mouse-6 Expression Beadchips. Duodenal and liver RNA samples from 3 Hfe−/− mice and 3 mice with dietary iron overload were used. As controls, RNA samples from the duodenum and liver of 4 wild-type mice (2 controls of the Hfe−/− mice and 2 controls of the mice with dietary iron overload) were used. All 10 samples were analyzed individually. The amplification of total RNA (300 ng), in vitro transcription, hybridization and scanning have been described before [27].

Data analysis

Array data were normalized with Chipster (v1.1.1) using the quantile normalization method. Quality control of the data included non-metric multidimensional scaling, dendrograms, hierarchical clustering, and 2-way clustering (heat maps). These analyses showed that data from one of the three duodenal samples from Hfe−/− mice were highly divergent from the other two. Thus, this sample was excluded from further analyses. The data were then filtered according to the SD of the probes. The percentage of data that did not pass through the filter was adjusted to 99.4%, implicating a SD value of almost 3. At this point, statistical analysis was performed using the empirical Bayes t-test for the comparison of 2 groups. Due to the small number of samples, the statistical results were considered as orientative and thus no filtering was applied to the data according to p-values. The remaining 280 probes were further filtered according to fold-change with ±1.4 as cut-off values for up- and down-regulated expression, respectively. The functional annotation tool DAVID (Database for Annotation, Visualization and Integrated Discovery) [43], [44] was used to identify enriched biological categories among the regulated genes as compared to all the genes present in Illumina's Sentrix Mouse-6 Expression Beadchip. The annotation groupings analyzed were: Gene Ontology biological process and molecular functions, SwissProt Protein Information Resources keywords, SwissProt comments, Kyoto Encyclopedia of Genes and Genomes and Biocarta pathways. Results were filtered to remove categories with EASE (expression analysis systematic explorer) scores greater than 0.05. Redundant categories with the same gene members were removed to yield a single representative category.

Quantitative Reverse-Transcriptase PCR

For this analysis, duodenal and liver RNA samples from 5 mice of each experimental group (Hfe−/− and dietary iron overload) and 4 mice from each control group (wild-type and normal diet) were used. Exceptionally, for the analysis of mRNA expression in the duodenum of Hfe−/− mice, only 4 samples were used. RNA samples (5 µg) were converted into first strand cDNA with a First Strand cDNA Synthesis kit (Fermentas, Burlington, Canada) using random hexamer primers. The relative expression levels of target genes in the duodenum and liver were assessed by Q-RT-PCR using the LightCycler detection system (Roche, Rotkreuz, Switzerland). The reaction setup, cycling program, standard curve method and primer pairs for Angptl4, Dnajb1 and Tfrc have been described before [27]. Mouse Hamp1 and Hamp2 primers have also been characterized previously [26]. The primer sets for the other target genes (Dataset S5) were designed using Primer3 [45], based on the complete cDNA sequences deposited in GenBank. The specificity of the primers was verified using NCBI Basic Local Alignment and Search Tool (BLAST) [46]. To avoid amplification of contaminating genomic DNA, both primers from each set were specific to different exons, when possible. Each cDNA sample was tested in triplicate. The mean and SD of the 3 crossing point (Cp) values were calculated for each sample and a SD cutoff level of 0.2 was set. Accordingly, when the SD of the triplicates of a sample was greater than 0.2, the most outlying replicate was excluded and the analysis was continued with the two remaining replicates. Using the standard curve method, the Cp values were then transformed by the LightCycler software into copy numbers. The expression value for each sample was the mean of the copy numbers for the sample's replicates. This value was normalized by dividing it by the geometric mean of the 4 internal control genes, an accurate normalization method [47]. The normalization factor was always considered as a value of 100 and the final result was expressed as relative mRNA expression level.

Statistical analyses

We performed statistical analyses of the microarray data using the empirical Bayes t-test for the comparison of 2 groups, and the p-values are shown in supplementary datasets S1-S4. For the Q-RT-PCR results, we used the Mann-Whitney test to evaluate differences in group values for Hfe−/− mice vs. wild-type mice and mice with dietary iron overload vs. untreated mice. Due to the small sample sizes, the statistical significance is only considered as orientative. Values are expressed as mean±SD.

Supporting Information

[Figure omitted. See PDF.]

Dataset S1.

List of genes differentially expressed in the liver of Hfe knockout mice

https://doi.org/10.1371/journal.pone.0007212.s001

(0.04 MB XLS)

Dataset S2.

List of genes differentially expressed in the liver of mice fed an iron-supplemented diet

https://doi.org/10.1371/journal.pone.0007212.s002

(0.05 MB XLS)

Dataset S3.

Genes whose expression was altered in the duodenum of Hfe knockout mice

https://doi.org/10.1371/journal.pone.0007212.s003

(0.05 MB XLS)

Dataset S4.

Genes whose expression was affected in the duodenum of mice fed an iron-supplemented diet

https://doi.org/10.1371/journal.pone.0007212.s004

(0.04 MB XLS)

Dataset S5.

Sequences of the primers used in the Q-RT-PCR experiments performed in this study

https://doi.org/10.1371/journal.pone.0007212.s005

(0.06 MB DOC)

Acknowledgments

We thank Mary Migas and Rosemary O'Neill for excellent technical assistance. We are grateful to the staff of the Finnish microarray core facility at Turku Centre of Biotechnology, especially Tiia Heinonen for performing the microarray hybridizations and scanning.

Author Contributions

Conceived and designed the experiments: AR REF RSB BRB SP. Performed the experiments: AR. Analyzed the data: AR TL. Contributed reagents/materials/analysis tools: REF RSB BRB SP. Wrote the paper: AR. Critically reviewed the manuscript and approved its final version: TL REF RSB BRB SP.

Citation: Rodriguez A, Luukkaala T, Fleming RE, Britton RS, Bacon BR, Parkkila S (2009) Global Transcriptional Response to Hfe Deficiency and Dietary Iron Overload in Mouse Liver and Duodenum. PLoS ONE4(9): e7212. https://doi.org/10.1371/journal.pone.0007212

References

1. Britton RS (1996) Metal-induced hepatotoxicity. Semin Liver Dis 16: 3–12.RS Britton1996Metal-induced hepatotoxicity.Semin Liver Dis16312

2. Parkkila S, Niemela O, Britton RS, Fleming RE, Waheed A, et al. (2001) Molecular aspects of iron absorption and HFE expression. Gastroenterology 121: 1489–1496.S. ParkkilaO. NiemelaRS BrittonRE FlemingA. Waheed2001Molecular aspects of iron absorption and HFE expression.Gastroenterology12114891496

3. Feder JN, Gnirke A, Thomas W, Tsuchihashi Z, Ruddy DA, et al. (1996) A novel MHC class I-like gene is mutated in patients with hereditary haemochromatosis. Nat Genet 13: 399–408.JN FederA. GnirkeW. ThomasZ. TsuchihashiDA Ruddy1996A novel MHC class I-like gene is mutated in patients with hereditary haemochromatosis.Nat Genet13399408

4. Fleming RE, Britton RS, Waheed A, Sly WS, Bacon BR (2005) Pathophysiology of hereditary hemochromatosis. Semin Liver Dis 25: 411–419.RE FlemingRS BrittonA. WaheedWS SlyBR Bacon2005Pathophysiology of hereditary hemochromatosis.Semin Liver Dis25411419

5. Papanikolaou G, Samuels ME, Ludwig EH, MacDonald ML, Franchini PL, et al. (2004) Mutations in HFE2 cause iron overload in chromosome 1q-linked juvenile hemochromatosis. Nat Genet 36: 77–82.G. PapanikolaouME SamuelsEH LudwigML MacDonaldPL Franchini2004Mutations in HFE2 cause iron overload in chromosome 1q-linked juvenile hemochromatosis.Nat Genet367782

6. Nemeth E, Ganz T (2006) Regulation of iron metabolism by hepcidin. Annu Rev Nutr 26: 323–342.E. NemethT. Ganz2006Regulation of iron metabolism by hepcidin.Annu Rev Nutr26323342

7. Abboud S, Haile DJ (2000) A novel mammalian iron-regulated protein involved in intracellular iron metabolism. J Biol Chem 275: 19906–19912.S. AbboudDJ Haile2000A novel mammalian iron-regulated protein involved in intracellular iron metabolism.J Biol Chem2751990619912

8. McKie AT, Marciani P, Rolfs A, Brennan K, Wehr K, et al. (2000) A novel duodenal iron-regulated transporter, IREG1, implicated in the basolateral transfer of iron to the circulation. Mol Cell 5: 299–309.AT McKieP. MarcianiA. RolfsK. BrennanK. Wehr2000A novel duodenal iron-regulated transporter, IREG1, implicated in the basolateral transfer of iron to the circulation.Mol Cell5299309

9. Gunshin H, Mackenzie B, Berger UV, Gunshin Y, Romero MF, et al. (1997) Cloning and characterization of a mammalian proton-coupled metal-ion transporter. Nature 388: 482–488.H. GunshinB. MackenzieUV BergerY. GunshinMF Romero1997Cloning and characterization of a mammalian proton-coupled metal-ion transporter.Nature388482488

10. Nemeth E, Tuttle MS, Powelson J, Vaughn MB, Donovan A, et al. (2004) Hepcidin regulates cellular iron efflux by binding to ferroportin and inducing its internalization. Science 306: 2090–2093.E. NemethMS TuttleJ. PowelsonMB VaughnA. Donovan2004Hepcidin regulates cellular iron efflux by binding to ferroportin and inducing its internalization.Science30620902093

11. Nicolas G, Chauvet C, Viatte L, Danan JL, Bigard X, et al. (2002) The gene encoding the iron regulatory peptide hepcidin is regulated by anemia, hypoxia, and inflammation. J Clin Invest 110: 1037–1044.G. NicolasC. ChauvetL. ViatteJL DananX. Bigard2002The gene encoding the iron regulatory peptide hepcidin is regulated by anemia, hypoxia, and inflammation.J Clin Invest11010371044

12. Pigeon C, Ilyin G, Courselaud B, Leroyer P, Turlin B, et al. (2001) A new mouse liver-specific gene, encoding a protein homologous to human antimicrobial peptide hepcidin, is overexpressed during iron overload. J Biol Chem 276: 7811–7819.C. PigeonG. IlyinB. CourselaudP. LeroyerB. Turlin2001A new mouse liver-specific gene, encoding a protein homologous to human antimicrobial peptide hepcidin, is overexpressed during iron overload.J Biol Chem27678117819

13. Zhou XY, Tomatsu S, Fleming RE, Parkkila S, Waheed A, et al. (1998) HFE gene knockout produces mouse model of hereditary hemochromatosis. Proc Natl Acad Sci U S A 95: 2492–2497.XY ZhouS. TomatsuRE FlemingS. ParkkilaA. Waheed1998HFE gene knockout produces mouse model of hereditary hemochromatosis.Proc Natl Acad Sci U S A9524922497

14. Ramm GA (2000) Animal models of iron overload based on excess exogenous iron. In: Barton JC, Edwards CQ, editors. Hemochromatosis: Genetics, Pathophysiology, Diagnosis and Treatment. New York: Cambridge University Press. pp. 494–507.GA Ramm2000Animal models of iron overload based on excess exogenous iron.JC BartonCQ EdwardsHemochromatosis: Genetics, Pathophysiology, Diagnosis and TreatmentNew YorkCambridge University Press494507

15. Bridle KR, Frazer DM, Wilkins SJ, Dixon JL, Purdie DM, et al. (2003) Disrupted hepcidin regulation in HFE-associated haemochromatosis and the liver as a regulator of body iron homoeostasis. Lancet 361: 669–673.KR BridleDM FrazerSJ WilkinsJL DixonDM Purdie2003Disrupted hepcidin regulation in HFE-associated haemochromatosis and the liver as a regulator of body iron homoeostasis.Lancet361669673

16. Ahmad KA, Ahmann JR, Migas MC, Waheed A, Britton RS, et al. (2002) Decreased liver hepcidin expression in the Hfe knockout mouse. Blood Cells Mol Dis 29: 361–366.KA AhmadJR AhmannMC MigasA. WaheedRS Britton2002Decreased liver hepcidin expression in the Hfe knockout mouse.Blood Cells Mol Dis29361366

17. Coppin H, Darnaud V, Kautz L, Meynard D, Aubry M, et al. (2007) Gene expression profiling of Hfe−/− liver and duodenum in mouse strains with differing susceptibilities to iron loading: identification of transcriptional regulatory targets of Hfe and potential hemochromatosis modifiers. Genome Biol 8: R221.H. CoppinV. DarnaudL. KautzD. MeynardM. Aubry2007Gene expression profiling of Hfe−/− liver and duodenum in mouse strains with differing susceptibilities to iron loading: identification of transcriptional regulatory targets of Hfe and potential hemochromatosis modifiers.Genome Biol8R221

18. Morey JS, Ryan JC, Van Dolah FM (2006) Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR. Biol Proced Online 8: 175–193.JS MoreyJC RyanFM Van Dolah2006Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR.Biol Proced Online8175193

19. Muckenthaler M, Roy CN, Custodio AO, Minana B, deGraaf J, et al. (2003) Regulatory defects in liver and intestine implicate abnormal hepcidin and Cybrd1 expression in mouse hemochromatosis. Nat Genet 34: 102–107.M. MuckenthalerCN RoyAO CustodioB. MinanaJ. deGraaf2003Regulatory defects in liver and intestine implicate abnormal hepcidin and Cybrd1 expression in mouse hemochromatosis.Nat Genet34102107

20. Abgueguen E, Toutain B, Bedrine H, Chicault C, Orhant M, et al. (2006) Differential expression of genes related to HFE and iron status in mouse duodenal epithelium. Mamm Genome 17: 430–450.E. AbgueguenB. ToutainH. BedrineC. ChicaultM. Orhant2006Differential expression of genes related to HFE and iron status in mouse duodenal epithelium.Mamm Genome17430450

21. Goetz DH, Holmes MA, Borregaard N, Bluhm ME, Raymond KN, et al. (2002) The neutrophil lipocalin NGAL is a bacteriostatic agent that interferes with siderophore-mediated iron acquisition. Mol Cell 10: 1033–1043.DH GoetzMA HolmesN. BorregaardME BluhmKN Raymond2002The neutrophil lipocalin NGAL is a bacteriostatic agent that interferes with siderophore-mediated iron acquisition.Mol Cell1010331043

22. Yang J, Mori K, Li JY, Barasch J (2003) Iron, lipocalin, and kidney epithelia. Am J Physiol Renal Physiol 285: F9–18.J. YangK. MoriJY LiJ. Barasch2003Iron, lipocalin, and kidney epithelia.Am J Physiol Renal Physiol285F918

23. Devireddy LR, Gazin C, Zhu X, Green MR (2005) A cell-surface receptor for lipocalin 24p3 selectively mediates apoptosis and iron uptake. Cell 123: 1293–1305.LR DevireddyC. GazinX. ZhuMR Green2005A cell-surface receptor for lipocalin 24p3 selectively mediates apoptosis and iron uptake.Cell12312931305

24. Mori K, Lee HT, Rapoport D, Drexler IR, Foster K, et al. (2005) Endocytic delivery of lipocalin-siderophore-iron complex rescues the kidney from ischemia-reperfusion injury. J Clin Invest 115: 610–621.K. MoriHT LeeD. RapoportIR DrexlerK. Foster2005Endocytic delivery of lipocalin-siderophore-iron complex rescues the kidney from ischemia-reperfusion injury.J Clin Invest115610621

25. Yan QW, Yang Q, Mody N, Graham TE, Hsu CH, et al. (2007) The adipokine lipocalin 2 is regulated by obesity and promotes insulin resistance. Diabetes 56: 2533–2540.QW YanQ. YangN. ModyTE GrahamCH Hsu2007The adipokine lipocalin 2 is regulated by obesity and promotes insulin resistance.Diabetes5625332540

26. Ilyin G, Courselaud B, Troadec MB, Pigeon C, Alizadeh M, et al. (2003) Comparative analysis of mouse hepcidin 1 and 2 genes: evidence for different patterns of expression and co-inducibility during iron overload. FEBS Lett 542: 22–26.G. IlyinB. CourselaudMB TroadecC. PigeonM. Alizadeh2003Comparative analysis of mouse hepcidin 1 and 2 genes: evidence for different patterns of expression and co-inducibility during iron overload.FEBS Lett5422226

27. Rodriguez A, Hilvo M, Kytomaki L, Fleming RE, Britton RS, et al. (2007) Effects of iron loading on muscle: genome-wide mRNA expression profiling in the mouse. BMC Genomics 8: 379.A. RodriguezM. HilvoL. KytomakiRE FlemingRS Britton2007Effects of iron loading on muscle: genome-wide mRNA expression profiling in the mouse.BMC Genomics8379

28. Norton JD (2000) ID helix-loop-helix proteins in cell growth, differentiation and tumorigenesis. J Cell Sci 113 (Pt 22): 3897–3905.JD Norton2000ID helix-loop-helix proteins in cell growth, differentiation and tumorigenesis.J Cell Sci113 (Pt 22)38973905

29. Damdinsuren B, Nagano H, Kondo M, Yamamoto H, Hiraoka N, et al. (2005) Expression of Id proteins in human hepatocellular carcinoma: relevance to tumor dedifferentiation. Int J Oncol 26: 319–327.B. DamdinsurenH. NaganoM. KondoH. YamamotoN. Hiraoka2005Expression of Id proteins in human hepatocellular carcinoma: relevance to tumor dedifferentiation.Int J Oncol26319327

30. Kautz L, Meynard D, Monnier A, Darnaud V, Bouvet R, et al. (2008) Iron regulates phosphorylation of Smad1/5/8 and gene expression of Bmp6, Smad7, Id1, and Atoh8 in the mouse liver. Blood 112: 1503–1509.L. KautzD. MeynardA. MonnierV. DarnaudR. Bouvet2008Iron regulates phosphorylation of Smad1/5/8 and gene expression of Bmp6, Smad7, Id1, and Atoh8 in the mouse liver.Blood11215031509

31. Andriopoulos B Jr, Corradini E, Xia Y, Faasse SA, Chen S, et al. (2009) BMP6 is a key endogenous regulator of hepcidin expression and iron metabolism. Nat Genet 41: 482–487.B. Andriopoulos JrE. CorradiniY. XiaSA FaasseS. Chen2009BMP6 is a key endogenous regulator of hepcidin expression and iron metabolism.Nat Genet41482487

32. Brown KE, Broadhurst KA, Mathahs MM, Weydert J (2007) Differential expression of stress-inducible proteins in chronic hepatic iron overload. Toxicol Appl Pharmacol 223: 180–186.KE BrownKA BroadhurstMM MathahsJ. Weydert2007Differential expression of stress-inducible proteins in chronic hepatic iron overload.Toxicol Appl Pharmacol223180186

33. Shayeghi M, Latunde-Dada GO, Oakhill JS, Laftah AH, Takeuchi K, et al. (2005) Identification of an intestinal heme transporter. Cell 122: 789–801.M. ShayeghiGO Latunde-DadaJS OakhillAH LaftahK. Takeuchi2005Identification of an intestinal heme transporter.Cell122789801

34. Qiu A, Jansen M, Sakaris A, Min SH, Chattopadhyay S, et al. (2006) Identification of an intestinal folate transporter and the molecular basis for hereditary folate malabsorption. Cell 127: 917–928.A. QiuM. JansenA. SakarisSH MinS. Chattopadhyay2006Identification of an intestinal folate transporter and the molecular basis for hereditary folate malabsorption.Cell127917928

35. Fleming RE, Sly WS (2002) Mechanisms of iron accumulation in hereditary hemochromatosis. Annu Rev Physiol 64: 663–680.RE FlemingWS Sly2002Mechanisms of iron accumulation in hereditary hemochromatosis.Annu Rev Physiol64663680

36. Trinder D, Olynyk JK, Sly WS, Morgan EH (2002) Iron uptake from plasma transferrin by the duodenum is impaired in the Hfe knockout mouse. Proc Natl Acad Sci U S A 99: 5622–5626.D. TrinderJK OlynykWS SlyEH Morgan2002Iron uptake from plasma transferrin by the duodenum is impaired in the Hfe knockout mouse.Proc Natl Acad Sci U S A9956225626

37. Zoller H, Pietrangelo A, Vogel W, Weiss G (1999) Duodenal metal-transporter (DMT-1, NRAMP-2) expression in patients with hereditary haemochromatosis. Lancet 353: 2120–2123.H. ZollerA. PietrangeloW. VogelG. Weiss1999Duodenal metal-transporter (DMT-1, NRAMP-2) expression in patients with hereditary haemochromatosis.Lancet35321202123

38. Fleming RE, Migas MC, Zhou X, Jiang J, Britton RS, et al. (1999) Mechanism of increased iron absorption in murine model of hereditary hemochromatosis: increased duodenal expression of the iron transporter DMT1. Proc Natl Acad Sci U S A 96: 3143–3148.RE FlemingMC MigasX. ZhouJ. JiangRS Britton1999Mechanism of increased iron absorption in murine model of hereditary hemochromatosis: increased duodenal expression of the iron transporter DMT1.Proc Natl Acad Sci U S A9631433148

39. Owen D, Kuhn LC (1987) Noncoding 3′ sequences of the transferrin receptor gene are required for mRNA regulation by iron. Embo J 6: 1287–1293.D. OwenLC Kuhn1987Noncoding 3′ sequences of the transferrin receptor gene are required for mRNA regulation by iron.Embo J612871293

40. Dupic F, Fruchon S, Bensaid M, Loreal O, Brissot P, et al. (2002) Duodenal mRNA expression of iron related genes in response to iron loading and iron deficiency in four strains of mice. Gut 51: 648–653.F. DupicS. FruchonM. BensaidO. LorealP. Brissot2002Duodenal mRNA expression of iron related genes in response to iron loading and iron deficiency in four strains of mice.Gut51648653

41. Edgar R, Domrachev M, Lash AE (2002) Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res 30: 207–210.R. EdgarM. DomrachevAE Lash2002Gene Expression Omnibus: NCBI gene expression and hybridization array data repository.Nucleic Acids Res30207210

42. GEO. Gene Expression Omnibus. GEO. Gene Expression Omnibus.http://www.ncbi.nlm.nih.gov/geo/. http://www.ncbi.nlm.nih.gov/geo/.

43. The DAVID Bioinformatic Resources. The DAVID Bioinformatic Resources.http://david.abcc.ncifcrf.gov. http://david.abcc.ncifcrf.gov.

44. Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, et al. (2003) DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol 4: P3.G. Dennis JrBT ShermanDA HosackJ. YangW. Gao2003DAVID: Database for Annotation, Visualization, and Integrated Discovery.Genome Biol4P3

45. Primer3. Primer3.http://frodo.wi.mit.edu/. http://frodo.wi.mit.edu/.

46. BLAST. Basic Local Alignment and Search Tool. BLAST. Basic Local Alignment and Search Tool.http://blast.ncbi.nlm.nih.gov/Blast.cgi. http://blast.ncbi.nlm.nih.gov/Blast.cgi.

47. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3: RESEARCH0034.J. VandesompeleK. De PreterF. PattynB. PoppeN. Van Roy2002Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes.Genome Biol3RESEARCH0034