- Research article
- Open Access
Transcriptome analysis for the identification of cellular markers related to trabecular meshwork differentiation
BMC Genomics volume 18, Article number: 383 (2017)
Development of primary open-angle glaucoma (POAG) is associated with the malfunctioning trabecular meshwork (TM). Cell therapy offers great potential for the treatment of POAG, but requires the generation of functional TM cells in vitro to replace the lost/dysfunctional cells. TM differentiation in vitro from various stem cell types must be monitored by the expression of specific markers. However, no single definitive marker of the TM has been identified.
To identify robust markers of TM differentiation, we performed global transcriptome profiling using high-density oligonucleotide microarray on ex vivo TM tissue and cultured TM progenitors. Corneal and scleral tissues were also used in the analysis. After removal of genes expressed in the cornea and sclera, 18 genes were identified that were differentially expressed in the TM relative to the other samples. CDH23, F5, KCNAB1, FGF9, SPP1, and HEY1 were selected among the genes highly expressed in the TM, together with BDNF which was repressed, compared to progenitors for further investigation. Expression analysis by qPCR verified the differential expression and immunofluorescence of the anterior segment confirmed strong expression in the TM.
Three independent cohort of expression studies have identified novel markers, fitting in identifying TM cells and in evaluating directed TM differentiation in vitro.
Primary open-angle glaucoma (POAG) is a leading cause of irreversible blindness worldwide. It is defined by a characteristic pattern of damage to the optic nerve head resulting in a progressive loss of vision. Elevated intraocular pressure (IOP) is an important risk factor for axonal damage at the optic nerve head [1,2,3]. In POAG patients, abnormalities in the trabecular meshwork (TM) lead to increased resistance to aqueous humour drainage and impedes outflow from the anterior chamber [4,5,6].
The TM located at the corneoscleral limbus is responsible for draining approximately 90% of the aqueous humour from the anterior chamber via conventional outflow pathway . Alterations in the composition and deposition of extracellular matrix (ECM) together with significant decline in the number of TM cells have been observed in the diseased tissue of POAG specimens [8,9,10,11,12,13]. Accumulation of the ECM components, fibronectin, elastin and fibrillar material change the ECM composition [8, 9]. Glycosaminoglycan profiling also found an increase in chondroitin sulfate and decline in hyaluronan in the diseased tissues [10, 11]. Other structural modifications include thickening and fusion of the trabeculae, and deposition of sheath-derived (SD) plaques [8, 12]. Furthermore, remnant cells in the meshwork seem to be under stress, with increased expression of stress proteins, crystallin alpha B and inducible nitric oxide synthase . Genome-wide expression studies found altered expression of certain genes with the pathology as well, such as those pertaining to inflammation, acute-phase response and G-protein signalling being elevated, while antioxidants and members of the solute carrier family were downregulated among others . Such aberrations in the TM are believed to interfere with its filtering mechanism and its ability to control IOP.
Recent work has established the therapeutic potential of TM cell replacement. TM stem cells, mesenchymal stem cells (MSC) and induced pluripotent stem cells (iPSC)-derived TM cells promoted TM regeneration and in some cases also alleviated the IOP in models of OAG [15,16,17,18]. To realize the potential of such therapies it is critical to establish means of generating functional TM cells in vitro; newly discovered stem cells of the TM are enabling this potential. We and others recently established methodologies to isolate and propagate stem cells from the TM [19, 20]. We established that these progenitors display the defining characteristics of mesenchymal stem cells (MSC) , and termed them TM-derived MSC (TM-MSC). In our quest for cell replacement therapy for POAG, we are interested in establishing an expandable source of functional TM cells from their progenitors. To this end, we are establishing methods to direct the differentiation of TM-MSC into functional TM tissue. An impediment to this effort is the lack of cellular markers that are required to monitor differentiation of the cells and to characterize their purity.
Since no single definitive marker is available to identify TM cells, a combination of genes, highly expressed in the TM, is utilized to assign TM lineage (Table 1). Expression analysis found high expression of aquaporin 1 (AQP1) in primary human TM cells . Gene expression profiling of TM cells followed by comparative analysis with Schlemm’s canal endothelium and subsequent reporter construct activity demonstrated the specificity of chitinase 3-like 1 (CHI3L1) to the TM in perfused anterior segments . Similarly, matrix gla protein (MGP) was identified to be highly expressed in the TM [24,25,26]. Myocilin (MYOC) is also reported to be highly expressed in the TM with its expression alterable by glucocorticoids and topography among other factors . In fact, this gene was first discovered and characterized in the TM as trabecular meshwork-inducible glucocorticoid response (TIGR) gene [28,29,30]. Matrix metallopeptidase 1 (MMP1) and ankyrin 3 (ANK3) are also among genes that have high expression in primary TM cells compared to cultured Schlemm’s canal endothelial cells and scleral fibroblasts as deduced by gene expression profiling . Another marker, ELAM1 is expressed in glaucomatous TM cells . These markers are commonly used collectively to distinguish the TM. In most instances, their specificity to the TM in the context of the anterior segment has not been demonstrated and it is likely that they are also highly expressed in the surrounding cell types. Moreover, complete molecular dissection of genes differentially expressed between the TM progenitors and mature TM cells has not been explored.
To identify suitable markers, genome-wide expression profiling of ex vivo TM tissue and cultured TM progenitors was performed to reveal genes differentially expressed between the mature and progenitor cells of the TM. To qualify as good differentiation markers, the genes have to be (1) differentially expressed between the TM progenitors and mature TM cells and (2) ideally specific to the TM in the anterior segment of the eye. From this study, we identified a panel of differentiation markers for the identification of mature TM cells as well as indicators of TM differentiation.
Research-grade corneoscleral tissues from deidentified human cadaver donors were procured from the Lions Eye Institute for Transplantation and Research (Tampa, FL, USA) and handled in accordance to the tenets of the Declaration of Helsinki. Informed written consent was attained from donors who were then registered on the local organ and tissue donor registry. For unregistered donors, written consent was acquired from the next of kin of the deceased donor. All donors were within 50–70 years of age and had no known record of chronic illnesses. Serological tests ensured that the samples were free of transmissible diseases. The corneoscleral specimens were harvested within 24 h of death, preserved in Optisol-GS at 4 °C and processed for the study within 13 days from preservation.
Extraction of tissues
The corneoscleral specimen was washed several times in phosphate-buffered saline (PBS, Gibco; Thermo Fisher Scientific, Carlsbad, CA, USA) and stained with 0.1% trypan blue (Sigma, St. Louis, MO, USA) for 1 min. The dye stains the TM and scleral spur intensely with limited coloration of the corneal endothelium. After removal of excess stain with PBS rinses, the tissues were harvested by a technique adapted from Tripathi and Tripathi  in a moist chamber under a stereomicroscope taking sterile precautions. The corneal and scleral samples were extracted by making parallel cuts anterior to the Schwalbe’s line and posterior to the scleral spur respectively. The tissues were later minced into smaller pieces and homogenized in Trizol reagent (Ambion; Thermo Fisher Scientific) for RNA isolation. For TM extraction, the surrounding corneal endothelium and scleral spur were peeled off, and the TM was gently dissected from the underlying tissue with a crescent knife. The harvested TM was further processed for the isolation of TM progenitors or homogenized in Trizol reagent for gene expression studies by passing the sample several times through a blunt 20-gauge needle fitted to a nuclease-free syringe.
Isolation and propagation of TM progenitors
TM-MSC culture was established from the harvested TM tissue as previously described . Briefly, the TM tissue was digested with 2 mg/ml type I colleganese (Worthington Biochemical Corporation, Lakewood, NJ, USA) overnight at 37 °C in a humidified atmosphere of 5% CO2 in the cell culture incubator. The tissue was further dissociated in 0.05% trypsin-EDTA and triturated into a single cell suspension. The cells were then plated in low glucose DMEM containing L-glutamine and 100 mg/l sodium pyruvate, 10% fetal bovine serum (FBS), 10 mM non-essential amino acids, 100 units/ml Pencillin and 100 μg/ml Streptomycin. FBS and all other reagents were purchased from Gibco. Medium was changed every 2–3 days. The cells were passaged with 0.05% trypsin-EDTA when a confluence of 80–90% was reached and seeded at a ratio of 1:4. Characterization of the TM-derived cells as MSC was performed by flow cytometric analyses for MSC markers with accordance to the ISCT criteria . Early passage cells at p3 – p4 were utilized for the study.
Total RNA was extracted from the Trizol/cell lysate followed by purification and DNase I treatment with the RNeasy miniprep kit (Qiagen, Venlo, Netherlands). Briefly, phase separation was performed by adding 200 μl chloroform per 1 ml of Trizol and after 3 min incubation, the mixture was centrifuged at 12 000 g for 15 min at 4 °C. The aqueous phase was then mixed with one volume of 70% ethanol and transferred to a RNeasy spin column with collection tube for the RNA to be purified and eluted according to the manufacturer’s protocol with an additional step of in-column DNase I treatment. RNA yield was determined on the NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific). RNA quality was confirmed with the Agilent RNA 6000 pico kit on the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA).
Tissue samples obtained from three donors and TM-MSC isolated from three other donors were utilized for the study. RNA harvested from the samples was amplified with the TotalPrep RNA amplification kit according to manufacturer’s protocol (Ambion). The resultant purified biotin-labelled complementary RNA (cRNA) was quantified and 750 ng cRNA was applied to the HumanRef-8 v3.0 expression beadchip (Illumina, San Diego, CA, USA) following the direct hybridization instructions. Cy3 conjugated to streptavidin was used to detect the hybridized cRNA. The chip was scanned on the BeadArray Reader and imaged with the BeadScan software (both from Illumina).
The raw data were background-adjusted (as calculated from the negative control probes) and converted into expression profile by GenomeStudio (Illumina), and further analyzed with GeneSpring GX v14.5 (Agilent Technologies). The data were normalized by percentile shift to 75th percentile and threshold raw signal was set as 10. Unsupervised hierachical clustering was used to visualize gene expression similarity. The individual microarrays for the respective sample set were averaged and pairwise analysis was performed in comparison to TM, so relative intensity is the fold change relative to TM. Differentially expressed genes were analysed by one-way ANOVA, with correction for Benjamini & Hochberg false discovery rate (FDR) test. Data were filtered using fold change cut-off of 2.0 and P adj ≤ 0.05 (q-value; defined as statistical significance). Entity lists generated from the analyses are available under Additional file 1. Our data have been deposited in NCBI's Gene Expression Omnibus (GEO) repository and accessible through GEO series accession number GSE87526.
Quantitative real-time polymerase chain reaction
cDNA was generated from RNA extracted from tissues derived from three more donors and three other independent TM-MSC lines using the high-capacity cDNA reverse transcription kit (Applied Biosystems, Foster City, CA, USA) according to prescribed protocol. Quantitative real-time polymerase chain reaction (qPCR) was performed by direct dye binding utilizing SYBR Green (Applied Biosystems). Primers were designed by Primer3 [34, 35] or PrimerBank  with an annealing temperature of 55 °C where possible (Additional file 2: Table S1). cDNA was analysed in duplicates on the 7900HT Fast Real-Time PCR system and analyzed using SDS 2.4 software (both from Applied Biosystems). Amplification of ACTB was performed for each cDNA for normalization of RNA content. Relative mRNA abundance was calculated by the ∆∆CT method of comparative quantification. All values were presented as mean ± standard error of the mean (mean ± SEM) from three biological replicates. Significance was assessed using Student’s t-test and significant differences were considered as those with P < 0.05.
The corneoscleral specimen was washed several times in PBS to remove detached cells. The excess conjunctiva was removed with an angular iris dissecting scissors and the specimen was cut in half. The cornea and sclera were trimmed and the sample comprising the corneoscleral limbus was embedded in the optimal cutting temperature (OCT) compound (Tissue-Tek; Sakura Finetek USA, Torrance, CA, USA) in a cryomold and frozen at −80 °C. The frozen tissue sample was cut into 8 μm sections at the cryostat (Leica, Wetzlar, Germany), mounted onto microscope slides (Thermo Fisher Scientific), air-dried for 15 min and stored at −80 °C. The cryosections were thawed at room temperature and rehydrated with water before proceeding to fixation and permeabilization.
The cryosections were washed in PBS and fixed in 4% paraformaldehyde (PFA) (Sigma) pH 7.4 for 10 min at room temperature. The tissue sections were then washed twice in PBS for 5 min and permeabilized in 0.25% Triton X-100 (Sigma) for 10 min, followed by three washes. To block non-specific binding of antibodies, the samples were incubated with 10% heat-inactivated goat serum and 0.3 M Glycine (Sigma) in PBST for 30 min. The sections were then incubated in primary antibody, diluted at optimized concentration in PBST containing 1% serum, overnight at 4 °C and washed thrice. Alexa Fluor 594-conjugated secondary antibody in 1% serum was added and incubated for 1 h in the dark at room temperature and washed three times. Nuclei were stained with 1 μg/ml 4′,6-diamidino-2-phenylindole (DAPI) (Invitrogen, Carlsbad, CA, USA) for 1 min followed by several washes. Finally, coverslip was mounted onto the glass slide with fluorescence mounting medium (Dako, Denmark) and affixed with nail polish. Images were acquired under the camera mode of the fluorescence microscope (Zeiss, Germany). Images were archived under the Zeiss AxioVision 4.8 image analysis software and processed in Adobe Photoshop (Adobe Systems, San Hose, CA, USA). Independent cohort of cornealscleral specimens were used for the study. Details of the antibodies used for the study are listed in Additional file 2: Table S2.
Expression patterns of known TM markers
Expression studies by other groups have identified several genes highly expressed in the TM, which have been used as markers of the TM. To assess the specificity of their expression in the TM relative to the peripheral tissues of cornea and sclera, we performed qPCR analysis with tissues derived from three donors. Surprisingly, expression of AQP1 was much higher in the adjacent tissues relative to TM (Fig. 1). The remaining markers, MGP, CHI3L1 and ELAM1 had greater expression in the TM relative to cornea, but generally showed higher expression in the sclera. Only MYOC was highly expressed in the TM compared to the peripheral tissues but not significant statistically compared to sclera (P ≥ 0.05). Given that previously identified TM markers are lacking in specificity, additional markers are needed.
Gene expression profiling of TM-MSC and TM
In this study we sought to identify robust differentiation markers of the TM. To this end, we searched for genes that display high differential expression between the TM-MSC and mature TM cells. Ex vivo TM, mostly comprised of differentiated TM cells, was utilized as a source of mature TM cells. To identify genes with differential expression, gene expression profiling by microarray technology was performed on TM tissue and TM-MSC derived from three independent donors of the same age group (50–70 years). Signal intensities were normalized by the method of percentile shift (75th) and differential genes that passed the fold change cut-off of 2.0 with statistical significance, as determined by one-way ANOVA followed by Benjamini-Hochberg multiple testing corrected P ≤ 0.05 were identified.
The resultant entity list had a considerable number of genes differentially expressed between the two sample sets, with 1651 and 1967 probes corresponding to those genes elevated in the TM and TM-MSC, respectively (Fig. 2a). MYOC and angiopoietin-like 7 (ANGPTL7) were most elevated in the TM, being over 1000-fold higher in the TM compared to TM-MSC. The top 20 genes elevated in the TM are listed in Table 2. These genes code for proteins expressed in the different cell compartments, such as intracellular proteins (PCP4, KRT13, KRT5, MYH11, HBB, FAM107A), cell surface proteins (CDH23, HLA-DRα, VCAM1, CXADR, TACSTD2), and secreted proteins (MYOC, ANGPTL7, CCL3L3, SERPINA3, APOD, FCGBP, PTGDS, C2ORF40). H19, a long non-coding RNA was also elevated.
Interestingly, HLA-DRα which we showed previously is not expressed in TM-MSC was highly expressed in the mature cells . THY-1 (CD90), a MSC marker present in TM-MSC, was repressed in the TM. TM-MSC expressed alpha-kinase 2 (ALPK2) and lysyl oxidase (LOX) at the greatest level compared to the TM, at 137- and 71-fold, respectively. Other genes highly repressed in the TM are listed in Table 3.
We analysed our data to assess whether previously reported TM markers were differentially expressed between the mature TM cells and TM-MSC. MYOC and MGP were robustly elevated in the TM at 1374- and 67-fold, respectively. Expression of AQP1 was about 4-fold higher in the TM. The remaining markers, CHI3L1, MMP1 and ELAM1, were not found in the list of differentially expressed genes. Further investigation was done to determine whether the differentially expressed markers are specific to the TM region of the eye.
Gene expression profiling of cornea and sclera
To find genes uniquely expressed in the TM in the context of the anterior segment of the eye, corneal and scleral tissues, harvested from three donor specimens, were included in the microarray study and pairwise analysis relative to the TM was performed as described earlier. The resulting gene expression data were processed into two entity lists: TM vs cornea and TM vs sclera (fold difference ≥ 2.0 and P adj ≤ 0.05). A limited number of genes were found to have differential expression: 47 genes in the TM vs cornea and 164 in the TM vs sclera data sets. Both gene lists were merged to determine genes that are differentially expressed in the cornea and sclera in common with respect to the TM (i.e. TM vs cornea/sclera) (Fig. 2a). Twenty genes had significantly higher expression in the TM and three were lower in comparison to the cornea and sclera (Table 4).
Surprisingly, none of the previously identified TM markers were present in the merged TM vs cornea/sclera gene list, indicating that they were not specific to the TM. Hence, our approach was modified to discover new markers from this gene set that are also differentially expressed between the mature TM cells and TM-MSC. This could identify genes that are both ‘TM-specific’ in the context of the anterior segment and robust enough to track TM differentiation.
Identification of trabecular meshwork differentiation markers
Unsupervised hierarchical clustering showed that the tissues clustered together, apart from TM-MSC, with the TM having greater extent of expression similarity to sclera than cornea (Fig. 2b). This confirms the divergent state of the differentiated cells in the tissues from the progenitor state of TM-MSC. Overlap of the three gene lists (TM vs TM-MSC, TM vs cornea and TM vs sclera) was performed to identify TM unique genes. From the comparative analyses, 18 genes were found to be differentially expressed in the TM tissue compared to TM-MSC, cornea and sclera (Fig. 2c). Specifically, 16 genes (namely CDH23, SLC16A10, SPP1, F5, KCNAB1, MGC33846, FGF9, TMEM178, HEY1, NEB, SORCS1, CD1D, ZBTB46, SYTL4, FAIM2 and BDNF) were elevated and two were repressed (ITPA and SLC7A2) in the TM tissue relative to the other tissues. Similar trend was observed in the TM vs TM-MSC gene set, except BDNF which was repressed in the TM. Details of the genes, fold change and regulation are presented in Table 5.
Six genes elevated in the TM were selected for further analysis based on largest divergence in expression from the cornea and sclera, and availability of antibodies. Cadherin-related 23 (CDH23), secreted phosphoprotein 1 (SPP1), coagulation factor 5 (F5), potassium voltage-gated channel, shaker-related subfamily, beta member 1 (KCNAB1), fibroblast growth factor 9 (FGF9), and hairy/enhancer-of-split related with YRPW motif 1 (HEY1) were chosen as TM differentiation markers (Table 5). CDH23 had the highest differential expression (420-fold) in the TM relative to TM-MSC. Its expression in the TM in comparison to the cornea and sclera were 9- and 5-fold, respectively. The other selected genes were elevated in the TM by fold differences of 35–100 compared to TM-MSC, and 5–50 relative to the cornea/sclera. As stated earlier, brain-derived neurotropic factor (BDNF) was the only gene repressed in the TM tissue relative to TM-MSC in the merged gene list and thus selected as a TM-MSC marker (i.e. negative differentiation marker). While BDNF was repressed by 13-fold in the TM in comparison to TM-MSC, its expression was higher in the TM tissue relative to the adjacent tissues.
Verification by qPCR was performed to confirm the microarray results. Expression of the shortlisted differentiation markers was assessed in tissues obtained from three new donors and three other TM-MSC lines. Deriving the samples from new biological replicates increases the reliability of the study by eliminating donor or batch specific variations.
qPCR result confirmed the repression of the differentiation markers in the cornea and sclera compared to the TM (Fig. 3a). F5 had the lowest expression in the two peripheral tissues followed by KCNAB1 with fold difference close to 0.01 relative to the TM (P < 0.0001). Remaining markers were also significantly repressed with fold changes between 0.1 and 0.3 (P < 0.05). This confirms higher expression of the selected markers in the TM relative to the cornea and sclera.
In comparison to TM-MSC, the positive differentiation markers were all expressed significantly higher in the TM, while BDNF was repressed (P < 0.05) (Fig. 3b). The positive markers were more than 100-fold higher, with F5 and CDH23 having strongest expression in the TM relative to TM-MSC. BDNF was substantially lower at 0.3-fold in the TM. The qPCR results from independent samples confirmed the differential expression of the identified markers in the two cell types. Thus, from the mRNA expression analysis, the markers are: 1) specific to the TM relative to the cornea and sclera; and 2) robustly differential in expression between the TM and TM-MSC.
Immunofluorescence of selected markers in the TM
To ensure the markers are highly expressed in the TM compared to the rest of the anterior segment, immunofluorescence was performed on corneoscleral tissue. Corneoscleral tissue sections were prepared from donor eyes by trimming off the excess cornea and sclera, embedded in OCT compound and cryosectioned transversely (Fig. 4). Frozen sections were fixed and permeablized, followed by immunostaining with commercially available antibodies to assess the expression of markers in the TM and surrounding tissues. Two markers, SPP1 and FGF9, stained brightly in the TM and were undetectable in the surrounding tissues (Fig. 5). Other positive markers of TM differentiation, KCNAB1, CDH23, F5 and HEY1 were clearly detected in the TM, though a weak signal was also detected in peripheral tissues, particularly in the sclera, and in some cases in the Schlemm’s canal. BDNF which was elevated in the TM relative to the cornea and sclera at the transcript level, was of stronger intensity in the TM compared to the surrounding tissues, but not as pronounced as the positive markers. As a TM-MSC marker, BDNF was expected to be expressed most intensely in the ‘insert portion’ of the TM, a region under the Schwalbe’s line of the cornea where TM progenitors are believed to reside. However, this was not apparent probably due to the small number of stem cells in this region. Corresponding staining in the TM-MSC showed negligible signal for all the TM markers, while BDNF expression was detected (Additional file 3: Figure S1).
The pathophysiological mechanisms causing POAG are not completely understood. However, the association of TM dysfunction with POAG, as characterized by structural changes in the TM, has been well established. Since, the structural and functional integrity of the TM is indispensable in maintaining normal IOP and healthy vision, cellular therapy is being explored to replace the lost or dysfunctional cells with healthy new ones [15, 16]. The capacity of progenitors of the TM to integrate into the TM tissue and differentiate has already been demonstrated .
Expression of a combination of genes highly expressed in the TM is generally accepted to identify the TM cells, however most currently known markers are also expressed in other ocular cells. In our attempt to find markers better suited as TM differentiation markers to track differentiation from their progenitors, we generated a list of genes differentially expressed in the TM-MSC and ex vivo TM tissue by comprehensive transcriptome profiling. To our knowledge, this is the first report dechipering molecular level differences between the mature TM cells and a pure population of its progenitors.
There are scattered reports on certain genes differentially expressed in the filtering compartment of the TM and its non-filtering stem cell niche. Through immunocharacterization in vivo, it was reported that ANK3 and HMFG1 are highly expressed in the niche, while CHI3L1 is low compared to the TM tissue . Contrary to this report, we found that ANK3 is upregulated in the TM compared to TM-MSC. In addition, HMFG1 and CHI3LI were not differentially expressed. Several other genes highly expressed in the TM (MYOC, AQP1, MGP and MMP1) were higher in the TM relative to TM-MSC. However, they have similar expression in the cornea and sclera, rendering them non-specific to the TM.
MYOC and ANGPTL7 were the most differentially expressed genes in the TM relative to TM-MSC. It had been shown previously by serial analysis of gene expression (SAGE) profiling that MYOC and ANGPTL7 are among the 40 most highly expressed genes in the TM . Besides the TM, MYOC expression was detected in the cornea, sclera, choroid, ciliary body, iris, lamina cribosa, retina and optic nerve as well as in the acellular vitreous and aqueous humour [38,39,40,41]. ANGPTL7 is highly expressed in the cornea and lens . Its expression in other ocular tissues is not well understood. To elucidate markers that are robust and specific to the TM, genes enriched in the TM tissue in relation to cornea, sclera and TM-MSC were identified through comprehensive transcriptome profiling and selected candidates were further evaluated. Relative mRNA levels, as measured by qPCR, confirmed highest expression of the transcripts in the TM. Immunocharacterization on corneoscleral sections also demonstrated high marker expression in the TM. Some markers displayed specificity to the TM while the rest had greater expression in the TM in the context of the anterior segment of the eye. In addition, BDNF which was repressed in the TM in relation to TM-MSC was identified as a negative indicator of TM differention.
The newly identified markers of TM differentiation (TM markers) include secreted proteins F5, SPP1 and FGF9, membrane proteins CDH23 and KCNAB1, and transcription factor HEY1. F5, a key cofactor for prothrombinase activity, is part of the coagulation process. The gene has been mapped to the glaucoma locus, GLC1A where MYOC also resides and is known to be downregulated in cultured human TM cells . Similarly SPP1, a phosphoprotein involved in matrix mineralization, is also downregulated in human TM cells in vitro and highly expressed in the TM tissue [24, 25, 43]. Primary TM cells in these studies will likely contain TM progenitors and repression of both proteins compared to the TM tissue support our data. SPP1 is also a cytokine that promotes interferon gamma (IFNG) and interleukin 12 (IL12) production, hence functioning in immune response [44, 45]. Other activities associated with SPP1 are migration, ECM remodelling, cell adhesion and survival . Interestingly, it is an inhibitor of calcification like MGP, another well established marker of the TM . FGF9 is involved in biological processes important for development, such as tissue repair, proliferation and differentiation of several cell types including induction of neural crest from which the TM was derived . Receptors for FGF signalling have been detected in cultured human TM cells and their activation stimulates proliferation and extracellular acidification . The association of these genes with TM function and maintenance is unknown.
CDH23 is a member of the cadherin superfamily that encodes calcium-dependent cell-cell glycoprotein. It plays an essential role in auditory function as mutations in the gene leads to deafness in humans and mice [49, 50]. As an integral part of stereocilia and hair bundle of cochlear hair cells [51, 52], it has a putative role in mechanotransduction. It is tempting to speculate that high expression of CDH23 in the TM may offer a similar mechanosensitivity to the tissue, allowing IOP homeostasis. It may also be involved in the structural integrity of the TM. KCNAB1 heterotetramerizes with alpha subunits and modulates the activity of the Kv channels [53,54,55]. It is involved in a range of functions including but not limited to regulating neurotransmitter release, heart rate, insulin secretion, smooth muscle contraction and cell volume. Transcriptional repressor, HEY1 which is activated by Notch signalling [56, 57], is important in development . It is essential for cardiovascular development [59, 60], and implicated in neurogenesis , somitogenesis and osteogenesis . The roles KCNAB1 and HEY1 play in the TM mechanism are yet to be known. BDNF which was recognized as a TM-MSC marker in this study is a neurotrophin involved in the regulation of survival, differentiation and stress response. It is known to be expressed by both cultured human TM cells and ex vivo TM tissue, with the former expressing it at a higher level . BDNF polymorphism has also been implicated with POAG progression , although its exact pathophysiological mechanism is unclear.
The screening of markers was performed on samples attained from new sets of donors for the various expression studies, thereby minimizing batch-specific and donor-specific biases. Moreover, ex vivo TM used in the study instead of cultured primary TM cells offers a more accurate molecular dissection of the mature TM cells. Marker gene expression was investigated in a commercially available TM cell line by qPCR analysis; four of the markers (F5, HEY1, FGF9 and KCNAB1) were found to be significantly elevated, while BDNF was repressed relative to TM-MSC (Additional file 4: Figure S2). We are attempting to establish conditions that induce TM differentiation from TM-MSC in vitro. From our preliminary results, we observed partial induction of F5, CDH23, SPP1, HEY1, KCNAB1 and FGF9, indicating that these markers may be useful to monitor the process of differentiation towards TM lineage (unpublished data). Thus, this panel of markers together with previously established TM markers could allow rigorous assessment of differentiation in transplantation studies targeting the TM as well as in the generation of mature TM cells in vitro for molecular applications.
Our study has identified a set of markers of TM differentiation. It is better to use the markers collectively rather than a single identifier to ensure robust differentiation. We have confirmed by independent means of expression analysis that these markers are: (1) differentially expressed between the TM and TM-MSC and 2) specific to the TM relative to other regions in the anterior segment. Upon successful differentiation of TM-MSC into TM lineage, the positive markers (TM markers) are expected to be upregulated, whereas BDNF being a negative indicator (TM-MSC marker) ideally should be downregulated.
Collaborative Normal-Tension Glaucoma Study Group. Comparison of glaucomatous progression between untreated patients with normal-tension glaucoma and patients with therapeutically reduced intraocular pressures. Am J Ophthalmol. 1998;126:487–97.
Collaborative Normal-Tension Glaucoma Study Group. The effectiveness of intraocular pressure reduction in the treatment of normal-tension glaucoma. Am J Ophthalmol. 1998;126:498–505.
The AGIS Investigators. The advanced glaucoma intervention study (AGIS): 7. the relationship between control of intraocular pressure and visual field deterioration. Am J Ophthalmol. 2000;130:429–40.
Grant WM. Clinical measurements of aqueous outflow. AMA Arch Ophthalmol. 1951;46:113–31.
Seiler T, Wollensak J. The resistance of the trabecular meshwork to aqueous humor outflow. Graefes Arch Clin Exp Ophthalmol. 1985;223:88–91.
Mäepea O, Bill A. The pressures in the episcleral veins, Schlemm’s canal and the trabecular meshwork in monkeys: effects of changes in intraocular pressure. Exp Eye Res. 1989;49:645–63.
Grant WM. Experimental aqueous perfusion in enucleated human eyes. Arch Ophthalmol. 1963;69:783–801.
Liitjen-Drecoll E, Futa R, Rohen JW. Ultrahistochemical studies on tangential sections of the trabecular meshwork in normal and glaucomatous eyes. Invest Ophthalmol Vis Sci. 1981;21:563–73.
Babizhayev MA, Brodskaya MW. Fibronectin detection in drainage outflow system of human eyes in ageing and progression of open-angle glaucoma. Mech Ageing Dev. 1989;47:145–57.
Knepper PA, Goossens W, Palmberg PF. Glycosaminoglycan stratification of the juxtacanalicular tissue in normal and primary open-angle glaucoma. Invest Ophthalmol Vis Sci. 1996;37:2414–25.
Knepper PA, Goossens W, Hvizd M, Palmbergf PF. Glycosaminoglycans of the human trabecular meshwork in primary open-angle glaucoma. Invest Ophthalmol Vis Sci. 1996;37:1360–7.
Gottanka J, Johnson DH, Martus P, Lütjen-Drecoll E. Severity of optic nerve damage in eyes with POAG is correlated with changes in the trabecular meshwork. J Glaucoma. 1997;6:123–32.
Alvarado J, Murphy C, Juster R. Trabecular meshwork cellularity in primary open-angle glaucoma and nonglaucomatous normals. Ophthalmology. 1984;91:564–79.
Liton PB, Luna C, Challa P, Epstein DL, Gonzalez P. Genome-wide expression profile of human trabecular meshwork cultured cells, nonglaucomatous and primary open angle glaucoma tissue. Mol Vis. 2006;12:774–90.
Du Y, Yun H, Yang E, Schuman JS. Stem cells from trabecular meshwork home to TM tissue in vivo. Invest Ophthalmol Vis Sci. 2013;54:1450–9.
Manuguerra-Gagné R, Boulos PR, Ammar A, Leblond FA, Krosl G, Pichette V, et al. Transplantation of mesenchymal stem cells promotes tissue regeneration in a glaucoma model through laser-induced paracrine factor secretion and progenitor cell recruitment. Stem Cells. 2013;31:1136–48.
Abu-Hassan DW, Li X, Ryan EI, Acott TS, Kelley MJ. Induced pluripotent stem cells restore function in a human cell loss model of open-angle glaucoma. Stem Cells. 2015;33:751–61.
Zhu W, Gramlich OW, Laboissonniere L, Jain A, Sheffield VC, Trimarchi JM, et al. Transplantation of iPSC-derived TM cells rescues glaucoma phenotypes in vivo. Proc Natl Acad Sci U S A. 2016;113:E3492–500.
Tay CY, Sathiyanathan P, Chu SWL, Stanton LW, Wong TT. Identification and characterization of mesenchymal stem cells derived from the trabecular meshwork of the human eye. Stem Cells Dev. 2012;21:1381–90.
Du Y, Roh DS, Mann MM, Funderburgh ML, Funderburgh JL, Schuman JS. Multipotent stem cells from trabecular meshwork become phagocytic TM cells. Invest Ophthalmol Vis Sci. 2012;53:1566–75.
Dominici M, Le Blanc K, Mueller I, Slaper-Cortenbach I, Marini F, Krause D, et al. Minimal criteria for defining multipotent mesenchymal stromal cells. The International Society for Cellular Therapy position statement. Cytotherapy. 2006;8:315–7.
Stamer WD, Seftor RE, Snyder RW, Regan JW. Cultured human trabecular meshwork cells express aquaporin-1 water channels. Curr Eye Res. 1995;14:1095–100.
Liton PB, Liu X, Stamer WD, Challa P, Epstein DL, Gonzalez P. Specific targeting of gene expression to a subset of human trabecular meshwork cells using the chitinase 3-like 1 promoter. Invest Ophthalmol Vis Sci. 2005;46:183–90.
Gonzalez P, Epstein DL, Borrás T. Characterization of gene expression in human trabecular meshwork using single-pass sequencing of 1060 clones. Invest Ophthalmol Vis Sci. 2000;41:3678–93.
Tomarev SI, Wistow G, Raymond V, Dubois S, Malyukova I. Gene expression profile of the human trabecular meshwork: NEIBank sequence tag analysis. Invest Ophthalmol Vis Sci. 2003;44:2588–96.
Gonzalez P, Caballero M, Liton PB, Stamer WD, Epstein DL. Expression analysis of the matrix GLA protein and VE-cadherin gene promoters in the outflow pathway. Invest Ophthalmol Vis Sci. 2004;45:1389–95.
Gasiorowski JZ, Russell P. Biological properties of trabecular meshwork cells. Exp Eye Res. 2009;88:671–5.
Polansky JR, Kurtz RM, Alvarado JA, Weinreb RN, Mitchell MD. Eicosanoid production and glucocorticoid regulatory mechanisms in cultured human trabecular meshwork cells. Prog Clin Biol Res. 1989;312:113–38.
Kubota R, Noda S, Wang Y, Minoshima S, Asakawa S, Kudoh J, et al. A novel myosin-like protein (myocilin) expressed in the connecting cilium of the photoreceptor: molecular cloning, tissue expression, and chromosomal mapping. Genomics. 1997;41:360–9.
Ortego J, Escribano J, Coca-Prados M. Cloning and characterization of subtracted cDNAs from a human ciliary body library encoding TIGR, a protein involved in juvenile open angle glaucoma with homology to myosin and olfactomedin. FEBS Lett. 1997;413:349–53.
Challa et al. IVOS 2003; 44:ARVO E-Abstract 3164.
Wang N, Chintala SK, Fini ME, Schuman JS. Activation of a tissue-specific stress response in the aqueous outflow pathway of the eye defines the glaucoma disease phenotype. Nat Med. 2001;7:304–9.
Tripathi RC, Tripathi BJ. Human trabecular endothelium, corneal endothelium, keratocytes, and scleral fibroblasts in primary cell culture. A comparative study of growth characteristics, morphology, and phagocytic activity by light and scanning electron microscopy. Exp Eye Res. 1982;35:611–24.
Koressaar T, Remm M. Enhancements and modifications of primer design program Primer3. Bioinformatics. 2007;23:1289–91.
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3--new capabilities and interfaces. Nucleic Acids Res. 2012;40:e115.
Wang X, Seed B. A PCR primer bank for quantitative gene expression analysis. Nucleic Acids Res. 2003;31:154e.
Liu Y, Munro D, Layfield D, Dellinger A, Walter J, Peterson K, et al. Serial analysis of gene expression (SAGE) in normal human trabecular meshwork. Mol Vis. 2011;17:885–93.
Fautsch MP, Johnson DH. Characterization of myocilin-myocilin interactions. Invest Ophthalmol Vis Sci. 2001;42:2324–31.
Karali A, Russell P, Stefani FH, Tamm ER. Localization of myocilin/trabecular meshwork--inducible glucocorticoid response protein in the human eye. Invest Ophthalmol Vis Sci. 2000;41:729–40.
Rao PV, Allingham RR, Epstein DL. TIGR/myocilin in human aqueous humor. Exp Eye Res. 2000;71:637–41.
Russell P, Tamm ER, Grehn FJ, Picht G, Johnson M. The presence and properties of myocilin in the aqueous humor. Invest Ophthalmol Vis Sci. 2001;42:983–6.
Niki D, Katsu K, Yokouchi Y. Ontogeny of angiopoietin-like protein 1, 2, 3, 4, 5, and 7 genes during chick embryonic development. Dev Growth Differ. 2009;51:821–32.
Gonzalez P, Zigler JSJ, Epstein DL, Borrás T. Identification and isolation of differentially expressed genes from very small tissue samples. Biotechniques. 1999;26:884–6. 888–92.
Denhardt DT, Giachelli CM, Rittling SR. Role of osteopontin in cellular signaling and toxicant injury. Annu Rev Pharmacol Toxicol. 2001;41:723–49.
Denhardt DT, Noda M, O’Regan AW, Pavlin D, Berman JS. Osteopontin as a means to cope with environmental insults: regulation of inflammation, tissue remodeling, and cell survival. J Clin Invest. 2001;107:1055–61.
Xue W, Comes N, Borrás T. Presence of an established calcification marker in trabecular meshwork tissue of glaucoma donors. Invest Ophthalmol Vis Sci. 2007;48:3184–94.
Yardley N, García-Castro MI. FGF signaling transforms non-neural ectoderm into neural crest. Dev Biol. 2012;372:166–77.
Wordinger RJ, Clark AF, Agarwal R, Lambert W, McNatt L, Wilson SE, et al. Cultured human trabecular meshwork cells express functional growth factor receptors. Invest Ophthalmol Vis Sci. 1998;39:1575–89.
Di Palma F, Pellegrino R, Noben-Trauth K. Genomic structure, alternative splice forms and normal and mutant alleles of cadherin 23 (Cdh23). Gene. 2001;281:31–41.
Bork JM, Peters LM, Riazuddin S, Bernstein SL, Ahmed ZM, Ness SL, et al. Usher syndrome 1D and nonsyndromic autosomal recessive deafness DFNB12 are caused by allelic mutations of the novel cadherin-like gene CDH23. Am J Hum Genet. 2001;68:26–37.
Siemens J, Lillo C, Dumont RA, Reynolds A, Williams DS, Gillespie PG, et al. Cadherin 23 is a component of the tip link in hair-cell stereocilia. Nature. 2004;428:950–5.
Söllner C, Rauch GJ, Siemens J, Geisler R, Schuster SC, Müller U, et al. Mutations in cadherin 23 affect tip links in zebrafish sensory hair cells. Nature. 2000;428:955–9.
Rettig J, Heinemann SH, Wunder F, Lorra C, Parcej DN, Dolly JO, et al. Inactivation properties of voltage-gated K+ channels altered by presence of beta-subunit. Nature. 1994;369:289–94.
Nakahira K, Shi G, Rhodes KJ, Trimmer JS. Selective interaction of voltage-gated K+ channel beta-subunits with alpha-subunits. J Biol Chem. 1996;271:7084–9.
Leicher T, Roeper J, Weber K, Wang X, Pongs O. Structural and functional characterization of human potassium channel subunit beta 1 (KCNA1B). Neuropharmacology. 1997;35:787–95.
Maier MM, Gessler M. Comparative analysis of the human and mouse Hey1 promoter: Hey genes are new Notch target genes. Biochem Biophys Res Commun. 2000;275:652–60.
Nakagawa O, McFadden DG, Nakagawa M, Yanagisawa H, Hu T, Srivastava D, et al. Members of the HRT family of basic helix-loop-helix proteins act as transcriptional repressors downstream of Notch signaling. Proc Natl Acad Sci U S A. 2000;97:13655–60.
Leimeister C, Externbrink A, Klamt B, Gessler M. Hey genes: a novel subfamily of hairy- and Enhancer of split related genes specifically expressed during mouse embryogenesis. Mech Dev. 1999;85:173–7.
Chin MT, Maemura K, Fukumoto S, Jain MK, Layne MD, Watanabe M, et al. Cardiovascular basic helix loop helix factor 1, a novel transcriptional repressor expressed preferentially in the developing and adult cardiovascular system. J Biol Chem. 2000;275:6381–7.
Kokubo H, Miyagawa-Tomita S, Nakazawa M, Saga Y, Johnson RL. Mouse hesr1 and hesr2 genes are redundantly required to mediate Notch signaling in the developing cardiovascular system. Dev Biol. 2005;278:301–9.
Sakamoto M, Hirata H, Ohtsuka T, Bessho Y, Kageyama R. The basic helix-loop-helix genes Hesr1/Hey1 and Hesr2/Hey2 regulate maintenance of neural precursor cells in the brain. J Biol Chem. 2003;278:44808–15.
De Jong DS, Vaes BL, Dechering KJ, Feijen A, Hendriks JM, Wehrens R, et al. Identification of novel regulators associated with early-phase osteoblast differentiation. J Bone Min Res. 2004;19:947–58.
Wordinger RJ, Lambert W, Agarwal R, Talati M, Clark AF. Human trabecular meshwork cells secrete neurotrophins and express neurotrophin receptors (Trk). Invest Ophthalmol Vis Sci. 2000;41:3833–41.
Nowak A, Szaflik JP, Gacek M, Przybylowska-Sygut K, Kamińska A, Szaflik J, et al. BDNF and HSP gene polymorphisms and their influence on the progression of primary open-angle glaucoma in a Polish population. Arch Med Sci. 2014;10:1206–13.
The authors thank Dr Jodhbir Singh Mehta and Dr Gary Peh from the Singapore Eye Research Institute for their generous provision of some samples utilized in this study. The authors also thank Dr Vibhor Kumar for his valuable input for the microarray analyses.
The work was supported by the Agency for Science, Technology and Research (A*STAR), Singapore.
Availability of data and materials
Our data has been deposited into Gene Expression Omnibus.
PS: Conception and design, collection and/or assembly of data, data analysis and interpretation, manuscript writing. CYT: Collection and/or assembly of data. LS: Conception and design, data analysis, financial support, manuscript writing. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All human samples were provided by the Lions Eye Institute for Transplantation and Research which attained the tissues following informed written consent from donors and/or the next of kin of the deceased donors. The work was performed according to institutional ethics guidelines.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Data on the differential gene lists TM vs TM-MSC, TM vs cornea and TM vs sclera. (XLSX 721 kb)
List of primer sequences used in the study. Table S2. Details of antibodies utilized for immunofluorescence. (DOCX 13 kb)
Immunofluorescence of TM markers on TM-MSC. (TIF 4091 kb)
qPCR analysis of identified marker genes in TM cells (TMC) relative to TM-MSC. *** P < 0.001; ** P < 0.01; * P < 0.05 (Student’s t-test). (TIF 678 kb)
About this article
Cite this article
Sathiyanathan, P., Tay, C.Y. & Stanton, L.W. Transcriptome analysis for the identification of cellular markers related to trabecular meshwork differentiation. BMC Genomics 18, 383 (2017). https://doi.org/10.1186/s12864-017-3758-7
- Human trabecular meshwork