Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Invest Ophthalmol Vis Sci. Author manuscript; available in PMC 2010 July 19.
Published in final edited form as:
PMCID: PMC2906229

Differential Gene Expression in the Pig Limbal Side Population: Implications for Stem Cell Cycling, Replication, and Survival



To define the molecular signature of limbal SP cells and identify signaling pathways associated with the phenotype of these putative stem cells.


Primary cultures of pig limbal epithelial cells stained with Hoechst 33342 were sorted by flow cytometry into SP and non-SP cells, and purified RNA was processed for microarray analysis with an oligonucleotide spotted array. Expressed transcripts for which SP and non-SP expressions differed by more that 1.5-fold in each paired set and by twofold overall were considered to be differentially expressed. Differential expression was validated by quantitative PCR and immunostaining. Data-mining methods were used to identify cellular processes that are either salient or depressed in the SP cells.


The microarray identified approximately 9000 distinct, expressed, and identifiable genes. Of those, 382 and 296 were either over- or underexpressed in the SP cells, respectively. Overrepresentation analysis indicated that SP cells are in a low metabolic and biosynthetic state. In addition, a pattern of elevated MXD1, MAXI2, DUSP5, p27/KIP1, and p57/KIP2 and decreased Cyclin D and CDK genes can be expected to slow intrinsic and mitogen-induced G1-to-S cell cycle transition. SP cells were also rich in genes associated with stem cell phenotype and genes providing protection against oxidative and/or xenobiotic damage.


Microarray analysis of pig limbal SP cells yielded a molecular signature underscoring a phenotype characterized by slow cycling and low metabolic activity. The results provide valuable insights for the preservation and/or replication of epithelial stem cells.

Stem cells are critical for the function of the ocular surface. In the limbocorneal system, stem cells are concentrated in the limbus.15 Limbal damage due to chemical or thermal injury, microbial infection, or autoimmune reactions results in limbal stem cell deficiencies reflected in corneal opacities, neovascularization, and/or generalized inflammation.68 Advances in the experimental research and clinical application of limbal transplantation for the treatment of limbal stem deficiencies have been rapid, progressing from the straightforward transplantation of contralateral biopsies to pre-expansion of the donor material in culture.913

The isolation and characterization of limbal stem cells should facilitate optimal enhancement of precursor cell proliferation during epithelial cell expansion in culture, thereby increasing the reconstructive capacity of expanded cell patches and reducing the size of the initial donor cell pool that needs to be excised from a healthy contralateral eye for successful transplantation.

During the past decade, putative stem cells have been isolated from multiple organs by using an established correlation between stemness and the ability to efflux large aromatic compounds, in particular Hoechst 33342 by the ABCG2/BCRP transporter.1418 Hoechst 33342-transporting cells are easily recognized and sorted by flow cytometry on the basis of their fluorescent emission characteristics and are presently referred as side population (SP) cells. Several laboratories, including ours, have isolated SP cells from mammalian CNJ and limbal epithelia and subjected these to stemness tests.1925 Using in vivo BrdU labeling and long-term (>2 months) chasing in young rabbits, we showed that in the epithelia (E) of both limbus and conjunctiva (CNJ), SP cells were highly enriched in cells that have been slow cycling in vivo and displayed other features associated with stem cells in vitro.21 In most organs examined, SP cells are quite rare; they typically amount to between 0.05% and 0.5% of the total cells in each studied tissue or organ. These low numbers are consistent with the experiments in bone marrow indicating that SP cells are between the most primitive, or basic, stem cell.26

We have recently completed a microarray-based study of differentially expressed genes, or molecular signature, of SP cells isolated from the human CNJE.27 The rarity of these cells poses unique challenges and their investigation, such as in differential gene expression studies, has been problematic. Preparations of cadaveric donor corneas typically yielded 250,000 to 500,000 limbal cells and fewer than 1,000 SP cells. Hence, with a view to delineating a molecular signature of limbal SP cells, we opted in this study to employ specimens obtained from the pig, a species from which we could simultaneously obtain a large number of fresh corneas from young animals and for which microarrays are commercially available.

Many of the genes differentially expressed in the limbal SP cells may underpin functional features that underscore the SP cell elsewhere and/or universal features of stem cells, such as in vivo slow cycling. Hence, to probe for the possible general relevance of genes that are differentially expressed genes in the limbal SP cells, we performed similar microarray measurements for the pig CNJE, a tissue that shares with the corneal epithelium a common environment, developmental origin, and PAX6 expression28,29 and used these results and recently obtained microarray data for the human CNJE27 for comparative assessments.


Tissue Procurement, Cell Isolation, and Culture

Fresh pig eyes with intact eyelids enucleated from killed 3-month-old pigs were delivered within 30 hours of excision by Pel Freeze (Rogers, AR). A human cornea from an unidentifiable adult Caucasian male was obtained from the National Disease Research Interchange (NDRI, Philadelphia, PA).

Pig corneas were excised and, after careful removal of all perilimbal conjunctival tissue, limbal strips were microdissected into quarter sections. The limbal strips and quartered conjunctivas were incubated for 16 to 20 hours at 4°C in 5 mg/mL Dispase (Roche, Nutley, NJ) dissolved in 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid-buffered (hb), penicillin-streptomycin–complemented Dulbecco’s modified Eagle’s medium and Ham’s F12 1:1 mix (D/F-12), with a side-to-side tilting motion. Most of the limbal epithelial sheets spontaneously fully or partially separated from the underlying stroma and are easily scooped out of the enzymatic solution. Conjunctival sheets were separated from the subepithelial matrix by gentle mechanical prodding. Epithelial sheets were trypsinized for 20 to 25 minutes at 37°C with orbital agitation. The dissociated cells were sieved through 100- and 40-µm filters, pelleted, resuspended in bicarbonate-buffered D/F12 complemented with 5% FBS, 10 ng/mL cholera toxin, 10 ng/mL epidermal growth factor (EGF), insulin, transferrin, selenium, and 0.5% dimethyl sulfoxide (SHEM), plated in 75-cm2 flasks at a density of 80,000 to 100,000 cells/cm2 and cultured for 16 hours at 37°C in a 5% CO2 incubator. The three limbal experiments performed used 30 corneas each and the duplicate conjunctival experiments were based on pools of four tissues. Unless stated otherwise, all reagents were purchased from Sigma (St. Louis, MO).

Flow Cytometry

After overnight culture, the medium was refreshed with prewarmed solution to remove floating cells (60%–70% of the total plated) and complemented with 5 µg/mL Hoechst 33342 for 1.5 hours. The treated cells were then released by quick trypsinization, spun down, and resuspended in ice-cold phenol red-free D/F12 complemented with 4% serum and 1 µg/mL propidium iodide (PI). Cell sorting was performed with a flow cytometer (INFLUX; Cytopeia, Seattle, WA). The source tube was maintained at 4°C and SP and non-SP (nSP) cells were sorted directly into 750 µL RNA isolation solution (Tri-Reagent LS; Molecular Research Center [MRC] Cincinnati, OH). We used a specific range of forward (FSC) and side (SSC) light-scattering levels to exclude the great majority of lymphocytes (cells with very low SSC) and to limit the amount of other complex nonepithelial cells present within the epithelial strata (cells with very high SSC such as melanocytes and dendritic cells).27 Fumitremorgin C (FTC; a generous gift from Susan Bates, National Cancer Institute, Bethesda, MD) was used to determine the involvement of ABCG2 in the observed SP cells.22

Microarray Processing

RNA isolation solutions containing the collected SP and nSP cells were adjusted to 1.0 mL with the addition of H2O and 1 µL poly acryl carrier (MRC). RNA was isolated according to the manufacturer’s instructions and dissolved in 5-µL water. RNA yields, determined using a quantitation reagent (RiboGreen; Molecular Probes, Eugene, OR), ranged from 2.0 to 2.4 pg RNA/cell. RNA integrity was evaluated with microarrays (2100 Bioanalyzer 6000 nanochips; Agilent Technologies, Santa Clara, CA). cDNA was synthesized with reverse transcriptase (SuperScript Choice; Invitrogen, Rockville, MD) and a modified oligo dT primer (Affymetrix, Santa Clara, CA). The initial cDNA was subjected to two cycles of in vitro transcription (IVT; ENZO BioArray HighYield Kit; Affymetrix) to generate amplified biotin-labeled cRNA. After electro-phoretic assessment of the suitability of this cRNA product for microarray processing, appropriately fragmented (restricted) cRNA was hybridized to the gene microarray (GeneChip Porcine Genome Array; probe sets interrogate approximately 23,256 transcripts; Affymetrix). Hybridized microarrays were stained with a streptavidin-phycoerythrin reagent and fluorescence images were captured with a laser scanner (G2500A; Agilent).

Data Analyses

Normalized signal intensities (SIs; relative measures of gene expression) and signal quality properties in the form of gene present (P) or not-present calls were extracted from the fluorescent images (Microarray Suite, ver. 5.0 [MAS 5.0]; Affymetrix) and collected in spreadsheet format (lists) for downstream processing. To identify transcripts that are differentially expressed for the SP/nSP comparison, we first derived an expressed transcript spreadsheet incorporating only those transcripts that received the MAS5 present rating in all three experiments in at least one cell type: SP or nSP. The Hierarchical Clustering Explorer ( provided in the public domain by the Human Computer Interaction Lab, University of Maryland, College Park, MD) were used to generate a heat map for these transcripts. From the list of expressed genes we derived two primary differential lists, one incorporating only those transcripts for which the SP/nSP expression (SI) ratio (R) is greater than 1.5 in each of the three SP/nSP comparisons conducted and a second one limited to those transcripts for which the same ratios were less than 1:1.5. Transcripts included in these primary lists were considered over- or underexpressed if the overall SP/nSP expression ratio was greater than 2.0 or smaller than 0.5, respectively. The combination of these two latter lists is taken to represent a SP molecular signature. The duplicate experiments from pig conjunctiva and four replicates of human conjunctival epithelium27 were processed in an identical manner.

Pig microarray transcript-to-gene conversions were based on multispecies gene homology annotations generously provided by Christopher K. Tuggle (Iowa State University, Ames, IA; Couture et al., manuscript in preparation). A shortened description of the annotation methodology is provided in Wang et al.30 The list was further enriched by adding annotations (from the NetAffx Analysis Center; Affymetrix; that were not included in the data from Iowa State University.

Annotations indicating moderate to strong homology to known genes from other species were used to derive HGNC (HUGO Gene Nomenclature Committee) gene symbols for the primary differential lists. The HGNC symbols lists were submitted to The Database for Annotation, Visualization and Integrated Discovery (DAVID; provided in the public domain by the National Institute of Allergy and Infectious Diseases (NIAID), Bethesda, MD) to identify overrepresented biological or molecular processes within the differentially expressed transcripts. DAVID analysis probes each gene list against a relevant population list, and calculates scores as probabilities and/or false-discovery rates (FDRs) reporting on the likelihood of overrepresentation within the ontological systems derived from the Gene Ontology (GO) consortia and other public genomic resources.

Real-Time PCR

Total RNA was prepared from 23,500 sorted pig limbal SP and nSP cells. Three quarters of this RNA was reverse transcribed (Omniscript enzyme; RT+ product; Qiagen). The remaining RNA was processed in parallel under identical conditions but with omission of the enzyme, to generate an RT product. Real-time PCR (40 synthesis cycles at 60°C) of the RT+ (triplicates, each using one twelfth of the total RT reaction) and RT (single tube using one fourth of the reaction) products were performed in a sequence-detection system (Prism 7900HT; Applied Biosystems, Foster City, CA; using the RT2 SYBR Green/ROX PCR Master Mix; SuperArrays, Frederick, MD), with PCR forward and reverse primers for pig β-actin (aagtactccgtgtggatcgg and cagtccgccta-gaagcattt; yields 131-bp amplicon), CXCR4 (ctgagaagcatgacggacaa and tcccaaagtaccagtttgcc; 113-bp amplicon), GPX2 (aacagcctcaagtacgtccg and gtcgtcataagggtagggca; 133-bp amplicon), and DUSP5 (ggttgaagcag-caagatggt and gagaggtgcaacgagaaagg; 113 bp amplicon). SP/nSP message ratios in the amplified material were calculated from cycle thresholds (Cts) according to the ΔΔCt method, by using β-actin Cts for normalization.27


Cultures in the printed glass slides were stained with FITC-conjugated anti-keratin14 (clone LL002; Abcam, Cambridge, MA). Cryosections of pig and human limbocorneal sections were stained using an antibody against aquaporin 3 (a generous gift from Dennis Brown, Harvard Medical School, Cambridge, MA).


SPs in Briefly Cultured Pig Limbal and Conjunctival Cells

After an overnight culture in SHEM between 20% and 25% of pig limbal cells harvested by trypsinization firmly adhered to the plastic substratum. Incubation of the attached cells with Hoechst 33342 led to well-defined side populations for both limbal and conjunctival cells (Fig. 1). Those side populations extended from the main Hoechst saturated G0/G1 group down to cells that fully excluded the dye. This degree of Hoechst exclusion is higher than we observed in earlier studies with cells incubated in suspension soon after trypsinization.22 This improvement may reflect a better metabolic status of the cells, on which Hoechst efflux is dependent, after recovery in culture. Limbal SP cells accounted for 0.89% (n = 5) of the live gated cells. As previously shown for human or rabbit ocular surface cells, the SP feature was largely abolished when the incubation included the ABCG2-specific inhibitor FTC.

Side populations in pig ocular surface cells. (A) Limbus. (B) Conjunctiva. SP and nSP cell regions are indicated in (A). Bottom insets: images indicate that cells were collected only from limited side (SSC)- and forward (FSC)-scatter ranges. Top insets ...

RNA and Microarray Quality Controls

The need to trypsinize, centrifuge, resuspend, and maintain the epithelial cells in cold storage for 20 to 40 minutes before sorting and lysis resulted in moderate DNA damage as assessed by biochip (2100 Bioanalyzer; Agilent) (data not shown). Although the damage was within acceptable limits, careful retrospective analysis of data integrity was warranted. In addition to RNA damage before amplification, normal early termination of the processivity of the reverse transcription and IVT enzymes invariably contribute to a gradual decrease in the concentration of available probe as a function of the probed region distance from the priming (polyA) 3′ end. If these factors occur unevenly across the samples, the SI ratio could display a dependence on probe location relative to that priming end.

The results for three actin transcripts, a housekeeping gene, one near the poly(A) end, one in the middle of the gene and one near the 5′ end of the ~1000-bp gene, showed that, whereas SI values decreased along the 3′ to 5′ direction, the SP/nSP expression ratio remained essentially unchanged, whether we used them as the ratio between the SP and nSP averages, or as the average of SP/nSP ratios for each individual paired experiment (Table 1). Ratios for transcripts of differentially expressed genes were also not demonstrably affected by the SI values of different representative transcripts, even when these values differed by as much as 50-fold (Fig. 2A). Finally, heat maps of expressed transcripts (Fig. 2B) showcased the similarity between nSP or SP lists and the difference between the two cell types.

Quality controls for microarray processing. Right: SP/nSP SI ratios calculated from different transcripts of selected overexpressed genes. Gene names are indicated on the x-axis. Note that calculated expression ratios are not dependent on the signal intensity ...
Relationship between Probe Locations Respective to the PolyA Site and Signal Intensities and SP/nSP SI Ratios for the β-Acting Gene

Differential Expression

The total transcript-present list included 11,300 probes of which, 518 and 447 were SP over- and underexpressed, respectively. By removing transcripts that presently cannot be linked to a known gene and based on a redundancy gene representation (i.e., more than one transcript probe per gene) of approximately 20%, we calculated that the total gene-present list, when fully curated, would include approximately 9000 distinct genes. Likewise, limiting the differentially expressed lists to transcripts that can be ascribed to a known (named) gene produced over- and underexpressed gene lists containing 382 and 296 genes, respectively, for a total of 678 genes or 7.5% of the total identified. This large percentile underscores the unique phenotype of limbal epithelial SP cells, as previously observed for the human conjunctival epithelium.27 The 40 genes with the greatest over- or underexpression ratios and average SI greater than 100 (we assumed that most low SI genes are less likely to be relevant for cell phenotype and function), are displayed in Table 2A and Table 2B and the complete differentially expressed gene list is provided in the Appendix, online at

SP Genes Expressed in the Pig Limbal Epithelium

The limbal data included in Table 2A were compared with data obtained in a similar fashion from two replicated SP/nSP microarray analyses of pig conjunctiva epithelium and with the data earlier data from human conjunctiva.27 The great majority of the 40 genes included in Table 2A were also overexpressed in SP cells of pig CNJ; only 8 did not comply with the overexpression criteria (Table 2A, RpCNJ column). Extending the correlation to human CNJ (Table 2A, RhCNJ column) provided a convenient means of reinforcing the identification of genes that may have an intrinsic relationship to either stem cells in general, epithelial stem cells at large, or to stem cells of the PAX-6 positive epithelial of the ocular surface. From Table 2 A it is possible to identify genes that may be overexpressed in stem cells of the ocular surface in either a species-independent (e.g., CXCR4, MXD1, GPX2, EPAS1, MGST1, GRTP1, LNX1, VNN2, and FBXO32) or species-dependent (e.g., RHBDD2, AQP3) manner, and others that may be unique to the corneal epithelium (e.g., ZAR1, EEA1).

Overall, gene underexpression was less pronounced than overexpression. No gene was strongly underexpressed, and whereas there were 104 genes with an SP/nSP ratio greater than 3, only 20 underexpressed genes possessed nSP/SP ratios greater than this substantial threshold.


Due to limited amounts of SP RNA available per experiment, validation by real-time PCR was performed only for three markedly overexpressed genes. The PCR-derived SP/nSP ratios for CXCR4, GPX2, and DUSP5 were 24-, 12- and 3.7- fold, respectively, higher in all three cases than the averages measured in the microarray experiment (Table 2A, lines 1, 4; and Table 3, respectively). Agarose gel electrophoresis of the real-time PCR reaction products confirmed that for all four genes used, a single product of the expected size was observed (see the Methods section).

Selected Categories of Genes Differentially Expressed in the Pig Limbal SP Cells

The scarcity of antibodies known to recognize pig protein homologues is a serious impediment for validation at the polypeptide level. Aquaporin 3 is highly overexpressed in the pig but not in the human (Table 2A). We identified an antibody that recognizes the pig and human proteins. Consistent with what would be expected for a stem cell biomarker, in the pig the antibody generated intense staining in isolated basal cells (Fig. 3A). In stark contrast, in human samples, the protein was strongly expressed in all basal and early suprabasal limbal (Fig. 3B) and conjunctival (not shown) cells.

Expression of aquaporin 3 in the limbus, In pig, strong expression occurs only in isolated basal cells. In human the expression occurs throughout the basal and early suprabasal cells.

DAVID Cluster and Gene Categories Analysis

Differentially expressed genes with SI averages greater than 100 were subjected to DAVID’s functional annotation clustering using DAVID’s human gene set as background. For the overexpressed genes, the largest annotation clusters were those associated with regulation of cellular and biological developmental processes (59 and 52 genes of the 215 recognized; enrichment scores of 2.68 and 2.33, respectively). The GO_TERMS with the most significant overrepresentation with FDRs under 2% are shown in Table 4. It is noteworthy that the most overrepresented terms are related to regulatory inhibition of cellular biological and metabolic activity. It also identifies transcription from the RNA polymerase II promoter (genes, EPAS1, MAFB, CHD4, SAP18, EGR1, LMCD1, ID1, KLF7, ID2, TNFRSF1A, CDKN1C, HES1, and MYD88). Conversely, DAVID’s overrepresentation analysis of the underexpressed gene set essentially showed that the SP cells were underactive in terms of the synthesis of ribosomes and cellular organelles, a feature consistent with an undifferentiated, generally quiescent cell.

Over- and Underrepresented GO Terms within the Differentially Expressed Genes

Topical Gene Analysis

Even though the differentially expressed gene list is quite diverse and includes many genes whose functions remain poorly defined, several features of expression were compatible with the concept that these SP cells form part of the limbal stem cell population (Table 3). The SP cells overexpresses CXCR4, a pleiotropic cytokine receptor which plays a variety of roles, such as being a critical component of hematopoietic stem cells homing in bone marrow,31,32 the homing of mesen-chymal stem cells to bone33 and epithelial cell proliferation.34 ALDH1A1 is a cytosolic enzyme that has been shown to be generally associated with stem cells in multiple organs or tissues, and its overexpression may increase overall oxidative protection.35 Overexpression of keratin-15 and -19 and reduced CJA1/connexin43 is a well-documented feature of hair follicle stem cells36 and limbal stem cells,37 respectively. The differential expression of several categories of genes involved in the control of MAPK phosphorylation, cell cycling, stem cell replication and survival mechanisms (Table 3) need to be considered.


SP cell populations are quite small, typically comprising between 0.05% and 0.5% of the total organ/tissue cell count. Sound microarray studies require, at minimum, nanogram quantities of RNA, amounts that are difficult to derive from either human or rodent sources. Given this context, to delineate the molecular signature of limbal SP cells we resorted to the pig, a species possessing appropriate ocular dimensions and abundantly accessible tissue for which commercial microarrays are available. Although sequencing of the pig genome is far from complete and the available microarray has a limited gene representation span (e.g., neither p63, a gene commonly included in limbal stem cell studies or ABCG2, the prospective SP-generating transporter is represented), our studies yielded a robust molecular signature that incorporates nearly 700 genes of the ~9000 recognized as limbal-expressed by the microarray. These results clearly show the unique nature of the SP cells.

Expression patterns that are consistent with known stem cell features can be gleaned from the overall differentially expressed gene list. Slow cycling and the implicit low responsiveness to growth stimulation under steady state conditions represent the most recognized features of epithelial stem cells. The differentially expressed gene list incorporates several expression patterns that will confer or may contribute to slow cycling or quiescence (Fig. 4). Multiple dual specificity phosphatases (DUSPs), the primary inhibitors of ERK, p38, andJNK, the terminal mitogen-activated protein kinases38 are overex-pressed (Table 3), and none is underexpressed (Appendix). Since ERK1/2, mediate growth factor-initiated activation of nuclear transcription factors that control the G1-to-S restriction (R) point, overexpression of DUSPs, in particular DUSP5, with its capacity for ERK dephosphorylation and nuclear sequestra-tion,39 can be predicted to have an inhibitory effect on the stimulation of proliferation by circulating mitogens. We have recently confirmed proliferative inhibition by this enzyme in a corneal cell line (Wolosin et al., unpublished data, 2009). DUSP14 has been directly shown to slow pancreatic β-cell proliferation.40

Schematic description of nuclear protein gene patterns that may affect proliferation of SP cells. Overexpression of multiple nuclear dual specificity phosphatases may simultaneously decrease proliferative and migratory responses to growth factors in the ...

The overexpression of the MAD proteins MXD1 and MXI1, through its effect on MYC, can also have a negative impact on proliferation (Table 3, Fig. 4). MAPK-activated MYC promotes transcriptional activities that facilitate the G1-to-S transition and affect multiple cell cycle stages.41 MYC transcriptional effects require its association with MYC associated factor X (MAX).42 MAX availability, however, is negatively controlled by its alternative association with MAD proteins. Hence, markedly augmented expression of MXD1/MAD1 and MXI1/MAD2 levels in the SP cells is likely to further dull the responsiveness of these cells to mitogen stimuli.

The third notable expression pattern with the potential to inhibit the rate of G1 to S transition in a mitogen-independent manner involves cyclin D, the protein that integrates signals controlling G1-to-S transition via retinoblastoma. SP cells overexpress the two main kinase inhibitors of this cyclin, p27/CIP1, and p57/KIP2, and conversely moderately underexpress cyclin D2 and CDK kinases that activate cell cycle progression after association with nonphosphorylated cyclin D isoforms (Table 3, Fig. 4).

HES1 and ID are dominant negative blockers for several differentiation-inducing helix-loop-helix transcription factors. Recently Nakamura et al.43 have shown that HES1 is a critical survival factor for the mouse limbal stem cell, in similarity to the effect of both HES1 and ID1 and in the hematopoietic system.44,45

In addition to slow cycling, geno- and cytoprotection are two functions that seem naturally associated with long-lived stem cells. The overexpressed gene list includes a large number of genes that have been shown to be involved in defense against free radicals, oxidative stress, and removal of aromatic DNA damage inducers or xenobiotic compounds (Table 3).

The potential significance of the results to human biology is supported by the results displayed in Table 2. More than 70% of the genes overexpressed in the limbal SP cells are also overexpressed in the conjunctival SP cells, and likewise most of the genes overexpressed in these two systems are also overexpressed in the human conjunctiva SP cells. However, as suggested by the case of aquaporin 3, whose overexpression appears to be a pig-specific phenomenon (Table 2, Fig. 3), such extrapolations necessitate a degree of caution and/or specific reconfirmation.

From a practical clinical perspective, with the proviso that further study solidifies the notion that SP cells are part of the limbal stem cell population, many of the features indicated by the gene expression patterns observed in these cells may prevent prompt expansion of stem and early precursor cells within total cultured populations for use in surgical limbal reconstructive procedures. It is the transient amplifying cells that are more likely to undergo rapid gains in these cultures. Modification of ex vivo cell behavior through either pharmacologic means or reversible ectopic gene or antisense expression may be useful to improve the fraction of stem/precursor cells within expanded limbal populations.

In summary, notwithstanding the present limitations imposed by incomplete pig genome sequence data, the differential global gene expression analysis outlined in this report demonstrated the feasibility and value of the microarray studies used in advancing limbal stem cell research.


The authors thank the Mount Sinai School of Medicine for the use of the Flow Cytometry, Microarray, and Real Time PCR shared facilities, without which the study could not have been completed.

Supported by National Institutes of Health Grants EY 014878, EY 015132, and EY01867 (Core Center Grant) and an Unrestricted grant from Research to Prevent Blindness, Inc.


Disclosure: M.A.M. Akinci, None; H. Turner, None; M. Taveras, None; J.M. Wolosin, None


1. Davanger M, Evensen A. Role of the pericorneal papillary structure in renewal of corneal epithelium. Nature. 1971;229:560–561. [PubMed]
2. Schermer A, Galvin S, Sun TT. Differentiation-related expression of a major 64 K corneal keratin in vivo and in culture suggests limbal location of corneal epithelial stem cells. J Cell Biol. 1986;103:49–62. [PMC free article] [PubMed]
3. Cotsarelis G, Cheng SZ, Dong G, Sun TT, Lavker RM. Existence of slow-cycling limbal epithelial basal cells that can be preferentially stimulated to proliferate: implications on epithelial stem cells. Cell. 1989;57:201–209. [PubMed]
4. Tseng SC. Concept and application of limbal stem cells. Eye. 1989;3:141–157. [PubMed]
5. Wolosin JM, Xiong X, Schutte M, Stegman Z, Tieng A. Stem cells and differentiation stages in the limbo-corneal epithelium. Prog Retin Eye Res. 2005;19:223–255. [PubMed]
6. Chen JJ, Tseng SC. Corneal epithelial wound healing in partial limbal deficiency. Invest Ophthalmol Vis Sci. 1990;31:1301–1314. [PubMed]
7. Puangsricharern V, Tseng SC. Cytologic evidence of corneal diseases with limbal stem cell deficiency. Ophthalmology. 1995;102:1476–1485. [PubMed]
8. Huang AJ, Tseng SC. Corneal epithelial wound healing in the absence of limbal epithelium. Invest Ophthalmol Vis Sci. 1991;32:96–105. [PubMed]
9. Kenyon KR, Tseng SC. Limbal autograft transplantation for ocular surface disorders. Ophthalmology. 1989;96:709–722. [PubMed]
10. Tsai RJ, Tseng SC. Human allograft limbal transplantation for corneal surface reconstruction. Cornea. 1994;13:389–400. [PubMed]
11. Tsubota K, Satake Y, Ohyama M, et al. Surgical reconstruction of the ocular surface in advanced ocular cicatricial pemphigoid and Stevens-Johnson syndrome. Am J Ophthalmol. 1996;122:38–52. [PubMed]
12. Pellegrini G, Traverso CE, Franzi AT, Zingirian M, Cancedda R, De Luca M. Long-term restoration of damaged corneal surfaces with autologous cultivated corneal epithelium. Lancet. 1996;349:990–993. [PubMed]
13. Hayashida Y, Nishida K, Yamato M, et al. Ocular surface reconstruction using autologous rabbit oral mucosal epithelial sheets fabricated ex vivo on a temperature-responsive culture surface. Invest Ophthalmol Vis Sci. 2005;46:1632–1639. [PubMed]
14. Goodell MA, Brose K, Paradis G, Conner AS, Mulligan RC. Isolation and functional properties of murine hematopoietic stem cells that are replicating in vivo. J Exp Med. 1996;183:1797–1806. [PMC free article] [PubMed]
15. Zhou S, Schuetz JD, Bunting KD, et al. The ABC transporter Bcrp1/ABCG2 is expressed in a wide variety of stem cells and is a molecular determinant of the side-population phenotype. Nat Med. 20012;7:1028–1034. [PubMed]
16. Zhou S, Morris JJ, Barnes Y, Lan L, Schuetz JD, Sorrentino BP. Bcrp1 gene expression is required for normal numbers of side population stem cells in mice, and confers relative protection to mitoxantrone in hematopoietic cells in vivo. Proc Natl Acad Sci U S A. 2002;99:12339–12344. [PubMed]
17. Ono M, Maruyama T, Masuda H, et al. Side population in human uterine myometrium displays phenotypic and functional characteristics of myometrial stem cells. Proc Natl Acad Sci U S A. 2007;104:18700–18705. [PubMed]
18. Hoshi N, Kusakabe T, Taylor BJ, Kimura S. Side population cells in the mouse thyroid exhibit stem/progenitor cell-like characteristics. Endocrinology. 2007;148:4251–4158. [PMC free article] [PubMed]
19. Wolosin JM, Tseng SCG, Budak MT, et al. The limbo-corneal epithelium and its stem cells: principles and applications. In: Greer VE, editor. Focus in Stem Cell Research. New York: Nova Science Publishers; 2004. pp. 59–92.
20. Paiva CS, Chen Z, Corrales RM, Pflugfelder SC, Li DQ. ABCG2 transporter identifies a population of clonogenic human limbal epithelial cells. Stem Cells. 2005;23:63–73. [PMC free article] [PubMed]
21. Budak MT, Alpdogan OS, Zhou M, Lavker RM, Akinci MA, Wolosin JM. Ocular surface epithelia contain ABCG2-dependent side population cells exhibiting features associated with stem cells. J Cell Sci. 2005;118:1715–1724. [PMC free article] [PubMed]
22. Epstein SP, Wolosin JM, Asbell PA. P63 expression levels in side population and low light scattering ocular surface epithelial cells. Trans Am Ophthalmol Soc. 2005;103:187–199. [PMC free article] [PubMed]
23. Wolosin JM. Cell markers and the side population phenotype in ocular surface epithelial stem cell characterization and isolation. Ocul Surf. 2006;4:10–23. [PubMed]
24. Umemoto T, Yamato M, Nishida K, Yang J, Tano Y, Okano T. Limbal epithelial side-population cells have stem cell-like properties, including quiescent state. Stem Cells. 2006;24:86–94. [PubMed]
25. Park KS, Lim CH, Min BM, et al. The side population cells in the rabbit limbus sensitively increased in response to the central cornea wounding. Invest Ophthalmol Vis Sci. 2006;47:892–900. [PubMed]
26. Goodell MA, Rosenzweig M, Kim H, et al. Dye efflux studies suggest that hematopoietic stem cells expressing low or undetectable levels of CD34 antigen exist in multiple species. Nat Med. 1997;3:1337–1345. [PubMed]
27. Akinci MAM, Turner H, Taveras M, et al. Molecular profiling of the of conjunctival epithelial stem side population cells identifies atypical cell surface markers and sources of a slow cycling phenotype. Invest Ophthalmol Vis Sci. Published online March 25, 2009. [PMC free article] [PubMed]
28. Davis JA, Reed RR. Role of Olf-1 and Pax-6 transcription factors in neurodevelopment. J Neurosci. 1996;16:5082–5094. [PubMed]
29. Koroma BM, Yang JM, Sundin OH. The Pax-6 homeobox gene is expressed throughout the corneal and conjunctival epithelia. Invest Ophthalmol Vis Sci. 1997;38:108–120. [PubMed]
30. Wang Y, Couture OP, Qu L, et al. Analysis of porcine transcriptional response to Salmonella enterica serovar Choleraesuis suggests novel targets of NFkappaB are activated in the mesenteric lymph node. BMC Genomics. 2008;9:437. [PMC free article] [PubMed]
31. Ratajczak MZ, Zuba-Surma E, Kucia M, Reca R, Wojakowski W, Ratajczak J. The pleiotropic effects of the SDF-1-CXCR4 axis in organogenesis, regeneration and tumorigenesis. Leukemia. 2006;20:1915–1924. [PubMed]
32. Chute JP. Stem cell homing. Curr Opin Hematol. 2006;13:399–406. [PubMed]
33. Kitaori T, Ito H, Schwarz EM, et al. Stromal cell-derived factor 1/CXCR4 signaling is critical for the recruitment of mesenchymal stem cells to the fracture site during skeletal repair in a mouse model. Arthritis Rheum. 2009;60:813–823. [PubMed]
34. Florin L, Maas-Szabowski N, Werner S, Szabowski A, Angel P. Increased keratinocyte proliferation by JUN-dependent expression of PTN and SDF-1 in fibroblasts. J Cell Sci. 2005;118:1981–1989. [PubMed]
35. Moreb JS. Aldehyde dehydrogenase as a marker for stem cells. Curr Stem Cell Res Ther. 2008;3:237–246. [PubMed]
36. Kloepper JE, Tiede S, Brinckmann J, et al. Immunophenotyping of the human bulge region: the quest to define useful in situ markers for human epithelial hair follicle stem cells and their niche. Exp Dermatol. 2008;7:592–609. [PubMed]
37. Matic M, Petrov IN, Chen S, Wang C, Dimitrijevich SD, Wolosin JM. Stem cells of the corneal epithelium lack connexins and metabolite transfer capacity. Differentiation. 1997;61:251–260. [PubMed]
38. Patterson KI, Brummer T, O’Brien PM, Daly RJ. Dual-specificity phosphatases: critical regulators with diverse cellular targets. Biochem J. 2009;418:475–489. [PubMed]
39. Mandl M, Slack DN, Keyse SM. Specific inactivation and nuclear anchoring of extracellular signal-regulated kinase 2 by the inducible dual-specificity protein phosphatase DUSP5. Mol Cell Biol. 2005;25:1830–1845. [PMC free article] [PubMed]
40. Klinger S, Poussin C, Debril MB, Dolci W, Halban PA, Thorens B. Increasing GLP-1-induced beta-cell proliferation by silencing the negative regulators of signaling cAMP response element modulator-alpha and DUSP14. Diabetes. 2008;57:584–593. [PubMed]
41. Mateyak MK, Obaya AJ, Sedivy JM. c-Myc regulates Cyclin D-Cdk4 and -Cdk6 activity but affects cell cycle progression at multiple independent points. Mol Cell Biol. 1999;19:4672–4683. [PMC free article] [PubMed]
42. Grandori C, Cowley SM, James LP, Eisenman RN. The Myc/Max/Mad network and the transcriptional control of cell behavior. Annu Rev Cell Dev Biol. 2000;16:653–699. [PubMed]
43. Nakamura T, Ohtsuka T, Sekiyama E, et al. Hes1 regulates corneal development and the function of corneal epithelial stem/progenitor cells. Stem Cells. 2008;26:1265–1274. [PubMed]
44. Kunisato A, Chiba S, Nakagami-Yamaguchi E, et al. HES-1 preserves purified hematopoietic stem cells ex vivo and accumulates side population cells in vivo. Blood. 2003;101:1777–1783. [PubMed]
45. Perry SS, Zhao Y, Nie L, Cochrane SW, Huang Z, Sun XH. Id1, but not Id3, directs long-term repopulating hematopoietic stem-cell maintenance. Blood. 2007;110:2351–2360. [PubMed]