Donor age and long-term culture do not negatively influence the stem potential of limbal fibroblast-like stem cells
Stem Cell Research & Therapyvolume 7, Article number: 83 (2016)
The Erratum to this article has been published in Stem Cell Research & Therapy 2016 7:106
In regenerative medicine the maintenance of stem cell properties is of crucial importance. Ageing is considered a cause of reduced stemness capability. The limbus is a stem niche of easy access and harbors two stem cell populations: epithelial stem cells and fibroblast-like stem cells. Our aim was to investigate whether donor age and/or long-term culture have any influence on stem cell marker expression and the profiles in the fibroblast-like stem cell population.
Fibroblast-like stem cells were isolated and digested from 25 limbus samples of normal human corneo-scleral rings and long-term cultures were obtained. SSEA4 expression and sphere-forming capability were evaluated; cytofluorimetric assay was performed to detect the immunophenotypes HLA-DR, CD45, and CD34 and the principle stem cell markers ABCG2, OCT3/4, and NANOG. Molecular expression of the principal mesenchymal stem cell genes was investigated by real-time PCR. Two-dimensional gel electrophoresis and mass spectrometric sequencing were performed and a stable proteomic profile was identified. The proteins detected were explored by gene ontology and STRING analysis. The data were reported as means ± SD, compared by Student’s unpaired t test and considering p < 0.05 as statistically significant.
The isolated cells did not display any hematopoietic surface marker (CD34 and CD45) and HLA-DR and they maintained these features in long-term culture. The expression of the stemness genes and the multilineage differentiation under in-vitro culture conditions proved to be well maintained. Proteomic analysis revealed a fibroblast-like stem cell profile of 164 proteins with higher expression levels. Eighty of these showed stable expression levels and were involved in maintenance of “the stem gene profile”; 84 were differentially expressed and were involved in structural activity.
The fibroblast-like limbal stem cells confirmed that they are a robust source of adult stem cells and that they have good plasticity, good proliferative capability, and long-term maintenance of stem cell properties, independently of donor age and long-term culture conditions. Our findings confirm that limbal fibroblast-like stem cells are highly promising for application in regenerative medicine and that in-vitro culture steps do not influence their stem cell properties. Moreover, the proteomic data enrich our knowledge of fibroblast-like stem cells.
The limbus, located at the junction of the cornea and conjunctiva of the ocular surface, is characterized by stromal invaginations (the palisades of Vogt). These structures give anatomical and functional properties protecting stem cells from insults and allowing constant renewal of the corneal epithelium [1–8].
Two stem cell populations have been described in the limbus: limbal epithelial stem cells (LESCs) and limbal stromal stem cells. Their similar and different features have not been defined unequivocally [9–12]. Polisetty et al.  proposed a role for limbal mesenchymal stromal cells (MSCs) in maintaining support of the stem cell niche. More recent data have validated the hypothesis of “niche stromal cells” based on the capacity of limbal stromal cells to be an efficient feeder layer for ex-vivo limbal epithelial cell expansion [8, 14, 15]. However, given the complexity of the limbal niche structure and its cellular components, which cellular type has the main role in normal tissue maintenance remains to be seen [16–18].
We previously identified a subpopulation of limbal stem cells, which we referred to as fibroblast-like limbal stem cells (f-LSCs) . We indicated a core set of attributes that uniquely characterize f-LSCs and that classify them as mesenchymal stem cells. In support of this, f-LSCs expressed stem cell surface antigens SSEA4 (stage-specific embryonic antigen-4), TRA 1–60, and TRA 1–81 and several nuclear transcription factors, such as OCT4 (octamer-binding transcription factor 4), NANOG (Homeobox protein NANOG), and SOX2 (SRY (sex determining region Y)-box 2), involved in self-renewal and maintenance of pluripotency of both embryonic and adult stem cells [20, 21]. The f-LSCs were positive for the limbal stem cell marker ABCG2 (ATP-binding cassette sub-family G member 2) and were negative for the LESC marker ΔNp63 (a splice variant of p63) [22–26]. More recently it has been hypothesized that age negatively affects stem cell number and potential, also referring to the LESC population [27, 28].
The current study was undertaken to explore whether ex-vivo expansion of f-LSCs could modify their stem molecular features. Our study aimed to increase knowledge in the field of limbal mesenchymal stem cell research regarding three facets: firstly, we aimed to discover whether long-term culture affects the gene expression and proteomic profile of f-LSCs; secondly, we proposed a preserved “proteomic stem cell pattern” in ageing and long-term culture conditions; and finally, we evaluated the maintenance of multilineage differentiation capability in vitro.
Establishment of limbal cell cultures
Human tissues were used in accordance with the Declaration of Helsinki and informed written consent was given by each patient. The study was approved by the Ethical Committee of the AOUP, University of Palermo (No. 09/2009).
Normal human corneo-scleral rings from donors aged between 24 and 74 years were obtained 2–3 hours post surgery in the Ophthalmology Department (AOUP, University of Palermo, Italy). The rings were kept in Hank’s Balanced Salt Solution (HBSS; PAA, Pashing, Austria) and then cut into small segments to facilitate isolation of the limbus from the sclera.
Enzymatic digestion and culture
The limbal segments were incubated with collagenase I (5 mg/ml; Sigma-Aldrich, St. Louis, MO, USA) overnight at 37 °C in a shaking bath. The following day, the digest was placed in a p60 dish culture (Corning, New York, USA) with the fibroblastic-maintenance medium (DMEM/F12 supplemented with 10 % embryonic stem cell-tested fetal bovine serum (EC-FBS; PAA), 1× ITS (5 μg/ml insulin, 5 μg/ml transferrin, 5 μg/ml selenium; PAA), and 20 ng/ml basic fibroblast growth factor (b-FGF; Preprotech, London, UK)) until cells reached confluence, changing the medium when necessary.
Subsequently, the f-LSC subculture was kept in the expansion medium (f-EM: DMEM/F12 supplemented with 5 % EC-FBS (PAA), 1× ITS (PAA), and 4 ng/ml b-FGF (Preprotech)).
The f-LSCs were subcultured and seeded at a density of 4 × 103 cells/cm2 and cell counts were performed at 24, 48, 72, 96, and 120 hours by optical microscope observation after staining with vital trypan blue. The doubling time was calculated online (http://www.doubling-time.com/compute.php; Roth V. 2006). For each sample of donor limbus, three sets of experiments were used for calculations.
The f-LSCs were placed in ultralow-attachment six-well plates (Corning) at a density of no more than 1 × 102 cells/cm2, and were cultured in f-EM without serum. Sphere formation was assessed by counting the number of spheres (cells > 3) under an optical microscope.
The spheres were transferred into a cell-culture chamber slide (Labtek II; Nunc, Waltham, MA, USA) and incubated in f-EM, under adhesion-condition culture at 37 °C in 5 % CO2 to allow their attachment. After 2 hours, cells were washed with phosphate-buffered saline (PBS) and fixed for 30 minutes in 2 % (wt/vol.) paraformaldehyde in PBS at room temperature, then incubated with 1 μg/ml PE-conjugated anti-human SSEA-4 monoclonal antibody (Miltenyi Biotec, Bergisch Gladbach, Germany) in PBS/bovine serum albumin (BSA; PAA) for 30 minutes, at room temperature. After incubation the cells were washed three times with PBS, counterstained with DAPI (Sigma Aldrich), and observed under a Zeiss Axiophot fluorescence microscope equipped with a Nikon DS-FI1 CCD camera.
Human bone marrow-derived MSCs (BM-MSCs; Lonza, Walkersville, MD, USA) were used as a stem marker positive control, HeLa (kindly obtained from Dr Salvatore Feo at University of Palermo) were used as an LESC marker positive control, and unstained f-LSCs were used as a negative control.
Sorting assay through magnetic isolation
The f-LSCs were magnetically separated for SSEA4 expression by MACS MicroBead Technology (Miltenyi Biotec) according to the manufacturer’s instructions.
Flow cytometry analysis
The cells were harvested and filtered through a 40-μm filter mesh and suspended at a concentration of 1 × 106 cells/ml. Then 100 μl of cell suspension containing 5 × 105 cells was used for each flow cytometric test.
Human anti-HLA-DR, human anti-CD34, and human anti-CD45 monoclonal antibodies were tested on f-LSCs and were detected with the appropriate secondary antibody (Table 1). The incubation conditions were in accordance with the manufacturers’ instructions.
Stem cell phenotypes
The cells were double-stained with human anti-SSEA4 and human anti-ABCG2 monoclonal antibody, both surface MSC markers. The f-LSCs were tested for SSEA4 and for the human nuclear markers ΔNp63 or OCT4 or NANOG monoclonal antibody, after permeabilization with PBS supplemented with 0.1 % saponin and 1 % BSA for 20 minutes.
The antibody dilution, incubation, and detection conditions are presented in Table 1.
All reaction mixtures were then acquired with a FACS Calibur flow cytometer (Becton-Dickinson, Franklin Lakes, NJ, USA) and analyzed with the CellQuest Pro software. BM-MSCs were used as a positive control for SSEA4, NANOG, ABCG2, and OCT4, and HeLa cells were used as a positive control for ΔNp63.
Analysis of cell cycle status of MSCs
Single cell suspensions of f-LSCs were obtained from two different culture passages: P4 and P30. For DNA content analysis, Nicoletti’s protocol was performed. Briefly, 1 × 106 cells were fixed in 70 % ethanol, rehydrated in PBS, and then resuspended in a DNA extraction buffer (with 0.2 M NaHPO4, 0.1 % Tritonx-100, pH 7.8). After staining with 1 μg/ml of propidium iodide (PI) for 5 minutes, the fluorescence intensity was determined by analysis on a FACS Calibur flow cytometer (Becton-Dickinson). Data acquisition was performed with CellQuest software (Becton Dickinson), and the percentages of phase G1, S, and G2 cells were calculated with the MODFIT-LT software program (Verity Software House, Inc. Augusta, Topsham, ME, USA).
RNA extraction, quantification, and retrotranscription
Total RNA was extracted and purified using E.Z.N.A. Total RNA Kit I (Omega Bio-Tek Inc., Norcross, GA, USA) according to the manufacturer’s instructions. RNA quantity and quality were assessed by Nano Drop 2000 (Thermo Scientific, Waltham, MA, USA), and 2 μg of f-LSC total RNA was reverse-transcribed to cDNA in a volume of 20 μl with Oligo dT primers (Applied Biosystems, Carlsbad, CA, USA) and the Reverse Transcriptase Rnase kit (Improm II; Promega, Madison, WI, USA).
Real-time quantitative PCR
Real-time quantitative PCR primers were purchased from Qiagen (QuantiTect® Primer Assays; Qiagen, Milan, Italy) and Eurofin MWG (Biotech, Bergish Gladbach, Germany) and are listed in Table 2. All reactions were performed using the Quantitect SYBR Green PCR Kit (Qiagen, Valencia, CA, USA) on the RotorGene Q Instrument (Qiagen, Valencia, CA, USA). Each cDNA sample was mixed with specific primer sets and PCR master mix. The PCR parameters included denaturation at 95 °C for 3 minutes, then 40 cycles at 95 °C for 20 seconds, annealing at 60 °C for 30 seconds, and elongation at 72 °C for 60 seconds. Reactions were performed at least in triplicate. The specificity of the amplified products was determined by the melting peak analysis. The relative quantification model with efficiency correction was applied to normalize the expression of the target gene to β-actin (used as the housekeeping gene) and to compare gene expression with BM-MSCs (used as a positive cell control), on Rest2009 software (Relative Expression Software Tool; Qiagen, Valencia, CA, USA) . The results were represented as histograms on GraphPad software (GraphPad Software, Inc., La Jolla, CA, USA).
The f-LSCs were scraped and incubated on ice for 30 minutes with M-RIPA buffer (50 mM Tris, pH 7.5, 0.1 % Nonidet P-40, 0.1 % deoxycholate, 150 mM NaCl, 4 mM EDTA) and a mixture of protease inhibitors (0.01 % aprotinin, 10 mM sodium pyrophosphate, 2 mM sodium orthovanadate, 1 mM PMSF). The whole cellular lysate was centrifuged at 12,000 rcf for 8 minutes to clear cell debris and the supernatant was dialyzed against ultrapure distilled water, lyophilized, and stored at –80 °C, as described previously [30, 31]. The protein concentration in the cellular extracts was determined using the Quick Start™ Bradford Protein Assay (BIO RAD, Segrate, Milan, Italy) according to the manufacturer’s instructions.
Two-dimensional gel electrophoresis (2D-IPG) was performed as described previously [30, 31]. Briefly, protein samples of f-LSCs were solubilized, and aliquots of 45 μg (analytical gels) or 1.5 mg (preparative gels) of total proteins were separately applied for isoelectrofocusing (IEF) using commercial sigmoidal IPG strips, 18 cm long with pH range 3.0–10; BIO RAD). The focused proteins were then separated on 9–16 % linear gradient polyacrylamide gels (SDS-PAGE) and were visualized by means of ammoniacal silver staining.
For image acquisition and data analysis, silver-stained gels were digitized using a computing densitometer and analyzed with ImageMaster 2D PLATINUM software (Amersham, Little Chalfont, Buckinghamshire, UK). Gel calibration was carried out using an internal standard and the support of the ExPaSy molecular biology server; the quantitative analysis of protein spots was performed in duplicate maps, and normalized as vol. % (integration of optical density over the spot area). The differential expression of proteins was evaluated when the difference in their values was ≥ 3 % volume. The labels correspond to the access number of the Swiss-Prot/TrEMBL database.
Mass spectrometric sequencing was performed with the Voyager DE-PRO (Applied Biosystems, Carlsbad, CA USA) mass spectrometer . Peptide mass fingerprinting was compared with the theoretical masses from the Swiss-Prot or NCBI sequence databases using Mascot (http://www.matrixscience.com/).
Gene ontology and network analysis
The genes expressing the invariant and variant proteins were analyzed for their enrichment in specific gene ontology annotation using the Gene Ontology Consortium website (http://geneontology.org/) (Additional file 1: Table S2) .
Network analysis was performed on the genes expressing the invariant and variant proteins using the STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) website (http://string-db.org/) .
Multilineage potential differentiation assays
Multilineage differentiation capability of f-LSCs at early and late passages (P4 and P20) was assessed. The following differentiation protocols were performed.
For osteogenic differentiation, 5 × 103/cm2 f-LSC cells were cultured in home-made differentiation medium consisting of DMEM supplemented with 15 % FBS, 10–4 mM dexamethasone (Sigma-Aldrich), 10 mM glycerophosphate (Sigma-Aldrich), and 0.05 mM ascorbic acid (Sigma-Aldrich) . After 21 days of culture in the differentiation medium, cells were stained with Alizarin red S (Sigma-Aldrich) to detect the calcium deposits. Briefly, the medium was removed and the cells were fixed with 4 % formaldehyde solution for 30 minutes, and after fixation were rinsed twice with distilled water and stained with 2 % Alizarin red S (pH 4.2) for 3 minutes. After incubation the cells were observed under a light optical microscope at 20–40× magnification.
Adipose differentiation was induced by seeding 5 × 103/cm2 of f-LSCs and culturing in home-made differentiation medium consisting of DMEM medium containing 10 % FBS, 0.5 nM 1-methyl-3-isobutylxanthine (IBMX; Sigma-Aldrich), 10–4 mM dexamethasone, 10 μg/ml insulin, and 100 μM indomethacin (Sigma-Aldrich) . To detect the presence of lipid droplets, cells were fixed with 4 % formaldehyde solution, rinsed twice in distilled water, stained with Oil Red O (Sigma-Aldrich), and observed under a light microscope at 20–40× magnification.
Chondrogenic differentiation was induced by culturing cell mass in serum-free DMEM (PAA) with 10 ng/ml TGFβ-3 (Preprotech), 0.05 mM ascorbic acid (Sigma-Aldrich), 2 mM sodium pyruvate (PAA), and 10–7 mM dexamethasone (Sigma-Aldrich) in a six-well culture plate. Briefly, micromass cultures were generated by seeding 5-μl droplets of a cell solution of 1.6 × 107 viable cells/ml in f-EM. After 2 hours of incubation at 37 °C, culture medium was replaced with chondrogenesis differentiation medium. After 21 days of culture the cell mass was evaluated morphologically and stained with Alcian blue to search for sulfated proteoglycan deposits.
All assays were performed in triplicate. The data were reported as means ± SD and compared using the appropriate version of Student’s unpaired t test. Test results were reported as two-tailed p values, where p < 0.05 was considered statistically significant.
Isolation of f-LSCs
We determined that approximately 3 × 104–5 × 104 cultured f-LSCs were isolated from each explanted human corneo-scleral ring after enzymatic digestion. The digests were kept in culture in f-EM (we referred to these culture passages as P0) and in 2 weeks a single cell population with an elongated shape was obtained (Fig. 1A.a).
Within 10 days the cells, kept in a low-adhesion culture condition, gave rise to floating spherical cell bodies (considered a hallmark of the stemness feature) with a clear and well-delineated border, which we referred to as “limbospheres” (Fig. 1A.b). When transferred to adhesion conditions the limbospheres spread out and a monolayer culture of fibroblast-like cells formed (Fig. 1A.c,d).
SSEA4 expression was investigated by immunofluorescence analysis both in f-LSC spheres (Fig. 1B.a–c) and in monolayer f-LSCs (Fig. 1B.f–h). BM-MSCs and HeLa cells were respectively used as positive control (Fig. 1B.d–i) and negative control (Fig. 1B.e–l). As reported in the literature, it did not seem possible to obtain spheres from cells that have no SSEA4 expression.
Flow cytometry demonstrated that SSEA4 was highly expressed (88.95 ± 7.8 %) in f-LSC monolayer cultures and did not significantly increase after sorting (90.76 ± 5.6 %). Immunophenotypes were analyzed in the presorting (total) population and the postsorting SSEA4+ population. Both populations showed almost no expression of CD34 (0.32 ± 0.01 % vs. 0.19 ± 0.02 %, respectively), CD45 (1.35 ± 0.7 % vs. 0.65 ± 0.2 %, respectively), and HLA-DR (0.25 ± 0.04 % vs. 0.16 ± 0.07 %, respectively) (Fig. 1C).
Since neither population displayed any significant differences in marker expression, the subsequent cell analyses were only performed on the total monolayer f-LSC population.
Proliferation in long-term culture
The distribution of f-LSCs in the various stages of the cell cycle was investigated at two different passages. The early passage P4 showed an increase in the percentage of cells arrested in the G1 phase compared with the late-passage P30 (76.15 ± 2.45 % vs. 89.90 ± 3.65 %, p < 0.05) and a decrease in the percentage of cells in the S phase and in the G2 phase (Fig. 2a, upper panel). Moreover, in P30 the mRNA level of CDKN1B (cyclin-dependent kinase inhibitor 1B, also known as P37) was increased up to 4.59-fold relative to P3 (Fig. 2a, left-lower panel). The culture doubling time was slightly extended in P30 with respect to P3 f-LSCs, reaching about 46.2 ± 5.37 hours vs. 39.6 ± 2.54 hours (Fig. 2a, right-lower panel).
Flow cytometry stem cell phenotype characterization in long-term culture
Phenotype characterization was performed at two different culture passages (P3 and P30). A collection of double-positive populations (dpp) was obtained for several stem markers. The data expressions detected with flow cytometry analysis are the following (P3 vs. P30): ABCG2+/SSEA4+ = 98.6 ± 7.1 % vs. 96.6 ± 4.7 % (p = NS); OCT4+/SSEA4+ = 92.1 ± 6.4 % vs. 88.3 ± 2.6 % (p = NS); and NANOG+/SSEA4+ = 95.42 ± 6.8 % vs. 80.52 ± 4.2 % (p = NS). The f-LSCs did not express Δnp63 (Δnp63+/SSEA4+ < 0.5 %). Figure 2b shows one representative experiment out of a total of 25 sets of experiments.
The positive controls are shown in Additional file 2: Figure S1.
Stem-related gene expression is not influenced by age and by long-term cultures
We compared the f-LSC stem molecular expression pattern in young donors (age < 45 years: N = 11, six males and five females; Group A) and older donors (age > 45 years: N = 14, eight males and six females; Group B). No substantial differences were found in the immune phenotype (Fig. 3a): both Groups A and B were negative for CD34 (0.21 ± 0.07 % vs. 0.19 ± 0.1 %, p = NS), CD45 (0.27 ± 0.1 % vs. 0.13 ± 0.08 %, p = NS), and HLA-DR (0.23 ± 0.12 % vs. 0.16 ± 0.04 %, p = NS) whereas they highly expressed SSEA4 (92.3 ± 5.1 % vs. 90.76 ± 4.61 %, p = NS). The expression analysis performed on three different culture passages (P3, P20, and P30) showed no substantial difference for all stem cell markers tested between f-LSCs of patients within the two different groups (Figs. 3b,c). Comparing the gene expression levels of the two groups, it was clear that there was no significant age-related difference. The f-LSC stem cell markers, NANOG, THY-1, and SOX-2, were firmly maintained, as well as ABCG2, thus confirming maintenance of the limbal nature. Only OCT4 relative expression showed a faint decrease that appeared to be correlated with long-term culture but not with age. Moreover, CK12 (cytokeratin 12, corneal epithelial marker) and ΔNp63 (LESC marker) were not expressed (or only weakly expressed; data not shown), thus confirming that the culture conditions were not suitable for epithelial differentiation. Our findings showed that cultured and expanded f-LSCs exhibited similar expression patterns both in long-term culture and soon after isolation, maintaining their own expression characteristics.
The f-LSC proteomic profile is not influenced by age and by long-term cultures
To further characterize the f-LSC phenotype, we performed a proteomic analysis. Two-dimensional gel electrophoresis assay resulted in a master gel that revealed 164 spots, 78 % of which had an average pixel density of about 0.3 (Fig. 4a). The list of assigned proteins is shown in Additional file 3: Table S1. The derived protein f-LSC profiles represented by the densitometry, expressed as a graph, are shown in Fig. 4b, which also shows the relative intensities of the single protein expression levels. In Additional file 4: Figure S2, panels A and B show the proteomic maps of f-LSCs from the same older donor (Group B) at two different culture passages (P4 vs. P20), panels C and D show the proteomic maps of f-LSCs from two different young donors of the same age (Group A) at the same passage culture, and panel E shows the proteomic map of P20 f-LSCs from one 64-year-old donor (Group B). By comparing the f-LSC P4/B and f-LSC P20/B proteomic profiles, we detected overlapping for about 90 % of the proteins expressed (Additional file 5: Figure S3A); to test the occurrence of possible individual variants we compared the proteomic profile of f-LSCs from different donors of the same age (Group A) and we found overlapping for about 83 % of the proteins (Additional file 5: Figure S3B). Moreover, comparing the proteomic profiles of Group A with Group B donor at the same culture passage, the proteomic profiles overlapped by about 68 % (Additional file 5: Figure S3C).
By comparative analysis between the proteomic profiles of all investigated conditions, a protein group with higher expression levels was detected. Eighty of the proteins maintained invariant expression levels, whereas 84 proteins proved to be differentially expressed (Additional file 1: Table S2). The Gene Ontology analysis showed that they were related to cell growth/maintenance and apoptosis, energy pathways and signal transduction, membrane protein metabolism, protein biosynthesis, folding and degradation, cytoskeleton, and motility (Additional file 6: Table S3 and S4). Both in the invariant protein group and the differential protein group the “catalytic activity” was the most abundant class (33.3 % vs. 52.9 %, respectively). The “structural molecule activity” was also well represented (21.6 % vs.13.7 %, unvaried vs. differential proteins) (Fig. 5a,b) [34–46].
The protein-interaction networks (PIN) for the unvaried protein groups revealed strong interactions between structural proteins (Fascin-1, cytoskeletal keratin 19), chaperones, including two heat shock proteins (HSP70, HSP90), DNA binding protein (i.e., elongation factors EEF1A, EEF1B2) and remodeling and proteosome complex (PSMA5,4,6). This profile is consistent with a stem phenotype that requires plasticity and good proliferative ability [47, 48].
Interestingly, the hypoxia-inducible factor-1 (HIF-1) pathway is a central node of differential protein group PIN with close correlation between superoxide dismutase 2 (SOD2), thioredoxin (TXN), peroxiredoxin 6 (PRDX6), and catepsin D (CTSD). It has been shown that HIF-1 is involved in self-renewal and maintenance of pluripotency, by promoting the expression of the putative stemness genes (OCT4, SOX2, NANOG, Klf4) [49–51].
Multilineage capability of f-LSCs
We tested the differentiation capability of f-LSCs at an early passage (P3) and a late passage (P10). Both were differentiated in vitro using osteoblastic, adipose, and chondrogenic home-made media, as described in Methods. Twenty-one days after osteoblastic induction, Alizarin red S revealed the presence of extracellular matrix mineralization in both f-LSC cultures suggesting osteoblastic differentiation (Fig. 6a). lcian blue positivity confirmed the presence of acid polysaccharides (such as glycosaminoglycans) suggesting chondrogenic differentiation 21 days after induction (Fig. 6b).
Moreover, oil red positivity confirmed the presence of neutral triglyceride and lipid vacuoles, suggesting adipogenic differentiation after 28 days of culture (Fig. 6c).
Ease of access and immune-privileged status make the eye an ideal organ for its potential application in regenerative medicine. In our study, we confirmed the ease of isolation of the subpopulation of f-LSCs, which make them ideal candidates for adult stem cell therapy. In spite of the great interest in clinical applications, no exhaustive information is yet available on f-LSCs, especially on their exact characteristics. The International Society for Cellular Therapy has defined criteria to identify MSCs . We adopted these criteria to characterize f-LSCs from the limbus in our study. In addition, this work provides convincing evidence that limbal stroma contains f-LSCs. The latter subpopulation resembles stem cells which reside in the stem cell niche, where they are maintained in an undifferentiated state. These observations were supported by Notara and Daniels  who showed that LESCs are discretely located in the basal layer of the corneal limbal epithelium, at the junction between the transparent cornea and the opaque sclera. The limbal palisades of Vogt have been proposed as the site of the LESC niche whereas f-LSCs seem to reside in the extracellular matrix [3, 4]. Corneal stromal stem cells have been located in the anterior stroma subadjacent to the basal side of the palisades of Vogt [54–60]. In our study we did not aim to define the exact localization of f-LSCs but confirmed that they are completely different from LESCs, as demonstrated by ΔNp63 and CK12 and CK19 negativity, and they can be easily isolated and grown.
In our research we showed that f-LSCs fulfilled the criteria for multipotency established today.
In particular, our group investigated the possible effects of long-term culture conditions and donor age on the f-LSC phenotype. The cells isolated did not display surface expression of any hematopoietic marker (CD34 and CD45) and HLA-DR. By contrast, they expressed a variety of stemness markers (SSEA4, ABCG2, OCT3/4, and NANOG).
They showed low adhesion growth capability as limbo-spheres and they should proliferate in vitro as adherent cells for long-term culture, maintaining positivity for ABCG2, OCT3/4, NANOG, SOX2, THY-1, and SSEA4 and negativity for CD34, CD45, and HLD-DR; and they have multilineage differentiation potential under in-vitro culture conditions. By comparing the cell cycle distribution and the proliferation curve in the early and late passages there was a weak increase in the percentage of cells in phase G1 with a concomitant decrease in the percentage of cells in phase S. This event is reflected as a slight elongation in the doubling time in late passage culture but was not analyzed in more detail.
The evaluation of molecular expression of f-LSCs isolated from different donors (with the same or different age) and at different culture passages revealed a strong stability of the “limbal stem molecular pattern”, not affected by long-term culture and donor age. No significant differences were found regarding specific osteogenic, adipogenic, or chondrogenic staining.
For the first time, we constructed a two-dimensional electrophoresis (2-DE) proteomic pattern of cultured f-LSCs. The derived pattern of resolved protein spots was highly consistent; we identified 164 proteins stably expressed in the different conditions, defining profiles enriched in proteins linked to cell plasticity and proliferative and self-renewal capability. Furthermore, overlapping protein expression profiles confirm the stability of the stem cell phenotype with higher expression of structural proteins and proteins involved in the stem molecular pathway.
f-LSCs represent a robust source of MSCs independently of donor age and in-vitro conditions; the stability of their proteomic pattern could be very promising, suggesting their potential use in regenerative medicine.
ABCG2, ATP-binding cassette sub-family G member 2; BM-MSC, bone marrow-derived mesenchymal stromal cell; f-LSC, fibroblast-like limbal stem cell; LESC, limbal epithelial stem cell; MSC, mesenchymal stromal cell; NANOG, Homeobox protein NANOG; OCT4, octamer-binding transcription factor 4; P, passage of culture; PBS, phosphate-buffered saline; PIN, protein-interaction networks; SOX2, SRY (sex determining region Y)-box 2; SSEA4, stage-specific embryonic antigen-4; ΔNp63, a splice variant of p63
Davanger M, Evensen A. Role of the pericorneal papillary structure in renewal of corneal epithelium. Nature. 1971;229:560–1.
Schermer A, Galvin S, Sun TT. Differentiation-related expression of a major 64 K corneal keratin in vivo and in culture. J Cell Biol. 1986;103:49–62.
Tseng SC. Concept and application of limbal stem cells. Eye. 1989;3:141–57.
Goldberg MF, Bron AJ. Limbal palisades of Vogt. Trans Am Ophthalmol Soc. 1982;80:155–71.
Schofield R. The stem cell system. Biomed Pharmacother. 1983;37:375–80.
Branch MJ, Hashmani K, Dhillon P, et al. Mesenchymal stem cells in the human corneal limbal stroma. Invest Ophthalmol Vis Sci. 2012;53:5109–16.
Li GG, Zhu YT, Xie HT, et al. Mesenchymal stem cells derived from human limbal niche cells. Invest Ophthalmol Vis Sci. 2012;53:5686–97.
Xie HT, Chen SY, Li GG, et al. Isolation and expansion of human limbal stromal niche cells. Invest Ophthalmol Vis Sci. 2012;53:279–86.
Dravida S, Pal R, Khanna A. The transdifferentiation potential of limbal fibroblast-like cells. Brain Res Dev Brain Res. 2005;160(2):239–51.
Chen Z, De Paiva CS, Luo L, et al. Characterization of putative stem cell phenotype in human limbal epithelia. Stem Cells. 2004;22:355–66.
Du Y, Funderburgh ML, Mann MM, et al. Multipotent stem cells in human corneal stroma. Stem Cells. 2005;23(9):1266–75.
Katikireddy KR, Dana R, Jurkunas UV. Differentiation potential of limbal fibroblasts and bone marrow mesenchymal stem cell to corneal epithelial cells. Stem Cells. 2014;32:717–29.
Polisetty N, Fatima A, Madhira SL, et al. Mesenchymal cells from limbal stroma of human eye. Mol Vis. 2008;14:431–42.
Bray LJ, Heazlewood CF, Atkinson K, et al. Evaluation of methods for cultivating limbal mesenchymal stromal cells. Cytotherapy. 2012;14:936–47.
Ainscough SL, Linn ML, Barnard Z, et al. Effects of fibroblast origin and phenotype on the proliferative potential of limbal epithelial progenitor cells. Exp Eye Res. 2011;92:10–9.
West JD, Dorà NJ, Collinson JM. Evaluating alternative stem cell hypotheses for adult corneal epithelial maintenance. World J Stem Cells. 2015;7(2):281–99.
Huang M, Wang B, Wan P, et al. Roles of limbal microvascular net and limbal stroma in regulating maintenance of limbal epithelial stem cells. Cell Tissue Res. 2015;359(2):547–63.
Soma T, Hayashi R, Sugiyama H, et al. Maintenance and distribution of epithelial stem/progenitor cells after corneal reconstruction using oral mucosal epithelial cell sheets. PLoS One. 2014;9(10):e110987.
Criscimanna A, Zito G, Taddeo A, et al. In vitro generation of pancreatic endocrine cells from human adult fibroblast-like limbal stem cells. Cell Transplant. 2012;21:73–90.
Takahashi K, Tanabe K, Ohnuki M, et al. Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell. 2007;131(5):861–72.
Yu J, Vodyanik MA, Smuga-Otto K, et al. Induced pluripotent stem cell lines derived from human somatic cells. Science. 2007;318(5858):1917–20.
De Paiva CS, Chen Z, Corrales RM, et al. ABCG2 transporter identifies a population of clonogenic human limbal epithelial cells. Stem Cells. 2005;23:63–73.
Kim M, Turnquist H, Jackson J, et al. The multidrug resistance transporter ABCG2 (breast cancer resistance protein 1) effluxes Hoechst 33342 and is overexpressed in hematopoietic stem cells. Clin Cancer Res. 2002;8:22–8.
Di Ioro E, Barbaro V, Ruzza A, et al. Isoforms of DeltaNp63 and the migration of ocular limbal epithelial cells. Stem Cells. 2005;23(1):63–73.
Pellegrini G, Dellambra E, Golisano O, et al. p63 identifies keratinocyte stem cells. Proc Natl Acad Sci U S A. 2001;98(6):3156–61.
Barbaro V, Testa A, Di Ioro E, et al. C/EBP delta regulates cell cycle and self-renewal of human limbal stem cells. J Cell Biol. 2007;177(6):1047–9.
Notara M, Shortt AJ, O’Callaghan AR, et al. The impact of age on the physical and cellular properties of the human limbal stem cell niche. AGE. 2013;35:289–300.
Chang CYA, McGhee JJ, Green CR, et al. Comparison of stem cell properties in cell populations isolated from human central and limbal corneal epithelium. Cornea. 2011;30(10):1155–62.
Livak KJ, Schimittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2–ΔΔCt method. Methods. 2001;25:402–8.
Pucci-Minafra I, Cancemi P, Albanese NN, et al. New protein cluster of breast cancer tissue proteomics using actin content as a cellularity indicator. J Proteomic Res. 2008;7:1412–8.
Pucci-Minafra I, Cancemi P, Fontana S, et al. Expanding the protein catalogue in the proteome reference map of human breast cancer cells. Proteomics. 2006;6(8):2609–26.
Ashburner M, Ball CA, Blake JA, et al. The Gene Ontology Consortium. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.
Szklarczyk D, Franceschini A, Wyder S, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43(Database issue):D447–52.
Perrini S, Laviola L, Cignarelli A, et al. Fat depot-related differences in gene expression, adiponectin secretion, and insulin action signalling in human adipocytes differentiated in vitro from precursor stromal cells. Diabetologia. 2008;51:155–64.
Sakaguchi M, Shingo T, Shimazaki T, et al. A carbohydrate-binding protein, Galectin-1, promotes proliferation of adult neural stem cells. Proc Natl Acad Sci U S A. 2006;103:7112–7.
Vasta GR. Galectins as pattern recognition receptors: structure, function, and evolution. Adv Exp Med Biol. 2012;946:21–36.
Kandri T, Lataillade JJ, Doucet C, et al. Proteomic study of galectin-1 expression in human mesenchymal stem cells. Stem Cell Dev. 2005;14:204–12.
Spike BT, Kelber JA, Booker E, et al. CRIPTO/GRP78 signaling maintains fetal and adult mammary stem cells ex vivo. Stem Cell Rep. 2014;2:427–39.
Luo S, Mao C, Lee B, et al. GRP78/BiP is required for cell proliferation and protecting the inner cell mass from apoptosis during early mouse embryonic development. Mol Cell Biol. 2006;26:5688–97.
Taipale M, Tucker G, Peng J, et al. A quantitative chaperone interaction network reveals the architecture of cellular protein homeostasis pathways. Cell. 2014;158:434–48.
Miharada K, Karlsson G, Rehn M, et al. Cripto regulates hematopoietic stem cells as a hypoxic-niche-related factor through cell surface receptor GRP78. Cell Stem Cell. 2011;9:330–44.
Witke W, Sutherland JD, Sharpe A, et al. Profilin I is essential for cell survival and cell division in early mouse development. Proc Natl Acad Sci U S A. 2001;98:3832–6.
Zheng J, Lu Z, Kocabas F, et al. Profilin 1 is essential for retention and metabolism of mouse hematopoietic stem cells in bone marrow. Blood. 2014;123:992–1001.
Chaerkady R, Kerr CL, Kandasamy K, et al. Comparative proteomics of human embryonic stem cells and embryonal carcinoma cells. Proteomics. 2010;10:1359–73.
Blanchoin L, Pollard TD, Mullins RD. Interactions of ADF/cofilin, ArP3/3 complex, capping protein and profilin in remodeling of branched actin filament networks. Curr Biol. 2000;10:1273–82.
Collet JF, Messens J. Structure, function and mechanism of thioredoxin proteins. Antioxid Redox Signal. 2010;13:1205–16.
Pipino C, Pierdomenico L, Di Tomo P, et al. Molecular and phenotypic characterization of human amniotic fluid-derived cells: a morphological and proteomic approach. Stem Cells Dev. 2015;24(12):1415-28. [Epub ahead of print]
Roche S, Delorme B, Oostendorp Robert AJ, et al. Comparative protoemic analysis of human mesenchym and embryonic stem cell: towars the definition of a mesenchymal stem cell proteomic signature. Proteomics. 2009;9:223–32.
Mimeault M, Batra SK. Hypoxia-inducing factors as master regulators of stemeness properties and altered metabolism of cancer- and metastasis-initiating cell. J Cell Mol Med. 2013;17:30–54.
De Miguel MP, Alcaina Y, Sainz de la Maza D, et al. Cell metabolism under microenvironmental low oxygen tension levels in stemness, proliferation and pluripotency. Curr Mol Med. 2015;15(4):343–59.
Saito S, Lin YC, Tsai MH. Emerging roles of hypoxia-inducible factors and reactive oxygen species in cancer and pluripotent stem cells. Kaohsiung J Med Sci. 2015;31:279–86.
Dominici M, Le Blanc K, Mueller I, et al. Minimal criteria for defining multipotent mesenchymal stromal cells. The International Society for Cellular Therapy position statement. Cytotherapy. 2006;8:315–7.
Notara M, Daniels JT. Biological principals and clinical potentials of limbal epithelial stem cells. Cell Tissue Res. 2008;331(1):135–43.
Pellegrini G, Traverso CE, Franzi AT, et al. Long-term restoration of damaged corneal surfaces with autologous cultivated corneal epithelium. Lancet. 1997;349(9057):990–3.
Sangwan VS, Matalia HP, Vemuganti GK, et al. Clinical outcome of autologous cultivated limbal epithelium transplantation. Indian J Ophthalmol. 2006;54(1):29–34.
Kolli S, Ahmad S, Lako M, et al. Successful clinical implementation of corneal epithelial stem cell therapy for treatment of unilateral limbal stem cell deficiency. Stem Cells. 2010;28(3):597–610.
Rama P, Matuska S, Paganoni G, et al. Limbal stem-cell therapy and long-term corneal regeneration. N Engl J Med. 2010;363(2):147–55.
Limb GA, Daniels JT. Ocular regeneration by stem cells: present status and future prospects. Br Med Bull Rev. 2008;85:47–61.
Amano S, Yamagami S, Mimura T, et al. Corneal stromal and endothelial cell precursors. Cornea. 2006;25:S73–7.
Pinnamaneni N, Funderburgh JL. Concise review: Stem cells in the corneal stroma. Stem Cells. 2012;30:1059–63.
The authors thank A Criscimanna and G Zito for their efforts in the first part of the limbus research area.
This work is dedicated to our former director Prof. Aldo Galluzzo, who prematurely died in 2011 when the Laboratory of Regenerative Medicine was created and research in this field started.
LT was responsible for conception and design, collection and assembly of data, data analysis and interpretation, and manuscript writing. RM was responsible for collection and assembly of data and revision of the manuscript. GC and SC were responsible for conception and design, provision of study material or patients, and revision of the manuscript. GP, MP, AC, and WA were responsible for data analysis and interpretation and drafting the manuscript. GDC was responsible for acquisition of data and revision of the manuscript. IP-M was responsible for manuscript writing and revision of the manuscript. CG was responsible for conception and design, data analysis and interpretation, manuscript writing and final approval of manuscript, manuscript drafting, revising critically for important intellectual content, and financial support. All authors agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. All authors read and approved the final manuscript.
Financial competing interests
The work was partially funded by PON 01_00829 2007/2013 with European Community funds administered by the Italian Ministry for the University and RIMEDRI PO FESR 2007/2013 (CG was Scientific Director for both).
The authors declare that they have applied for a patent concerning the possible use of f-LSCs in type 1 diabetes.
Nonfinancial competing interests
The authors declare that they have no nonfinancial competing interests in relation to this manuscript.
Table S2 Catalogue of the unvaried and differential expressed proteins in f-LSCs samples. (ZIP 110 kb)
Figure S1 Positive Control of cytofluorimetric assay: Human bone marrow mesechymal stem cells (BM-MSCs) stained positive for SSEA4 (A), ABCG2 (B), OCT4 ( C ) and NANOG (D); HeLa cells stained positive for ΔNp63 (E). A cell suspension of BM-MSCs and HeLa was obtained by trypsinization and centrigugation. The cell pellet was resuspended in 100μl of PBS to test SSEA4 or ABCG2. The samples were resuspended in PBS supplemented with 0.1% saponin and 1% BSA for 20 minutes for permeabilization cells and finally stained for OCT4 or NANOG or ΔNp63. (TIF 686 kb)
Table S1 Catalogue of the protein spots identified in the 2D-IPG prototype map of f-LSCs. Table reports the following information: Protein names, accession numbers (AC) and abbreviated names correspond to the nomenclature used in the Swiss-Prot database. The experimental values of pI and MW for every isoelectric spot were calculated with ImageMaster 2D Platinum system; the theoretical values represent the predicted MW and pI for each identified protein according to Swiss-Prot and TrEMBL database. Identification methods: 1, MALDI-TOF; 2, N-terminal sequencing by automated Edman degradation and 3, Western Blotting. (DOC 198 kb)
A) The proteomic map of P4 f-LSCs from >45 years donor; B) The proteomic map of P20 f-LSCs from >45 years donor f-LSCs; C, D) The proteomic maps of P10 f-LSCs from 2 donors < 45 years; E) The proteomic map of P20 f-LSCs from one 64 years old donor. (TIF 2807 kb)
Figure S3 Figure.5: A) Proteomic profile analysis of P4 f-LSCs>45 years vs. f-LSC P20 f-LSCs >45 years, overlapped for about 90% of proteins expressed; B) f-LSCs from 2 different donors of the same age (<45 years) overlapped for about 83% of proteins; C) proteomic profiles of the same cutlure passage f-LSCs from 2 donors of different ages (< 45 years vs. >45 years) overlapped for about 68%. P= passage. (TIF 25509 kb)
Table S3. The genes expressing the invariant proteins were analyzed for their enrichment in specific gene ontology annotation. Table S4. The genes expressing the differenzial proteins were analyzed for their enrichment in specific gene ontology annotation. (XLS 16 kb)