Next Article in Journal
Oxidants, Antioxidants and Thiol Redox Switches in the Control of Regulated Cell Death Pathways
Next Article in Special Issue
High Glucose Concentrations Affect Band 3 Protein in Human Erythrocytes
Previous Article in Journal
Polyphenols and Sesquiterpene Lactones from Artichoke Heads: Modulation of Starch Digestion, Gut Bioaccessibility, and Bioavailability following In Vitro Digestion and Large Intestine Fermentation
Previous Article in Special Issue
Management of Traumatic Brain Injury: From Present to Future
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of A2E-Induced Oxidative Stress on Retinal Epithelial Cells: New Insights on Differential Gene Response and Retinal Dystrophies

1
Department of Biomedical and Dental Sciences and Morphofunctional Imaging, Division of Medical Biotechnologies and Preventive Medicine, University of Messina, 98125 Messina, Italy
2
Department of Biomolecular Strategies, Genetics and Avant-Garde Therapies, I.E.ME.S.T., 90139 Palermo, Italy
3
Department of Chemical, Biological, Pharmaceutical and Environmental Sciences, University of Messina, 98125 Messina, Italy
*
Authors to whom correspondence should be addressed.
Antioxidants 2020, 9(4), 307; https://doi.org/10.3390/antiox9040307
Submission received: 20 March 2020 / Revised: 7 April 2020 / Accepted: 8 April 2020 / Published: 10 April 2020
(This article belongs to the Special Issue Cellular Oxidative Stress)

Abstract

:
Oxidative stress represents one of the principal inductors of lifestyle-related and genetic diseases. Among them, inherited retinal dystrophies, such as age-related macular degeneration and retinitis pigmentosa, are well known to be susceptible to oxidative stress. To better understand how high reactive oxygen species levels may be involved in retinal dystrophies onset and progression, we performed a whole RNA-Seq experiment. It consisted of a comparison of transcriptomes’ profiles among human retinal pigment epithelium cells exposed to the oxidant agent N-retinylidene-N-retinylethanolamine (A2E), considering two time points (3h and 6h) after the basal one. The treatment with A2E determined relevant differences in gene expression and splicing events, involving several new pathways probably related to retinal degeneration. We found 10 different clusters of pathways involving differentially expressed and differentially alternative spliced genes and highlighted the sub- pathways which could depict a more detailed scenario determined by the oxidative-stress-induced condition. In particular, regulation and/or alterations of angiogenesis, extracellular matrix integrity, isoprenoid-mediated reactions, physiological or pathological autophagy, cell-death induction and retinal cell rescue represented the most dysregulated pathways. Our results could represent an important step towards discovery of unclear molecular mechanisms linking oxidative stress and etiopathogenesis of retinal dystrophies.

1. Introduction

Oxidative stress, recently defined as “a state where oxidative forces exceed the antioxidant systems due to loss of the balance between them”, represents one of the principal inductors of lifestyle-related and genetic diseases [1]. Among them, retinal dystrophies and, in particular, age-related macular degeneration (AMD) and the subgroup of retinitis pigmentosa (RP) are well known to be susceptible to oxidative stress [2,3]. Nowadays AMD is considered the principal cause of visual disability among patients over 50 years [4]. The typical cellular sign of early AMD is represented by drusen, characteristic macular pigmentary deposits, associated with intermediate vision loss [5]. A “dry” and a “wet” form of AMD are currently known, the first more diffused but the second responsible for 90% of acute blindness due to AMD [6]. As previously mentioned, the most relevant risk factor associated with AMD etiopathogenesis is represented by high levels of oxidative stress damaging the macula, generally induced by production of advanced glycation end products (AGE) and exposition to environmental factors [7]. The majority of such effects are also exerted by dysregulation of vascular endothelial growth factor (VEGF), impairment proteins and organelles clearance and glial cell dysfunctions [8]. RP consists of a very heterogeneous group of inherited eye disorders characterized by progressive vision loss [9]. RP has an incidence of 1:4000 people worldwide and represents the most prevalent form of photoreceptor-related diseases [10]. Primary symptoms can already occur in childhood or adolescence and, generally, consist of night blindness and gradual reduction of the visual field due to progressive death of rods. Total blindness, instead, is a feature typical of the late stage of the disease, following degeneration of macula photoreceptors [11]. Rods represent about 95% of all photoreceptors, and oxidative metabolism of fatty acids is their main energy source [12]. Main causes of rod death are genetic mutations, and more than 80 RP-causative genes have been already identified (https://sph.uth.edu/RetNet/sum-dis.htm#B-diseases), even if a relevant number of them are still unknown [13]. Conversely, cone degeneration is usually a late event frequently resulting from cytotoxic effects of high oxygen levels in the retina after rod reduction. Thus, oxidative damage is considered the first cause of cone apoptosis and progressive vision loss [14]. Interestingly, AMD and RP can also arise due to mutations in genes expressed in other retinal cell types, such as MERTK [15], RLBP1 [16] and RPE65 [17] expressed in retinal pigment epithelium (RPE). Originally, only a trophic function was hypothesized for RPE cells. Nowadays, it is well known that RPE is a monolayer of neural-crista-derived pigmented epithelial cells interacting with Bruch’s membrane and choriocapillaris on the basolateral side and with the outer segments of the photoreceptors on the apical one [17]. RPE plays many vital roles for photoreceptor cells and the most fascinating certainly is the protection from oxidative stress [18]. Recent studies confirmed a high level of reactive oxygen species (ROS) in RPE, and fatty acids are one of their molecular targets. If oxidized, they may impair transduction pathways and gene expression [19]. Although fatty acid oxidation was already confirmed to cause macular degeneration, oxidative stress mechanism in RP development requires further clarifications [20]. Therefore, to better understand how high ROS levels may lead to retinal dystrophies onset and progression, we performed a comparison of transcriptomes’ profiles among human RPE cells exposed to the oxidant agent N-retinylidene-N-retinylethanolamine (A2E). A2E is a toxic bis-retinoid that derives from the condensation and the oxidation of the trans-retinal [21]. Throughout life, A2E and other complex lipids accumulate and form lipofuscin in the RPE, ultimately determining photoreceptor death [22]. Additionally, A2E photo-oxidation is able to generate singlet oxygen, a highly reactive molecule that contributes to the increase of level of toxic metabolites such as epoxides and endoperoxides [23]. Moreover, lipid peroxidation is also responsible for A2E cleavage that releases cytotoxic reactive aldehydes. These reactive species could affect the lipid membranes fluidity and can damage the DNA and the cellular proteins [24]. On these bases, it was hypothesized that lipofuscin, A2E and its oxidized metabolites could accumulate if the cellular antioxidant system is unable to fight the oxidative-damaged lipids in rods and cones. Consequently, accumulation and photo-oxidation of A2E could lead to several retinal dystrophies, like already established for age-related macular degeneration (AMD) [25]. Various studies based on ARPE-19 cell line present pharmacological solution to retinal dystrophies. The same cell line is also used to identify new candidate compounds able to protect RPE against A2E oxidation [26]. As highlighted in the manuscript, the treatment with such a compound determines relevant differences in gene expression and splicing events, involving several new pathways probably related to retinal degeneration.

2. Materials and Methods

2.1. Cell Culture

Human RPE-derived Cells (H-RPE—Human Retinal Pigment Epithelial Cells, Clonetics™, Lonza, Walkersville, USA) were grown in T-75 flasks containing RtEGM™ Retinal Pigment Epithelial Cell Growth Medium BulletKit® (Clonetics™, Lonza, Walkersville, MD, USA) with 2% v/v fetal bovine serum (FBS), 1% of penicillin/streptomycin and incubated at 37 °C with 5% CO2. H-RPE cells were then plated into 96-well plates (4 × 104 cells/well) and cultured for 24 h to reach confluence before treatment. Subsequently, A2E was added to a final concentration of 20 μM for 6 h before rinsing with medium. Control cell groups were incubated without A2E. Confluent cultures were transferred to PBS supplemented with calcium, magnesium and glucose (PBS–CMG) and then exposed to blue light emitted by a tungsten-halogen source (470 ± 20 nm; 0.4 mW/mm2) for 30 min to induce phototoxicity of A2E and incubated at 37 °C. The 1–3 generation of subcultured RPE cells were used in this study.

2.2. MTT Assay

Cell viability was determined by mitochondrial-dependent reduction of methylthiazolyldiphenyl-tetrazolium bromide (MTT) (Sigma-Aldrich, St. Louis, MO, USA) to formazan-insoluble crystals. Briefly, 10 μL of 5 mg/mL of MTT in PBS was added to the cells following the A2E treatment. After incubation at 37 °C for 2 h, 100 μL of 10% SDS in 0.01 mol/L HCl was added to dissolve the crystals and incubated for 16 h. The absorbance was measured in a Dynatech microplate reader at 570 nm. Results were expressed as percentage of viable cells normalized with control conditions in the absence of A2E.

2.3. Total RNA Sequencing

Total RNA was extracted by TRIzolTM Reagent (InvitrogenTM, ThermoFisher Scientific, Waltham, MA, USA), following manufacturer’s protocol, and quantified at Qubit 2.0 fluorimeter by Qubit® RNA assay kit (InvitrogenTM, ThermoFisher Scientific, Waltham, MA, USA). The RNA-seq samples consisted of 3 factor groups, represented by Human RPE cells, before the treatment with A2E and at 2 following different time points of 3 h and 6 h, respectively. For each group, 3 biological replicates were considered, for a total of 9 samples. The selection of 3 h and 6 h time points was based on previous experiments realized by our research group (unpublished data), confirmed by outcomes from MTT assay in this work. Such results highlighted that in wider time intervals, the death rate of oxidative stressed cells might be so high as to invalidate the following expression analysis. Libraries were generated using 1 µg of total RNA by the TruSeq Stranded Total RNA Sample Prep Kit with Ribo-Zero H/M/R (Illumina, San Diego, CA, USA), according to manufacturer’s protocols. Sequencing runs were performed on an HiSeq 2500 Sequencer (Illumina, San Diego, CA, USA), using the HiSeq SBS Kit v4 (Illumina, San Diego, CA, USA). The experiment was repeated thrice.

2.4. Quality Validation and Read Mapping

Sequence reads were generated from RPE-specific cDNA libraries on the Illumina HiSeq 2500 Sequencer. Obtained raw sequences were filtered to remove low-quality reads (average per base Phred score <30) and adaptor sequences. The quality of analyzed data was checked using FastQC (v.0.11.9) (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and QualiMap (v.2.2.1) [27], while trimming was realized by Trimmomatic (v.0.39). Filtered data were then mapped by CLC Genomics Workbench v.20.0 (https://digitalinsights.qiagen.com/products-overview/analysis-and-visualization/qiagen-clc-genomics-workbench/) against the Homo sapiens genome hg38 and the RNA database v.91, on Ensembl database. RNA-seq analysis was conducted using the following settings: quality trim limit = 0.01, ambiguity trim maximum value = 2. Map to annotated reference was as follows: mismatch cost = 2, insertion and deletion costs = 3, minimum length fraction and minimum similarity fraction = 0.8, maximum number of hits for a read = 10, strand-specific = both, expression value = TPM. Raw data are available upon request.

2.5. Gene Expression and Statistical Analysis

Mapped reads were quantified by alignment-dependent and alignment-independent methods. The first approach uses the expectation-maximization (EM) algorithm [28] in order to determine expressions even in cases where the majority of reads map equally well to multiple genes or transcripts. Once the algorithm has converged, every non-uniquely mapping read was assigned randomly to a particular transcript according to the abundances of the transcripts within the same mapping. The transcript per million reads (TPM) values were, then, computed from the counts assigned to each transcript. The second method has a higher accuracy for the point expression estimation and also allows the user to bootstrap the expression quantification to get an estimate of the technical variability. This approach was applied by the Salmon tool [29] using the transcript fasta files downloaded from GENCODE v.32 (gencode.v32.transcripts.fa). Salmon was run with the following settings for the RNA-seq data: quant –index/index --libType U --unmatedReads /single.fastq --incompatPrior ‘0.0’ --biasSpeedSamp ’5 ’ --fldMax ‘1000’ --fldMean ‘250’ --fldSD ‘25’ --forgettingFactor ‘0.65’ --maxReadOcc ‘100’ --numBiasSamples ‘2000000’ --numAuxModelSamples ‘5000000’ --numPreAuxModelSamples ‘5000’ --numGibbsSamples ‘0’ --numBootstraps ‘0’ --thinningFactor ‘16’ --sigDigits ‘3’ --vbPrior ‘1e-05’ -o/output. Once obtained, Salmon outputs were imported using tximport R package version 1.10.0 and lengthScaledTPM method [30] to generate read counts and Transcripts Per Million (TPMs). Low expressed transcripts and genes were filtered based on the data mean–variance trend analysis. The expected decreasing trend between data mean and variance was observed when expressed transcripts were determined as which had 1 of the 9 samples with count per million reads (CPM) 1, which provided an optimal filter of low expression. A gene was expressed if any of its transcripts with the above criteria were expressed. The trimmed mean of M-values (TMM) method was used to normalize the gene and transcript read counts to l o g 2 -CPM [31]. The principal component analysis (PCA) plot showed that the RNA-seq data did not have distinct batch effects, so it was possible to proceed with downstream analyses.

2.6. DE, DAS and DTU Analysis

Limma R package was used for differential expression analyses [32]. General linear models were established to compare gene and transcript expression changes at the different conditions of experimental design, setting the contrast groups as 0 h.untreated versus 3 h.treated, 0 h.untreated versus 6 h.treated, 3 h.treated versus 6 h.treated, 0 h.untreated versus (3 h.treated + 6 h.treated)/2. For differentially expressed (DE) genes/transcripts, the l o g 2 fold change ( L 2 F C ) of gene/transcript abundance were calculated based on contrast groups and significance of expression changes were determined using the t-test [33]. P-values of multiple testing were adjusted with BH to correct false discovery rate (FDR) [34]. A gene/transcript was significantly DE in a contrast group if it had adjusted p-value <0.01 and L 2 F C 1. At the alternative splicing level, differential transcript usage (DTU) transcripts were determined by comparing the L 2 F C of a transcript to the weighted average of L 2 F C s (weights were based on their standard deviation) of all remaining transcripts in the same gene. A transcript was determined as significant DTU if it had adjusted p-value <0.01 and Δ PS 0.1. For differentially alternative spliced (DAS) genes, L 2 F C of each individual transcript was compared to gene level L 2 F C , which was calculated as the weighted average of L 2 F C s of all transcripts of the gene. Then p-values of individual transcript comparisons were summarized to a single gene level p-value with F-test. A gene was significantly DAS in a contrast group if it had an adjusted p-value <0.01 and any of its transcript had a Δ Percent Spliced ( Δ PS) ratio 0.1. Finally, time points (0 h, 3 h, 6 h) in groups (untreated, treated) were used for time-series trend analysis. Natural Cubic Spline method with degree of freedom was used to generate time-series trend (Figure S1).

2.7. Gene-Enrichment and Functional Pathway Analysis

The up- and down-regulated genes were analyzed by the Database for Annotation, Visualization and Integrated Discovery (DAVID) 6.8 [35]. This tool is based on more than 40 annotation categories, including GO terms to protein–protein interactions, from disease associations to gene functional summaries, and many others. In DAVID annotation system, EASE Score, a modified Fisher Exact P-Value, is adopted to measure the gene-enrichment in annotation terms. The EASE score provides a conservative adjustment to the Fisher exact probability that weights significance in favor of themes supported by more genes. The EASE score is calculated by penalizing (removing) one gene within the given category from the list and calculating the resulting Fisher exact probability for that category.

2.8. Selection of Single-Pathway “Master genes” and Selection of Retinitis Pigmentosa Candidate Genes by ToppGene Prioritization

In order to highlight new candidate genes involved into retinitis pigmentosa, based on oxidative-related candidate pathways, we chose the 15 most altered genes for each one. Firstly, we considered them for each time point; then, we chose the commons in all time points. Subsequently, chosen genes underwent prioritization by ToppGene (https://toppgene.cchmc.org), a web-based tool able to classify a selected group of candidate genes from a large set of genes correlated with a pathology, giving each one a score. The score is based on the intersection of data from various databases of annotations related to cellular and physiological functions, analyzing complex networks shared between genes already known to cause the disease (training genes) and candidate genes (test genes). Training genes were obtained from RetNet online database.

2.9. Data Validation by qRT-PCR.

Ten most dysregulated mRNAs from candidate genes previously identified were selected and validated by quantitative Real-Time-Polymerase Chain Reaction (qRT-PCR), in order to validate RNA-Seq data. Reverse transcription was carried out according to the manufacturer’s protocol of GoScript™ Reverse Transcription System (Promega, Madison, WI, USA). The obtained cDNA was subjected to RT-PCR in the ABI 7500 fast sequence detection system (Applied Biosystems, Foster City, CA, USA), using BRYT-Green-based PCR reaction. PCR amplification was performed in a total reaction mixture of 20 μL containing 20 ng cDNA, 10 μL 2 × GoTaq1qPCR Master Mix (Promega, Madison, WI, USA) and 0.2 μM of each primer. PCR was run with the standard thermal cycle conditions using the two-step qRT-PCR method: an initial denaturation at 95 °C for 30 s, followed by 40 cycles of 30 s at 95 °C and 30 s at 60 °C. Each reaction was replicated six times, considering all analyzed time points (18 samples), and the average threshold cycle (Ct) was calculated for each replicate. The expression of mRNAs was calculated relative to expression level of endogenous control β-actin, and the relative expression of genes was calculated using the 2−ΔΔCt method [36]. The results were shown as the mean ± SEM (Standard Error of Mean). Statistical significance was determined by analysis of variance between groups (ANOVA), followed by Bonferroni post-hoc test. Finally, a linear regression analysis was performed to check the correlation of the FC of the gene expression ratios between qRT-PCR and RNA-Seq. The statistical analyses were all performed using IBM SPSS 26.0 software (https://www.ibm.com/analytics/us/en/technology/spss/). Adjusted p-values <0.05 were considered statistically significant. The research was approved by the Scientific Ethics Committee of the Azienda Ospedaliera Universitaria—Policlinico “G. Martino” Messina.

3. Results

3.1. MTT Cell Viability Assay Results

The MTT cell viability assay showed a relevant and different trend in RPE-treated cells versus control. The addition of A2E to cultures led to a dose-dependent increase in cell death percentage (Figure 1).

3.2. Sequencing Analysis and Mapping Statistics

RNA sequencing carried out on Illumina HiSeq 2500 yielded a total of 96,346,180 quality reads (mean mapping quality = 29) with a percentage of 67.8% uniquely mapped. A total of 16,173 genes and 69,653 transcripts were identified out of 60,609 and 227,462 reference counterparts, respectively, considering the whole human transcriptome. The annotated reference assembly v.32 (GRCh38.p13) was downloaded from GeneCode FTP server (ftp://ftp.ebi.ac.uk/pub/databases/gencode/Gencode_human/). All previous mapping statistics were based on average values calculated for all three replicates in each time point. Details are available in Figure S2.

3.3. Analysis of Gene Expression Profile of RPE Cells

As previously cited, in our transcriptome study, 16,173 genes (Table S1) and 69,653 transcripts (Table S2) were expressed with log2-CPM values ≥1 considered as the average value per sample. We first analyzed differential expression at the gene level (DE) (Figure 2a and Table S3). Considering the stringent criteria we have chosen, we identified 2432 genes that were significantly differentially expressed in response to A2E treatment. Of these, 59.7% resulted as up-regulated, while 40.3% down-regulated, showing a decreasing time-dependent trend for up-regulated ones and an increasing trend for down-regulated ones. Afterward, the analysis of transcript-level data allowed us to identify genes that were DAS between the contrast groups (Figure 2b and Table S4). We detected 5119 DAS genes, of which 1101 were also DE genes (regulated by both transcription and AS) and 4108 were regulated by AS only. Therefore, considering all 6540 genes that showed significantly altered levels of differential gene and/or transcript-level expression, 78.3% were differentially alternatively spliced, highlighting a consistent response to the induced oxidative stress. Furthermore, to identify the specific transcripts that characterize a gene as DAS, a DTU analysis was performed (Figure 2b and Table S5). Globally, around 12% (8587) of expressed transcripts were classified as DTU. The next step consisted in the analysis of early response to the induced oxidative stress in comparison with the response to late treatment with A2E. About 40% (858) and 55% (2457) of DE and DAS genes, respectively, resulted as common to both time point observations, while the residual was unique to 3 h and 6 h (Figure 3). Consequently, it is evident how changes in gene-level expression and alternative splicing occurred throughout the whole period, either transiently (occurring after 3 h and returning to initial level after 6 h), occurring later (only after 6 h from treatment) or enduring throughout the whole period. Such results probably suggest the different responses of RPE cells as a growing resistance to oxidative stress stimuli in a time-dependent manner.

3.4. DE and DAS Genes Highlighted Different Functionality Patterns

Functional enrichment analyses of DE and DAS genes showed relevant differences, as already supposed by the overlap of only 1011 genes (Figure 2). The most significantly enriched terms for DE genes were linked to vascular events and very heterogeneous responses to oxidative stress. As evidenced by hierarchical clustering of total gene expression levels of DE genes, adaptive, transient and late expression profiles followed the A2E induced stress, and an analog response was highlighted by transcript expression profiles of individual DE genes (Figure 4). Functional annotation of individual DE genes clusters was associated with methylation of RNA (cluster 1), various stress response (cluster 2), microtubule assembly and activity (cluster 3), angiogenesis (cluster 6), lipid biosynthesis (cluster 7), focal adhesions (cluster 8) and respiratory electron transport (cluster 10). For DAS genes, the most enriched functional terms deal with alteration of cellular proliferation and cell death. Hierarchical clustering of the DTU transcripts and expression profiles of individual DAS genes also evidenced heterogeneous response profiles. Functional annotation of DTU clusters of genes revealed enrichment of terms related to misfolded and damaged protein removal (clusters 1 and 3), autophagy (cluster 2), lipid biosynthesis (cluster 6), regulation of DNA damage response (cluster 7), induction of cell death (cluster 8) and translation regulation (cluster 10) (Figure 5).

3.5. Early Cellular Response to Induced Stress Mainly Involves Pre-mRNA Splicing and Glycolysis-Related DE and DAS Genes

The applied statistical model used in our analyses permitted us to determine precisely at which time point each DE and DAS gene first showed a significant change, together with the magnitude and course of that change. The top 15 significant genes or transcripts belonging to the considered groups were deeply investigated, allowing the identification of new candidate genes and specific pathways possibly involved in retinal dystrophy etiopathogenesis (Figure 6). After the first 3 h of treatment, only one gene (HNRNPA3P6) reached a relevant DE status, and it is involved in cytoplasmic trafficking of RNA and pre-mRNA splicing. Main DAS genes that showed early important differences (CTSH and GPI), instead, were related to glycolysis, and the same biochemical pathway was the most correlated to initial DTU transcripts expression dysregulation (ENST00000615999.4, ENST00000588991.7, ENST00000586425.2, ENST00000525807.5 and ENST00000550050.5).

3.6. Late RPE Cell Response to A2E Treatment Could Impair Bioenergetic Specific Reactions, Extracellular Matrix Integrity and Neurotransmission-Related DE and DAS Genes

While the number of significant and relevant DE genes was limited in the early stage of treatment, it grew over time, reaching a climax at 6 h (Figure 6). At this time point, many of the DE genes (ACTG1, CCN2, RPL19, RPL3, P4HB, RPS11, FTL, CAPZB and RNA5-8SN2) (Figure 7 and Figure 8) resulted as linked to new particular pathways, as iron metabolism, plasma lipoproteins assembly and F-actin capping. Furthermore, several over-expressed DE genes at 6 h (TTC8, ARL3, REEP6, GUCA1B and PDE6G) resulted as associated with alternative splicing of retina-preferred gene transcripts. About DAS genes, instead, the most significant ones (ACADVL, GPI and LTBP3) resulted as correlated to already-involved pathways (e.g., glycolysis). However, more interestingly, the same genes were together with DTU transcripts (ENST00000586425.2 and ENST00000588991.7), but even more interestingly, with other biochemical activities regarding oxidative processes in mitochondria (e.g., fatty acids reactions) and low conductance of potassium channels (as evidenced by DE transcripts ENST00000577650.5, ENST00000451956.1, ENST00000556690.5 and ENST00000551173.5).

3.7. The Transcriptome Comparison between Untreated (Time Zero) and Treated (3 h + 6 h) Rpe Cells Revealed the Possible Impairment of Retinal Cells Crosstalk and Synapses, Leading to Rescue or Cell Death

Many DE and DAS genes only showed significant differences only if analyzed during the whole considered treatment period. Among these DE genes, five (CCN2, ACTG1, UTP14C, TMSB4XP6 and TMEM189-UBE2V1) were linked to extracellular matrix constituent secretion and cellular junctions, as well as to misfolded protein ubiquitination. Furthermore, the number of DAS genes and DTU transcripts changed during treatment was very high. The first ones (ACADVL, GPI, HNRNPA1, CD81, CD63, CTSB, CTSH, LTBP3, CAPNS1 and PAX6) were enriched in terms related to neuropeptide catabolic processes and extracellular vesicles in the crosstalk of cells, while the second ones (ENST00000586425.2, ENST00000550050.5, ENST00000588991.7, ENST00000527343.5, ENST00000584364.5, ENST00000263645.9, ENST00000615999.4 and ENST00000518154.5) mainly regard dendrite regeneration. Therefore, the overall scenario could reveal a dual response by stressed cells: on one hand, the alteration of retinal cell crosstalk and synapses could lead to various forms of cell death (e.g., autophagy). On the other hand, RPE tries to survive, increasing regeneration of capable parts (e.g., dendrites), by synaptic plasticity. Detailed info on pathways enrichment are available in Table S6 and Table S7).

3.8. The Most Significant DAS Genes Represented the Main Retinal Dystrophy Candidate Genes

The ToppGene prioritization analysis on known causative retinal dystrophy genes that intersected with most significant DE and DAS genes previously described revealed 19 candidate genes (Bonferroni corrected p-value <0.05). Of these, seven showed a strong association with the training genes (Bonferroni corrected p-value <0.01). Most of the 19 significant obtained candidate genes, included three of most statistically associated (PAX6, CTSH and HNRPA1) belong to DAS genes, further highlighting the influence of oxidative stress on alternative splicing (Table 1). Details of ToppGene results are available in Table S8.

3.9. qRT-PCR Validation

To validate the authenticity and reproducibility of the RNA-Seq results, the 15 selected mRNAs were validated by qRT-PCR analysis, and obtained expression profiles were similar to the results of transcriptome analysis (Table S9). Moreover, the ANOVA statistics, conducted to compare the means among multiple groups, highlighted high significance (p-values < 0.05). The linear regression analysis showed a significantly positive correlation of the relationship between gene expression ratios of qRT-PCR and RNA-Seq for all selected time points (Figure S3), confirming our transcriptomic data validity.

4. Discussion

Retinal dystrophies like age-related macular degeneration and, particularly, retinitis pigmentosa represent a very heterogeneous group of ocular pathologies characterized by a very complex pattern of environmental and genetic causes. One of the most challenging aspects regards the incomplete knowledge of all causative genes and their involved biochemical and molecular pathways, leading to a huge group of orphan forms [37]. Gene mutations or dysfunctional processes not only in the retina but also in RPE could cause inherited retinal degeneration, age-related macular degeneration and other retinal diseases [38]. Such a feature highlights the relevant role of RPE, a high metabolic demand monolayer of pigmented cells that plays fundamental functions for both rods and cones, such as metabolite transport and photoreceptor excitability, regulation of visual cycle, secretion of growth factors, phagocytosis of photoreceptor outer segments (POSs) and oxidative stress defense. Regarding the latter point, oxidative stress represents one of the major lethal mechanisms responsible for age-related RPE damages [39]. Many studies have demonstrated that accumulation of lipid deposit called lipofuscin generates reactive oxygen species through phototoxicity in RPE cells [40,41,42]. Oxidative stress triggered by photo-oxidation of bis-retinoid A2E, a lipofuscin constituent, is well known to be a progression factor of age-related macular degeneration and also in genetic macular degeneration syndromes such as Stargardt disease [43], but very little is known about A2E involvement in retinitis pigmentosa.
In this study, we treated RPE cells with A2E during a follow-up of two time points (3 h and 6 h) after exposure and compared them to untreated time zero controls. The main purpose of our experiment was the discovery of new pathways potentially involved in retinal dystrophies development, with the further detection of new candidate genes that could be associated or causative of such ocular diseases, emerging from the expression analysis in such altered conditions.
Starting from an initial average value of 16,173 detected genes per sample, about 2432 showed changes in their expression level and 119 were differentially alternative spliced with transiently, late or enduring fluctuations. Selected altered DE and DAS genes were then functionally and statistically analyzed and clustered into final 10 candidate “macro-pathways”, showing a very variable time of exposure-related trends. Considering the average fold-change of each constituting genes and their reciprocal connections, we revealed a more detailed functional network. Such connections could help to depict several causative/associative clusters, underlining a more complex pattern of possible retinal dystrophies etiologies.
From analyses of DE genes related pathways, it emerged that, after an early increase of apoptosis processes, the programmed cell death decreased in both considered time points following A2E exposure, probably due to activation of rescue systems and to a limited percentage of survived cells. An opposite trend was shown by steroid receptor and nucleoside transport activities, which evidenced a huge up-regulation of involved genes after 3 h and 6 h. Such results could reflect late alterations in RPE antioxidant and anti-inflammatory abilities [44], as well as inhibition of photoreceptor outer segment (POS) phagocytosis and impairment of ion currents in retinal cells [45,46]. Furthermore, a very interesting result was the sinusoidal trend involving isoprenoid pathway, related to cholesterol-dependent homeostasis of POS [47] and angiogenesis [48]. After 3 h from treatment, it could be possible that disc bulk membranes increased trying to improve phototransduction by residual photoreceptors, despite the decreased choriocapillaris viability due to the reduced vascular endothelial growth factor (VEGF) [49]. The situation reversed at 6 h observation, when the discs turnover was drastically decreased and the VEGF level increased, which could contribute to the subretinal neovascularization already characterized in wet age macular degeneration [50].
Additionally, the extensive alternative splicing information identified from DAS-related pathway analysis highlighted a much higher degree of complexity of regulation in response to A2E-induced oxidative stress, which has been significantly underestimated by analysis of DE genes only. In particular, speed and extent of the oxidative stress-induced AS suggested that AS, together with the transcriptional response, is a major driver of transcriptome reprogramming for RPE cell death and their attempts to survive. From 3 h and up to 6 h from treatment, an impairment of intracellular traffic related to Rab proteins, already reported in choroideremia [51], was observed along with the alteration of autophagy and accumulation of proteins and damaged organelles. These events are typical of AMD [52] and are also enforced by inactivation of chaperone genes [53,54]. This scenario could reflect a strong reduction of macroautophagy (a catabolic cell survival system) and of a hybrid autophagy–phagocytosis-degradative pathway called LC3-associated phagocytosis (LAP) [55], which plays a critical role in visual pigment regeneration, as well as the complete degradation of phagosomes [56]. The other resulted pathway “big cluster”, instead, highlighted a global up-regulation of DAS-involved genes up to the end of the last considered time point. In particular, we demonstrated the dynamic contribution of AS by changes in multiple different mechanisms of transcription and translation. Photosensitization of A2E stimulates oxidative DNA damage, such as the production of 8-oxo-guanines in telomeres, leading to their possible damage [57]. Thus, the resulted alteration of telomerase RNA localization to Cajal body could accelerate the RPE senescence [58], and this process could be further increased by the reduction of FGFR1 signaling and the consequent POS phagocytosis decrease [59]. Moreover, DNA damage response could increase the activity of miRNAs involved in [60,61] and cell death genes transcription by TP53 [62], determining a possible role in retina degeneration. In the meantime, as already discussed, retinal cells could try to fight against induced stress, and one resistance mechanism could be represented by the improved maturation of the large ribosomal unit (LSU) [63] and by the increased polyribosome activity, showing a fundamental role in translation regulation of many retinal genes [64,65]. Detailed analysis of DE and DAS gene-involved pathways, with their possible impact on retinal dystrophies etiopathogenesis, is reported in Table 2.

5. Conclusions

We realized a whole RNA-seq experiment on RPE cells treated with A2E, considering two time points (3 h and 6 h) after the basal one. We found 10 different clusters of pathways involving DE and DAS genes, with many highlighted sub-pathways, which could depict a more detailed scenario determined by induced oxidative stress. Regulation and/or alterations of angiogenesis, extracellular matrix integrity, isoprenoid-mediated reactions, physiological or pathological autophagy, cell death induction and retinal cell rescue represented the most dysregulated pathway, probably involved in retinal degeneration. Assembly of splicing and transcriptional networks from analyzed data will further define the contribution of AS, as an extra level of regulation, and the interplay and coordination of the transcriptional and AS responses. However, it is fundamental to highlight several limitations of our study: RPE-cultured cells were not in contact with photoreceptors’ outer segments and did not perform any phagocytosis. Moreover, the short-term response (at 3 h and 6 h) detected in vitro do not surely reflect what happens in vivo. Finally, even if it was important to underline the importance to realize an in vivo experiment to confirm what observed in RPE cells, our results could represent an important step towards discovery of unclear molecular mechanisms involved in etiopathogenesis of retinal dystrophies.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-3921/9/4/307/s1, Figure S1: Time-series trend analysis, Figure S2: Gene expression data statistics, Figure S3: qRT-PCR validation of ten most differentially expressed genes, Table S1: Gene expressed in RPE cells transcriptome analysis, Table S2: Gene transcripts expressed in RPE cells transcriptome analysis, Table S3: DE genes profile of RPE Cells, Table S4: DAS genes profile of RPE Cells, Table S5: DTU transcripts profile of RPE Cells, Table S6: Significant DE vs DAS gene enrichment by DAVID, Table S7: Significant DE vs DAS gene enrichment by ClueGO, Table S8: ToppGene prioritization detailed results, Table S9: Expression profiles of 10 among most DE genes.

Author Contributions

Conceptualization, L.D.; methodology, L.D.; software, L.D. and C.S.; validation, C.S. and S.A.; formal analysis, L.D.; investigation, L.D.; resources, C.S.; data curation, L.D. and C.S.; writing—original draft preparation, L.D.; writing—review and editing, C.S., R.D. and C.R.; visualization, C.R. and R.D.; supervision, A.S.; project administration, A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sies, H. Oxidative stress: A concept in redox biology and medicine. Redox Biol. 2015, 4, 180–183. [Google Scholar] [CrossRef] [Green Version]
  2. Donato, L.; Scimone, C.; Nicocia, G.; D’Angelo, R.; Sidoti, A. Role of oxidative stress in Retinitis pigmentosa: New involved pathways by an RNA-Seq analysis. Cell Cycle 2019, 18, 84–104. [Google Scholar] [CrossRef] [Green Version]
  3. Datta, S.; Cano, M.; Ebrahimi, K.; Wang, L.; Handa, J.T. The impact of oxidative stress and inflammation on RPE degeneration in non-neovascular AMD. Prog. Retin. Eye Res. 2017, 60, 201–218. [Google Scholar] [CrossRef] [PubMed]
  4. Randolph, S.A. Age-related macular degeneration. Workplace Health Saf. 2014, 62, 352. [Google Scholar] [CrossRef] [PubMed]
  5. Gheorghe, A.; Mahdi, L.; Musat, O. Age-Related Macular Degeneration. Rom. J. Ophthalmol. 2015, 59, 74–77. [Google Scholar] [PubMed]
  6. Lim, L.S.; Mitchell, P.; Seddon, J.M.; Holz, F.G.; Wong, T.Y. Age-related macular degeneration. Lancet 2012, 379, 1728–1738. [Google Scholar] [CrossRef]
  7. Kivinen, N.; Koskela, A.; Kauppinen, A.; Kaarniranta, K. Pathogenesis of age-related macular degeneration—dialogue between autophagy and inflammasomes. Duodecim 2017, 133, 641–646. [Google Scholar] [PubMed]
  8. Xu, B.; Zhang, H.; Zhu, M.; Le, Y.Z. Critical Role of Trophic Factors in Protecting Muller Glia: Implications to Neuroprotection in Age-Related Macular Degeneration, Diabetic Retinopathy, and Anti-VEGF Therapies. Adv. Exp. Med. Biol. 2019, 1185, 469–473. [Google Scholar] [CrossRef]
  9. Natarajan, S. Decoding retinitis pigmentosa. Indian J. Ophthalmol. 2013, 61, 91–94. [Google Scholar] [CrossRef]
  10. Verbakel, S.K.; van Huet, R.A.C.; Boon, C.J.F.; den Hollander, A.I.; Collin, R.W.J.; Klaver, C.C.W.; Hoyng, C.B.; Roepman, R.; Klevering, B.J. Non-syndromic retinitis pigmentosa. Prog. Retin. Eye Res. 2018, 66, 157–186. [Google Scholar] [CrossRef]
  11. Fahim, A. Retinitis pigmentosa: Recent advances and future directions in diagnosis and management. Curr. Opin. Pediatr. 2018, 30, 725–733. [Google Scholar] [CrossRef] [PubMed]
  12. Agbaga, M.P.; Merriman, D.K.; Brush, R.S.; Lydic, T.A.; Conley, S.M.; Naash, M.I.; Jackson, S.; Woods, A.S.; Reid, G.E.; Busik, J.V.; et al. Differential composition of DHA and very-long-chain PUFAs in rod and cone photoreceptors. J. Lipid Res. 2018, 59, 1586–1596. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Birtel, J.; Gliem, M.; Oishi, A.; Muller, P.L.; Herrmann, P.; Holz, F.G.; Mangold, E.; Knapp, M.; Bolz, H.J.; Charbel Issa, P. Genetic testing in patients with retinitis pigmentosa: Features of unsolved cases. Clin. Exp. Ophthalmol. 2019, 47, 779–786. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Campochiaro, P.A.; Mir, T.A. The mechanism of cone cell death in Retinitis Pigmentosa. Prog. Retin. Eye Res. 2018, 62, 24–37. [Google Scholar] [CrossRef]
  15. Audo, I.; Mohand-Said, S.; Boulanger-Scemama, E.; Zanlonghi, X.; Condroyer, C.; Demontant, V.; Boyard, F.; Antonio, A.; Mejecase, C.; El Shamieh, S.; et al. MERTK mutation update in inherited retinal diseases. Hum. Mutat. 2018, 39, 887–913. [Google Scholar] [CrossRef]
  16. Scimone, C.; Donato, L.; Esposito, T.; Rinaldi, C.; D’Angelo, R.; Sidoti, A. A novel RLBP1 gene geographical area-related mutation present in a young patient with retinitis punctata albescens. Hum. Genom. 2017, 11, 18. [Google Scholar] [CrossRef]
  17. Miraldi Utz, V.; Coussa, R.G.; Antaki, F.; Traboulsi, E.I. Gene therapy for RPE65-related retinal disease. Ophthalmic Genet. 2018, 39, 671–677. [Google Scholar] [CrossRef]
  18. Sparrow, J.R.; Hicks, D.; Hamel, C.P. The retinal pigment epithelium in health and disease. Curr. Mol. Med. 2010, 10, 802–823. [Google Scholar] [CrossRef]
  19. Nowak, J.Z. Oxidative stress, polyunsaturated fatty acids-derived oxidation products and bisretinoids as potential inducers of CNS diseases: Focus on age-related macular degeneration. Pharmacol. Rep. 2013, 65, 288–304. [Google Scholar] [CrossRef]
  20. Donato, L.; Scimone, C.; Rinaldi, C.; D’Angelo, R.; Sidoti, A. Non-coding RNAome of RPE cells under oxidative stress suggests unknown regulative aspects of Retinitis pigmentosa etiopathogenesis. Sci. Rep. 2018, 8, 16638. [Google Scholar] [CrossRef]
  21. Crouch, R.K.; Koutalos, Y.; Kono, M.; Schey, K.; Ablonczy, Z. A2E and Lipofuscin. Prog. Mol. Biol. Transl. Sci. 2015, 134, 449–463. [Google Scholar] [CrossRef] [PubMed]
  22. Sparrow, J.R.; Fishkin, N.; Zhou, J.; Cai, B.; Jang, Y.P.; Krane, S.; Itagaki, Y.; Nakanishi, K. A2E, a byproduct of the visual cycle. Vis. Res. 2003, 43, 2983–2990. [Google Scholar] [CrossRef] [Green Version]
  23. Parmar, V.M.; Parmar, T.; Arai, E.; Perusek, L.; Maeda, A. A2E-associated cell death and inflammation in retinal pigmented epithelial cells from human induced pluripotent stem cells. Stem Cell Res. 2018, 27, 95–104. [Google Scholar] [CrossRef] [PubMed]
  24. Dontsov, A.; Koromyslova, A.; Ostrovsky, M.; Sakina, N. Lipofuscins prepared by modification of photoreceptor cells via glycation or lipid peroxidation show the similar phototoxicity. World J. Exp. Med. 2016, 6, 63–71. [Google Scholar] [CrossRef]
  25. Alaimo, A.; Linares, G.G.; Bujjamer, J.M.; Gorojod, R.M.; Alcon, S.P.; Martinez, J.H.; Baldessari, A.; Grecco, H.E.; Kotler, M.L. Toxicity of blue led light and A2E is associated to mitochondrial dynamics impairment in ARPE-19 cells: Implications for age-related macular degeneration. Arch. Toxicol. 2019, 93, 1401–1415. [Google Scholar] [CrossRef] [PubMed]
  26. Alaimo, A.; Di Santo, M.C.; Dominguez Rubio, A.P.; Chaufan, G.; Garcia Linares, G.; Perez, O.E. Toxic effects of A2E in human ARPE-19 cells were prevented by resveratrol: A potential nutritional bioactive for age-related macular degeneration treatment. Arch. Toxicol. 2019, 94, 553–572. [Google Scholar] [CrossRef]
  27. Okonechnikov, K.; Conesa, A.; Garcia-Alcalde, F. Qualimap 2: Advanced multi-sample quality control for high-throughput sequencing data. Bioinformatics 2016, 32, 292–294. [Google Scholar] [CrossRef]
  28. Li, B.; Ruotti, V.; Stewart, R.M.; Thomson, J.A.; Dewey, C.N. RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics 2010, 26, 493–500. [Google Scholar] [CrossRef] [Green Version]
  29. Patro, R.; Duggal, G.; Love, M.I.; Irizarry, R.A.; Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 2017, 14, 417–419. [Google Scholar] [CrossRef] [Green Version]
  30. Soneson, C.; Matthes, K.L.; Nowicka, M.; Law, C.W.; Robinson, M.D. Isoform prefiltering improves performance of count-based methods for analysis of differential transcript usage. Genome Biol. 2016, 17, 12. [Google Scholar] [CrossRef] [Green Version]
  31. Pereira, M.B.; Wallroth, M.; Jonsson, V.; Kristiansson, E. Comparison of normalization methods for the analysis of metagenomic gene abundance data. BMC Genom. 2018, 19, 274. [Google Scholar] [CrossRef] [Green Version]
  32. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
  33. Guo, W.T.N.; Stephen, G.; Milne, I.; Calixto, C.; Waugh, R.; Brown, J.W.; Zhang, R. 3D RNA-seq—A powerful and flexible tool for rapid and accurate differential expression and alternative splicing analysis of RNA-seq data for biologists. bioRxiv 2019, 656686. [Google Scholar] [CrossRef]
  34. Ge, Y.; Sealfon, S.C.; Speed, T.P. Some Step-down Procedures Controlling the False Discovery Rate under Dependence. Stat. Sin. 2008, 18, 881–904. [Google Scholar] [PubMed]
  35. Jiao, X.; Sherman, B.T.; Huang da, W.; Stephens, R.; Baseler, M.W.; Lane, H.C.; Lempicki, R.A. DAVID-WS: A stateful web service to facilitate gene/protein list analysis. Bioinformatics 2012, 28, 1805–1806. [Google Scholar] [CrossRef] [Green Version]
  36. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  37. Nivison-Smith, L.; Milston, R.; Madigan, M.; Kalloniatis, M. Age-related macular degeneration: Linking clinical presentation to pathology. Optom. Vis. Sci. 2014, 91, 832–848. [Google Scholar] [CrossRef] [PubMed]
  38. Vancura, P.; Csicsely, E.; Leiser, A.; Iuvone, P.M.; Spessert, R. Rhythmic Regulation of Photoreceptor and RPE Genes Important for Vision and Genetically Associated With Severe Retinal Diseases. Investig. Ophthalmol. Vis. Sci. 2018, 59, 3789–3799. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Fuhrmann, S.; Zou, C.; Levine, E.M. Retinal pigment epithelium development, plasticity, and tissue homeostasis. Exp. Eye Res. 2014, 123, 141–150. [Google Scholar] [CrossRef] [Green Version]
  40. D’Angelo, R.; Donato, L.; Venza, I.; Scimone, C.; Aragona, P.; Sidoti, A. Possible protective role of the ABCA4 gene c.1268A>G missense variant in Stargardt disease and syndromic retinitis pigmentosa in a Sicilian family: Preliminary data. Int. J. Mol. Med. 2017, 39, 1011–1020. [Google Scholar] [CrossRef] [Green Version]
  41. Donato, L.; Scimone, C.; Rinaldi, C.; Aragona, P.; Briuglia, S.; D’Ascola, A.; D’Angelo, R.; Sidoti, A. Stargardt Phenotype Associated With Two ELOVL4 Promoter Variants and ELOVL4 Downregulation: New Possible Perspective to Etiopathogenesis? Investig. Ophthalmol. Vis. Sci. 2018, 59, 843–857. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Wiktor, A.; Sarna, M.; Wnuk, D.; Sarna, T. Lipofuscin-mediated photodynamic stress induces adverse changes in nanomechanical properties of retinal pigment epithelium cells. Sci. Rep. 2018, 8, 17929. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Moiseyev, G.; Nikolaeva, O.; Chen, Y.; Farjo, K.; Takahashi, Y.; Ma, J.X. Inhibition of the visual cycle by A2E through direct interaction with RPE65 and implications in Stargardt disease. Proc. Natl. Acad. Sci. USA 2010, 107, 17551–17556. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Kaarniranta, K.; Machalinska, A.; Vereb, Z.; Salminen, A.; Petrovski, G.; Kauppinen, A. Estrogen signalling in the pathogenesis of age-related macular degeneration. Curr. Eye Res. 2015, 40, 226–233. [Google Scholar] [CrossRef]
  45. Akanuma, S.; Soutome, T.; Hisada, E.; Tachikawa, M.; Kubo, Y.; Hosoya, K. Na+-independent nucleoside transporters regulate adenosine and hypoxanthine levels in Muller cells and the inner blood-retinal barrier. Investig. Ophthalmol. Vis. Sci. 2013, 54, 1469–1477. [Google Scholar] [CrossRef] [Green Version]
  46. Dos Santos-Rodrigues, A.; Pereira, M.R.; Brito, R.; de Oliveira, N.A.; Paes-de-Carvalho, R. Adenosine transporters and receptors: Key elements for retinal function and neuroprotection. Vitam. Horm. 2015, 98, 487–523. [Google Scholar] [CrossRef]
  47. Fliesler, S.J.; Keller, R.K. Isoprenoid metabolism in the vertebrate retina. Int. J. Biochem. Cell Biol. 1997, 29, 877–894. [Google Scholar] [CrossRef]
  48. Roysommuti, S.; Thaeomor, A.; Khimsuksri, S.; Lerdweeraphon, W.; Wyss, J.M. Perinatal taurine imbalance alters the interplay of renin-angiotensin system and estrogen on glucose-insulin regulation in adult female rats. Adv. Exp. Med. Biol. 2013, 776, 67–80. [Google Scholar] [CrossRef]
  49. Okabe, K.; Kobayashi, S.; Yamada, T.; Kurihara, T.; Tai-Nagara, I.; Miyamoto, T.; Mukouyama, Y.S.; Sato, T.N.; Suda, T.; Ema, M.; et al. Neurons limit angiogenesis by titrating VEGF in retina. Cell 2014, 159, 584–596. [Google Scholar] [CrossRef] [Green Version]
  50. Joyal, J.S.; Gantner, M.L.; Smith, L.E.H. Retinal energy demands control vascular supply of the retina in development and disease: The role of neuronal lipid and glucose metabolism. Prog. Retin. Eye Res. 2018, 64, 131–156. [Google Scholar] [CrossRef]
  51. Preising, M.; Ayuso, C. Rab escort protein 1 (REP1) in intracellular traffic: A functional and pathophysiological overview. Ophthalmic Genet. 2004, 25, 101–110. [Google Scholar] [CrossRef] [PubMed]
  52. Kauppinen, A.; Paterno, J.J.; Blasiak, J.; Salminen, A.; Kaarniranta, K. Inflammation and its role in age-related macular degeneration. Cell Mol. Life Sci. 2016, 73, 1765–1786. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Piippo, N.; Korhonen, E.; Hytti, M.; Skottman, H.; Kinnunen, K.; Josifovska, N.; Petrovski, G.; Kaarniranta, K.; Kauppinen, A. Hsp90 inhibition as a means to inhibit activation of the NLRP3 inflammasome. Sci. Rep. 2018, 8, 6720. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Subrizi, A.; Toropainen, E.; Ramsay, E.; Airaksinen, A.J.; Kaarniranta, K.; Urtti, A. Oxidative stress protection by exogenous delivery of rhHsp70 chaperone to the retinal pigment epithelium (RPE), a possible therapeutic strategy against RPE degeneration. Pharm. Res. 2015, 32, 211–221. [Google Scholar] [CrossRef]
  55. Dhingra, A.; Alexander, D.; Reyes-Reveles, J.; Sharp, R.; Boesze-Battaglia, K. Microtubule-Associated Protein 1 Light Chain 3 (LC3) Isoforms in RPE and Retina. Adv. Exp. Med. Biol. 2018, 1074, 609–616. [Google Scholar] [CrossRef]
  56. Frost, L.S.; Lopes, V.S.; Bragin, A.; Reyes-Reveles, J.; Brancato, J.; Cohen, A.; Mitchell, C.H.; Williams, D.S.; Boesze-Battaglia, K. The Contribution of Melanoregulin to Microtubule-Associated Protein 1 Light Chain 3 (LC3) Associated Phagocytosis in Retinal Pigment Epithelium. Mol. Neurobiol. 2015, 52, 1135–1151. [Google Scholar] [CrossRef] [Green Version]
  57. Sparrow, J.R.; Zhou, J.; Cai, B. DNA is a target of the photodynamic effects elicited in A2E-laden RPE by blue-light illumination. Investig. Ophthalmol. Vis. Sci. 2003, 44, 2245–2251. [Google Scholar] [CrossRef] [Green Version]
  58. Wang, J.; Feng, Y.; Han, P.; Wang, F.; Luo, X.; Liang, J.; Sun, X.; Ye, J.; Lu, Y.; Sun, X. Photosensitization of A2E triggers telomere dysfunction and accelerates retinal pigment epithelium senescence. Cell Death Dis. 2018, 9, 178. [Google Scholar] [CrossRef] [Green Version]
  59. Cornish, E.E.; Natoli, R.C.; Hendrickson, A.; Provis, J.M. Differential distribution of fibroblast growth factor receptors (FGFRs) on foveal cones: FGFR-4 is an early marker of cone photoreceptors. Mol. Vis. 2004, 10, 1–14. [Google Scholar]
  60. Donato, L.; Bramanti, P.; Scimone, C.; Rinaldi, C.; Giorgianni, F.; Beranova-Giorgianni, S.; Koirala, D.; D’Angelo, R.; Sidoti, A. miRNAexpression profile of retinal pigment epithelial cells under oxidative stress conditions. FEBS Open Bio 2018, 8, 219–233. [Google Scholar] [CrossRef] [Green Version]
  61. Olejniczak, M.; Kotowska-Zimmer, A.; Krzyzosiak, W. Stress-induced changes in miRNA biogenesis and functioning. Cell Mol. Life Sci. 2018, 75, 177–191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Westlund, B.S.; Cai, B.; Zhou, J.; Sparrow, J.R. Involvement of c-Abl, p53 and the MAP kinase JNK in the cell death program initiated in A2E-laden ARPE-19 cells by exposure to blue light. Apoptosis 2009, 14, 31–41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. McCann, K.L.; Charette, J.M.; Vincent, N.G.; Baserga, S.J. A protein interaction map of the LSU processome. Genes Dev. 2015, 29, 862–875. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Marchesi, N.; Thongon, N.; Pascale, A.; Provenzani, A.; Koskela, A.; Korhonen, E.; Smedowski, A.; Govoni, S.; Kauppinen, A.; Kaarniranta, K.; et al. Autophagy Stimulus Promotes Early HuR Protein Activation and p62/SQSTM1 Protein Synthesis in ARPE-19 Cells by Triggering Erk1/2, p38(MAPK), and JNK Kinase Pathways. Oxid. Med. Cell Longev. 2018, 2018, 4956080. [Google Scholar] [CrossRef] [Green Version]
  65. McGrew, D.A.; Hedstrom, L. Towards a pathological mechanism for IMPDH1-linked retinitis pigmentosa. Adv. Exp. Med. Biol. 2012, 723, 539–545. [Google Scholar] [CrossRef]
Figure 1. MTT determination of A2E treatment in retinal pigment epithelium (RPE) cells. Cell death was assessed at considered time points (3 h and 6 h) in A2E treated samples compared to basal untreated group. Results are shown as mean ± standard error of mean (n = 3). p-value < 0.05.
Figure 1. MTT determination of A2E treatment in retinal pigment epithelium (RPE) cells. Cell death was assessed at considered time points (3 h and 6 h) in A2E treated samples compared to basal untreated group. Results are shown as mean ± standard error of mean (n = 3). p-value < 0.05.
Antioxidants 09 00307 g001
Figure 2. Summary figure of expressed genes and significant DE, DE + DAS, DE + DTU, DAS and DTU genes from analysis of the considered time points of RPE cell data. (A) Number of genes regulated only by transcription (DE), only by alternative splicing (DAS) and by both transcription and alternative splicing (DE + DAS); (B) number of transcripts regulated only by transcription (DE), only by alternative splicing (DAS) and by both transcription and alternative splicing (DE + DAS). DTU = Differential Transcript Usage. AS = Alternative Splicing.
Figure 2. Summary figure of expressed genes and significant DE, DE + DAS, DE + DTU, DAS and DTU genes from analysis of the considered time points of RPE cell data. (A) Number of genes regulated only by transcription (DE), only by alternative splicing (DAS) and by both transcription and alternative splicing (DE + DAS); (B) number of transcripts regulated only by transcription (DE), only by alternative splicing (DAS) and by both transcription and alternative splicing (DE + DAS). DTU = Differential Transcript Usage. AS = Alternative Splicing.
Antioxidants 09 00307 g002
Figure 3. Comparison of the gene lists generated during differential expression analyses. Euler diagrams of DE (A) and DAS (B) genes identified during expression analyses in considered conditions of experimental design, setting the contrast groups as 0 h.untreated versus 3 h.treated, 0 h.untreated versus 6 h.treated, 3 h.treated versus 6 h.treated, 0 h.untreated versus (3 h.treated + 6 h.treated)/2.
Figure 3. Comparison of the gene lists generated during differential expression analyses. Euler diagrams of DE (A) and DAS (B) genes identified during expression analyses in considered conditions of experimental design, setting the contrast groups as 0 h.untreated versus 3 h.treated, 0 h.untreated versus 6 h.treated, 3 h.treated versus 6 h.treated, 0 h.untreated versus (3 h.treated + 6 h.treated)/2.
Antioxidants 09 00307 g003
Figure 4. Hierarchical clustering and heatmap of DE genes and Key GO terms. DE genes show segregation into 10 coexpressed clusters, of which the main ones related to oxidative stress were highlighted (circled) and linked to GO specific terms. Full results of GO enrichment analyses are shown in Tables S6 and S7. The z-score scale represents mean-subtracted regularized log-transformed transcripts per million (TPMs).
Figure 4. Hierarchical clustering and heatmap of DE genes and Key GO terms. DE genes show segregation into 10 coexpressed clusters, of which the main ones related to oxidative stress were highlighted (circled) and linked to GO specific terms. Full results of GO enrichment analyses are shown in Tables S6 and S7. The z-score scale represents mean-subtracted regularized log-transformed transcripts per million (TPMs).
Antioxidants 09 00307 g004
Figure 5. Hierarchical clustering and heatmap of DE genes + DTU transcripts from DAS genes and Key GO terms. DE genes and DTU transcripts from DAS genes show segregation into 10 coexpressed clusters, of which the main ones related to oxidative stress were highlighted (circled) and linked to GO specific terms. Full results of GO enrichment analyses are shown in Tables S6 and S7. The z-score scale represents mean-subtracted regularized log-transformed TPMs.
Figure 5. Hierarchical clustering and heatmap of DE genes + DTU transcripts from DAS genes and Key GO terms. DE genes and DTU transcripts from DAS genes show segregation into 10 coexpressed clusters, of which the main ones related to oxidative stress were highlighted (circled) and linked to GO specific terms. Full results of GO enrichment analyses are shown in Tables S6 and S7. The z-score scale represents mean-subtracted regularized log-transformed TPMs.
Antioxidants 09 00307 g005
Figure 6. Volcano plots of significant DE and DAS genes and of DE and DTU transcripts. Volcano plots of significant (adjusted p-value < 0.01) DE genes (A), DE transcripts (B), DAS genes (C) and DTU transcripts (D). The low expressed genes and transcripts were filtered. The top 15 considered elements with the smallest corrected p-values and bigger fold-changes are highlighted, and different colors refer to different contrast groups. DE genes: log2FC vs. −log10(FDR) at gene level; DAS genes: maximum ΔPS of transcript in a gene vs. −log10(FDR) at gene level; DE transcripts: log2FC vs. −log10(FDR) at gene level at transcript level and DTU transcripts: ΔPS vs. −log10(FDR) at transcript level.
Figure 6. Volcano plots of significant DE and DAS genes and of DE and DTU transcripts. Volcano plots of significant (adjusted p-value < 0.01) DE genes (A), DE transcripts (B), DAS genes (C) and DTU transcripts (D). The low expressed genes and transcripts were filtered. The top 15 considered elements with the smallest corrected p-values and bigger fold-changes are highlighted, and different colors refer to different contrast groups. DE genes: log2FC vs. −log10(FDR) at gene level; DAS genes: maximum ΔPS of transcript in a gene vs. −log10(FDR) at gene level; DE transcripts: log2FC vs. −log10(FDR) at gene level at transcript level and DTU transcripts: ΔPS vs. −log10(FDR) at transcript level.
Antioxidants 09 00307 g006
Figure 7. Expression profiles of most important DE genes/transcripts; first cluster associated with newly discovered pathways linked to oxidative stress. Detailed gene/transcript expression profiles across the time course of the first cluster of most important DE genes (ACTG1, CAPZB, CCN2, FTL and RPS11) linked to newly identified candidate pathways linked to oxidative stress.
Figure 7. Expression profiles of most important DE genes/transcripts; first cluster associated with newly discovered pathways linked to oxidative stress. Detailed gene/transcript expression profiles across the time course of the first cluster of most important DE genes (ACTG1, CAPZB, CCN2, FTL and RPS11) linked to newly identified candidate pathways linked to oxidative stress.
Antioxidants 09 00307 g007
Figure 8. Expression profiles of most important DE genes/transcripts secondo cluster associated with newly discovered pathways linked to oxidative stress. Detailed gene/transcript expression profiles across the time course of the second cluster of most important DE genes (P4HB, RNA5-8SN2, RPL3 and RPL19) linked to newly identified candidate pathways linked to oxidative stress.
Figure 8. Expression profiles of most important DE genes/transcripts secondo cluster associated with newly discovered pathways linked to oxidative stress. Detailed gene/transcript expression profiles across the time course of the second cluster of most important DE genes (P4HB, RNA5-8SN2, RPL3 and RPL19) linked to newly identified candidate pathways linked to oxidative stress.
Antioxidants 09 00307 g008
Table 1. ToppGene prioritization analysis results. The ToppGene prioritization analysis on known causative retinal dystrophy genes that intersected with most significant DE and DAS genes obtained during the whole transcriptome analysis revealed 19 candidate genes (Bonferroni corrected p-value < 0.05), of which seven had a strong association with training genes (Bonferroni-corrected p-value < 0.01). All shown p-values are Bonferroni-corrected.
Table 1. ToppGene prioritization analysis results. The ToppGene prioritization analysis on known causative retinal dystrophy genes that intersected with most significant DE and DAS genes obtained during the whole transcriptome analysis revealed 19 candidate genes (Bonferroni corrected p-value < 0.05), of which seven had a strong association with training genes (Bonferroni-corrected p-value < 0.01). All shown p-values are Bonferroni-corrected.
RankGene SymbolGeneIdGO: Mol. Func. ScoreGO: Mol. Func. p-ValueGO: Bio. Proc. ScoreGO: Bio. Proc. p-ValueGO: Cell. Comp. ScoreGO: Cell. Comp. p-ValueHuman Pheno. ScoreHuman Pheno. p-ValuePathway ScorePathway p-ValuePubmed ScorePubmed p-ValueDisease ScoreDisease p-ValueAverage ScoreOverall p-Value
1PAX650800.04220.09311.00.01020.28140.08071.00.01350.00.5031.00.0010.9997.75 × 10−100.61761 88 × 10−12
2ACTG1710.81430.01060.9960.02780.96860.01270.9990.01350.00.5030.3300.0320.7300.0060.72675 15 × 10−11
3TGFBI70450.63290.02880.9820.04620.46830.05340.9880.03310.00.5030.2040.0730.8190.0030.62260.002
4CCN214900.43180.04910.9990.02090.43620.05440.9990.02530.00.5030.2040.0730.5510.0110.55290.003
5CTSH15120.27960.06770.9990.01850.99860.00360.9760.03580.00.5030.00.5360.4450.0160.53020.003
6GNAI227710.92630.00290.9620.05730.95000.01660.9990.01950.00.5030.0950.0730.00.5120.59510.004
7LTBP240530.59180.03750.8130.10490.03500.22540.9990.01350.00.5030.2040.0730.6360.0090.51680.007
8HNRNPA131780.04210.09310.9100.07720.54810.04990.9990.02140.00.5030.8010.0060.00.5120.49130.011
9GPI28210.04210.09310.9650.05630.95630.01560.9860.03350.00.5030.4900.0190.00.5120.52490.011
10CD819750.67020.02130.9990.02090.30570.07820.9990.02320.00.5030.0950.0730.00.5120.46910.011
11FTL25120.26220.07660.5870.14030.23700.09790.9990.02050.00.5030.0950.0730.3310.0210.41510.012
12ITGAV36850.47320.04760.9990.01560.98850.0094−1.00.00.00.5030.0760.0730.00.5120.45630.014
13CAPZB8320.50610.04260.9980.02420.99340.0079−1.00.00.00.5030.2050.0730.00.5120.49190.015
14LTBP340540.04210.09310.9820.04580.03500.22540.9990.01780.00.5030.00.5360.6360.0090.38640.018
15CTSB15080.29960.06420.9150.07530.38620.06000.1830.05530.00.5030.3430.0320.00.5120.37940.026
16P4HB50340.29960.06420.8510.09600.23700.09790.9410.04000.00.5030.3510.0320.00.5120.43170.031
17MATN241470.47470.04640.9870.04080.37960.0606−1.00.00.00.5030.2040.0730.00.5120.38340.042
18TMEM189-UBE2V1387522−1.00.0−1.00.0−1.00.0−1.00.00.00.5030.0590.0730.5790.0110.27460.050
19ACADVL370.62720.02900.9030.07920.24220.08570.9800.03480.00.5030.00.5360.00.5120.40140.053
20CD1519770.30150.06390.7760.11130.44890.05420.9940.03130.00.5030.00.5360.00.5120.37550.065
21CAPNS18260.04210.09310.8190.10390.14280.1432−1.00.00.00.5030.0760.0730.00.5120.22770.103
22CD639670.00.56520.9860.04160.33900.0613−1.00.00.00.5030.0950.0730.00.5120.27350.104
23RNA5-8SN21098642810.32500.0559−1.00.00.24660.0846−1.00.0−1.00.00.00.536−1.00.00.19050.112
24RPL361220.32500.05590.9520.06090.28140.0807−1.00.00.00.5030.00.5360.00.5120.26370.117
25RPS1162050.32500.05590.8770.08610.28140.0807−1.00.00.00.5030.00.5360.00.5120.25710.132
26RPL1961430.32500.05590.8770.08610.28140.0807−1.00.00.00.5030.00.5360.00.5120.25710.132
27SLC16A391230.00.56520.2120.20390.30570.0782−1.00.00.00.5030.0650.0730.00.5120.17030.192
28UTP14C97240.00.56520.2200.20370.28140.0807−1.00.00.00.5030.00.5360.00.5120.12510.353
29TMSB4XP67120−1.00.0−1.00.0−1.00.0−1.00.0−1.00.00.00.536−1.00.00.01.0
Table 2. Detailed analysis of DE and DAS gene-involved pathways, with their possible impact on retinal dystrophies etiopathogenesis. DE and DAS genes were dysregulated during the whole analysis and showed several fluctuations during observed time points, suggesting that changes in gene-level expression and alternative splicing occurred throughout the whole period, either transiently (occurring after 3 h and returning to initial level after 6 h), occurring later (only after 6 h from treatment) or enduring throughout the whole period.
Table 2. Detailed analysis of DE and DAS gene-involved pathways, with their possible impact on retinal dystrophies etiopathogenesis. DE and DAS genes were dysregulated during the whole analysis and showed several fluctuations during observed time points, suggesting that changes in gene-level expression and alternative splicing occurred throughout the whole period, either transiently (occurring after 3 h and returning to initial level after 6 h), occurring later (only after 6 h from treatment) or enduring throughout the whole period.
De Gene-Involved PathwaysExpression ChangesDas Gene-Involved Pathways
RNA methyltransferase3 h and 6 h = DOWN-REGULATEDEndosomal sorting complex required for transport (ESCRT) and RAB geranylgeranylation
TRAF6 mediated NF-kB activation and activation of IKK by MEKK1Phagopore assembly site membrane; C-terminal protein lipidation; protein localization to microtubule cytoskeleton; regulation of TNFR1 signaling; TNF signaling
Condensed chromosome outer kinetochore and kinesin complexMethylation; activation of chaperone genes; nucleotide-sugar biosynthetic process
Transport of nucleoside and free purine and pyrimidine; histone pre-mRNA DCP binding; formation of AT-AC complex; respiratory electron transport3 h and 6 h = UP-REGULATEDNegative regulation of FGFR1 signaling
Regulation of telomerase RNA localization to Cajal body; maturation of LSU-rRNA; miRNAs involved in DNA damage response
Mithocondrial intermembrane space and TP53 regulates transcription of cell death genes
Polysomal ribosome
Mevalonate pathway and cholesterol biosynthesis3 h = UP-REGULATED
6 h = DOWN-REGULATED
Cholesterol biosynthesis
CBL binds and ubiquinates Sprouty and MAP kinase phosphatase activity3 h = DOWN-REGULATED
6 h = UP-REGULATED
/
PTK2/SRC-1 phosphorylates BCAR1

Share and Cite

MDPI and ACS Style

Donato, L.; D’Angelo, R.; Alibrandi, S.; Rinaldi, C.; Sidoti, A.; Scimone, C. Effects of A2E-Induced Oxidative Stress on Retinal Epithelial Cells: New Insights on Differential Gene Response and Retinal Dystrophies. Antioxidants 2020, 9, 307. https://doi.org/10.3390/antiox9040307

AMA Style

Donato L, D’Angelo R, Alibrandi S, Rinaldi C, Sidoti A, Scimone C. Effects of A2E-Induced Oxidative Stress on Retinal Epithelial Cells: New Insights on Differential Gene Response and Retinal Dystrophies. Antioxidants. 2020; 9(4):307. https://doi.org/10.3390/antiox9040307

Chicago/Turabian Style

Donato, Luigi, Rosalia D’Angelo, Simona Alibrandi, Carmela Rinaldi, Antonina Sidoti, and Concetta Scimone. 2020. "Effects of A2E-Induced Oxidative Stress on Retinal Epithelial Cells: New Insights on Differential Gene Response and Retinal Dystrophies" Antioxidants 9, no. 4: 307. https://doi.org/10.3390/antiox9040307

APA Style

Donato, L., D’Angelo, R., Alibrandi, S., Rinaldi, C., Sidoti, A., & Scimone, C. (2020). Effects of A2E-Induced Oxidative Stress on Retinal Epithelial Cells: New Insights on Differential Gene Response and Retinal Dystrophies. Antioxidants, 9(4), 307. https://doi.org/10.3390/antiox9040307

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop