Gene array identification of Ipf1/Pdx1-/-regulated genes in pancreatic progenitor cells
© Svensson et al; licensee BioMed Central Ltd. 2007
Received: 15 March 2007
Accepted: 23 November 2007
Published: 23 November 2007
The homeodomain transcription factor IPF1/PDX1 exerts a dual role in the pancreas; Ipf1/Pdx1 global null mutants fail to develop a pancreas whereas conditional inactivation of Ipf1/Pdx1 in β-cells leads to impaired β-cell function and diabetes. Although several putative target genes have been linked to the β-cell function of Ipf1/Pdx1, relatively little is known with respect to genes regulated by IPF1/PDX1 in early pancreatic progenitor cells.
Microarray analyses identified a total of 111 genes that were differentially expressed in e10.5 pancreatic buds of Ipf1/Pdx1-/- embryos. The expression of one of these, Spondin 1, which encodes an extracellular matrix protein, has not previously been described in the pancreas. Quantitative real-time RT-PCR analyses and immunohistochemical analyses also revealed that the expression of FgfR2IIIb, that encodes the receptor for FGF10, was down-regulated in Ipf1/Pdx1-/- pancreatic progenitor cells.
This microarray analysis has identified a number of candidate genes that are differentially expressed in Ipf1/Pdx1-/- pancreatic buds. Several of the differentially expressed genes were known to be important for pancreatic progenitor cell proliferation and differentiation whereas others have not previously been associated with pancreatic development.
The pancreas is an endodermally derived organ that forms from a ventral and a dorsal evagination of the foregut epithelium. These two evaginations, the dorsal and ventral pancreatic buds, subsequently grow, branch and differentiate into distinct pancreatic cell types . The homeodomain transcription factor Insulin Promoter Factor 1/Pancreatic and Duodenal homeobox 1 (IPF1/PDX1) is one of the earliest markers of the developing pancreas. IPF1/PDX1 is expressed already at ~10 somites stage at the regions of the dorsal and ventral gut endoderm from which the pancreatic buds evaginate . IPF1/PDX1 expression remains high in pancreatic epithelial cells until ~e10.5 after which it is down-regulated  and remains low in proliferating pancreatic epithelial cells. Strong IPF1/PDX1 expression reappears in the differentiating β-cells as they emerge at ~e13  and high level of IPF1/PDX1 expression is maintained in adult β-cells where IPF1/PDX1 controls the expression of several key β-cell genes, including the insulin gene, thereby ensuring normal β-cell function and glucose homeostasis [4, 5].
Loss of Ipf1/Pdx1 gene function in mice and humans results in pancreatic agenesis demonstrating a key role for the Ipf1/Pdx1 gene in pancreatic development [6–8]. Ipf1/Pdx1 is, however, not required for the initiation of the pancreatic program and the initial stages of pancreas development, i.e. the formation of the pancreatic buds, still occurs in Ipf1/Pdx1-/- mice [7, 9]. Although the pancreatic program is initiated in Ipf1/Pdx1 deficient embryos, the subsequent growth of the embryonic pancreas is arrested, resulting in pancreas agenesis [6, 7, 9]. Recombination experiments between pancreatic epithelium and pancreatic mesenchyme have demonstrated that the pancreatic developmental defect observed in Ipf1/Pdx1-/- embryos is confined to the epithelial cells . Thus, pancreatic mesenchyme isolated from Ipf1/Pdx1-/- e10.5 dorsal pancreatic buds could support the growth of wt e10.5 dorsal pancreatic epithelium whereas the reverse combination failed to grow . These data provide evidence for a cell-autonomous role for Ipf1/Pdx1 in early pancreatic progenitor cells. To date, no direct or indirect Ipf1/Pdx1 downstream genes have, however, been identified that can explain the pancreatic phenotype observed in Ipf1/Pdx1-/- embryos.
To identify Ipf1/Pdx1 target genes in early pancreatic progenitor cells and to begin to unravel the molecular mechanisms that lead to the attenuation of pancreatic growth in Ipf1/Pdx1-/- mice we performed microarray analyses on cDNA prepared from Ipf1/Pdx1-/- e10.5 buds and stage matched littermate wildtype controls. We have identified genes that are differentially expressed in Ipf1/Pdx1-/- pancreatic buds and a subset of these was chosen for further expression analysis by quantitative real-time (qRT) RT-PCR, in situ hybridization and immunohistochemistry. In agreement with the pancreatic developmental defect observed in Ipf1/Pdx1-/- embryos, several of the differentially expressed genes identified in this study encode factors linked to pancreatic progenitor cell proliferation and differentiation.
Gene expression changes in Ipf1/Pdx1-/-pancreatic buds
Top ranked 73 down-regulated genes in Ipf1/Pdx1-/- vs. Ipf1/Pdx1 +/+ e10.5 dorsal pancreatic buds
Genbank acc. No.
Extracellular matrix proteins
secreted phosphoprotein 1
-9.8 U *
spondin 1, (f-spondin) extracellular matrix protein
-3.8 K **
solute carrier family 38, member 5
solute carrier family 2 (facilitated glucose transporter), member 3
ATP-binding cassette, sub-family G (WHITE), member 1
carbohydrate sulfotransferase 2
v-erb-b2 erythroblastic leukemia viral oncogene homolog 3 (avian)
-3.2 U *
microsomal glutathione S-transferase 1
-3.2 K *
Lines homolog 2 (Drosophila)
basic leucine zipper nuclear factor
paired box gene 6
-5.8 K *
Nkx transcription factor related, locus 1 (Drosophila)
-4.2 U ***
nuclear factor, erythroid derived 2,-like 1
nuclear respiratory factor 1
tumor-associated calcium signal transducer 1
schwannomin interacting protein 1
glucocorticoid modulatory element binding protein 2
high mobility group nucleosomal binding domain 3
polo-like kinase 2 (Drosophila)
mitogen activated protein kinase 13
protein tyrosine phosphatase, receptor type, G
integral membrane protein 2B
RIKEN cDNA 1810057C19 gene
-6.1 U *
cytidine 5'-triphosphate synthase 2
nuclear prelamin A recognition factor-like
RIKEN cDNA 9130213B05 gene
retinol binding protein 4, plasma
gulonolactone (L-) oxidase
serine (or cysteine) proteinase inhibitor, clade B, member 1a
ESTs, transcribed locus
tumor protein D52
periostin, osteoblast specific factor
RIKEN cDNA 1110003E01 gene
transmembrane protein 34
methionine aminopeptidase-like 1
THUMP domain containing 3
WAP four-disulfide core domain 2
keratin complex 2, basic, gene 7
vitamin K epoxide reductase complex, subunit 1
Mid1 interacting protein 1 (gastrulation specific G12-like (zebrafish))
proprotein convertase subtilisin/kexin type 2
aldo-keto reductase family 1, member C13
serine protease inhibitor, Kunitz type 2
SPARC related modular calcium binding 1
RIKEN cDNA 2810457I06 gene
RIKEN cDNA 4921517J23 gene
RIKEN cDNA 4833432M17 gene
CDP-diacylglycerol synthase (phosphatidate cytidylyltransferase) 2
hyperparathyroidism 2 homolog (human)
ELAV (embryonic lethal, abnormal vision, Drosophila)-like 1
ESTs, transcribed locus
H19 fetal liver mRNA
expressed sequence AI427122
RIKEN cDNA 1700022L09 gene
proviral integration site 2
Keratin complex 2, basic, gene 8
Top ranked 38 up-regulated genes in Ipf1/Pdx1-/- vs. Ipf1/Pdx1+/+ e10.5 dorsal pancreatic buds
Genbank acc. No
secreted frizzled-related sequence protein 2
lin-9 homolog (C. elegans)
secreted frizzled-related sequence protein 1
cell division cycle associated 5
E26 avian leukemia oncogene 1, 5' domain
RIKEN cDNA 2310057D15 gene
cDNA sequence BC002059
microtubule-associated protein, RP/EB family, member 1
neuron specific gene family member 1
RIKEN cDNA 6330403K07 gene
deiodinase, iodothyronine type III, antisense
expressed sequence AI449753
tripartite motif protein 6
RIKEN cDNA C230004F18 gene
RIKEN cDNA 2700007P21 gene
chemokine-like factor super family 3
plastin 3 (T-isoform)
DNA-damage inducible transcript 3
metastasis associated 3
mesoderm specific transcript
platelet derived growth factor receptor, alpha polypeptide
cofactor required for Sp1 transcriptional activation, subunit 2
RIKEN cDNA D130019J16 gene
vacuolar protein sorting 24 (yeast)
RIKEN cDNA 9930012K11 gene
Thyroid hormone receptor associated protein 3
methionine-tRNA synthetase 2 (mitochondrial)
chorionic somatomammotropin hormone 2
RIKEN cDNA 0710005M24 gene
heat shock 27 kDa protein 8
transmembrane protein 15
DNA segment, Chr 6, ERATO Doi 253, expressed
RIKEN cDNA 2010301N04 gene
neurofilament, light polypeptide
Altered expression of genes associated with cell proliferation and differentiation
Given the pancreatic growth arrest phenotype of Ipf1/Pdx1-/- mice we selected a total of 10 differentially expressed genes identified by the microarray (9 down-regulated and 1 up-regulated) that encode factors linked to growth and differentiation for further expression analyses and validation by quantitative real-time RT-PCR. The selected genes include genes previously associated with pancreatic progenitor cell proliferation and differentiation, including NK6 transcription factor related, locus 1 (Nkx6.1), v-erb-b2 erythroblastic leukemia viral oncogene homolog 3 (ErbB3), Secreted phosphoprotein 1 (Spp1) and Paired box gene 6 (Pax6), as well as genes encoding factors that previously have not been linked to pancreas development, such as Spondin 1 (Spon1), Microsomal glutathione S-transferase 1 (Mgst1), Reticulon 1 (Rtn1), Basic leucine zipper nuclear factor 1 (Blzf1), Schwannomin interacting protein 1 (Schip1), and Secreted frizzled related protein1 (Sfrp1).
In contrast, three genes, Blzf1, a gene up-regulated by retinoids [32, 33], Schip1, which is involved in cell cycle regulation [34, 35], and Sfrp1, which encode a secreted antagonist to Wnt signaling , which all were identified as differentially expressed in the microarray analysis could not be validated when analyzed by quantitative real-time RT-PCR and were excluded from further analyses (data not shown). Taken together, these data show that the expression of several genes implicated in growth, differentiation and adhesion, some of which have previously been associated with pancreatic development, is impaired in Ipf1/Pdx1-/- pancreatic progenitor cells.
Ptf1a/p48 and FgfR2IIIb expression is perturbed in Ipf1/Pdx1-/-pancreatic buds
The expression of two genes, pancreas specific transcription factor 1a/p48 (Ptf1a/p48) and fibroblast growth factor receptor 2IIIb (FgfR2IIIb), that previously have been implicated to have a role in pancreatic development [37–39] but which were not spotted on the microarray chips were also analyzed by quantitative real-time RT-PCR. FgfR2IIIb is expressed in the developing pancreas and in FgfR2IIIb null mutant mice pancreatic growth is reduced and the branching and morphogenesis of the pancreatic ductal epithelium is impaired . Quantitative real-time RT-PCR analysis revealed that FgfR2IIIb expression was reduced by ~75% in e10.5 Ipf1/Pdx1-/- dorsal pancreatic buds compared to that of stage matched Ipf1/Pdx1+/+ dorsal pancreatic buds (Fig. 1). Immunohistochemical analyses confirmed the reduced expression of FGFR2IIIb in e10.5 pancreatic progenitor cells of Ipf1/Pdx1-/- deficient mice (Fig. 2e, f).
PTF1a/p48 has been implicated in commitment and proliferation of pancreatic progenitor cells and Ptf1a/p48-/- null mutant mice fail to develop a pancreas [37, 38]. Quantitative real-time RT-PCR analyses showed that Ptf1a/p48 expression was reduced by ~91% in Ipf1/Pdx1-/- e10.5 dorsal pancreatic buds as compared to wild-type littermates (Fig. 1). In situ hybridization analyses confirmed the reduced expression of Ptf1a/p48 in dorsal pancreatic buds of Ipf1/Pdx1-/- embryos compared to wildtype littermates (Fig. 2g, h). These results suggest that in early pancreatic progenitor cells both FgfR2IIIb and Ptf1a/p48 are direct or indirect Ipf1/Pdx1 downstream target genes.
Spon1expression in pancreatic progenitor cells overlaps temporally with the phase of proliferation
Spon1 encodes the extracellular matrix protein F-spondin, also known as vascular smooth muscle cell growth-promoting factor. F-spondin stimulates the proliferation and growth of vascular SMC during ovarian growth and development . In addition, F-spondin has been reported to promote neurite outgrowth from spinal chord, hippocampal and sensory neurons [23, 25, 27] as well as the differentiation of neural progenitor cells . Spon1 has not previously been described to be expressed in the developing pancreas. Thus, we next determined the temporal and spatial expression of Spon1 in the developing pancreas by in situ hybridization.
In Ipf1/Pdx1 deficient embryos, pancreatic development becomes arrested around e10 [6, 7, 9] but the molecular mechanism(s) underlying this early developmental arrest has remained largely unknown. To begin to identify genes regulated by Ipf1/Pdx1 in early pancreatic progenitor cells that could help to explain the developmental arrest observed in Ipf1/Pdx1 -/- embryos we embarked on a gene array profiling and quantitative real-time RT-PCR approach. Given the early developmental arrest and hence limiting amount of tissue, i.e. dorsal buds from e10.5 Ipf1/Pdx1 +/+ and Ipf1/Pdx1 -/- litter pairs, that could be retrieved for further cDNA synthesis we first performed a pilot study using the 15 K cDNA chip that was hybridized twice with cDNA prepared from dorsal buds of a matched Ipf1/Pdx1 +/+ and Ipf1/Pdx1 -/- e10.5 litter pair.
This pilot study was followed by two hybridizations of the 15 K cDNA chip with cDNA prepared from another pair of matched Ipf1/Pdx1 +/+ and Ipf1/Pdx1 -/- e10.5 dorsal buds. The 21 K cDNA chip was originally hybridized with cDNA from 4 pairs of Ipf1/Pdx1 +/+ and Ipf1/Pdx1 -/- e10.5 dorsal buds in a standard manner. However due to technical problems in the hybridization step one of the arrays had to be excluded from the experiment resulting in a 3 arrays from the 21 K experiment. In summary a total of 4 independent pairs of Ipf1/Pdx1 +/+ and Ipf1/Pdx1 -/- e10.5 dorsal buds were analyzed on four 15 K arrays and three 21 K arrays.
Hence, there were some limitations in our microarray study. We had a total of four biological replicates per condition: the 15 K experiment used only 2 biological replicates per condition and the 21 K experiment used 3 biological replicates per condition. This limitation was due to challenges in retrieving sufficient tissue material for subsequent cDNA synthesis. We analyzed separately the data from the two arrays, so the results of these analyses are based on a small number of biological replicates. Furthermore, the experimental design was unbalanced in various ways (different types of technical replication in the assays using the 15 K arrays, unbalanced dye assignments, etc.). Thus, the list of differentially expressed genes identified solely by our microarray experiments might be sub-optimal in terms of a bona fide candidate list. However, it provided us with sufficient information for selecting potentially interesting Ipf1/Pdx1 regulated genes that we followed-up with quantitative real-time RT-PCR, in situ hybridization and immunohistochemistry.
Here we show that the expression of FgfR2IIIb is reduced in Ipf1/Pdx1-/- e10.5 pancreatic buds, thus linking IPF1/PDX1 to FGF signaling in pancreatic progenitor cells. FGF signaling has been shown to be important both for normal pancreas development and maintenance of pancreatic function in adults [5, 40, 41]. Mice with a targeted disruption of Fgf10 show an early block in pancreas development that is reminiscent of that observed in Ipf1/Pdx1 null mutant mice. Alike the Ipf1/Pdx1-/- mice the pancreatic program is initiated in Fgf10-/- mice, i.e. the evagination of the pancreatic buds occur, but the subsequent proliferation and branching of the pancreatic epithelium is perturbed . Moreover, mice lacking the FGF10 receptor, FgfR2IIIb, display a pancreatic growth defect with impaired branching and morphogenesis of the ductal epithelium. . Moreover, forced expression of Fgf10 in the developing pancreatic epithelium results in pancreatic hyperplasia, adding further evidence for a role for FGF10-FGFR2IIIb signaling in pancreatic progenitor cell proliferation [40, 41].
Taken together, the almost identical pancreatic phenotype observed in mice lacking either Ipf1/Pdx1, the receptor FgfR2IIIb, or Fgf10, and the reduced pancreatic expression of FgfR2IIIb in Ipf1/Pdx1-/- embryos provide evidence for a role for IPF1/PDX1 in regulating FGF signaling in the developing pancreatic epithelium. Thus, Ipf1/Pdx1 appears to regulate pancreatic growth by ensuring high level expression of FGFR2IIIb in pancreatic progenitor cells. Interestingly, Ipf1/Pdx1 has also been shown to regulate the expression of another FGF-receptor, FGFR1c, in adult β-cells thereby ensuring normal β-cell function and glucose homeostasis .
Epidermal growth factor (EGF) signaling has also been linked to pancreatic growth and development [1, 19, 43]. Several EGF receptors, including ErbB3, have been shown to be expressed in the developing pancreatic epithelium . In addition, ErbB3 has been described to be expressed in the pancreatic mesenchyme [18, 19]. Mice lacking a functional ErbB1/Egfr gene have mildly perturbed pancreatic growth  and mice lacking a functional ErbB3 gene present with pancreatic hypoplasia [18, 44]. Additional support for a role of the EGF family in pancreas development originates from in vitro studies in which EGF was shown to stimulate the growth of e13.5 rat pancreatic epithelium that had been depleted of mesenchyme . Thus, in addition to the reduced expression of FgfR2IIIb in progenitor cells of Ipf1/Pdx1 deficient mice the reduced expression of ErbB3 is likely to contribute to the perturbed pancreatic development that these mice present with.
Ptf1a/p48 expression is severely reduced in Ipf1/Pdx1 null mutant mice, providing evidence that Ptf1a/p48, directly or indirectly, is a downstream target gene of Ipf1/Pdx1. In support of a role for Ipf1/Pdx1 in regulating Ptf1a/p48 expression, inactivation of the Ipf1/Pdx1 gene in the pancreatic epithelium at progressively later stages also results in reduced expression of Ptf1a/p48 . In contrast, Kawaguchi et al. have shown that transgenic expression of Ipf1/Pdx1 under the control of the Ptf1a/p48 promoter in an Ipf1/Pdx1-/- null background partially rescued the pancreatic developmental arrest , suggesting that Ptf1a/p48 expression is independent of Ipf1/Pdx1. However, the Ptf1a-Ipf1/Pdx1 transgenic rescue approach leaves open the possibility that the partially rescued phenotype is the result of leaky transgene expression due to positional integration effects and/or transgene copy number. Moreover, since the hsp68 minimal promoter was coupled to the Ptf1a/p48 promoter in the transgenic rescue approach, the transgene expression may differ from that of the endogenous Ptf1a/p48 promoter.
Alike Ipf1/Pdx1, Nkx6.1 is expressed in pancreatic progenitors from embryonic day 9.5 and thereafter the expression of Nkx6.1 also becomes restricted to mature β-cells. Nkx6.1 has previously been identified as a downstream target gene of IPF1/PDX1; β-cell specific inactivation of the Ipf1/Pdx1 gene resulted in a loss of Nkx6.1 expression in β-cells [4, 11] and in agreement with previous findings [4, 11] we show that Nkx6.1 expression is abolished also in pancreatic progenitor cells of Ipf1/Pdx1 null mutant mice. These data provide strong evidence that Ipf1/Pdx1 act upstream of Nkx6.1 in both pancreatic progenitor cells and adult β-cells.
F-spondin, encoded by Spon1, is a 90 kDa extracellular matrix protein that has been reported to play an important role in proliferation and growth of vascular SMC during ovarian growth and development . F-spondin has also been showed to promote neurite outgrowth from spinal chord, hippocampal, and sensory neurons [23, 25, 27]. Other reports describe a role for F-spondin in the development of the nervous system where it is predominately expressed during early stages of development but down-regulated during neural differentiation . F-spondin has also been suggested to promote differentiation of neural precursor cells to nerve cell-like cells . Interestingly, a F-spondin related protein, R-spondin, has been shown to function as a growth factor; injection of R-spondin in mice resulted in an increased proliferation of intestinal crypt epithelial cells in both the small intestine and colon . Similar to the transient expression of Spon1 in the developing nervous system, Spon1 is transiently expressed, from e9.5 to e15.5 in progenitor cells of the developing pancreas. Thus, Spon1 expression in the developing pancreas overlaps with the phase of pancreatic progenitor cell proliferation. Whether F-spondin is involved in regulation of proliferation and/or differentiation of pancreatic progenitor cells will have to be determined by functional analyses in vitro and/or in vivo.
Taken together, our data show that Ipf1/Pdx1, directly or indirectly, regulates the expression of genes, such as Ptf1a/p48, FGFR2IIIb, and ErbB3, implicated in pancreatic growth and differentiation, providing evidence that the perturbed expressions of these factors provide a mechanistic explanation for the pancreatic developmental arrest that Ipf1/Pdx1-/- mice present with. In addition we show that Spon1, which encodes an extracellular matrix protein, is expressed in the developing pancreas and that Spon1 expression is down-regulated in Ipf1/Pdx1-/- embryos.
The animal studies were approved by the Institutional Animal Care and Use Committee of Umeå University and conducted in accordance with the Guidelines for the care and use of Laboratory Animals. The generation of Ipf1/Pdx1-/- mice has been previously described . Ipf1/Pdx1 wildtype and null mutant embryos were obtained from our local breeding stock.
Isolation of embryos and tissue
Ipf1/Pdx1+/- mice were mated and embryonic day 0.5 (e0.5) was determined by the presence of morning vaginal plug. Pregnant females were sacrificed at e10.5 and embryos were transferred to ice-cold Lebovitz L-15 medium (Gibco). Dorsal pancreatic buds were dissected out in L-15 medium and snap frozen in liquid nitrogen for later RNA isolation. Embryos sacrificed for immunohistochemistry and in situ hybridization were fixated in 4% paraformaldehyde (SIGMA) in 0.1 M PBS for 2 hours and overnight respectively at +4°C followed by dehydration in 20% sucrose overnight at 4°C and mounted in Tissue-Tek (Sakura) and stored at -80°C.
Genomic DNA was extracted from embryonic tails by proteinase K (Boehringer) digestion followed by 2-propanol precipitation. Genotyping was performed by PCR using two internal Ipf1/Pdx1 primers; 5'-AGAAGCTGGCCACTAGCCTCT-3' and 5'-CTGTGGGCAACAAGGGAGTT-3', and two internal Neo primers; 5'-AGAAGCTATTCGGCTATGA-3' and 5'-TTTCCACCATGATATTCG-3'. The PCR reactions result in a 150 bp and a 548 bp product respectively if target sequence is present.
Microarray, hybridization and scanning
The cDNA microarrays used in this study were of two types. The first contain approximately 15,000 clones printed in duplicate, in two identical fields, yielding 30,000 probes. The set of probes represents approximately 12,000 unique genes obtained through large-scale, in-house EST sequencing of three cDNA libraries. The first library originates from a neural tissue stem cell compartment (lateral ventricular wall), the second from neurospheres (neural stem cells cultured in vitro) and the third from a hematopoietic stem cell line expressing the Lhx2 gene . The second cDNA array used in this study contains 20,600 clones derived from two different clone sets: a 15,000 mouse cDNA set from National Institute of Aging (NIH) and a 5,400 cDNA clone set obtained from Research Genetics. Microarray slides were hybridized at 42°C in a water bath for 16–20 hours (15 K array) or at 37°C in a Genetac hybridization station (Genomic Solutions) for 15 hours (21 K array). The 15 K arrays were scanned using the G2565BA DNA microarray scanner (Agilent Technologies) and images analyzed in the GenePix5.0 software (Axon Instruments, CA, USA). The 21 K arrays were scanned using the Scanarray 4000 XL microarray scanner (Perkin Elmer) at three settings (laser power/PMT): 80/80, 90/90 and 100/100 (in that order), where the numbers were percentages of maximum values. Images were analyzed using the ScanarrayExpress software (Perkin Elmer).
RNA isolation, amplification and labeling
Total RNA was isolated from individual pancreatic buds using Nucleospin RNA II kit (Clontech) according to manufacturer's guidelines. cDNA was generated and amplified using SuperSMART™ PCR cDNA synthesis kit (Clontech). 2 μg of each cDNA preparation was labeled with the fluorescent dyes Cy3 or Cy5 using BD Atlas™ SMART™ Fluorescent Probe Amplification kit (BD Biosciences Clontech) according to manufacturer's recommendations. Samples hybridized to the 15 K cDNA microarrays were labeled and analyzed as follows: In the first set, the labeled samples (one wildtype bud cDNA labeled with Cy3, one Ipf1/Pdx1 null mutant pancreatic bud cDNA labeled with Cy5) were split in two, and analyzed as duplicates on two separate microarrays. In the second set, the samples (one wildtype bud cDNA and one Ipf1/Pdx1 null mutant pancreatic bud cDNA) were split before labeling, and analyzed as a dye-swapped replicate on two separate microarrays.
Samples hybridized to the 21 K cDNA microarrays were labeled and analyzed as follows: two wildtype and one Ipf1/Pdx1 null mutant pancreatic bud were individually labeled with Cy3 and two Ipf1/Pdx1 null mutants and one wildtype pancreatic bud were individually labeled with Cy5. Samples were pooled pair wise and hybridized to three 21 K cDNA microarrays.
One wt cDNA sample and one Ipf1/Pdx1 null mutant cDNA sample hybridized to the 21 K array were identical to second set of samples hybridized to the 15 K array. Due to the limited amount of material, the rest of the samples were hybridized to one array each (15 K or 21 K array).
In total, cDNA generated from pancreatic buds from 4 individual Ipf1/Pdx1 null mutant and 4 wildtype littermates respectively were individually labeled, pooled pair wise and hybridized to four 15 K cDNA microarrays and to three 21 K cDNA microarrays.
Microarray data processing and statistical analysis
The 15 K and 21 K experiments were analyzed separately as described below. Given some selection criteria, discussed below, 111 genes were considered differentially expressed. 73 of these genes were down-regulated (Table 1) and 38 up-regulated (Table 2) in Ipf1/Pdx1 null mutants. Additional information and all array data are reported to ArrayExpress  with the following accession numbers. [ArrayExpress: E-TABM-280] and [ArrayExpress: E-MEXP-1047].
The 21 K microarray analyses were carried out using ScanArrayExpress 2.1 (PerkinElmer) and the in-house S-Plus library UmeaSAMED 1.39. The data analysis included: combined image analysis extracting median intensity values, complete filtration, local background correction, RLS 100, removed C-spots, MA-print-tip normalization with span = 0.5, censoring with minimum value 8, and the B-test, see  for detailed descriptions of the procedures. Genes with top-ranked B-statistics (top 100) and at least two-fold regulated (i.e. the absolute value of the average log-ratio 1), were selected for further studies.
The 15 K microarray analyses were carried out using R (1.7), Bioconductor, LIMMA (1.3), and the aroma package [50–52]. Several filtering steps, normalization and statistical analyses were done, as detailed in ArrayExpress. The foreground intensities without background subtraction were extracted using the median intensity values. A feature was filtered if it was flagged as "Not Found" in GenePix, if both channels were saturated, the percentage of foreground pixels above the local background plus two of its standard deviations were below 70 for both channels, or the signal-to-noise ratio for both channels were below 3. Data was normalized using MA-print-tip normalization . The intensities for duplicated features were averaged after normalization. The empirical Bayes moderated t-test implemented in the LIMMA package (i.e. the B-test) was used  to identify differentially expressed genes . The four 15 K assays were analyzed in two separate sub-experiments, each including the two assays derived from the same biological replicates. Genes that were at least two-fold up-regulated (or down-regulated) in both sub-experiments and had B-values larger than 10 in sub-experiment 1 and 2 in sub-experiment 2 were selected for further studies.
Quantitative real- time RT-PCR
Gene expression levels were measured using the ABI PRISM 7000 Sequence Detection System (Applied Biosystems) according to the manufactures recommendations. Expression of the β-2-microglobulin (b2M and TBP) was used to normalize expression levels. Primer sequences were the following: B2M: 5'-GCTATCCAGAAAACCCCTCAAA-3', 5'-CTGTGTTACGTAGCAGTTCAGTATGTTC-3'; Blzf1: 5'-ACAGCCCAGCTTTCTGAACAG-3', 5'-ACTCATCTGCCATTACCCTGCTT'3'; ErbB3: 5'-CTCTCCCCACCAGGGTTAGAG-3', 5'-TGACGAAAGGGTCCCTTCC-3'; Fgfr2IIIb: 5'-CAGTAAATACGGGCCTGATGG-3', 5'-GAACAGAGCCAGCACTTCTGC-3'; Foxa2: 5'-GAGCCATCCGACTGGAGCA-3', 5'-GGAATGAGCCCGTCGCTAG-3'; Glu: 5'-CCGCCGTGCCCAAGA-3', 5'-CATCATGACGTTTGGCAATGTT-3'; Ipf1: 5'-TAGGACTCTTTCCTGGGACCAA-3', 5'-AATAAAAAGGGTACAAACTTGAGCGT-3'; Mgst1: 5'-CGCATTCCAGAGGATAACCAA-3', 5'-CAAAGCCAGCGCAGTCTTCT-3'; Nkx6.1: 5'-TCAGGTCAAGGTCTGGTTCCA-3', 5'GGTCTCCGAGTCCTGCTT-3'; P48: 5'-TCATCTGCCATCGAGGCAC-3', 5'-ACCATAATCCGGGTCACTGG-3'; Pax6: 5'-GCAAACAACCTGCCTATGCAA-3', 5'-GGCAGCATGCACGAGTATGA-3'; Reticulon: 5'-CCTGGTCCAGGACCTGGTG-3', 5'-AGCCACATGAGGACTGCAAAT-3'; Schip1: 5'-AAAAAGTCTCCCGTCGCTGA-3', 5'-TTCGCTTATGTGAGGCATGTGT-3'; Sfrp1: 5'-CCTCCATGCGACAACGAGTT-3', 5'-TTGATTTTCATCCTCAGTGCAAAC-3'; Spondin1: 5'-GAAGCAGCTACCGAGTGACACTT-3', 5'-CTTTGAGGGCAATTAACGTGAAA-3'; Spp1: 5'-TCAAAGTCTAGGAGTTTCCAGGTTTC-3', 5'-TGTGGCATCAGGATACTGTTCATC-3'; TBP: 5'-GAATTGTACCGCAGCTTCAAAA-3', 5'-AGTGCAATGGTCTTTAGGTCAAGTT-3'.
Immunohistochemistry was performed as described . Primary antibodies used were: rabbit anti-IPF1/PDX1 , rabbit Nkx6.1 , rabbit anti-FgfR2 (Research Diagnostics, Inc) and glucagon anti-guinea pig (Electra-Box Diagnostika Ab). Secondary antibodies used were: Cy3 anti-rabbit (Jacksson) and Alexa 488 anti-guinea pig (Molecular Probes).
In situ hybridization was performed as described . Probes used were: Ptf1a/p48 (provided by O. Hagenbüchle), Sox9 (covering region 1365–2305 bp of image clone 5320371) and Spon1 (covering region 1237–4035 bp of image clone 4221758).
We thank members of our laboratory for technical instructions, suggestions and helpful discussions. This work was supported by grants from Wallenberg Consortium North, Kempe foundation, the Swedish Research Council and the Juvenile Diabetes Research Foundation (to H.E).
- Edlund H: Transcribing pancreas. Diabetes. 1998, 47: 1817-23. 10.2337/diabetes.47.12.1817.View ArticlePubMedGoogle Scholar
- Jonsson J, Ahlgren U, Edlund T, Edlund H: IPF1, a homeodomain protein with a dual function in pancreas development. Int J Dev Biol. 1995, 39: 789-98.PubMedGoogle Scholar
- Ohlsson H, Karlsson K, Edlund T: IPF1, a homeodomain-containing transactivator of the insulin gene. Embo J. 1993, 12: 4251-9.PubMed CentralPubMedGoogle Scholar
- Ahlgren U, Jonsson J, Jonsson L, Simu K, Edlund H: beta-cell-specific inactivation of the mouse Ipf1/Pdx1 gene results in loss of the beta-cell phenotype and maturity onset diabetes. Genes Dev. 1998, 12: 1763-8.PubMed CentralView ArticlePubMedGoogle Scholar
- Hart AW, Baeza N, Apelqvist A, Edlund H: Attenuation of FGF signalling in mouse beta-cells leads to diabetes. Nature. 2000, 408: 864-8. 10.1038/35048589.View ArticlePubMedGoogle Scholar
- Jonsson J, Carlsson L, Edlund T, Edlund H: Insulin-promoter-factor 1 is required for pancreas development in mice. Nature. 1994, 371: 606-9. 10.1038/371606a0.View ArticlePubMedGoogle Scholar
- Offield MF, Jetton TL, Labosky PA, Ray M, Stein RW, Magnuson MA, Hogan BL, Wright CV: PDX-1 is required for pancreatic outgrowth and differentiation of the rostral duodenum. Development. 1996, 122: 983-95.PubMedGoogle Scholar
- Stoffers DA, Zinkin NT, Stanojevic V, Clarke WL, Habener JF: Pancreatic agenesis attributable to a single nucleotide deletion in the human IPF1 gene coding sequence. Nat Genet. 1997, 15: 106-10. 10.1038/ng0197-106.View ArticlePubMedGoogle Scholar
- Ahlgren U, Jonsson J, Edlund H: The morphogenesis of the pancreatic mesenchyme is uncoupled from that of the pancreatic epithelium in IPF1/PDX1-deficient mice. Development. 1996, 122: 1409-16.PubMedGoogle Scholar
- Williams C, Wirta V, Meletis K, Wikstrom L, Carlsson L, Frisen J, Lundeberg J: Catalog of gene expression in adult neural stem cells and their in vivo microenvironment. Exp Cell Res. 2006, 312: 1798-812. 10.1016/j.yexcr.2006.02.012.View ArticlePubMedGoogle Scholar
- Pedersen JK, Nelson SB, Jorgensen MC, Henseleit KD, Fujitani Y, Wright CV, Sander M, Serup P: Endodermal expression of Nkx6 genes depends differentially on Pdx1. Dev Biol. 2005, 288: 487-501. 10.1016/j.ydbio.2005.10.001.View ArticlePubMedGoogle Scholar
- Jensen J, Serup P, Karlsen C, Nielsen TF, Madsen OD: mRNA profiling of rat islet tumors reveals nkx 6.1 as a beta-cell-specific homeodomain transcription factor. J Biol Chem. 1996, 271: 18749-58. 10.1074/jbc.271.31.18749.View ArticlePubMedGoogle Scholar
- Oster A, Jensen J, Edlund H, Larsson LI: Homeobox gene product Nkx 6.1 immunoreactivity in nuclei of endocrine cells of rat and mouse stomach. J Histochem Cytochem. 1998, 46: 717-21.View ArticlePubMedGoogle Scholar
- Oster A, Jensen J, Serup P, Galante P, Madsen OD, Larsson LI: Rat endocrine pancreatic development in relation to two homeobox gene products (Pdx-1 and Nkx 6.1). J Histochem Cytochem. 1998, 46: 707-15.View ArticlePubMedGoogle Scholar
- Rudnick A, Ling TY, Odagiri H, Rutter WJ, German MS: Pancreatic beta cells express a diverse set of homeobox genes. Proc Natl Acad Sci USA. 1994, 91: 12203-7. 10.1073/pnas.91.25.12203.PubMed CentralView ArticlePubMedGoogle Scholar
- Sander M, Neubuser A, Kalamaras J, Ee HC, Martin GR, German MS: Genetic analysis reveals that PAX6 is required for normal transcription of pancreatic hormone genes and islet development. Genes Dev. 1997, 11: 1662-73. 10.1101/gad.11.13.1662.View ArticlePubMedGoogle Scholar
- St-Onge L, Sosa-Pineda B, Chowdhury K, Mansouri A, Gruss P: Pax6 is required for differentiation of glucagon-producing alpha-cells in mouse pancreas. Nature. 1997, 387: 406-9. 10.1038/387406a0.View ArticlePubMedGoogle Scholar
- Erickson SL, O'Shea KS, Ghaboosi N, Loverro L, Frantz G, Bauer M, Lu LH, Moore MW: ErbB3 is required for normal cerebellar and cardiac development: a comparison with ErbB2-and heregulin-deficient mice. Development. 1997, 124: 4999-5011.PubMedGoogle Scholar
- Kritzik MR, Krahl T, Good A, Gu D, Lai C, Fox H, Sarvetnick N: Expression of ErbB receptors during pancreatic islet development and regrowth. J Endocrinol. 2000, 165: 67-77. 10.1677/joe.0.1650067.View ArticlePubMedGoogle Scholar
- Kilic G, Wang J, Sosa-Pineda B: Osteopontin is a novel marker of pancreatic ductal tissues and of undifferentiated pancreatic precursors in mice. Dev Dyn. 2006, 235: 1659-67. 10.1002/dvdy.20729.PubMed CentralView ArticlePubMedGoogle Scholar
- Shijubo N, Uede T, Kon S, Nagata M, Abe S: Vascular endothelial growth factor and osteopontin in tumor biology. Crit Rev Oncog. 2000, 11: 135-46.PubMedGoogle Scholar
- Sulzbacher I, Birner P, Trieb K, Lang S, Chott A: Expression of osteopontin and vascular endothelial growth factor in benign and malignant bone tumors. Virchows Arch. 2002, 441: 345-9. 10.1007/s00428-002-0671-4.View ArticlePubMedGoogle Scholar
- Burstyn-Cohen T, Frumkin A, Xu YT, Scherer SS, Klar A: Accumulation of F-spondin in injured peripheral nerve promotes the outgrowth of sensory axons. J Neurosci. 1998, 18: 8875-85.PubMedGoogle Scholar
- Burstyn-Cohen T, Tzarfaty V, Frumkin A, Feinstein Y, Stoeckli E, Klar A: F-Spondin is required for accurate pathfinding of commissural axons at the floor plate. Neuron. 1999, 23: 233-46. 10.1016/S0896-6273(00)80776-X.View ArticlePubMedGoogle Scholar
- Feinstein Y, Borrell V, Garcia C, Burstyn-Cohen T, Tzarfaty V, Frumkin A, Nose A, Okamoto H, Higashijima S, Soriano E, et al: F-spondin and mindin: two structurally and functionally related genes expressed in the hippocampus that promote outgrowth of embryonic hippocampal neurons. Development. 1999, 126: 3637-48.PubMedGoogle Scholar
- Feinstein Y, Klar A: The neuronal class 2 TSR proteins F-spondin and Mindin: a small family with divergent biological activities. Int J Biochem Cell Biol. 2004, 36: 975-80. 10.1016/j.biocel.2004.01.002.View ArticlePubMedGoogle Scholar
- Klar A, Baldassare M, Jessell TM: F-spondin: a gene expressed at high levels in the floor plate encodes a secreted protein that promotes neural cell adhesion and neurite extension. Cell. 1992, 69: 95-110. 10.1016/0092-8674(92)90121-R.View ArticlePubMedGoogle Scholar
- Miyamoto K, Morishita Y, Yamazaki M, Minamino N, Kangawa K, Matsuo H, Mizutani T, Yamada K, Minegishi T: Isolation and characterization of vascular smooth muscle cell growth promoting factor from bovine ovarian follicular fluid and its cDNA cloning from bovine and human ovary. Arch Biochem Biophys. 2001, 390: 93-100. 10.1006/abbi.2001.2367.View ArticlePubMedGoogle Scholar
- Schubert D, Klar A, Park M, Dargusch R, Fischer WH: F-spondin promotes nerve precursor differentiation. J Neurochem. 2006, 96: 444-53. 10.1111/j.1471-4159.2005.03563.x.View ArticlePubMedGoogle Scholar
- Busenlehner LS, Codreanu SG, Holm PJ, Bhakat P, Hebert H, Morgenstern R, Armstrong RN: Stress sensor triggers conformational response of the integral membrane protein microsomal glutathione transferase 1. Biochemistry. 2004, 43: 11145-52. 10.1021/bi048716k.View ArticlePubMedGoogle Scholar
- Steiner P, Kulangara K, Sarria JC, Glauser L, Regazzi R, Hirling H: Reticulon 1-C/neuroendocrine-specific protein-C interacts with SNARE proteins. J Neurochem. 2004, 89: 569-80. 10.1111/j.1471-4159.2004.02345.x.View ArticlePubMedGoogle Scholar
- Tong JH, Duprez E, Lanotte M: JEM-1, a novel nuclear co-factor: localisation and functional interaction with AP-1. Leukemia. 1999, 13: 1982-92. 10.1038/sj/leu/2401560.View ArticlePubMedGoogle Scholar
- Tong JH, Fant X, Duprez E, Benoit G, Uphoff CC, Drexler HG, Pla JC, Lofvenberg E, Lanotte M: Expression patterns of the JEM-1 gene in normal and tumor cells: ubiquity contrasting with a faint, but retinoid-induced, mRNA expression in promyelocytic NB4 cells. Leukemia. 1998, 12: 1733-40. 10.1038/sj.leu.2401173.View ArticlePubMedGoogle Scholar
- Goutebroze L, Brault E, Muchardt C, Camonis J, Thomas G: Cloning and characterization of SCHIP-1, a novel protein interacting specifically with spliced isoforms and naturally occurring mutant NF2 proteins. Mol Cell Biol. 2000, 20: 1699-712. 10.1128/MCB.20.5.1699-1712.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Muranen T, Gronholm M, Renkema GH, Carpen O: Cell cycle-dependent nucleocytoplasmic shuttling of the neurofibromatosis 2 tumour suppressor merlin. Oncogene. 2005, 24: 1150-8. 10.1038/sj.onc.1208283.View ArticlePubMedGoogle Scholar
- Kawano Y, Kypta R: Secreted antagonists of the Wnt signalling pathway. J Cell Sci. 2003, 116: 2627-34. 10.1242/jcs.00623.View ArticlePubMedGoogle Scholar
- Kawaguchi Y, Cooper B, Gannon M, Ray M, MacDonald RJ, Wright CV: The role of the transcriptional regulator Ptf1a in converting intestinal to pancreatic progenitors. Nat Genet. 2002, 32: 128-34. 10.1038/ng959.View ArticlePubMedGoogle Scholar
- Krapp A, Knofler M, Ledermann B, Burki K, Berney C, Zoerkler N, Hagenbuchle O, Wellauer PK: The bHLH protein PTF1-p48 is essential for the formation of the exocrine and the correct spatial organization of the endocrine pancreas. Genes Dev. 1998, 12: 3752-63.PubMed CentralView ArticlePubMedGoogle Scholar
- Pulkkinen MA, Spencer-Dene B, Dickson C, Otonkoski T: The IIIb isoform of fibroblast growth factor receptor 2 is required for proper growth and branching of pancreatic ductal epithelium but not for differentiation of exocrine or endocrine cells. Mech Dev. 2003, 120: 167-75. 10.1016/S0925-4773(02)00440-9.View ArticlePubMedGoogle Scholar
- Hart A, Papadopoulou S, Edlund H: Fgf10 maintains notch activation, stimulates proliferation, and blocks differentiation of pancreatic epithelial cells. Dev Dyn. 2003, 228: 185-93. 10.1002/dvdy.10368.View ArticlePubMedGoogle Scholar
- Norgaard GA, Jensen JN, Jensen J: FGF10 signaling maintains the pancreatic progenitor cell state revealing a novel role of Notch in organ development. Dev Biol. 2003, 264: 323-38. 10.1016/j.ydbio.2003.08.013.View ArticlePubMedGoogle Scholar
- Bhushan A, Itoh N, Kato S, Thiery JP, Czernichow P, Bellusci S, Scharfmann R: Fgf10 is essential for maintaining the proliferative capacity of epithelial progenitor cells during early pancreatic organogenesis. Development. 2001, 128: 5109-17.PubMedGoogle Scholar
- Edlund H: Pancreas: how to get there from the gut?. Curr Opin Cell Biol. 1999, 11: 663-8. 10.1016/S0955-0674(99)00033-2.View ArticlePubMedGoogle Scholar
- Miettinen PJ, Huotari M, Koivisto T, Ustinov J, Palgi J, Rasilainen S, Lehtonen E, Keski-Oja J, Otonkoski T: Impaired migration and delayed differentiation of pancreatic islet cells in mice lacking EGF-receptors. Development. 2000, 127: 2617-27.PubMedGoogle Scholar
- Cras-Meneur C, Elghazi L, Czernichow P, Scharfmann R: Epidermal growth factor increases undifferentiated pancreatic embryonic cells in vitro: a balance between proliferation and differentiation. Diabetes. 2001, 50: 1571-9. 10.2337/diabetes.50.7.1571.View ArticlePubMedGoogle Scholar
- Hale MA, Kagami H, Shi L, Holland AM, Elsasser HP, Hammer RE, MacDonald RJ: The homeodomain protein PDX1 is required at mid-pancreatic development for the formation of the exocrine pancreas. Dev Biol. 2005, 286: 225-37. 10.1016/j.ydbio.2005.07.026.View ArticlePubMedGoogle Scholar
- Kim KA, Kakitani M, Zhao J, Oshima T, Tang T, Binnerts M, Liu Y, Boyle B, Park E, Emtage P, et al: Mitogenic influence of human R-spondin1 on the intestinal epithelium. Science. 2005, 309: 1256-9. 10.1126/science.1112521.View ArticlePubMedGoogle Scholar
- ArrayExpress. [http://www.ebi.ac.uk/arrayexpress/#ae-main]
- Ryden P, Andersson H, Landfors M, Naslund L, Hartmanova B, Noppa L, Sjostedt A: Evaluation of microarray data normalization procedures using spike-in experiments. BMC Bioinformatics. 2006, 7: 300-10.1186/1471-2105-7-300.PubMed CentralView ArticlePubMedGoogle Scholar
- Bengtsson H: aroma – An R Object-oriented Microarray Analysis environment. Mathematical Sciences. 2004, 18:Google Scholar
- Gentleman RIaR: R: A Language for Data Analysis and Graphics. Journal of Computational and Graphical Statistics. 1996, 5: 299-314. 10.2307/1390807.Google Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-PubMedGoogle Scholar
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30: e15-10.1093/nar/30.4.e15.PubMed CentralView ArticlePubMedGoogle Scholar
- Lonnstedt I, Britton T: Hierarchical Bayes models for cDNA microarray gene expression. Biostatistics. 2005, 6: 279-91. 10.1093/biostatistics/kxi009.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.