Molecular Vision 2012; 18:2961-2975
Received 17 July 2012 | Accepted 13 December 2012 | Published 15 December 2012
The first two authors contributed equally to this work.
1Department of Ophthalmology, the Second Affiliated Hospital of Harbin Medical University, Harbin, China; 2Department of Ophthalmology, Inner Mongolia Forestry general Hospital, Yakeshi, Inner Mongolia, China
Correspondence to: Shaomin Peng, Department of Ophthalmology, the Second Affiliated Hospital of Harbin Medical University, Harbin, 150086, China; Phone: +86 0451-86605851; Fax: +86 0451-86605757; email: email@example.com
Purpose: The molecular mechanisms associated with human retinal pigment epithelium (RPE) development constitute the basis for cell replacement therapy for the treatment of retinal degenerative diseases. In the current study, gene expression analysis of the human fetal RPE during development was performed and was compared with the human native RPE.
Methods: Microdissection of the human RPE at three time points (13 weeks and 16 weeks of gestation and in mature adult eyes) was performed, and total RNA was isolated. Equal amounts of RNA were pooled from two or three independent donor eyes for each time point in each group. Gene expression was analyzed by hybridization to microarray chips. Validation was accomplished by comparing the microarray expression proﬁles with quantitative real-time reverse transcriptase-polymerase chain reaction (qRT2-PCR) analysis of selected genes and by comparing selected expression proﬁles with predicted proﬁles based on previous studies.
Results: Of the 45,033 probe sets on the microarray, 30,736 were detected. A total of 3,498 differentially expressed genes could be clustered into eight patterns of expression that were statistically significant. Analysis of the expression patterns of genes coding for key functions (pigment synthesis, visual cycle, phagocytosis, adherens and tight junctions, and transcellular transport) indicated that the human RPE achieves a high degree of maturity during early pregnancy. Compared with 154 signature genes in the RPE, 148 candidate genes were identified in this study, including 53 downregulated genes and 5 upregulated genes. The qRT2-PCR results showed similar expression trends to those obtained by microarray analysis at the three time points.
Conclusions: This study demonstrated gene expression profiles in the human RPE during normal development. These ﬁndings indicate that the human RPE has different expression patterns than those of other animals. The results of this study may be helpful in furthering the understanding of the developmental processes occurring in humans and of the differentiation of RPE cells derived from human embryonic stem cells and from human induced pluripotent stem cells.
The retinal pigment epithelium (RPE) cells in human eyes form a quiescent, polarized epithelial monolayer located between the neural retina and the vascular choroid, and these cells serve to support and maintain the photoreceptor cells and other outer retinal cells via multiple mechanisms, including formation of the blood-retinal barrier, absorption of stray light, supply of nutrients to the neural retina, and regeneration of visual pigment, as well as the uptake and recycling of the shed outer segments of photoreceptors . Because of its important function in supporting photoreceptors, dysfunction and loss of the RPE leads to photoreceptor degeneration or apoptosis. Substantial evidence supports the notion that the dysfunction and death of RPE cells play a critical role in the pathogenesis of age-related macular degeneration (AMD) [2-4], which is the leading cause of blindness among the elderly in the developed world. As the population continues to age, the number of people in the United States with advanced AMD is expected to exceed 2.9 million by 2020 .
The emerging strategy of cell replacement therapy has provided a new approach to the treatment of AMD. Various types of dissociated RPE cells, such as immortalized adult RPE cell lines, human fetal RPE cells, and RPE cells derived from human embryonic stem cells (hESCs) and human induced pluripotent stem cells (hiPSCs), have been transplanted into the subretinal space in animal models of retinal degeneration caused by dysfunction of the RPE [6-13]. Moreover, Steven Schwartz reported on the first set of phase I clinical trials, which are ongoing, in which two patients were treated with RPE cells derived from hESCs . Although many of these studies have demonstrated the protection of photoreceptors and even improvement in visual function after transplantation, the visual protection and improvement have been found to be different due to the differential sources of cells; for example, fetal RPE cells may yield better results than adult RPE cells . Nevertheless, most studies have shown that transplanted cells die within 2 weeks to several months and that long-term survival is not achieved [12,13,16]. Thus, we are still faced with major obstacles to such cell-based therapies in clinical application. In addition, Liao’s studies showed that only 42 genes among 108 selected RPE signature genes are commonly shared by human fetal RPE, hESC-RPE, and hiPSC-RPE cells . Moreover, rather than a polygonal monolayer structure, which is characteristic of a functional RPE that is formed after transplantation, it is common to find clumps of cells that fail to integrate with the host monolayer.
Therefore, the culturing of
In the present study, we analyzed thousands of genes expressed in the different stages of the fetal and adult RPE using microarray and quantitative real-time reverse transcriptase–polymerase chain reaction (qRT2-PCR) analysis. These RPE gene expression profiles, combined with a bioinformatics analysis, enabled us to better understand the changes in gene expression and the intrinsic program of gene regulation that accompany human RPE development. Additionally, we expect that this research will contribute in the future to the optimization of cell replacement therapy and the differentiation of hESCs and hiPSCs.
This research followed the tenets of the Declaration of Helsinki on the use of human material for research and written informed consent was obtained before sample collection. Fetal human eyes (13 and 16 weeks of gestation) and adult eyes (40–60 years old), obtained from Harbin Medical University (Harbin, China), were dissected to obtain the posterior eye cups. After the retinas were carefully peeled off, the eye cups were treated with 2% dispase (Invitrogen-Gibco, Grand Island, NY) in PBS at 37 °C for 20 min. The RPE sheets were then gently peeled from the choroid with fine forceps under a dissecting microscope, and the sheets were stored in RNAlater (Qiagen, Valencia, CA).
Equal amounts of RNA were pooled from the three independent replicates used for each time point in each group. RNA was extracted using TRIzol (Invitrogen, Carlsbad, CA) and the RNeasy kit (Qiagen), according to the manufacturer’s instructions. The concentration and quality of RNA were determined using a bioanalyzer (Agilent Technologies, Palo Alto, CA) and/or a nanodrop ND-1000 spectrophotometer (Thermo Fisher Scientific Inc., Wilmington, DE). All of the samples had A260 /A280 ratios of total RNA of 2.0, and the ratio of the 28S to 18S rRNA bands was greater than 1.8.
Approximately 5 μg of total RNA in each sample set (pooled separately from
Raw data were extracted as pair files using NimbleScan software (version 2.5). The NimbleScan software’s implementation of
robust multichip averaging offers quantile normalization and background correction. The three gene summary files were imported
into Agilent GeneSpring software (version 11.0) for further analysis. Genes that had values greater than or equal to the lower
cutoff—100.0 in all samples (“All Targets Value”)—were chosen for data analysis. Differentially expressed genes were identified
using a combination of t-statistics, false discovery rates, and fold-change screening. Spots with twofold or more differences
in intensity between any two groups were marked as differentially expressed genes. The gene expression patterns were analyzed
using Short Time-series Expression Miner (STEM, version 1.2.2b) software. Gene ontology (GO) analysis
Seven genes that exhibited significant differential expression profiles based on cDNA microarray analysis were chosen for further validation by two-step qRT2–PCR analysis. Total RNA was extracted using an RNeasy Mini Kit (Qiagen). First-strand cDNA was synthesized from a 2 µg aliquot of the total RNA samples using the oligo-dT primer and MMLV reverse transcriptase (Epicenter, Madison, WI). Real-time RT–PCR was performed with a Rotor-Gene 3000 cycler (Corbett Research) and SybrGreen I (Invitrogen) according to the manufacturer’s protocol. The running conditions were as follows: 5 min at 95 °C, followed by 40 cycles of 10 s at 95 °C (denaturation), 15 s at 59 °C (annealing), and 20 s at 72 °C (extension). The experiments were performed in triplicate with a minimum of two biologic repetitions. Glyceraldehyde-3-phosphate dehydrogenase was used as a control to normalize the data. A total of 40 cycles were run, and the C(t) values were used to evaluate the relative expression levels of the tested genes. Standard curves were generated for each gene. To determine the specificity of DNA quantification, a melting curve was generated for each gene. Specific primers were designed with Primer software (version 5.0, Integrated DNA Technologies, Skokie, IL). The primers used for qRT2-PCR are listed in Table 1.
To generate gene expression profiles during human RPE development, three time points were chosen: 13 weeks and 16 weeks
Gene cluster analysis was performed using STEM version 1.2.2b software. STEM implements a novel clustering method that is based on a set of distinct and representative short temporal expression profiles, where each probe in the data set is assigned to the profile with the closest match. The genes expressed during RPE development could be classified into 16 clusters, and 8 were judged to be statistically significant. These 8 clusters included 3,498 probe sets. As shown in Figure 1, clusters 2 and 3 and clusters 12 and 13 represented the genes with consistently decreased expression and consistently increased expression throughout development, respectively. Clusters 7 and 9 represented the genes that were expressed at low levels at 13 weeks of gestation and in adulthood but that were upregulated in the 16-week-old fetus. In contrast to cluster 6, cluster 8 represented genes that were expressed at high levels in the 13-week-old fetus and in adulthood but were downregulated in the 16-week-old fetus. Genes with consistently decreasing expression during retinal development were more common than genes with consistently increasing expression. Gene ontology analysis revealed that these differentially expressed genes were associated with, for example, differentiation, signal transduction, protein translation, metabolism, DNA binding, and transcription.
With the aim of better characterizing the expression patterns of genes with critical functions within the retinal pigment epithelium, we gathered information on genes grouped by function, namely pigment synthesis, visual cycle, phagocytosis, adherens and tight junctions, and transcellular transport. Genes in these five groups were analyzed, and the results are shown in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, and Figure 9.
I. Pigment synthesis-- The pigment synthesis pathway is involved in eumelanin synthesis within melanosomes and melanosomal transporters. Ten genes were identified as participating in RPE pigment synthesis (Figure 2). All of the genes in this group exhibited significant downregulation during development. The gene with the largest decrease was TYR (140X), which promotes the first two rate-limiting steps of melanin synthesis.
II. Visual cycle-- This pathway serves to regenerate or better reisomerize the all-trans-retinal complex into the 11-cis-retinal complex. The reisomerization occurs in the RPE. Several reviews have described the chemistry and enzymology of the visual cycle in detail [23,24]. As shown in Figure 3, several genes in the visual cycle were stably expressed during development, including ABCA4, RBP1, RBP2, RBP5, RLBP1 (CRALBP), CRABP1, and RDH11. While other genes in this group showed minimal changes between 13- and 16-week fetuses, a significant decrease in expression was observed between the 16-week-old fetuses and adults, such as RPE65, LRAT, ALDH1A2, RDH5, and RDH14. The gene with the largest decrease in expression was RPE65 (60X), which participates in a rate-limiting step of the visual cycle. IRBP could not be detected in our microarray.
III. Phagocytosis-- Some of the molecular mechanisms responsible for the phagocytosis of photoreceptor outer segments by the RPE have been examined recently, and several genes implicated in the three stages of RPE phagocytosis (recognition, ingestion, and digestion) have been identified . All of the genes that have been reported to be associated with phagocytosis were present in 13-week-old fetuses. The majority of genes exhibited constant expression levels throughout development (Figure 4). The expression of MERTK and ITGAV (integrin αν) decreased during development; however, because the magnitude of the changes was not very great, these genes did not appear in our differentially expressed gene list. The only differentially expressed gene was MYO7A, which is responsible for the rapid relocation of phagosomes from the apical processes back to the base of the RPE cells , and its expression decreased during development.
IV. Adherens and tight junctions-- Adherens junctions were present before tight junctions and initial cell-to-cell contacts began to form. The cadherin and catenin families control the formation, maintenance, and function of adherens junctions . The expression of Cadherin-7, 10, 12, and 24 was significantly downregulated throughout development, while other cadherin family members exhibited stable expression (Figure 5). None of the catenin family members exhibited differential expression during development.
Tight junctions are composed of a growing list of proteins that can be grouped into transmembrane, adaptor, and effector proteins
 (Figure 6). Many of the typical transmembrane proteins, such as F11R, JAM2, JAM3, CRB3, occludin, and claudins 1–5, 7, 11, 12, 15,
and 19, were present in our microarray. Claudins play an important role in determining tight junction selectivity and permeability
, and they undergo constant remodeling during development . In our microarray, claudin-1, 2, 3, 5, and 19 exhibited
Genes coding for adaptor proteins such as ZO-1, ZO-2, ZO-3, TJAP1, MPDZ, MAGI1, MAGI3, ACTB, and CTNNA1 (α-catenin) were present in the fetuses and showed no changes in expression levels compared with those present in the adults. These adaptor proteins are important components of the tight junction complex structure for two reasons: first, they anchor the junction to the cytoskeleton, and second, they may help to regulate the overall function of the structure . The expression of MAGI2, PAR-3 (PARD3), and AF-6 (MLLT4) was significantly downregulated from 13 weeks of gestation to adulthood. Effector proteins include different types of signaling components, such as GTP-binding proteins, protein kinases, and phosphatases . Among these effector proteins, cell cycle regulators such as CDK4 and cyclin D1 showed significant downregulation during development, suggesting that the cell cycle rate was slowed. Another gene exhibiting higher expression in 13-week-old fetuses compared with adults was RAB3B, which is responsible for polarized vesicle transport.
V. Transcellular transport-- To satisfy the retina’s large requirement for glucose, the RPE expresses high levels of glucose transporters (Figure 7). The current study demonstrated that the housekeeping transporter GLUT3 (SLC2A3) was stably expressed, but that GLUT1 (SLC2A1) was downregulated during development. There was also stable expression of the mRNA of GLUT 4 (SLC2A4), GLUT 8 (SLC2A8), and GLUT 9 (SLC2A9), whereas the expression of GLUT 10 (SLC2A10) continuously decreased from the 13-week-old fetus to adulthood.
Na-K-ATPase, which is located in the apical membrane, provides the energy for transepithelial transport. Functional Na-K-ATPase is a heterodimeric protein consisting of an α-subunit and a β-subunit . Four isoforms of the α-subunit (α1, α2, α3, and α4) and three isoforms of the β-subunit (β1, β2, and β3) have been identified and were examined in this study, but only the mRNA expression of the β1-subunit decreased during development (Figure 8). The precise role of the β-subunit is still not known, but there is some evidence that the β-subunit may be important for cell-to-cell adhesion .
The RPE establishes a polarized distribution of ion pumps, channels, and transporters to mediate the active transcellular
transport of ions and organic solutes. Although many transporters have already been identified as highly enriched in the human
RPE , a genome-wide analysis has not yet been performed. Figure 9 shows some of the genes that are differentially expressed during development. The microarray data revealed that
A recent comprehensive study comparing the gene expression profiles of the fetal and adult RPE with somatic tissues identified 154 signature genes that are unique to the RPE . Our cross-reference analysis showed that 148 of the candidate signature genes found through their Affymetrix platform could be identified in our NimbleGen microarray. Among the 148 genes, 53 were downregulated, while only 5 were upregulated from the 13-week-old fetus to adulthood. This finding is consistent with the general gene expression pattern that the majority of expressed genes in RPE exhibited, i.e., a trend of downregulation during development. (Fifty-eight differentially expressed genes are listed in Table 2 and Table 3.) There were no differentially expressed genes present between 13 and 16 weeks of gestation. Gene ontology analysis indicated that 53 downregulated genes were involved in transport, eye development, cell cycle, metabolism, protein translation, signal transduction, visual perception, stress responses, inflammation, transcription regulation, pigment synthesis, cell adhesion and matrix biogenesis. The five upregulated genes were AHR, growth arrest specific 1 (GAS1), GRAMD3, LGALS8, and LSR. Gene ontology analysis revealed that these genes were associated with the cell cycle, signal transduction, embryonic development, and stress responses. It is worth mentioning that GAS1 is an important transcriptome for the negative regulation of the S phase mitotic cell cycle, and the disruption of GAS1 results in ectopic RPE proliferation and transdifferentiation in the retina . The expression of GAS1 increased in the RPE during development, and it was expected that, subsequent to the initial establishment of the RPE in the optic vesicles, proliferation in the presumptive RPE would cease, leading to the formation of a single layer of cuboidal cells. Our results further highlight the important role of upregulated GAS1 in normal RPE development. Among these 148 signature genes, we found that 90 exhibited minimal changes in expression that varied less than twofold. These genes included many common housekeeping genes, and they were not clustered into any of the expression patterns. Although these genes were less interesting from the standpoint of differential expression analysis, they are a key validation component of the data.
Seven genes that exhibited high degrees
Although the RPE is fundamentally important to retinal development and function and it is a critical focus of retinal degenerative
diseases and therapeutic interventions, relatively little is known about the molecular mechanisms of human RPE development.
This is the first study to demonstrate gene expression profiles in the human RPE during normal development. Approximately
4 weeks into human embryogenesis, the neural ectoderm grows out from the diencephalons to form the optic vesicle; it then
The RPE is an unusual epithelium, and it performs numerous functions that are essential for visual function. To better characterize the expression patterns of genes with critical functions within the RPE, we also performed a more focused analysis of the expressed genes grouped by function. All of the genes associated with pigment synthesis exhibited consistently decreasing expression between the fetuses and adults. This finding is in line with previous studies , which found that the majority of the melanin pigment in the RPE is synthesized prenatally within a few weeks of embryonic development, and it is stored in the melanosomes for the entire lifespan of the individual. Age, various types of light damage, and oxygen free radicals can contribute to the reduction of melanosomes in adults. It is known that melanin can exert a protective function in RPE cells in several ways. First, melanin shields these cells from the damage induced by sunlight and ultraviolet radiation . Second, melanin can counteract the oxidative stress caused in the RPE by free radicals derived from lipid peroxidation products  and by accumulated iron . Third, the high binding capacity of melanin to metal ions  and exogenous chemicals  also lends support to the protective role of melanin in the eye. A reduced capability to absorb light energy and increases in oxidative stress with aging are also important factors in the cascade of events leading to AMD, which may be the reason that AMD is more common in elderly people. The expression levels of genes involved in the visual cycle exhibited a downregulation trend. This finding is somewhat surprising because, with the formation of the photoreceptor outer segments, one would expect upregulation of visual cycle protein expression, as has been found in chickens and rodents [20,49]. By analyzing the expression levels of genes associated with phagocytosis, we found that there was also no expected upregulation trend with the formation of the outer segments and the lengthening of the RPE microvilli, but there was instead a slight decrease with development. To further validate the expression patterns of genes involved in the visual cycle and in phagocytosis in the human RPE, we compared our data with previous experimental results. Miller’s laboratory identified a set of genes that are unique to the RPE using a microarray technique . Among these genes, three (RPE65, RBP1, and RDH11) are involved in the visual cycle and four (ITGAV, VDP, LAMP2, and GULP1) are involved in phagocytosis. RPE65, RBP1, ITGAV, and GULP1 exhibited higher expression levels in the native fetal RPE compared with the native adult RPE. The expression levels of RDH11, VDP, and LAMP2 remained relatively unchanged between native fetuses and native adults. It was reported that the outer segments of the cones and rods are not fully developed until the 30th week of gestation . All of these results support the hypothesis that the expression of genes associated with the visual cycle and phagocytosis is not coordinated with the morphological maturation of the outer photoreceptor segments.
The apical junctional complex, composed of adherens and tight junctions, is an integrated signaling complex that helps to regulate cell shape, polarity, and proliferation . Most of the proteins found in all epithelial cells were presented in our microarray. These proteins form a complicated interacting network and are often regulated by posttranslational modifications. Therefore, it is difficult to correlate the expression levels of a handful of mRNAs with biologic processes. However, there were several genes of interest. Claudin-19 showed minimal changes in expression between 13 and 16 weeks of gestation but significantly decreased expression between 16-week-old fetuses and adults. As claudin-19 is the predominant claudin in the human RPE , changes in the expression level of claudin-19 may indicate changes in the selectivity and permeability of tight junctions. It is speculated that a relatively static extracellular environment is needed for embryogenesis because the development of the neural retina occurs slowly and the visual neurophysiology is inactive during embryogenesis. The higher expression of claudin-19 in the fetal RPE is thought to result in the formation of more “tight” junctions and a more powerful outer blood-retinal barrier, with greater selectivity and lower permeability. This hypothesis was supported by our previous study, in which we found that the mRNA of claudin-19 was more prominent than other claudins in cultured human RPE cells isolated from fetal eyes.
The TERs of hfRPE were approximately 1,000 and 400 Ω×cm2 in serum medium and serum-free medium, respectively [32,51], which is much higher than the value of 206 Ω×cm2 that was estimated for the adult human RPE in vivo . The permeability of the barrier increased 10 fold with the knockdown of claudin-19 by small interfering RNA . There are two reasons for the decreasing expression of claudin-19 in adults: on one hand, increased barrier permeability is needed with increasing activity of the visual neurophysiology and the increased requirements of nutrients and ions; on the other hand, proinflammatory cytokines increase junction permeability via decreases in the expression of claudin-19 . As another example, the TER in the ARPE-19 cell line, which was transformed from adult cells without the expression of claudin-19, was found to be approximately 50 Ω×cm2 . The expression of JAM-C, AF-6, and PAR-3 was consistently downregulated in the 13-week-old fetus. As described in the literature, JAM, AF-6, and the PAR-3/PAR-6/aPKC complex play major roles in the initial formation of adherens junctions and contribute to cell polarity , and the expression of JAM and AF-6 decreased substantially after the formation of adherens junctions. Thus, it is reasonable to believe that adherens junctions and cell polarity are established before 13 weeks of gestation in the human RPE.
The expression of many membrane pumps, channels, and transporters is up- or downregulated during human RPE development. A
few observations are notable
The results presented here demonstrate that, although the human RPE follows the same pattern of differentiation as that in other vertebrates at the macroscopic level, there are notable differences at the molecular level. In chicks and rodents, the expression of the visual cycle proteins is coordinated with the formation of the outer photoreceptor segments and displays a gradual polarization of the RPE, but the human RPE appears to achieve a high degree of maturity as early as 13 weeks of gestation (corresponding to the early phase of human retinal development). These results are supported by the observation that human RPE cells isolated from 16-week-old fetuses are already highly polarized, with a high TER, and the visual cycle is already well established . The different developmental pattern of the human RPE may be a mechanism to protect the photoreceptors from toxins during the long period of gestation.
Recently, the molecular signature of the human RPE was identified, and 154 RPE signature genes were validated by microarray and qRT2-PCR analysis, both in the RPE and in an independent set of 78 tissues . We cross-referenced our data with this reported data set for two reasons: first, we wanted to validate the accuracy of our microarray data; and second, this well defined data set enabled easier identification of important regulators during RPE development. Using their identified signature gene list, we found 58 genes showing twofold or more differences in expression at the three examined time points. These genes participate in a wide spectrum of activities, such as transport, eye development, the cell cycle, metabolism, signal transduction, visual perception, and pigment synthesis, suggesting that these 58 genes are the core developmental genes in the human RPE.
The present study not only offers information on the trends of changes in gene expression levels during human RPE development, but has also identified a set of genes that are crucial for developmental regulation. These results may be helpful in furthering our understanding of the developmental processes of the human RPE.
The authors are grateful to all donors and their families. This work was supported by National Natural Science Foundation of China (Grant No.30772381).