<title>Introduction</title> <p>The retinal pigment epithelium (RPE) is a monolayer of hexagonally packed, highly pigmented, polarized cells located on the posterior wall of the eye, whose apical membranes are intimately associated with the outer segments of photoreceptor cells of the neural retina (NR). The RPE forms the outer blood-retinal barrier and carries out important physiologic and protective functions necessary for visual processing in rod and cone cells, such as retinoid metabolism, phagocytosis of discarded rod and cone outer segments, absorption of stray light to preserve visual acuity, control of water and ion flow between the neural retina and choroid, and protection of the neural retina from oxidative damage [##REF##15221858##1##–##REF##15987797##3##]. These functions require maintenance of intimate association of the RPE with the NR, which if disrupted leads to severe ocular pathologies. In situ, the RPE is both non-proliferative and non-migratory, yet these cells have been shown to exhibit a high degree of plasticity in vitro. The plasticity in function and phenotype of the RPE can be recapitulated in vivo when damage occurs to the retina in the form of a retinal detachment or tear. Therefore retinal detachments and tears require surgical repair of the RPE/NR interaction [##REF##15814473##4##]. However, as many as 10% of all rhegmatogenous retinal detachments, in which a retinal tear occurs, result in failure due to the occurrence of proliferative vitreoretinopathy (PVR).</p> <p>PVR is characterized by the formation of epiretinal membranes, comprised in part of dedifferentiated RPE cells that undergo epithelial-mesenchymal transformation (EMT) and contribute to this fibroplastic response [##UREF##0##2##,##REF##7525361##5##,##REF##9716190##6##]. While the cause of PVR remains unknown, one aspect of this disease includes changes in the expression of a variety of genes that regulate RPE cell phenotype. Identification of the transcription factors that maintain RPE cells in a differentiated non-proliferative and non-migratory state could provide potential therapeutic targets. Studies of EMT in a variety of cell types in development and disease have begun to identity such factors [##REF##14644200##7##, ####REF##12631389##8##, ##REF##15713751##9####15713751##9##]. More recently, high-throughput technology such as microarray analysis has identified changes in gene expression in RPE cells undergoing EMT in vitro, including genes involved in DNA synthesis and repair, cell cycle, intracellular signaling, and cell adhesion [##REF##12483023##10##,##REF##11820701##11##]. However, the mechanisms by which these many changes are controlled during EMT, including the transcriptional regulators that may coordinate this process, remain to be elucidated.</p> <p>The purpose of the present study was to identify transcription factors that regulate EMT in RPE cells. Genes were identified that are upregulated or down-regulated during RPE EMT, and the genomics tool Promoter Analysis and Interaction Network Toolset (PAINT v 3.3) [##REF##14583114##12##] was then used to generate models of the promoter regions of these genes including predictions of those TFs that regulate their expression. We then tested the validity of these models using RT–PCR to analyze expression of the TFs in differentiated and undifferentiated RPE cells and indeed identified several TFs that are differentially expressed between these two RPE cell states. The results of these studies indicate that a combination of computational and biochemical approaches can be successfully applied to analyze these complex events.</p>
<title>Methods</title> <title>Computational analysis of gene regulatory regions</title> <p>Identification and analysis of the regulatory regions of genes expressed in RPE cells was performed using the Promoter Analysis and Interaction Network Toolset version 3.3 (PAINT) program [##REF##14583114##12##]. The target gene set analyzed comprised sixty genes differentially expressed during EMT of RPE cells, selected from previous literature reports and additional genes under study in our laboratory [##REF##12483023##10##,##REF##11820701##11##] (##TAB##0##Table 1##). This set was divided into two subsets representing those preferentially expressed in the undifferentiated or differentiated state, and further refined by functional gene ontology [##REF##10802651##13##]. Homologous genes in humans, mice, rats, and chickens were identified using the Ensembl database and the Ensembl gene ID for each was used in the PAINT program [##REF##15608235##14##,##REF##11752248##15##]. For each gene, PAINT identified the putative transcription start site (TSS) and subsequent promoter analysis to identify transcription response elements (TREs) was performed to 5,000 bp upstream of the TSS, with exclusion of complimentary strand analysis and 1.0 core similarity threshold [##REF##14583114##12##]. In addition, sequence data corresponding to the first intron of each gene was retrieved from Ensembl and also entered into PAINT, as FASTA formatted sequences, for further analysis. TREs within both the promoters as well as first introns were identified using the MATCH/TRANSFAC database [##REF##12520026##16##]. These TREs were entered into the Feasibility Network Builder module of PAINT (FeasNet Builder), which constructed a candidate interaction matrix (CIM), a graphic representation of the occurrence of these TREs within the gene set. Enrichment analysis was performed using PAINT to compute the Fisher's exact test p-values indicating relative over-representation of TREs within the selected gene set as compared to the larger background gene set, the 588 genes present on the Human Atlas Array (Cat. No. 7740–1; Clontech, Palo Alto, CA) used to originally identify the set of sixty differentially regulated RPE genes [##REF##12483023##10##,##REF##11820701##11##]. In each of the analyses, multiple testing correction was applied using a false discovery rate (FDR) estimate [##UREF##1##17##]. In all cases, the multiple testing corrections were not informative as they did not result in any over-represented TREs, since the FDR was consistently above 90% for all TREs. Therefore, we followed a discovery approach and chose a threshold of p&lt;0.1 on the Fisher's exact test p-value to identify those TREs to be included in further filtering as described below. Models of RPE gene interaction networks, based upon results of the FeaseNet analysis, were graphically generated using GraphViz [##UREF##2##18##].</p> <title>Generation of models for gene regulatory regions and selection of targets for biochemical analysis</title> <p>Global regulatory models for gene sets coordinately expressed in RPE cells, as well as for individual genes differentially expressed in differentiated or undifferentiated RPE cells, were developed based upon PAINT-derived computational data and were constructed by comparing phylogenetically conserved transcriptional regulation across the four species: human, mouse, rat, and chicken. To establish criteria for selection of specific transcription factors for further analysis, we assigned values to each TRE based upon their frequency of detection across the coordinately expressed gene sets and their evolutionary conservation. An evolutionary conservation factor (ECF) of 1 to 3 points was assigned to each TF, where 1 indicates presence on human genes, 2 indicates human and mouse or rat genes, and 3 indicates a presence on human and chicken as well as either mouse or rat genes. In addition, TREs were scored according to a frequency ratio (FR) derived from the ratio of the percent occurrence of a given TRE in a specific gene subset divided by its frequency of occurrence in the background gene set, for the human genome data. Ultimate TRE selection was based upon criteria filters of a combined ECF score of 2 or greater along with a FR greater than 3.</p> <title>Cell culture</title> <p>Primary cultures of chick embryo RPE cells (cRPE) were established from RPE tissues obtained from fertile White Leghorn chicken eggs maintained in a humidified atmosphere at 37 °C until embryonic day 10, corresponding to Hamburger and Hamilton stage 36 [##REF##1304821##19##,##REF##1628748##20##]. Eyes were enucleated, the anterior segment and vitreous were removed, and the posterior eyecup divided in half and incubated in HBSG (HEPES buffered saline with glucose containing 1 mg/ml glucose, 10 mM HEPES, pH 7.4, 3 mM KCl, and 0.15 mM NaCl). The neural retina was removed and the eyecup was incubated in HBSG containing 20 mM EDTA (ethylene diamine tetraacetic acid) for 1 h. The eye cups were then rinsed for an additional 30 min in HBSG and the RPE was dissected from the choroid. RPE tissue was collected by gentle centrifugation and resuspended in MEM (Eagle's minimum essential medium; Sigma-Aldrich, St. Louis, MO) supplemented with 10% FBS (Fetal Bovine Serum; Sigma-Aldrich), 0.22% sodium bicarbonate and 1% penicillin, streptomycin and amphotericin (Gibco/Invitrogen, Carlsbad, CA). The RPE tissue was mechanically dissociated into single cells and plated onto plastic six well tissue culture plates at a density of 1 eye/well. The plates were previously coated with 1.1 μg/cm<sup>2</sup> mouse laminin (Sigma-Aldrich) 12 h before cell plating. The cultures were maintained in a humidified atmosphere of 5% CO<sub>2</sub> at 37 °C for 48 h at which time all non-adherent cells were removed by rinsing with fresh medium [##REF##3148924##21##].</p> <p>The human RPE-derived cell line ARPE-19 [##REF##8698076##22##] was maintained cultured in DMEM-F12 (Catalog number D8900; Dulbecco's modified eagles minimum nutrient mixture F12 Ham; Sigma-Aldrich) supplemented with 5% FBS, 2 mM L-glutamine, 0.348% sodium bicarbonate and 1% Antibiotic-Antimycotic in T25 culture flasks in a humidified atmosphere of 5% CO<sub>2</sub>. To maintain cells in an undifferentiated state, they were passaged before obtaining confluence. To obtain differentiated cells, cells were grown to confluence and then maintained in DMEM-F12 as above except that the serum was reduced to 1% [##REF##12657613##23##]. These cultures reach confluence 2–3 weeks after passaging and differentiate within 4–6 weeks, though the cultures can be kept in a differentiated state for extended culture periods. After 4–6 weeks in culture the cells exhibit hexagonal packing of pigmented, polarized epithelia and expression of CRALBP and RPE65 typical of morphological and biochemical markers, respectively, of RPE cells in vivo [##REF##8698076##22##].</p> <title>RNA extraction and reverse transcriptase polymerase chain reaction amplification</title> <p>Total RNA was extracted from all cell types using the Micro-to-Midi Total RNA Purification System (Invitrogen, Carlsbad, CA) following manufacturers specifications for isolation of RNA from animal tissues. Purity and concentration of RNA from each sample was assessed on a spectrophotometer at 260 nm wavelength (A<sub>260</sub>). All primers used in this study are listed in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-1.pdf">Appendix 1</ext-link> and were designed using the web-based tool <ext-link ext-link-type="uri" xlink:href="http://www.genefisher.de">GeneFisher</ext-link> Primer Design program, and all ranged between 20 and 22 nucleotides in length, with melting temperatures of 50–65 °C and G-C content between 40%–60% [##REF##8877506##24##]. Primer specificity was determined by the nucleotide-nucleotide basic local alignment search tool (BLASTn) set to the specific species genome, with acceptance of a primer pair based upon expect-value of less than 1 (e-value &lt;1) [##REF##9254694##25##]. RT–PCR was performed using the SuperScript III One-Step RT–PCR System with Platinum Taq (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. The reaction parameters were (a) cDNA synthesis, 1 cycle at 55 °C for 30 min, (b) denaturation, 1 cycle at 94 °C for 2 min, (c) amplification, 40 cycles at 94 °C for 15 s, (d) melting at primer-specific temperatures for 30 s, (e) extension, 68 °C for 1 min, and (f) final extension, 1 cycle at 68 °C for 5 min. RT–PCR products were analyzed by agarose gel electrophoresis on 1.25% gels containing ethidium bromide (1 μg/ml). The resultant bands were visualized and recorded under ultraviolet light using the Kodak 1D photo system.</p>
<title>Results</title> <title>Transcriptional regulatory network analysis of differentially expressed retinal pigment epithelium genes reveals gene set-specific clusters of transcription response elements</title> <p>To identify target genes for the analysis of transcription response elements by PAINT, a gene set containing two clusters representing those genes whose expression is preferentially associated with the undifferentiated versus differentiated state of the ARPE-19 cells was established (##TAB##0##Table 1##). Thus each gene cluster would be predicted to be preferentially associated with one or more TREs involved in their respective coordinate regulation of expression. The genes within the differentiated and undifferentiated clusters were separately analyzed using PAINT to identify and statistically analyze the occurrence of TREs within the promoter region, including 5,000 bp of sequence upstream from the transcription start site, as well as the entire first intron sequence. The complete results of this analysis by the FeasNet Builder module of PAINT are listed in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-2.pdf">Appendix 2</ext-link> for the promoter regions and <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-3.pdf">Appendix 3</ext-link> for the first introns. The analysis identified those TREs that are statistically over-represented on the promoters and first introns of genes within either of the two gene clusters, as compared to their occurrence with the larger background set of promoter and intron regions. TREs with a p-value of less than 0.1 were deemed significant and included for further analysis. These data were then converted into graphic representations of the TRE network or CIM, using the FeasNet Viewer module of PAINT, as illustrated in ##FIG##0##Figure 1## for the promoter regions and in ##FIG##1##Figure 2## for the first introns. Inspection of the matrix patterns for the respective differentiated and undifferentiated gene clusters indicated that distinct subsets of TREs were associated with each of these two clusters, suggesting that these TREs could represent components of the networks regulating coordinated expression of genes within each cluster. ##FIG##0##Figure 1## and ##FIG##1##Figure 2## represent a subset of the complete CIM, and include only those TREs with a p-value of less than 0.1. The complete CIMs representing the results of the analysis of the full 60 gene set are tabulated in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-2.pdf">Appendix 2</ext-link> and <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-3.pdf">Appendix 3</ext-link>, as well as illustrated in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-4.pdf">Appendix 4</ext-link> and <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-5.pdf">Appendix 5</ext-link>. The results of the CIM analysis indicate that several specific TREs are statistically overrepresented within the regulatory regions of genes within each gene subset, and furthermore that several of these TREs differ between the two subsets. These distinctly represented TREs are thus candidates for further analysis regarding their potential regulatory role in coordinating gene expression during RPE cell differentiation.</p> <p>To further visualize the potential interrelationships of the candidate TREs and genes within each cluster, the data generated by the FeasNet Viewer module was graphically presented using GraphViz to generate a regulatory network diagram for the human gene set, as illustrated in ##FIG##2##Figure 3## for the promoter regions and ##FIG##3##Figure 4## for the first introns. These results identify nodes within each network and further identify those TREs with the potential to coordinately regulate subsets of genes within each cluster. The TREs that comprise the nodes of this visualization represent all those TREs which are overrepresented among the genes of each cluster, and were found to be associated with as few as one gene within each cluster to as many as 28 genes within a cluster. Those TREs exhibiting high levels of interconnectivity are among the best candidates for coordinate regulation of RPE genes.</p> <p>To identify those TREs of particular significance in the differential regulation of RPE genes during EMT, we next performed an analysis comparing all the overrepresented TREs from the CIM analysis with respect to their relative frequency of occurrence between the differentiated, undifferentiated and background gene clusters. While the previous analyses identified TREs overrepresented in either or both of the two gene clusters, the frequency analysis further distinguished those TREs preferentially associated with either the differentiated or undifferentiated gene cluster. The results, which are shown in ##FIG##4##Figure 5## for the promoter regions and ##FIG##5##Figure 6## for the intron regions, indicate that a select subgroup of the TREs can be assigned as potential regulators of one of each of the two cell states. For example, in the promoter regions, Hand1:E47 and COMP1 show increased relative frequency among down-regulated genes, whereas Oct-1 and SREBP-1 do so among upregulated genes. In the first intron regions, CDP and IRF show increased relative frequency among down-regulated genes, whereas GATA-3 and MAZR do so among upregulated genes.</p> <title>Phylogenetic comparison of PAINT analyses identifies evolutionarily conserved transcriptional regulatory elements across divergent species</title> <p>The above studies of gene regulatory networks regulating RPE cell differentiation were focused on human genome sequence analysis. However, evolutionary conservation analysis of regulatory regions can also be useful for identifying functionally important sites. Therefore, to further identify those elements of these hypothetical networks that would most likely be of functional significance due to their evolutionary conservation, we performed parallel PAINT analyses similar to those described above for human genomic elements for three additional model species: mouse (<italic>Mus musculus</italic>), rat (<italic>Rattus norvegicus</italic>), and chicken (<italic>Gallus gallus</italic>). The results of the FeasNet Builder analysis in the three additional species gene sets indicate that, as observed for the analysis of human genes, several overrepresented TREs are again identified (<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-6.pdf">Appendix 6</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-8.pdf">Appendix 8</ext-link>, and <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-10.pdf">Appendix 10</ext-link> for the promoter regions and in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-7.pdf">Appendix 7</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-9.pdf">Appendix 9</ext-link>, and <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-11.pdf">Appendix 11</ext-link> for the first introns) for each species. The graphical candidate interaction matrix for each of these three species is shown in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-12.pdf">Appendix 12</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-13.pdf">Appendix 13</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-14.pdf">Appendix 14</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-15.pdf">Appendix 15</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-16.pdf">Appendix 16</ext-link>, and <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-17.pdf">Appendix 17</ext-link>, while the GraphViz output of the FeasNet Viewer analyses illustrating the regulatory network diagrams for these three species is shown in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-18.pdf">Appendix 18</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-19.pdf">Appendix 19</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-20.pdf">Appendix 20</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-21.pdf">Appendix 21</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-22.pdf">Appendix 22</ext-link>, and <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-23.pdf">Appendix 23</ext-link>, which both again indicate the presence of a complex regulatory network within each of these species. As for the human genome, we also compared the frequency of occurrence of each identified TRE between each gene cluster and the background gene set for each of the three model species. These frequency analyses are shown in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-24.pdf">Appendix 24</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-25.pdf">Appendix 25</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-26.pdf">Appendix 26</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-27.pdf">Appendix 27</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-28.pdf">Appendix 28</ext-link>, and <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-29.pdf">Appendix 29</ext-link>. As for the human genome analysis discussed above, these results of the CIM, GraphViz and frequency analyses identify TREs for each of the three species which are also candidates for regulation of gene expression during RPE differentiation. Comparison of the results from the human as well as model system analyses identified both phylogenetically conserved as well as species-specific TREs, which were incorporated into criteria tables and gene regulatory region models as described below (##TAB##1##Table 2## and ##FIG##6##Figure 7##).</p> <title>Compilation of human and cross-species transcriptional regulatory network analyses generates global as well as gene-specific regulatory models</title> <p>To generate a global working model of gene expression regulation in the RPE, we compared the above results across the four species analyzed and identified those TREs of highest evolutionary conservation, as well as highest frequency of occurrence, by calculation of an evolutionary conservation factor (ECF) score, as well as a frequency ratio (FR) score, for each TRE identified in the computational analysis, as described in Methods. The results of these calculations are shown in ##TAB##1##Table 2##, with the table containing only those results for TREs that passed the three criteria of a p-value &gt;0.1, an ECF greater or equal to 2, and an FR greater or equal to 3. The full data set of calculated ECFs and FRs for all the TREs identified in the computational analysis is shown in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-30.pdf">Appendix 30</ext-link>, for the promoter regions and <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-31.pdf">Appendix 31</ext-link> for the first intron regions. These calculations provide a scheme for ranking TREs, and by extension their corresponding TFs, for inclusion into models of RPE gene regulation, as well as for further studies probing their expression and function in RPE differentiation.</p> <p>The TREs passing all these criteria, and hence representing potentially evolutionarily conserved regulators of RPE gene expression, were incorporated into archetypal gene models for the coordinated up- and down-regulation of genes during RPE cell differentiation, as shown in ##FIG##6##Figure 7##. Although evolutionary conservation is a potential indicator of conserved function, in recognition of the likelihood that species-specific aspects to coordinated gene regulation also exist, we also developed an additional model with modified ECF criteria focused on regulation in the human genome, such that incorporation of TREs into this second model was dependent on the TRE occurring within the human gene set and any one additional species, while maintaining the same p-value and FR criteria filters. The resulting models are shown in ##FIG##7##Figure 8##, and share many elements of the trans-species model, in some cases eliminating TREs not found in the human analysis, while adding some others due to the less stringent ECF criteria. Finally, additional gene-specific regulation models were constructed for selected genes representing paired members of multigene families, or other markers, with reciprocal expression during RPE cell differentiation. These paired models were constructed for the cell adhesion proteins N- and R-cadherin, for the lactate transporters MCT3 and −4, and for α-smooth muscle actin and RPE65 (##FIG##8##Figure 9##). These gene regulation models identify those TREs most likely to coordinate expression of specific genes as well as broader sets of up- and down-regulated genes during RPE cell differentiation, and provided targets for validation of these models through direct biochemical analysis as described below [##REF##8698076##22##,##REF##12657613##23##,##UREF##3##26##].</p> <title>Validation of gene models by identification of transcription factors whose expression is dependent upon the state of retinal pigment epithelium cell differentiation</title> <p>The above gene regulatory models identified TREs that could serve as elements of the gene regulatory network during RPE cell differentiation. For these TREs to play a role in regulation of their associated genes, the activity of the corresponding transcription factors would be expected to be dynamically regulated at appropriate times to effect such control. One common mechanism of regulation for TF activity is at the transcriptional level resulting in differential steady-state mRNA expression levels. Thus, to test the computationally derived models, we used RT–PCR to assay for the presence of mRNA encoding transcription factors (TFs) predicted to play a regulatory role in the RPE. The results obtained from comparison of mRNA extracted from undifferentiated and differentiated RPE cells, using both the human ARPE-19 cell line and primary embryonic chicken RPE cells, are shown in ##FIG##9##Figure 10##, ##FIG##10##Figure 11##, ##FIG##11##Figure 12##. The analyses of ARPE-19 cells were performed using both undifferentiated cells that been in culture for one month, exhibiting morphology of fusiform, unpigmented mesenchymal cells, as well as differentiated ARPE-19 cells exhibiting morphology of polygonal, darkly pigmented, epithelial cells (##FIG##9##Figure 10A,B##). RT–PCR was first performed on each respective cell population targeting mRNAs for α-SMA and RPE65, respective markers of the mesenchymal undifferentiated and epithelial differentiated state of RPE cells [##REF##8698076##22##]. These results indeed demonstrated that mRNA encoding α-SMA, but not RPE65, was expressed in undifferentiated cells, whereas RPE65 mRNA was readily detected among differentiated cells with a reduced level of α-SMA mRNA (##FIG##9##Figure 10C,D##). Semi-quantitative RT–PCR was used to further distinguish these levels of α-SMA mRNA, which more clearly distinguished these two cell states (##FIG##9##Figure 10E##). Having verified that these two markers were distinctly expressed between these two cell populations, we then used RT–PCR to determine the levels of mRNA encoding the specific transcription factors previously identified by the computational analysis, which were expected to fall into three categories, exhibiting either quantitative, qualitative, or no differences between the two test cell populations. Of these TFs, mRNAs encoding four were found to be reciprocally expressed in differentiated versus undifferentiated ARPE-19 cells, with Oct-1 and TFE3 detected only in differentiated cells, and Core Binding Factor and NKX3A detected only in undifferentiated cells (##FIG##10##Figure 11##). mRNAs encoding additional transcription factors, including GATA-1, IRF-1, and SMAD3, were detected in both cell states (##FIG##10##Figure 11A##). Semi-quantitative RT–PCR was then used to further analyze differences in the expression patterns of these factors, with GATA-1 detected at higher levels in differentiated cells, whereas IRF-1 and SMAD3 were detected at higher levels in undifferentiated cells (##FIG##10##Figure 11B##).</p> <p>Similar analyses were performed using freshly isolated chicken RPE tissues and primary chick RPE cultures. When cultured, chick RPE cells re-enter the cell cycle and de-differentiate [##REF##9066129##27##]. RNA was prepared from freshly isolated RPE cells as well as from cells cultured for five days in vitro, and both cell populations were probed for TFs corresponding to those TREs identified in association with both the global as well as gene-specific regulatory models. The results indicated that, similar to ARPE19 cells, and consistent with model predictions, the primary chick RPE cells also revealed reciprocal TF expression between differentiated and undifferentiated cells. mRNAs encoding AML-1 and HNF-1 were detected only in the differentiated chicken RPE cells, whereas mRNAs encoding HNF-3 and SREBP-1 were detected only in the undifferentiated cells, while two additional mRNAs encoding TFs DP-1 and TFII-I were detected in both cell populations (##FIG##11##Figure 12##). Finally, as an adjunct to the PAINT-derived analyses, an additional series of RT–PCR reactions were performed to determine whether other TFs, not identified through PAINT, but known from prior studies to be involved in EMT of epithelial cells other than RPE, were expressed in chick or ARPE-19 cells. As shown in ##FIG##12##Figure 13##, RT–PCR amplification of mRNA from undifferentiated and differentiated primary chick RPE cells generated similar levels of amplicons for Slug, Twist and SIP1, whereas Snail was detected at higher levels in undifferentiated cells, and LEF1 was detected only in differentiated cells. When similar analyses were performed with total RNA isolated from differentiated or undifferentiated ARPE-19 cells, Slug, Snail, Twist, and SIP1 were not detected, while SMAD2 was detected at equal levels in both samples. Interestingly, while LEF-1 was also detected in both cell populations, a distinct additional amplicon was detected in differentiated RPE cells, indicating that differential splicing of this gene transcript occurs during the course of RPE differentiation. Our inability to detect expression of certain classical mediators of EMT such as Snail, Slug, Twist or SIP1 in ARPE-19 cells, while we were able to detect them in primary cultures of embryonic chick RPE cells, may be related to the different stages of development represented by these two model systems (embryonic versus adult), to the unique properties of RPE cells compared to other epithelial cell types that may exhibit species-specific differences, or to some specific phenotypic property of ARPE-19 cells that arose during their derivation [##REF##16801545##28##,##REF##17000672##29##]. Overall, the results of the RT–PCR analyses indicate that the computational biology approach was successful at identifying transcription factors whose expression is regulated during RPE cell differentiation, and thus may play a role in control of differential gene expression and modulation of RPE cell phenotype.</p>
<title>Discussion</title> <p>The results of the present study have permitted the construction of several hypothetical models for regulation of genes in RPE cells during EMT, each generated using a different set of theoretical boundaries and statistical criteria. While the computational approach using the PAINT toolset has been previously applied to other cell types [##REF##14583114##12##,##REF##15252187##30##], to the best of our knowledge the present work represents the first application to the analysis of RPE cell differentiation. A strength of the models developed here is that they make strong predictions of which TFs would be expected to be differentially acting during phenotypic changes in RPE cells, predictions which were successfully tested and positively borne out by the RT–PCR analyses in the present studies. These results form the basis for design of future studies that will be directed at testing the function of these various TFs in regulating RPE cell phenotype. These experiments are guided by the integration of the experimental results into a comprehensive model for RPE gene regulation (##FIG##13##Figure 14##), which indicates for each TRE included in the final model, the various criteria filters that led to its inclusion, including evolutionary conservation, frequency of occurrence, position in a gene regulatory network node, and generation of a positive amplicon in RT–PCR validation assays. Two TFs, Oct1 and HNF1, although not previously identified with respect to EMT, pass all of these four criteria, may play a unique role in this context in RPE cells, and thus are identified as excellent candidates for direct functional analysis in future studies.</p> <p>While not previously analyzed in the context of RPE cells, the transcription factors identified in the present study by the PAINT and RT–PCR analyses can be categorized with respect to other cell types into three groups that include (1) TFs not previously associated with EMT; (2) TFs that, while not previously known to directly affect EMT, have been shown to regulate cellular processes that are components of EMT, and (3) TFs previously shown to directly affect EMT principally in other cell types. The first group includes the TF GATA-1, which of all the TFs identified in this study is the only one not directly linked to an EMT-related process. This factor is expressed in cells of the erythroid lineage and is essential for proper erythroid development, but its potential role in regulation of epithelial cell phenotype remains to be determined [##REF##15252187##30##, ####REF##8321207##31##, ##REF##8206977##32####8206977##32##]. The second group encompasses the factors such as Oct-1, HNF-1, NKX3A, IRF-1, SREBP-1, and Core Binding Factor, which have not been specifically linked to EMT, yet regulate processes such as cell migration, cell adhesion and metabolic pathways associated with EMT. Oct-1 and HNF-1 act as important regulators of development processes such as neural tube development [##REF##1365913##33##, ####REF##9232592##34##, ##REF##12860991##35##, ##REF##2407588##36####2407588##36##]. NKX3A, a homolog of NKX2–5 that functions to activate N-cadherin expression in cardiac development, may function in a similar manner by activating N-cadherin expression, which has been shown to be highly expressed in metastatic cancer cells [##REF##9537602##37##, ####REF##14645532##38##, ##REF##15925411##39##, ##REF##11137288##40##, ##REF##9226374##41####9226374##41##] and is upregulated during RPE de-differentiation. Interferon regulatory factor-1 (IRF-1) is another factor that falls into this category, in that it plays an important tumor suppressive role in a wide variety of human neoplasias [##REF##10524230##42##,##REF##11846974##43##]. Sterol regulatory element binding protein-1 (SREBP-1) is known to affect expression of lipogenic genes in the liver, which is of interest insofar as cells undergoing EMTs possess altered fatty acid and glucose/insulin metabolism [##REF##11343966##44##]. Previous work has shown a switch to aerobic glycolysis when cells begin to migrate in the initial stages of EMT, an action that may be mediated by SREBP-1 [##REF##15073099##45##]. Core binding factor (CBF) may also be indirectly involved in EMT, in that it interacts with members of the TGF-β signaling factor to influence cell growth and differentiation [##REF##15864279##46##]. Finally, SMAD3 and TFE3 constitute the last group and have previously been strongly implicated in the signaling pathways associated with TGF-β induced EMT, whereby they activate LEF-1 transcription, a major EMT inducer [##REF##15118077##47##, ####REF##10973944##48##, ##REF##10707962##49##, ##REF##10082515##50##, ##REF##16253118##51####16253118##51##]. These two TFs are thus also identified as excellent candidates for further analysis in RPE cells, since they were identified through the PAINT analysis and have also been previously identified as regulators of genes associated with EMT in several cell types [##REF##15689496##52##, ####REF##14617750##53##, ##REF##9765209##54##, ##REF##16123809##55##, ##REF##16737956##56####16737956##56##], and in one report in RPE cells [##REF##15273699##57##]. Thus, while this study has identified several novel potential regulators of the RPE, the concordance between certain results of the present study and prior reports provides further validation of a combined in silico computational approach as an adjunct to in vivo as well as in vitro biochemical and cell biologic studies.</p> <p>One potential limitation of the present approach is indicated by the apparent lack of identification by the PAINT analysis of some TFs that have been previously associated with EMT in other cell types and that may play a role in RPE as well. These include additional downstream mediators of TGF-β signaling pathways besides SMAD3 such as Snail, Slug, Twist, SMAD2, SIP1, β-catenin, and LEF-1 [##REF##16123809##55##,##REF##15273699##57##, ####REF##16010672##58##, ##REF##16568083##59##, ##REF##15294153##60##, ##REF##16232121##61####16232121##61##]. As one approach to addressing this, we performed RT–PCR assays to determine the presence of mRNAs corresponding to these TFs, and did detect several of these in RPE cells, although only LEF-1 was indicated to be differentially expressed between differentiated and undifferentiated cells. Interestingly, differential splicing of LEF-1 mRNA, as detected here, has been reported previously to generate several isoforms that may confer distinct functionality on this TF [##REF##10756202##62##]. The primers used in our RT–PCR flank the third through sixth exons, the third of which encodes a premature stop codon that can generate a form of LEF-1 retaining its β-catenin binding site but lacking DNA binding domain and nuclear localization signal [##REF##10756202##62##]. Given the key role of LEF-1 in the TGFβ signaling pathway, this may indicate one possible means through which modulation of such pathways occurs during RPE differentiation. For any TRE to be analyzed via PAINT, its sequences must be available in TRANSFAC database, and our manual inspection of this database revealed that no sequences are available for TREs corresponding to Slug in any species, and Snail, Twist and SIP1 sequences are available only for the mouse genome, whereas only SMAD-2 and LEF-1 sequences are available for all species analyzed in this study. Of these TREs for which at least partial sequence data was available, although some were indeed detected by PAINT in some genes within the clusters, only LEF-1 was enriched with a p-value &lt;0.1, but had a low ECF value.</p> <p>A second limitation of the present study is the inherent variability observed in the phenotype of ARPE-19 cells. Several reports have indicated that there is a degree of variability, depending on the culture conditions such as serum concentration and growth substrate, as well as differences between ARPE-19 cells and native human RPE [##REF##14982971##63##,##REF##16234463##64##]. While we acknowledge that this inherent variability exists, the ARPE-19 cells used in the present study were from undifferentiated and well differentiated cell cultures, respectively, as defined by both morphological as well as biochemical criteria.</p> <p>In situ, the RPE is a monolayer of morphologically and functionally polarized non-proliferative and non-migratory cells whose unique properties are essential to the proper development and function of the retina. However, these cells are known to exhibit a high degree of plasticity in phenotype and function both in vitro and in vivo [##REF##15221858##1##,##REF##15987797##3##]. Delineating the mechanisms underlying this plasticity is essential to understanding the conditions under which RPE cells undergo these changes, and is critical to developing preventive and therapeutic interventions for conditions in which RPE plasticity may lead to retinal diseases such as proliferative vitreoretinopathy (PVR) [##REF##9716190##6##]. Current therapeutic techniques used to treat retinal detachments and their complications are limited to invasive surgical procedures aimed at physically re-attaching the sensory portion of the retina to the underlying RPE, and removing epiretinal membranes, such as laser- or cryo-therapy, supplemented by pneumatic retinopexy, scleral buckling or vitrectomy. Presently, PVR occurs as a complication in up to 10% of surgical retinal detachment repairs, making it the most common post-surgical complication associated with these procedures [##REF##15814473##4##]. Development of non-surgical or adjunct treatments for PVR will require a better understanding of the underlying biology of the genetic and epigenetic mechanisms regulating RPE cell phenotype and underlying the plasticity exhibited by RPE cells. Since this plasticity likely reflects changes in the expression of a wide variety of gene products, and thus ultimately the coordinated function of several transcription factors, the present study was designed to apply the tools of computational biology to identify transcription factors whose function could modulate changes in RPE cell phenotype. The TFs identified in this study thus become excellent candidates for further analysis of their role in this process.</p> <p>In conclusion, we have predicted and experimentally verified the differential expression of several transcription factors including Oct-1, HNF-1, SMAD3, TFE, Core binding factor, GATA-1, IRF, NKX3A, SREBP-1, and LEF-1 that may be of importance in the regulation of genes during EMT of RPE cells, as determined first by computational analysis and modeling, and then tested by direct RT–PCR analysis. The results clearly indicate that several of these TFs are differentially regulated during RPE differentiation and thus may play a role in epithelial-mesenchymal transformations of RPE cells in both developmental and disease processes. These TFs are thus excellent targets for further studies directed at testing their role as regulators of RPE cell phenotype, and consequently may also provide future targets for therapeutic interventions in cases of PVR and other related disorders of the eye.</p>
<p>This is an open-access article distributed under the terms of the
        Creative Commons Attribution License, which permits unrestricted use,
        distribution, and reproduction in any medium, provided the original
        work is properly cited.</p> <title>Purpose</title> <p>Phenotypic transformation of retinal pigment epithelial (RPE) cells contributes to the onset and progression of ocular proliferative disorders such as proliferative vitreoretinopathy (PVR). The formation of epiretinal membranes in PVR may involve an epithelial-mesenchymal transformation (EMT) of RPE cells as part of an aberrant wound healing response. While the underlying mechanism remains unclear, this likely involves changes in RPE cell gene expression under the control of specific transcription factors (TFs). Thus, the purpose of the present study was to identify TFs that may play a role in this process.</p> <title>Methods</title> <p>Regulatory regions of genes that are differentially regulated during phenotypic transformation of ARPE-19 cells, a human RPE cell line, were subjected to computational analysis using the promoter analysis and interaction network toolset (PAINT). The PAINT analysis was used to identify transcription response elements (TREs) statistically overrepresented in the promoter and first intron regions of two reciprocally regulated RPE gene clusters, across four species including the human genome. These TREs were then used to construct transcriptional regulatory network models of the two RPE gene clusters. The validity of these models was then tested using RT–PCR to detect differential expression of the corresponding TF mRNAs during RPE differentiation in both undifferentiated and differentiated ARPE-19 and primary chicken RPE cell cultures.</p> <title>Results</title> <p>The computational analysis resulted in the successful identification of specific transcription response elements (TREs) and their cognate TFs that are candidates for serving as nodes in a transcriptional regulatory network regulating EMT in RPE cells. The models predicted TFs whose differential expression during RPE EMT was successfully verified by reverse transcriptase polymerase chain reaction (RT–PCR) analysis, including Oct-1, hepatocyte nuclear factor 1 (HNF-1), similar to mothers against decapentaplegic 3 (SMAD3), transcription factor E (TFE), core binding factor, erythroid transcription factor-1 (GATA-1), interferon regulatory factor-1 (IRF), natural killer homeobox 3A (NKX3A), Sterol regulatory element binding protein-1 (SREBP-1), and lymphocyte enhancer factor-1 (LEF-1).</p> <title>Conclusions</title> <p>These studies successfully applied computational modeling and biochemical verification to identify biologically relevant transcription factors that are likely to regulate RPE cell phenotype and pathological changes in RPE in response to diseases or trauma. These TFs may provide potential therapeutic targets for the prevention and treatment of ocular proliferative disorders such as PVR.</p>
<title>Appendix 1. Primer sequences used for ARPE-19 and Chicken Primary RPE cells RT-PCR amplification of selected transcription factor mRNAs.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-1.pdf">Appendix 1</ext-link>.” This will initiate the download of a pdf archive that contains the file. Oligonucleotide primers for RT-PCR amplification were designed and synthesized as described in Methods. Sequences are indicated for both forward and reverse primers used with either human ARPE-19 or chick primary RPE cells. The size of the predicted amplicon is indicated in base pairs.</p> <title>Appendix 2: Statistically over-represented TREs derived from the PAINT analysis of the human promoter regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-2.pdf">Appendix 2</ext-link>.” This will initiate the download of a pdf archive that contains the file. Represented for each TRE identified in the promoter sequences is the cluster in which the TRE was determined to be enriched, the TRANSFAC identifier and the p-value of the TRE in both gene sets.</p> <title>Appendix 3: Statistically over-represented TREs derived from the PAINT analysis of the human first intron regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-3.pdf">Appendix 3</ext-link>.” This will initiate the download of a pdf archive that contains the file. Represented for each TRE identified in the first intron sequences is the cluster in which the TRE was determined to be enriched, the TRANSFAC identifier and the p-value of the TRE in both gene sets.</p> <title>Appendix 4. Complete Candidate Interaction Matrix (CIM) for transcription response elements from PAINT analysis of human gene promoters.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-4.pdf">Appendix 4</ext-link>.” This will initiate the download of a pdf archive that contains the file. The RPE gene set was analyzed by PAINT and a graphic CIM was generated as described in Methods. The y-axis lists the Ensembl Gene identifiers for each gene and the x-axis lists the TRANSFAC TRE identifiers for each TRE found at least once in the promoter region. Genes listed along the y-axis are divided into two clusters that are either up-regulated (blue) or down-regulated (green) during EMT of RPE cells. TREs listed along the x-axis are clustered according to related occurrence pattern calculated using Jaccard's coefficient. The elements within the matrix are color-coded based upon the p-value of each TRE found in the regulatory regions of the genes. A red dot represents a TRE that is statistically significant and therefore over-represented in our gene set, while a blue dot signifies an under-represented TRE and a grey dot stands for a TRE with no statistical significance in our gene list. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 5. Complete Candidate Interaction Matrix (CIM) for transcription response elements from PAINT analysis of human gene first introns.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-5.pdf">Appendix 5</ext-link>.” This will initiate the download of a pdf archive that contains the file. The RPE gene set was analyzed by PAINT and a graphic CIM was generated as described in Methods. The y-axis lists the Ensembl Gene identifiers for each gene and the x-axis lists the TRANSFAC TRE identifiers for each TRE found at least once in the first intron region. Genes listed along the y-axis are divided into two clusters that are either up-regulated (blue) or down-regulated (green) during EMT of RPE cells. TREs listed along the x-axis are clustered according to related occurrence pattern calculated using Jaccard's coefficient. The elements within the matrix are color-coded based upon the p-value of each TRE found in the regulatory regions of the genes. A red dot represents a TRE that is statistically significant and therefore over-represented in our gene set, while a blue dot signifies an under-represented TRE and a grey dot stands for a TRE with no statistical significance in our gene list. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 6: Statistically over-represented TREs derived from the PAINT analysis of the mouse promoter regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-6.pdf">Appendix 6</ext-link>.” This will initiate the download of a pdf archive that contains the file. Represented for each TRE identified in the promoter sequences is the cluster in which the TRE was determined to be enriched, the TRANSFAC identifier and the p-value of the TRE in both gene sets.</p> <title>Appendix 7: Statistically over-represented TREs derived from the PAINT analysis of the mouse first intron regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-7.pdf">Appendix 7</ext-link>.” This will initiate the download of a pdf archive that contains the file. Represented for each TRE identified in the first intron sequences is the cluster in which the TRE was determined to be enriched, the TRANSFAC identifier and the p-value of the TRE in both gene sets.</p> <title>Appendix 8: Statistically over-represented TREs derived from the PAINT analysis of the rat promoter regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-8.pdf">Appendix 8</ext-link>.” This will initiate the download of a pdf archive that contains the file. Represented for each TRE identified in the promoter sequences is the cluster in which the TRE was determined to be enriched, the TRANSFAC identifier and the p-value of the TRE in both gene sets.</p> <title>Appendix 9: Statistically over-represented TREs derived from the PAINT analysis of the rat first intron regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-9.pdf">Appendix 9</ext-link>.” This will initiate the download of a pdf archive that contains the file. Represented for each TRE identified in the first intron sequences is the cluster in which the TRE was determined to be enriched, the TRANSFAC identifier and the p-value of the TRE in both gene sets.</p> <title>Appendix 10: Statistically over-represented TREs derived from the PAINT analysis of the chicken promoter regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-10.pdf">Appendix 10</ext-link>.” This will initiate the download of a pdf archive that contains the file. Represented for each TRE identified in the promoter sequences is the cluster in which the TRE was determined to be enriched, the TRANSFAC identifier and the p-value of the TRE in both gene sets.</p> <title>Appendix 11: Statistically over-represented TREs derived from the PAINT analysis of the chicken first intron regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-11.pdf">Appendix 11</ext-link>.” This will initiate the download of a pdf archive that contains the file. Represented for each TRE identified in the first intron sequences is the cluster in which the TRE was determined to be enriched, the TRANSFAC identifier and the p-value of the TRE in both gene sets.</p> <title>Appendix 12. Complete Candidate Interaction Matrix (CIM) for transcription response elements from PAINT analysis of mouse gene promoters.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-12.pdf">Appendix 12</ext-link>.” This will initiate the download of a pdf archive that contains the file. The RPE gene set was analyzed by PAINT and a graphic CIM was generated as described in Methods. The y-axis lists the Ensembl Gene identifiers for each gene and the x-axis lists the TRANSFAC TRE identifiers for each TRE found at least once in the promoter region. Genes listed along the y-axis are divided into two clusters that are either up-regulated (blue) or down-regulated (green) during EMT of RPE cells. TREs listed along the x-axis are clustered according to related occurrence pattern calculated using Jaccard's coefficient. The elements within the matrix are color-coded based upon the p-value of each TRE found in the regulatory regions of the genes. A red dot represents a TRE that is statistically significant and therefore over-represented in our gene set, while a blue dot signifies an under-represented TRE and a grey dot stands for a TRE with no statistical significance in our gene list. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 13. Complete Candidate Interaction Matrix (CIM) for transcription response elements from PAINT analysis of mouse gene first intron.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-13.pdf">Appendix 13</ext-link>.” This will initiate the download of a pdf archive that contains the file. The RPE gene set was analyzed by PAINT and a graphic CIM was generated as described in Methods. The y-axis lists the Ensembl Gene identifiers for each gene and the x-axis lists the TRANSFAC TRE identifiers for each TRE found at least once in the first intron region. Genes listed along the y-axis are divided into two clusters that are either up-regulated (blue) or down-regulated (green) during EMT of RPE cells. TREs listed along the x-axis are clustered according to related occurrence pattern calculated using Jaccard's coefficient. The elements within the matrix are color-coded based upon the p-value of each TRE found in the regulatory regions of the genes. A red dot represents a TRE that is statistically significant and therefore over-represented in our gene set, while a blue dot signifies an under-represented TRE and a grey dot stands for a TRE with no statistical significance in our gene list. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 14. Complete Candidate Interaction Matrix (CIM) for transcription response elements from PAINT analysis of rat gene promoters.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-14.pdf">Appendix 14</ext-link>.” This will initiate the download of a pdf archive that contains the file. The RPE gene set was analyzed by PAINT and a graphic CIM was generated as described in Methods. The y-axis lists the Ensembl Gene identifiers for each gene and the x-axis lists the TRANSFAC TRE identifiers for each TRE found at least once in the promoter region. Genes listed along the y-axis are divided into two clusters that are either up-regulated (blue) or down-regulated (green) during EMT of RPE cells. TREs listed along the x-axis are clustered according to related occurrence pattern calculated using Jaccard's coefficient. The elements within the matrix are color-coded based upon the p-value of each TRE found in the regulatory regions of the genes. A red dot represents a TRE that is statistically significant and therefore over-represented in our gene set, while a blue dot signifies an under-represented TRE and a grey dot stands for a TRE with no statistical significance in our gene list. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 15. Complete Candidate Interaction Matrix (CIM) for transcription response elements from PAINT analysis of rat gene first introns.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-15.pdf">Appendix 15</ext-link>.” This will initiate the download of a pdf archive that contains the file. The RPE gene set was analyzed by PAINT and a graphic CIM was generated as described in Methods. The y-axis lists the Ensembl Gene identifiers for each gene and the x-axis lists the TRANSFAC TRE identifiers for each TRE found at least once in the first introns region. Genes listed along the y-axis are divided into two clusters that are either up-regulated (blue) or down-regulated (green) during EMT of RPE cells. TREs listed along the x-axis are clustered according to related occurrence pattern calculated using Jaccard's coefficient. The elements within the matrix are color-coded based upon the p-value of each TRE found in the regulatory regions of the genes. A red dot represents a TRE that is statistically significant and therefore over-represented in our gene set, while a blue dot signifies an under-represented TRE and a grey dot stands for a TRE with no statistical significance in our gene list. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 16. Complete Candidate Interaction Matrix (CIM) for transcription response elements from PAINT analysis of chicken gene promoters.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-16.pdf">Appendix 16</ext-link>.” This will initiate the download of a pdf archive that contains the file. The RPE gene set was analyzed by PAINT and a graphic CIM was generated as described in Methods. The y-axis lists the Ensembl Gene identifiers for each gene and the x-axis lists the TRANSFAC TRE identifiers for each TRE found at least once in the promoter region. Genes listed along the y-axis are divided into two clusters that are either up-regulated (blue) or down-regulated (green) during EMT of RPE cells. TREs listed along the x-axis are clustered according to related occurrence pattern calculated using Jaccard's coefficient. The elements within the matrix are color-coded based upon the p-value of each TRE found in the regulatory regions of the genes. A red dot represents a TRE that is statistically significant and therefore over-represented in our gene set, while a blue dot signifies an under-represented TRE and a grey dot stands for a TRE with no statistical significance in our gene list. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 17. Complete Candidate Interaction Matrix (CIM) for transcription response elements from PAINT analysis of chicken gene first intron.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-17.pdf">Appendix 17</ext-link>.” This will initiate the download of a pdf archive that contains the file. The RPE gene set was analyzed by PAINT and a graphic CIM was generated as described in Methods. The y-axis lists the Ensembl Gene identifiers for each gene and the x-axis lists the TRANSFAC TRE identifiers for each TRE found at least once in the first intron region. Genes listed along the y-axis are divided into two clusters that are either up-regulated (blue) or down-regulated (green) during EMT of RPE cells. TREs listed along the x-axis are clustered according to related occurrence pattern calculated using Jaccard's coefficient. The elements within the matrix are color-coded based upon the p-value of each TRE found in the regulatory regions of the genes. A red dot represents a TRE that is statistically significant and therefore over-represented in our gene set, while a blue dot signifies an under-represented TRE and a grey dot stands for a TRE with no statistical significance in our gene list. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 18. Transcriptional regulatory network diagram for TREs associated with mouse promoter regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-18.pdf">Appendix 18</ext-link>.” This will initiate the download of a pdf archive that contains the file. The graphical representation was derived using GraphViz as described in Methods. The ellipses and diamonds represent individual genes divided into up-regulated (blue) and down-regulated (green) clusters. The boxes represent TREs, with arrows indicating gene-TRE associations. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 19. Transcriptional regulatory network diagram for TREs associated with rat promoter regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-19.pdf">Appendix 19</ext-link>.” This will initiate the download of a pdf archive that contains the file. The graphical representation was derived using GraphViz as described in Methods. The ellipses and diamonds represent individual genes divided into up-regulated (blue) and down-regulated (green) clusters. The boxes represent TREs, with arrows indicating gene-TRE associations. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 20. Transcriptional regulatory network diagram for TREs associated with chicken promoter regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-20.pdf">Appendix 20</ext-link>.” This will initiate the download of a pdf archive that contains the file. The graphical representation was derived using GraphViz as described in Methods. The ellipses and diamonds represent individual genes divided into up-regulated (blue) and down-regulated (green) clusters. The boxes represent TREs, with arrows indicating gene-TRE associations. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 21. Transcriptional regulatory network diagram for TREs associated with mouse first intron regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-21.pdf">Appendix 21</ext-link>.” This will initiate the download of a pdf archive that contains the file. The graphical representation was derived using GraphViz as described in Methods. The ellipses and diamonds represent individual genes divided into up-regulated (blue) and down-regulated (green) clusters. The boxes represent TREs, with arrows indicating gene-TRE associations. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 22. Transcriptional regulatory network diagram for TREs associated with rat first intron regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-22.pdf">Appendix 22</ext-link>.” This will initiate the download of a pdf archive that contains the file. The graphical representation was derived using GraphViz as described in Methods. The ellipses and diamonds represent individual genes divided into up-regulated (blue) and down-regulated (green) clusters. The boxes represent TREs, with arrows indicating gene-TRE associations. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 23. Transcriptional regulatory network diagram for TREs associated with chicken first intron regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-23.pdf">Appendix 23</ext-link>.” This will initiate the download of a pdf archive that contains the file. The graphical representation was derived using GraphViz as described in Methods. The ellipses and diamonds represent individual genes divided into up-regulated (blue) and down-regulated (green) clusters. The boxes represent TREs, with arrows indicating gene-TRE associations. Note: This is a large-format figure and should be viewed at enhanced magnification to read the text in specific sections of the figure.</p> <title>Appendix 24. Frequency analysis of TRE representation in mouse promoter regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-24.pdf">Appendix 24</ext-link>.” This will initiate the download of a pdf archive that contains the file. Frequency of occurrence of each TRE in the mouse gene promoter regions was determined from the PAINT analysis as described in Methods. The y-axis indicates the frequency of each TRE among the up-regulated (blue) or down-regulated (green) gene clusters as well as among the full background gene set (black).</p> <title>Appendix 25. Frequency analysis of TRE representation in mouse first intron regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-25.pdf">Appendix 25</ext-link>.” This will initiate the download of a pdf archive that contains the file. Frequency of occurrence of each TRE in the mouse gene first intron regions was determined from the PAINT analysis as described in Methods. The y-axis indicates the frequency of each TRE among the up-regulated (blue) or down-regulated (green) gene clusters as well as among the full background gene set (black).</p> <title>Appendix 26. Frequency analysis of TRE representation in rat promoter regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-26.pdf">Appendix 26</ext-link>.” This will initiate the download of a pdf archive that contains the file. Frequency of occurrence of each TRE in the rat gene promoter regions was determined from the PAINT analysis as described in Methods. The y-axis indicates the frequency of each TRE among the up-regulated (blue) or down-regulated (green) gene clusters as well as among the full background gene set (black).</p> <title>Appendix 27. Frequency analysis of TRE representation in rat first intron regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-27.pdf">Appendix 27</ext-link>.” This will initiate the download of a pdf archive that contains the file. Frequency of occurrence of each TRE in the rat gene first intron regions was determined from the PAINT analysis as described in Methods. The y-axis indicates the frequency of each TRE among the up-regulated (blue) or down-regulated (green) gene clusters as well as among the full background gene set (black).</p> <title>Appendix 28. Frequency analysis of TRE representation in chicken promoter regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-28.pdf">Appendix 28</ext-link>.” This will initiate the download of a pdf archive that contains the file. Frequency of occurrence of each TRE in the chicken gene promoter regions was determined from the PAINT analysis as described in Methods. The y-axis indicates the frequency of each TRE among the up-regulated (blue) or down-regulated (green) gene clusters as well as among the full background gene set (black).</p> <title>Appendix 29. Frequency analysis of TRE representation in chicken first intron regions.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-29.pdf">Appendix 29</ext-link>.” This will initiate the download of a pdf archive that contains the file. Frequency of occurrence of each TRE in the chicken gene first intron regions was determined from the PAINT analysis as described in Methods. The y-axis indicates the frequency of each TRE among the up-regulated (blue) or down-regulated (green) gene clusters as well as among the full background gene set (black).</p> <title>Appendix 30: Criteria table for transcription response element inclusion in gene promoter region regulatory network models.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-30.pdf">Appendix 30</ext-link>.” This will initiate the download of a pdf archive that contains the file. Represented for each TRE identified in the promoter regions of each gene are the TRANSFAC identifier name, the p-value, the ECF (evolutionary conservation factor), and the FR (frequency ratio) as defined in Methods. TREs included in this table are those with p-value &lt;0.1, ECF =2, and FR &gt;3.</p> <title>Appendix 31: Criteria table for transcription response element inclusion in gene first intron regulatory network models.</title> <p>To access the data, click or select the words “<ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-31.pdf">Appendix 31</ext-link>.” This will initiate the download of a pdf archive that contains the file. Represented for each TRE identified in the first intron regions of each gene are the TRANSFAC identifier name, the p-value, the ECF (evolutionary conservation factor), and the FR (frequency ratio) as defined in Methods. TREs included in this table are those with p-value &lt;0.1, ECF =2, and FR &gt;3.</p> <title>Acknowledgments</title> <p>Grant support: NIH RO1EY06658 (G.B.G., C.H.P.), NIHT32ES07282 (G.B.G., C.H.P.), DARPA F30602–01–2–0578 (R.V., P.C., G.E.G.), NIH RO1EY012042 (N.J.P.), Thomas Jefferson University-Wills Eye Hospital Pilot Program Research Award 080–01000-A92201 (G.B.G., N.J.P., C.H.P.), and a Foerderer Fellowship from Thomas Jefferson University College of Graduate Studies (C.H.P.).</p> <p>Portions of this work have been previously presented in either poster or paper presentation format at the following meetings: 2004 Thomas Jefferson University chapter of Sigma-Xi Scientific Research Society, Annual Student Research Symposium, 2004 Annual Meeting of The Association for Research in Vision and Ophthalmology, 27th Annual Meeting of the Philadelphia chapter of the Society for Neuroscience, 44th Annual Meeting of The American Society for Cell Biology, 2005 Thomas Jefferson University chapter of Sigma-Xi Scientific Research Society, Annual Student Research Symposium, 2005 Gordon Conference on Cell Contact and Adhesion, 3rd Annual Retreat of the Greater Philadelphia Bioinformatics Alliance, 2006 Lennox K. Black International Prize for Excellence in Medicine, 4th Symposium: Molecular Basis of Neurologic Disease, 2007 Thomas Jefferson University chapter of Sigma-Xi Scientific Research Society, Annual Student Research Symposium, Greater Philadelphia Bioinformatics Alliance Spring Retreat 2007, 2007 Annual Meeting of The Association for Research in Vision and Ophthalmology.</p>
<fig id="f1" fig-type="figure" position="float"><label>Figure 1</label><caption><p>Candidate interaction matrix for statistically enriched transcription response elements from promoter analysis and interaction network toolset analysis of human gene promoters. The retinal pigment epithelium (RPE) gene set was analyzed by promoter analysis and interaction network toolset (PAINT) and a graphic candidate interaction matrix (CIM) was generated as described in Methods. The y-axis lists the Ensembl Gene identifiers for each gene and the x-axis lists the TRANSFAC identifiers for each transcription response element (TRE) found at least once in the promoter region of one or more genes. Genes listed along the y-axis are divided into two clusters that are either upregulated (blue) or down-regulated (green) during epithelial-mesenchymal transformation (EMT) of RPE cells. TREs listed along the x-axis are clustered according to related occurrence pattern calculated using Jaccard's coefficient. The elements within the matrix are color-coded based on the p-value of each TRE found in the regulatory regions of the genes. A red dot represents a TRE that is statistically significant and therefore over-represented in our gene set, while a blue dot signifies an under-represented TRE and a gray dot stands for a TRE with no statistical significance in our gene list. This figure represents the subset of enriched TREs for the human genome; the full CIMs for human and other genomes analyzed are shown in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-4.pdf">Appendix 4</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-5.pdf">Appendix 5</ext-link>, <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-12.pdf">Appendix 12</ext-link>, and <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-13.pdf">Appendix</ext-link><ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a999/app-13.pdf">13</ext-link>.</p></caption></fig> <fig id="f2" fig-type="figure" position="float"><label>Figure 2</label><caption><p>Candidate interaction matrix for statistically enriched transcription response elements from promoter analysis and interaction network toolset analysis of human gene first introns. The retinal pigment epithelium (RPE) gene set was analyzed by promoter analysis and interaction network toolset (PAINT) and a graphic candidate interaction matrix (CIM) was generated as described in Methods. The y-axis lists the Ensembl Gene identifiers for each gene and the x-axis lists the TRANSFAC identifiers for each transcription response element (TRE) found at least once in the first intron region. Genes listed along the y-axis are divided into two clusters that are either upregulated (blue) or down-regulated (green) during EMT of RPE cells. TREs listed along the x-axis are clustered according to related occurrence pattern calculated using Jaccard's coefficient. The elements within the matrix are color-coded based on the p-value of each TRE found in the regulatory regions of the genes. A red dot represents a TRE that is statistically significant and therefore over-represented in our gene set, while a blue dot signifies an under-represented TRE and a gray dot stands for a TRE with no statistical significance in our gene list. This figure represents the subset of enriched TREs for the human genome; the full CIMs for human and other genomes analyzed are shown in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-14.pdf">Appendix 14</ext-link> through <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-17.pdf">Appendix 17</ext-link>.</p></caption></fig> <fig id="f3" fig-type="figure" position="float"><label>Figure 3</label><caption><p>Transcriptional regulatory network diagram for transcription response element associated with human promoter regions. The graphical representation was derived using GraphViz as described in Methods. The ellipses and diamonds represent individual genes divided into upregulated (blue) and down-regulated (green) clusters. The boxes represent TREs, with arrows indicating gene-TRE associations. Corresponding network diagrams for mouse, rat and chick promoter regions are shown in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-18.pdf">Appendix 18</ext-link> through <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-20.pdf">Appendix 20</ext-link>.</p></caption></fig> <fig id="f4" fig-type="figure" position="float"><label>Figure 4</label><caption><p>Transcriptional regulatory network diagram for transcription response element associated with human first intron regions. The graphical representation was derived using GraphViz as described in Methods. The ellipses and diamonds represent individual genes divided into upregulated (blue) and down-regulated (green) clusters. The boxes represent TREs, with arrows indicating gene-TRE associations. Corresponding network diagrams for mouse, rat and chick first intron regions are shown in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-21.pdf">Appendix 21</ext-link> through <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-23.pdf">Appendix 23</ext-link>.</p></caption></fig> <fig id="f5" fig-type="figure" position="float"><label>Figure 5</label><caption><p>Frequency analysis of transcription response element representation in human promoter regions. Frequency of occurrence of each transcription response element (TRE) in the human gene promoter regions was determined from the promoter analysis and interaction network toolset (PAINT) analysis as described in Methods. The y-axis indicates the frequency of each TRE among the upregulated (blue) or down-regulated (green) gene clusters as well as among the full background gene set (black). The x-axis indicates the over-represented TREs, ordered by increasing p-value. Corresponding frequency analyses for mouse, rat and chick promoter regions are shown in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-24.pdf">Appendix 24</ext-link> through <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-26.pdf">Appendix 26</ext-link>.</p></caption></fig> <fig id="f6" fig-type="figure" position="float"><label>Figure 6</label><caption><p>Frequency analysis of transcription response element representation in human first intron regions. Frequency of occurrence of each transcription response element (TRE) in the human gene first intron regions was determined from the PAINT analysis as described in Methods. The y-axis indicates the frequency of each TRE among the upregulated (blue) or down-regulated (green) gene clusters as well as among the full background gene set (black). The x-axis indicates the over-represented TREs, ordered by increasing p-value. Corresponding frequency analyses for mouse, rat and chick first intron regions are shown in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-27.pdf">Appendix 27</ext-link> through <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-29.pdf">Appendix 29</ext-link>.</p></caption></fig> <fig id="f7" fig-type="figure" position="float"><label>Figure 7</label><caption><p>Archetypal cross-species gene regulatory region models of undifferentiated and differentiated gene clusters. These models incorporate transcription response element (TREs) that were found to be over-represented in results of both the human and chicken, as well as either the mouse or rat, promoter analysis and interaction network toolset (PAINT) analysis. TSS represents transcriptional start site.</p></caption></fig> <fig id="f8" fig-type="figure" position="float"><label>Figure 8</label><caption><p>Human gene regulatory region models of undifferentiated and differentiated gene clusters. These models incorporate transcription response elements (TREs) that were found to be over-represented in results of the human, and either the chicken, mouse or rat, promoter analysis and interaction network toolset (PAINT) analysis. TSS represents transcriptional start site.</p></caption></fig> <fig id="f9" fig-type="figure" position="float"><label>Figure 9</label><caption><p>Gene regulatory region models for specific reciprocally-regulated gene pairs. Models for the paired genes that are reciprocally regulated during EMT of RPE cells models including N- and R-cadherin (<bold>A</bold>), α-SMA and RPE-65 (<bold>B</bold>), and MCT-3 and −4 (<bold>C</bold>). Models were constructed by including only those TREs that are over-represented in both the human and chicken, as well as either the mouse or rat, PAINT analysis.</p></caption></fig> <fig id="f10" fig-type="figure" position="float"><label>Figure 10</label><caption><p>Reverse transcriptase polymerase chain reaction analysis of markers during retinal pigment epithelium cell differentiation. mRNA was isolated from undifferentiated (<bold>A</bold>) or differentiated (<bold>B</bold>) ARPE-19 cells and subjected to reverse transcriptase polymerase chain reaction amplification to detect SMA (<bold>C</bold>) or RPE65 (<bold>D</bold>) as described in Methods. Phase-contrast micrographs represent undifferentiated (<bold>A</bold>) or differentiated (<bold>B</bold>) ARPE-19 cells after one week (<bold>A</bold>) or 52 weeks (<bold>B</bold>) of culture. <bold>C</bold> and <bold>D</bold> represent RT–PCR amplification of mRNA samples isolated from ARPE-19 cells maintained in culture for 1, 2, 5, 7, 10, 12, and 52 weeks, using primers to detect mRNA for either αSMA (<bold>C</bold>) or RPE65 (<bold>D</bold>). <bold>E</bold> represents RT–PCR amplification for a series of 25, 30, or 35 cycles to detect αSMA using mRNA isolated from ARPE-19 cells that are undifferentiated (U) differentiated (<bold>D</bold>). The first lane in <bold>C-E</bold> represents a DNA standard ladder of 300, 400, and 500 bp.</p></caption></fig> <fig id="f11" fig-type="figure" position="float"><label>Figure 11</label><caption><p>Reverse transcriptase polymerase chain reaction amplification of transcription response element mRNAs during ARPE-19 retinal pigment epithelium cell differentiation RNA was isolated from undifferentiated and differentiated ARPE-19 cells and subjected to RT–PCR analysis to detected transcription response element (TRE) mRNAs as described in Methods. In <bold>A</bold>, all reactions were performed for 40 cycles, where lane 1 represents DNA standards, lane 15 represents the positive control primers for GADPH, and lane 16 is the negative control with no mRNA template. The intervening lanes in <bold>A</bold> represent primers specific for the following TFs: 2 Core binding factor, 3 E2F1, 4 Evi-1, 5 GATA1, 6 HNF-1, 7 IRF-1, 8 Nkx2–5, 9 NKX3A, 10 Oct-1, 11 SMAD3, 12 SREBP-1, 13 TFE3, 14 v-Myb. In <bold>B</bold>, semi-quantitative RT–PCR was also done for either 30, 35, or 40 cycles as indicated using primers specific for GATA-1, IRF-1, or SMAD3. The first lane in <bold>A</bold> represents a standard DNA ladder at 300, 400, and 500 bp, while in <bold>B</bold> the DNA standards are at 400 and 500 bp.</p></caption></fig> <fig id="f12" fig-type="figure" position="float"><label>Figure 12</label><caption><p>Reverse transcriptase polymerase chain reaction amplification of transcription response elements mRNAs during chick embryo retinal pigment epithelium cell differentiation RNA was isolated from undifferentiated cultured chick embryo retinal pigment epithelium (RPE) cells or differentiated fresh RPE tissue and subjected to RT–PCR analysis to detect transcription response elements (TREs) mRNAs as described in Methods. Lanes 1 and 11 represent DNA standard ladders at 300, 400, 500, 600, and 700 bp, lane 10 represents the positive control for GADPH, and the remaining lanes represent primers specific for the following TFs: 2 FoxD3, 3 AML-1, 4 HNF-3α, 5 HNF-1, 6 E2F1, 7 DP1, 8 TFII-I, 9 SREBP-1.</p></caption></fig> <fig id="f13" fig-type="figure" position="float"><label>Figure 13</label><caption><p>Reverse transcriptase polymerase chain reaction amplification of EMT-associated transcription response elements mRNAs during ARPE-19 and chick embryo retinal pigment epithelium cell differentiation RNA was isolated from undifferentiated and differentiated ARPE-19 or chick embryo retinal pigment epithelium (cRPE) cells and subjected to RT–PCR analysis to detect transcription response elements (TRE) mRNAs as described in Methods. Lanes 1 and 16 represent DNA standard ladders at 300, 400, and 500 bp, lanes 8 and 15 represent the positive controls for GADPH, and the remaining lanes represent primers specific for the following TFs: 2 and 9, Slug; 3 and 10, Snail; 4 and 11, Twist; 5 and 12, SIP-1; 6 and 13, SMAD-2; 7 and 14, LEF-1.</p></caption></fig> <fig id="f14" fig-type="figure" position="float"><label>Figure 14</label><caption><p>Comprehensive cross-species models for reciprocal regulation of genes during retinal pigment epithelium cell differentiation Models for regulatory regions, including promoters and first introns, of differentiated and undifferentiated gene clusters, were constructed as described in the text. transcription response elements (TREs) inclusion criteria indicated in these models are frequency ratio, evolutionary conservation factor, and RT–PCR detection of mRNA expression. Symbols representing TREs indicates passage of corresponding criteria filter, as indicated in key, by TRE, filled boxes signify that a TRE has passed all criteria filters.</p></caption></fig>
<table-wrap id="t1" position="float"><label>Table 1</label><caption><title>Gene ontology list for target sequences used in PAINT analysis</title></caption><table frame="hsides" rules="groups"><col width="134" span="1"/><col width="100" span="1"/><col width="80" span="1"/><col width="81" span="1"/><col width="72" span="1"/><col width="80" span="1"/><thead><tr><th valign="bottom" align="left" scope="col" rowspan="1" colspan="1"><bold>Gene Ontology</bold></th><th valign="bottom" align="left" scope="col" rowspan="1" colspan="1"><bold>Gene name</bold></th><th valign="bottom" align="left" scope="col" rowspan="1" colspan="1"><bold>Cluster</bold></th><th valign="bottom" align="left" scope="col" rowspan="1" colspan="1"><bold>Gene Ontology</bold></th><th valign="bottom" align="left" scope="col" rowspan="1" colspan="1"><bold>Gene name</bold></th><th valign="bottom" align="left" scope="col" rowspan="1" colspan="1"><bold>Cluster</bold></th></tr></thead><tbody><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Growth factor binding<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">IGFBP-1<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1">Cell cycle effectors<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">CDC25A<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">IGFBP-3<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Cyclin H<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">EMAP II<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Cdc2-related protein kinase<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">IL-13<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">G(1)/G(S)/G(2) beta2<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"><hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"><hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Cell death<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Caspase 4<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1">Cell adhesion<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">I-CAM-1<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">ICE LAP3<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">V-CAM-1<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"><hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"><hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">E-selectin<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Cell organization<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">alpha-SMA<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Integrin alpha4<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">ZO1<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Integrin beta5<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Thymosin beta10<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">N-cadherin<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"><hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"><hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Fibronectin<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Signal transduction<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">STK-2<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">CD44 antigen<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">CAK<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Integrin alpha5<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">MAPK3<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Integrin beta4<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">PDGF-B<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">R-cadherin<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">CD33 antigen<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">RGS19IP1<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1">Cell metabolism<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">CLK-1<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">FADK2<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">STK-1<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Tyk2<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Tyrosinase<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">RAB5A<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">PEDF<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">PKCgamma17<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">RPE65<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">IGF-1<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">TYRP2<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">VEGF<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">MCP-1<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1">Ion binding<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">CAM IV<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Neuromodulin<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">SPARC/Osteonectin<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">FADK2<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">PAK-C alpha<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1">Transcription factors<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">CRE BP-1<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"><hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">FGFR-1<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">TFAP 2<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"><hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">FGFR-3<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">HDGF<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1">Proliferation/differentiation<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">STAT6<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"><hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"><hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"><hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">GDF-1<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Inter/intra cellular transport<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">MAL<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">MCT-4<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Undifferentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1"> <hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">MCT-3<hr/></td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated<hr/></td></tr><tr><td valign="bottom" align="left" scope="row" rowspan="1" colspan="1"> </td><td valign="bottom" align="left" rowspan="1" colspan="1"> </td><td valign="bottom" align="left" rowspan="1" colspan="1"> </td><td valign="bottom" align="left" rowspan="1" colspan="1"> </td><td valign="bottom" align="left" rowspan="1" colspan="1">Bestrophin</td><td valign="bottom" align="left" rowspan="1" colspan="1">Differentiated</td></tr></tbody></table></table-wrap> <table-wrap id="t2" position="float"><label>Table 2</label><caption><title>Criteria table for transcription response element inclusion in gene regulatory network models</title></caption><table frame="hsides" rules="groups"><col width="153" span="1"/><col width="99" span="1"/><col width="81" span="1"/><col width="72" span="1"/><thead><tr><th valign="top" align="center" scope="col" rowspan="1" colspan="1"><bold>TRE</bold></th><th valign="top" align="center" scope="col" rowspan="1" colspan="1"><bold>p-value</bold></th><th valign="top" align="center" scope="col" rowspan="1" colspan="1"><bold>Evolutionary conservation factor</bold></th><th valign="top" align="center" scope="col" rowspan="1" colspan="1"><bold>Frequency ratio</bold></th></tr></thead><tbody><tr><td colspan="4" valign="top" align="left" scope="col" rowspan="1">Promoter<hr/></td></tr><tr><td valign="top" align="center" scope="row" rowspan="1" colspan="1">E2F-1/V$E2F1_Q3_01<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">0.00095<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">3.725<hr/></td></tr><tr><td valign="top" align="center" scope="row" rowspan="1" colspan="1">Oct-1/V$OCT1_02<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">0.00801<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">3<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">4.25<hr/></td></tr><tr><td valign="top" align="center" scope="row" rowspan="1" colspan="1">Nkx2–5/V$NKX25_02<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">0.05404<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2.3<hr/></td></tr><tr><td valign="top" align="center" scope="row" rowspan="1" colspan="1">Poly A/V$LDSPOLYA_B<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">0.07291<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2.16<hr/></td></tr><tr><td valign="top" align="center" scope="row" rowspan="1" colspan="1">IRF/V$IRF_Q6_01<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">0.08969<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2.33<hr/></td></tr><tr><td valign="top" align="center" scope="row" rowspan="1" colspan="1">COMP1/V$COMP1_01<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">0.03569<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2.52<hr/></td></tr><tr><td colspan="4" valign="top" align="left" scope="col" rowspan="1"><hr/></td></tr><tr><td colspan="4" valign="top" align="left" scope="col" rowspan="1">First intron<hr/></td></tr><tr><td colspan="4" valign="top" align="left" scope="col" rowspan="1"><hr/></td></tr><tr><td valign="top" align="center" scope="row" rowspan="1" colspan="1">Evi-1/V$EVI1_05<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">0.00003<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">13.86<hr/></td></tr><tr><td valign="top" align="center" scope="row" rowspan="1" colspan="1">Oct-1/V$OCT1_02<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">0.00043<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">6.32<hr/></td></tr><tr><td valign="top" align="center" scope="row" rowspan="1" colspan="1">GR/V$GRE_C<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">0.00045<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">14.6<hr/></td></tr><tr><td valign="top" align="center" scope="row" rowspan="1" colspan="1">SMAD-3/V$SMAD3_Q6<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">0.00071<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">3<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">3.18<hr/></td></tr><tr><td valign="top" align="center" scope="row" rowspan="1" colspan="1">HNF-1/V$HNF1_01<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">0.00918<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">2<hr/></td><td valign="top" align="center" rowspan="1" colspan="1">7.35<hr/></td></tr><tr><td valign="top" align="center" scope="row" rowspan="1" colspan="1">AP-1/V$AP1_C</td><td valign="top" align="center" rowspan="1" colspan="1">0.01893</td><td valign="top" align="center" rowspan="1" colspan="1">2</td><td valign="top" align="center" rowspan="1" colspan="1">4.87</td></tr></tbody></table></table-wrap>
<table-wrap-foot><p>The three columns indicate the gene name, gene ontology cluster assignment, and inclusion into either the differentiated or undifferentiated retinal pigment epithelium expression cluster. The genes were organized within the ontology clusters according to The Gene Ontology Consortium (2000).</p></table-wrap-foot> <table-wrap-foot><p>Represented for each transcription response elements (TRE) identified in the promoter and intron regions of each gene are the TRANSFAC identifier name, the p-value, the evolutionary conservation factor (ECF), and the frequency ratio (FR) as defined in Methods. TREs included in this table are those with p-value &lt;0.1, ECF=2, and FR&gt;3. The full table of values for all TREs analyzed is contained in <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-30.pdf">Appendix 30</ext-link> (promoters) and <ext-link ext-link-type="uri" xlink:href="http://www.molvis.org/molvis/v14/a170/app-31.pdf">Appendix 31</ext-link> (introns).</p></table-wrap-foot>
<graphic xlink:href="mv-v14-1414-f1"/> <graphic xlink:href="mv-v14-1414-f2"/> <graphic xlink:href="mv-v14-1414-f3"/> <graphic xlink:href="mv-v14-1414-f4"/> <graphic xlink:href="mv-v14-1414-f5"/> <graphic xlink:href="mv-v14-1414-f6"/> <graphic xlink:href="mv-v14-1414-f7"/> <graphic xlink:href="mv-v14-1414-f8"/> <graphic xlink:href="mv-v14-1414-f9"/> <graphic xlink:href="mv-v14-1414-f10"/> <graphic xlink:href="mv-v14-1414-f11"/> <graphic xlink:href="mv-v14-1414-f12"/> <graphic xlink:href="mv-v14-1414-f13"/> <graphic xlink:href="mv-v14-1414-f14"/>
[{"label": ["2"], "citation": ["Grunwald G. Structure and Function of the Retinal Pigment Epithelium.In: Tasman W, Jaeger EA, editors. Foundations of Clinical Ophthalmology. Philadelphia(PA): Lippincott, Williams and Wilkins; 2004. p. 1\u201321."]}, {"label": ["17"], "surname": ["Benjamini", "Hochberg"], "given-names": ["Y", "Y"], "article-title": ["Controlling the false discovery rate: a practical and powerful approach to multiple testing. "], "source": ["J R Statist Soc B"], "year": ["1995"], "volume": ["57"], "fpage": ["289"], "lpage": ["300"]}, {"label": ["18"], "surname": ["Gansner", "North"], "given-names": ["ER", "SC"], "article-title": ["An open graph visualization system and its applications to software engineering."], "source": ["Softw Pract Exper"], "year": ["1999"], "volume": ["00"], "issue": ["S1"], "fpage": ["1"], "lpage": ["5"]}, {"label": ["26"], "surname": ["Grunwald"], "given-names": ["GB"], "article-title": ["Cadherin cell adhesion molecules in retinal development and pathology."], "source": ["Progress in Retinal and Eye Research"], "year": ["1996"], "volume": ["15"], "fpage": ["363"], "lpage": ["92"]}]
Mol Vis. 2008 Aug 4; 14:1414-1428