- Research article
- Open Access
- Open Peer Review
Distinct gene subsets in pterygia formation and recurrence: dissecting complex biological phenomenon using genome wide expression data
BMC Medical Genomics volume 2, Article number: 14 (2009)
Pterygium is a common ocular surface disease characterized by fibrovascular invasion of the cornea and is sight-threatening due to astigmatism, tear film disturbance, or occlusion of the visual axis. However, the mechanisms for formation and post-surgical recurrence of pterygium are not understood, and a valid animal model does not exist. Here, we investigated the possible mechanisms of pterygium pathogenesis and recurrence.
First we performed a genome wide expression analysis (human Affymetrix Genechip, >22000 genes) with principal component analysis and clustering techniques, and validated expression of key molecules with PCR. The controls for this study were the un-involved conjunctival tissue of the same eye obtained during the surgical resection of the lesions. Interesting molecules were further investigated with immunohistochemistry, Western blots, and comparison with tear proteins from pterygium patients.
Principal component analysis in pterygium indicated a signature of matrix-related structural proteins, including fibronectin-1 (both splice-forms), collagen-1A2, keratin-12 and small proline rich protein-1. Immunofluorescence showed strong expression of keratin-6A in all layers, especially the superficial layers, of pterygium epithelium, but absent in the control, with up-regulation and nuclear accumulation of the cell adhesion molecule CD24 in the pterygium epithelium. Western blot shows increased protein expression of beta-microseminoprotein, a protein up-regulated in human cutaneous squamous cell carcinoma. Gene products of 22 up-regulated genes in pterygium have also been found by us in human tears using nano-electrospray-liquid chromatography/mass spectrometry after pterygium surgery. Recurrent disease was associated with up-regulation of sialophorin, a negative regulator of cell adhesion, and never in mitosis a-5, known to be involved in cell motility.
Aberrant wound healing is therefore a key process in this disease, and strategies in wound remodeling may be appropriate in halting pterygium or its recurrence. For patients demonstrating a profile of 'recurrence', it may be necessary to manage as a poorer prognostic case and perhaps, more adjunctive treatment after resection of the primary lesion.
Global gene expression has been used successfully to elicit biological behavior in different soft tissue tumors. Pterygium as a human disease, noted to be more prevalent than 20% of some populations,[2, 3] is of immense biological interest for a few reasons.
First, the pathogenesis of this condition is hotly debated. Hypothesis driven approaches have not resolved the relative importance of competing mechanisms for this disease. Theories that have been proposed include inflammatory influence, degeneration of connective tissue, genetic instability, angiogenesis, redox-related toxicity, cellular proliferation, aberration of apoptosis, exuberant wound healing, altered lipid metabolism, mast cell infiltration. and stem cell dysfunction. Conventional approaches to disease mechanism, by virtue of their narrow focus, were not helpful to assess relative contribution of widely heterogenous processes. Furthermore, a fundamental issue about the diseased tissue remains un-resolved in this context: the origin of the epithelium overlying pterygial lesions, though suspected to be conjunctival in origin, is not entirely certain.
Second, the spectrum of tumor size and behavior is tremendous, ranging from the inconspicuous lesion barely encroaching on the peripheral cornea, to the rapidly growing, menacing tumors that obscure the visual axes and threaten vision. In some contexts, additional molecular events. may drive the original tumor to behave differently from the original lesion. It is also intriguing that unlike frank malignancies, these lesions, however fast growing, do not erode through the full thickness cornea, highlighting the presence of distinct processes from those manifested in malignant tumors.
Third, pterygium is a mixed soft tissue tumor that is strongly associated with non-ionising ultraviolet radiation. Unlike the case of cutaneous melanomas, the cell type that responded to the environmental trigger may not be epithelial but rather, the fibro-vascular component – a truly unusual phenomenon in human tumor biology because ultraviolet radiation-induced tumors encountered in humans are epithelial in origin and ultraviolet radiation effects on keloids, if any, are generally inhibitory.
Lastly, the inability of researchers to discover an analogous tumor in animals or reconstitute the disease in organ cultures imply that pterygium involves specie-specific mechanisms mediated by in-vivo cell-cell or cell-matrix interactions. The molecules involved in this disease may have therefore arisen from divergent evolution in the human ocular surface.
In view of the controversies on the multiple mechanisms of pterygia formation, we advocate an unbiased, global gene expression approach to decipher the dominant gene expression patterns that may underlie specific molecular events of interest. Recently, the genome-wide microarray data in whole tissue pterygium as well as microarray data limited to cultured pterygial fibroblasts, have been published. We have performed a study using gene microarray on recurrent pterygium, primary pterygium and un-involved conjunctiva, examined the gene subsets which were differentially regulated, compared our findings with these previous microarray studies,[16, 17] and discussed some of the findings in the light of a complementary proteomics approach.
In clinical practice, the treatment of this condition is surgical excision. However, some cases aggressively recur after surgery. Conjunctival auto-grafting as an adjunctive procedure may reduce recurrence, though the explanation for this is not entirely clear. For these reasons, we speculate that a non-hypothesis driven approach may also be useful to discover gene expression signatures that predispose lesions to a more aggressive phenotype. Such information will benefit clinicians, who can appropriately anticipate otherwise 'un-expected' biological behavior in their treatment of this disorder.
Antibodies against MUC5AC, keratin 6 and CD24 were from Labvision, Neomarkers, Fremont, CA, USA, anti-MSMB antibodies were from US Biologicals, Massachusetts, MA, USA, anti-CEACAM5 and anti-NR4A2 antibodies were from Abcam, Cambridge, UK.
Samples used for the Study
The procurement and use of both human tissues in this study was in compliance with the tenets of the Declaration of Helsinki. The study was approved by the Institutional Review Board of Singapore Eye Research Institute. Written informed consent was obtained from donors after explanation of the nature and possible consequences of the study. Human tissues samples were obtained from patients diagnosed with primary pterygium and came from different races, Chinese, Malay and Indian. All patients underwent pterygium excision in conjunction with the use of an upper bulbar conjunctival free autograft placed over the site of the original lesion. All pterygia specimens used were nasal pterygia. The whole pterygium tissue and a small portion of the conjunctival patch (approx. 1 × 3 mm) from the superotemporal conjunctiva were collected and were rapidly frozen in liquid nitrogen after removal and stored at -150°C. For immunohistology purpose, the samples were collected on ice and embedded in Optimal Cutting Temperature (OCT, Sakura, USA) in the laboratory.
The procedure for the microarray experiment has previously been published. All microarray chips and related protocols and equipment for the processing of these chips were from Affymetrix Inc., Santa Clara, CA. The human genome GeneChip U133A consisting of more than 22000 probe sets was used for this study.
A group of 8 primary pterygium samples, harvested from 4 males (aged 40 to 50 years old) and 4 females (aged 50 to 60 years old) was used in this experiment. Another group of 4 control conjunctiva tissues, each pooled from 4 individual samples, was used as controls to the diseased tissues. Considering that only a tiny piece of the conjunctiva tissue could be obtained from each patient during the surgery, pooling of 4 conjunctiva tissues was necessary to obtain enough starting material.
Total RNA was extracted using TRIzol Reagent (Invitrogen, CA) and purified with RNeasy Mini Kit (Qiagen, Valencia, CA) according to the manufacturer's instructions. Five micrograms of each purified RNA sample were prepared according to the Affymetrix standard protocol. Fifteen micrograms of biotin-labelled cRNA using BioArray RNA Transcript Labelling Kit (ENZO Life Sciences, NY) were fragmented and the appropriate volume injected separately into the probe array chips.
The transcripts were hybridized onto the immobilized oligonucleotide sequence on array for 16 hours at 45°C under 60 rpm rotation using GeneChip Hybridization Oven 640. Washing and Strepavidin-staining steps were performed using Affymetrix Fluidics Station 450. The chips were scanned using GeneChip Scanner3000 and the image data were further analyzed using Microarray Suite v.5.0. Data analysis included pre-processing and normalization before identification of differentially expressed genes. Pre-processing adjusted for non-specific binding and background noise, whereas normalization removed systematic variation in the data due to effects other than biological differences. We used the Robust Multi-array Average (RMA) model  to extract the gene expression signals from probe intensities without taking the mismatch probe signals into consideration. Cross-array normalization was performed using the intensity-based log ratio median method  with the first array as the reference. Gene-level normalization was performed by normalizing all samples to the median of the expression level of the control (un-involved conjunctiva) samples. The data were annotated using gene annotation headings, bioprocesses and molecular functions from the NetsAffx database http://www.affymetrix.com/analysis/index.affx.
Analysis of Microarray Data
The data reported in this study have been deposited in NCBI's Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/ with GEO series accession number GSE2513). Data were visualised and explored using the GeneSpring GX 7.3 platform (Agilent Technology, Redwood City, CA).
For selection of differentially expressed genes, the modified t-statistic (SAM)  with 100% of the standard deviation percentile as the fudge constant was used. The threshold for significantly changed genes was set at a false discovery rate (FDR) of 5%. After delineating a list of differentially expressed genes, we performed various types of analysis to search for a pattern of global gene expression. The methods ranged from K-means clustering, to the construction of hierarchical dendrograms on a subset of genes, as well as the visualisation of gene expression pattern summarised by principal component analysis. In additional to analysis of differentially expressed genes, we employed the Gene Set Enrichment Analysis (GSEA) method for identification of those pathways that were more affected in pterygium tissues. The goal of GSEA is to determine whether members of a gene set tend to occur toward the top or bottom of the list, in which case the gene set is correlated with the phenotypic class distinction. In contrast to the methods for extraction of differentially expressed genes, GSEA considers the collective up- or down-regulation of a gene set rather than individual genes. The GSEA software package downloaded from the Board Institute's website was used for the identification of activated or deactivated pathways in primary pterygial and conjunctival tissues. Microarray data for 8 primary pterygial samples and 4 samples of conjunctival tissue were used. Various datasets from publicly accessible databases (Biocarta, STKE, PubMed and KEGG) were dissected into over 300 gene sets. The two phenotypic classes used were pterygium and uninvolved conjunctiva tissues.
As an independent method of analyzing signaling pathways in pterygium, we exported the list of genes significantly up-regulated or down-regulated by 2 fold into the Pathway Studio 5.0 software (Ariadne Genomics Inc, Rockville, MD). The option 'Find all shortest paths between selected entities' was used on the up-regulated and down-regulated genes sequentially. The number of connectivities was limited to 2 and the genes mapped into known pathways displayed on a chart, and the relevant relationships listed in a table. The purpose of performing this analysis was to identify important upstream regulators of the differentially expressed genes as well as downstream effectors, in an unbiased fashion, based on known biological knowledge.
Methods used to extract and analyse tear proteins from patients who had undergone surgery for pterygium have already been described. Wherever possible, the fold change data for genes corresponding to detected tear proteins were tabulated.
Real Time Reverse Transcription Polymerase Chain Reaction
Six pairs of independent pterygium and uninvolved conjunctiva tissues were used in the real time reverse transcription polymerase chain reaction (qPCR) experiment. One microgram of the each total RNA preparation was reverse-transcribed to single stranded cDNAs using an oligo-dT primer with Superscript II Rnase H Reverse Transcriptase (Invitrogen, USA). Primer pairs specific to each gene used in the qPCR are listed in Additional file 1. qPCR was performed using SYBR Green PCR Master Mix (Applied Biosystems) in an ABI Prism 7700 Sequence Detection System (Applied Biosystems). The thermal cycling conditions were as follows: 95°C for 10 min, 45 cycles at 95°C for 30 s and 60°C for 1 min. Data obtained from qPCR were analyzed using the comparative CT method as previously described by Livak. Paired t-statistics with p-value < 0.05 was used to determine whether the qPCR results in pterygial tissue were significantly different from those in un-involved conjunctiva.
A separate group of paired pterygial and conjunctival tissue samples from three patients were used for immunohistochemistry. Sections were cut at 5 μm thicknesses from blocks of freshly frozen pterygium and matched un-involved conjunctiva embedded in the OCT. The post-fixed sections were incubated with the primary antibody at 4°C overnight. Specific antibodies against the following were used: MUC5AC, MSMB, CD24, CEACAM5, keratin 6 and NR4A2. After 3 washes with PBS, secondary antibody used was either fluorescein isothiocyanate (FITC)- or Rhodamine- conjugated anti IgG (Santa Cruz, USA), incubated with the sections for 40 minutes at room temperature. After 3 final washes each section was mounted with in a fluorescence mounting medium (DAKO, Denmark). Images were captured using a 40× Achrostigmat lens on an Axioplan2 microscope equipped with an AxioCam MR camera (Carl Zeiss, Germany).
Global Gene Expression Profiling and analysis of pterygium versus conjunctiva
Gene expression analysis
With gene level normalization using un-involved conjunctival samples as reference, gene expression analysis showed that a total of 105 probe sets were significantly changed (p < 0.05) by at least 2 fold in primary pterygial tissue from uninvolved conjunctiva. Among these, 60 probe sets were up-regulated and 45 probe sets were down-regulated (Table 1). Without performing gene level normalization, 114 unique genes out of 156 probe sets were significantly changed in pterygium. (Figure 1)
Gene Ontology Analysis
All significantly changed genes were categorized based on biological functions using categories from the GeneOntology database (Figure 2). Among the up-regulated genes, genes coding for cell adhesion (8), extracellular matrix (ECM) (7) and structural proteins (8) accounted for 38% of up-regulated genes. Examples of up-regulated cell adhesion molecules include CADH11, MCAM, FN1, CEACAM5, ECM component include COL1A2, COL1A1, COL4A2, COL6A3 and structural proteins, ACTA2, KRT12, KRT6, KRT23, SPRR1B and SPRR3. Several of these genes also encode for proteins involved in wound healing, including COL6A3, CD24, FN1, KRT6a, TFF1, SPRR1b, and MUC5AC. Another group of up-regulated genes is the 'mitogenic proteins' group exemplified by MSMB, CEACAM5 and CD24.
Some immune response and transcription factors were down-regulated in primary pterygium. Examples of the former include immunoglobulin heavy chain genes, whereas examples of the latter include ATF3, BTG2, ERG1, FOS, FOSB, JUN, NR4A1 and NR4A2. Other down-regulated genes include those encoding for transport proteins, ie., HBB and HBD, and those involved in stress-response, including DUSP1 and GADD45B.
When the composition of the list of significantly changed genes was studied in terms of certain gene ontology categories, we observed that there is an over-representation of such genes compared to their expected frequency, considering the proportion of the probes representing these genes present on the chip (Table 2). This suggests that structural molecules and molecules with transporter activity play a role in pterygium.
Analysis of enriched pathways
Gene set enrichment analysis (GSEA) was performed using gene expression data from the conjunctival and primary pterygium groups (Figure 3). Significantly enriched pathways in pterygium were related to ECM formation, tissue invasion, angiogenesis and lipid metabolism. Candidate genes found to be involved in the ECM pathway were cadherins (CDH5 and CDH3), integrins (eg., ITGb1), cystatin C (CST3) and intracellular adhesion molecule 1(ICAM1). Angiogenesis related genes include akt1, akt2, tyrosine kinase 2 (Tyk2), phosphinositide-3-kinases (PIK3C2A, PIK3CD, PIK3R1, PIK3R2, PIK3R4) and janus kinases (JAK1, JAK2). The genes involved in 'tissue invasion' pathways include plasminogen activator and receptor (PLAU, PLAUR), matrix metalloproteinase (MMP1,2,9), tissue inhibitor of metalloproteinase 1 (TIMP1, TIMP3), serine (or cysteine) proteinase inhibitor clade B member 5 (SERPINB5), s100 calcium binding protein (S100A4), secreted phosphoprotein1 (spp1) and mucin (MUC1). The lipid formation pathway was also significantly enriched in pterygium tissues, for example, the sterol regulatory element binding transcription factor1 (SREBF1), dual specific phosphatase 14 (DUSP14), hexokinase2 (HK2) and the solute carrier family 2 member 4 (SLC2A4).
The only significantly enriched pathway in un-involved conjunctiva was related to apoptosis. The genes in this pathway include the calpains (CAPN1, CAPN2, CAPN3, CAPNS1) and bcl-antagonist of cell death (BAD2).
Pathway analysis to identify upstream and downstream signals
Up-regulated genes and down-regulated genes were analysed in pathway studio. Figure 4 summarised signaling network that resulted from this analysis. Primary pterygium may be characterized by stress-induced down-regulation of transcription factors (Egr1, Jun and Fos), with defective wound healing as the major process responsible for the disease phenotype. Scrutiny of the networks in Figure 4 illustrates a few other relevant processes: up-regulation of PECAM1 and down-regulation of TFPI-2 may contribute to aberrant vascularisation, whereas down-regulation of DUSP1 and up-regulation of GP75 may contribute to an abnormal response to oxidative stress. The identities and biological/molecular functions of these relationships are shown in Additional files 2 and 3.
To extract some useful knowledge or patterns from global gene expression data, we performed a few cluster analyses. The K-means cluster method was computed using 5 clusters on the 156 gene probes differentially expressed in pterygium compared to controls, over 1000 iterations. Table 3 shows the final characteristics of the clusters. Two clusters contained genes that were down-regulated in pterygium, whereas three other clusters showed genes that were up-regulated. Visual examination of the role of these genes revealed biological functions that are consistent with the above-mentioned pathway and gene ontology studies.
To identify co-expressed gene sets or the similarity of samples, we further performed hierarchical clustering analysis on the data for the 156 gene probes. Two closely related clusters of genes up-regulated in conjunctiva was detected (Figure 5A), and similarly, 2 clusters of up-regulated genes in pterygium were detected (Figure 5B). Scrutiny of the composition of these clusters yields a number of cytoskeletal proteins, immunoglobulins, cancer markers and transcription factors.
We further performed principal component analysis on the entire set of genes in the chip, as well as several selected categories of Gene Ontology previously implicated in the pathogenesis of pterygium (Figure 6) for the following two purposes: 1. To show the compactness of the samples (represented by the distance between circles) in each condition and 2. To reveal the possible existence of a plane which can separate the 2 conditions. The 3D scatter plots of the first three prinicipal components show that the expression of genes coding for ECM component (Figure 6J) was able to differentiate the conjunctiva and pterygium samples better than those in the antioxidant (Figure 6F) and apoptosis (Figure 6G) categories. When all the genes were included in the analysis (Figure 6A), a clear plane between the pterygial and conjunctival samples was not so evident. This may be due to the inclusion of many genes that were not significantly regulated but whose expression data contained extensive noise. For this analysis, we deliberately included processes like 'reproduction' (Figure 6H) being unlikely to play a significant role in pterygium pathogenesis, acted as a negative control.
Differences between primary and recurrent pterygium
Recurrent pterygium [see Additional file 4] differs from primary cases as it has a different morphology and often worse prognosis. Table 4 shows the genes that were differentially regulated in recurrent compared to primary pterygium. Additional file 5 shows that a large number of genes in recurrent pterygium were down-regulated or up-regulated relative to primary pterygium and conjunctiva tissues. Since there were a rather large number of differentially regulated genes, only those with at least 2 fold change between primary and recurrent pterygia are shown in table 4. It is interesting that the up-regulated genes are not the same as those up-regulated in primary pterygium compared to conjunctiva. Examples include stearoyl-CoA desaturase 5 that converts fatty acids into monounsaturated forms in the endoplasmic reticulum and involved in dyslipidemia, the ubiquinol-cytochrome c reductase involved in the mitochondrial respiratory chain, the DNA recombination repair protein RAD51, the neuronal thread protein (AD7C-NTP), which is a extracellular protein involved in apoptosis, and the gene for Mediterranean fever (MEFV) which is involved in regulation of transcription and inflammatory response. Involvement of the gene for sialophorin (SPN or CD43) and NEK5 (Never in mitosis gene a) are consistent with increased cell migration because SPN is a negative regulator of cell adhesion and NEK5 is involved in microtubule function and motility. Additional file 6 shows possible mechanisms of pterygium recurrence involving these mediators.
Validation of gene microarray data
Relative quantitative real time polymerase chain reaction (qPCR) was performed as an independent laboratory approach to validate the transcriptional level changes in the microarray experiment. These results (Figure 7) show that for all genes selected for validation with qPCR, the direction and magnitude of changes were consistent with the results obtained from the microarray analysis.
We then addressed the question whether protein levels were also affected in pterygium. Several proteins were chosen arbitrarily to address the protein expression levels as well as the localization in pterygium tissues. Our results (Figure 8) show that protein expression levels for keratin 6, CEACAM5, CD24, MSMB and MUC5AC were increased, whereas NR4A2 and IGFBP3 were reduced in pterygium, consistent with the transcript changes detected by the microarray analysis. Furthermore, the immunofluorescent staining yielded interesting information about the localization of these proteins in pterygium. Immunofluorescent staining (Figure 8A) shows that keratin 6a stained strongly in the superficial layer in pterygial epithelium but not in un-involved conjunctiva. Similarly, CEACAM5 was detected in the squamous layer of pterygial epithelium but was not detectable in conjunctiva. There was up-regulation and nuclear accummulation of the cell adhesion molecule CD24 in the pterygium epithelium, but CD24 was not detectable in conjunctiva. MSMB, on the other hand, was detected prominently in the basal epithelial layer of pterygium, with some staining also in the superficial stromal cells adjacent to the basal epithelia. Lastly, MUC5AC protein was detected in both pterygial and conjunctival tissue sections predominantly in goblet cells. However, the intensity of MUC5AC staining in pterygium was stronger than that in conjunctiva. IGFBP3, present in un-involved conjunctival stroma and epithelium, and NR4A2, present in conjunctival epithelium, were reduced in pterygium (Figure 8B).
Comparison of microarray data with tear proteins
Another report from our center has previously documented the analysis of tear proteins in tears of patients after surgery for pterygium. Briefly, this study employed reverse-phase high-pressure liquid chromatography followed by tryptic digestion and characterisation of proteins using nanoLC-nano-ESI-MS/MS. Eleven genes corresponding to the detected tear proteins were significantly up-regulated in pterygium compared to un-involved conjunctival tissue (Table 5). Examples of the up-regulated genes were those encoding for prolactin induced protein I and the S100 A8 proteins. Another eleven genes corresponding to the detected tear proteins were found to be down-regulated in pterygium. These include lipophilin C, ribonuclease 4, complement C3 and histone 1b. Twenty-nine genes corresponding to detected tear proteins were not significantly up or down-regulated in pterygium. Since there were no controls in this study, it is difficult to interpret some of these findings. Furthermore it was not possible to compare the microarray results with every instance of the protein data because more than one type of probe sets may contribute to the synthesis of the protein and not all gene probe sets have signal-noise ratio high enough for this analysis.
A global gene expression analysis of pterygium showed distinct differences between primary pterygium and uninvolved conjunctiva. Several pathways were significantly affected in pterygium. These were: increase in the production of extracellular matrix, structural proteins, mitotic proteins, and protein involved in tissue invasion.
Recurrent pterygia demonstrated a different signature composed of other perturbed genes. For example, the COL4A6 (AL031177) and the RAB6B (BC002510) were down- and up-regulated respectively in recurrent pterygia compared to primary pterygia. The former encodes for one of the 6 subunits of collagen IV, a major component of the basement membrane, whereas the latter is a RAS family oncogene. The finding suggests that recurrence is a distinct biological phenomenon from the formation of primary pterygium, even though in general, we did not detect obvious microscopic changes between primary and recurrent pterygia.
Microarray analysis of gene expression is a useful approach for understanding the molecular mechanism of disease. We used both qPCR and immunohistochemistry to show that both RNA level and protein levels of specific mediators were dysregulated, validating the results from the gene microarray approach. Some 22 genes corresponding to tear proteins detected in pterygium subjects were significantly up or down-regulated in pterygium relative to conjunctival tissue. Many of the processes discovered in this study are highly novel and were not previously associated with pterygium, for example, immunohistochemistry show that cell adhesion molecules (ie., CD24) may be increased but abnormally localized in the nuclei in pterygium epithelium and therefore cell adhesion properties may be disturbed. Another example is the evidence for Goblet cell dysfunction. In pterygium, there was elevation of transcript and protein expression of MUC5AC and MSMB. MUC5AC is a well-known marker of conjunctival mucous-secreting Goblet cells, and MSMB is a cutaneous squamous cell carcinoma related protein that is associated with respiratory tract Goblet cells.
Comparison with previous studies
The data from a previous study utilized only 2 primary and 1 recurrent pterygia, with no stratification of data between the clinical sub-types. This study highlighted the up-regulation of 29 genes common to primary and recurrent pterygium. Our results supported the dysregulation of 19 of these genes in pterygium (Table 6), for example, TRAP100, MIP-4, RBP-1, MAP-17 and PECAM1. One apparent discrepancy was that the CLIC2 was up-regulated in our study, but significantly down-regulated in John-Aryankalayil et al. Discrepancies reflect differences in study methodology, for example, this study used an older microarray chip HG_U95Av2, containing considerably less probes than the U133A; and the earlier study may have differences in probe, chip or gene-level normalizations which were not reported. Differences in normalization may account for the differences in the folds of change in specific genes between different studies. In John-Aryankalayil et al's study, although a few up-regulated genes from different functional categories were tabulated, there was no attempt to evaluate the patterns of expression or relative contribution of broad biological processes. Unlike this study, we did not restrict to listing differentially expressed genes by fold change.
Previous studies[28, 29] have shown that gene expression profiles in tumor and wound response were similar, lending support to our opinion that wound response is the major theme in global pterygial gene expression.
Possible mechanisms of pterygia formation
We believe that our data support a predominantly wound healing pattern of gene expression in pterygium. Genes encoding for extracellular matrix, structural and adhesion molecules including wound healing related proteins, collagen subtypes, keratin 6A and fibronectin were significantly up-regulated in pterygium. Other up-regulated proteins like small proline rich protein 1B (SPRR1B), CD24, S100 calcium binding protein, SPARC, TFF1, SPRR1B and SERPINB13 also govern the wound healing process and cornification of epithelium, again, supporting the hypothesis of aberrant wound response in pterygium. Fibronectin alternatively splice transcript with EDA, up-regulated in our pterygial specimens, was shown to be up-regulated in wound healing. In order to further validate this finding, semi-quantitative PCR was performed on paired pterygial and conjunctival tissue specimens from 3 patients to examine the level of fibronectin transcript with EDA. The results show that pterygium expressed higher levels of fibronectin transcript with EDA compared to uninvolved conjunctiva (Figure 8C).
Our data also suggest that in pterygia formation, there may be increased cellular proliferation and motility on one hand, and reduced cell death on the other. Genes up-regulated in pterygium include fibronectin (FN1), CEACAM5 (CEA), CD24, SPARC, MSMB and TFF1. CEACAM5 (CEA), a common cancer marker, was up-regulated 3.8-fold in pterygium. SPARC was known to have anti-adhesive property, and induced cancer cell motility. The transcripts coding for microseminoprotein (MSMB) or PSP94 was up-regulated 3.6-fold, whereas transcripts for calcium binding protein S100A8 was up-regulated 5.7-fold in pterygium. Interestingly, in our analysis of tears from patients with pterygium, S100 proteins were also detected. Trefoil factors promote restitution of epithelial cells and are abundantly secreted onto the mucosal surface rapidly after mucosal injury. In contrast to the above, genes coding for apoptosis (TGM2, IGFBP3 and DUSP1) were down-regulated in pterygium, reinforcing the over-proliferative tendency in pterygium.
The stress-inducible transcription regulator genes including ATF3, BTG2, EGR1, ERG2, FOS, JUN, NR4A1 and NR4A2 were surprisingly, down-regulated in pterygium relative to uninvolved conjunctiva. One has to bear in mind however, that our tissue specimens may represent chronic UV stimulation. The expression of these transcription factors may very well be elevated shortly after the commencement of the initial stimulus.
The transcript level of immunoglobulin subunit in pterygium was relatively lower as compared to uninvolved conjunctiva. This does not support a role for B cell mediated immunity in pterygium formation. However, the co-expression of some immunoglobulin genes (Figure 4) may have some biological significance, such as an immunological basis for the inappropriate wounding response.
Possible mechanisms of pterygia recurrence
One hypothesis is that a low level of inflammation may stimulate epithelial and fibroblast cell migration due to factors such as SPN [see Additional file 3A]. This allows continuation of inflammatory factors and cytokines to be produced, potentiating a cycle [see Additional file 3B]. Factors such as SPN are unique to recurrence only, since they were up-regulated compared to primary pterygium cases.
In an alternative hypothesis, the formation of an exuberant scar or a 'go' signal, may or may not be coupled with a reduction of 'stop' or inhibitory signals. One such 'stop signal' may be the MSMB protein. Any up-regulation of this anti-metastatic gene in pterygium could give an inhibiting signal to slow the growth rate in pterygium. When pterygium is surgically removed, this inhibiting factor may also be reduced, encouraging the remnants of fibrous tissue to proliferate more aggressively, attaining the phenotype of recurrent pterygium. This may also explain the success of conjunctival autographs as adjunctive procedure after excision of the lesion in the prevention of recurrence. The graft tissue may 'replenish' any surgically-induced loss of 'stop' signals. In support of this hypothesis, the level of MSMB was significantly depressed in recurrent compared to primary pterygia.
Strengths and Limitations
The strengths of this study include the use of a variety of analytical approaches to interpret the microarray data. The tissue specimens have been harvested in a center which manages a high volume of such disease. Samples have been processed in a very standardized fashion, obtained from patients with accurate clinical diagnoses.
One limitation of the study is that we did not evaluate p53 and related genes such as the Vasoendothelial growth factor (VEGF). Pterygia specimens contained a mixture of cell types whereas un-involved conjunctiva consisted largely of epithelial tissue. Therefore, any observations concerning differential gene expression could be related to differences in composition of cells as well as to pathological processes.
We did not notice regional differences in the intensity of CD24 immunofluorescent staining along the length of the pterygium epithelium on longitudinal section [see Additional file 7]. However, we cannot conclusively state that there is no difference in the pattern of CD24 expression between the head and the body of the pterygium. Since the limits of the pterygium body was defined clinically by the surgeon, we were uncertain if the bulk of the pterygium body was available for immunohistology.
Since we did not obtain temporally located pterygia in this study, our findings may not be applicable to this less frequently occurring type of pterygia.
Specimens analysed were at the later stage of pterygium formation when surgical removal was needed. No data at the early stage of disease were available for investigation.
In addition to these limitations, there are inherent limitations with hieuristic and even more deterministic methods of analyzing global gene expression. For example, the results of the k-means clustering depend arbitrarily on the initial number of clusters and the number of iterations. However, we attempt to reduce the effects of these shortcomings by employing a variety of different data-mining techniques and interpreting results as a whole and not in isolation. We drew conclusions based on findings that are consistent between different procedures and refrained from over-interpretation of isolated anomalies.
If changes in global gene expression in cases of recurrent pterygium do not occur until a very late time point, it may not be possible to predict recurrence of the condition by analysing the tissue obtained during the time of first excision. However, the gene expression data may still be useful for understanding the biology of recurrence. The tear protein analysis also has limitations. It is difficult to determine whether the detected proteins were physiological, related to the surgical trauma, or arising from the disease. Nevertheless the 22 significantly up-regulated genes that have products in human tears may be potential biomarkers for the disease.
This study illustrates that an unbiased global gene expression approach is useful to address the disease mechanisms in a controversial condition. Further studies on more specimens may allow the use of a sub-set of genes to prognosticate lesions. Non-surgical treatment for pterygium currently does not exist. Since wound healing and matrix dysregulation is a major theme, known modulators of wound healing may be explored in pterygium in a targeted way, using topical application. Other novel processes require further evaluation. Studies in the regulation of cell adhesion by CD24 and mucous processing pathways are required to understand pterygium formation. Studies are also required to elicit the origin of the overlying epithelium in pterygium tissue.
Nevertheless, our study suggests that after resection of the initial tumor, clinicians can modify treatment in selected patients based on objective criteria after analysing tumor gene expression or tumor proteins.
Based on differential gene profiling in pterygium, it can be concluded that an aberrant wound healing process is the major pathogenetic process in pterygial formation. The other secondary processes such as cornification and attempted barrier formation may be compensatory in nature. This contrasts sharply with frank malignant neoplasms, where gene expression signatures are dominanted by increased proliferation, reduced apoptosis, cell cycling anomalies and genomic instability. The existence of an aggressive phenotype to account for post excision recurrence may be related to an imbalance of growth signals rather than due to mere prolongation of the stimulus that initiated primary pterygial formation.
(complementary) deoxyribonucleic acid
false discovery rate
Gene expression omnibus
Gene set enrichment analysis
optimal cutting temperature
- (q) PCR:
(quantitative) polymerase chain reaction
Robust Multi-array Average
(complementary) ribonucleic acid.
West RB, Nuyten DS, Subramanian S, Nielsen TO, Corless CL, et al: Determination of stromal signatures in breast carcinoma. PLoS Biol. 2005, 3: e187-10.1371/journal.pbio.0030187.
Luthra R, Nemesure BB, Wu SY, Xie SH, Leske MC: Frequency and risk factors for pterygium in the Barbados Eye Study. Arch Ophthalmol. 2001, 119: 1827-1832.
Wu K, He M, Xu J, Li S: Pterygium in aged population in Doumen County, China. Yan Ke Xue Bao. 2002, 18: 181-184.
Di Girolamo N, Chui J, Coroneo MT, Wakefield D: Pathogenesis of pterygia: role of cytokines, growth factors, and matrix metalloproteinases. Prog Retin Eye Res. 2004, 23: 195-228. 10.1016/j.preteyeres.2004.02.002.
Dushku N, John MK, Schultz GS, Reid TW: Pterygia pathogenesis: corneal invasion by matrix metalloproteinase expressing altered limbal epithelial basal cells. Arch Ophthalmol. 2001, 119: 695-706.
Perra MT, Maxia C, Corbu A, Minerba L, Demurtas P, et al: Oxidative stress in pterygium: relationship between p53 and 8-hydroxydeoxyguanosine. Mol Vis. 2006, 12: 1136-1142.
Aspiotis M, Tsanou E, Gorezis S, Ioachim E, Skyrlas A, et al: Angiogenesis in pterygium: study of microvessel density, vascular endothelial growth factor, and thrombospondin-1. Eye. 2007, 21: 1095-1101. 10.1038/sj.eye.6702495.
Kim YJ, Park ES, Song KY, Park SC, Kim JC: Glutathione transferase (class pi) and tissue transglutaminase (Tgase C) expression in pterygia. Korean J Ophthalmol. 1998, 12: 6-13.
Tan DT, Liu YP, Sun L: Flow cytometry measurements of DNA content in primary and recurrent pterygia. Invest Ophthalmol Vis Sci. 2000, 41: 1684-1686.
Tan DT, Tang WY, Liu YP, Goh HS, Smith DR: Apoptosis and apoptosis related gene expression in normal conjunctiva and pterygium. Br J Ophthalmol. 2000, 84: 212-216. 10.1136/bjo.84.2.212.
Song YS, Ryu YH, Choi SR, Kim JC: The involvement of adult stem cells originated from bone marrow in the pathogenesis of pterygia. Yonsei Med J. 2005, 46: 687-692.
Peiretti E, Dessi S, Putzolu M, Fossarello M: Hyperexpression of low-density lipoprotein receptors and hydroxy-methylglutaryl-coenzyme A-reductase in human pinguecula and primary pterygium. Invest Ophthalmol Vis Sci. 2004, 45: 3982-3985. 10.1167/iovs.04-0176.
Nakagami T, Murakami A, Okisaka S, Ebihara N: Mast cells in pterygium: number and phenotype. Jpn J Ophthalmol. 1999, 43: 75-79. 10.1016/S0021-5155(98)00069-0.
Coroneo MT, Di Girolamo N, Wakefield D: The pathogenesis of pterygia. Curr Opin Ophthalmol. 1999, 10: 282-288. 10.1097/00055735-199908000-00011.
Ramalho FS, Maestri C, Ramalho LN, Ribeiro-Silva A, Romao E: Expression of p63 and p16 in primary and recurrent pterygia. Graefes Arch Clin Exp Ophthalmol. 2006, 244: 1310-1314. 10.1007/s00417-006-0287-5.
John-Aryankalayil M, Dushku N, Jaworski CJ, Cox CA, Schultz G, et al: Microarray and protein analysis of human pterygium. Mol Vis. 2006, 12: 55-64.
Solomon A, Grueterich M, Li DQ, Meller D, Lee SB, et al: Overexpression of Insulin-like growth factor-binding protein-2 in pterygium body fibroblasts. Invest Ophthalmol Vis Sci. 2003, 44: 573-580. 10.1167/iovs.01-1185.
Zhou L, Beuerman RW, Foo Y, Liu S, Ang LP, et al: Characterisation of human tear proteins using high-resolution mass spectrometry. Ann Acad Med Singapore. 2006, 35: 400-407.
Wong YW, Chew J, Yang H, Tan DT, Beuerman R: Expression of insulin-like growth factor binding protein-3 in pterygium tissue. Br J Ophthalmol. 2006, 90: 769-772. 10.1136/bjo.2005.087486.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, et al: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264. 10.1093/biostatistics/4.2.249.
Yang H, Haddad H, Tomas C, Alsaker K, Papoutsakis ET: A segmental nearest neighbor normalization and gene identification method gives superior results for DNA-array analysis. Proc Natl Acad Sci USA. 2003, 100: 1122-1127. 10.1073/pnas.0237337100.
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, et al: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.
Livak KJ, Flood SJ, Marmaro J, Giusti W, Deetz K: Oligonucleotides with fluorescent dyes at opposite ends provide a quenched probe system useful for detecting PCR product and nucleic acid hybridization. PCR Methods Appl. 1995, 4: 357-362.
Lafontaine PO, Arnal M, Buron N, Solary E, Bron AM, et al: Trefoil factor family mRNA and protein expression in pterygium. Int J Oncol. 2005, 27: 997-1003.
Haider AS, Peters SB, Kaporis H, Cardinale I, Fei J, et al: Genomic analysis defines a cancer-specific gene expression signature for human squamous cell carcinoma and distinguishes malignant hyperproliferation from benign hyperplasia. J Invest Dermatol. 2006, 126: 869-881. 10.1038/sj.jid.5700157.
Weiber H, Andersson C, Murne A, Rannevik G, Lindstrom C, et al: Beta microseminoprotein is not a prostate-specific protein. Its identification in mucous glands and secretions. Am J Pathol. 1990, 137: 593-603.
Chang HY, Sneddon JB, Alizadeh AA, Sood R, West RB, et al: Gene expression signature of fibroblast serum response predicts human cancer progression: similarities between tumors and wounds. PLoS Biol. 2004, 2: E7-10.1371/journal.pbio.0020007.
Iyer VR, Eisen MB, Ross DT, Schuler G, Moore T, et al: The transcriptional program in the response of human fibroblasts to serum. Science. 1999, 283: 83-87. 10.1126/science.283.5398.83.
Muro AF, Chauhan AK, Gajovic S, Iaconcig A, Porro F, et al: Regulated splicing of the fibronectin EDA exon is essential for proper skin wound healing and normal lifespan. J Cell Biol. 2003, 162: 149-160. 10.1083/jcb.200212079.
Campo McKnight DA, Sosnoski DM, Koblinski JE, Gay CV: Roles of osteonectin in the migration of breast cancer cells into bone. J Cell Biochem. 2006, 97: 288-302. 10.1002/jcb.20644.
Liu AY, Bradner RC, Vessella RL: Decreased expression of prostatic secretory protein PSP94 in prostate cancer. Cancer Lett. 1993, 74: 91-99. 10.1016/0304-3835(93)90049-F.
Taupin D, Podolsky DK: Trefoil factors: initiators of mucosal healing. Nat Rev Mol Cell Biol. 2003, 4: 721-732. 10.1038/nrm1203.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/2/14/prepub
We would like to thank Tracy Chai for performing the microarray and immunofluorescent staining experiments and Dr Wong Yong Wee for contributing to the analysis in the paper. The study is funded by Singapore Biomedical Research Council (BMRC) 03/1/35/19/231 and National Medical Research Council (NMRC) IBG. BMRC has approved the submission of this manuscript for publication.
The authors declare that they have no competing interests.
LT drafted and revised the manuscript, performed the analysis and coordinated study. JC carried out the molecular studies and microarray experiments, HY carried out the analysis on GSEA and provided help with bioinformatics, LPKA conceived of the study, coordinated study, revised the manuscript. DTHT critically reviewed manuscript, provided patient materials and RB initiated the study and formulated the design, revised the manuscript and approved the submission.
Electronic supplementary material
Additional File 2: Primary and recurrent pterygium specific genes. Heat map showing pterygium specific and recurrence specific genes. Note that in recurrent pterygium a larger number of genes were relatively down-regulated (green) compared to primary pterygium and un-involved conjunctiva, and a smaller but still large number of genes were up-regulated (red) compared to primary pterygium and un-involved conjunctiva. (DOC 82 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Tong, L., Chew, J., Yang, H. et al. Distinct gene subsets in pterygia formation and recurrence: dissecting complex biological phenomenon using genome wide expression data. BMC Med Genomics 2, 14 (2009) doi:10.1186/1755-8794-2-14
- Global Gene Expression
- Cutaneous Squamous Cell Carcinoma
- S100 Calcium Binding Protein
- Conjunctival Tissue