<title>Background</title> <p>An understanding of the molecular and cellular regulation of pancreas development is emerging [##REF##12094230##1##, ####REF##10842072##2##, ##REF##10615124##3##, ##REF##9836511##4##, ##REF##10600706##5####10600706##5##]. Expression of the transcription factor <italic>Pdx1 </italic>is essential for pancreas development and is initiated at Theiler stage (TS) 13 in the region of gut endoderm destined to become the pancreas [##REF##9334273##6##, ####REF##11577200##7##, ##REF##14699589##8####14699589##8##]. At TS14, the foregut endoderm evaginates to form the dorsal pancreas bud [##REF##9334273##6##,##REF##9620856##9##,##REF##10725233##10##]. The ventral bud appears somewhat later (TS17-TS20). Expression of <italic>Ptf1a</italic>, another critical regulatory factor, is detected at this stage and is essential for the generation of both exocrine and endocrine cell types [##REF##14736742##11##, ####REF##15031114##12##, ##REF##12185368##13####12185368##13##]. The 'secondary transition', from TS20 to TS22, marks the differentiation of pancreas precursors into endocrine and exocrine cell types. The notch signaling pathway plays a critical role in this process through the lateral inhibition of neighboring cells [##REF##10842072##2##,##REF##10615124##3##,##REF##10476967##14##,##REF##10868931##15##]. Subsequently, endocrine progenitors express the essential basic helix-loop-helix transcription factor <italic>Neurog3 </italic>[##REF##11973276##16##, ####REF##10677506##17##, ##REF##10825208##18####10825208##18##]. In response to <italic>Neurog3 </italic>expression, endocrine precursor cells express a number of transcriptional regulators, including <italic>B2</italic>/<italic>NeuroD</italic>, <italic>Pax6</italic>, <italic>Isl1</italic>, <italic>Nkx2-2</italic>, <italic>Nkx6-1</italic>, and others, that play roles in the differentiation and maturation of the various endocrine cells types [##REF##14699589##8##,##REF##11776473##19##]. By TS24 the majority of cell fates are established and remodeling of the pancreas begins with initially scattered endocrine cells formed at duct tips starting to migrate. At TS26, isletogenesis occurs as endocrine cells fuse and form recognizable 'islets', while acinar cells gain their mature ultrastructure. Pancreas development continues postnatally, with β-cells gaining the ability to sense glucose levels and respond with pulsatile insulin release.</p> <p>Analysis of the transcriptomes of precursor cells present at different stages of pancreas development is expected to further facilitate a definition of the genetic cascades essential for endocrine and exocrine differentiation. Towards this end a number of microarray expression profiling studies have been carried out on the developing pancreas [##REF##12789652##20##, ####REF##12592642##21##, ##REF##14660441##22##, ##REF##12636919##23##, ##REF##12829622##24##, ##REF##12086925##25##, ##REF##18071024##26####18071024##26##]. Serial analysis of gene expression (SAGE), like microarrays, provides a quantitative analysis of gene expression profiles. A major advantage of SAGE, however, is that the data are digital, making it easily shared amongst investigators and compared across different experiments and tissues.</p> <p>In this study we describe the construction and analyses of ten SAGE libraries from TS17 to TS26 (embryonic days 10.5-18.5) mouse pancreases as well as from adult islets and ducts. <italic>Pdx1 </italic>enhanced green fluorescent protein (EGFP) and <italic>Neurog3 </italic>EGFP reporter strains [##REF##14660441##22##] were employed to allow fluorescence activated cell sorting (FACS) purification of pancreatic and endocrine progenitor cell populations, respectively, at early stages of mouse pancreas development. To our knowledge we are the first group to generate SAGE libraries from embryonic pancreas tissues. In sum, we sequenced over 2 million SAGE tags representing over 200,000 tag types, providing a truly comprehensive view of pancreas development. To validate our results, we assessed the temporal expression profiles of 44 genes by quantitative real-time PCR (qRT-PCR) and categorized the TS22 pancreas staining patterns of 601 genes in the <italic>GenePaint </italic>database [##REF##17953485##27##,##REF##14681479##28##], providing insight into the expression profiles of hundreds of transcripts previously not described in the pancreas. We then used the libraries to construct a network of predicted transcription factor interactions describing β-cell development, and validated selected linkages in this network using chromatin immunoprecipitation followed by qPCR (ChIP-qPCR) to detect enrichment of binding sites. Taken together, we anticipate these data will act as a framework for future studies on the regulatory networks driving pancreas development and function.</p>
<title>Materials and methods</title> <title>Mouse maintenance and islet isolation</title> <p>All mice were bred and maintained at the British Columbia Cancer Research Centre animal facility according to the guidelines of the Canadian Council on Animal Care. All protocols were approved by the University of British Columbia Animal Care Committee. Mice were housed in micro-isolator units, provided with Purina mouse food and autoclaved water <italic>ad libidum</italic>, and were maintained at 20°C ± 2°C under a light/dark cycle (light, 05:00-19:00; dark, 19:00-05:00). Males were mated overnight with up to three females and females were checked for plugs before 9:30 the following morning. Plugged mice were considered to be 0.5 days post coitum.</p> <title>Generation of SAGE libraries</title> <p>Embryos, obtained from timed pregnant female <italic>Pdx1 </italic>EGFP and <italic>Ngn3 </italic>EGFP mice mated with C57Bl/6J mice, were staged using Theiler criteria. Embryos were collected in ice cold phosphate buffered saline (PBS) solution. The developing pancreas was isolated by micro-dissection while the embryo was submerged in PBS. For TS17-TS22 <italic>Pdx1 </italic>EGFP and <italic>Ngn3 </italic>EGFP FACS sorted libraries, a single cell suspension was generated by Trypsin/EDTA digestion. The sample was prepared for sorting by staining with propidium iodide and passing through a 40 μm filter with a final resuspension in FACS buffer (PBS containing 2% fetal bovine serum and 2 mM EDTA). EGFP positive or negative cells were sorted and collected as appropriate into a microcentrifuge tube containing Trizol<sup>® </sup>using BD FACSVantage instruments (operated by the Terry Fox Laboratory at the BC Cancer Research Centre). Between three and eight litters of mice were used to generate each library, depending on the stage. Islets and ducts from at least 10 separate mice were purified from 8-10-week-old ICR males by collagenase digestion and gradient centrifugation as previously described [##REF##15618358##55##] and placed in Trizol<sup>®</sup>. Total RNA was extracted using Trizol<sup>® </sup>and RNA quality was assessed using a Bioanalyser (Agilent, Santa Clara, CA, USA). For the TS22 and TS26 whole pancreas as well as the islet and duct libraries roughly 10 μg of total RNA was used to construct each SAGE library. For the FACS sorted libraries a SAGElite approach [##REF##17135571##56##] was used with a minimum of 20 ng of total RNA. All libraries were constructed as previously described [##REF##16352711##36##].</p> <title>SAGE data analysis</title> <p>SAGE data were analyzed using DiscoverySpace4 [##REF##17210078##57##]. All SAGE libraries were generated as part of the Mouse Atlas of Gene Expression project [##REF##16352711##36##] and filtered for sequence quality so that each tag had a 95% or greater probability of being correct, using the PHRED score quality assessment software [##REF##9521922##58##]. Tag to gene mapping was performed using the mouse Refseq, MGC, and Ensembl databases using the DiscoverySpace program. Tags were considered sense position matches if they mapped in the sense orientation to the gene and antisense matches if they mapped in the opposite orientation. A tag was considered unambiguous if it matched a single sense position gene in all of the databases, and ambiguous if it mapped to multiple genes in a sense position regardless of the mapping position.</p> <p>The specificity of tags was determined by first obtaining the counts for the tags in 205 different Mouse Atlas libraries. From this the mean of the tag counts in all the libraries (M<sub>a</sub>) was determined and compared to tag count in the library of interest (C<sub>i</sub>) to obtain the mean ratio (M<sub>r</sub>). The total number of libraries the tag was found in, or library count, was next determined (L<sub>a</sub>) as was the total counts of the tag in the library under analysis (C<sub>i</sub>). The specificity (S) was then calculated as:</p> <p></p> <p>Thus, tags with a high mean ratio that appear in relatively few libraries and are expressed more abundantly will have the highest specificities.</p> <p>In short, the data were first normalized and FOM analysis for the <italic>K</italic>-means algorithm with Euclidean distance was performed on normalized data using 5 iterations and a maximum of 20 clusters. For this we used the multiexperiment viewer from the Institute for Genomic Research [##REF##12613259##59##]. From this a cluster number of 14 was chosen and used in a PoissonC-based clustering strategy specifically designed for SAGE data [##REF##15239836##39##]. GO term comparisons using EASE [##REF##14519205##60##] and KEGG pathway enrichments were calculated using Webgestalt [##REF##15980575##61##], comparing against the Webgestalt mouse database using a hypergeometric test to calculate <italic>p</italic>-values and a minimum of two genes required in a gene set. Heatmaps were generated using relative tag abundances using the multiexperiment viewer from the Institute for Genomic Research [##REF##12613259##59##].</p> <title><italic>GenePaint </italic>analyses</title> <p>Images of <italic>in situ </italic>hybridization staining patterns for whole embryo sagittal sections were obtained from the <italic>GenePaint </italic>website [##REF##14681479##28##,##UREF##6##62##] for 923 genes identified in our clusters. Additional data file 4 lists the gene names and <italic>GenePaint </italic>probe IDs for all of these genes. When appropriate, higher magnification images of the area of the embryo containing the pancreas were obtained and the pancreas outlined. The brightness and contrast of some of the images were altered using Photoshop to better assess the staining pattern. Genes were classified as showing trunk, tip, epithelial, mesenchmymal or vasculature staining, based in part on the similarity of their staining pattern to staining patterns for <italic>Ins2</italic>, <italic>Gcg</italic>, <italic>Neurog3</italic>, <italic>Ela3</italic>, <italic>Pnliprp1</italic>, <italic>Foxa2</italic>, <italic>Sox9</italic>, <italic>Prrx1</italic>, or <italic>Ets1</italic>.</p> <title>Quantitative real-time PCR</title> <p>Probes for <italic>Ins1</italic>, <italic>Ins2</italic>, <italic>Pdx1</italic>, <italic>Neurog3</italic>, <italic>Ptf1a</italic>, <italic>Nkx2-2</italic>, <italic>Amy</italic>, and <italic>β-actin </italic>were purchased from Applied Biosystems (Foster City, CA, USA). All other primers were designed using Primer3 [##UREF##7##63##] and spanned introns where possible (Additional data file 5). Amplicons were between 80 and 120 bp for efficient amplification. Primer efficiencies were determined using a dilution series of TS24 cDNA. Only primer pairs with an efficiency greater than 0.8 were used in subsequent analyses. An ABI 7500 real-time PCR system (Applied Biosystems) and SYBR<sup>® </sup>Green supermix (Applied Biosystems) or Universal PCR Master Mix (Applied Biosystems) was used for all reactions. Triplicate cDNAs were obtained by reverse transcription of 1 μg of total RNA from newly isolated pancreas tissue for each reverse transcription. Generated cDNA (10 ng) was used in each reaction with all reactions done in duplicate. Samples were normalized to GAPDH, and the fold increase compared to TS23 as calculated using 2<sup>-ΔΔCt</sup> 
[##REF##11328886##64##]. To determine significance, the non-parametric Mann Whitney test was used to compare the ΔCt values of the reference and the sample using Prism4 (GraphPad Software, La Jolla, CA, USA). SOTA analysis was performed using the multiexperiment viewer from the Institute for Genomic Research [##REF##12613259##59##] using Euclidean distance, 4 max. cycles and default parameters.</p> <title>Network generation</title> <p>To generate networks, transcription factors with evidence of expression in the relevant cell types in addition to being in our SAGE data (in at least two libraries to eliminate possible artifacts), from the literature or in the EPConDB/T1Dbase [##UREF##5##48##] were identified. The binding site preferences of the identified transcription factors were then obtained from Transfac or the literature (Additional data file 6). These sites were then used to search the CisRED database [##REF##16381958##49##] to identify conserved motifs in the proximal promoters of factors that were co-expressed or expressed in daughter cells. This allowed us to create predicted transcription factor interaction lists, with literature reported interactions layered on these predictions [##REF##14699589##8##,##REF##15298336##41##]. In cases where two factors with the same binding site were co-expressed - for example, Neurog3 and Neurod1 - both are shown to be linked to the target site. The resulting network was then visualized using Biotapestry [##REF##15907831##65##].</p> <title>Chromatin immunoprecipitation</title> <p>Min6 cells were homogenized in 1% formaldehyde and incubated at room temperature for 10 minutes, prior to the addition of 0.125 M glycine followed by incubation at room temperature for 5 minutes. Cells were pelleted and resuspended in 5 volumes ChIP cell lysis buffer (10 mM Tris-Cl, pH 8.0, 10 mM NaCl, 3 mM MgCl<sub>2</sub>, 0.5% NP-40) containing protease inhibitor cocktail (Roche, Laval, QC, Canada). Cells were then rehomogenized, incubated on ice for 5 minutes and repelleted. The pellet was resuspended in 1.5 volumes ChIP nuclear lysis buffer (1% SDS, 5 mM EDTA, 50 mM Tris-Cl, pH 8.1) supplemented with protease inhibitor cocktail tablets (Roche). Cells were passed through a 26.5 gauge needle prior to sonication in a water-ice bath (Sonicator 3000, Misonix, Farmingdale, NY, USA) for 6 cycles of 30 s on, 40 s off, and 60 μg of chromatin was pre-cleared with 250 μl Protein G beads (100 μl, Active Motif, Carlsbad, CA, USA), Protein Inhibitor Cocktail (0.5 μl, Active Motif) and supplemented with 7.5 μl ChIP dilution buffer (0.1% SDS, 10% Triton X-100, 1.67 M NaCl, 167 M Tris-Cl, pH 8.1). Beads were precipitated and 3 μg of rabbit α-Foxa2 (M-20; Santa Cruz, Santa Cruz, CA, USA), rabbit α-Pdx1 (Millipore, Billerica, MA, USA) antibody or normal rabbit IgG (Santa Cruz) was added to supernatants. Fresh Protein G beads were also blocked with 1 mg/mg bovine serum albumin and 0.1 mg/mg salmon sperm DNA in ChIP dilution buffer. Following overnight incubation, rocking at 4°C, the samples were incubated with the blocked beads for 4 h, rocking at 4°C. The beads were then precipitated and washed in low salt buffer (0.1% SDS, 1% Triton X-100, 2 mM EDTA, 20 mM Tris-Cl, pH 8.1, 150 mM NaCl), high salt buffer (low salt buffer with 500 mM NaCl), lithium chloride buffer (0.25 M LiCl, 1% NP-40, 1% deoxycholate, 1 mM EDTA, 10 mM Tris-Cl, pH 8.1) and 2× TE buffer (10 mM Trix-HCl, ph 8.0, 1 mM EDTA). Beads were resuspended in 125 μl elution buffer (1% SDS, 0.1 M NHCO<sub>3</sub>) and rotated for 1 h at room temperature. NaCl (0.192 M) was added to reverse crosslinks and samples were incubated overnight at 65°C. Samples were then incubated with Proteinase K (Invitrogen, Carlsbad, CA, USA) and RNaseA (Sigma, Oakville, ON, Canada) for 1 h at 50°C. DNA was purified by two rounds of phenol-chloroform extraction and ethanol precipitation and resuspended in 50 μl dH<sub>2</sub>O. To detect enrichment, we used qRT-PCR, as described above, using the primers listed in Additional data file 5.</p> <title>Accession numbers</title> <p>The GenBank [##UREF##2##33##] accession numbers for the genes discussed are: <italic>Pdx1 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_008814">NM_008814</ext-link>), <italic>Ptf1a </italic>(<ext-link ext-link-type="gen" xlink:href="NM_018809">NM_018809</ext-link>), <italic>Neurog3 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_009719">NM_009719</ext-link>), <italic>Neurod1 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_010894">NM_010894</ext-link>), <italic>Pax6 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_013627">NM_013627</ext-link>), <italic>Isl1 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_021459">NM_021459</ext-link>), <italic>Nkx2-2 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_010919">NM_010919</ext-link>), <italic>Nkx6-1 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_144955">NM_144955</ext-link>), <italic>Myt3 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_173868">NM_173868</ext-link>), <italic>Foxa2 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_010446">NM_010446</ext-link>), <italic>Mafa </italic>(<ext-link ext-link-type="gen" xlink:href="NM_194350">NM_194350</ext-link>), <italic>Pax4 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_011038">NM_011038</ext-link>), <italic>Bhlhb8 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_010800">NM_010800</ext-link>), <italic>Gcg </italic>(<ext-link ext-link-type="gen" xlink:href="NM_008100">NM_008100</ext-link>), <italic>Iapp </italic>(<ext-link ext-link-type="gen" xlink:href="NM_010491">NM_010491</ext-link>), <italic>Ins1 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_008386">NM_008386</ext-link>), <italic>Ins2 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_008387">NM_008387</ext-link>), <italic>Pnlip </italic>(<ext-link ext-link-type="gen" xlink:href="NM_026925">NM_026925</ext-link>), <italic>Sdha </italic>(<ext-link ext-link-type="gen" xlink:href="NM_023281">NM_023281</ext-link>), <italic>HbS1L </italic>(<ext-link ext-link-type="gen" xlink:href="NM_001042593">NM_001042593</ext-link>), <italic>B2m </italic>(<ext-link ext-link-type="gen" xlink:href="NM_009735">NM_009735</ext-link>), <italic>Plunc </italic>(<ext-link ext-link-type="gen" xlink:href="NM_011126">NM_011126</ext-link>), <italic>Cldn13 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_020504">NM_020504</ext-link>), <italic>Pomc </italic>(<ext-link ext-link-type="gen" xlink:href="NM_008895">NM_008895</ext-link>), <italic>Prm2 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_008933">NM_008933</ext-link>), <italic>Alb </italic>(<ext-link ext-link-type="gen" xlink:href="NM_009654">NM_009654</ext-link>), <italic>Amy1 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_007446">NM_007446</ext-link>), <italic>Gast </italic>(<ext-link ext-link-type="gen" xlink:href="NM_010257">NM_010257</ext-link>), <italic>Tcf12 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_011544">NM_011544</ext-link>), <italic>Zfp326 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_018759">NM_018759</ext-link>), <italic>Meox1 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_010791">NM_010791</ext-link>), <italic>Zfp446 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_175558">NM_175558</ext-link>), <italic>Rnf6 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_028774">NM_028774</ext-link>), <italic>Son </italic>(<ext-link ext-link-type="gen" xlink:href="ENSMUSG00000022961">ENSMUSG00000022961</ext-link>), <italic>Nr1d1 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_145434">NM_145434</ext-link>), <italic>Bcl6b </italic>(<ext-link ext-link-type="gen" xlink:href="NM_007528">NM_007528</ext-link>), <italic>Arx </italic>(<ext-link ext-link-type="gen" xlink:href="ENSMUSG00000035277">ENSMUSG00000035277</ext-link>), <italic>Mfng </italic>(<ext-link ext-link-type="gen" xlink:href="NM_008595">NM_008595</ext-link>), <italic>Hes1 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_008235">NM_008235</ext-link>), <italic>Notch1 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_008714">NM_008714</ext-link>), <italic>β-actin </italic>(<ext-link ext-link-type="gen" xlink:href="NM_007393">NM_007393</ext-link>), <italic>Foxa2 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_010446">NM_010446</ext-link>), <italic>Myt1 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_008665">NM_008665</ext-link>), <italic>Myt1L </italic>(<ext-link ext-link-type="gen" xlink:href="NM_001093775">NM_001093775</ext-link>), <italic>Ela3 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_026419">NM_026419</ext-link>), <italic>Pnliprp1 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_018874">NM_018874</ext-link>), <italic>Sox9 </italic>(<ext-link ext-link-type="gen" xlink:href="NM_011448">NM_011448</ext-link>), <italic>Prrx1 </italic>(<ext-link ext-link-type="gen" xlink:href="BC092372">BC092372</ext-link>), <italic>Ets1 </italic>(<ext-link ext-link-type="gen" xlink:href="BC114351">BC114351</ext-link>).</p>
<title>Results</title> <title>Validating the biological significance of the pancreas SAGE libraries</title> <p>In order to gain further insights into pancreas development and to provide a complementary analysis to available microarray data, we generated ten SAGE libraries from the mouse pancreas tissues by sequencing a total of 2,266,558 tags (Table ##TAB##0##1##). These libraries are publicly available at the Mouse Atlas [##UREF##0##29##] or CGAP SAGE websites [##UREF##1##30##] and can be analyzed using tools available through these sites. A total of 208,412 different tag types were detected in these libraries after stringent quality selection.</p> <p>To confirm that the libraries accurately represent the cell types intended (Table ##TAB##0##1##), we assessed the distribution of tags in the libraries for genes with well-characterized expression profiles in pancreas development. Figure ##FIG##0##1## shows that transcription factors expressed in pancreas progenitor epithelial cells, such as <italic>Pdx1 </italic>and <italic>Nkx2-2</italic>, can be found in our TS17-TS19 <italic>Pdx1 </italic>EGFP+ libraries. Tags for these genes were also found frequently in the <italic>Neurog3 </italic>EGFP+ libraries. This is in agreement with the known expression of these factors. For example, <italic>Pdx1 </italic>is expressed in essentially all pancreas epithelial cells prior to the secondary transition while its expression after the secondary transition is abundant only in β-cells and β-cell precursors [##REF##14699589##8##]. Prior to the secondary transition <italic>Neurog3 </italic>expression is quite low; however, at the start of the secondary transition its expression increases dramatically [##REF##15604203##31##] and is subsequently lost quickly thereafter. This is precisely what we see in our data - low <italic>Neurog3 </italic>levels in the <italic>Pdx1 </italic>EGFP+ libraries, high expression in the <italic>Neurog3 </italic>EGFP+ libraries and diminishing expression in the TS22 and TS26 whole pancreas libraries, with no expression in the <italic>Neurog3 </italic>EGFP- or the adult islet or duct libraries. <italic>Neurod1</italic>, <italic>Isl1</italic>, <italic>Pax6 </italic>and <italic>Pax4 </italic>expression occurs subsequent to <italic>Neurog3</italic>, but unlike <italic>Neurog3 </italic>their expression is maintained in endocrine cell types [##REF##14699589##8##]. In our data it is clear that the expression of all of these genes is most abundant in the <italic>Neurog3 </italic>EGFP+ libraries, or the islet library, as would be predicted. <italic>Ptf1a </italic>and <italic>Bhlhb8 </italic>(<italic>Mist1</italic>) are two transcription factors known to drive exocrine cell development. <italic>Ptf1a </italic>was found only in the TS22 whole pancreas library, and while low levels of <italic>Bhlhb8 </italic>were noted in the TS22 <italic>Neurog3 </italic>EGFP+ library, much higher levels were found in the duct cell library. Markers of mature exocrine cells showed peak expression in the TS26 whole pancreas or adult duct libraries, with moderate expression also in the islet library, suggesting a low level of exocrine cell contamination in this library. <italic>Glucagon </italic>expression peaked in the <italic>Neurog3 </italic>EGFP+ libraries, which is not surprising as Glucagon-positive cells are relatively abundant at these time points compared to in the adult islet. <italic>Iapp</italic>, <italic>Ins1 </italic>and <italic>Ins2 </italic>were all most abundant in the islet library, as was expected. The expression of these genes was also noted in the duct library, suggesting some level of islet cell contamination in this library. In sum, the expression profiles of these selected markers in our data match predictions based on their known expression profiles, indicating that our libraries accurately reflect the cell types and stages intended.</p> <title>Count and specificity thresholds</title> <p>In SAGE data, tags with very low counts (especially those present as singletons) are enriched in error tags and their counts have little statistical power. It is useful, therefore, to use a minimum tag count threshold. To determine what count level to threshold our data at, in order to maximize the comprehensiveness of the data, while at the same time ensuring a high level of reliability, we assessed how different tag count thresholds affected the number of tags that mapped to known pancreas expressed transcripts or expressed sequence tags (ESTs). This analysis revealed that a threshold of a minimum raw count of 4 provided a good compromise between the number of tags kept and the percentage of tags that mapped to known pancreas expressed transcripts or ESTs (Additional data file 1). Additionally, in comparisons using Audic and Claverie statistics [##REF##9331369##32##], tags with a count of 4 were statistically different from 0 at <italic>p </italic>≤ 0.05. From the 10 pancreas SAGE libraries, 16,233 tags met this threshold. Of these, 70% (11,656) mapped to known transcripts using the Refseq [##UREF##2##33##], Ensembl transcript [##UREF##3##34##], and MGC [##UREF##4##35##] databases with 85% (9,918) of these mapped unambiguously in the sense direction. These 9,918 unambiguously mapped sense tags represented 7,911 different genes, suggesting that many of the genes have alternative transcript termination sites, although this remains to be validated. A further 11% (1,817) of tags mapped only to the genome and possibly represent novel genes, leaving 17% (2,760) of tags we were unable to map. These results suggest the comprehensive nature of our data and suggest that our libraries are potentially a rich source of novel pancreas expressed transcripts.</p> <p>It was of particular interest to us to identify genes with pancreas specific functions, rather than genes with ubiquitous roles in development or cellular function. We wanted, therefore, to institute a further threshold based on the specificity of the tags to the pancreas libraries. For this, we obtained the counts for the 11,735 tags that mapped unambiguously to a specific transcript or mapped uniquely to the genome in a total of 205 different SAGE libraries [##REF##16352711##36##], including the libraries created here. Next, we calculated the specificities (S values) of each of these tags to each of the 205 libraries by dividing the ratio of the tag count in the library of interest versus its mean count in all the other libraries, multiplied by the log of its count in the library of interest, by the number of libraries the tag was found in. Tags were then ranked on their maximum specificity in any one of the pancreas libraries. Table ##TAB##1##2## lists the 25 most specific tags identified in the pancreas libraries. As expected, tags that map to markers of mature pancreas cell types (that is, <italic>Ins1</italic>, <italic>Ins2</italic>, <italic>Pnlip</italic>) were very high on the list.</p> <p>To validate that these rankings accurately reflect the level of restriction of a gene's expression pattern, we compared our results with TS22 whole embryo <italic>in situ </italic>hybridization staining patterns using the <italic>GenePaint </italic>database [##REF##17953485##27##,##REF##14681479##28##]. We did this with sets of transcripts with high (S &gt; 0.1, representing 5% of the genes), medium (0.001 &gt; S &lt; 0.1, representing 25% of the genes), and low (S &lt; 0.001, representing 70% of the genes) S values. Figure ##FIG##1##2## indicates that the calculated S values correlated extremely well with the relative restriction of the staining seen in the TS22 whole embryo sections. Genes with high S values showed staining specifically in the pancreas, genes with medium S values showed staining in the pancreas and a limited number of other tissues, and genes with low S values showed broad staining throughout the embryo. Additionally, our metric met biological expectation and genes with known pancreas specificity (<italic>Ins1 </italic>S = 27.9, <italic>Ins2 </italic>S = 62.7, <italic>Gcg </italic>S = 10.985, and so on) had very high S values, while housekeeping genes (<italic>Sdha </italic>S = 0.0006, <italic>HbS1L </italic>S = 0.0002, <italic>B2m </italic>S = 0.0005) had very low S values. Meanwhile, genes with restricted expression to other tissues either did not meet our count threshold (<italic>Plunc</italic>, <italic>Cldn13</italic>, <italic>Pomc</italic>, <italic>Prm2</italic>, and so on) [##REF##17519037##37##] or had very low S values (<italic>Alb </italic>S = 0.0007). Together, these observations provided confidence in our specificity metric and we set a threshold of a minimum S of 0.002, as this value occurs roughly at the inflection point between medium and high S values in the plot of S value versus cumulative tag types represented (Figure ##FIG##1##2##). In sum, 2,536 (approximately 20%) tags met this threshold.</p> <title>SAGE tag clustering</title> <p>We next wanted to group the tags based on their differential expression during pancreas development so as to segregate them based on their potential functional significance to the different stages and cell types represented by our libraries. First, a FOM analysis for the <italic>K</italic>-means algorithm with Euclidean distance was performed on normalized data, essentially as described [##REF##16524483##38##]. Based on these results we performed a 14-cluster analysis using the PoissonC algorithm [##REF##15239836##39##] with subsequent hand curation to finalize the clusters (Figure ##FIG##2##3## and Additional data file 2).</p> <p>A summary of the clusters (Table ##TAB##2##3##) revealed that tags for genes with similar known pancreas function cluster together. For example, genes essential to endocrine cell specification were predominantly found in cluster 5, pancreatic enzyme genes in clusters 11 and 12, and islet hormone genes in cluster 13. The clusters also showed differential enrichment for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway terms (Table ##TAB##2##3##). Of interest, the clusters also had distinctively different median specificities, with cluster 5 containing genes with the highest median S, followed by cluster 13. These two clusters are enriched in genes in the mature onset diabetes of the young KEGG pathway and contain many endocrine specific factors, and this reflects the specialized nature of these cells. Cluster 14 had the lowest median S and the flattest expression profile of the clusters. In sum, these data suggested that the clusters represented biologically distinct gene sets.</p> <title>Validation of SAGE tag clusters</title> <p>To validate the identified clusters, we first compared our data to lists of genes determined to be enriched in pancreatic progenitors, endocrine cells, or islets using Affymetrix microarray analysis of <italic>Pdx1 </italic>EGFP+ and <italic>Neurog3 </italic>EGFP+ cells and islet tissues, similar to those used here [##REF##14660441##22##]. There were 107 genes present in both genes sets and the representation of each enrichment group from the array analysis in our clusters calculated (Additional data file 3). Of the 29 genes identified as enriched in pancreatic progenitors in the microarray analysis, we identified 13 of these in clusters 1-3 or cluster 9 that show peak expression early in pancreas development. Another 11 were found in clusters 10 and 11 that show peak expression in the TS26 whole pancreas library or the duct library, stages and tissue types that were not used in the array analyses. Of 24 genes identified in the array study as enriched in endocrine cells, 19 were found in cluster 5, with 2 more in cluster 4, both of which show peak expression in the <italic>Neurog3 </italic>EGFP+ libraries here. Of the genes identified as islet enriched in the array studies, 16 of 54 were classified as such in our study; a further 20 were found in clusters 11 and 12 that have peak expression in the ducts, again a tissue not represented in the array studies; and a further 10 were found in clusters 5 or 8 that show peak expression in the <italic>Neurog3 </italic>EGFP+ libraries and islet library, respectively. Overall, the two data sets compare well and the majority of genes were identified as enriched in the same cell populations, although the differences in the tissues used in each study, specifically our inclusion of developing whole pancreas and adult duct libraries, did cause differences in some of the results.</p> <p>To further confirm that our clusters accurately group genes with similar temporal expression profiles, we analyzed the expression of 44 genes through pancreas development using qRT-PCR. Selected targets included <italic>Ins2</italic>, <italic>Nkx2-2</italic>, <italic>Pdx1</italic>, <italic>Neurog3</italic>, <italic>Amy1</italic>, and <italic>Ptf1a</italic>, which all have well established expression profiles as references. We then used a self-organizing tree algorithm (SOTA) clustering analysis to group the obtained temporal expression profiles for these genes. This allowed us to determine if groupings similar to those found in the SAGE data cluster analysis were observed. In our SOTA analysis, genes with four distinct expression profiles were identified (Figure ##FIG##3##4##): one group with peak expression in the islet sample, one with peak expression in the TS26 whole pancreas, one with peak expression from TS21-TS26, and one with peak expression in the ducts sample. All of the genes in the SOTA groups containing <italic>Ins2</italic>, <italic>Mafa</italic>, <italic>Pdx1</italic>, and <italic>Nkx2-2</italic>, which are markers of the endocrine lineage, were from clusters 1, 4, 5, and 13. Three of the six genes in the SOTA group with peak expression at TS26 were from clusters 4 and 5, although each of these showed relatively high expression in either the TS22 or TS26 whole pancreas libraries. Of the genes in the SOTA group with peak expression from TS21-TS26, one was from cluster 3, two were from cluster 5 and one was from cluster 9. Clusters 3 and 9 are enriched in mesenchymal factors (see below). Since no mesenchymal cells should be present in the islet and duct samples, it makes sense for these genes to have this expression profile. Two genes from cluster 5 were in this SOTA group, including <italic>Neurog3</italic>, which is known to be developmentally restricted in expression, and <italic>Gast</italic>, likely reflecting the relative number of Gastrin-producing cells in the different samples. Of the 11 genes in the SOTA group with peak expression in the ducts sample, 4 were from clusters 7 and 12, while the rest were found in the other clusters, although significantly excluding clusters 13 and 8. All of the genes in this group had counts in the duct library, despite being in clusters with peak expression in other libraries, although they all had, in general, low overall tag counts.</p> <title><italic>GenePaint </italic>analysis</title> <p>Taken together, the data suggested that the generated clusters represent transcript sets with distinct roles in pancreas development. To further confirm this, we assessed whether the transcripts identified in each of the SAGE tag clusters had spatial expression profiles consistent with these roles using the <italic>GenePaint </italic>database [##REF##17953485##27##,##REF##14681479##28##]. For each of the 923 genes present in our clusters and in the <italic>GenePaint </italic>database, we analyzed the <italic>in situ </italic>hybridization staining pattern in the pancreas from TS22 whole embryo sections. In sum, 601 of the genes showed informative staining, and these were categorized based on their staining patterns into one of five expression domains found in the pancreas [##REF##17609113##40##] (Figure ##FIG##4##5##). For the remaining 316 genes, either the probes did not show stain in any sections or sections with pancreas were not present in the database. Regardless, we identified 88 genes expressed in the tips of epithelial branches that at E14.5 primarily contain exocrine progenitor cell types. A further 81 genes were identified as expressed in the trunk of the epithelial branches that contains endocrine and ductal progenitor cells; 221 genes were identified as expressed throughout the epithelium; and a further 51 were found only in the mesenchyme, and 42 in the vasculature. For a full categorization of the genes see Additional data file 4. There were 124 (13%) genes identified in our SAGE data that were not detected in the pancreas at the time point assessed. The average tag count for these genes was only 6.8 while for detected genes it was 24, suggesting this is, in part, due to the low expression levels of these genes. Moreover, the number of genes not detected was highest in clusters 1, 2 and 13, which include genes that show low relative expression at TS22.</p> <p>Analysis of the representation of each of the five staining profiles in the 14 identified SAGE tag clusters (Figure ##FIG##5##6##) revealed that some clusters are far more predictive of specific expression profiles than others. For example, clusters 9 and 12 were not at all predictive of a given staining pattern. In contrast, 73% of the genes analyzed in cluster 4 showed pan-epithelial staining, while 59% of those in cluster 8 showed trunk staining and a further 24% showed pan-epithelial staining. These data suggest that several of the clusters represent genes with distinct spatial expression profiles. Significantly, these profiles are consistent with the known roles of genes within the clusters. For example, cluster 5, which is enriched in genes involved in mature onset diabetes of the young, contains tags with peak expression in the <italic>Neurog3 </italic>EGFP+ libraries and genes in this cluster predominately show pan-epithelial or trunk expression. It is apparent from these results, in combination with the median profiles of the clusters (Figure ##FIG##2##3##), that clusters 1, 2, and 4 represent genes appropriately expressed spatially and temporally to be functionally significant in pancreatic progenitor cells, while genes in clusters 4, 5, and 8 are likely functionally significant in endocrine progenitors, genes from clusters 6, 10, and 11 in exocrine progenitors, genes from cluster 3 in mesenchymal cells, and genes from clusters 8 and 13 in adult islet cells. Taken together, these analyses allowed us to identify lists of transcripts, many of which have not previously been characterized in pancreas development, that are appropriately expressed spatially and temporally to play functionally significant roles in each of the major phases of pancreas development.</p> <title>Identification of a transcriptional cascade in endocrine pancreas development</title> <p>Transcriptional regulators, in concert with signaling factors, provide the genetic instructions that drive endocrine pancreas development. Based on our analyses, genes in the endocrine lineage (that is, from pancreatic progenitor cells through endocrine progenitors to adult islet cells) can be identified based on their presence in clusters 1, 2, 4, 5, 8 and 13. We identified 58 tags for transcriptional regulators in these clusters representing 43 different factors. Eliminating factors only expressed in one library, and for which we could not find additional support for their expression, left 38 different factors for which there is good evidence of their expression in the endocrine pancreas lineage. Eight of these genes were represented by multiple tag types, with four separate tags mapping to <italic>Neurod1 </italic>and <italic>Isl1</italic>. Figure ##FIG##6##7a## shows a heatmap of the expression of these factors in endocrine cell development as detected in our SAGE data. Many of these factors have well-established roles in pancreas development, including <italic>Ngn3</italic>, <italic>Pax4</italic>, <italic>Nkx2-2</italic>, <italic>Pdx1</italic>, <italic>Isl1</italic>, and <italic>Neurod1</italic>. However, several factors with uncharacterized roles were also identified, including: <italic>Tcf12</italic>, <italic>Zfp326 </italic>and <italic>Meox1</italic>, which were most abundant in <italic>Pdx1 </italic>EGFP+ cells; <italic>Zfp446</italic>, <italic>Rnf6</italic>, and <italic>Son</italic>, which were most abundant in <italic>Neurog3 </italic>EGFP+ cells; and <italic>Nr1d1</italic>, <italic>Myt3</italic>, and <italic>Bcl6b</italic>, which were most abundant in islet cells.</p> <p>The transcriptional circuits essential to β-cell development are only beginning to be elucidated [##REF##14699589##8##,##REF##15298336##41##] and the majority of regulatory interactions driving pancreas development remain unknown. Expression profiling has been used extensively to generate networks with co-expressed genes 'linked' using the assumption that co-expression implies co-regulation [##REF##16651540##42##, ####REF##15749820##43##, ##REF##16698233##44##, ##REF##15778709##45####15778709##45##]. The promoters of co-expressed genes are often then analyzed for common sequence elements and the presence of known binding sites. One advancement on this technique is to include analysis of sequence conservation under the elements to better predict 'real' sites [##REF##15084257##46##]. Such methods have been successfully applied to mammalian cells to identify networks active in B-cells and during macrophage activation [##REF##16698233##44##,##REF##15778709##45##,##REF##17534062##47##]. Therefore, to similarly gain further insight into β-cell development, we utilized our transcription factor cascade as it represents sets of co-expressed transcription factors.</p> <p>For this, we first eliminated factors for which we could not find direct evidence of expression in the cell types in the β-cell lineage, that is, pancreas progenitors, endocrine progenitors, and mature β-cells, using the <italic>GenePaint in situ </italic>data, literature, or the EPConDB/T1Dbase [##UREF##5##48##], or as mentioned above were found in only one library. We then searched the proximal promoters of each of the remaining transcription factors, using the CisRED database of highly conserved sequence elements in the proximal promoters of genes from multiple species [##REF##16381958##49##], for binding sites of co-expressed transcription factors or of factors expressed in the parent cell type. From this we generated predicted transcription factor interaction lists. Overall, our predictions identified 28% of transcription factor interactions reported in the literature [##REF##14699589##8##], with some binding sites less predictive than others. Randomly generated interaction lists using the same transcription factors identified only 3% of literature reported sites, indicating that our interaction predictions are significantly enriched for true sites (<italic>p </italic>= 5.7 × 10<sup>-6 </sup>using Fisher's exact test). These predicted interactions, as well as literature reported interactions, were then used to build a transcriptional network describing endocrine pancreas development (Figure ##FIG##6##7b##). In sum, a total of 217 novel predicted interactions are shown in this network, suggesting that it is a rich source of hypotheses for future validation. To validate this and to gain an understanding of the false positive rate associated with our interaction predictions, we utilized ChIP-qPCR to validate predicted targets of Foxa2 and Pdx1 in the β-cell portion of the network using Min6 cells. For Foxa2 we found high levels of enrichment for predicted sites in the promoters of <italic>Pdx1</italic>, <italic>Myt1</italic>, <italic>Myt3</italic>, and <italic>Nkx2-2</italic>, but not for predicted sites in the promoters of <italic>Foxa2</italic>, <italic>Neurod1</italic>, or <italic>Nkx6-1 </italic>(Figure ##FIG##6##7c##). For Pdx1 modest levels of enrichment were achieved for sites in the promoters of <italic>Ins1</italic>, <italic>Pdx1</italic>, <italic>Myt1</italic>, <italic>Neurod1</italic>, <italic>Nkx2-2</italic>, and <italic>Nkx6-1</italic>, although no enrichment was seen for a predicted site in the <italic>Foxa2 </italic>promoter (Figure ##FIG##6##7d##). No enrichment was seen at either of two sites that were not predicted to contain a binding site for either Foxa2 or Pdx1. In total, 10 out of 14 (71%) predicted interactions proved valid, with 4 of 7 (57%) validating for Foxa2 and 6 of 7 (86%) for Pdx1.</p>
<title>Discussion</title> <p>Taken together, our data suggest that the pancreas libraries generated and analyzed in this study accurately reflect the cell types intended and are highly comprehensive. Our analysis of the tag distribution of genes with well characterized expression profiles in the libraries confirmed their known expression. As well, 79% of the genes we identified, after instituting a count threshold, have independent validation of pancreas expression via EST information. Of the genes we analyzed using the <italic>GenePaint </italic>database, 20% were not detected, comparable with the expected detection rate of this technique, with this number as low as 13% for clusters with expression at TS22, the time point analyzed. It has been suggested that a SAGE library of 120,000-160,000 tags is roughly equivalent in sensitivity to a microarray [##REF##15475240##50##]. Six of the ten libraries generated here were sequenced to greater than 300,000 tags per library with the remaining four sequenced to approximately 100,000 tags per library. It is clear that these libraries are a rich source of information on the expression profiles of a large diversity of known and unknown transcripts throughout pancreas development. Interestingly, multiple tags were identified for many genes with critical roles in pancreas development and function, including <italic>Ins1</italic>, <italic>Ins2</italic>, <italic>Isl1</italic>, <italic>Neurod1</italic>, and <italic>Pax6</italic>; the biological relevance of this remains to be elucidated. Additionally, we identified 1,049 tags that meet both our count and specificity thresholds and mapped only to the genome. Although some of these tags are likely to be error tags, they represent a rich source of information on completely novel transcripts.</p> <p>To fully appreciate the data generated in this study it is important to consider the cell types the libraries represent. Two of the libraries generated use a transgenic mouse strain that expresses EGFP under control of the <italic>Pdx1 </italic>promoter, while four libraries were generated using a transgenic strain expressing EGFP under control of the <italic>Neurog3 </italic>promoter [##REF##14660441##22##]. Prior to the secondary transition, <italic>Pdx1 </italic>marks all pancreas epithelial cells, with the exception of rare glucagon-expressing cells. At the stages used here, TS17 and TS19, the collected cells therefore represent a relatively pure population of pancreas epithelial progenitor cells, which will eventually give rise to all subsequent endocrine and exocrine cell types in the pancreas. <italic>Pdx1 </italic>is also expressed in the duodenum at these time points and some level of duodenum specific transcripts may also be present, although care was taken in the dissections to minimize this. <italic>Neurog3 </italic>is extremely specific in its expression and in the pancreas is exclusively expressed in endocrine precursors, the cells that give rise to all of the hormone-producing cell types in the pancreas. <italic>Neurog3 </italic>expression itself is transient and decreases substantially after TS23 [##REF##10677506##17##]. However, the EGFP protein is relatively stable and continues to mark endocrine cells as late as TS26 [##REF##18071024##26##]. Thus, the <italic>Neurog3 </italic>EGFP+ libraries collected here, from TS20 to TS22, represent a mixture of endocrine progenitor cell-types, predominantly α and β cell precursors, in various stages of differentiation and maturation. The <italic>Neurog3 </italic>EGFP- TS20 library contains a variety of pancreas cell types, including mesenchymal cells, as well as epithelial and exocrine progenitors. Similarly, the TS22 and TS26 whole pancreas libraries are composed of various cell types, at TS22 predominantly mesenchymal cells and various epithelial cell types, including endocrine and exocrine progenitors, while at TS26 a greater relative abundance of exocrine cells is found. The Islet and Duct libraries represent fully mature cell types. These libraries were constructed from hand-picked islets and ducts isolated by collagenase digestion and gradient centrifugation and, as such, contain a minimum of contaminating exocrine cells. The Islet library is predominantly composed of β-cells that compose approximately 80% of a rodent islet [##REF##14687913##51##], with α-cells composing much of the remaining cells, with smaller numbers of δ, ε, and PP cells also present. For the Duct library, epithelial cells from duct tips were collected, which have been proposed to harbor pancreas stem cells [##REF##10830272##52##], although transcripts from contaminating exocrine and islet cells were inevitably identified as well.</p> <p>One of the major strengths of SAGE data is the ability to easily compare data amongst different libraries. We took advantage of this strength and were able to institute a specificity-based threshold for our data by obtaining the counts for tags that met our count threshold in 205 different SAGE libraries. In general, metrics for specificity assess the level of expression in the tissue of interest versus the total expression. One study used a metric based on Shannon entropy to account for bias from differing levels of expression in different tissues [##REF##15833120##53##]. This metric is not, however, directly applicable to SAGE data due to the large number of tags with no counts in many samples. Here we use a metric that accounts for the relative expression of a tag, the number of libraries it is found in, as well as its absolute expression level. Thus, tags expressed at high-count levels in the library of interest and with low levels in very few other libraries have the highest specificity, while tags with low counts in numerous samples have the lowest specificity. We show that this metric performs well and that it accurately identifies genes with relative levels of restriction compared with available <italic>GenePaint </italic>data, and meets biological expectation.</p> <p>The use of a specificity threshold to assess data from transcriptome-based approaches to pancreas development is novel and previous studies have relied on the identification of genes of interest based on their differential expression in pancreas development. However, differential expression can occur for many reasons and many of the tags that did not meet our specificity threshold showed large changes in count levels between the libraries (data not shown), suggesting that assessing differential expression alone does not screen out many widely expressed genes and that the addition of a specificity threshold is helpful in producing gene lists with the highest functional relevance. For example, of 30 genes identified as being the most significantly differentially expressed from TS21-TS25 in <italic>Neurog3 </italic>EGFP+ cells using the PancChip 6 [##REF##18071024##26##], 10 did not make our count cut-off and 15 did not make our specificity cut-off. Additionally, with a few exceptions - tags for <italic>Mafa </italic>were found at a count of only one in the islet and TS22 <italic>Neurog3 </italic>EGFP+ library and, thus, missed our count threshold - almost all of the transcription factors known to play significant roles in pancreas development were present on our lists. Thus, by instituting our count and specificity thresholds we feel we have generated a tag list more enriched for genes of biological significance to pancreas development and function than previously accomplished.</p> <p>Our data in general confirm the accuracy of our SAGE tag clustering. In our qRT-PCR validations, 10 of the 44 genes assessed had a maximum tag count of only 10 and our qRT-PCR data did not match the SAGE data as well, for these genes, as it did using genes with higher tag counts. Thus, some caution in interpreting the expression profiles of transcripts with low tag counts is necessary. Regardless, it was clear that tags with even very low counts were still indicative of actual expression. Clearly, deeper sequencing, particularly of the libraries at 100,000 tag depth, would improve the statistical power of the observed expression differences and provide a more robust detection of rare transcripts.</p> <p>Five distinct domains of expression have been identified in the developing pancreas [##REF##17609113##40##]. Although, several factors expressed in endocrine cells, which are found in the epithelial trunk, have been identified, relatively few factors are known to be expressed in the remaining expression domains [##REF##14699589##8##,##REF##17609113##40##]. Here we identify and categorize into the 5 expression domains the expression profiles of over 400 genes in the developing pancreas. In sum, we provide spatial and temporal expression data on hundreds of transcripts not previously characterized in pancreas development. Combining our SAGE tag clustering and <italic>GenePaint </italic>data provides a clear indication as to the stages of pancreas development at which many of these genes are likely functionally significant. This, and studies on retinal development [##REF##15226823##54##], demonstrate the power of combining SAGE tag clustering with large scale <italic>in situ </italic>hybridization data.</p> <p>Further, we utilized predicted and literature reported binding data to gain insight into the cascade of transcription factors that drives endocrine pancreas development. As already mentioned, our predicted interactions identified 28% of literature reported sites. The CisRED database used to generate our interaction lists identifies conserved sequence elements only in proximal promoters (approximately 1 kb upstream). Of course, not all regulatory elements are in proximal promoters and this is likely a major reason a higher percentage of known interactions was not identified. Future experimentation to locate and analyze regulatory elements across the genome will better predict all possible interactions. Optimally false positive error rates in such networks are calculated using comparisons with comprehensive binding data for multiple factors. However, this type of data is not available for the relevant transcription factors in pancreas tissues. We therefore attempted to gain an understanding of our false positive rate by utilizing ChIP-qPCR to validate predicted binding sites for Foxa2 and Pdx1. Although limited in scale, these validation experiments give us a reasonable estimation of our false positive rate as approximately 29%. However, this error rate varied significantly from 14% for Pdx1 to 46% for Foxa2. In part, this is due to our reliance on available binding data for the different transcription factors in the network, and thus the reliability of the predicted interactions is highly dependent on the quality of the available binding data, which varies significantly. It is also possible that some of the interactions we tested do not occur in the cell line used, but do in fact occur <italic>in vivo</italic>, although this remains to be tested. Regardless, it is a given that better binding data for a greater number of transcription factors are essential to the construction of more accurate networks in this fashion. Despite this, it is clearly of great interest to begin to elucidate the relevance of many of these predicted interactions in pancreas development and function.</p>
<title>Conclusion</title> <p>Previous studies have used microarrays to assess expression profiles in the developing pancreas. Here we utilized SAGE to provide a more comprehensive and complementary analysis, and to yield further insights into signals critical in regulating pancreas development. By instituting a specificity threshold for our data we have isolated a list of tags highly enriched for genes with significant roles in pancreas development and function. Equivalent analyses have not been performed in previous array studies, and highlight one of the strengths of SAGE as such comparisons are relatively straightforward with this technique. Our data represent a highly comprehensive analysis of pancreas development. Our classification of the <italic>in situ </italic>staining patterns of over 400 genes, in combination with our cluster analysis, allowed us to identify gene sets describing each of the major stages of pancreas development represented by our libraries. These data provide evidence for the pancreas enriched expression of hundreds of factors that have not previously been implicated in pancreatic development, suggesting that our data are a rich source of genes with novel roles in pancreas development. Further, from our data we constructed a predictive transcriptional network driving endocrine pancreas development that predicts many novel interactions that are of great interest for future analysis. In sum, these data provide insight into the regulatory networks driving pancreas development and function and provide a framework for future studies.</p>
<p>This is an open access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <p>The expression profile of different developmental stages of the murine pancreas and predictions of transcription factor interactions, provides a framework for pancreas regulatory networks and development.</p> <title>Background</title> <p>Despite recent advances, the transcriptional hierarchy driving pancreas organogenesis remains largely unknown, in part due to the paucity of comprehensive analyses. To address this deficit we generated ten SAGE libraries from the developing murine pancreas spanning Theiler stages 17-26, making use of available <italic>Pdx1 </italic>enhanced green fluorescent protein (EGFP) and <italic>Neurog3 </italic>EGFP reporter strains, as well as tissue from adult islets and ducts.</p> <title>Results</title> <p>We used a specificity metric to identify 2,536 tags with pancreas-enriched expression compared to 195 other mouse SAGE libraries. We subsequently grouped co-expressed transcripts with differential expression during pancreas development using <italic>K</italic>-means clustering. We validated the clusters first using quantitative real time PCR and then by analyzing the Theiler stage 22 pancreas <italic>in situ </italic>hybridization staining patterns of over 600 of the identified genes using the <italic>GenePaint </italic>database. These were then categorized into one of the five expression domains within the developing pancreas. Based on these results we identified a cascade of transcriptional regulators expressed in the endocrine pancreas lineage and, from this, we developed a predictive regulatory network describing beta-cell development.</p> <title>Conclusion</title> <p>Taken together, this work provides evidence that the SAGE libraries generated here are a valuable resource for continuing to elucidate the molecular mechanisms regulating pancreas development. Furthermore, our studies provide a comprehensive analysis of pancreas development, and insights into the regulatory networks driving this process are revealed.</p>
<title>Abbreviations</title> <p>ChIP, chromatin immunoprecipitation; EGFP, enhanced green fluorescent protein; EST, expressed sequence tag; FACS, fluorescence activated cell sorting; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; PBS, phosphate buffered saline; qRT-PCR, quantitative real time PCR; S, specificity; SAGE, Serial analysis of gene expression; SOTA, self-organizing tree algorithm; TS, Theiler stage.</p> <title>Authors' contributions</title> <p>BH was involved in all aspects of the project, including design, tissue collection, data analysis, qRT-PCR validations, <italic>GenePaint </italic>analysis, and so on. BZ assisted with qRT-PCR and <italic>in situ </italic>studies, and performed the <italic>Foxa2 </italic>ChIP experiment. JW performed tissue collection. TA performed tissue collection. MB performed the <italic>Pdx1 </italic>ChIP experiment. PH is a senior author. SJ is a senior author. MM is a senior author. CH is a senior author and was involved in experimental design and planning.</p> <title>Additional data files</title> <p>The following additional data files are available. Additional data file ##SUPPL##0##1## is a figure showing a comparison of the number and percentage of tags that map to known pancreas ESTs or transcripts as a function of minimum tag count. Additional data file ##SUPPL##1##2## is an Excel spreadsheet showing the tag counts for the 2,536 tags used in our SAGE tag <italic>K</italic>-means cluster analysis. Additional data file ##SUPPL##2##3## is a figure showing the percentage of genes identified as <italic>Pdx1 </italic>EGFP+, Neurog EGFP+, or Islet enriched by Gu <italic>et al</italic>. [##REF##14660441##22##] in each of the identified <italic>K</italic>-means cluster groups. Additional data file ##SUPPL##3##4## is a table showing the <italic>GenePaint</italic>-based staining classification results. Additional data file ##SUPPL##4##5## is a table listing the genes and primer sequences used in qRT-PCR validation studies. Additional data file ##SUPPL##5##6## is a table listing the obtained binding information for identified transcription factors.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>The authors would like to acknowledge members of the BC Cancer Agency Genome Sciences Centre sequencing, bioinformatics, and SAGE library construction teams. In addition, we are particularly grateful to Douglas Melton and Guoqiang Gu for providing the <italic>Pdx1 </italic>EGFP and <italic>Ngn3 </italic>EGFP mouse strains, as well as to Bruce Verchere and Galina Smirnova for their assistance in islet and duct collection. We would also like to thank Genome BC for their assistance in project management and the staff of the Animal Resource Center at the BBCRC. MAM, PAH, SJ, and CDH are scholars of the Michael Smith Foundation for Health Research. PAH is also a CIHR New Investigator. Funding was provided by the Genome Canada contracts 'Mouse Atlas of Development' and 'MORGEN' and the National Institute of Health (USA), with infrastructure support provided by the BC Cancer Foundation.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p>Heatmap of SAGE tag counts for genes with known expression profiles in pancreas development. Tags for genes with well characterized expression profiles in pancreas development were identified and their normalized counts obtained in each of the ten SAGE libraries created. A heatmap, generated using the multi-experiment viewer as described in the Materials and methods, of these results is shown based on the counts of the tags per hundred thousand (TPH). SAGE tags used include: TACACGTTCTGACAACT (<italic>Nkx2-2</italic>); AAGTGGAAAAAAGAGGA (<italic>Pdx1</italic>); TAGTTTTAACAGAAAAC (<italic>Foxa2</italic>); ACCTTCACACCAAACAT (<italic>Hnf4a</italic>); AATGCAGAGGAGGACTC (<italic>Neurod1</italic>); CAGGGTTTCTGAGCTTC (<italic>Neurog3</italic>); TCATTTGACTTTTTTTT (<italic>Isl1</italic>); GATTTAAGAGTTTTATC (<italic>Pax6</italic>); CAGCAGGACGGACTCAG (<italic>Pax4</italic>); CAGTCCATCAACGACGC (<italic>Ptf1a</italic>); AGAAACAGCAGGGCCTG (<italic>Bhlhb8</italic>); GACCACACTGTCAAACA (<italic>Cpa1</italic>); CCCTGGGTTCAGGAGAT (<italic>Ctrb1</italic>); TTGCGCTTCCTGGTGTT (<italic>Ela1</italic>); ACCACCTGGTAACCGTA (<italic>Gcg</italic>); GCCGGGCCCTGGGGAAG (<italic>Ghrl</italic>); CTAAGAATTGCTTTAAA (<italic>Iapp</italic>); GCCCTGTTGGTGCACTT (<italic>Ins1</italic>); TCCCGCCGTGAAGTGGA (<italic>Ins2</italic>). The libraries shown include: <italic>Pdx1 </italic>EGFP+ TS17 (P+ TS17); <italic>Pdx1 </italic>EGFP+ TS19 (P+ TS19); <italic>Neurog3 </italic>EGFP- TS20 (N- TS20); <italic>Neurog3 </italic>EGFP+ TS20 (N+ TS20); <italic>Neurog3 </italic>EGFP+ TS21 (N+ TS21); <italic>Neurog3 </italic>EGFP+ TS22 (N+ TS22); whole pancreas TS22 (WTS22); whole pancreas TS26 (WTS26); adult isolated ducts (Ducts); adult isolated islets (Islets).</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p>Specificity threshold accurately predicts spatial expression restriction. A plot of specificity (S) versus cumulative tag types represented shows the distribution of tags into tags with high (S &gt; 0.1; top), medium (0.001 &gt; S &lt; 0.1, middle), and low (S &lt; 0.001, bottom) S values. Representative <italic>in situ </italic>hybridization staining patterns from TS22 whole embryo saggital sections obtained from <italic>GenePaint </italic>are shown for each specificity group. Relevant <italic>GenePaint </italic>probe IDs can be found in Additional data file 4. Arrows indicate the location of the pancreas (p).</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p>Median plots of identified SAGE tag <italic>K</italic>-means cluster analysis using 14 clusters. We clustered 2, 536 SAGE tags with a count greater than 4 in one of the SAGE libraries and with a minimum specificity of 0.002 and that map unambiguously to a specific transcript or genome location into 14 clusters using <italic>K</italic>-means clustering using a PoissonC algorithm as described in the Materials and methods. The median normalized tag counts for the tags in each of the clusters is shown plotted against the indicated SAGE libraries. The libraries shown include: <italic>Pdx1 </italic>EGFP+ TS17 (P+ TS17); <italic>Pdx1 </italic>EGFP+ TS19 (P+ TS19); <italic>Neurog3 </italic>EGFP- TS20 (N- TS20); <italic>Neurog3 </italic>EGFP+ TS20 (N+ TS20); <italic>Neurog3 </italic>EGFP+ TS21 (N+ TS21); <italic>Neurog3 </italic>EGFP+ TS22 (N+ TS22); whole pancreas TS22 (WTS22); whole pancreas TS26 (WTS26); adult isolated ducts (Ducts); adult isolated islets (Islets). A full list of the tags, the cluster they belong to, and their counts in each of the libraries is shown in Additional data file 2.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p>SOTA clustering of temporal expression profiles from qRT-PCR analysis of 44 genes in pancreas development. qRT-PCR was used to determine the relative expression levels of the indicated genes during pancreas development at the TSs indicated. The relative level of expression of each gene was normalized and a SOTA analysis used to group the genes. Heatmaps of the relative expression levels of the genes in the SOTA groups, including the SOTA centroid, with peak expression in <bold>(a) </bold>the islets, <bold>(b) </bold>the TS26 developing pancreas, <bold>(c) </bold>the TS21-TS26 developing pancreas, or <bold>(d) </bold>the ducts are shown. The data shown are averages of the results obtained from pancreases from three separate litters (pancreases from an individual litter were pooled) or islet/duct collections with triplicate reactions from the separate RNA extractions.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p>Representative <italic>in situ </italic>staining patterns for genes expressed in each of the identified expression profiles. Representative genes for each of the identified spatial expression profiles, including genes with known and previously un-described, or novel, staining profiles in pancreas development, are shown. For this, images of <italic>in situ </italic>hybridization staining patterns for whole embryo sagittal sections were obtained from the <italic>GenePaint </italic>website and magnified to show the pancreas (outlined in red). Relevant <italic>GenePaint </italic>probe IDs can be found in Additional data file 4.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p>The percent association of genes in each <italic>K</italic>-means cluster with the five expression domains in the pancreas. The <italic>in situ </italic>hybridization staining profiles of 605 genes with informative stain in TS22 pancreas tissue were classified into the groups shown using the <italic>GenePaint </italic>database. Additional data file 4 lists the full categorization of each of these genes. The percentage of genes with each staining profile in each of the SAGE tag <italic>K</italic>-means cluster is shown.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p>A cascade of transcription factors expressed in endocrine cell development. Tags that mapped unambiguously to transcription factors and met our count and specificity thresholds were identified and their normalized counts obtained in each of the SAGE libraries that represent the development of the endocrine lineage. <bold>(a) </bold>A heatmap, generated using the multi-experiment viewer as described in the Materials and methods, of these results is shown. Tags are organized by the order of their expression in the libraries. The libraries shown include: <italic>Pdx1 </italic>EGFP+ TS17 (P+ TS17); <italic>Pdx1 </italic>EGFP+ TS19 (P+ TS19); <italic>Neurog3 </italic>EGFP+ TS20 (N+ TS20); <italic>Neurog3 </italic>EGFP+ TS21 (N+ TS21); <italic>Neurog3 </italic>EGFP+ TS22 (N+ TS22); adult isolated islets (Islets). <bold>(b) </bold>A transcription factor network in β-cell development. The network was generated as described in the materials and methods using Biotapestry to visualize the network [65]. Ins1 and Ins2 are shown as the final products expressed in mature β-cells. Literature reported interactions [8,41] are included with arrow heads indicating positive regulation and perpendicular lines repression. Interactions in pancreatic progenitors (PP), endocrine progenitors (EP) and mature β-cells are shown. <bold>(c) </bold>The fold enrichment of predicted Foxa2 binding sites, as detected by qPCR, in the promoters of <italic>Pdx1</italic>, <italic>Myt1</italic>, <italic>Myt3</italic>, <italic>Neurod1</italic>, <italic>Nkx2-2</italic>, <italic>Nkx6-1</italic>, and of two sites not predicted to contain a Foxa2 site (<italic>NegF1 </italic>and <italic>NegF2</italic>) in ChIP experiments using an anti-Foxa2 antibody compared to IgG. <bold>(d) </bold>The fold enrichment of predicted Pdx1 binding sites, as detected by qPCR, in the promoters of <italic>Foxa2</italic>, <italic>Ins1</italic>, <italic>Myt1</italic>, <italic>Neurod1</italic>, <italic>Nkx2-2</italic>, <italic>Nkx6-1</italic>, and of two sites not predicted to contain a Pdx1 site (<italic>NegP1 </italic>and <italic>NegP2</italic>) in ChIP experiments using an anti-Pdx1 antibody compared to IgG. The data shown are averages of the results obtained from three separate ChIP experiments each with duplicate reactions. Error bars indicate the standard deviation of the averaged replicates.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Summary of pancreas SAGE libraries generated</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Accession</td><td align="left">Stage</td><td align="left">Tissue subtype</td><td align="left">Cell types represented</td><td align="left">Library type</td><td align="center">Tags sequenced*</td><td align="center">Tag types</td></tr></thead><tbody><tr><td align="left">SM161/SM244</td><td align="left">TS17</td><td align="left"><italic>Pdx1 </italic>EGFP+<sup>†</sup></td><td align="left">All pancreas epithelial cells with the exception of rare Glucagon-positive cells</td><td align="left">Long SAGElite</td><td align="center">306,588</td><td align="center">44,491</td></tr><tr><td align="left">SM231</td><td align="left">TS19</td><td align="left"><italic>Pdx1 </italic>EGFP+</td><td align="left">All pancreas epithelial cells with the exception of rare Glucagon-positive cells</td><td align="left">Long SAGElite</td><td align="center">317,716</td><td align="center">49,572</td></tr><tr><td align="left">SM162/SM245</td><td align="left">TS20</td><td align="left"><italic>Ngn3 </italic>EGFP-<sup>†</sup></td><td align="left">A mixture of pancreas cell types composed predominantly of mesenchymal cells and pancreas epithelial progenitors as well as those destined to become exocrine cell types</td><td align="left">Long SAGElite</td><td align="center">308,745</td><td align="center">47,695</td></tr><tr><td align="left">SM243/SM160</td><td align="left">TS20</td><td align="left"><italic>Ngn3 </italic>EGFP+</td><td align="left">All endocrine progenitor cells as well as endocrine cells at various stages of maturation</td><td align="left">Long SAGElite</td><td align="center">320,473</td><td align="center">51,847</td></tr><tr><td align="left">SM225/SM249</td><td align="left">TS21</td><td align="left"><italic>Ngn3 </italic>EGFP+</td><td align="left">All endocrine progenitor cells as well as endocrine cells at various stages of maturation</td><td align="left">Long SAGElite</td><td align="center">313,503</td><td align="center">58,864</td></tr><tr><td align="left">SM232</td><td align="left">TS22</td><td align="left"><italic>Ngn3 </italic>EGFP+</td><td align="left">All endocrine progenitor cells as well as endocrine cells at various stages of maturation</td><td align="left">Long SAGElite</td><td align="center">301,222</td><td align="center">37,726</td></tr><tr><td align="left">SM223</td><td align="left">TS22</td><td align="left">Whole</td><td align="left">A mixture of pancreas cell types composed predominantly of pancreas epithelial cells differentiating into exocrine cell types with some endocrine cells and mesenchymal cells</td><td align="left">Long SAGE</td><td align="center">98,189</td><td align="center">13,676</td></tr><tr><td align="left">SM016</td><td align="left">TS26</td><td align="left">Whole</td><td align="left">A mixture of pancreas cell types composed predominantly of pancreas epithelial cells differentiating into exocrine cell types with some endocrine cells and mesenchymal cells</td><td align="left">Long SAGE</td><td align="center">81,130</td><td align="center">17,963</td></tr><tr><td align="left">SM102</td><td align="left">DPN70</td><td align="left">Isolated ducts</td><td align="left">Hand picked adult ducts isolated by collagenase treatment and gradient centrifugation</td><td align="left">Long SAGE</td><td align="center">119,024</td><td align="center">23,528</td></tr><tr><td align="left">SM017</td><td align="left">DPN70</td><td align="left">Isolated islets</td><td align="left">Hand picked adult islets isolated by collagenase treatment and gradient centrifugation composed of each of the major endocrine cell types</td><td align="left">Long SAGE</td><td align="center">99,968</td><td align="center">16,039</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Top 25 most specific transcripts in the pancreas SAGE libraries</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Tag</td><td align="left">Accession/location</td><td align="left">Symbol</td><td align="center"><italic>Pdx1</italic>-GFP+ (TS17)</td><td align="center"><italic>Pdx1</italic>-GFP+ (TS19)</td><td align="center"><italic>Neurog3</italic>-GFP- (TS20)</td><td align="center"><italic>Neurog3</italic>-GFP+ (TS20)</td><td align="center"><italic>Neurog3</italic>-GFP+ (TS 21)</td><td align="center"><italic>Neurog3</italic>-GFP+ (TS22)</td><td align="center">Whole (TS22)</td><td align="center">Whole (TS26)</td><td align="center">Ducts</td><td align="center">Islets</td><td align="center">MaxS<sup>†</sup></td></tr></thead><tbody><tr><td align="left">TCCCGCCGTGAAGTGGA</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_008387">NM_008387</ext-link></td><td align="left"><italic>Ins2</italic></td><td align="center">0*</td><td align="center">0.31</td><td align="center">0</td><td align="center">13.11</td><td align="center">57.43</td><td align="center">3,298.9</td><td align="center">4.07</td><td align="center">139.28</td><td align="center">1,422.4</td><td align="center">2,2471.19</td><td align="center">62.72</td></tr><tr><td align="left">TTCTGTCTGGGCTTCCT</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_023333">NM_023333</ext-link></td><td align="left"><italic>2210010C04Rik</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">77.65</td><td align="center">651.97</td><td align="center">109.03</td><td align="center">33.43</td></tr><tr><td align="left">GCCCTGTTGGTGCACTT</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_008386">NM_008386</ext-link></td><td align="left"><italic>Ins1</italic></td><td align="center">0</td><td align="center">2.83</td><td align="center">0</td><td align="center">45.56</td><td align="center">19.59</td><td align="center">839.25</td><td align="center">6.11</td><td align="center">9.86</td><td align="center">207.52</td><td align="center">3,116</td><td align="center">27.90</td></tr><tr><td align="left">TTAGGAGGCTGCTGCTG</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_026925">NM_026925</ext-link></td><td align="left"><italic>Pnlip</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">1,760.99</td><td align="center">116.04</td><td align="center">18.10</td></tr><tr><td align="left">CCCTGGGTTCAGGAGAT</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_025583">NM_025583</ext-link></td><td align="left"><italic>Ctrb1</italic></td><td align="center">0</td><td align="center">0.31</td><td align="center">0</td><td align="center">0</td><td align="center">31.21</td><td align="center">74.36</td><td align="center">17.31</td><td align="center">3,162.83</td><td align="center">1,443.41</td><td align="center">385.12</td><td align="center">18.05</td></tr><tr><td align="left">GCCCTGTGGATGCGCTT</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_008387">NM_008387</ext-link></td><td align="left"><italic>Ins2</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0.33</td><td align="center">16.27</td><td align="center">0</td><td align="center">0</td><td align="center">15.96</td><td align="center">432.14</td><td align="center">17.58</td></tr><tr><td align="left">GTGTGCGCTGGTGGCGA</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_007919">NM_007919</ext-link></td><td align="left"><italic>Ela2</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">69.03</td><td align="center">181.48</td><td align="center">4</td><td align="center">11.75</td></tr><tr><td align="left">GCATCGTGAGCTTCGGC</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_007919">NM_007919</ext-link></td><td align="left"><italic>Ela2</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">2.32</td><td align="center">0</td><td align="center">1,329.96</td><td align="center">2,680.13</td><td align="center">1,156.37</td><td align="center">11.24</td></tr><tr><td align="left">GTGTGCGCCGGCGGCGA</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_026419">NM_026419</ext-link></td><td align="left"><italic>Ela3</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">1</td><td align="center">1.02</td><td align="center">636.02</td><td align="center">369.67</td><td align="center">23.01</td><td align="center">11.14</td></tr><tr><td align="left">ACCACCTGGTAACCGTA</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_008100">NM_008100</ext-link></td><td align="left"><italic>Gcg</italic></td><td align="center">7.5</td><td align="center">63.26</td><td align="center">0.65</td><td align="center">2,554.97</td><td align="center">1,952.71</td><td align="center">550.42</td><td align="center">34.63</td><td align="center">25.88</td><td align="center">124.34</td><td align="center">326.1</td><td align="center">10.99</td></tr><tr><td align="left">AAAGTATGCAAATAGCT</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_026918">NM_026918</ext-link></td><td align="left"><italic>1810010M01Rik</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">194.75</td><td align="center">934.27</td><td align="center">459.15</td><td align="center">9.90</td></tr><tr><td align="left">CAGACTAAGTACCCATA</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_009885">NM_009885</ext-link></td><td align="left"><italic>Cel</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0.66</td><td align="center">1</td><td align="center">0</td><td align="center">750.65</td><td align="center">375.55</td><td align="center">16.01</td><td align="center">8.81</td></tr><tr><td align="left">TTTTACTTCTAAGAGTC</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_021331">NM_021331</ext-link></td><td align="left"><italic>G6pc2</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0.31</td><td align="center">0</td><td align="center">3.32</td><td align="center">0</td><td align="center">0</td><td align="center">5.88</td><td align="center">221.07</td><td align="center">7.74</td></tr><tr><td align="left">CCCGGGTGCAAGAAGAA</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_018874">NM_018874</ext-link></td><td align="left"><italic>Pnliprp1</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">5.93</td><td align="center">12.62</td><td align="center">18.26</td><td align="center">16.3</td><td align="center">1,135.22</td><td align="center">250.37</td><td align="center">8</td><td align="center">7.40</td></tr><tr><td align="left">TCCCTTCAACCTTAGAC</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_011271">NM_011271</ext-link></td><td align="left"><italic>Rnase1</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0.33</td><td align="center">0</td><td align="center">221.87</td><td align="center">1,249.33</td><td align="center">170.05</td><td align="center">6.48</td></tr><tr><td align="left">TTAAACCAGAGTTCATA</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_023333">NM_023333</ext-link></td><td align="left"><italic>2210010C04Rik</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">10.08</td><td align="center">0</td><td align="center">5.66</td></tr><tr><td align="left">GCCTACAACTAAACTGT</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_023182">NM_023182</ext-link></td><td align="left"><italic>Ctrl</italic></td><td align="center">0</td><td align="center">0.31</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">27.12</td><td align="center">491.5</td><td align="center">195.06</td><td align="center">5.46</td></tr><tr><td align="left">GCACCAAGTACACATAT</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_029706">NM_029706</ext-link></td><td align="left"><italic>Cpb1</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">303.22</td><td align="center">209.2</td><td align="center">21.01</td><td align="center">5.11</td></tr><tr><td align="left">TTGCGCTTCCTGGTGTT</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_033612">NM_033612</ext-link></td><td align="left"><italic>Ela1</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">8.4</td><td align="center">0</td><td align="center">4.93</td></tr><tr><td align="left">TGGGAGTGGAGGATGCC</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_026925">NM_026925</ext-link></td><td align="left"><italic>Pnlip</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">29.41</td><td align="center">9</td><td align="center">4.83</td></tr><tr><td align="left">TTCCAAGTGGAGGAGGT</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_018874">NM_018874</ext-link></td><td align="left"><italic>Pnliprp1</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0.31</td><td align="center">0</td><td align="center">0</td><td align="center">10.18</td><td align="center">163.93</td><td align="center">36.97</td><td align="center">1</td><td align="center">4.78</td></tr><tr><td align="left">CTAAGAATTGCTTTAAA</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_010491">NM_010491</ext-link></td><td align="left"><italic>Iapp</italic></td><td align="center">0</td><td align="center">0.31</td><td align="center">0</td><td align="center">3.43</td><td align="center">6.64</td><td align="center">49.8</td><td align="center">0</td><td align="center">2.47</td><td align="center">25.21</td><td align="center">170.05</td><td align="center">4.50</td></tr><tr><td align="left">CAGTCCATCAACGACGC</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_018809">NM_018809</ext-link></td><td align="left"><italic>Ptf1a</italic></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">7.13</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">4.36</td></tr><tr><td align="left">CAAAGAATGCAATCTGA</td><td align="left"><ext-link ext-link-type="gen" xlink:href="nt_039700">nt_039700</ext-link></td><td/><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">7.13</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">4.36</td></tr><tr><td align="left">CTTGCAGTCTGAGTTCG</td><td align="left"><ext-link ext-link-type="gen" xlink:href="nt_039413">nt_039413</ext-link></td><td/><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">7.13</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">4.36</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Summary of SAGE tag <italic>K</italic>-means cluster data</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Cluster</td><td align="center">Number of tags in the cluster</td><td align="center">Number of genes in the cluster</td><td align="center">Number of genome maps in the cluster</td><td align="center">Number assessed by <italic>GenePaint</italic>*</td><td align="center">Number assessed by QPCR</td><td align="center">Median S<sup>†</sup></td><td align="left">Previously characterized genes in the cluster</td><td align="left">Selected GO categories and KEGG pathways enriched in the cluster<sup>‡</sup></td></tr></thead><tbody><tr><td align="left">1</td><td align="center">154</td><td align="center">85</td><td align="center">61</td><td align="center">37</td><td align="center">3</td><td align="center">0.0079</td><td align="left"><italic>Nkx6-2</italic></td><td align="left">Transcriptional activator activity <italic>p </italic>= 0.02; development <italic>p </italic>= 0.049</td></tr><tr><td align="left">2</td><td align="center">49</td><td align="center">29</td><td align="center">19</td><td align="center">15</td><td align="center">2</td><td align="center">0.0037</td><td/><td align="left">Metabolism <italic>p </italic>= 0.01; cell organization and biogenesis <italic>p </italic>= 0.035</td></tr><tr><td align="left">3</td><td align="center">58</td><td align="center">40</td><td align="center">15</td><td align="center">14</td><td align="center">2</td><td align="center">0.0044</td><td/><td align="left">Receptor activity <italic>p </italic>= 0.028; development <italic>p </italic>= 0.030</td></tr><tr><td align="left">4</td><td align="center">292</td><td align="center">115</td><td align="center">175</td><td align="center">45</td><td align="center">4</td><td align="center">0.00895</td><td align="left"><italic>Hes6</italic>, <italic>Pdx1</italic>, <italic>Sox9</italic></td><td align="left">Regulation of transcription <italic>p </italic>= 0.027; maturity onset diabetes of the young <italic>p </italic>= 0.002</td></tr><tr><td align="left">5</td><td align="center">1,008</td><td align="center">427</td><td align="center">542</td><td align="center">175</td><td align="center">13</td><td align="center">0.03555</td><td align="left"><italic>Arx</italic>, <italic>Gcg</italic>, <italic>Ghrl</italic>, <italic>Iapp</italic>, <italic>Isl1</italic>, <italic>Nkx2-2</italic>, <italic>Myt1</italic>, <italic>Neurog3</italic>, <italic>Neurod1</italic>, <italic>Pax4</italic>, <italic>Pax6</italic>, <italic>Pou3f4</italic>, <italic>Pyy</italic></td><td align="left">Secretory pathway <italic>p </italic>&lt; 0.001; hormone activity <italic>p </italic>= 0.049; maturity onset diabetes of the young <italic>p </italic>&lt; 0.001</td></tr><tr><td align="left">6</td><td align="center">60</td><td align="center">41</td><td align="center">17</td><td align="center">16</td><td align="center">1</td><td align="center">0.00465</td><td/><td/></tr><tr><td align="left">7</td><td align="center">21</td><td align="center">11</td><td align="center">10</td><td align="center">7</td><td align="center">1</td><td align="center">0.008</td><td/><td/></tr><tr><td align="left">8</td><td align="center">78</td><td align="center">46</td><td align="center">28</td><td align="center">29</td><td align="center">5</td><td align="center">0.012</td><td align="left"><italic>Pax6</italic></td><td align="left">Eye morphogenesis <italic>p </italic>= 0.020; type II diabetes mellitus <italic>p </italic>= 0.001</td></tr><tr><td align="left">9</td><td align="center">23</td><td align="center">16</td><td align="center">6</td><td align="center">10</td><td align="center">2</td><td align="center">0.0041</td><td/><td align="left">Cell proliferation <italic>p </italic>= 0.028</td></tr><tr><td align="left">10</td><td align="center">401</td><td align="center">281</td><td align="center">107</td><td align="center">122</td><td align="center">4</td><td align="center">0.0158</td><td align="left"><italic>Id2</italic></td><td align="left">Response to endogenous stimulus <italic>p </italic>= 0.021</td></tr><tr><td align="left">11</td><td align="center">76</td><td align="center">57</td><td align="center">10</td><td align="center">23</td><td align="center">1</td><td align="center">0.00555</td><td align="left"><italic>Amy1</italic>, <italic>Cel</italic>, <italic>Clps</italic>, <italic>Ela1</italic>, <italic>Pnliprp2</italic>, <italic>Reg1</italic></td><td align="left">Protein catabolism <italic>p </italic>= 0.002</td></tr><tr><td align="left">12</td><td align="center">154</td><td align="center">122</td><td align="center">13</td><td align="center">56</td><td align="center">3</td><td align="center">0.0074</td><td align="left"><italic>Ela1</italic>, <italic>Pnlip</italic>, <italic>Reg3d</italic></td><td align="left">Growth factor binding <italic>p </italic>= 0.005; carboxypeptidase activity <italic>p </italic>= 0.013; regulation of cell growth <italic>p </italic>= 0.027</td></tr><tr><td align="left">13</td><td align="center">136</td><td align="center">84</td><td align="center">42</td><td align="center">30</td><td align="center">3</td><td align="center">0.01835</td><td align="left"><italic>Iapp</italic>, <italic>Ins1</italic>, <italic>Ins2</italic></td><td align="left">Secretion <italic>p </italic>= 0.03; maturity onset diabetes of the young <italic>p </italic>&lt; 0.001; type II diabetes mellitus <italic>p </italic>&lt; 0.001; type I diabetes mellitus <italic>p </italic>= 0.003</td></tr><tr><td align="left">14</td><td align="center">56</td><td align="center">47</td><td align="center">4</td><td align="center">22</td><td align="center">0</td><td align="center">0.00335</td><td/><td align="left">Protein metabolism <italic>p </italic>= 0.020</td></tr></tbody></table></table-wrap>
<disp-formula>S = M<sub>r </sub>log<sub>1.3</sub>(C<sub>i</sub>)/L<sub>a</sub></disp-formula>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional data file 1</title><p>Number and percentage of tags that map to known pancreas ESTs or transcripts as a function of minimum tag count.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional data file 2</title><p>Tag counts for the 2,536 tags used in our SAGE tag <italic>K</italic>-means cluster analysis.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional data file 3</title><p>Percentage of genes identified as <italic>Pdx1 </italic>EGFP+, Neurog EGFP+, or Islet enriched by Gu <italic>et al</italic>. [##REF##14660441##22##] in each of the identified <italic>K</italic>-means cluster groups.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional data file 4</title><p><italic>GenePaint</italic>-based staining classification results.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional data file 5</title><p>Genes and primer sequences used in qRT-PCR validation studies.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional data file 6</title><p>Binding information for identified transcription factors.</p></caption></supplementary-material>
<table-wrap-foot><p>*After 95% quality cutoffs for all tags. <sup>†</sup>The <italic>Pdx1 </italic>EGFP and <italic>Ngn3 </italic>EGFP transgenic strains were obtained from Douglas Melton as described in Gu <italic>et al</italic>. [22]. DPN, days post natal.</p></table-wrap-foot> <table-wrap-foot><p>*Tag counts are shown as tags per 100,000. This indicates the total number of times a given SAGE tag appears in the library per 100,000 tags and is used to normalize for libraries of varying size. <sup>†</sup>S is the specificity of the tag. Specificity is calculated as described in the Materials and methods. The maximum S in any one of the libraries created here is indicated.</p></table-wrap-foot> <table-wrap-foot><p>*Refers to the number of genes analyzed by <italic>in situ </italic>hybridization using <italic>GenePaint </italic>[62] on TS22 whole embryo cryo-sections that gave informative staining. <sup>†</sup>S is the specificity of the tag. Specificity is calculated as described in the Materials and methods. <sup>‡</sup>GO term enrichments and <italic>p</italic>-values were calculated using EASE while KEGG pathway enrichments and <italic>p</italic>-values using Webgestalt as described in the Materials and methods.</p></table-wrap-foot>
<graphic xlink:href="gb-2008-9-6-r99-1"/> <graphic xlink:href="gb-2008-9-6-r99-2"/> <graphic xlink:href="gb-2008-9-6-r99-3"/> <graphic xlink:href="gb-2008-9-6-r99-4"/> <graphic xlink:href="gb-2008-9-6-r99-5"/> <graphic xlink:href="gb-2008-9-6-r99-6"/> <graphic xlink:href="gb-2008-9-6-r99-7"/>
<media xlink:href="gb-2008-9-6-r99-S1.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r99-S2.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r99-S3.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r99-S4.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r99-S5.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r99-S6.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"article-title": ["The Mouse Atlas Website"]}, {"article-title": ["CGAP SAGE Website"]}, {"article-title": ["The Refseq Database"]}, {"article-title": ["The Ensembl Database"]}, {"article-title": ["The Mammalian Gene Collection"]}, {"article-title": ["EPConDB/T1DBase"]}, {"article-title": ["GenePaint Database"]}, {"article-title": ["Primer3"]}]
Genome Biol. 2008 Jun 14; 9(6):R99