Introduction: Despite druggable events to be present in 80% of neuroblastomapatients within the Princess Maxima Center precision medicine program 'iTHER', clinical uptake of treatment recommendations has been low, and the clinical impact for individual patients remains hard to predict. This stresses the need for a method integrating genomics and transcriptomics with functional approaches into therapeutic decision making. Methods: We aimed to launch an online repository integrating genomics and transcriptomics with high-throughput drug screening (HTS) of nineteen commonly used neuroblastoma cell lines and fifteen neuroblastoma patient-derived organoids (NBL-PDOs). Cell lines, NBL-PDOs and their parental tumors were characterized utilizing (lc)WGS, WES and RNAseq. Cells were exposed to similar to 200 compounds. Results were transferred to the R2 visualization platform. Results: A powerful reference set of cell lines is available, reflecting distinct known pharmacologic vulnerabilities. HTS identified additional therapeutic vulnerabilities, such as a striking correlation between a positive mesenchymal signature and sensitivity to BCL2-inhibitor venetoclax. Finally, we explored personalized drug sensitivities within iTHER, demonstrating HTS can support genomic and transcriptomic results, thereby strengthening the rationale for clinical uptake. Conclusion: We established a dynamic publicly available dataset with detailed genomic, transcriptomic, and pharmacological annotation of classical neuroblastoma cell lines as well as novel sharable NBL-PDOs, representing the heterogeneous landscape of neuroblastoma. We anticipate that in vitro drug screening will be complementary to genomic-guided precision medicine by supporting clinical decision making, thereby improving prognosis for all neuroblastoma patients in the future.

BACKGROUND:Pediatric cancer is the leading cause of disease-related death in children and the need for better therapeutic options remains urgent. Due to the limited number of patients, target and drug development for pediatrics is often supplemented by data from studies focused on adult cancers. Recent evidence shows that pediatric cancers possess different vulnerabilities that should be explored independently from adult cancers.METHODS:Using the publicly available Genomics of Drug Sensitivity in Cancer database, we explore therapeutic targets and biomarkers specific to the pediatric solid malignancies Ewing sarcoma, medulloblastoma, neuroblastoma, osteosarcoma, and rhabdomyosarcoma. Results are validated using cell viability assays and high-throughput drug screens are used to identify synergistic combinations.RESULTS:Using published drug screening data, PARP is identified as a drug target of interest across multiple different pediatric malignancies. We validate these findings, and we show that efficacy can be improved when combined with conventional chemotherapeutics, namely topoisomerase inhibitors. Additionally, using gene set enrichment analysis, we identify ribosome biogenesis as a potential biomarker for PARP inhibition in pediatric cancer cell lines.CONCLUSION:Collectively, our results provide evidence to support the further development of PARP inhibition and the combination with TOP1 inhibition as a therapeutic approach in solid pediatric malignancies. Additionally, we propose ribosome biogenesis as a component to PARP inhibitor sensitivity that should be further investigated to help maximize the potential utility of PARP inhibition and combinations across pediatric solid malignancies.

BACKGROUND:Immunotherapy in high-risk neuroblastoma (HR-NBL) does not live up to its full potential due to inadequate (adaptive) immune engagement caused by the extensive immunomodulatory capacity of HR-NBL. We aimed to tackle one of the most notable immunomodulatory processes in neuroblastoma (NBL), absence of major histocompatibility complex class I (MHC-I) surface expression, a process greatly limiting cytotoxic T cell engagement. We and others have previously shown that MHC-I expression can be induced by cytokine-driven immune modulation. Here, we aimed to identify tolerable pharmacological repurposing strategies to upregulate MHC-I expression and therewith enhance T cell immunogenicity in NBL.METHODS:Drug repurposing libraries were screened to identify compounds enhancing MHC-I surface expression in NBL cells using high-throughput flow cytometry analyses optimized for adherent cells. The effect of positive hits was confirmed in a panel of NBL cell lines and patient-derived organoids. Compound-treated NBL cell lines and organoids were cocultured with preferentially expressed antigen of melanoma (PRAME)-reactive tumor-specific T cells and healthy-donor natural killer (NK) cells to determine the in vitro effect on T cell and NK cell cytotoxicity. Additional immunomodulatory effects of histone deacetylase inhibitors (HDACi) were identified by transcriptome and translatome analysis of treated organoids.RESULTS:Drug library screening revealed MHC-I upregulation by inhibitor of apoptosis inhibitor (IAPi)- and HDACi drug classes. The effect of IAPi was limited due to repression of nuclear factor kappa B (NFκB) pathway activity in NBL, while the MHC-I-modulating effect of HDACi was widely translatable to a panel of NBL cell lines and patient-derived organoids. Pretreatment of NBL cells with the HDACi entinostat enhanced the cytotoxic capacity of tumor-specific T cells against NBL in vitro, which coincided with increased expression of additional players regulating T cell cytotoxicity (eg, TAP1/2 and immunoproteasome subunits). Moreover, MICA and MICB, important in NK cell cytotoxicity, were also increased by entinostat exposure. Intriguingly, this increase in immunogenicity was accompanied by a shift toward a more mesenchymal NBL cell lineage.CONCLUSIONS:This study indicates the potential of combining (immuno)therapy with HDACi to enhance both T cell-driven and NKcell-driven immune responses in patients with HR-NBL.

ABSTRACTHomozygous inactivation of the CDKN2A locus is one of the most common genomic aberrations in human cancer. The locus codes for two unrelated and distinctly regulated proteins: p14ARF and p16INK4a, which inhibit MDM2 and CDK4/6, respectively. Loss of CDKN2A is also a recurrent event in relapsed neuroblastoma, a childhood tumour that arises from neural crest cells. To examine the consequences of the loss of the two distinct gene transcripts in neuroblastoma, we used the CRISPR-Cas9 system to knockout p14, p16 and p14+p16 in SY5Y cells. RNA sequencing of the transcriptome revealed a striking shift towards an immature Schwann cell precursor-like phenotype with mesenchymal characteristics, specifically in the p16 and p14+p16 knockouts. High-throughput drug screening of p16 and p14+p16 knockout clones identified a large increase in sensitivity to EGFR inhibitors. On protein level, we were able to confirm that EGFR pathway activation is higher in p14+p16 knockout cells and that treatment with the EGFR inhibitor afatinib resulted in higher levels of apoptosis. Afatinib also reduced tumour growth in vivo in xenografts transplanted with p14+p16 knockout SY5Y cells. Overall, our study suggests that CDKN2A deletion in neuroblastoma relates to a phenotypic shift towards a more progenitor like state and increases sensitivity to EGFR inhibitors.

Neuroblastomas are childhood tumors with frequent fatal relapses after induction treatment, which is related to tumor evolution with additional genomic events. Our whole-genome sequencing data analysis revealed a high frequency of somatic cytosine > adenine (C > A) substitutions in primary neuroblastoma tumors, which was associated with poor survival. We showed that increased levels of C > A substitutions correlate with copy number loss (CNL) of OGG1 or MUTYH. Both genes encode DNA glycosylases that recognize 8-oxo-guanine (8-oxoG) lesions as a first step of 8-oxoG repair. Tumor organoid models with CNL of OGG1 or MUTYH show increased 8-oxoG levels compared to wild-type cells. We used CRISPR-Cas9 genome editing to create knockout clones of MUTYH and OGG1 in neuroblastoma cells. Whole-genome sequencing of single-cell OGG1 and MUTYH knockout clones identified an increased accumulation of C > A substitutions. Mutational signature analysis of these OGG1 and MUTYH knockout clones revealed enrichment for C > A signatures 18 and 36, respectively. Clustering analysis showed that the knockout clones group together with tumors containing OGG1 or MUTYH CNL. In conclusion, we demonstrate that defects in 8-oxoG repair cause accumulation of C > A substitutions in neuroblastoma, which contributes to mutagenesis and tumor evolution.

MLL-rearranged acute lymphoblastic leukemia (ALL) represents a highly aggressive ALL subtype, characterized by aberrant DNA methylation patterns. DNA methyltransferase inhibitors, such as decitabine have previously been demonstrated to be effective in eradicating MLL-rearranged ALL cells in vitro. Here, we assessed the in vivo anti-leukemic potential of low-dose DNA methyltransferase inhibitor decitabine using a xenograft mouse model of human MLL-rearranged ALL. Furthermore, we explored whether prolonged exposure to low-dose decitabine could chemo-sensitize MLL-rearranged ALL cells toward conventional chemotherapy as well as other known epigenetic-based and anti-neoplastic compounds. Our data reveal that decitabine prolonged survival in xenograft mice of MLL-rearranged ALL by 8.5 days (P = .0181), but eventually was insufficient to prevent leukemia out-growth, based on the examination of the MLLAF4 cell line SEM. Furthermore, we observe that prolonged pretreatment of low-dose decitabine mildly sensitized toward the conventional drugs prednisolone, vincristine, daunorubicin, asparaginase, and cytarabine in a panel of MLL-rearranged cell lines. Additionally, we assessed synergistic effects of decitabine with other epigenetic-based or anticancer drugs using high-throughput drug library screens. Validation of the top hits, including histone deacetylase inhibitor panobinostat, BCL2 inhibitor Venetoclax, MEK inhibitor pimasertib, and receptor tyrosine kinase foretinib, revealed additive and moderate synergistic effects for the combination of each drug together with decitabine in a cell line-dependent manner.

Chromosome 17q gains are almost invariably present in high-risk neuroblastoma cases. Here, we perform an integrative epigenomics search for dosage-sensitive transcription factors on 17q marked by H3K27ac defined super-enhancers and identify TBX2 as top candidate gene. We show that TBX2 is a constituent of the recently established core regulatory circuitry in neuroblastoma with features of a cell identity transcription factor, driving proliferation through activation of p21-DREAM repressed FOXM1 target genes. Combined MYCN/TBX2 knockdown enforces cell growth arrest suggesting that TBX2 enhances MYCN sustained activation of FOXM1 targets. Targeting transcriptional addiction by combined CDK7 and BET bromodomain inhibition shows synergistic effects on cell viability with strong repressive effects on CRC gene expression and p53 pathway response as well as several genes implicated in transcriptional regulation. In conclusion, we provide insight into the role of the TBX2 CRC gene in transcriptional dependency of neuroblastoma cells warranting clinical trials using BET and CDK7 inhibitors.

Abstract BACKGROUND: The anti-apoptotic protein B cell lymphoma/leukaemia 2 (BCL-2) is highly expressed in the majority of all neuroblastomas. In previous preclinical studies, we have shown that treatment of BCL-2-dependent neuroblastoma with BCL-2 inhibitors leads to programmed cell death. These results have contributed to the initiation of a phase I trial to study the safety and pharmacokinetics of Venetoclax in children with relapsed or refractory neuroblastoma. AIM: The current study aims to identify and validate targeted combination strategies to prevent or overcome neuroblastoma resistance to Venetoclax. METHODS: Targeted drug candidates for combination treatment with Venetoclax were identified by high-throughput drug screening of non-resistant and Venetoclax-resistant BCL-2-dependent neuroblastoma cell lines. Top hits were subsequently tested more extensively in vitro and in vivo to validate the screening results. RESULTS: High-throughput drug screens identified the MDM2 inhibitor Idasanutlin as one of the strongest re-sensitizers for Venetoclax in Venetoclax-resistant BCL-2-dependent neuroblastoma cells with wild-type p53. Subsequent in vitro validation showed that Idasanutlin induced cyclin dependent kinase inhibitor 1 (p21)-mediated growth arrest in non-resistant neuroblastoma cells, but induced a BCL-2-associated X protein (BAX)-mediated apoptotic response in Venetoclax-resistant neuroblastoma cells in the presence of Venetoclax. In vivo combination of Venetoclax with Idasanultin resulted in a remarkably improved anticancer effect compared to single agent therapy, with very good partial and complete responses in BCL-2-dependent neuroblastoma xenografts. Combination screens additionally revealed that the clinically studied targeted inhibitors CUDC-907 (PI3K/HDAC), Prexasertib (CHK1) and Talazoparib (PARP1/2) improved the efficacy of Venetoclax more potently than ALK inhibitors in BCL-2-dependent neuroblastoma cells harboring ALK mutations. CONCLUSION: Our findings suggest that the clinical use of Venetoclax for the treatment of children with BCL-2-dependent neuroblastoma tumors can be improved by combination therapy with targeted inhibitors. The presence of additional neuroblastoma driving genomic events is not always predictive of the optimal combination strategy for BCL-2-dependent neuroblastoma subgroups. These findings should guide the design of combination Phase 2 trials of Venetoclax with other targeted compounds such as Idasanutlin, CUDC-907, Prexasertib and Talazoparib. Citation Format: M. Emmy M. Dolman, Laurel T. Bate-Eya, Lindy Vernooij, Bianca Koopmans, Lindy K. Alles, Anke H. Essing, Daphne Lelieveld, David A. Egan, Mark Kerstjens, Ronald W. Stam, Huib N. Caron, Jan J. Molenaar. Preclinical identification of Venetoclax combination strategies with Idasanutlin, CUDC-907, Prexasertib or Talazoparib for neuroblastoma treatment [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2630.

Olfactory ensheathing cells (OECs) have neuro-restorative properties in animal models for spinal cord injury, stroke, and amyotrophic lateral sclerosis. Here we used a multistep screening approach to discover genes specifically contributing to the regeneration-promoting properties of OECs. Microarray screening of the injured olfactory pathway and of cultured OECs identified 102 genes that were subsequently functionally characterized in cocultures of OECs and primary dorsal root ganglion (DRG) neurons. Selective siRNA-mediated knockdown of 16 genes in OECs (ADAMTS1, BM385941, FZD1, GFRA1, LEPRE1, NCAM1, NID2, NRP1, MSLN, RND1, S100A9, SCARB2, SERPINI1, SERPINF1, TGFB2, and VAV1) significantly reduced outgrowth of cocultured DRG neurons, indicating that endogenous expression of these genes in OECs supports neurite extension of DRG neurons. In a gain-of-function screen for 18 genes, six (CX3CL1, FZD1, LEPRE1, S100A9, SCARB2, and SERPINI1) enhanced and one (TIMP2) inhibited neurite growth. The most potent hit in both the loss- and gain-of-function screens was SCARB2, a protein that promotes cholesterol secretion. Transplants of fibroblasts that were genetically modified to overexpress SCARB2 significantly increased the number of regenerating DRG axons that grew toward the center of a spinal cord lesion in rats. We conclude that expression of SCARB2 enhances regenerative sprouting and that SCARB2 contributes to OEC-mediated neuronal repair.

In humans, the Crumbs homologue-1 (CRB1) gene is mutated in progressive types of autosomal recessive retinitis pigmentosa and Leber congenital amaurosis. The severity of the phenotype due to human CRB1 or mouse Crb1 mutations is dependent on the genetic background. Mice on C57BL/6J background with Crb1 mutations show late onset of retinal spotting phenotype or no phenotype. Recently, we showed that conditional deletion of mouse Crb2 in the retina results in early retinal disorganization leading to severe and progressive retinal degeneration with concomitant visual loss that mimics retinitis pigmentosa due to mutations in the CRB1 gene. Recent studies in the fruit fly and zebrafish suggest roles of the Crumbs (CRB) complex members in the regulation of cellular signalling pathways including the Notch1, mechanistic target of rapamycin complex 1 (mTORC1) and the Hippo pathway. Here, we demonstrate that mice backcrossed to C57BL/6J background with loss of CRB2 in the retina show a progressive disorganization and degeneration phenotype during late retinal development. We used microarray gene profiling to study the transcriptome of retinas lacking CRB2 during late retinal development. Unexpectedly, the retinas of newborn mice lacking CRB2 showed no changes in the transcriptome during retinal development. These findings suggest that loss of CRB2 in the developing retina results in retinal disorganization and subsequent degeneration without major changes in the transcriptome of the retina. These mice might be an interesting model to study the onset of retinal degeneration upon loss of CRB proteins.

PURPOSE: The ciliary body (CB) of the human eye consists of the non-pigmented (NPE) and pigmented (PE) neuro-epithelia. We investigated the gene expression of NPE and PE, to shed light on the molecular mechanisms underlying the most important functions of the CB. We also developed molecular signatures for the NPE and PE and studied possible new clues for glaucoma. METHODS: We isolated NPE and PE cells from seven healthy human donor eyes using laser dissection microscopy. Next, we performed RNA isolation, amplification, labeling and hybridization against 44×k Agilent microarrays. For microarray conformations, we used a literature study, RT-PCRs, and immunohistochemical stainings. We analyzed the gene expression data with R and with the knowledge database Ingenuity. RESULTS: The gene expression profiles and functional annotations of the NPE and PE were highly similar. We found that the most important functionalities of the NPE and PE were related to developmental processes, neural nature of the tissue, endocrine and metabolic signaling, and immunological functions. In total 1576 genes differed statistically significantly between NPE and PE. From these genes, at least 3 were cell-specific for the NPE and 143 for the PE. Finally, we observed high expression in the (N)PE of 35 genes previously implicated in molecular mechanisms related to glaucoma. CONCLUSION: Our gene expression analysis suggested that the NPE and PE of the CB were quite similar. Nonetheless, cell-type specific differences were found. The molecular machineries of the human NPE and PE are involved in a range of neuro-endocrinological, developmental and immunological functions, and perhaps glaucoma.

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.

Using the Braak staging for neurofibrillary changes as an objective indicator of the progression of Alzheimer's disease, we have performed a systematic search for global gene expression changes in the prefrontal cortex during the course of Alzheimer's disease. In the prefrontal cortex, senile plaques and neurofibrillary changes start to appear around Braak stage III, allowing for the detection of changes in gene expression before, during and after the onset of Alzheimer's disease neuropathology. Two distinct patterns of tightly co-regulated groups of genes were observed: (i) an increase in expression in early Braak stages, followed by a decline in expression in later stages (the UPDOWN clusters; containing 865 genes) and (ii) a decrease in expression in early Braak stages, followed by an increase in expression in later stages (the DOWNUP clusters; containing 983 genes). The most profound changes in gene expression were detected between Braak stages II and III, just before or at the onset of plaque pathology and neurofibrillary changes in the prefrontal cortex. We also observed an increase in intracellular beta amyloid staining from Braak stages I to III and a clear decrease in Braak stages IV to VI. These data suggest a link between specific gene expression clusters and Alzheimer's disease-associated neuropathology in the prefrontal cortex. Gene ontology over-representation and functional gene network analyses indicate an increase in synaptic activity and changes in plasticity during the very early pre-symptomatic stage of the disease. In later Braak stages, the decreased expression of these genes suggests a reduction in synaptic activity that coincides with the appearance of plaque pathology and neurofibrillary changes and the clinical diagnosis of mild cognitive impairment. The interaction of the ApoE genotype with the expression levels of the genes in the UPDOWN and DOWNUP clusters demonstrates that the accelerating role of ApoE-ε4 in the progression of Alzheimer's disease is reflected in the temporal changes in gene expression presented here. Since the UPDOWN cluster contains several genes involved in amyloid precursor protein processing and beta amyloid clearance that increase in expression in parallel with increased intracellular beta amyloid load, just before the onset of plaque pathology in the prefrontal cortex, we hypothesize that the temporally orchestrated increase in genes involved in synaptic activity represents a coping mechanism against increased soluble beta amyloid levels. As these gene expression changes occur before the appearance of Alzheimer's disease-associated neuropathology, they provide an excellent starting point for the identification of new targets for the development of therapeutic strategies aimed at the prevention of Alzheimer's disease.

Background To determine level, variability and functional annotation of gene expression of the human retinal pigment epithelium (RPE), the key tissue involved in retinal diseases like age-related macular degeneration and retinitis pigmentosa. Macular RPE cells from six selected healthy human donor eyes (aged 63–78 years) were laser dissected and used for 22k microarray studies (Agilent technologies). Data were analyzed with Rosetta Resolver, the web tool DAVID and Ingenuity software. Results In total, we identified 19,746 array entries with significant expression in the RPE. Gene expression was analyzed according to expression levels, interindividual variability and functionality. A group of highly (n = 2,194) expressed RPE genes showed an overrepresentation of genes of the oxidative phosphorylation, ATP synthesis and ribosome pathways. In the group of moderately expressed genes (n = 8,776) genes of the phosphatidylinositol signaling system and aminosugars metabolism were overrepresented. As expected, the top 10 percent (n = 2,194) of genes with the highest interindividual differences in expression showed functional overrepresentation of the complement cascade, essential in inflammation in age-related macular degeneration, and other signaling pathways. Surprisingly, this same category also includes the genes involved in Bruch's membrane (BM) composition. Among the top 10 percent of genes with low interindividual differences, there was an overrepresentation of genes involved in local glycosaminoglycan turnover. Conclusion Our study expands current knowledge of the RPE transcriptome by assigning new genes, and adding data about expression level and interindividual variation. Functional annotation suggests that the RPE has high levels of protein synthesis, strong energy demands, and is exposed to high levels of oxidative stress and a variable degree of inflammation. Our data sheds new light on the molecular composition of BM, adjacent to the RPE, and is useful for candidate retinal disease gene identification or gene dose-dependent therapeutic studies.

Olfactory ensheathing glia (OEG) are a specialized type of glia that support the growth of primary olfactory axons from the neuroepithelium in the nasal cavity to the brain. Transplantation of OEG in the injured spinal cord promotes sprouting of injured axons and results in reduced cavity formation, enhanced axonal and tissue sparing, remyelination, and angiogenesis. Gene expression analysis may help to identify the molecular mechanisms underlying the ability of OEG to recreate an environment that supports regeneration in the central nervous system. Here, we compared the transcriptome of cultured OEG (cOEG) with the transcriptomes of cultured Schwann cells (cSCs) and of OEG directly obtained from their natural environment (nOEG), the olfactory nerve layer of adult rats. Functional data mining by Gene Ontology (GO)‐analysis revealed a number of overrepresented GO‐classes associated with tissue repair. These classes include “response to wounding,” “blood vessel development,” “cell adhesion,” and GO‐classes related to the extracellular matrix and were overrepresented in the set of differentially expressed genes between both comparisons. The current screening approach combined with GO‐analysis has identified distinct molecular properties of OEG that may underlie their efficacy and interaction with host tissue after implantation in the injured spinal cord. These observations can form the basis for studies on the function of novel target molecules for therapeutic intervention after neurotrauma. © 2008 Wiley‐Liss, Inc.

PURPOSE:To describe gene expression differences between healthy, young human retinal pigment epithelium (RPE) cells from the macular area and RPE cells from two locations in the retinal periphery.METHODS:RPE cells from six human donor eyes, ages 17-36, without histopathological abnormalities, were dissected by laser and isolated from cryosections. Total RNA was isolated, amplified, and hybridized to a custom made oligonucleotide array containing 22,000 genes. Bioinformatic analysis was carried out using the computer programs Rosetta Resolver and the webtools EASE/David and GoStat. Confirmatory real time PCR (RT-PCR) and immunohistochemistry were performed according to standard protocols.RESULTS:Microarray and statistical analysis yielded 438 genes that were differentially expressed between macular RPE, and at least one out of two peripheral RPE locations. Out of these genes, 33 that showed fold changes of four, or higher, were selected for RT-PCR confirmation. For 17 genes (51%), a significant differential expression was found, while 11 additional genes (33%) showed a similar trend. Immuno-staining of one target (WFDC1) confirmed its differential expression on the protein level. Functional annotation and overrepresentation analysis independently defined extracellular matrix (ECM) genes as a statistically overrepresented class of genes in our RPE dataset. In total, 33 ECM genes were differentially expressed between macular and peripheral RPE regions. A subset of proteins corresponding to these genes is known to be present in Bruch's membrane.CONCLUSIONS:Our data showed that consistent topographical gene expression differences in the human RPE constitute around 1-5% of the RPE transcriptome. These changes may underlie topographical differences in RPE physiology, and pathology and may reflect local differences in the molecular composition and turnover of Bruch's membrane.