Defining the human macula transcriptome and candidate retinal disease genes using EyeSAGE

Recent Achievements in the Heterogeneity of Mammalian and Human Retinal Pigment Epithelium: In Search of a Stem Cell

Retinal pigment epithelium (RPE) cells are important fundamentally for the development and function of the retina. In this regard, the study of the morphological and molecular properties of RPE cells, as well as their regenerative capabilities, is of particular importance for biomedicine. However, these studies are complicated by the fact that, despite the external morphological similarity of RPE cells, the RPE is a population of heterogeneous cells, the molecular genetic properties of which have begun to be revealed by sequencing methods only in recent years. This review carries out an analysis of the data from morphological and molecular genetic studies of the heterogeneity of RPE cells in mammals and humans, which reveals the individual differences in the subpopulations of RPE cells and the possible specificity of their functions. Particular attention is paid to discussing the properties of "stemness," proliferation, and plasticity in the RPE, which may be useful for uncovering the mechanisms of retinal diseases associated with pathologies of the RPE and finding new ways of treating them.

Transcriptomic and Chromatin Accessibility Analysis of the Human Macular and Peripheral Retinal Pigment Epithelium at the Single-Cell Level

Some human retinal diseases are characterized by pathology that is restricted to specific cell types and to specific regions of the eye. Several disease entities either selectively affect or spare the macula, the retina region at the center of the posterior pole. Photoreceptor cells in the macula are involved in high-acuity vision and require metabolic support from non-neuronal cell types. Some macular diseases involve the retinal pigment epithelium (RPE), an epithelial cell layer with several metabolic-support functions essential for the overlying photoreceptors. In the current study, the ways in which RPE confers region-specific disease susceptibility were determined by examining heterogeneity within RPE tissue from human donors. RPE nuclei from the macular and peripheral retina were profiled using joint single-nucleus RNA and ATAC sequencing. The expression of several genes differed between macular and peripheral RPE. Region-specific ATAC peaks were found, suggesting regulatory elements used exclusively by macular or peripheral RPE. Across anatomic regions, subpopulations of RPE were identified that appeared to have differential levels of expression of visual cycle genes. Finally, loci associated with age-related macular degeneration were examined for a better understanding of RPE-specific disease phenotypes. These findings showed variations in the regulation of gene expression in the human RPE by region and subpopulation, and provide a source for a better understanding of the molecular basis of macular disease.

Single-cell RNA sequencing in vision research: Insights into human retinal health and disease

Gene expression provides valuable insight into cell function. As such, vision researchers have frequently employed gene expression studies to better understand retinal physiology and disease. With the advent of single-cell RNA sequencing, expression experiments provide an unparalleled resolution of information. Instead of studying aggregated gene expression across all cells in a heterogenous tissue, single-cell technology maps RNA to an individual cell, which facilitates grouping of retinal and choroidal cell types for further study. Single-cell RNA sequencing has been quickly adopted by both basic and translational vision researchers, and single-cell level gene expression has been studied in the visual systems of animal models, retinal organoids, and primary human retina, RPE, and choroid. These experiments have generated detailed atlases of gene expression and identified new retinal cell types. Likewise, single-cell RNA sequencing investigations have characterized how gene expression changes in the setting of many retinal diseases, including how choroidal endothelial cells are altered in age-related macular degeneration. In addition, this technology has allowed vision researchers to discover drivers of retinal development and model rare retinal diseases with induced pluripotent stem cells. In this review, we will overview the growing number of single-cell RNA sequencing studies in the field of vision research. We will summarize experimental considerations for designing single-cell RNA sequencing experiments and highlight important advancements in retinal, RPE, choroidal, and retinal organoid biology driven by this technology. Finally, we generalize these findings to genes involved in retinal degeneration and outline the future of single-cell expression experiments in studying retinal disease.

The Pathophysiological Significance of Fibulin-3

Fibulin-3 (also known as EGF-containing fibulin extracellular matrix protein 1 (EFEMP1)) is a secreted extracellular matrix glycoprotein, encoded by the EFEMP1 gene that belongs to the eight-membered fibulin protein family. It has emerged as a functionally unique member of this family, with a diverse array of pathophysiological associations predominantly centered on its role as a modulator of extracellular matrix (ECM) biology. Fibulin-3 is widely expressed in the human body, especially in elastic-fibre-rich tissues and ocular structures, and interacts with enzymatic ECM regulators, including tissue inhibitor of metalloproteinase-3 (TIMP-3). A point mutation in EFEMP1 causes an inherited early-onset form of macular degeneration called Malattia Leventinese/Doyne honeycomb retinal dystrophy (ML/DHRD). EFEMP1 genetic variants have also been associated in genome-wide association studies with numerous complex inherited phenotypes, both physiological (namely, developmental anthropometric traits) and pathological (many of which involve abnormalities of connective tissue function). Furthermore, EFEMP1 expression changes are implicated in the progression of numerous types of cancer, an area in which fibulin-3 has putative significance as a therapeutic target. Here we discuss the potential mechanistic roles of fibulin-3 in these pathologies and highlight how it may contribute to the development, structural integrity, and emergent functionality of the ECM and connective tissues across a range of anatomical locations. Its myriad of aetiological roles positions fibulin-3 as a molecule of interest across numerous research fields and may inform our future understanding and therapeutic approach to many human diseases in clinical settings.

Spectacle: An interactive resource for ocular single-cell RNA sequencing data analysis

Single-cell RNA sequencing has revolutionized ocular gene expression studies. This technology has enabled researchers to identify expression signatures for rare cell types and characterize how gene expression changes across biological conditions, such as topographic region or disease status. However, sharing single-cell RNA sequencing results remains a major obstacle, particular for individuals without a computational background. To address these limitations, we developed Spectacle, an interactive web-based resource for exploring previously published single-cell RNA sequencing data from ocular studies. Spectacle is powered by a locally developed R package, cellcuratoR, which utilizes the Shiny framework in R to generate interactive visualizations for single-cell expression data. Spectacle contains five pre-processed ocular single-cell RNA sequencing data sets and is accessible via the web at OcularGeneExpression.org/singlecell. With Spectacle, users can interactively identify which cell types express a gene of interest, detect transcriptomic subpopulations within a cell type, and perform highly flexible differential expression analyses. The freely-available Spectacle system reduces the bioinformatic barrier for interacting with rich single-cell RNA sequencing studies from ocular tissues, making it easy to quickly identify cell types that express a gene of interest.

Bone Morphogenetic Protein (BMP)4 But Not BMP2 Disrupts the Barrier Integrity of Retinal Pigment Epithelia and Induces Their Migration: A Potential Role in Neovascular Age-Related Macular Degeneration

Disruption of retinal pigment epithelial (RPE) barrier integrity and RPE migration are hallmark features in neovascular age-related macular degeneration (nAMD), but the underlying causes and pathophysiology are not completely well-defined. Herein, we aimed to evaluate the effect of bone morphogenetic proteins (BMPs) on the barrier function and migration of RPE. In particular, we investigated the role of BMP2 and BMP4 in these processes as our analysis of RNA-sequencing (seq) data from human donor eyes demonstrated that they are highly differentially expressed BMP members in macular RPE/choroid versus macular retina. We used electrical cell-substrate impedance sensing (ECIS) system to monitor precisely in real time the barrier integrity and migration of ARPE-19 after treatment with various concentrations of BMP2 or BMP4. Immunofluorescence was also used to assess the changes in the expression and the organization of the key tight junction protein, zona occludens (ZO)-1, in ARPE-19 cells under BMP2 or BMP4 treatment. This was followed by measuring the activity of matrix metalloproteinases (MMPs). Finally, RNA-seq and ELISA were used to determine the local and circulating levels of BMP2 and BMP4 in retinas and serum samples from nAMD donors. Our ECIS results showed that BMP4 but not BMP2 decreased the transcellular electrical resistance (TER) of ARPE-19 and increased their migration in comparison with control (vehicle-treated cells). Furthermore, immunofluorescence showed a disorganization of ZO-1 in BMP4-treated ARPE-19 not in BMP2-treated cells or vehicle-treated controls. This effect of BMP4 was associated with significant increases in the activity of MMPs, specifically MMP2. Lastly, these results were corroborated by additional findings that circulating but not local BMP4 levels were significantly higher in nAMD donor samples compared to controls. Collectively, our results demonstrated unreported effects of BMP4 on inducing RPE dysfunction and suggest that BMP4 but not BMP2 may represent a potential therapeutic target in nAMD.

Statistical driver genes as a means to uncover missing heritability for age-related macular degeneration

Background

Age-related macular degeneration (AMD) is a progressive retinal disease contributing to blindness worldwide. Multiple estimates for AMD heritability (h²) exist; however, a substantial proportion of h² is not attributable to known genomic loci. The International AMD Genomics Consortium (IAMDGC) gathered the largest dataset of advanced AMD (ADV) cases and controls available and identified 34 loci containing 52 independent risk variants defining known AMD h². To better define AMD heterogeneity, we used Pathway Analysis by Randomization Incorporating Structure (PARIS) on the IAMDGC data and identified 8 statistical driver genes (SDGs), including 2 novel SDGs not discovered by the IAMDGC. We chose to further investigate these pathway-based risk genes and determine their contribution to ADV h², as well as the differential ADV subtype h².

Methods

We performed genomic-relatedness-based restricted maximum-likelihood (GREML) analyses on ADV, geographic atrophy (GA), and choroidal neovascularization (CNV) subtypes to investigate the h² of genotyped variants on the full DNA array chip, 34 risk loci (n = 2758 common variants), 52 variants from the IAMDGC 2016 GWAS, and the 8 SDGs, specifically the novel 2 SDGs, PPARA and PLCG2.

Results

Via GREML, full chip h² was 44.05% for ADV, 46.37% for GA, and 62.03% for CNV. The lead 52 variants’ h² (ADV: 14.52%, GA: 8.02%, CNV: 13.62%) and 34 loci h² (ADV: 13.73%, GA: 8.81%, CNV: 12.89%) indicate that known variants contribute ~ 14% to ADV h². SDG variants account for a small percentage of ADV, GA, and CNV heritability, but estimates based on the combination of SDGs and the 34 known loci are similar to those calculated for known loci alone. We identified modest epistatic interactions among variants in the 2 SDGs and the 52 IAMDGC variants, including modest interactions between variants in PPARA and PLCG2.

Conclusions

Pathway analyses, which leverage biological relationships among genes in a pathway, may be useful in identifying additional loci that contribute to the heritability of complex disorders in a non-additive manner. Heritability analyses of these loci, especially amongst disease subtypes, may provide clues to the importance of specific genes to the genetic architecture of AMD.

A multiethnic genome-wide analysis of 44,039 individuals identifies 41 new loci associated with central corneal thickness

Central corneal thickness (CCT) is one of the most heritable human traits, with broad-sense heritability estimates ranging between 0.68 to 0.95. Despite the high heritability and numerous previous association studies, only 8.5% of CCT variance is currently explained. Here, we report the results of a multiethnic meta-analysis of available genome-wide association studies in which we find association between CCT and 98 genomic loci, of which 41 are novel. Among these loci, 20 were significantly associated with keratoconus, and one (RAPSN rs3740685) was significantly associated with glaucoma after Bonferroni correction. Two-sample Mendelian randomization analysis suggests that thinner CCT does not causally increase the risk of primary open-angle glaucoma. This large CCT study explains up to 14.2% of CCT variance and increases substantially our understanding of the etiology of CCT variation. This may open new avenues of investigation into human ocular traits and their relationship to the risk of vision disorders.

Amelotin is expressed in retinal pigment epithelium and localizes to hydroxyapatite deposits in dry age-related macular degeneration

Deposition of hydroxyapatite (HAP) basal to the retinal pigment epithelium (RPE) is linked to the progression of age-related macular degeneration (AMD). Serum-deprivation of RPE cells in culture mimics some features of AMD. We now show that serum-deprivation also leads to the induction of amelotin (AMTN), a protein involved in hydroxyapatite mineralization in enamel. HAP is formed in our culture model and is blocked by siRNA inhibition of AMTN expression. In situ hybridization and immunofluorescence imaging of human eye tissue show that AMTN is expressed in RPE of donor eyes with geographic atrophy ("dry" AMD) in regions with soft drusen containing HAP spherules or nodules. AMTN is not found in hard drusen, normal RPE, or donor eyes diagnosed with wet AMD. These findings suggest that AMTN is involved in formation of HAP spherules or nodules in AMD, and as such provides a new therapeutic target for slowing disease progression.

Genome-wide association meta-analysis of corneal curvature identifies novel loci and shared genetic influences across axial length and refractive error

Corneal curvature, a highly heritable trait, is a key clinical endophenotype for myopia - a major cause of visual impairment and blindness in the world. Here we present a trans-ethnic meta-analysis of corneal curvature GWAS in 44,042 individuals of Caucasian and Asian with replication in 88,218 UK Biobank data. We identified 47 loci (of which 26 are novel), with population-specific signals as well as shared signals across ethnicities. Some identified variants showed precise scaling in corneal curvature and eye elongation (i.e. axial length) to maintain eyes in emmetropia (i.e. HDAC11/FBLN2 rs2630445, RBP3 rs11204213); others exhibited association with myopia with little pleiotropic effects on eye elongation. Implicated genes are involved in extracellular matrix organization, developmental process for body and eye, connective tissue cartilage and glycosylation protein activities. Our study provides insights into population-specific novel genes for corneal curvature, and their pleiotropic effect in regulating eye size or conferring susceptibility to myopia.

The Utah Protocol for Postmortem Eye Phenotyping and Molecular Biochemical Analysis

Purpose

Current understanding of local disease pathophysiology in AMD is limited. Analysis of the human disease-affected tissue is most informative, as gene expression, expressed quantitative trait loci, microenvironmental, and epigenetic changes can be tissue, cell type, and location specific. Development of a novel translational treatment and prevention strategies particularly for earlier forms of AMD are needed, although access to human ocular tissue analysis is challenging. We present a standardized protocol to study rapidly processed postmortem donor eyes for molecular biochemical and genomic studies.

Methods

We partnered with the Utah Lions Eye Bank to obtain donor human eyes, blood, and vitreous, within 6 hours postmortem. Phenotypic analysis was performed using spectral-domain optical coherence tomography (SD-OCT) and color fundus photography. Macular and extramacular tissues were immediately isolated, and the neural retina and retinal pigment epithelium/choroid from each specimen were separated and preserved. Ocular disease phenotype was analyzed using clinically relevant grading criteria by a group of four ophthalmologists incorporating data from SD-OCT retinal images, fundus photographs, and medical records.

Results

The use of multimodal imaging leads to greater resolution of retinal pathology, allowing greater phenotypic rigor for both interobserver phenotype and known clinical diagnoses. Further, our analysis resulted in excellent quality RNA, which demonstrated appropriate tissue segregation.

Conclusions

The Utah protocol is a standardized methodology for analysis of disease mechanisms in AMD. It uniquely allows for simultaneous rigorous phenotypic, molecular biochemical, and genomic analysis of both systemic and local tissues. This better enables the development of disease biomarkers and therapeutic interventions.

A multiethnic genome-wide association study of primary open-angle glaucoma identifies novel risk loci

Primary open-angle glaucoma (POAG) is a leading cause of irreversible vision loss, yet much of the genetic risk remains unaccounted for, especially in African-Americans who have a higher risk for developing POAG. We conduct a multiethnic genome-wide association study (GWAS) of POAG in the GERA cohort, with replication in the UK Biobank (UKB), and vice versa, GWAS in UKB with replication in GERA. We identify 24 loci (P < 5.0 × 10^-8), including 14 novel, of which 9 replicate (near FMNL2, PDE7B, TMTC2, IKZF2, CADM2, DGKG, ANKH, EXOC2, and LMX1B). Functional studies support intraocular pressure-related influences of FMNL2 and LMX1B, with certain Lmx1b mutations causing high IOP and glaucoma resembling POAG in mice. The newly identified loci increase the proportion of variance explained in each GERA race/ethnicity group, with the largest gain in African-Americans (0.5-3.1%). A meta-analysis combining GERA and UKB identifies 24 additional loci. Our study provides important insights into glaucoma pathogenesis.

Complete Transcriptome Profiling of Normal and Age-Related Macular Degeneration Eye Tissues Reveals Dysregulation of Anti-Sense Transcription

Age-related macular degeneration (AMD) predominantly affects the retina and retinal pigment epithelium in the posterior eye. While there are numerous studies investigating the non-coding transcriptome of retina and RPE, few significant differences between AMD and normal tissues have been reported. Strand specific RNA sequencing of both peripheral retina (PR) and RPE-Choroid-Sclera (PRCS), in both AMD and matched normal controls were generated. The transcriptome analysis reveals a highly significant and consistent impact on anti-sense transcription as well as moderate changes in the regulation of non-coding (sense) RNA. Hundreds of genes that do not express anti-sense transcripts in normal PR and PRCS demonstrate significant anti-sense expression in AMD in all patient samples. Several pathways are highly enriched in the upregulated anti-sense transcripts—in particular the EIF2 signaling pathway. These results call for a deeper exploration into anti-sense and noncoding RNA regulation in AMD and their potential as therapeutic targets.

A large multi-ethnic genome-wide association study identifies novel genetic loci for intraocular pressure

Elevated intraocular pressure (IOP) is a major risk factor for glaucoma, a leading cause of blindness. IOP heritability has been estimated to up to 67%, and to date only 11 IOP loci have been reported, accounting for 1.5% of IOP variability. Here, we conduct a genome-wide association study of IOP in 69,756 untreated individuals of European, Latino, Asian, and African ancestry. Multiple longitudinal IOP measurements were collected through electronic health records and, in total, 356,987 measurements were included. We identify 47 genome-wide significant IOP-associated loci (P < 5 × 10⁻⁸); of the 40 novel loci, 14 replicate at Bonferroni significance in an external genome-wide association study analysis of 37,930 individuals of European and Asian descent. We further examine their effect on the risk of glaucoma within our discovery sample. Using longitudinal IOP measurements from electronic health records improves our power to identify new variants, which together explain 3.7% of IOP variation.

Intraocular pressure (IOP) is a major risk factor for glaucoma. Here, Choquet and co-authors perform a multi-ethnic genome-wide association study of repeat IOP measurements in 69,756 individuals and identify 40 novel loci, 36 of which show directionally consistent effects in glaucoma.

Isolation and Comparative Transcriptome Analysis of Human Fetal and iPSC-Derived Cone Photoreceptor Cells

Loss of cone photoreceptors, crucial for daylight vision, has the greatest impact on sight in retinal degeneration. Transplantation of stem cell-derived L/M-opsin cones, which form 90% of the human cone population, could provide a feasible therapy to restore vision. However, transcriptomic similarities between fetal and stem cell-derived cones remain to be defined, in addition to development of cone cell purification strategies. Here, we report an analysis of the human L/M-opsin cone photoreceptor transcriptome using an AAV2/9.pR2.1:GFP reporter. This led to the identification of a cone-enriched gene signature, which we used to demonstrate similar gene expression between fetal and stem cell-derived cones. We then defined a cluster of differentiation marker combination that, when used for cell sorting, significantly enriches for cone photoreceptors from the fetal retina and stem cell-derived retinal organoids, respectively. These data may facilitate more efficient isolation of human stem cell-derived cones for use in clinical transplantation studies.

Genetic association study of exfoliation syndrome identifies a protective rare variant at LOXL1 and five new susceptibility loci

Exfoliation syndrome (XFS) is the most common known risk factor for secondary glaucoma and a major cause of blindness worldwide. Variants in two genes, LOXL1 and CACNA1A, have previously been associated with XFS. To further elucidate the genetic basis of XFS, we collected a global sample of XFS cases to refine the association at LOXL1, which previously showed inconsistent results across populations, and to identify new variants associated with XFS. We identified a rare protective allele at LOXL1 (p.Phe407, odds ratio (OR) = 25, P = 2.9 × 10^-14) through deep resequencing of XFS cases and controls from nine countries. A genome-wide association study (GWAS) of XFS cases and controls from 24 countries followed by replication in 18 countries identified seven genome-wide significant loci (P < 5 × 10^-8). We identified association signals at 13q12 (POMP), 11q23.3 (TMEM136), 6p21 (AGPAT1), 3p24 (RBMS3) and 5q23 (near SEMA6A). These findings provide biological insights into the pathology of XFS and highlight a potential role for naturally occurring rare LOXL1 variants in disease biology.

Transcriptome Analysis on Monocytes from Patients with Neovascular Age-Related Macular Degeneration

Mononuclear phagocytes (MPs), including monocytes/macrophages, play complex roles in age-related macular degeneration (AMD) pathogenesis. We reported altered gene-expression signature in peripheral blood mononuclear cells from AMD patients, and a chemokine receptor signature on AMD monocytes. To obtain comprehensive understanding of MP involvement, particularly in peripheral circulation in AMD, we performed global gene expression analysis in monocytes. We separated monocytes from treatment-naïve neovascular AMD (nvAMD) patients (n = 14) and age-matched controls (n = 15), and performed microarray and bioinformatics analysis. Quantitative real-time PCR was performed on other sets of nvAMD (n = 25), atrophic AMD (n = 21), and controls (n = 28) for validation. This validated microarray genes (like TMEM176A/B and FOSB) tested, including differences between nvAMD and atrophic AMD. We identified 2,165 differentially-expressed genes (P < 0.05), including 79 genes with log2 fold change ≥1.5 between nvAMD and controls. Functional annotation using DAVID and TANGO demonstrated immune response alterations in AMD monocytes (FDR-P <0.05), validated by randomized data comparison (P < 0.0001). GSEA, ISMARA, and MEME analysis found immune enrichment and specific involved microRNAs. Enrichment of differentially-expressed genes in monocytes was found in retina via SAGE data-mining. These genes were enriched in non-classical vs. classical monocyte subsets (P < 0.05). Therefore, global gene expression analysis in AMD monocytes reveals an altered immune-related signature, further implicating systemic MP activation in AMD.

Looking into the future: Using induced pluripotent stem cells to build two and three dimensional ocular tissue for cell therapy and disease modeling

Retinal degenerative diseases are the leading cause of irreversible vision loss in developed countries. In many cases the diseases originate in the homeostatic unit in the back of the eye that contains the retina, retinal pigment epithelium (RPE) and the choriocapillaris. RPE is a central and a critical component of this homeostatic unit, maintaining photoreceptor function and survival on the apical side and choriocapillaris health on the basal side. In diseases like age-related macular degeneration (AMD), it is thought that RPE dysfunctions cause disease-initiating events and as the RPE degenerates photoreceptors begin to die and patients start loosing vision. Patient-specific induced pluripotent stem (iPS) cell-derived RPE provides direct access to a patient's genetics and allow the possibility of identifying the initiating events of RPE-associated degenerative diseases. Furthermore, iPS cell-derived RPE cells are being tested as a potential cell replacement in disease stages with RPE atrophy. In this article we summarize the recent progress in the field of iPS cell-derived RPE "disease modeling" and cell therapies and also discuss the possibilities of developing a model of the entire homeostatic unit to aid in studying disease processes in the future. This article is part of a Special Issue entitled SI: PSC and the brain.

Meta-analysis of gene-environment-wide association scans accounting for education level identifies additional loci for refractive error

Myopia is the most common human eye disorder and it results from complex genetic and environmental causes. The rapidly increasing prevalence of myopia poses a major public health challenge. Here, the CREAM consortium performs a joint meta-analysis to test single-nucleotide polymorphism (SNP) main effects and SNP × education interaction effects on refractive error in 40,036 adults from 25 studies of European ancestry and 10,315 adults from 9 studies of Asian ancestry. In European ancestry individuals, we identify six novel loci (FAM150B-ACP1, LINC00340, FBN1, DIS3L-MAP2K1, ARID2-SNAT1 and SLC14A2) associated with refractive error. In Asian populations, three genome-wide significant loci AREG, GABRR1 and PDE10A also exhibit strong interactions with education (P<8.5 × 10(-5)), whereas the interactions are less evident in Europeans. The discovery of these loci represents an important advance in understanding how gene and environment interactions contribute to the heterogeneity of myopia.

Vision from next generation sequencing: multi-dimensional genome-wide analysis for producing gene regulatory networks underlying retinal development, aging and disease

Genomics and genetics have invaded all aspects of biology and medicine, opening uncharted territory for scientific exploration. The definition of "gene" itself has become ambiguous, and the central dogma is continuously being revised and expanded. Computational biology and computational medicine are no longer intellectual domains of the chosen few. Next generation sequencing (NGS) technology, together with novel methods of pattern recognition and network analyses, has revolutionized the way we think about fundamental biological mechanisms and cellular pathways. In this review, we discuss NGS-based genome-wide approaches that can provide deeper insights into retinal development, aging and disease pathogenesis. We first focus on gene regulatory networks (GRNs) that govern the differentiation of retinal photoreceptors and modulate adaptive response during aging. Then, we discuss NGS technology in the context of retinal disease and develop a vision for therapies based on network biology. We should emphasize that basic strategies for network construction and analyses can be transported to any tissue or cell type. We believe that specific and uniform guidelines are required for generation of genome, transcriptome and epigenome data to facilitate comparative analysis and integration of multi-dimensional data sets, and for constructing networks underlying complex biological processes. As cellular homeostasis and organismal survival are dependent on gene-gene and gene-environment interactions, we believe that network-based biology will provide the foundation for deciphering disease mechanisms and discovering novel drug targets for retinal neurodegenerative diseases.

Transcriptome of the human retina, retinal pigmented epithelium and choroid

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 high-throughput 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 high-throughput techniques.

ARHGEF12 influences the risk of glaucoma by increasing intraocular pressure

Primary open-angle glaucoma (POAG) is a blinding disease. Two important risk factors for this disease are a positive family history and elevated intraocular pressure (IOP), which is also highly heritable. Genes found to date associated with IOP and POAG are ABCA1, CAV1/CAV2, GAS7 and TMCO1. However, these genes explain only a small part of the heritability of IOP and POAG. We performed a genome-wide association study of IOP in the population-based Rotterdam Study I and Rotterdam Study II using single nucleotide polymorphisms (SNPs) imputed to 1000 Genomes. In this discovery cohort (n = 8105), we identified a new locus associated with IOP. The most significantly associated SNP was rs58073046 (β = 0.44, P-value = 1.87 × 10(-8), minor allele frequency = 0.12), within the gene ARHGEF12. Independent replication in five population-based studies (n = 7471) resulted in an effect size in the same direction that was significantly associated (β = 0.16, P-value = 0.04). The SNP was also significantly associated with POAG in two independent case-control studies [n = 1225 cases and n = 4117 controls; odds ratio (OR) = 1.53, P-value = 1.99 × 10(-8)], especially with high-tension glaucoma (OR = 1.66, P-value = 2.81 × 10(-9); for normal-tension glaucoma OR = 1.29, P-value = 4.23 × 10(-2)). ARHGEF12 plays an important role in the RhoA/RhoA kinase pathway, which has been implicated in IOP regulation. Furthermore, it binds to ABCA1 and links the ABCA1, CAV1/CAV2 and GAS7 pathway to Mendelian POAG genes (MYOC, OPTN, WDR36). In conclusion, this study identified a novel association between IOP and ARHGEF12.

Meta-analysis of Genome-Wide Association Studies Identifies Novel Loci Associated With Optic Disc Morphology

Primary open-angle glaucoma is the most common optic neuropathy and an important cause of irreversible blindness worldwide. The optic nerve head or optic disc is divided in two parts: a central cup (without nerve fibers) surrounded by the neuroretinal rim (containing axons of the retinal ganglion cells). The International Glaucoma Genetics Consortium conducted a meta-analysis of genome-wide association studies consisting of 17,248 individuals of European ancestry and 6,841 individuals of Asian ancestry. The outcomes of the genome-wide association studies were disc area and cup area. These specific measurements describe optic nerve morphology in another way than the vertical cup-disc ratio, which is a clinically used measurement, and may shed light on new glaucoma mechanisms. We identified 10 new loci associated with disc area (CDC42BPA, F5, DIRC3, RARB, ABI3BP, DCAF4L2, ELP4, TMTC2, NR2F2, and HORMAD2) and another 10 new loci associated with cup area (DHRS3, TRIB2, EFEMP1, FLNB, FAM101, DDHD1, ASB7, KPNB1, BCAS3, and TRIOBP). The new genes participate in a number of pathways and future work is likely to identify more functions related to the pathogenesis of glaucoma.

Transcriptomic analysis across nasal, temporal, and macular regions of human neural retina and RPE/choroid by RNA-Seq

Proper spatial differentiation of retinal cell types is necessary for normal human vision. Many retinal diseases, such as Best disease and male germ cell associated kinase (MAK)-associated retinitis pigmentosa, preferentially affect distinct topographic regions of the retina. While much is known about the distribution of cell types in the retina, the distribution of molecular components across the posterior pole of the eye has not been well-studied. To investigate regional difference in molecular composition of ocular tissues, we assessed differential gene expression across the temporal, macular, and nasal retina and retinal pigment epithelium (RPE)/choroid of human eyes using RNA-Seq. RNA from temporal, macular, and nasal retina and RPE/choroid from four human donor eyes was extracted, poly-A selected, fragmented, and sequenced as 100 bp read pairs. Digital read files were mapped to the human genome and analyzed for differential expression using the Tuxedo software suite. Retina and RPE/choroid samples were clearly distinguishable at the transcriptome level. Numerous transcription factors were differentially expressed between regions of the retina and RPE/choroid. Photoreceptor-specific genes were enriched in the peripheral samples, while ganglion cell and amacrine cell genes were enriched in the macula. Within the RPE/choroid, RPE-specific genes were upregulated at the periphery while endothelium associated genes were upregulated in the macula. Consistent with previous studies, BEST1 expression was lower in macular than extramacular regions. The MAK gene was expressed at lower levels in macula than in extramacular regions, but did not exhibit a significant difference between nasal and temporal retina. The regional molecular distinction is greatest between macula and periphery and decreases between different peripheral regions within a tissue. Datasets such as these can be used to prioritize candidate genes for possible involvement in retinal diseases with regional phenotypes.

Differential Gene Expression in Age-Related Macular Degeneration

Gene expression is the first step in ascribing function between an associated gene and disease. Understanding how variation in a gene influences expression, particularly in tissues affected by the disease, may help elucidate what influences the phenotypic outcome of that disease. Previous studies of the genetics of age-related macular degeneration (AMD) have identified several risk factors, but have not yet bridged the gap between gene association and identifying a specific mechanism or function that is involved in the pathogenesis of AMD. Advances in genomic technologies, such as RNA sequencing (RNA-seq), single cell RNA-seq, bilsulfite sequencing, and/or whole genome methylation, will be powerful tools for identifying genes/pathways that are differentially expressed in those with AMD versus those without AMD. These technologies should advance the field of AMD research so that appropriate preventive and therapeutic targets can be developed.

Generation of a foveomacular transcriptome

PURPOSE

Organizing molecular biologic data is a growing challenge since the rate of data accumulation is steadily increasing. Information relevant to a particular biologic query can be difficult to extract from the comprehensive databases currently available. We present a data collection and organization model designed to ameliorate these problems and applied it to generate an expressed sequence tag (EST)-based foveomacular transcriptome.

METHODS

Using Perl, MySQL, EST libraries, screening, and human foveomacular gene expression as a model system, we generated a foveomacular transcriptome database enriched for molecularly relevant data.

RESULTS

Using foveomacula as a gene expression model tissue, we identified and organized 6,056 genes expressed in that tissue. Of those identified genes, 3,480 had not been previously described as expressed in the foveomacula. Internal experimental controls as well as comparison of our data set to published data sets suggest we do not yet have a complete description of the foveomacula transcriptome.

CONCLUSIONS

We present an organizational method designed to amplify the utility of data pertinent to a specific research interest. Our method is generic enough to be applicable to a variety of conditions yet focused enough to allow for specialized study.

Drusen and photoreceptor abnormalities in African-Americans with intermediate non-neovascular age-related macular degeneration

PURPOSE/AIM

To investigate the relationship of drusen and photoreceptor abnormalities in African-American (AA) patients with intermediate non-neovascular age-related macular degeneration (AMD).

MATERIALS AND METHODS

AA patients with intermediate AMD (n = 11; age 52-77 years) were studied with spectral-domain optical coherence tomography. Macular location and characteristics of large drusen (≥125 µm) were determined. Thickness of photoreceptor laminae was quantified overlying drusen and in other macular regions. A patient with advanced AMD (age 87) was included to illustrate the disease spectrum.

RESULTS

In this AA patient cohort, the spectrum of changes known to occur in AMD, including large drusen, sub-retinal drusenoid deposits and geographic atrophy, were identified. In intermediate AMD eyes (n = 17), there were 183 large drusen, the majority of which were pericentral in location. Overlying the drusen there was significant thinning of the photoreceptor outer nuclear layer (termed ONL(+)) as well as the inner and outer segments (IS + OS). The reductions in IS + OS thickness were directly related to ONL(+) thickness. In a fraction (∼8%) of paradrusen locations with normal lamination sampled within ∼280 µm of peak drusen height, ONL(+) was significantly thickened compared to age and retinal-location-matched normal values. Topographical maps of the macula confirmed ONL thickening in regions neighboring and distant to large drusen.

CONCLUSIONS

We confirm there is a pericentral distribution of drusen across AA-AMD maculae rather than the central localization in Caucasian AMD. Reductions in the photoreceptor laminae overlying drusen are evident. ONL(+) thickening in some macular areas of AA-AMD eyes may be an early phenotypic marker for photoreceptor stress.

Comprehensive analysis of gene expression in human retina and supporting tissues

Understanding the influence of gene expression on the molecular mechanisms underpinning human phenotypic diversity is fundamental to being able to predict health outcomes and treat disease. We have carried out whole transcriptome expression analysis on a series of eight normal human postmortem eyes by RNA sequencing. Here we present data showing that ∼80% of the transcriptome is expressed in the posterior layers of the eye and that there is significant differential expression not only between the layers of the posterior part of the eye but also between locations of a tissue layer. These differences in expression also extend to alternative splicing and splicing factors. Differentially expressed genes are enriched for genes associated with psychiatric, immune and cardiovascular disorders. Enrichment categories for gene ontology included ion transport, synaptic transmission and visual and sensory perception. Lastly, allele-specific expression was found to be significant forCFH,C3 andCFB, which are known risk genes for age-related macular degeneration. These expression differences should be useful in determining the underlying biology of associations with common diseases of the human retina, retinal pigment epithelium and choroid and in guiding the analysis of the genomic regions involved in the control of normal gene expression.

Altered gene expression in dry age-related macular degeneration suggests early loss of choroidal endothelial cells

PURPOSE

Age-related macular degeneration (AMD) is a major cause of blindness in developed countries. The molecular pathogenesis of early events in AMD is poorly understood. We investigated differential gene expression in samples of human retinal pigment epithelium (RPE) and choroid from early AMD and control maculas with exon-based arrays.

METHODS

Gene expression levels in nine human donor eyes with early AMD and nine control human donor eyes were assessed using Affymetrix Human Exon ST 1.0 arrays. Two controls did not pass quality control and were removed. Differentially expressed genes were annotated using the Database for Annotation, Visualization and Integrated Discovery (DAVID), and gene set enrichment analysis (GSEA) was performed on RPE-specific and endothelium-associated gene sets. The complement factor H (CFH) genotype was also assessed, and differential expression was analyzed regarding high AMD risk (YH/HH) and low AMD risk (YY) genotypes.

RESULTS

Seventy-five genes were identified as differentially expressed (raw p value <0.01; ≥50% fold change, mean log2 expression level in AMD or control ≥ median of all average gene expression values); however, no genes were significant (adj. p value <0.01) after correction for multiple hypothesis testing. Of 52 genes with decreased expression in AMD (fold change <0.5; raw p value <0.01), 18 genes were identified by DAVID analysis as associated with vision or neurologic processes. The GSEA of the RPE-associated and endothelium-associated genes revealed a significant decrease in genes typically expressed by endothelial cells in the early AMD group compared to controls, consistent with previous histologic and proteomic studies. Analysis of the CFH genotype indicated decreased expression of ADAMTS9 in eyes with high-risk genotypes (fold change = -2.61; raw p value=0.0008).

CONCLUSIONS

GSEA results suggest that RPE transcripts are preserved or elevated in early AMD, concomitant with loss of endothelial cell marker expression. These results are consistent with the notion that choroidal endothelial cell dropout or dedifferentiation occurs early in the pathogenesis of AMD.

Gene expression profile in human trabecular meshwork from patients with primary open-angle glaucoma

PURPOSE

To identify the specific genes in human trabecular meshwork (TM) related to POAG.

METHODS

Primary open-angle glaucoma TM specimens were obtained from routine trabeculectomy surgery. Nonglaucomatous control TM specimens were dissected from donor eyes using the same approach as a standard trabeculectomy. All cases were screened for myocilin (MYOC) mutations. Total RNA was extracted, labeled, and hybridized to Illumina HumanWG-6 BeadChips. Expression data were normalized and analyzed using the R package limma in Bioconductor. Pathway analyses were performed using DAVID Bioinformatics Resources.

RESULTS

Our study included surgical TM specimens from 15 cases and 13 controls. One case was identified with a heterozygous Q368X MYOC mutation. If TMs were available from both eyes in an individual, the expression data were combined for analysis. The following three comparisons were performed for differential analyses: (1) MYOC POAG case versus 14 non-MYOC POAG cases, (2) MYOC POAG case versus 13 controls, and (3) 14 non-MYOC POAG cases versus 13 controls. Limited by one MYOC case in comparisons 1 and 2, expression changes were reported comparing the fold changes but without P values. Comparison 3 identified 483 genes, including 36 components of TM exosomes. Gene ontology analysis identified several enriched functional clusters, including cell adhesion, extracellular matrix, and secretion.

CONCLUSIONS

This is the largest TM expression study of POAG cases and controls performed to date and represents the first report of TM expression in a patient having POAG with a Q368X MYOC mutation. Our data suggest the potential role of endocytic and exosome pathways in the pathogenesis of POAG.

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

BACKGROUND

The retina is a complex tissue comprised of multiple cell types that is affected by a diverse set of diseases that are important causes of vision loss. Characterizing the transcripts, both annotated and novel, that are expressed in a given tissue has become vital for understanding the mechanisms underlying the pathology of disease.

RESULTS

We sequenced RNA prepared from three normal human retinas and characterized the retinal transcriptome at an unprecedented level due to the increased depth of sampling provided by the RNA-seq approach. We used a non-redundant reference transcriptome from all of the empirically-determined human reference tracks to identify annotated and novel sequences expressed in the retina. We detected 79,915 novel alternative splicing events, including 29,887 novel exons, 21,757 3' and 5' alternate splice sites, and 28,271 exon skipping events. We also identified 116 potential novel genes. These data represent a significant addition to the annotated human transcriptome. For example, the novel exons detected increase the number of identified exons by 3%. Using a high-throughput RNA capture approach to validate 14,696 of these novel transcriptome features we found that 99% of the putative novel events can be reproducibly detected. Further, 15-36% of the novel splicing events maintain an open reading frame, suggesting they produce novel protein products.

CONCLUSIONS

To our knowledge, this is the first application of RNA capture to perform large-scale validation of novel transcriptome features. In total, these analyses provide extensive detail about a previously uncharacterized level of transcript diversity in the human retina.

A Prevalence of Imprinted Genes within the Total Transcriptomes of Human Tissues and Cells

Genomic imprinting is an epigenetic phenomenon that causes a differential expression of paternally and maternally inherited alleles of a subset of genes (the so-called imprinted genes). Imprinted genes are distributed throughout the genome and it is predicted that about 1% of the human genes may be imprinted. It is recognized that the allelic expression of imprinted genes varies between tissues and developmental stages. The current study represents the first attempt to estimate a prevalence of imprinted genes within the total human transcriptome. In silico analysis of the normalized expression profiles of a comprehensive panel of 173 established and candidate human imprinted genes was performed, in 492 publicly available SAGE libraries. The latter represent human cell and tissue samples in a variety of physiological and pathological conditions. Variations in the prevalence of imprinted genes within the total transcriptomes (ranging from 0.08% to 4.36%) and expression profiles of the individual imprinted genes are assessed. This paper thus provides a useful reference on the size of the imprinted transcriptome and expression of the individual imprinted genes.

Human RNA integrity after postmortem retinal tissue recovery

PURPOSE

To assess the parameters for postmortem retinal tissue recovery and processing that affect the quality of RNA extracted from the retina/retinal pigment epithelium (RPE) complex.

METHODS

RNA was extracted from retina/RPE samples. The RNA quality was determined based on qualitative/quantitative measurements made with a Bioanalyzer (Agilent) and on the expression of a long retinal gene (RPE65). After a pilot analysis on rats, ocular RNA was extracted from human donor eyeballs (group A) explanted according to conventional procedures for cornea transplantation. In a second experiment, another group of human donor eyeballs (group B) were processed in a much shorter time. The postmortem interval (T) comprised two periods: T1, the time between a donor's death and enucleation, and T2, the time between eyeball explantation and immersion of the excised retina/RPE sample in preservative solution (T = T1 + T2).

RESULTS

A short T2 was correlated with good quality of RNA extracted from the retina/RPE complex (p = 0.043) and successful expression of a tissue-specific gene (p = 0.007). No other parameter appeared to influence RNA quality.

CONCLUSIONS

The time between eyeball explantation and immersion of the retina/RPE sample in preservative solution was the chief parameter affecting the quality of RNA extracted from the retina/RPE complex.

Screening of candidate genes for primary open angle glaucoma

PURPOSE

Primary open-angle glaucoma (POAG) is one of the leading causes of irreversible blindness in the world. To make progress in understanding POAG, it is necessary to identify more POAG-causing genes.

METHODS

Using haplotype analysis, we found that mutational region is located on chromosome 2 in two families. Furthermore, we screened 11 candidate genes on chromosome 2 by protein-protein interaction (PPI) analysis, including mutS homolog 6 (MSH6), mutS homolog 2 (MSH2), v-rel reticuloendotheliosis viral oncogene homolog (REL), endothelial PAS domain protein 1 (EPAS1), vaccinia related kinase 2 (VRK2), F-box protein 11 (FBXO11), EGF containing fibulin-like extracellular matrix protein 1 (EFEMP1), reticulon 4 (RTN4), RAB1A, member RAS oncogene family (RAB1A), ARP2 actin-related protein 2 homolog (ACTR2), and calmodulin 2 (phosphorylase kinase, delta; CALM2). These 11 genes are all predicted to be related to trabecular meshwork changes and progressive loss of retinal ganglion cells in POAG patients.

RESULTS

According to our study, FBXO11 and VRK2 may interact with tumor protein p53 to regulate mitochondrial membrane permeability, mitochondrial membrane organization, and apoptosis. MSH2 is responsible for repairing DNA mismatches and RTN4 is for neuronal regeneration. Therefore, they are supposed to play a negative role in cellular process in POAG. CALM2 may be involved in retinal ganglion cell death and oxidative damage to cell communication.

CONCLUSIONS

The results demonstrate that the genes above may be associated with pathogenesis of POAG.

Transcriptome changes in age-related macular degeneration

Age-related macular degeneration (AMD) is a debilitating, common cause of visual impairment. While the last decade has seen great progress in understanding the pathophysiology of AMD, the molecular changes that occur in eyes with AMD are still poorly understood. In the current issue of Genome Medicine, Newman and colleagues present the first systematic transcriptional profile analysis of AMD-affected tissues, providing a comprehensive set of expression data for different regions (macula versus periphery), tissues (retina versus retinal pigment epithelium (RPE)/choroid), and disease state (control versus early or advanced AMD). Their findings will serve as a foundation for additional systems-level research into the pathogenesis of this blinding disease.

Chemistry and biology of vision

Visual perception in humans occurs through absorption of electromagnetic radiation from 400 to 780 nm by photoreceptors in the retina. A photon of visible light carries a sufficient amount of energy to cause, when absorbed, a cis,trans-geometric isomerization of the 11-cis-retinal chromophore, a vitamin A derivative bound to rhodopsin and cone opsins of retinal photoreceptors. The unique biochemistry of these complexes allows us to reliably and reproducibly collect continuous visual information about our environment. Moreover, other nonconventional retinal opsins such as the circadian rhythm regulator melanopsin also initiate light-activated signaling based on similar photochemistry.

A dominant mutation in RPE65 identified by whole-exome sequencing causes retinitis pigmentosa with choroidal involvement

Linkage testing using Affymetrix 6.0 SNP Arrays mapped the disease locus in TCD-G, an Irish family with autosomal dominant retinitis pigmentosa (adRP), to an 8.8 Mb region on 1p31. Of 50 known genes in the region, 11 candidates, including RPE65 and PDE4B, were sequenced using di-deoxy capillary electrophoresis. Simultaneously, a subset of family members was analyzed using Agilent SureSelect All Exome capture, followed by sequencing on an Illumina GAIIx platform. Candidate gene and exome sequencing resulted in the identification of an Asp477Gly mutation in exon 13 of the RPE65 gene tracking with the disease in TCD-G. All coding exons of genes not sequenced to sufficient depth by next generation sequencing were sequenced by di-deoxy sequencing. No other potential disease-causing variants were found to segregate with disease in TCD-G. The Asp477Gly mutation was not present in Irish controls, but was found in a second Irish family provisionally diagnosed with choroideremia, bringing the combined maximum two-point LOD score to 5.3. Mutations in RPE65 are a known cause of recessive Leber congenital amaurosis (LCA) and recessive RP, but no dominant mutations have been reported. Protein modeling suggests that the Asp477Gly mutation may destabilize protein folding, and mutant RPE65 protein migrates marginally faster on SDS-PAGE, compared with wild type. Gene therapy for LCA patients with RPE65 mutations has shown great promise, raising the possibility of related therapies for dominant-acting mutations in this gene.

Serial analysis of gene expression (SAGE) in normal human trabecular meshwork

PURPOSE

To identify the genes expressed in normal human trabecular meshwork tissue, a tissue critical to the pathogenesis of glaucoma.

METHODS

Total RNA was extracted from human trabecular meshwork (HTM) harvested from 3 different donors. Extracted RNA was used to synthesize individual SAGE (serial analysis of gene expression) libraries using the I-SAGE Long kit from Invitrogen. Libraries were analyzed using SAGE 2000 software to extract the 17 base pair sequence tags. The extracted sequence tags were mapped to the genome using SAGE Genie map.

RESULTS

A total of 298,834 SAGE tags were identified from all HTM libraries (96,842, 88,126, and 113,866 tags, respectively). Collectively, there were 107,325 unique tags. There were 10,329 unique tags with a minimum of 2 counts from a single library. These tags were mapped to known unique Unigene clusters. Approximately 29% of the tags (orphan tags) did not map to a known Unigene cluster. Thirteen percent of the tags mapped to at least 2 Unigene clusters. Sequence tags from many glaucoma-related genes, including myocilin, optineurin, and WD repeat domain 36, were identified.

CONCLUSIONS

This is the first time SAGE analysis has been used to characterize the gene expression profile in normal HTM. SAGE analysis provides an unbiased sampling of gene expression of the target tissue. These data will provide new and valuable information to improve understanding of the biology of human aqueous outflow.

Dissection of chromosome 16p12 linkage peak suggests a possible role for CACNG3 variants in age-related macular degeneration susceptibility

PURPOSE

Age-related macular degeneration (AMD) is a complex disorder of the retina, characterized by drusen, geographic atrophy, and choroidal neovascularization. Cigarette smoking and the genetic variants CFH Y402H, ARMS2 A69S, CFB R32Q, and C3 R102G have been strongly and consistently associated with AMD. Multiple linkage studies have found evidence suggestive of another AMD locus on chromosome 16p12 but the gene responsible has yet to be identified.

METHODS

In the initial phase of the study, single-nucleotide polymorphisms (SNPs) across chromosome 16 were examined for linkage and/or association in 575 Caucasian individuals from 148 multiplex and 77 singleton families. Additional variants were tested in an independent dataset of unrelated cases and controls. According to these results, in combination with gene expression data and biological knowledge, five genes were selected for further study: CACNG3, HS3ST4, IL4R, Q7Z6F8, and ITGAM.

RESULTS

After genotyping additional tagging SNPs across each gene, the strongest evidence for linkage and association was found within CACNG3 (rs757200 nonparametric LOD* = 3.3, APL (association in the presence of linkage) P = 0.06, and rs2238498 MQLS (modified quasi-likelihood score) P = 0.006 in the families; rs2283550 P = 1.3 × 10(-6), and rs4787924 P = 0.002 in the case-control dataset). After adjusting for known AMD risk factors, rs2283550 remained strongly associated (P = 2.4 × 10(-4)). Furthermore, the association signal at rs4787924 was replicated in an independent dataset (P = 0.035) and in a joint analysis of all the data (P = 0.001).

CONCLUSIONS

These results suggest that CACNG3 is the best candidate for an AMD risk gene within the 16p12 linkage peak. More studies are needed to confirm this association and clarify the role of the gene in AMD pathogenesis.

Building the developmental oculome: systems biology in vertebrate eye development and disease

The vertebrate eye is a sophisticated multicomponent organ that has been actively studied for over a century, resulting in the identification of the major embryonic and molecular events involved in its complex developmental program. Data gathered so far provides sufficient information to construct a rudimentary network of the various signaling molecules, transcription factors, and their targets for several key stages of this process. With the advent of genomic technologies, there has been a rapid expansion in our ability to collect and process biological information, and the use of systems-level approaches to study specific aspects of vertebrate eye development has already commenced. This is beginning to result in the definition of the dynamic developmental networks that operate in ocular tissues, and the interactions of such networks between coordinately developing ocular tissues. Such an integrative understanding of the eye by a comprehensive systems-level analysis can be termed the 'oculome', and that of serial developmental stages of the eye as it transits from its initiation to a fully formed functional organ represents the 'developmental oculome'. Construction of the developmental oculome will allow novel mechanistic insights that are essential for organ regeneration-based therapeutic applications, and the generation of computational models for eye disease states to predict the effects of drugs. This review discusses our present understanding of two of the individual components of the developing vertebrate eye--the lens and retina--at both the molecular and systems levels, and outlines the directions and tools required for construction of the developmental oculome.

Identification of disease-causing mutations in autosomal dominant retinitis pigmentosa (adRP) using next-generation DNA sequencing

PURPOSE

To determine whether massively parallel next-generation DNA sequencing offers rapid and efficient detection of disease-causing mutations in patients with monogenic inherited diseases. Retinitis pigmentosa (RP) is a challenging application for this technology because it is a monogenic disease in individuals and families but is highly heterogeneous in patient populations. RP has multiple patterns of inheritance, with mutations in many genes for each inheritance pattern and numerous, distinct, disease-causing mutations at each locus; further, many RP genes have not been identified yet.

METHODS

Next-generation sequencing was used to identify mutations in pairs of affected individuals from 21 families with autosomal dominant RP, selected from a cohort of families without mutations in "common" RP genes. One thousand amplicons targeting 249,267 unique bases of 46 candidate genes were sequenced with the 454GS FLX Titanium (Roche Diagnostics, Indianapolis, IN) and GAIIx (Illumina/Solexa, San Diego, CA) platforms.

RESULTS

An average sequence depth of 70× and 125× was obtained for the 454GS FLX and GAIIx platforms, respectively. More than 9000 sequence variants were identified and analyzed, to assess the likelihood of pathogenicity. One hundred twelve of these were selected as likely candidates and tested for segregation with traditional di-deoxy capillary electrophoresis sequencing of additional family members and control subjects. Five disease-causing mutations (24%) were identified in the 21 families.

CONCLUSION

This project demonstrates that next-generation sequencing is an effective approach for detecting novel, rare mutations causing heterogeneous monogenic disorders such as RP. With the addition of this technology, disease-causing mutations can now be identified in 65% of autosomal dominant RP cases.

A new strategy to identify and annotate human RPE-specific gene expression

Background

To identify and functionally annotate cell type-specific gene expression in the human retinal pigment epithelium (RPE), a key tissue involved in age-related macular degeneration and retinitis pigmentosa.

Methodology

RPE, photoreceptor and choroidal cells were isolated from selected freshly frozen healthy human donor eyes using laser microdissection. RNA isolation, amplification and hybridization to 44 k microarrays was carried out according to Agilent specifications. Bioinformatics was carried out using Rosetta Resolver, David and Ingenuity software.

Principal Findings

Our previous 22 k analysis of the RPE transcriptome showed that the RPE has high levels of protein synthesis, strong energy demands, is exposed to high levels of oxidative stress and a variable degree of inflammation. We currently use a complementary new strategy aimed at the identification and functional annotation of RPE-specific expressed transcripts. This strategy takes advantage of the multilayered cellular structure of the retina and overcomes a number of limitations of previous studies. In triplicate, we compared the transcriptomes of RPE, photoreceptor and choroidal cells and we deduced RPE specific expression. We identified at least 114 entries with RPE-specific gene expression. Thirty-nine of these 114 genes also show high expression in the RPE, comparison with the literature showed that 85% of these 39 were previously identified to be expressed in the RPE. In the group of 114 RPE specific genes there was an overrepresentation of genes involved in (membrane) transport, vision and ophthalmic disease. More fundamentally, we found RPE-specific involvement in the RAR-activation, retinol metabolism and GABA receptor signaling pathways.

Conclusions

In this study we provide a further specification and understanding of the RPE transcriptome by identifying and analyzing genes that are specifically expressed in the RPE.

Targeted high-throughput DNA sequencing for gene discovery in retinitis pigmentosa

The causes of retinitis pigmentosa (RP) are highly heterogeneous, with mutations in more than 60 genes known to cause syndromic and non-syndromic forms of disease. The prevalence of detectable mutations in known genes ranges from 25 to 85%, depending on mode of inheritance. For example, the likelihood of detecting a disease-causing mutation in known genes in patients with autosomal dominant RP (adRP) is 60% in Americans and less in other populations. Thus many RP genes are still unknown or mutations lie outside of commonly tested regions. Furthermore, current screening strategies can be costly and time-consuming.We are developing targeted high-throughput DNA sequencing to address these problems. In this approach, a microarray with oligonucleotides targeted to hundreds of genes is used to capture sheared human DNA, and the sequence of the eluted DNA is determined by ultra-high-throughput sequencing using next-generation DNA sequencing technology. The first capture array we have designed contains 62 full-length retinal disease genes, including introns and promoter regions, and an additional 531 genes limited to exons and flanking sequences. The full-length genes include all genes known to cause at least 1% of RP or other inherited retinal diseases. All of the genes listed in the RetNet database are included on the capture array as well as many additional retinal-expressed genes. After validation studies, the first DNA's tested will be from 89 unrelated adRP families in which the prevalent RP genes have been excluded. This approach should identify new RP genes and will substantially reduce the cost per patient.

Differential gene expression in the developing human macula: microarray analysis using rare tissue samples

The macula is a unique and important region in the primate retina that achieves high resolution and color vision in the central visual field. We recently reported data obtained from microarray analysis of gene expression in the macula of the human fetal retina (Kozulin et al., Mol Vis 15:45–59, 1). In this paper, we describe the preliminary analyses undertaken to visualize differences and verify comparability of the replicates used in that study, report the differential expression of other gene families obtained from the analysis, and show the reproducibility of our findings in several gene families by quantitative real-time PCR.

Time-dependent reduction of glutamine synthetase in retina of diabetic rats

Diabetic Retinopathy (DR) is one of the most common complications of diabetes and a major cause of blindness worldwide. We studied the transcriptome of the diabetic retina using Series Analysis of Gene Expression (SAGE) technology and observed a 45.6% reduction in transcript levels of glutamine synthetase (GS) in streptozotocin-induced diabetic rats compared with normal rats. RT-PCR and colorimetric enzyme activity assays revealed significant differences in GS mRNA expression and enzyme activity as early as the first month of diabetes development, with a progressive decrease in GS mRNA level and enzyme activity over a 12-month period. Northern blot analysis indicated a linear correlation between the reduction in GS expression and the time course of diabetic retinopathy (r = 0.802, p < 0.0001), which was validated by real-time RT-PCR (r = 0.731, p < 0.001). Our results implicate GS as a possible biomarker for evaluating the severity of developed diabetic retinopathy over the time course of diabetes progression.

Application of serial analysis of gene expression to the study of human genetic disease

Sequence tag analysis using serial analysis of gene expression (SAGE) is a powerful strategy for the quantitative analysis of gene expression in human genetic disorders. SAGE facilitates the measurement of mRNA transcripts and generates a non-biased gene expression profile of normal and pathological disease tissue. In addition, the SAGE technique has the capacity of detecting the expression of novel transcripts allowing for the identification of previously uncharacterised genes, thus providing a unique advantage over the traditional microarray-based approach for expression profiling. The technique has been successful in providing pathological gene expression profiles in a number of common genetic disorders including diabetes, cardiovascular disease, Parkinson disease and Down syndrome. When combined with next generation sequencing platforms, SAGE has the potential to become a more powerful and sensitive technique making it more amenable for diagnostic use. This review will therefore discuss the application of SAGE to several common genetic disorders and will further evaluate its potential use in diagnosing human genetic disease.

Expression of the diabetes risk gene wolframin (WFS1) in the human retina

Wolfram syndrome 1 (WFS1, OMIM 222300), a rare genetic disorder characterized by optic nerve atrophy, deafness, diabetes insipidus and diabetes mellitus, is caused by mutations of WFS1, encoding WFS1/wolframin. Non-syndromic WFS1 variants are associated with the risk of diabetes mellitus due to altered function of wolframin in pancreatic islet cells, expanding the importance of wolframin. This study extends a previous report for the monkey retina, using immunohistochemistry to localize wolframin on cryostat and paraffin sections of human retina. In addition, the human retinal pigment epithelial (RPE) cell line termed ARPE-19 and retinas from both pigmented and albino mice were studied to assess wolframin localization. In the human retina, wolframin was expressed in retinal ganglion cells, optic axons and the proximal optic nerve. Wolframin expression in the human retinal pigment epithelium (RPE) was confirmed with intense cytoplasmic labeling in ARPE-19 cells. Strong labeling of the RPE was also found in the albino mouse retina. Cryostat sections of the mouse retina showed a more extended pattern of wolframin labeling, including the inner nuclear layer (INL) and photoreceptor inner segments, confirming the recent report of Kawano et al. [Kawano, J., Tanizawa, Y., Shinoda, K., 2008. Wolfram syndrome 1 (Wfs1) gene expression in the normal mouse visual system. J. Comp. Neurol. 510, 1-23]. Absence of these cells in the human specimens despite the use of human-specific antibodies to wolframin may be related to delayed fixation. Loss of wolframin function in RGCs and the unmyelinated portion of retinal axons could explain optic nerve atrophy in Wolfram Syndrome 1.

Genomewide linkage scans for ocular refraction and meta-analysis of four populations in the Myopia Family Study

PURPOSE

Genomewide linkage scans were performed in Caucasian (CAUC) and Old Order Amish (OOA) families to identify genomic regions containing genes responsible for refractive error control. We also performed a meta-analysis by combining these results with our previous linkage results from Ashkenazi Jewish (ASHK) and African American (AFRAM) families.

METHODS

Two hundred seventy-one CAUC and 411 OOA participants (36 and 61 families, respectively) were recruited to participate in the Myopia Family Study. Recruitment criteria were designed to enrich the sample for multiplex myopic families. Genomewide, model-free, multipoint linkage analyses were performed separately for each population by using >370 microsatellite markers. Empirical significance levels were determined via gene-dropping simulations. A meta-analysis was performed by combining linkage results from the CAUC, OOA, AFRAM, and ASHK samples, and results were compared to previously reported loci for myopia and refraction.

RESULTS

Suggestive evidence of linkage was found at 12q24 (LOD = 4.583, P = 0.00037) and 4q21 (LOD = 2.72, P = 0.0028) in the CAUC sample and at 5qter (LOD = 3.271, P = 0.0014) in the OOA. Meta-analysis linkage results were largely driven by population-specific signals from ASHK and AFRAM families. The meta-analysis showed suggestive evidence of linkage to 4q21-22 (meta-P = 0.00214) adjacent to the previously reported MYP9 and MYP11 loci.

CONCLUSIONS

The results showed suggestive evidence of linkage of ocular refraction to 12q24 and 4q21 in CAUC and to 5qter in OOA families. The meta-analysis supports the view that several genes play a role in refractive development across populations. In MFS families, four broad genomic regions (on 1p, 4q, 7p, and 12q) most likely contain genes that influence ocular refraction.