Transcriptome of the human retina, retinal pigmented epithelium and choroid

YGENO-08716; No. of pages: 12; 4C: Genomics xxx (2015) xxx–xxx

Contents lists available at ScienceDirect

Genomics journal homepage: www.elsevier.com/locate/ygeno

a

1 0

a r t i c l e

11 12 13 14

Article history: Received 28 August 2014 Accepted 24 January 2015 Available online xxxx

15 16 17 32 18 33 19 20 21 22 23 24

Keywords: Age-related macular degeneration cDNA microarray Glaucoma Retina Retinal pigmented epithelium Choroid RNA-Seq Serial analysis of gene expression Transcriptome

Center for Applied Genomics, The Children's Hospital of Philadelphia, Philadelphia, Pa 19104, USA Department of Ophthalmology, University of Pennsylvania, Philadelphia, Pa 19104, USA Department of Biostatistics and Epidemiology, University of Pennsylvania, Philadelphia, Pa 19104, USA d Department of Ophthalmology, University of Alabama School of Medicine, Birmingham, Al 35294, USA c

a b s t r a c t

F

The retina and its adjacent supporting tissues – retinal pigmented epithelium (RPE) and choroid – are critical structures in human eyes required for normal visual perception. Abnormal changes in these layers have been implicated in diseases such as age-related macular degeneration and glaucoma. With the advent of highthroughput methods, such as serial analysis of gene expression, cDNA microarray, and RNA sequencing, there is unprecedented opportunity to facilitate our understanding of the normal retina, RPE, and choroid. This information can be used to identify dysfunction in age-related macular degeneration and glaucoma. In this review, we describe the current status in our understanding of these transcriptomes through the use of highthroughput techniques. © 2015 Published by Elsevier Inc.

4.

5.

R

Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Structure of the human eye: retina and choroid . . . . . . . . . . . . . . Transcriptome analysis of the retina and RPE/choroid . . . . . . . . . . . 3.1. High-throughput technologies: cDNA microarray, SAGE and RNA-Seq 3.2. RNA qualities and methods of case ascertainment . . . . . . . . . Exploring the human retina transcriptome . . . . . . . . . . . . . . . . 4.1. Number of genes detected in retina . . . . . . . . . . . . . . . . 4.2. Effect of aging on gene expression . . . . . . . . . . . . . . . . 4.3. Effect of retina anatomic location on gene expression . . . . . . . . 4.4. RNA-Seq assesses alternative splicing . . . . . . . . . . . . . . . 4.5. Comparison between RNA-Seq and Microarray . . . . . . . . . . . Exploring the human RPE/choroid transcriptome . . . . . . . . . . . . . 5.1. Number of genes detectable by microarray . . . . . . . . . . . . . 5.2. Effect of anatomic location on gene expression . . . . . . . . . . .

N C O

1. 2. 3.

R

Contents

U

40 41 42 43 44 45 46 47 48 49 50 51 52 53

P

i n f o

37 35 34 36 39 38

R O

b

O

6 7 8 9

D

5

Lifeng Tian a, Krista L. Kazmierkiewicz b, Anita S. Bowman b, Mingyao Li c, Christine A. Curcio d, Dwight E. Stambolian b,⁎

E

4Q6

T

3

Transcriptome of the human retina, retinal pigmented epithelium and choroid

C

2Q5

Review

E

1

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

. . . . . . . . . . . . . .

0 0 0 0 0 0 0 0 0 0 0 0 0 0

Abbreviations: AMD, age-related macular degeneration; eQTL, expression quantitative trait loci; GWAS, genome-wide association studies; FPKM, Fragments Per Kilobase of gene per Million mapped fragment; MAQC, MicroArray Quality Control; MR, macular retina; PLIER, Probe Logarithmic Intensity Error; PR, peripheral retina; RPE, retinal pigment epithelium; RNASeq, RNA sequencing; SAGE, serial analysis of gene expression. ⁎ Corresponding author at: 313B Stellar-Chance Labs, 422 Curie Blvd., Philadelphia, PA 19104, USA. Fax: +1 215 573 6728. E-mail addresses: [email protected] (L. Tian), [email protected] (K.L. Kazmierkiewicz), [email protected] (A.S. Bowman), [email protected] (M. Li), [email protected] (C.A. Curcio), [email protected] (D.E. Stambolian).

http://dx.doi.org/10.1016/j.ygeno.2015.01.008 0888-7543/© 2015 Published by Elsevier Inc.

Please cite this article as: L. Tian, et al., Genomics (2015), http://dx.doi.org/10.1016/j.ygeno.2015.01.008

25 26 27 28 29 30 31

2

54 55 56 57 58 59 60 61 62 63 64 65

L. Tian et al. / Genomics xxx (2015) xxx–xxx

6.

Transcriptome differences between human retina and RPE/choroid . . . . 6.1. Comparison of differential expression . . . . . . . . . . . . . 6.2. Comparison of pathway analysis . . . . . . . . . . . . . . . 7. The transcriptome in AMD and glaucoma . . . . . . . . . . . . . . . 7.1. Age-related macular degeneration . . . . . . . . . . . . . . . 7.2. Glaucoma . . . . . . . . . . . . . . . . . . . . . . . . . . 8. Bioinformatic databases for human retina and RPE/choroid transcriptomes 9. Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . Competing interests . . . . . . . . . . . . . . . . . . . . . . . . . . . Authors' contributions . . . . . . . . . . . . . . . . . . . . . . . . . . Acknowledgments . . . . . . . . . . . . . . . . . . . . . . . . . . . . References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

. . . . . . . . . . . .

0 0 0 0 0 0 0 0 0 0 0 0

F

66

109 110

The retina is a specialized neural tissue lining the back of the eye that is responsible for vision. The retina originates as an outgrowth of the brain during ontogenesis and is thus part of the central nervous system. There are 21 chorioretinal layers and clinically important spaces, many of which can be visualized non-invasively at high-resolution (Fig. 1) [7]. The macula has the highest overall density of neurons and is responsible for acute central vision (diameter = 6 mm = 21° of visual angle). The

90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106

111 112 113 114 115

C

88 89

E

86 87

R

84 85

R

82 83

O

80 81

C

78 79

N

76 77

U

74 75

3 . Transcriptome analysis of the retina and RPE/choroid

150

3.1. High-throughput technologies: cDNA microarray, SAGE and RNA-Seq

151

A cDNA microarray consists of immobilized probes complementary to known transcripts on a solid substrate [12]. Isolated RNA is labeled with fluorescent dyes and hybridized to the cDNA microarrays, washed, and scanned with a laser scanner. The amount of fluorescent dye intensity is a measure of gene expression. In early versions of cDNA microarrays, biases and artifacts produced inconsistent results among the same samples. In 2006, quality control standards were developed by the MicroArray Quality Control (MAQC) to address these issues [13]. Still, microarrays are unable to identify RNA editing events or novel isoforms, and they cannot accurately measure absolute expression levels due to hybridization and background variation [14]. Serial analysis of gene expression (SAGE) can perform a quantitative analysis of transcripts without requiring prior knowledge of the

152 153

O

2 . Structure of the human eye: retina and choroid

72 73

R O

108

70 71

P

107

In the past decade, multiple common genetic variants have been identified for prevalent diseases through genome-wide association studies (GWAS). Translating the results of GWAS to new treatments for disease still faces several challenges including: 1) Single nucleotide polymorphisms (SNPs) significant in GWAS are not causative and may be in linkage disequilibrium with causative SNPs that are distantly located; 2) associated genes harboring these SNPs have not been evaluated in tissues affected by disease; and 3) gene expression in cells central to the disease have not been evaluated for regulatory effects of these SNPs. While accumulating scientific evidence suggests that regulatory changes contribute to the phenotypes interrogated by GWAS, progress in defining these regulatory changes has been slow. This problem is apparent in the eye, a highly specialized structure composed of multiple layers and cells derived from mesoderm and neuroectoderm that has produced multiple GWAS hits outside coding regions for age-related macular degeneration (AMD), glaucoma, and refractive error, three widely prevalent eye conditions. Here, we focus on AMD and glaucoma, the most prevalent causes of untreatable vision loss in the elderly. The transcriptome is a complex mixture of multiple isoforms for known genes and non-coding RNAs that will require multiple approaches for characterization. These multiple transcript types are heritable, based on evidence that gene expression is heritable between monozygotic twins [1], Centre d'Etude du Polymorphisme Humain families [2,3], and Icelandic families [4]. The transcriptome is also cell/ tissue specific, requiring cell specific characterization. For example, there is lack of significant sharing of cis-expression quantitative trait loci (eQTLs) between adipose tissue and blood [4] or between brain cortical tissue and blood mononuclear cells [5], making it necessary to characterize the transcriptome of each cell type. Vision scientists have long recognized the importance of characterizing the transcriptome through their use of microarrays to evaluate expression. This technique persists as a mainstay among vision scientists even though microarrays have been shown to have biases in hybridization strength, as well as the potential for cross hybridization to probes with similar sequences [6]. Herein, we provide a brief overview of the eye structure with an emphasis on its posterior layers, and follow with a review of SAGE, cDNA microarray, and RNA sequencing (RNA-Seq) studies of the protein-coding transcriptomes of human retina, retinal pigmented epithelium (RPE) and choroid.

116 117

D

68 69

nasal retina (close to the nose, seeing the temporal visual field) has 1.4–3-fold more neurons, depending on the cell type, than temporal retina (close to the temple, seeing the nasal visual field.) The cellular composition of the neurosensory retina is highly organized with a variety of neurons and two glial types. Major neuron types have been historically classified based on morphology, laminar distribution, topography, brain connections, physiological responses to light stimuli, transmitter pharmacology, and molecular signatures, such as immunoreactivity to specific antibodies. These cells likely differ by their RNA transcriptional programs as well. A recent classification of retinal cell types indicates three cone photoreceptors, one rod photoreceptor, two horizontal cells, 13 bipolar cells, ≥29 amacrine cells, and ≥20 ganglion cells. Only the ganglion cells project to the brain (via the optic nerve). Ganglion cells are preferentially affected in glaucoma, a prevalent cause of vision loss for which high intraocular pressure is an important risk factor. Glial cells important for retinal function are Müller cells that span all the neuronal layers and astrocytes among ganglion cell axons en route to the optic nerve head. The RPE and choroidal vasculature constitute the photoreceptor support system, which is affected in AMD, a major cause of vision loss in the elderly of European descent worldwide. Polarized RPE has demanding dual roles serving photoreceptors apically and choroid basolaterally. Distinctive extracellular lesions that differentially confer risk for AMD progression distributes on both aspects of this key cell layer [8–10]. The retina has two vascular beds with different properties and propensity for disease. The retinal circulation serves the inner retina and is within the blood–retina barrier. The choroid serves the photoreceptors and RPE, and it is part of the systemic circulation. The RPE maintains the outer limit of the blood–retina barrier with junctional complexes. The choroid is distinguished by the highest blood flow in the body, especially under the macula, and it thins markedly with aging [11]. Choroidal cells include vascular and lymphatic endothelia, smooth muscle cells, fibroblasts, melanocytes, mast cells, autonomic neuronal ganglia, and resident and transient cells of monocyte lineage.

E

1 . Introduction

T

67

Please cite this article as: L. Tian, et al., Genomics (2015), http://dx.doi.org/10.1016/j.ygeno.2015.01.008

118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149

154 155 156 157 158 159 160 161 162 163 164

3

R O

O

F

L. Tian et al. / Genomics xxx (2015) xxx–xxx

173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194

E

R

171 172

R

169 170

N C O

167 168

transcript sequence. Compared to microarrays, SAGE represents an unbiased, comprehensive representation of the transcriptome [15,16]. A short tag i.e., 9–14 bases of a gene transcript is linked with 20–50 other such tags in a cloned DNA fragment. Sequencing of such a clone provides the sequence of 20–50 tags, and a set of a few thousand such clones represents a library of SAGE tags [17]. Most individual SAGE tags can be assigned to specific genes by alignment. It can identify low abundance transcripts and detect relatively small differences in their expression. However, it is limited to simultaneously analyzing a small number of expression profiles. RNA-Seq can measure absolute gene expression levels and identify the exact sequence of transcripts belonging to mRNA and non-coding RNA [6]. RNA-Seq can provide a measure of gene expression as well as the underlying genetic variants that influence gene expression. This additional information provided by RNA-Seq compared to microarrays and SAGE has made RNA-Seq the current method of choice for studying the transcriptome. The major advantages of RNA-Seq include:

U

165 166

C

T

E

D

P

Fig. 1. Chorioretinal tissue layers in human eye. A. Cross-section of a human eye. Green bracket shows extent of the posterior pole that is illustrated in panel B, including optic nerve head and macula with fovea. Light shines into this eye from the right. Schematic from David Fisher designs. B. In vivo high-resolution imaging of a normal human retina (black bar) and choroid (white bar) shows tissue layers of the 6 mm diameter macula and adjacent optic nerve clearly. Light shines into this eye from above. The neurosensory retina contains multiple bands of alternating high and low reflectivity that coincide in part with the anatomical layers. Blue bar delimits layers occupied by photoreceptors. Arrowheads indicate a dip at the fovea, shown by histology in panel C. The choroidal vasculature contains lumens of large vessels. It is bounded externally by the sclera, which has a more homogeneous reflectivity. Swept-source optical coherence tomography scan, a 58 year-old male, courtesy of R.F. Spaide, MD. C. Ex vivo high-resolution histology of the fovea, which is located in the center of macula and is responsible for high acuity vision. The neurosensory retina including thickened photoreceptor layers and choroid are indicated by black, blue, and white arrows, respectively. A foveal pit is created by inner retinal neurons, Müller cells, and accompanying retinal vasculature being swept to the side of the visual axis. This facilitates light capture by tightly packed inner segments of cone photoreceptors on the opposite aspect of the retina. The photoreceptors and Müller cells (glial cells) have long fibers that together make a thick layer that is 14% of total retinal thickness in this region (Henle fiber layer). The RPE is a simple cuboidal epithelium that sits on Bruch's membrane. Bruch's membrane is the inner wall of the choroid and functions as a vessel wall between choriocapillaries and RPE, as well as substrate for RPE attachment. RPE maintains health of the photoreceptors and the choroid. RPE is visible clinically by its melanosomes and lipofuscin, which is autofluorescent due to long-lasting vitamin A by-products. The choroidal vasculature contains fenestrated capillaries adjacent to the RPE and larger arteries and veins near the sclera. A 63-year-old female recovered b3 h after death, osmium post-fixation, epoxy embedding, 1 μm thick section, and toluidine blue stain. Figure was prepared by J.D. Messinger, DC. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

• Identification of known and previously unknown transcripts as opposed to array-based methodologies that depend on annotation of known transcripts. • Accurate quantification of transcript levels by counting reads that map to transcripts as opposed to inferring levels from hybridization intensity signals as in microarrays. Background noise levels are essentially zero in RNA-Seq, making it easier to assay rare transcripts provided sequencing depth is sufficient. • The counting methodology used to ascertain expression levels in RNA-Seq provides a more accurate assessment of expression levels throughout the dynamic range than can be obtained through arraybased methods.

• Integration of RNA-Seq and whole-genome DNA sequencing data will 195 allow determination of RNA-DNA differences and allele-specific 196 expression within individual samples. 197 198

However, the cost of RNA-Seq, compared with microarray, is significantly higher, especially when deep depth of sequencing coverage is desired. Another challenge to the downstream data analysis of RNASeq is the unevenness of the depth of coverage across the expressed regions in the genome [14].

199

3.2. RNA qualities and methods of case ascertainment

204

For studies utilizing human eye tissues obtained from eye banks, the interval between donor death and RNA extraction for subsequent analysis should be considered when assessing the quality of published data. This interval can vary widely between studies. RNA deteriorates by deadenylation of the 3′ poly-A tail, decapping of the 5′ end, and degradation by 5′ and 3′ exonucleases, with a major role being played by the 5′ exonuclease. An analysis of RNA stability through a time course simulation of tissue recovery using eye bank conditions was performed by Malik et al. using pig eyes [18], which are of similar size to human eyes. At different time intervals after animal sacrifice, the authors isolated retina and RPE from the eyes and assessed the recovered RNA. They determined that RNA for housekeeping and RPE-specific genes deteriorated starting at 5 h. Immersion into a cold RNA stabilization medium (RNAlater) delayed degradation onset at least 24 h. Retina, with many more cells than RPE (N100 M vs. 4–6 M), had more stable RNA than RPE. However, RNA specific to potentially vulnerable minority cell populations was not investigated in this study. The possible reason for more rapid degradation of RNA from the RPE, a cell central to AMD,

205

Please cite this article as: L. Tian, et al., Genomics (2015), http://dx.doi.org/10.1016/j.ygeno.2015.01.008

200 201 202 203

206 207 208 209 210 211 212 213 214 215 216 217 218 Q8 219 220 221 222

236 237 238 239 240 241 242 243 244 245

E

234 235

Author

Samples

t1:4

Yoshida et al. [27]

Five human retinas

Hornan et al. [31]

• 2 young (ages 13–14) • 3 elderly (ages 62–72) • 24+ h to processing Human neural retinas (ages 40–77)

Cai et al. [28]

• 12 pooled fovea macula • 2 pooled macula • 10 peripheral • Time to processing not given Twelve human retinas

t1:6

R

t1:5

R

t1:3

Platform

Results

Microarray

17 (1.4%) genes were expressed at higher levels in retinas from young, whereas 7 (0.6%) genes were designated as elderly-dominant. Northern blot analysis and qRT-PCR results confirmed the changes in expression in 8 of 10 genes examined.

T

232 233

C

230 231

E

229

F

Table 1 Studies of retinal transcriptome.

227 228

O

t1:1 t1:2

225 226

clinical records are the preferred method of case ascertainment, as optic nerve damage may not be visible by macroscopic inspection until advanced stages (C.A. Girkin, MD, MSPH, personal communication, 5/9/14). In Tables 1 and 2, we summarize the studies of retina and RPE/choroid transcriptome reviewed in this manuscript, including sample size and post-mortem delay to preservation. Finally, studies should be scrutinized for their dissection techniques, and particularly how they name samples of different regions. Tissue punches covering differing extents of retina can yield variable results, especially in the macula, which contains the singularity of densely packed cone photoreceptors. The diameter of punches and their point of centration should be specified with dimensions (e.g., a punch 8 mm in diameter centered at 13 mm temporal to the fovea), as well as names (“periphery”), because names are variably defined. “Foveal” punches intended to capture only or primarily cone photoreceptors should be ≤1 mm in diameter. An annulus remaining after removal of the fovea (as just defined) from the 6-mm diameter macula is most accurately described as perifovea. The perifovea has a 9:1 rod:cone ratio in a young adult and 6:1 in an older adult. Outside the perifovea, the rod:cone ratio ranges from 30:1 just outside the macula to 10:1 at the retinal edge. Photoreceptors are the majority cell population and dictate the organization and cellular content of both the RPE/choroid and inner retinal layers. Therefore, rod:cone ratios should be factored into experimental design.

D

246

could be due to the high metabolism of native RPE and mechanical damage during dissection of this single cell layer [18]. In our studies, we examined both RNA integrity number at the time of processing after overnight shipping in cold stabilization medium and assessed expression of genes specific for numerically sparse cell types with good success [19]. In addition to RNA quality, studies utilizing diseased donor eyes should be checked for methods of case ascertainment, typically done with family report, clinical records, post-mortem inspection and/or histology of the fundus, or a combination. Clinical records of eye donors can vary considerably in source and in length of time before donor death. Pathology grading systems exist for donor eyes based on the color fundus photography systems used for large epidemiologic and clinical studies [7,20–22]. Through recent clinical experience with optical coherence tomography (OCT, Fig. 1B), a high-resolution imaging technology that provides detailed cross-sectional views of retina, it is now known that important AMD pathology is not captured well by color photography. Even less detail is apparent with color photography of partly opacified post-mortem retina. Identifying and quantifying these newly discovered AMD-relevant endophenotypes (e.g., subretinal drusenoid deposits [9] and outer retinal tubulation [23]) will be important for quantitative assessment of gene expression in tissues. Ex vivo OCT of short post-mortem eyes, if validated by histology, holds great promise for the molecular pathology laboratory [7,24,25]. For glaucoma,

R O

223 224

L. Tian et al. / Genomics xxx (2015) xxx–xxx

P

4

Microarray

qRT-PCR results of eight genes enriched in fovea macula and four genes enriched in peripheral retina

Microarray

The transcriptome of the human retina is affected by age and topographic location.

O

• 6 young (ages 18–43) ○ 6 macula vs. 6 peripheral

• 48 h to processing 10 tissues from 6 human eyes

Microarray

Ocular Tissue Database (OTDB) for normal eye gene expression is provided.

Young et al. [30]

• stored in RNAlater, time to processing not given Fifteen human eyes

Microarray

Significant gene expression fold changes (N1.5) were found in adult versus fetal retina.

Sharon et al. [17]

• 9 human fetal eyes (four male and five female) • 6 adult normal eyes (three male and three female) • b6 h to processing, stored in RNAlater, shipped overnight 2 human retinas (ages 88 and 44)

SAGE

An average of 13,779 UniGene transcripts was detected in retina.

RNA-Seq

75% of known exons were detected. Novel exons increase the number of exons by 3%.

U

t1:8

Wagner et al. [35]

N

t1:7

C

• 6 elderly (ages 72–79) ○ 6 macula vs. 6 peripheral

t1:9

• time to processing not given 3 normal human retinas (ages 42, 44, 46, one male, two female)

t1:10

Farkas et al. [29]

t1:11

Li et al. [19]

• time to processing not given 8 normal human retinas

RNA-Seq

80% of the transcriptome is expressed in retina. There is significant differential expression between retinal layers and locations of a tissue layer.

Whitmore et al. [32]

• b6 h to processing, stored in RNAlater, shipped overnight, processed upon arrival 4 normal human retinas

RNA-Seq

No genes were differentially expressed between temporal and nasal retina.

t1:12

• b6 h to processing, stored in liquid nitrogen • nasal, temporal, and macular regions

Please cite this article as: L. Tian, et al., Genomics (2015), http://dx.doi.org/10.1016/j.ygeno.2015.01.008

247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270

L. Tian et al. / Genomics xxx (2015) xxx–xxx t2:1 t2:2

5

Table 2 Studies of RPE/choroid transcriptome.

t2:3

Author

Samples

Platform

Results

t2:4

Bowes Rickman et al. [42]

17 human eyes

SAGE

Radeke et al. [39]

• ages 28–66 • b6 h to processing, stored in RNAlater 15 human eyes

Three differential expression lists were reported: macula retina enriched, periphery retina enriched, and RPE enriched.

Microarray

76 genes were differentially expressed between macular and peripheral regions. 29 genes selected for validation via qPCR and 21 were confirmed.

Van Soest et al. [40]

• from each eye, one macular and one peripheral punch of the RPE/choroid complex was taken • 0.5 to 7.4 h to processing, stored in RNAlater, stored up to a week 6 human eyes

Microarray

438 genes were differentially expressed between macular and peripheral regions. Validation of 33 genes using RT-PCR resulted in an 84% correlation.

Booij et al. [36]

• ages 17–36 • 2 females, 4 males • from each eye, one macular and two peripheral fragments were isolated • 24+ h to processing 6 human eyes

Microarray

Attempted to identify the expression profile of macular RPE and demonstrated the lack of knowledge regarding the RPE transcriptome. Functional analysis revealed enrichment for pathways such as oxidative phosphorylation and ATP synthesis.

Booij et al. [37]

• ages 63–78 • macular RPE isolated • 30+ h to processing 5 human eyes

Microarray

Strunnikova et al. [38]

• ages 63–78 • isolated RPE, choroid, and photoreceptor layers • 16–22 h then frozen 4 adult Caucasian (ages 64–89)

Microarray

154 RPE signature genes were identified. Using qRT-PCR, 48 genes were highly expressed both in vivo and in vitro. Demonstrated that culturing cells can change the expression, but did not affect enrichment of signature genes

Li et al. [19]

• 4 fetal native RPE samples • 4 fetal native choroid samples • b24 h to processing • fetal RPE primary cultures • ARPE-19 cell line 8 normal human RPE/choroid/sclera

RNA-Seq

926 genes were differentially expressed between macular and peripheral RPE/choroid/sclera.

Whitmore et al. [32]

• b6 h to processing, stored in RNAlater, shipped overnight, processed upon arrival 4 normal human RPE/choroid

RNA-Seq

81 genes showed increased expression in the nasal RPE/choroid and 39 genes showed decreased expression.

t2:5 Q1

• normal eyes (ages 45–65 and 75–87) ○ 9 males, 1 female

t2:10

t2:11

R

R

• b6 h to processing, stored in liquid nitrogen • nasal, temporal, and macular regions

4 . Exploring the human retina transcriptome

272

4.1. Number of genes detected in retina

273 274

There are over 60 functionally different neuronal cell types present in the retina. It is estimated that a complete transcriptome of the mammalian retina might consist of approximately 25,000 transcripts [26]. Early attempts to quantitate the retina transcriptome were done with focused arrays. Yoshida et al. used an array containing 2400 probes for human neuronal genes to assay RNA from two young and three elderly post-mortem retinas [27]. Fifty-two percent of the 2400 gene probes reportedly hybridized strongly with the retinal cDNA. Cai et al. assessed differential gene expression of young maculas and old maculas, and young and old peripheral retinas from 12 healthy donor eyes on an Affymetrix Human Genome U133 plus 2.0 array, containing 54,600 gene probes [28]. On average, 49% of the probes were detected in human retinas. There was no significant difference in the number of genes detected between human young macula (26,686 ± 319), old macula (26,956 ± 275), young peripheral retina (27,122 ± 108), or old peripheral retina (26,533 ± 490). Sharon et al. used the SAGE technique to assess the genes expressed in the human peripheral and macular retina from two individuals [17]. They detected an average of 13,779 UniGene transcripts from two peripheral retinas and one macular retina.

279 280 281 282 283 284 285 286 287 288 289 290 291 292

U

277 278

N C O

271

275 276

D

E

t2:9 Q3

T

O

R O

P

t2:8

C

t2:7

E

t2:6 Q2

F

• AMD eyes (ages 78–92) ○ 3 males, 2 females

Identified 114 RPE-specific genes. Selected 39 for validation. 85% were validated through literature and s-QPCR confirmed RPE expression for the remaining genes.

Recently, RNA-Seq has revealed that the retinal transcriptome is more complex than previously reported. Farkas et al. reported RNA-Seq results from three normal human retinas [29]. They were able to detect 75% to 82% of all exons annotated in eight annotation reference tracks (UCSC, RefSeq, CCDS, Vega, Ensembl, Aceview, Gencode, and LincRNAs) at an average read depth of five or greater. Li et al. performed a more comprehensive analysis of eight normal human retinas [19]. They estimated the expression levels of 23,569 RefSeq protein-coding genes using the Fragments Per Kilobase of gene per Million mapped fragments (FPKM) metric. With coverage depth ranging from 66 to 133 million paired-end reads per sample, they detected the expression (i.e., FPKM N 0) of the majority of known protein-coding genes. The average number of expressed RefSeq genes was 18,318 (78%) in macular retina (MR), and 18,850 (80%) in peripheral retina (PR). This is a significant increase over the SAGE results reported by Sharon et al. [17]. The reported RNA-Seq results reflect the complexity of the retinal transcriptome that was not appreciated with cDNA microarrays and SAGE.

293 294

4.2. Effect of aging on gene expression

311

295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310

Yoshida et al. identified 17 (1.4%) genes expressed at higher levels in 312 young retinas and 7 (0.6%) genes expressed at higher levels in older 313 retinas [27]. Northern blot analysis and qRT-PCR confirmed the changes 314

Please cite this article as: L. Tian, et al., Genomics (2015), http://dx.doi.org/10.1016/j.ygeno.2015.01.008

369

Although Cai et al. found no significant difference in expression between regions of eyes of different age groups, 105 out of 26,700 gene probes were differentially expressed between macula and periphery if all age groups are combined [28]. They found 24 gene probes expressed at higher levels in the peripheral than the macular retina. A second study by Hornan et al. compared the expression profile of macular vs. peripheral retina using 2–4 mm retinal punches from ten peripheral and two macular retinas and found 8 genes enriched in macula and 4 genes enriched in peripheral retina [31]. The eight fovea maculaenriched genes are: OPN1LW, HDAC9, NPIP, SRGAP3, SRGAP2, FLJ, CRP, and BHM2. Cai et al. also reported two genes, HDAC9 and SRGAP2, to be enriched in the macula [28]. Hornan et al. reported 4 peripheral retina enriched genes — RP2, YAP, RHO, and SAG [31]. One gene, RHO, was also enriched in peripheral retina by Cai et al. [28]. In a recent transcriptomic analysis of human neural retina and RPE/choroid by RNA-Seq, Whitmore et al. reported that 1) no genes were differentially expressed between temporal and nasal retina; 2) 30 genes showed increased expression in nasal retina and 128 genes showed decreased expression when comparing nasal retina vs. macular retina; and 3) 37 genes showed increased expression in temporal retina and 323 genes showed decreased expression when comparing temporal retina vs. macular retina [32]. Li et al. showed that 1) 2051 genes were differentially expressed between macular and peripheral retina; and 2) 60–70% of them had higher expression levels in macula than in periphery [19].

370

4.4. RNA-Seq assesses alternative splicing

371

Farkas et al. performed RNA-Seq from three normal human retinas and detected 79,915 novel alternative splicing events, including 29,887 novel exons, 21,757 3′ and 5′ alternative splice sites, and 28,271 exon skipping events [29]. They also identified 116 potential novel genes in the human retina.

340 341

346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368

372 373 374 375

379 380 381 382 383 384 385 386 387 388 389 390 391

5. Exploring the human RPE/choroid transcriptome

392

The RPE is a simple cuboidal epithelium apposed to the neurosensory retina (Fig. 1C). It is best understood for its function in providing vitamin A derivatives for phototransduction and the daily phagocytosis of spent photoreceptor outer segment tips. However, the RPE also transfers from the systemic vasculature all the nutrients required by photoreceptors, acts as a waste-management system, maintains the choriocapillaris vascular endothelium and Bruch's membrane, regulates the outer retinal component of the blood–retina barrier, and orchestrates the secretion of cytokines, among many other functions. The RPE is involved in many macular diseases especially AMD.

393 394

T

338 339

C

336 337

E

334 335

R

332 333

R

330 331

O

328 329

C

326 327

N

324 325

U

322 323

377 378

F

345

321

Previous studies have shown RNA-Seq to be superior at quantifying transcript levels through read counts in comparison to microarrays that infer levels from hybridization intensity signals [6,33]. Background noise levels are essentially zero in RNA-Seq, which enables one to assay rare transcripts provided sequencing depth is sufficient [34]. The counting method of RNA-Seq also provides a more accurate assessment of expression levels throughout the dynamic range compared to arraybased methods. A comparison of microarray data from Wagner et al. [35] with RNA-Seq data by Li et al. [19] confirms that RNA-Seq is more sensitive for low abundance transcripts. Li et al. found 9864 genes with a FPKM value lower than 1 [19]; whereas, only 94 genes having a Plier number lower than 1 were found using microarray [35] (Fig. 2). The difference in sensitivity between microarray and RNA-Seq for rare transcripts makes RNA-Seq the preferred method for transcriptome evaluation.

O

4.3. Effect of retina anatomic location on gene expression

319 320

376

R O

344

317 318

4.5. Comparison between RNA-Seq and Microarray

E

342 343

in expression in 8 of 10 genes examined and included increases in IFN-responsive transcription factor subunit (ISGF3G), creatine kinase B (CKB), and pancreatic amylase (AMY2A), and a decrease in TGF-beta receptor interacting protein 1 (TRIP1), LPS-induced TNF-alpha factor (PIG7), alpha-1 (E)-catenin (CTNNA1), ubiquitin hydrolase (USP9X), GABA receptor beta-3 subunit (GABRB3), and alpha-1 Type VII collagen (COL7A1). Cai et al. also reported changes in the retina transcriptome with age [28]. Macular and peripheral retina RNA from six young and six older human donor eyes were analyzed by microarray. They found 85 genes expressed at higher levels in young macula and 55 genes expressed at higher levels in older maculas. Fifty-two genes were more highly expressed in young peripheral retina and 34 genes more highly expressed in older peripheral retina. Although there was evidence of differential expression in this study, none of the differences were statistically significant (F-test, p N 0.05). The studies of Cai et al. [28] and Yoshida et al. [27] have little overlap for differentially expressed genes. Even for genes identified in both, the directions of enrichment might be different; for example, ENO3 was differentially expressed in both studies, but Yoshida et al. [27] reported enrichment in young retina while Cai et al. [28] found the gene enriched in the older macula. Another study by Young et al. compared gene expression patterns in normal human fetal (N = 9) and adult retinas (N = 6) using Illumina HumanHT-12 v4 Expression BeadChips [30]. A significant number of differentially expressed genes (N = 1185, fold change ≥ 1.5) were found between adult and fetal retina. These differentially expressed genes were involved in development, cell death/growth, cancer functions, and signaling pathways.

P

315 316

L. Tian et al. / Genomics xxx (2015) xxx–xxx

D

6

395 396 397 398 399 400 401 402

5.1. Number of genes detectable by microarray

403

Booij et al. used microarray to assess the expression profile of the RPE from six eyes of normal donors between the ages of 63 and 78 [36]. Using six 22 k custom arrays on laser dissected macular RPE cells they found that the highly expressed genes in the RPE represented functional categories in oxidative phosphorylation, ribosome pathway, and ATP synthesis. A subsequent study, carried out by the same group, assessed the purity of the RPE samples through microarray analysis of RPE, adjacent photoreceptors, and choroid cell layers [37]. The justification for this design is based on the difficulty of isolating RPE cells free of contamination from the retina and/or choroid. Genes with 2.5 × higher expression in the RPE relative to either the retina or choroid was limited to 114 RPE-specific genes. Those genes with high expression levels in two or more tissues were not analyzed further. Thirty-nine of the top highly expressed genes in the RPE were selected for validation and confirmed by a literature review and semi-quantitative PCR. Pathway analysis of the 114 genes revealed an overrepresentation of genes in categories such as membrane transport, vision, and glycoproteins. Three canonical pathways RARactivation, retinol metabolism, and GABA receptor signaling were identified. In 2010, Strunnikova et al. reported a microarray study of adult RPE, fetal RPE, and the RPE-derived cell line ARPE-19 [38]. A gene was designated as RPE-specific when its RPE cell expression was 10-fold higher than its mean relative expression in 18 non-ocular cell lines. One hundred fifty-four genes fulfilled this criterion and were labeled as signature RPE genes. Subsequent validation via qRT-PCR of these 154 genes confirmed 100 genes were RPE signature genes. This standard chosen by these authors is good for capturing large effects. However, it is not complete, as low abundance genes known to be expressed in RPE by other methods were not included in the set.

404 405

Please cite this article as: L. Tian, et al., Genomics (2015), http://dx.doi.org/10.1016/j.ygeno.2015.01.008

406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434

7

R O

O

F

L. Tian et al. / Genomics xxx (2015) xxx–xxx

Fig. 2. Distributions of gene expression values from microarray (Plier) and RNA-Seq (FPKM). X axis shows the log2 value of the Plier number; Y axis shows a histogram of the number of genes.

449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476

C

447 448

E

445 446

R

443 444

R

441 442

N C O

439 440

U

437 438

P

To uncover expression differences present in the macula that confer vulnerability to macular diseases, Radeke et al. assessed gene expression for the macula and extramacular regions of RPE/choroid from middle-aged (45–65 years old) normal eyes, older (75–87 years old) normal eyes, and AMD afflicted eyes [39]. Using the Agilent Human 1A microarray platform for the AMD samples and the Agilent Whole Human Genome array for the middle-aged and older normal samples, differentially expressed genes were determined between macular and extramacular regions of the RPE/choroid complex. Of the 201 genes identified as differentially expressed in the RPE/choroid, 76 genes remained as differentially expressed after retinal contaminants were removed. Twenty-nine genes were chosen for validation by qRT-PCR, and 21 of the 29 genes were confirmed. The differentially expressed genes grouped into categories including inflammatory response or regulation, angiogenesis, and those associated with extracellular matrix (ECM). A second study by van Soest et al. that sought to identify differential expression between locations was done on normal human macula RPE/choroid and peripheral RPE with a 22 K custom oligonucleotide microarray [40]. The custom array contained RPE and neural retina transcripts identified through public databases and other experimental data from total RNA hybridizations of RPE/choroid, retina, and brain samples. They reported 438 genes differentially expressed between macular and peripheral RPE. Thirty-three were chosen for validation using RT-PCR and 84% were validated. Whitmore et al. demonstrated that 1) in the nasal vs. macular RPE/ choroid comparison, 81 genes showed increased expression in nasal RPE/choroid and 39 genes showed decreased expression; 2) in the temporal vs. macular RPE/choroid comparison, 70 genes were increased in temporal RPE/choroid and 44 genes were decreased; and 3) in the comparison of temporal vs. nasal RPE/choroid, 3 genes were increased in the nasal RPE/choroid and 11 were decreased [32]. It is important to relate gene expression to gene products where possible, because the former is not necessarily a predictor of the latter. Skeie and Mahajan performed proteomic analysis of RPE/choroid in three anatomic regions called fovea, macula, and periphery from three normal human eyes [41]. They defined the foveal region as 4 mm in diameter, thus including both cone-dominant and rod-dominant subregions of the macula (see above). In RPE/choroid, they identified 1) 4204 unique proteins in peripheral, 4595 in macular, and 4409 in foveal punches; 2) high expression of 66 proteins in the foveal RPE/

D

436

choroid only and 251 in the foveal and macular RPE/choroid compared to the peripheral RPE/choroid. Gene ontology analysis showed that proteins with the highest expression in the periphery were enriched in the category of metabolites/energy. Proteins with the highest expression in the macula and fovea were enriched in the response to stimulus, reproduction, homeostasis, or immune system process categories, whereas these categories were absent in the periphery.

E

5.2. Effect of anatomic location on gene expression

T

435

477 478 479 480 481 482 483

6. Transcriptome differences between human retina and RPE/choroid 484 The first attempt to catalogue the genes of retina and RPE used post-mortem retina and RPE from a 44-year-old male and 88-year-old female [17]. Using SAGE analysis, the catalogue of expressed tags in both retina and RPE was completed without measurements of differential expression. In a second study, Bowes Rickman et al. did a more extensive SAGE analysis of macula and peripheral retina genes from five patients and a second SAGE analysis of macula and peripheral RPE/choroid from another five patients [42]. This data was combined with that of Sharon et al. [17] to identify genes enriched in macular retina, peripheral retina, and RPE. They reported three differential expression lists; macula retina enriched, periphery retina enriched, and RPE enriched. The macular retina enriched genes included 270 tags that corresponded to 209 genes. The periphery retina enriched list included 200 tags corresponding to 173 genes and the RPE enriched list contained 143 tags or 123 genes. Using a different approach, Booij et al. set out to find differentially expressed genes between the macula regions of the retina, RPE and choroid from three eyes aged from 63–78 years old [37]. The authors performed two microarrays to compare retina vs. RPE expression and RPE vs. choroid expression. The intersection of the RPE enriched genes from both assays was used as the RPE-Specific genes. Even though they combined RPE enriched genes in both microarrays to find RPEspecific enriched genes, results from the separate microarray experiments were briefly analyzed. From the 33,712 features on the array, 1904 probes were found to have 2.5-fold higher expression in the macular RPE compared to the macular retina, corresponding to 1792 genes. In the RPE vs. choroid assay, 1127 probes were expressed 2.5-fold higher in RPE than choroid, corresponding to 1066 genes.

485

6.1. Comparison of differential expression

513

486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512

We compared the enriched genes from all three above studies to our 514 RNA-Seq results [19]. For macular retina, 65/148 genes reported 515

Please cite this article as: L. Tian, et al., Genomics (2015), http://dx.doi.org/10.1016/j.ygeno.2015.01.008

544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579

C

542 543

E

540 541 Q9

R

535 536

R

533 534

O

531 532

C

529 530

N

527 528

U

525 526

583

In an effort to identify the underlying genetic architecture of AMD and glaucoma, multiple GWAS have been done on both diseases. In addition, direct approaches such as RNA-Seq or microarray transcriptome analysis of the diseased eye have been initiated.

584 585

7.1. Age-related macular degeneration

588

F

Bowes Rickman et al. [42] reported that individual macula-enriched genes were involved in pathways of metabolite precursors and energy substrates including glucose, alcohol, and monosaccharides, but the only pathway found to be enriched was glycolysis. Our list of macula genes from RNA-Seq gave a slightly different result, being significantly enriched in psychiatric, neurological, and chemical dependency disease classes. Most of the pathways enriched were signaling pathways (MAPK signaling, Wnt), as well as glycosphingolipid biosynthesis. The biological processes enriched by RNA-Seq in macular retina included neuron generation/differentiation, ion transport, cell signaling, and vision. Both macular and peripheral retina is dominated by rod photoreceptors, which are very sensitive to light and contribute minimally to color vision [43]. Bowes Rickman et al. reported overrepresentation of genes involved in detection, response, perception of light stimulus, cognition, and eye/structural development for peripheral retina [42]. This overlaps with the results of our RNA-Seq dataset. The list of genes enriched in RPE reported by Bowes Rickman et al. [42] provided no significant pathways, most likely because of the small number of reported genes. On the other hand, Booij et al. reported that RPE enriched genes were represented in cardiovascular, vision, metabolic, and immune diseases [37]. The biological processes for these genes include cell adhesion, vasculature development, immune response, response to wounding, and cell motility. Genes from our RNA-Seq dataset gave similar functional enrichment and biological processes as Booij et al. [37]. The RNA-Seq enriched pathways were similar to the results of Booij et al. [37] with the addition of glycosaminoglycan degradation. In summary, our RNA-Seq dataset demonstrated that 1) for the comparison in retina, the differentially expressed genes that had higher expression levels in macula were significantly enriched for psychiatric and neurological disorders, for chemical dependency in disease category, and for ion transport, transmission of nerve impulse, and synaptic transmission in GO; genes that had higher expression levels in periphery were significantly enriched for translational elongation, sensory perception of light stimulus, and visual perception in GO, and for ribosomes in KEGG; and 2) The overall pattern was different for the comparison in the RPE/choroid layer. Differentially expressed genes that had higher expression levels in macula were significantly enriched for muscle contraction and multicellular organismal process in GO and for dilated cardiomyopathy, hypertrophic cardiomyopathy, and calcium muscle contraction in KEGG; genes that had higher

523 524

7. The transcriptome in AMD and glaucoma

The AMD Gene Consortium recently reported a meta-GWAS analysis of 17,000 patients and 60,000 controls and found 19 loci associated with advanced AMD [44]. Newman et al. completed transcriptome profiling of age-matched normal and AMD donor eyes with the Agilent Whole Human Genome in situ oligonucleotide array [45]. They identified a number of differentially expressed genes in both the RPE/choroid and retina using a permuted p-value b 0.1 and fold change ≥ 1.5 as cutoffs. We compared the genes found near the 19 AMD disease loci by the AMD Gene Consortium with the genes found differentially expressed between diseased and normal eyes by Newman et al. [45] and found an overlap of 24 genes (Table 3). To assess the relationship among these 24 genes, we performed a DAVID analysis. Concentrating on known disease phenotypes, we found an enrichment of terms related to age-related macular degeneration, cardiovascular disease, type 2 diabetes, Alzheimer's disease, and other neurological diseases (Fig. 3). A second smaller cluster was enriched for cancer terms (specifically, breast and colorectal cancer). Functional enrichment analysis of these 24 genes showed enrichment of terms for immune response, regulation of immune processes, and regulation of apoptosis, as well as negative regulation of apoptosis (Figs. 4, 5). The overlap of these susceptibility loci and the identified differentially expressed genes may be coincidental. However, it may also suggest there is yet another layer of complexity to AMD. Cellular immune and inflammatory response appears to be integral to the AMD phenotype. The pathway analysis also reflects a link between AMD and the genetic response that causes inflammation (Fig. 5). While a definitive transcriptome of the AMD disease state has not yet been established, several groups are trying to define its expression profile and underlying genetic cause(s). Whitmore et al. assessed gene expression in the RPE/choroid in nine early AMD and seven normal post-mortem samples using the Affymetrix GeneChip Human Exon 1.0 ST array [46]. They found 75 genes differentially expressed between AMD and normal RPE/choroid with 52 down-regulated and 23 up-regulated in AMD. A DAVID analysis of the down-regulated genes in the AMD samples revealed enrichment in terms associated with vision, sensory perception, and the plasma membrane. We compared the set of differentially expressed genes from normal vs. diseased RPE/choroid samples as reported by Newman et al. [45] and Whitmore et al. [46] and found an overlap of only 15 genes. There are important differences between the two studies. First, the AMD donors selected by Newman et al. included both pre-AMD as well as late AMD [45]. Whitmore et al. [46] focused on early AMD (defined as eyes without GA or CNV) as determined by clinical ophthalmoscopy [44]. We did a DAVID functional annotation analysis of the 15 overlapping genes and found enrichment for terms such as differentiation, developmental protein, glycoprotein, and transport. A number of these genes are also transcriptional-factor targets. This suggests that there are early gene expression changes in early AMD that could serve as potential biomarkers. Whitmore et al. analyzed the top 35 genes reported as RPE-specific by Booij et al. [37] and found 13/35 genes to be increased in AMD samples [46].

O

539

522

R O

6.2. Comparison of pathway analysis

520 521

P

538

518 519

expression levels in peripheral RPE/choroid were enriched for develop- 580 mental in disease category and for biological adhesion, cell adhesion, 581 and cell development in GO [19]. 582

T

537

enriched by Bowes Rickman et al. [42] were enriched in our RNA-Seq gene list (44% agreement). For the peripheral retina, 75/135 genes reported by Bowes Rickman et al. [42] were enriched in both studies (55.6% agreement). Comparing the RPE enriched genes between the studies by Booij et al. [37] and Bowes Rickman et al. [42], we found a small overlap of five genes. When comparing microarray results from Booij et al. [37] with our RNA-Seq data, we found 1021/1358 genes enriched in our RNA-Seq (75.2% agreement). When comparing our RNA-Seq data to the RPE data of Bowes Rickman et al. [42], we find agreement for only ten genes (9.7% agreement). This discrepancy in the degree of overlap between our data and previous studies could be due to the small number of RPE enriched genes in the SAGE experiments, or because the SAGE analysis required RPE enriched genes to fulfill stricter criteria, a requirement not used in the microarray or RNA-Seq studies. Whitmore et al. [32] also compared their dataset with Li et al. [19] for 91 genes implicated by GWAS in eye diseases. They found that the expression values roughly follow the diagonal, and FPKM values were consistently higher for genes in their RPE/choroid samples compared to the RPE/choroid/sclera samples of Li et al. [19]; the authors attributed this finding to enrichment with RPE/choroid in their samples that lacked sclera.

D

516 517

L. Tian et al. / Genomics xxx (2015) xxx–xxx

E

8

Please cite this article as: L. Tian, et al., Genomics (2015), http://dx.doi.org/10.1016/j.ygeno.2015.01.008

586 587

589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639

L. Tian et al. / Genomics xxx (2015) xxx–xxx Table 3 The 24 overlapping genes found between GWAS signals and the Newman et al. study. The related SNP ID from the GWAS study is listed as well as the direction of expression in the tissue.

t3:5

Gene symbol

t3:6 t3:7 t3:8 t3:9 t3:10 t3:11 t3:12 t3:13 t3:14 t3:15 t3:16 t3:17 t3:18 t3:19

Up-regulated in AMD retina DDR1 Discoidin domain receptor tyrosine kinase 1 GTF2H4 General transcription factor IIH, polypeptide 4 STK19 Serine/threonine kinase 19 CFI Complement factor I C3 Complement component 3 APOE Apolipoprotein E VEGFA Vascular endothelial growth factor A C4B Complement component 4B (Chido blood group) APOC1 apolipoprotein C-I CCDC109B Coiled-coil domain containing 109B B3GALTL Beta 1,3-galactosyltransferase-like

t3:20 t3:21 t3:22 t3:23 t3:24 t3:25 t3:26 t3:27

DDR1 IER3 GTF2H4 APOE VEGFA APOC1

t3:28 t3:29 t3:30 t3:31

SYN3 HSPH1 TNXB EHMT2

t3:32 t3:33

Down-regulated in AMD RPE

t3:34 t3:35 t3:36

CFB TNFSF14 TNFRSF10A

t3:37 t3:38

STK19 ADAMTS9

t3:39 t3:40 t3:41 t3:42 t3:43 t3:44

VAV1 NLRC5 HSPH1 TUBB TNXB C2

Related SNP rs3130783 rs3130783 rs429608 rs4698775 rs2230199 rs4420638 rs943080 rs429608 rs4420638 rs4698775 rs9542236

Glaucoma is a neurodegenerative disease resulting in damage to the optic nerve. The trabecular meshwork (TM) in the anterior segment of the eye is a series of canals that drain aqueous fluid. These canals narrow with age causing an increase in intraocular pressure (IOP) that is associated with glaucoma. The change in gene expression associated with the narrowing of TM canals has been studied to identify potential causes of glaucoma. Liu et al. analyzed the expression profile of TM from 13 normal and 15 glaucomatous eye samples, ages 40 to 86 years old, with Illumina Human WG-6 BeadChips [47]. A total of 483 differentially expressed genes existed between normal and glaucomatous TM. Gene ontology analysis revealed enrichment for terms such as signal peptide, secreted, glycoprotein, extracellular region, cell adhesion, phosphate transport, and calcium ion binding. Also reported in this paper was one Primary Open Angle Glaucoma (POAG) case with a Myocilin (MYOC) mutation, a gene known to be associated with POAG. Differential expression between the POAG eye with MYOC mutation and the POAG eyes without MYOC mutation revealed 55 genes with fold changes ≥ 2-fold. Myocilin was not found to be differentially expressed between the MYOC and non-MYOC POAG cases. We compared the genes located in the POAG GWAS loci [48] with the 483 genes differentially expressed in the TM [47] and found no overlap.

641 642

R O

rs3130783 rs3130783 rs3130783 rs4420638 rs943080 rs4420638

Down-regulated in AMD retina

P

rs5749482 rs9542236 rs429608 rs429608

rs429608 rs6795735

rs2230199 rs1864163 rs9542236 rs3130783 rs429608 rs429608

U

N C O

R

R

It is becoming common practice to make expression information available via databases. Such resources, when maintained, can be of significant value for storing and querying genomic data. Examples of genomic databases that are commonly used include NCBI, dbSNP, Genome Browser, OMIM, GEO, UniProt, and others. Lesser known are the databases for a specific tissue or cell type. The RetinaCentral database (http://www.retinacentral.org), a repository of information on the mammalian retinome, was created to collect and store information on genes and proteins expressed in the retina. It encompasses data from 27 studies that employed a range of techniques including cDNA sequencing, microarray, and SAGE [49]. The web interface allows querying and filtering based on chromosomal location,

E

T

rs429608 rs2230199 rs13278062

E

Complement factor B Tumor necrosis factor superfamily, member 14 Tumor necrosis factor receptor superfamily, member 10a Serine/threonine kinase 19 ADAM metallopeptidase with thrombospondin type 1 motif, 9 vav 1 guanine nucleotide exchange factor NLR family, CARD domain containing 5 Heat shock 105 kDa/110 kDa protein 1 Tubulin, beta class I Tenascin XB Complement component 2

643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662

8. Bioinformatic databases for human retina and RPE/choroid 663 transcriptomes 664

D

Synapsin 3 Heat shock 105 kDa/110 kDa protein 1 Tenascin XB Euchromatic histone-lysine N-methyltransferase 2

640

O

Up-regulated in AMD RPE Discoidin domain receptor tyrosine kinase 1 Immediate early response 3 General transcription factor IIH, polypeptide 4 Apolipoprotein E Vascular endothelial growth factor A Apolipoprotein C-I

7.2. Glaucoma

F

Gene name

C

t3:1 t3:2 t3:3 t3:4 Q4

9

Fig. 3. Disease term enrichment of 24 genes near AMD disease loci. X axis shows the disease terms; Y axis shows the number of genes enriched in the specific disease term.

Please cite this article as: L. Tian, et al., Genomics (2015), http://dx.doi.org/10.1016/j.ygeno.2015.01.008

665 666 667 668 669 670 671 672 673 674 675 676

L. Tian et al. / Genomics xxx (2015) xxx–xxx

R O

O

F

10

Fig. 4. Functional enrichment analysis of 24 genes near AMD disease loci. X axis shows the biological function categories; Y axis shows the number of genes enriched in the specific function category.

691 692 693 694 695 696 697 698 699 700 701

702 703

9. Conclusions

731

The transcriptome encompasses all RNA molecules within a cell or a population of cells. It represents the link between an organism's genome and its phenotype. Tools for studying RNAs have been available for decades and include Northern blots, reverse-transcription PCR (RT-PCR), and expressed sequence tags (ESTs). The transition to a rapid and high-throughput quantification of the transcriptome became possible with the development of serial analysis of gene expression (SAGE), gene expression cDNA microarrays, and RNA-Seq. Recently, RNA-Seq has shown a more complex retinal transcriptome than previously reported. Farkas et al. [29] and our group [19] performed comprehensive analysis of the retina and RPE gene expression

732 733

E

D

P

RETINOBASE (http://www.lbgi.fr/wikili/index.php/RetinoBase) is a database designed for the analysis and visualization of microarray gene expression profiles of the retina summarizing datasets from 28 different microarray experiments. Five different model systems are represented in these datasets including drosophila, zebra fish, rat, mouse, and human [52]. It appears that this database is no longer active. EyeSAGE (http://neibank.nei.nih.gov/EyeSAGE/index.shtml) is a relational database that provides a comprehensive picture of gene expression in the human macula and includes transcription profiles of the retina and RPE/choroid [42]. The database was generated from SAGE (serial analysis of gene expression) libraries, both long and short SAGE, and also includes some microarray data. Validation was carried out for several known retina-specific genes to corroborate differential expression that was detected by either SAGE or microarray. The database is available for download in Microsoft Access format on NEIBank's website. In addition, the NEIBank's EyeBrowse can be used to visualize expression levels for various genes or candidate disease regions from the EyeSAGE datasets. Data from an Affymetrix GeneChip Human Exon 1.0 ST analysis of ocular tissues is available [35]. Based on six post-mortem eyes, the tissue profiling includes retina, optic nerve head, optic nerve, ciliary body, TM, sclera, lens, cornea, RPE/choroid, and iris. Normalization was done using the Probe Logarithmic Intensity Error (PLIER) method. The web interface allows users to query both transcript and exon level PLIER expression estimates from various ocular tissues. Finally, RetNet (https://sph.uth.edu/retnet/) is a database based on published literature was created for the purpose of providing genes and loci associated with retinal diseases [53]. It is a collection of accepted retinal disease genes/loci found in the literature.

T

C

E

689 690

R

687 688

R

685 686

O

683 684

C

681 682

N

679 680

category (such as retinal localization or retina/RPE-specific), study, and phenotype. Although last updated in 2005, it still serves as a useful reference for retina/RPE gene expression profiling and the selection of candidate genes for further study. Alternative splicing adds another layer of complexity to the human transcriptome and interest in isoforms has grown substantially. The Human Alternative Splicing Database (HASDB) is a database of tissuespecific alternative splicing for human genes and includes 46 different human tissues including retina [50]. Bayesian statistics are used to create a scoring function of tissue specificity (TS), as well as a measure of robustness (rTS) to assess stability of the TS value. Extensive validation was performed for tissue specificity. Although the HASDB itself can no longer be found online, HASDB data are incorporated into a successor database known as the Alternative Splicing Annotation Project (ASAP). Katsanis et al. have published an excellent resource detailing the design of several computational tools for the manipulation of dbEST to predict clusters in specific tissues [51]. In addition to providing a mechanism to manipulate dbEST to find ESTs specific to a tissue, the study focused on the retina and by association positional candidates for disorders of the retina. They combined information from the identified retina-specific ESTs with the “genome-mapping locations for retinopathy linked loci” to identify positional candidates. This information was incorporated into a database called RETBASE. While this repository was useful, it is no longer available online.

U

677 678

Fig. 5. Pathway analysis of 24 genes near AMD disease loci. X axis lists top 3 pathways; Y axis shows the number of genes enriched in the specific pathway.

Please cite this article as: L. Tian, et al., Genomics (2015), http://dx.doi.org/10.1016/j.ygeno.2015.01.008

704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730

734 735 736 737 738 739 740 741 742

L. Tian et al. / Genomics xxx (2015) xxx–xxx

751

The authors declare that they have no competing interests.

752

Authors' contributions

753 754 755

All authors contributed to the manuscript. LF, KLK, and ASB wrote the first draft and all other authors amended the manuscript. All authors read and approved the final manuscript.

756

Acknowledgments

[14] [15] [16]

[17]

757 758

We thank our laboratory members for helpful discussion of the manuscript. CAC acknowledges support from: NEI EY06109; Arnold and 759 Mabel Beckman Initiative for Macular Research; EyeSight Foundation of 760 Q10 Alabama, Research to Prevent Blindness. DES acknowledges support 761 from: NEI EY023164 and EY020483; Arnold and Mabel Beckman 762 Initiative for Macular Research.

[1] V.G. Cheung, L.K. Conlin, T.M. Weber, M. Arcaro, K.Y. Jen, M. Morley, R.S. Spielman, Natural variation in human gene expression assessed in lymphoblastoid cells, Nat. Genet. 33 (3) (2003) 422–425. [2] E.E. Schadt, S.A. Monks, T.A. Drake, A.J. Lusis, N. Che, V. Colinayo, T.G. Ruff, S.B. Milligan, J.R. Lamb, G. Cavet, P.S. Linsley, M. Mao, R.B. Stoughton, S.H. Friend, Genetics of gene expression surveyed in maize, mouse and man, Nature 422 (6929) (2003) 297–302. [3] B.E. Stranger, A.C. Nica, M.S. Forrest, A. Dimas, C.P. Bird, C. Beazley, C.E. Ingle, M. Dunning, P. Flicek, D. Koller, S. Montgomery, S. Tavare, P. Deloukas, E.T. Dermitzakis, Population genomics of human gene expression, Nat. Genet. 39 (10) (2007) 1217–1224. [4] V. Emilsson, G. Thorleifsson, B. Zhang, A.S. Leonardson, F. Zink, J. Zhu, S. Carlson, A. Helgason, G.B. Walters, S. Gunnarsdottir, M. Mouy, V. Steinthorsdottir, G.H. Eiriksdottir, G. Bjornsdottir, I. Reynisdottir, D. Gudbjartsson, A. Helgadottir, A. Jonasdottir, A. Jonasdottir, U. Styrkarsdottir, S. Gretarsdottir, K.P. Magnusson, H. Stefansson, R. Fossdal, K. Kristjansson, H.G. Gislason, T. Stefansson, B.G. Leifsson, U. Thorsteinsdottir, J.R. Lamb, J.R. Gulcher, M.L. Reitman, A. Kong, E.E. Schadt, K. Stefansson, Genetics of gene expression and its effect on disease, Nature 452 (7186) (2008) 423–428. [5] E.L. Heinzen, D. Ge, K.D. Cronin, J.M. Maia, K.V. Shianna, W.N. Gabriel, K.A. WelshBohmer, C.M. Hulette, T.N. Denny, D.B. Goldstein, Tissue-specific genetic control of splicing: implications for the study of complex traits, PLoS Biol. 6 (12) (2008) e1. [6] J.C. Marioni, C.E. Mason, S.M. Mane, M. Stephens, Y. Gilad, RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays, Genome Res. 18 (9) (2008) 1509–1517. [7] C.A. Curcio, J.D. Messinger, K.R. Sloan, A. Mitra, G. McGwin, R.F. Spaide, Human chorioretinal layer thicknesses measured in macula-wide, high-resolution histologic sections, Invest. Ophthalmol. Vis. Sci. 52 (7) (2011) 3943–3954. [8] S.A. Zweifel, Y. Imamura, T.C. Spaide, T. Fujiwara, R.F. Spaide, Prevalence and significance of subretinal drusenoid deposits (reticular pseudodrusen) in agerelated macular degeneration, Ophthalmology 117 (9) (2010) 1775–1781. [9] S.A. Zweifel, R.F. Spaide, C.A. Curcio, G. Malek, Y. Imamura, Reticular pseudodrusen are subretinal drusenoid deposits, Ophthalmology 117 (2) (2010) 303-12.e1. [10] C.A. Curcio, J.D. Messinger, K.R. Sloan, G. McGwin, N.E. Medeiros, R.F. Spaide, Subretinal drusenoid deposits in non-neovascular age-related macular degeneration: morphology, prevalence, topography, and biogenesis model, Retina 33 (2) (2013) 265–276. [11] R.S. Ramrattan, T.L. van der Schaft, C.M. Mooy, W.C. de Bruijn, P.G. Mulder, P.T. de Jong, Morphometric analysis of Bruch's membrane, the choriocapillaris, and the choroid in aging, Invest. Ophthalmol. Vis. Sci. 35 (6) (1994) 2857–2864. [12] M. Schena, D. Shalon, R.W. Davis, P.O. Brown, Quantitative monitoring of gene expression patterns with a complementary DNA microarray, Science 270 (5235) (1995) 467–470. [13] MAQC Consortium, L. Shi, L.H. Reid, W.D. Jones, R. Shippy, J.A. Warrington, S.C. Baker, P.J. Collins, F. de Longueville, E.S. Kawasaki, K.Y. Lee, Y. Luo, Y.A. Sun, J.C. Willey, R.A. Setterquist, G.M. Fischer, W. Tong, Y.P. Dragan, D.J. Dix, F.W. Frueh, F.M. Goodsaid, D. Herman, R.V. Jensen, C.D. Johnson, E.K. Lobenhofer, R.K. Puri, U. Schrf, J. ThierryMieg, C. Wang, M. Wilson, P.K. Wolber, L. Zhang, S. Amur, W. Bao, C.C. Barbacioru,

D

764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 Q11 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811

[20]

C

E

R

R

N C O

U

[21]

E

References

[19]

T

763

[18]

F

Competing interests

747 748

O

750

745 746

A.B. Lucas, V. Bertholet, C. Boysen, B. Bromley, D. Brown, A. Brunner, R. Canales, X.M. Cao, T.A. Cebula, J.J. Chen, J. Cheng, T.M. Chu, E. Chudin, J. Corson, J.C. Corton, L.J. Croner, C. Davies, T.S. Davison, G. Delenstarr, X. Deng, D. Dorris, A.C. Eklund, X.H. Fan, H. Fang, S. Fulmer-Smentek, J.C. Fuscoe, K. Gallagher, W. Ge, L. Guo, X. Guo, J. Hager, P.K. Haje, J. Han, T. Han, H.C. Harbottle, S.C. Harris, E. Hatchwell, C.A. Hauser, S. Hester, H. Hong, P. Hurban, S.A. Jackson, H. Ji, C.R. Knight, W.P. Kuo, J.E. LeClerc, S. Levy, Q.Z. Li, C. Liu, Y. Liu, M.J. Lombardi, Y. Ma, S.R. Magnuson, B. Maqsodi, T. McDaniel, N. Mei, O. Myklebost, B. Ning, N. Novoradovskaya, M.S. Orr, T.W. Osborn, A. Papallo, T.A. Patterson, R.G. Perkins, E.H. Peters, R. Peterson, K.L. Philips, P.S. Pine, L. Pusztai, F. Qian, H. Ren, M. Rosen, B.A. Rosenzweig, R.R. Samaha, M. Schena, G.P. Schroth, S. Shchegrova, D.D. Smith, F. Staedtler, Z. Su, H. Sun, Z. Szallasi, Z. Tezak, D. Thierry-Mieg, K.L. Thompson, I. Tikhonova, Y. Turpaz, B. Vallanat, C. Van, S.J. Walker, S.J. Wang, Y. Wang, R. Wolfinger, A. Wong, J. Wu, C. Xiao, Q. Xie, J. Xu, W. Yang, L. Zhang, S. Zhong, Y. Zong, W. Slikker Jr., The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements, Nat. Biotechnol. 24 (9) (2006) 1151–1161. J.H. Malone, B. Oliver, Microarrays, deep sequencing and the true measure of the transcriptome, BMC Biol. 9 (2011) (34-7007-9-34). V.E. Velculescu, L. Zhang, B. Vogelstein, K.W. Kinzler, Serial analysis of gene expression, Science 270 (5235) (1995) 484–487. L. Zhang, W. Zhou, V.E. Velculescu, S.E. Kern, R.H. Hruban, S.R. Hamilton, B. Vogelstein, K.W. Kinzler, Gene expression profiles in normal and cancer cells, Science 276 (5316) (1997) 1268–1272. D. Sharon, S. Blackshaw, C.L. Cepko, T.P. Dryja, Profile of the genes expressed in the human peripheral retina, macula, and retinal pigment epithelium determined through serial analysis of gene expression (SAGE), Proc. Natl. Acad. Sci. U. S. A. 99 (1) (2002) 315–320. K.J. Malik, C.D. Chen, T.W. Olsen, Stability of RNA from the retina and retinal pigment epithelium in a porcine model simulating human eye bank conditions, Invest. Ophthalmol. Vis. Sci. 44 (6) (2003) 2730–2735. M. Li, C. Jia, K.L. Kazmierkiewicz, A.S. Bowman, L. Tian, Y. Liu, N.A. Gupta, H.V. Gudiseva, S.S. Yee, M. Kim, T. Dentchev, J.A. Kimble, J.S. Parker, J.D. Messinger, H. Hakonarson, C.A. Curcio, D. Stambolian, Comprehensive analysis of gene expression in human retina and supporting tissues, Hum. Mol. Genet. (2014). Age-Related Eye Disease Study Research Group, The Age-Related Eye Disease Study system for classifying age-related macular degeneration from stereoscopic color fundus photographs: the Age-Related Eye Disease Study Report Number 6, Am J. Ophthalmol. 132 (5) (2001) 668–681. C.A. Curcio, N.E. Medeiros, C.L. Millican, The Alabama age-related macular degeneration grading system for donor eyes, Invest. Ophthalmol. Vis. Sci. 39 (7) (1998) 1085–1096. T.W. Olsen, X. Feng, The Minnesota grading system of eye bank eyes for age-related macular degeneration, Invest. Ophthalmol. Vis. Sci. 45 (12) (2004) 4484–4490. S.A. Zweifel, M. Engelbert, K. Laud, R. Margolis, R.F. Spaide, K.B. Freund, Outer retinal tubulation: a novel optical coherence tomography finding, Arch. Ophthalmol. 127 (12) (2009) 1596–1602. N. Bagheri, B.A. Bell, V.L. Bonilha, J.G. Hollyfield, Imaging human postmortem eyes with SLO and OCT, Adv. Exp. Med. Biol. 723 (2012) 479–488. N.H. Brown, A.F. Koreishi, M. McCall, J.A. Izatt, C.B. Rickman, C.A. Toth, Developing SDOCT to assess donor human eyes prior to tissue sectioning for research, Graefes Arch. Clin. Exp. Ophthalmol. 247 (8) (2009) 1069–1080. A. Swaroop, D.J. Zack, Transcriptome analysis of the retina, Genome Biol. 3 (8) (2002) (REVIEWS1022). S. Yoshida, B.M. Yashar, S. Hiriyanna, A. Swaroop, Microarray analysis of gene expression in the aging human retina, Invest. Ophthalmol. Vis. Sci. 43 (8) (2002) 2554–2560. H. Cai, M.A. Fields, R. Hoshino, L.V. Priore, Effects of aging and anatomic location on gene expression in human retina, Front. Aging Neurosci. 4 (2012) 8. M.H. Farkas, G.R. Grant, J.A. White, M.E. Sousa, M.B. Consugar, E.A. Pierce, Transcriptome analyses of the human retina identify unprecedented transcript diversity and 3.5 Mb of novel transcribed sequence via significant alternative splicing and novel genes, BMC Genomics 14 (1) (2013) 486. T.L. Young, F. Hawthorne, S. Feng, X. Luo, E. St Germain, M. Wang, R. Metlapally, Whole genome expression profiling of normal human fetal and adult ocular tissues, Exp. Eye Res. 116 (2013) 265–278. D.M. Hornan, S.N. Peirson, A.J. Hardcastle, R.S. Molday, M.E. Cheetham, A.R. Webster, Novel retinal and cone photoreceptor transcripts revealed by human macular expression profiling, Invest. Ophthalmol. Vis. Sci. 48 (12) (2007) 5388–5396. S.S. Whitmore, A.H. Wagner, A.P. DeLuca, A.V. Drack, E.M. Stone, B.A. Tucker, S. Zeng, T.A. Braun, R.F. Mullins, T.E. Scheetz, Transcriptomic analysis across nasal, temporal, and macular regions of human neural retina and RPE/choroid by RNA-Seq, Exp. Eye Res. 129 (2014) 93–106. A. Mortazavi, B.A. Williams, K. McCue, L. Schaeffer, B. Wold, Mapping and quantifying mammalian transcriptomes by RNA-Seq, Nat. Methods 5 (7) (2008) 621–628. A. Sirbu, G. Kerr, M. Crane, H.J. Ruskin, RNA-Seq vs dual- and single-channel microarray data: sensitivity analysis for differential expression and clustering, PLoS One 7 (12) (2012) e50986. A.H. Wagner, V.N. Anand, W.H. Wang, J.E. Chatterton, D. Sun, A.R. Shepard, N. Jacobson, I.H. Pang, A.P. Deluca, T.L. Casavant, T.E. Scheetz, R.F. Mullins, T.A. Braun, A.F. Clark, Exon-level expression profiling of ocular tissues, Exp. Eye Res. 111 (2013) 105–111. J.C. Booij, S. van Soest, S.M. Swagemakers, A.H. Essing, A.J. Verkerk, P.J. van der Spek, T.G. Gorgels, A.A. Bergen, Functional annotation of the human retinal pigment epithelium transcriptome, BMC Genomics 10 (2009) (164-2164-10-164). J.C. Booij, J.B. ten Brink, S.M. Swagemakers, A.J. Verkerk, A.H. Essing, P.J. van der Spek, A.A. Bergen, A new strategy to identify and annotate human RPE-specific gene expression, PLoS One 5 (5) (2010) e9341.

R O

749

profiles using RNA-Seq. Both studies revealed large differences in expression between these layers of the posterior eye and a smaller, but significant number of differences between macula and periphery of each layer. With even larger scale RNA-Seq datasets, it will be possible to assess the impact of gene expression on molecular networks that will lead to an understanding of how RNA/protein variation results in disease.

P

743 744

11

[22] [23]

[24] [25]

[26] [27]

[28] [29]

[30]

[31]

[32]

[33] [34]

[35]

[36]

[37]

Please cite this article as: L. Tian, et al., Genomics (2015), http://dx.doi.org/10.1016/j.ygeno.2015.01.008

812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 Q12 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897

[49]

[50] [51]

[52]

[53]

F

[48]

O

[47]

R O

[46]

P

[45]

A.R. Webster, H.E. Wichmann, T.W. Winkler, T.Y. Wong, A.F. Wright, D. Zelenika, M. Zhang, L. Zhao, K. Zhang, M.L. Klein, G.S. Hageman, G.M. Lathrop, K. Stefansson, R. Allikmets, P.N. Baird, M.B. Gorin, J.J. Wang, C.C. Klaver, J.M. Seddon, M.A. Pericak-Vance, S.K. Iyengar, J.R. Yates, A. Swaroop, B.H. Weber, M. Kubo, M.M. Deangelis, T. Leveillard, U. Thorsteinsdottir, J.L. Haines, L.A. Farrer, I.M. Heid, G.R. Abecasis, AMD Gene Consortium, Seven new loci associated with age-related macular degeneration, Nat. Genet. 45 (4) (2013) 433–439 (439e1-2). A.M. Newman, N.B. Gallo, L.S. Hancox, N.J. Miller, C.M. Radeke, M.A. Maloney, J.B. Cooper, G.S. Hageman, D.H. Anderson, L.V. Johnson, M.J. Radeke, Systems-level analysis of age-related macular degeneration reveals global biomarkers and phenotype-specific functional networks, Genome Med. 4 (2) (2012) 16. S.S. Whitmore, T.A. Braun, J.M. Skeie, C.M. Haas, E.H. Sohn, E.M. Stone, T.E. Scheetz, R.F. Mullins, Altered gene expression in dry age-related macular degeneration suggests early loss of choroidal endothelial cells, Mol. Vis. 19 (2013) 2274–2297. Y. Liu, R.R. Allingham, X. Qin, D. Layfield, A.E. Dellinger, J. Gibson, J. Wheeler, A.E. Ashley-Koch, W.D. Stamer, M.A. Hauser, Gene expression profile in human trabecular meshwork from patients with primary open-angle glaucoma, Invest. Ophthalmol. Vis. Sci. 54 (9) (2013) 6382–6389. J.N. Cooke Bailey, L. Sobrin, M.A. Pericak-Vance, J.L. Haines, C.J. Hammond, J.L. Wiggs, Advances in the genomics of common eye diseases, Hum. Mol. Genet. 22 (R1) (2013) R59–R65. H.L. Schulz, T. Goetz, J. Kaschkoetoe, B.H. Weber, The retinome — defining a reference transcriptome of the adult mammalian retina/retinal pigment epithelium, BMC Genomics 5 (1) (2004) 50. Q. Xu, B. Modrek, C. Lee, Genome-wide detection of tissue-specific alternative splicing in the human transcriptome, Nucleic Acids Res. 30 (17) (2002) 3754–3766. N. Katsanis, K.C. Worley, G. Gonzalez, S.J. Ansley, J.R. Lupski, A computational/ functional genomics approach for the enrichment of the retinal transcriptome and the identification of positional candidate retinopathy genes, Proc. Natl. Acad. Sci. U. S. A. 99 (22) (2002) 14326–14331. R.K. Kalathur, N. Gagniere, G. Berthommier, L. Poidevin, W. Raffelsberger, R. Ripp, T. Leveillard, O. Poch, RETINOBASE: a web database, data mining and analysis platform for gene expression data on retina, BMC Genomics (2008) 9 (2082164-9-208). S.P. Daiger, B.F. Rossiter, J. Greenberg, A. Christoffels, W. Hide, Data services and software for identifying genes and mutations causing retinal degeneration, Invest. Ophthalmol. Vis. Sci. 12 (39) (1997) 295.

N

C

O

R

R

E

C

T

E

[38] N.V. Strunnikova, A. Maminishkis, J.J. Barb, F. Wang, C. Zhi, Y. Sergeev, W. Chen, A.O. Edwards, D. Stambolian, G. Abecasis, A. Swaroop, P.J. Munson, S.S. Miller, Transcriptome analysis and molecular signature of human retinal pigment epithelium, Hum. Mol. Genet. 19 (12) (2010) 2468–2486. [39] M.J. Radeke, K.E. Peterson, L.V. Johnson, D.H. Anderson, Disease susceptibility of the human macula: differential gene transcription in the retinal pigmented epithelium/ choroid, Exp. Eye Res. 85 (3) (2007) 366–380. [40] S.S. van Soest, G.M. de Wit, A.H. Essing, J.B. ten Brink, W. Kamphuis, P.T. de Jong, A.A. Bergen, Comparison of human retinal pigment epithelium gene expression in macula and periphery highlights potential topographic differences in Bruch's membrane, Mol. Vis. 13 (2007) 1608–1617. [41] J.M. Skeie, V.B. Mahajan, Proteomic landscape of the human choroid-retinal pigment epithelial complex, JAMA Ophthalmol. 132 (11) (2014) 1271–1281. [42] C. Bowes Rickman, J.N. Ebright, Z.J. Zavodni, L. Yu, T. Wang, S.P. Daiger, G. Wistow, K. Boon, M.A. Hauser, Defining the human macula transcriptome and candidate retinal disease genes using EyeSAGE, Invest. Ophthalmol. Vis. Sci. 47 (6) (2006) 2305–2316. [43] C.A. Curcio, K.R. Sloan, R.E. Kalina, A.E. Hendrickson, Human photoreceptor topography, J. Comp. Neurol. 292 (4) (1990) 497–523. [44] L.G. Fritsche, W. Chen, M. Schu, B.L. Yaspan, Y. Yu, G. Thorleifsson, D.J. Zack, S. Arakawa, V. Cipriani, S. Ripke, R.P. Igo Jr., G.H. Buitendijk, X. Sim, D.E. Weeks, R.H. Guymer, J.E. Merriam, P.J. Francis, G. Hannum, A. Agarwal, A.M. Armbrecht, I. Audo, T. Aung, G.R. Barile, M. Benchaboune, A.C. Bird, P.N. Bishop, K.E. Branham, M. Brooks, A.J. Brucker, W.H. Cade, M.S. Cain, P.A. Campochiaro, C.C. Chan, C.Y. Cheng, E.Y. Chew, K.A. Chin, I. Chowers, D.G. Clayton, R. Cojocaru, Y.P. Conley, B.K. Cornes, M.J. Daly, B. Dhillon, A.O. Edwards, E. Evangelou, J. Fagerness, H.A. Ferreyra, J.S. Friedman, A. Geirsdottir, R.J. George, C. Gieger, N. Gupta, S.A. Hagstrom, S.P. Harding, C. Haritoglou, J.R. Heckenlively, F.G. Holz, G. Hughes, J.P. Ioannidis, T. Ishibashi, P. Joseph, G. Jun, Y. Kamatani, N. Katsanis, C.N. Keilhauer, J.C. Khan, I.K. Kim, Y. Kiyohara, B.E. Klein, R. Klein, J.L. Kovach, I. Kozak, C.J. Lee, K.E. Lee, P. Lichtner, A.J. Lotery, T. Meitinger, P. Mitchell, S. Mohand-Said, A.T. Moore, D.J. Morgan, M.A. Morrison, C.E. Myers, A.C. Naj, Y. Nakamura, Y. Okada, A. Orlin, M.C. Ortube, M.I. Othman, C. Pappas, K.H. Park, G.J. Pauer, N.S. Peachey, O. Poch, R.R. Priya, R. Reynolds, A.J. Richardson, R. Ripp, G. Rudolph, E. Ryu, J.A. Sahel, D.A. Schaumberg, H.P. Scholl, S.G. Schwartz, W.K. Scott, H. Shahid, H. Sigurdsson, G. Silvestri, T.A. Sivakumaran, R.T. Smith, L. Sobrin, E.H. Souied, D.E. Stambolian, H. Stefansson, G.M. Sturgill-Short, A. Takahashi, N. Tosakulwong, B.J. Truitt, E.E. Tsironi, A.G. Uitterlinden, C.M. van Duijn, L. Vijaya, J.R. Vingerling, E.N. Vithana,

U

898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 Q13 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934

L. Tian et al. / Genomics xxx (2015) xxx–xxx

D

12

Please cite this article as: L. Tian, et al., Genomics (2015), http://dx.doi.org/10.1016/j.ygeno.2015.01.008

935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971