Next Article in Journal
Navigating the Multiverse of Antisense RNAs: The Transcription- and RNA-Dependent Dimension
Next Article in Special Issue
Reassessing the Abundance of miRNAs in the Human Pancreas and Rodent Cell Lines and Its Implication
Previous Article in Journal
Identification of Differentially Expressed Intronic Transcripts in Osteosarcoma
Previous Article in Special Issue
In Silico Identification and Characterization of circRNAs as Potential Virulence-Related miRNA/siRNA Sponges from Entamoeba histolytica and Encystment-Related circRNAs from Entamoeba invadens
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Oxidative Stress and Its Modulation by Ladostigil Alter the Expression of Abundant Long Non-Coding RNAs in SH-SY5Y Cells

Department of Biological Chemistry, The Institute of Life Sciences, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
*
Author to whom correspondence should be addressed.
Non-Coding RNA 2022, 8(6), 72; https://doi.org/10.3390/ncrna8060072
Submission received: 25 September 2022 / Revised: 11 October 2022 / Accepted: 19 October 2022 / Published: 25 October 2022
(This article belongs to the Collection Feature Papers in Non-Coding RNA)

Abstract

:
Neurodegenerative disorders, brain injury, and the decline in cognitive function with aging are accompanied by a reduced capacity of cells in the brain to cope with oxidative stress and inflammation. In this study, we focused on the response to oxidative stress in SH-SY5Y, a human neuroblastoma cell line. We monitored the viability of the cells in the presence of oxidative stress. Such stress was induced by hydrogen peroxide or by Sin1 (3-morpholinosydnonimine) that generates reactive oxygen and nitrogen species (ROS and RNS). Both stressors caused significant cell death. Our results from the RNA-seq experiments show that SH-SY5Y cells treated with Sin1 for 24 h resulted in 94 differently expressed long non-coding RNAs (lncRNAs), including many abundant ones. Among the abundant lncRNAs that were upregulated by exposing the cells to Sin1 were those implicated in redox homeostasis, energy metabolism, and neurodegenerative diseases (e.g., MALAT1, MIAT, GABPB1-AS1, NEAT1, MIAT, GABPB1-AS1, and HAND2-AS1). Another group of abundant lncRNAs that were significantly altered under oxidative stress included cancer-related SNHG family members. We tested the impact of ladostigil, a bifunctional reagent with antioxidant and anti-inflammatory properties, on the lncRNA expression levels. Ladostigil was previously shown to enhance learning and memory in the brains of elderly rats. In SH-SY5Y cells, several lncRNAs involved in transcription regulation and the chromatin structure were significantly induced by ladostigil. We anticipate that these poorly studied lncRNAs may act as enhancers (eRNA), regulating transcription and splicing, and in competition for miRNA binding (ceRNA). We found that the induction of abundant lncRNAs, such as MALAT1, NEAT-1, MIAT, and SHNG12, by the Sin1 oxidative stress paradigm specifies only the undifferentiated cell state. We conclude that a global alteration in the lncRNA profiles upon stress in SH-SY5Y may shift cell homeostasis and is an attractive in vitro system to characterize drugs that impact the redox state of the cells and their viability.

1. Introduction

A decline in cognitive function occurs along with brain aging. In humans and rodents, it includes structural changes and extensive neuronal cell death [1]. The young brain successfully copes with stress (e.g., oxidative stress), a capacity that drops in the aging brain [2,3]. Failure in the stress response leads to pathological processes, as documented in brain injury, ischemia and reperfusion (I/R), as well as neurodegenerative diseases [4,5].
Gradual changes in the gene expression profiles occur throughout the progression of normal aging and neurodegenerative diseases [6,7]. In previous studies, ladostigil, a bifunctional reagent with antioxidant and anti-inflammatory activities, was shown as a novel neuroprotective drug [8]. Old rats that were chronically treated with ladostigil showed improved performance in memory and learning tasks relative to untreated old rats [9,10,11]. Nevertheless, the direct effect of oxidative stress at a cellular level is masked by the communication and feedback circuits of different cell types, including the microglia, astrocytes, and endothelial cells with neurons. To overcome this difficulty, the antioxidant function of ladostigil was studied using a neuroblastoma cell line [8,12]. Upon oxidative stress, an imbalance in mitochondrial redox was monitored, which was attenuated by ladostigil [12].
Emerging evidence suggests the involvement of long non-coding RNAs (lncRNAs) in neurodevelopmental, metabolic, and inflammatory diseases [13]. Furthermore, specific lncRNAs are implicated in all major neurodegenerative diseases, including Alzheimer’s disease (AD) [14,15]. Specifically, the highly expressed lncRNAs BACE1-AS, 51A, and BC200 were implicated in the production and accumulation of Aβ and were implicated in AD pathology in mice with AD pathology; SOX21-AS1 expression determines the level of neuronal oxidative stress and suppresses neuronal apoptosis [16]. Other lncRNAs have been linked to Parkinson’s disease (PD) [17,18]. For example, HOTAIR induces the apoptosis of dopaminergic neurons and MALAT1 (also called NEAT2) is associated with α-synuclein expression and PD pathogenesis [19]. The role of lncRNAs in governing the cellular oxidative levels was investigated beyond the context of neurodegenerative diseases and brain injury. Examples include oxidative stress that impacts pancreatic functioning [20], cardiac myocytes [21], lung epithelia [16], and more. Multiple layers of regulation are involved in lncRNA modulation, including chromatin remodeling [20], the regulation of gene expression enhancers [22], competition for miRNA binding sites on mRNAs [23], the recruitment of transcription factors to promoters, alterations in splicing, interference with expression as antisense transcripts, molecular decoys, and more [24,25]. Moreover, the involvement of circRNAs in manipulating oxidative stress has been established in cell systems and in the pathogenesis of diseases of the CNS [26]. Studies have shown a correlation between oxidative stress and lncRNA expression [27]. However, the mechanism by which lncRNAs can change in response to drugs that alter cellular homeostasis is mostly unknown.
In this study, we studied SH-SY5Y cells, a neuroblastoma-derived cell line that carries some dopaminergic neuronal features [28]. We tested the expression profiles of the cells under varying conditions by RNA-seq [12]. We exposed the SH-SY5Y cells to oxidative stress and monitored their viability and gene expression changes [12]. We applied conditions that resembled abrupt (by H2O2) and prolonged (by Sin1; 3-morpholinosydnonimine) oxidative stress and monitored the beneficial effect of ladostigil on the viability of SH-SY5Y cells under stress conditions. We focused on the mild upregulation of abundant lncRNAs following treatment with Sin1 and determined the impact of ladostigil under oxidative stress conditions on lncRNA expression. Our investigation of the gene expression levels along with the attenuation of cell death by ladostigil suggests the beneficial effect of lncRNAs on cell viability while coping with oxidative stress.

2. Results

2.1. Ladostigil Suppresses Cell Death

Figure 1A shows the results of flow cytometry to determine the fraction of necrotic cells (each experiment covered 50,000 cells). Propidium iodide (PI) was used as a marker for permeable membranes of dead cells. After 5 h of exposure to 40 μM H2O2, we saw a significant increase in PI-positive cells. In the presence of ladostigil, the fraction of PI-positive cells was attenuated (8.0% to 5.4%; p-value = 0.02). The attenuation in cell death was also detected 24 h later (Supplementary Figure S1). After treatment with ladostigil, the fraction of PI-positive cells in cells exposed to 40 μM and 80 μM H2O2 was reduced (Supplementary Figure S1). Therefore, we conclude that ladostigil protects cells against non-reversible cell death induced by acute and short-lived stress of H2O2.

2.2. Apoptotic Signal in SH-SY5Y Cells Exposed to Sin1

To better mimic a prolonged and steady increase in oxidation, we exposed the cells to Sin1 at varying concentrations. Sin1 in living cells acts as a peroxynitrite donor (ONOO) and generates nitric oxide (NO) and superoxide (O2) [29,30]. In addition to the effect of Sin1 on the redox state of cells, it was reported as an inhibitor of high-affinity choline transporters and trafficking in SH-SY5Y cells [31]. We monitored cell viability (using MTT) by increasing the levels of Sin1 (0 to 350 mM). In parallel, we performed the same experiment but with the preincubation of ladostigil (5.4 μM, 2 h prior to Sin1 exposure). Calculating the IC50 for each of the graphs resulted in the IC50 for Sin1 being 247.0 μM, while the ICD50 of ladostigil in the presence of Sin1 changed to 350.4 μM. These results allow the beneficial effect of ladostigil in coping with Sin1-induced cellular stress to be quantified (Figure 1B). The effect of prolonged oxidative stress produced by high levels of Sin1 on cell viability was assessed by flow cytometry using PI (Figure 1C). Cell viability was tested at 5 h upon increasing the concentration of Sin1. Under these harsh conditions, we also stained SH-SY5Y cells with Annexin-V, a sensitive marker for early apoptosis. Staining cells with Annexin-V showed that the fraction of apoptotic cells was insensitive to increased levels of Sin1 (Figure 1D). We concluded that, under a high concentration of Sin1, cell death was increased. However, the fraction that underwent apoptosis remained stable (~10%), supporting an increase in necrotic cells, rather than apoptosis. Recall that, at later time points, separating cells by PI and Annexin-V positives is no longer valid, as the membranes of the cells become fully permeable (discussed in [32]).

2.3. LncRNA Expression Levels in SH-SY5Y Cells

Sin1 at a low concentration (50 μM) produces mild but long-lasting stress due to the ongoing accumulation of nitric oxide (NO) and superoxide (O2). We focused on the ncRNA levels using RNA-seq (RNA length cutoff at >200 nt) 24 h after treatment. Figure 2A shows that the cells’ transcriptome comprised 70% mRNAs of coding genes (13,626 coding mRNAs out of a total of 19,475 transcripts), with the rest of the transcripts being partitioned among many ncRNA biotypes. The dominant class was lncRNA (long ncRNA, 63.2% of all ncRNAs), followed by pseudogenes (28.8%). However, the expression levels measured by TMM for most of these transcripts were very low (<4 TMM for biological triplicates, see the Section 4). By considering only transcripts that satisfied this threshold, we focused on 780 ncRNAs expressing transcripts that accounted for 6.7% of all transcripts (total of 11,588, Figure 2B). The rest of the analysis will consider only these ncRNAs (Supplementary Table S1). A subdivision according to biotypes confirmed that most ncRNA transcripts belonged to the lncRNA group (73%), followed by the unprocessed pseudogenes group (19%, Figure 2B).
Figure 2C shows a sorted list of high-expressing lncRNAs in SH-SY5Y cells (above TMM of 100). A manual classification of these highly expressed transcripts into groups revealed that about 16% of these transcripts belonged to mitochondrial ncRNAs (e.g., MT-RNR1 and MT-RNR2 for 12S and 16S rRNAs of the mitochondrial ribosome, respectively), and an additional 10% were transcripts involved in X-chromosomal inhibition and chromatin compaction (e.g., XIST, FTX, and JPX, collectively termed XCI, Figure 2C). Several of the other abundant lncRNAs (Figure 2C) were implicated in the regulation of oxidative stress, particularly in the context of the CNS and neurons (e.g., MALAT1, GAS5, and NEAT1) [33].

2.4. Chronic Oxidative Stress Upregulated a Collection of LncRNAs

Studies have shown a dependency between acute oxidative stress and the upregulation of lncRNA expression [27,34]. We tested the profile of the cells’ ncRNA levels 24 h after the Sin1 stimulus. Figure 3A shows the number of differentially expressed (DE) ncRNA transcripts partitioned into ncRNA types. We found that the strongest effect induced by Sin1 was the upregulation of ncRNAs. The fraction of ncRNAs among DE downregulated genes was only 5.6%, while among the upregulated transcripts, it accounted for 17.1% (p-value = 0.012). Moreover, among the downregulated genes, the ratio of lncRNAs to all other molecular ncRNA biotypes was 2.1, but it was 5.1 for the upregulated genes (Figure 3A). We conclude that the fraction of the expressed ncRNAs, and specifically the lncRNA class, significantly increased under Sin1-induced stress. While many coding genes were induced by Sin1 [12], the ncRNA expression fold change was modest (bounded by <3-fold, Figure 3B).
Table 1 shows the list of lncRNAs following Sin1 treatment (listed genes with >0.5% of total expression, see the Section 4). All 18 Sin1-upregulated ncRNAs (Table 1) were statistically significant (p-value FDR < 1.0 × 10−3; Supplementary Table S2). The highest-expressed lncRNA was MALAT1 (Figure 2C), which was also upregulated by Sin1 (Table 1, p-value FDR < 2.5 × 10−4).
Table 1 presents the Sin1-induced lncRNAs that met the threshold of absolute expression and that were statistically significant. Among them were NEAT1 (upregulated 1.49-fold, FDR p-value 3.5 × 10−5) and HAND2-AS1 (upregulated 1.33-fold, FDR p-value 6.7 × 10−25), both implicated in oxidative stress modulation. Notably, among the Sin1-upregulated genes, many were annotated as antisense (21% of all transcripts upregulated by Sin1; Supplementary Table S3). Table 1 also lists lncRNAs that were downregulated by Sin1 (five transcripts, Supplementary Table S3). Note that, in comparison with the upregulated lncRNAs, the absolute expression levels of these genes were low (<40 TMM).

2.5. Cancer-Related SNHG Family Members Are Differentially Expressed by Sin1

SNHGs (small nucleolar RNA host genes) belong to a large family of lncRNAs whose expression greatly fluctuates among cell types and conditions [35]. We observed 17 SNHGs expressed transcripts, which accounted for 3.3% of all highly expressed lncRNAs (Supplementary Table S4, 780 genes).
Figure 4A presents the expression levels (log10TMM) of the identified SNHGs. Four of the SNHGs (SNHG5, SNHG14, SNHG16, and SNHG29) were highly expressed, with >100 TMM each. The fold change in expression following Sin1 treatment was small to modest. SNHG3, SNHG4, SNHG12, and SNHG15 were significantly downregulated by Sin1. SNHG29 and SNHG5, which were highly expressed lncRNAs (TMM 328.3 and 131.4, respectively), were slightly upregulated by 10.1% and 8.1%, respectively (FDR p-value of 1.6 × 10−4 and 0.029; Figure 4, Supplementary Table S4). Our observations are consistent with the proposed role of SNHGs in miRNA competition and as molecular decoys in SH-SY5Y [36].

2.6. Impact of Ladostigil on ncRNA Profiles

To test the potential of ladostigil to cope with prolonged mild stress by Sin1 (50 μM, 24 h), we applied ladostigil (5.4 mM) and measured the effect on gene expression 24 h after treatment. Under such mild conditions of oxidation, the level of cell death was minimal (below 10%) [12]. As a control, we incubated cells with ladostigil for 2 h, which had negligible effects on cell viability (24 h post-incubation, Supplementary Table S1).
Figure 5A shows four of the most significantly upregulated lncRNAs by Sin1 (p-value FDR ranged from 1 × 10−23 to 1 × 10−51), but were unchanged in the presence of ladostigil (24 h). Among these transcripts were TSPEAR-AS1 and TSPEAR-AS2. TRPEAR is a low-expressing protein that participates in the regulation of the Notch signaling pathway. While information on antisense transcripts is scarce, samples of kidney renal clear cell carcinoma with higher TSPEAR-AS1 expression showed poor survival when compared with patients with low expression (hazard ratio of 1.95; p-value of 8.8 × 10−6) [37]. The other examples, AC106881.1 (UNC5C antisense RNA 1; Figure 5A) and BX005019.1, are exclusively expressed in brain tissues. BX005019.1 is located within the MIR137 host gene, a main player in the differentiation and maturation of the nervous system [38].
We tested whether ladostigil reduces the damage induced by oxidative stress via lncRNA expression. To this end, we considered the expression profile following Sin1 as a reference, and applied a strict statistical threshold. Only 8 genes out of 19,475 annotated genes were identified as significantly induced by ladostigil, with 4 coding genes and 4 lncRNAs (Figure 5B). RPARP-AS1 (also called C10orf95 Antisense RNA 1) overlapped with the production of a short (213 aa) protein of unknown function (Figure 5C). Possible regulation of the splicing machinery was proposed for AC105339.2, which overlapped with the FSD2 transcript (fibronectin type III and SPRY domain containing 2; also called A1L4K1). SNHG21 and SCARNA15 (Figure 5C) were two other overlapping genes that were implicated in the redox homeostasis of cancer cells and stress conditions [39].
We further inspected the genomic organization of these four ncRNAs upregulated by ladostigil in cells exposed to Sin1 (24 h). We used GeneHancer (GH) [40] to test the likelihood of the ladostigil-upregulated ncRNAs acting as enhancers (eRNAs). GeneHancer integrates information from numerous independent databases.
Table 2 presents the prediction confidence for the four ladostigil-upregulated ncRNAs as enhancers in view of their potential gene target regulation. Elite enhancers were defined by high confidence scores for enhancers according to multiple independent resources. All of the genes listed were classified as elite enhancers. Note that only 7% of the >100k analyzed genes and transcripts in the GH database were marked as elite enhancers. GH double elite was assigned to two of the four transcripts. Under stress, ladostigil induction of lncRNAs may impact global processes, as reflected by eRNAs that act as enhancers.
It is an accepted view that each lncRNA molecule is involved in multiple cellular functions. For example, we propose a function for AC105339.2 as an enhancer (eRNA, Table 2) and also as an antisense (Figure 5). Moreover, the human RNA interactome has a wealth of experimental evidence for protein, miRNA, and lncRNA interactions [41]. The notion of multi-function activities for lncRNAs was validated for several abundant ones. Both GAS5 and the mitochondrial pseudogene MTND4P12 were among the most abundant lncRNAs in SH-SY5Y (Figure 2C), and both were significantly upregulated in the presence of ladostigil (Supplementary Table S2). GAS5 is implicated in human cancer as a tumor suppressor that is modulated by miRNAs. In SH-SY5Y cells exposed to Sin-1, the level of GAS5 is already high (TMM = 188) and is slightly upregulated by ladostigil (p-value 5 × 10−5). We report on 45 MBSs on GAS5 (mapped to 35 unique miRNAs, Supplementary Table S5) using a stringent seed-matching scheme [42]. Interestingly, the GAS5 MBSs contained several highly expressed miRNAs, such as miR-10b (21.8% of all expressed miRNAs in undifferentiated SH-SY5Y), miR-21, miR-148a, miR-26a, miR-222, and more. We expect that the increase in GAS5 transcripts by ladostigil will directly bind miRNAs and consequently impact cell survival through ceRNA activity (Supplementary Table S5). A similar activity of ceRNA can be performed by MTND4P12 (MT-ND4 pseudogene 12), which has MBSs for 10 different miRNAs, among which miR-140 and miR-24 are relatively abundant in SH-SY5Y cells.

2.7. A Shift in Abundant ncRNAs upon Stress Is Cell-State-Dependent

SH-SY5Y cells can undergo in vitro differentiation with retinoic acid (RA), and as a result, develop a predominantly mature dopaminergic-like neuronal phenotype [43]. We performed semi-quantitative RT-PCR to test the impact of gene expression change for several of the major ncRNAs that were implicated in the response to oxidation stress according to the cell state (undifferentiated vs. RA differentiated). We assessed the effect on gene expression following the exposure of cells to Sin1 only and to Sin1 in the presence of ladostigil (5.4 μM; Figure 6). The results from the gel analysis of the amplicons of the lncRNAs indicate that the differentiated cells were signified by reduced expressions of MALAT-1 and MIAT, but not NETA-1. In accordance with the results from the RNA-seq analysis, the expression levels of MALAT-1, MIAT, and NEAT-1 were induced by Sin1, and ladostigil partially reverted the elevated expression of these major lncRNAs (no change was observed in β-actin as an internal control; Figure 6A). The maximal upregulation in expression was measured in MIAT (lncRNA myocardial infarction-associated transcript). The MIAT was shown to maintain human lens epithelial cells (HLECs), where its expression is critical to the viability, proliferation, and migration of the cells in response to oxidative stress [27]. Overall, Sin1 induced many of the abundant ncRNAs in the undifferentiated cells, while ladostigil attenuated the effect. In contrast, the effect of Sin1 was quite modest and no induction was observed in undifferentiated cells.

3. Discussion

In this work, we examined undifferentiated SH-SY5Y cells under oxidative stress conditions by inspecting the lncRNA expression profiles [44,45,46]. In brain injury, ischemia, neuroinflammation, and other pathophysiological processes, changes in the amount and composition of ncRNAs (e.g., miRNAs and lncRNAs) are frequently observed [47,48]. Several discoveries regarding ncRNAs emerged from the analysis of SH-SY5Y under continuous oxidative stress. Firstly, although ncRNAs accounted for 30% of all transcripts detected, many were lowly expressed (only 6.7% with TMM > 4, Figure 2), reflecting the high sensitivity of the sequencing methodology [49]. Secondly, following sin1, most lncRNAs were upregulated, with only a few being downregulated (Figure 3, Fisher exact test, p-value < 1 × 10−5). Finally, the absolute fold change (FC) of the Sin1-upregulated lncRNAs was low (Table 1). We argue that, even with such a low increase or decrease in the overall expression (Figure 3B bounded by 3.0-fold), the dysregulation of abundant lncRNAs substantially impacts cells’ states.
Despite the growing interest in using lncRNAs for prognosis [50], knowledge of lncRNA binding specificity is limited, with only a few highly expressed lncRNAs having been studied in depth. Even under the mild oxidation stimulation that was introduced to undifferentiated SH-SY5Y cells (50 μM Sin1, 24 h), tens of lncRNAs were up-regulated, a few of them belonging to the highly expressed lncRNAs. Many of the highly expressed lncRNAs have been linked to neurodegenerative, inflammatory, metabolic, and systemic diseases [33,51]. It is likely that the same lncRNAs may act on multiple layers of regulation, such as chromatin, transcription, biological machines (e.g., spliceosomes), recruiting or stabilizing proteins, and competing for binding sites [52]. Inspecting the genomic location supports the notion of multiple functions that can be carried out by a specific lncRNA [53]. For example, GABPB1-AS1 (Table 1) suppresses the transcription of GABPB1, a gene product that acts upon damaged membranes via peroxidases. As a result, in both health and disease, it regulates cellular antioxidant capacity and cell viability [15]. Several of the most abundant lncRNAs that were affected by changing conditions (e.g., MALAT1, XIST, GAS5, and NEAT1) [24] have been proposed to act via multiple modes of action [52]. NEAT1 lncRNA, in particular, was initially linked to chromatin organization but has since been shown to play a role in regulating oxidative stress, pain control, and inflammation [54,55]. In neurons, NEAT1 could reverse the damage caused by superoxide and counteract the H2O2-induced neuronal damage. The dysregulation of HAND2-AS1 can alter energy metabolism at the organismal level [56]. MALAT1, the most abundant lncRNA in SH-SY5Y, is involved in a variety of cell biology processes, including autophagy and ischemia-induced apoptosis [57,58]. MALAT1 and NEAT1 play a neuroprotective role in AD and PD, respectively [59,60]. MIAT, NEAT1, and MALAT1 are coordinately changed in schizophrenia [61]. Other lncRNAs have been linked to autism, amyotrophic lateral sclerosis (ALS), and other diseases.
The direct role of miRNAs in maintaining SH-SY5Y cell homeostasis has been extensively studied [62]. Several abundant ncRNAs play an active role in such sponging mechanisms. For example, MALAT1 promotes oxygen–glucose deprivation/reoxygenation (OGD/R)-induced neuronal injury through miRNA competition [59,60]. Similarly, NORAD promotes cell survival by sponging miRNAs [23]. In such cases, even a small change in the accessibility of miRNAs alters the ratio of free to bound cellular miRNAs, leading to a revised cell state by the available miRNA pool [63].
A special case concerns the SNHG genes (Figure 4) that were originally studied in cancer and are characterized by a rich repertoire of miRNA binding sites (MBS). As a result, several of the SNHGs may collaboratively titrate out miRNAs and thus release their bound targets, a mechanism known as ceRNA [25,64]. Of these, SNHG12 was significantly downregulated in the presence of Sin1 in undifferentiated SY-SY5Y cells. In these cells, SNHG12 participates in the unfolded protein response (UPR), apoptosis, inflammation, and oxidative stress [65]. It has been proposed that several of the SNHG members carry anti-inflammatory roles under oxygen–glucose deprivation/reoxygenation (OGD/R) conditions by targeting specific miRNAs [33,66]. As a general theme, lncRNAs, such as SHNGs, act via the sponging of specific miRNAs (e.g., [67,68]).
Among the abundant neuronal lncRNAs that were altered by Sin1 treatment, MALAT-1, MIAT, NEAT1, and MIR137HG (Table 1) contain miR-106a MBSs [69]. This also applies to members of the SNHG lncRNAs (Supplementary Table S4). Specifically, we identified 17 members of the SNHGs that were expressed in SH-SY5Y (Figure 4). Note that MBSs for miR-106a were found in nine of them (SNHG1, SNHG3, SNHG5-7, SNHG10-12, and SNHG14, Supplementary Table S4). Many of these targets were confirmed by miR-CLIP with miR-106a as bait [70]. An inverse link was established between the level of available miR-106a and Treg anti-inflammatory response (e.g., allergy, inflammatory bowel disease (IBD), and LPS-induced macrophages) [71]. In SH-SY5Y, we anticipate that the collaborative sponging of miR-106a by abundant lncRNAs, including SNHGs, underlies the beneficial impact of ladostigil.
Interestingly, among the ncRNA-induced genes following Sin1 is MIR137HG (Table 1), a precursor of hsa-miR-137 and hsa-miR-2682, two highly expressed miRNAs of the CNS and especially the brain cortex. In genome-wide association studies, the MIR137HG gene has been linked to schizophrenia and other psychiatric disorders [72,73]. The function of has-miR-137 and has-miR-2682 in cell proliferation was reported in the context of cancer, where SNHG3 acts as a direct sponge of miR-2682 [74]. Thus, it is likely that the suppression of SNHG3 (Figure 4) and the induction of MIR137HG (1.5 fold, p-value 4.5 × 10−23, Table 1) together increased the accessibility of miR-137 and miR-2682.
An attractive hypothesis for the involvement of abundant ncRNAs under long-lasting stress conditions concerns regulation by the antioxidation response. We recently showed that oxidative stress in SH-SY5Y undifferentiated cells led to an alteration in ER stress and the redox state [12]. Nrf2 (gene name: NFE2L2), a master transcription factor (TF) of redox homeostasis, is upregulated by Sin1, whereas ladostigil inhibits NFE2L2 expression [12]. Furthermore, Nrf2 has been identified as a regulatory hub for lncRNAs [75] and miRNAs [76,77]. The capacity of ladostigil to improve cell viability (Figure 1) could reflect a cross-talk of MALAT1 and Nrf2 in cytotoxicity. We identified MALAT1, the highest expressed lncRNA in undifferentiated SH-SY5Y cells, to be further upregulated by Sin1 (1.24-fold). MALAT1 has been shown to inhibit apoptosis via Nrf2 activation. This is mediated via the suppression of lipid peroxidation and damaged DNA [60].
In this work, we focused on undifferentiated SH-SY5Y cells exposed to abrupt (by H2O2) and long-term oxidative stress (by Sin1). However, exposing SH-SY5Y cells to various combinations of growth factors and morphogens (e.g., RA) drives their differentiation toward mature dopaminergic-related neurons [78,79]. After exposing SH-SY5Y to RA for several days, profound changes in mitochondrial metabolism were observed [80]. While differentiated SH-SY5Y cells maintained the same mitochondrial number, they were distinguished by elevated membrane potentials. Overall, the differentiated cells were more resistant to cytotoxicity and mitochondrial dysfunction. In accordance with these results, we show that RA-differentiated SH-SY5Y coped with Sin1 where no induction of abundant ncRNAs was recorded (Figure 6). Importantly, alternative differentiation protocols that were presented (e.g., the two-step RA and BDNF combination) drove the cells to become more vulnerable to energetic stress, despite the increased ATP production [28]. We anticipate that inspecting specific changes in the lncRNA profile is useful for the study of functional neurons in health and disease [51].
In summary, we present evidence for reduced cytotoxicity in undifferentiated SH-SY5Y cells that were treated by ladostigil. We found that ladostigil shifted the IC50 for cell viability of Sin1 from ~250 μM to ~350 μM, presumably by its antioxidative activity. There were numerous mechanisms by which lncRNAs responded to oxidative stress, including the regulation of transcription by enhancers, competing with miRNAs, modulating transcription factors, and molecular decoys. Among the hundreds of lncRNAs identified in the cells (Supplementary Table S2), tens of abundant ones were moderately upregulated. We show that RA-differentiated cells were quite robust to Sin1 and ladostigil. We anticipate that the oxidative stress-driven shifts in the lncRNA expression levels participated in the cells’ responses to external stress, which distinguished their differentiated state. We propose that even a modest dysregulation in the lncRNA levels may affect cell homeostasis. We propose lncRNAs as cellular indicators for assessing the effects of novel drugs (e.g., ladostigil) under chronic stress conditions that characterize the aging brain and neurodegenerative diseases.

4. Materials and Methods

4.1. Materials

All reagents were purchased from Sigma Aldrich (Burlington, MA, USA) unless otherwise stated. Ladostigil, (6-(N-ethyl, N-methyl) carbamyloxy)-N propargyl-1(R)-aminoindan tartrate, was a gift from Spero BioPharma (Cambridge, MA, USA). Sin1 (3-(4-5 Morpholinyl) sydnone imine hydrochloride) was obtained from Sigma-Aldrich (Cat-M5793). Media products MEM and F12 (ratio 1:1), heat-inactivated fetal calf serum (10% FCS, Cat. C8056, Sigma-Aldrich), and L-Alanyl-L Glutamine were obtained from Biological Industries IMBH (Beit-Haemek, Israel). All tissue culture materials were purchased from IMBH, (Beit-Haemek, Israel). H2O2 (30%, Cat. H1009, Sigma-Aldrich, Rehovot, Israel) was diluted in cold water and the solutions were used within 2 h. Cell viability was determined using MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide; Cat M5655, Sigma-Aldrich, Rehovot, Israel).

4.2. SH-SY5Y Cell Culture

Human neuroblastoma-derived SH-SY5Y cells derived from the original SK-N-SH clone were obtained from the ATCC (American Type Culture Collection, Manassas, VA, USA) [78]. Cells were cultured in Minimum Essential Media (MEM and F12 ratio 1:1, 4.5 g/l glucose) with 10% fetal calf serum (FCS) and 1:10 L-Alanyl-L-Glutamine. Cells were incubated at 37 °C under a humidified atmosphere of 5% CO2. Ladostigil was added to the culture medium 2 h prior to the activation of oxidative stress, and analysis was performed. Cells were tested 24 h later (i.e., 26 h after the addition of ladostigil), unless otherwise stated.
SH-SY5Y cells were differentiated by retinoic acid (RA) according to the published protocol [81]. Briefly, 1.5 × 104 cells/cm2 were seeded in medium supplemented with 10% fetal bovine serum (FBS, Cat. F7524, Sigma-Aldrich). After 24 h, the cells were introduced to a fresh medium with only 1% FBS, supplemented with 10 µM of RA (all-trans-RA, Cat. R625, Sigma-Aldrich), and maintained for an additional 7 days. On day 4, the same medium was replaced with fresh medium. On day 7, the cells were harvested and used for the preparation of the cDNA for the RT-PCR experiments.

4.3. Cell Viability Assay

SH-SY5Y cells were cultured at a density of 2 × 104 cells per well in 96-well plates with 200 µL of medium. An MTT assay was conducted using colorimetric measurement for cell viability, with viable cells producing a dark blue formazan product [82]. The cells were treated in the absence or presence of ladostigil (at 5.4 mM and 54 mM). The drug is soluble in water. The MTT solution in phosphate-buffered saline (PBS, pH 7.2) was prepared as a working stock of 5 mg/mL. After 24 h, the culture medium was supplemented with 10 µL of concentrated MTT per well. Absorption was determined in an ELISA reader at λ = 535 nm, using absorption at λ = 635 nm as a baseline. Cell viability was expressed as the percentage of untreated cells. Each experimental condition was repeated 8 times. The results of the survival assays are presented as the mean ± SEM. When appropriate, p-values < 0.05 were calculated and considered statistically significant. We report on genes that met the FDR-adjusted p-value of ≤0.05 (q-value).

4.4. Flow Cytometry

SH-SY5Y cells were cultured in 6-well or 12-well plates to reach a 70–80% confluence level (37 °C and 5% CO2) before analysis by fluorescence-activated single-cell sorting (FACS). Propidium iodide (PI, Cat. R37108, Thermo-Fisher (Waltham, MA, USA) was used as a marker for dead cells. PI and Annexin V-FITC were purchased from MBL (MEBCYTO-Apoptosis Kit; Woburn, MA, USA). Early apoptosis was detected using conjugated Annexin V [83]. The fluorescence of Annexin-V indicates the presence of phosphatidylserine (PS) at the outer leaflet of the plasma membrane. FACS analysis was conducted with 50,000 cells and the fraction of dead cells stained with PI was determined. Cells were stained with Annexin V according to the manufacturer’s MBL protocol. Cells in the early apoptosis phase are defined as Annexin V-positive and PI-negative. Discriminative staining with PI and Annexin-V by FACS was discussed in [32,84].

4.5. RNA Sequencing

SH-SY5Y cells were pre-incubated for 2 h with ladostigil at a concentration of 5.4 μM. The cells were exposed to Sin1 (t = 0) and harvested at 10 h or 24 h following the addition of ladostigil. The total RNA was extracted using an RNeasy Plus Universal Mini Kit (QIAGEN, GmbH, Hilden, Germany) according to the manufacturer’s protocol. Briefly, a mini spin column was used and centrifuged for 15 s at ≥8000× g and room temperature. Samples with an RNA integrity number (RIN) of >9.0, as measured by an Agilent 2100 Bioanalyzer, were considered for further analysis. Total RNA samples (1 μg RNA) were enriched for mRNAs by the pull-down of poly(A) RNA and libraries were prepared using a KAPA Stranded mRNA-Seq Kit (QIAGEN, GmbH, Hilden, Germany) according to the manufacturer’s protocol. As the RNA-seq experiment was performed using a protocol for the removal of ribosomal RNAs, restricting molecule length to a minimum of >200 nucleotides, a sequencing depth to many of the lowly expressing transcripts could be achieved. RNA-seq libraries were sequenced using Illumina NextSeq 500 to generate 85 bp single-end reads. RNA-seq data files were deposited in ArrayExpress [85] under accession number E-MTAB-10450.

4.6. Reverse Transcription Polymerase–Chain Reaction (RT-PCR)

Total RNA was extracted with Trizol (Thermo-Fisher, Waltham, MA, USA)), and RT was performed using a Ready-To-Go first-strand synthesis kit (Amersham Pharmacia Biotech, Cat. 27925901, Uppsala, Sweden) according to the manufacturer’s instructions. RNA was reverse-transcribed into cDNA (1 μg) and used in the PCR reaction. The primers used for the PCR reactions are listed in Table 3. The PCR conditions consisted of denaturation at 95 °C for 2 min and 35 cycles (10 s at 95 °C, 15 s at 60 °C, and 15 min at 72 °C for extension). The amplicon from the β-actin was used as the internal control for normalization. The PCR products were separated on 1.5–2% agarose gel and stained with ethidium bromide, followed by densitometry measurement (using ImageJ).

4.7. Bioinformatic Analysis and Statistics

All next-generation sequencing (NGS) data underwent quality control using FastQC [86] and were processed using Trimmomatic [87]. All genomic loci were annotated using GENCODE version 37 and aligned to GRCh38 using a STAR aligner [88] and default parameters. All statistical tests were performed using R-base functions. Figures were generated using the R package ggplot2. The experiments contained a minimum of three biological replicates. The trimmed mean of M-values (TMM) was used for the normalization of the RNA read counts. A TMM value of > 4 was used as the lower cutoff value for the gene expression levels. For quantifying relative expression, the sum of TMM values from all lncRNA (> 4) was considered as 100% and only ncRNAs representing > 0.5% of that total TMM are reported. Differential expression analysis was performed using edgeR [89]. Unless otherwise stated, genes with a false discovery rate (FDR)-adjusted p-value of ≤ 0.05 and an absolute log fold-change (FC) above 0.3 were considered as significantly differentially expressed. Data on the miRNA–lncRNA–protein interactome were downloaded from the experimentally verified RAID v2.0 database [69].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ncrna8060072/s1, Table S1: RNA-seq TMM results of triplicate tests for SH-SY5Y cells (NT, lado, NT/lado and Sin1, Sin1/ladostigil. Table S2: RNA-seq-expressed (DE) genes; Sin1 treated vs. Sin1/Ladostigil treatment (24 h). Table S3: ncRNA collection of 780 DE significant genes for Sin1 vs. NT (TMM > 4; Source of Figure 2 and Figure 3). Table S4: SNHG representatives (Source of Figure 4), and their predicted and validated miRNA binding sites. Table S5: Conserved MBSs on ncRNA GAS5 and MTND4P12. Figure S1: Flow cytometry analysis of cells under oxidative stress in response to treatment with ladostigil.

Author Contributions

K.Z. led the research for all experiments, data collection, analysis and visualization. E.G. and T.E. carried out the cellular and drug calibration experiments. M.L. contributed in analysis, conceptualization, mentoring, visualization and writing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Israel Science Foundation (ISF), grant 2753/20 (M.L) and the ZC4H2 Associated Rare Disorders (ZARD) Pilot Grant (fellowship to K.Z.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data used in this study are available as processed data in the Supplementary Tables. The original data from RNA-seq are publicly available from ArrayExpress under accession number E-MTAB-10450.

Acknowledgments

We would like to deeply thank Marta Weinstock for her constant support and supply of ladostigil. We thank Elyad Lezmi for his contribution to the data processing and initial analysis, Naomi van Wijk for reading the manuscript, and the rest of Linial’s lab members for their useful comments and fruitful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ADAlzheimer’s disease
CNScentral nervous system
DEdifferentially expressed
ERendoplasmic reticulum
FACSfluorescence-activated cell sorting
FCfold change
FCSfetal calf serum
FDRfalse discovery rate
H2O2hydrogen peroxide
MBSsmiRNA binding sites
MEMminimum essential medium
NTnot treated
Nrf2nuclear factor erythroid 2-related factor 2
O2superoxide anion
ONOOperoxynitrite
PDParkinson’s disease
PIpropidium iodide
I/Rischemia and reperfusion
RNSreactive nitrogen species
ROSreactive oxygen species
TFtranscription factor
TMMtrimmed mean of M-values.

References

  1. Head, E.; Liu, J.; Hagen, T.; Muggenburg, B.; Milgram, N.; Ames, B.; Cotman, C. Oxidative damage increases with age in a canine model of human brain aging. J. Neurochem. 2002, 82, 375–381. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Kumar, H.; Lim, H.-W.; More, S.V.; Kim, B.-W.; Koppula, S.; Kim, I.S.; Choi, D.-K. The role of free radicals in the aging brain and Parkinson’s disease: Convergence and parallelism. Int. J. Mol. Sci. 2012, 13, 10478–10504. [Google Scholar] [CrossRef] [Green Version]
  3. Stefanatos, R.; Sanz, A. The role of mitochondrial ROS in the aging brain. FEBS Lett. 2018, 592, 743–758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Cobb, C.A.; Cole, M.P. Oxidative and nitrative stress in neurodegeneration. Neurobiol. Dis. 2015, 84, 4–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Gandhi, S.; Abramov, A.Y. Mechanism of oxidative stress in neurodegeneration. Oxidative Med. Cell. Longev. 2012, 2012, 428010. [Google Scholar] [CrossRef] [Green Version]
  6. Loring, J.; Wen, X.; Lee, J.; Seilhamer, J.; Somogyi, R. A gene expression profile of Alzheimer’s disease. DNA Cell Biol. 2001, 20, 683–695. [Google Scholar] [CrossRef]
  7. Cooper-Knock, J.; Kirby, J.; Ferraiuolo, L.; Heath, P.R.; Rattray, M.; Shaw, P.J. Gene expression profiling in human neurodegenerative disease. Nat. Rev. Neurol. 2012, 8, 518–530. [Google Scholar] [CrossRef]
  8. Weinreb, O.; Amit, T.; Bar-Am, O.; Youdim, M.B. Ladostigil: A novel multimodal neuroprotective drug with cholinesterase and brain-selective monoamine oxidase inhibitory activities for Alzheimer’s disease treatment. Curr. Drug Targets 2012, 13, 483–494. [Google Scholar] [CrossRef]
  9. Panarsky, R.; Luques, L.; Weinstock, M. Anti-inflammatory effects of ladostigil and its metabolites in aged rat brain and in microglial cells. J. Neuroimmune Pharmacol. 2012, 7, 488–498. [Google Scholar] [CrossRef] [PubMed]
  10. Shoham, S.; Linial, M.; Weinstock, M. Age-Induced Spatial Memory Deficits in Rats Are Correlated with Specific Brain Region Alterations in Microglial Morphology and Gene Expression. J. Neuroimmune Pharm. 2019, 14, 251–262. [Google Scholar] [CrossRef]
  11. Linial, M.; Stern, A.; Weinstock, M. Effect of ladostigil treatment of aging rats on gene expression in four brain areas associated with regulation of memory. Neuropharmacology 2020, 177, 108229. [Google Scholar] [CrossRef] [PubMed]
  12. Zohar, K.; Lezmi, E.; Eliyahu, T.; Linial, M. Ladostigil Attenuates Induced Oxidative Stress in Human Neuroblast-like SH-SY5Y Cells. Biomedicines 2021, 9, 1251. [Google Scholar] [CrossRef] [PubMed]
  13. Geisler, S.; Coller, J. RNA in unexpected places: Long non-coding RNA functions in diverse cellular contexts. Nat. Rev. Mol. Cell Biol. 2013, 14, 699–712. [Google Scholar] [CrossRef] [Green Version]
  14. Ang, C.E.; Trevino, A.E.; Chang, H.Y. Diverse lncRNA mechanisms in brain development and disease. Curr. Opin. Genet. Dev. 2020, 65, 42–46. [Google Scholar] [CrossRef]
  15. Wang, X.; Shen, C.; Zhu, J.; Shen, G.; Li, Z.; Dong, J. Long Noncoding RNAs in the Regulation of Oxidative Stress. Oxidative Med. Cell. Longev. 2019, 2019, 1318795. [Google Scholar] [CrossRef]
  16. Zhang, D.; Lee, H.; Haspel, J.A.; Jin, Y. Long noncoding RNA FOXD3-AS1 regulates oxidative stress-induced apoptosis via sponging microRNA-150. FASEB J. 2017, 31, 4472–4481. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Liu, W.; Zhang, Q.; Zhang, J.; Pan, W.; Zhao, J.; Xu, Y. Long non-coding RNA MALAT1 contributes to cell apoptosis by sponging miR-124 in Parkinson disease. Cell Biosci. 2017, 7, 1–9. [Google Scholar] [CrossRef] [Green Version]
  18. Zhou, Y.; Xu, Z.; Yu, Y.; Cao, J.; Qiao, Y.; Qiao, H.; Suo, G. Comprehensive analysis of the lncRNA-associated ceRNA network identifies neuroinflammation biomarkers for Alzheimer’s disease. Mol. Omics. 2019, 15, 459–469. [Google Scholar] [CrossRef]
  19. Yao, J.; Wang, X.Q.; Li, Y.J.; Shan, K.; Yang, H.; Wang, Y.N.Z.; Yao, M.D.; Liu, C.; Li, X.M.; Shen, Y. Long non-coding RNA MALAT 1 regulates retinal neurodegeneration through CREB signaling. EMBO Mol. Med. 2016, 8, 346–362. [Google Scholar] [CrossRef]
  20. Saeedi Borujeni, M.J.; Esfandiary, E.; Baradaran, A.; Valiani, A.; Ghanadian, M.; Codoner-Franch, P.; Basirat, R.; Alonso-Iglesias, E.; Mirzaei, H.; Yazdani, A. Molecular aspects of pancreatic beta-cell dysfunction: Oxidative stress, microRNA, and long noncoding RNA. J. Cell Physiol. 2019, 234, 8411–8425. [Google Scholar] [CrossRef]
  21. Meng, K.; Jiao, J.; Zhu, R.-R.; Wang, B.-Y.; Mao, X.-B.; Zhong, Y.-C.; Zhu, Z.-F.; Yu, K.-W.; Ding, Y.; Xu, W.-B. The long noncoding RNA hotair regulates oxidative stress and cardiac myocyte apoptosis during ischemia-reperfusion injury. Oxidative Med. Cell. Longev. 2020, 2020, 1645249. [Google Scholar] [CrossRef] [PubMed]
  22. Chen, H.; Du, G.; Song, X.; Li, L. Non-coding transcripts from enhancers: New insights into enhancer activity and gene expression regulation. Genom. Proteom. Bioinform. 2017, 15, 201–207. [Google Scholar] [CrossRef]
  23. Zhou, X.; Wang, Z.; Xu, B.; Ji, N.; Meng, P.; Gu, L.; Li, Y. Long non-coding RNA NORAD protects against cerebral ischemia/reperfusion injury induced brain damage, cell apoptosis, oxidative stress and inflammation by regulating miR-30a-5p/YWHAG. Bioengineered 2021, 12, 9174–9188. [Google Scholar] [CrossRef] [PubMed]
  24. Moran, V.A.; Perera, R.J.; Khalil, A.M. Emerging functional and mechanistic paradigms of mammalian long non-coding RNAs. Nucleic Acids Res. 2012, 40, 6391–6400. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Yang, M.; Lu, H.; Liu, J.; Wu, S.; Kim, P.; Zhou, X. lncRNAfunc: A knowledgebase of lncRNA function in human cancer. Nucleic Acids Res. 2022, 50, D1295–D1306. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, Y.; He, W.; Ibrahim, S.A.; He, Q.; Jin, J. Circular RNAs: Novel Players in the Oxidative Stress-Mediated Pathologies, Biomarkers, and Therapeutic Targets. Oxidative Med. Cell. Longev. 2021, 2021, 6634601. [Google Scholar] [CrossRef]
  27. Kim, C.; Kang, D.; Lee, E.K.; Lee, J.S. Long noncoding RNAs and RNA-binding proteins in oxidative stress, cellular senescence, and age-related diseases. Oxidative Med. Cell. Longev. 2017, 2017, 2062384. [Google Scholar] [CrossRef] [Green Version]
  28. Forster, J.; Köglsberger, S.; Trefois, C.; Boyd, O.; Baumuratov, A.; Buck, L.; Balling, R.; Antony, P. Characterization of differentiated SH-SY5Y as neuronal screening model reveals increased oxidative vulnerability. J. Biomol. Screen. 2016, 21, 496–509. [Google Scholar] [CrossRef] [Green Version]
  29. Singh, R.J.; Hogg, N.; Joseph, J.; Konorev, E.; Kalyanaraman, B. The peroxynitrite generator, SIN-1, becomes a nitric oxide donor in the presence of electron acceptors. Arch. Biochem. Biophys. 1999, 361, 331–339. [Google Scholar] [CrossRef]
  30. Mathy-Hartert, M.; Mouithys-Mickalad, A.; Kohnen, S.; Deby-Dupont, G.; Lamy, M.; Hans, P. Effects of propofol on endothelial cells subjected to a peroxynitrite donor (SIN-1). Anaesthesia 2000, 55, 1066–1071. [Google Scholar] [CrossRef]
  31. Cuddy, L.K.; Gordon, A.C.; Black, S.A.; Jaworski, E.; Ferguson, S.S.; Rylett, R.J. Peroxynitrite donor SIN-1 alters high-affinity choline transporter activity by modifying its intracellular trafficking. J. Neurosci. 2012, 32, 5573–5584. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Rieger, A.M.; Nelson, K.L.; Konowalchuk, J.D.; Barreda, D.R. Modified annexin V/propidium iodide apoptosis assay for accurate assessment of cell death. J. Vis. Exp. 2011, 50, e2597. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Tripathi, S.; Shree, B.; Mohapatra, S.; Swati; Basu, A.; Sharma, V. The Expanding Regulatory Mechanisms and Cellular Functions of Long Non-coding RNAs (lncRNAs) in Neuroinflammation. Mol. Neurobiol. 2021, 58, 2916–2939. [Google Scholar] [CrossRef]
  34. Giannakakis, A.; Zhang, J.; Jenjaroenpun, P.; Nama, S.; Zainolabidin, N.; Aau, M.Y.; Yarmishyn, A.A.; Vaz, C.; Ivshina, A.V.; Grinchuk, O.V. Contrasting expression patterns of coding and noncoding parts of the human genome upon oxidative stress. Sci. Rep. 2015, 5, 9737. [Google Scholar] [CrossRef] [Green Version]
  35. Zimta, A.A.; Tigu, A.B.; Braicu, C.; Stefan, C.; Ionescu, C.; Berindan-Neagoe, I. An Emerging Class of Long Non-coding RNA With Oncogenic Role Arises From the snoRNA Host Genes. Front Oncol. 2020, 10, 389. [Google Scholar] [CrossRef]
  36. Zhang, S.; Zhu, T.; Li, Q.; Sun, G.; Sun, X. Long Non-Coding RNA-Mediated Competing Endogenous RNA Networks in Ischemic Stroke: Molecular Mechanisms, Therapeutic Implications, and Challenges. Front Pharm. 2021, 12, 765075. [Google Scholar] [CrossRef]
  37. Nagy, A.; Munkacsy, G.; Gyorffy, B. Pancancer survival analysis of cancer hallmark genes. Sci. Rep. 2021, 11, 6047. [Google Scholar] [CrossRef]
  38. Mahmoudi, E.; Cairns, M.J. MiR-137: An important player in neural development and neoplastic transformation. Mol. Psychiatry 2017, 22, 44–55. [Google Scholar] [CrossRef] [Green Version]
  39. Beneventi, G.; Munita, R.; Cao Thi Ngoc, P.; Madej, M.; Ciesla, M.; Muthukumar, S.; Krogh, N.; Nielsen, H.; Swaminathan, V.; Bellodi, C. The small Cajal body-specific RNA 15 (SCARNA15) directs p53 and redox homeostasis via selective splicing in cancer cells. NAR Cancer 2021, 3, zcab026. [Google Scholar] [CrossRef]
  40. Fishilevich, S.; Nudel, R.; Rappaport, N.; Hadar, R.; Plaschkes, I.; Iny Stein, T.; Rosen, N.; Kohn, A.; Twik, M.; Safran, M.; et al. GeneHancer: Genome-wide integration of enhancers and target genes in GeneCards. Database 2017, 2017. [Google Scholar] [CrossRef]
  41. Lang, B.; Armaos, A.; Tartaglia, G.G. RNAct: Protein-RNA interaction predictions for model organisms with supporting experimental data. Nucleic Acids Res. 2019, 47, D601–D606. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Jeggari, A.; Marks, D.S.; Larsson, E. miRcode: A map of putative microRNA target sites in the long non-coding transcriptome. Bioinformatics 2012, 28, 2062–2063. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Korecka, J.A.; van Kesteren, R.E.; Blaas, E.; Spitzer, S.O.; Kamstra, J.H.; Smit, A.B.; Swaab, D.F.; Verhaagen, J.; Bossers, K. Phenotypic characterization of retinoic acid differentiated SH-SY5Y cells by transcriptional profiling. PLoS ONE 2013, 8, e63862. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Weinberg, M.S.; Wood, M.J. Short non-coding RNA biology and neurodegenerative disorders: Novel disease targets and therapeutics. Hum. Mol. Genet. 2009, 18, R27–R39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Riva, P.; Ratti, A.; Venturin, M. The long non-coding RNAs in neurodegenerative diseases: Novel mechanisms of pathogenesis. Curr. Alzheimer Res. 2016, 13, 1219–1231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Wu, Y.Y.; Kuo, H.C. Functional roles and networks of non-coding RNAs in the pathogenesis of neurodegenerative diseases. J. Biomed. Sci. 2020, 27, 49. [Google Scholar] [CrossRef] [Green Version]
  47. Lyu, Y.; Bai, L.; Qin, C. Long noncoding RNAs in neurodevelopment and Parkinson’s disease. Anim. Model Exp. Med. 2019, 2, 239–251. [Google Scholar] [CrossRef] [Green Version]
  48. Fan, Y.; Li, J.; Yang, Q.; Gong, C.; Gao, H.; Mao, Z.; Yuan, X.; Zhu, S.; Xue, Z. Dysregulated Long Non-coding RNAs in Parkinson’s Disease Contribute to the Apoptosis of Human Neuroblastoma Cells. Front Neurosci. 2019, 13, 1320. [Google Scholar] [CrossRef] [Green Version]
  49. Jathar, S.; Kumar, V.; Srivastava, J.; Tripathi, V. Technological Developments in lncRNA Biology. In Long Non Coding RNA Biology; Rao, M.R.S., Ed.; Springer: Singapore, 2017; pp. 283–323. [Google Scholar]
  50. Isin, M.; Ozgur, E.; Cetin, G.; Erten, N.; Aktan, M.; Gezer, U.; Dalay, N. Investigation of circulating lncRNAs in B-cell neoplasms. Clin. Chim. Acta 2014, 431, 255–259. [Google Scholar] [CrossRef]
  51. Taft, R.J.; Pang, K.C.; Mercer, T.R.; Dinger, M.; Mattick, J.S. Non-coding RNAs: Regulators of disease. J. Pathol. 2010, 220, 126–139. [Google Scholar] [CrossRef]
  52. Wang, W.-T.; Ye, H.; Wei, P.-P.; Han, B.-W.; He, B.; Chen, Z.-H.; Chen, Y.-Q. LncRNAs H19 and HULC, activated by oxidative stress, promote cell migration and invasion in cholangiocarcinoma through a ceRNA manner. J. Hematol. Oncol. 2016, 9, 117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Zeni, P.F.; Mraz, M. LncRNAs in adaptive immunity: Role in physiological and pathological conditions. RNA Biol. 2021, 18, 619–632. [Google Scholar] [CrossRef]
  54. Yan, W.; Chen, Z.-Y.; Chen, J.-Q.; Chen, H.-M. LncRNA NEAT1 promotes autophagy in MPTP-induced Parkinson’s disease through stabilizing PINK1 protein. Biochem. Biophys. Res. Commun. 2018, 496, 1019–1024. [Google Scholar] [CrossRef] [PubMed]
  55. Butler, A.A.; Johnston, D.R.; Kaur, S.; Lubin, F.D. Long noncoding RNA NEAT1 mediates neuronal histone methylation and age-related memory impairment. Sci. Signal. 2019, 12, eaaw9277. [Google Scholar] [CrossRef]
  56. Kang, Y.; Zhu, X.; Xu, Y.; Tang, Q.; Huang, Z.; Zhao, Z.; Lu, J.; Song, G.; Xu, H.; Deng, C. Energy stress-induced lncRNA HAND2-AS1 represses HIF1α-mediated energy metabolism and inhibits osteosarcoma progression. Am. J. Cancer Res. 2018, 8, 526. [Google Scholar]
  57. Wu, Q.; Yi, X. Down-regulation of Long Noncoding RNA MALAT1 Protects Hippocampal Neurons Against Excessive Autophagy and Apoptosis via the PI3K/Akt Signaling Pathway in Rats with Epilepsy. J. Mol. Neurosci. 2018, 65, 234–245. [Google Scholar] [CrossRef]
  58. Zeng, R.; Zhang, R.; Song, X.; Ni, L.; Lai, Z.; Liu, C.; Ye, W. The long non-coding RNA MALAT1 activates Nrf2 signaling to protect human umbilical vein endothelial cells from hydrogen peroxide. Biochem. Biophys. Res. Commun. 2018, 495, 2532–2538. [Google Scholar] [CrossRef]
  59. Li, L.; Xu, Y.; Zhao, M.; Gao, Z. Neuro-protective roles of long non-coding RNA MALAT1 in Alzheimer’s disease with the involvement of the microRNA-30b/CNR1 network and the following PI3K/AKT activation. Exp. Mol. Pathol. 2020, 117, 104545. [Google Scholar] [CrossRef]
  60. Jia, Y.; Yi, L.; Li, Q.; Liu, T.; Yang, S. LncRNA MALAT1 aggravates oxygen-glucose deprivation/reoxygenation-induced neuronal endoplasmic reticulum stress and apoptosis via the miR-195a-5p/HMGA1 axis. Biol. Res. 2021, 54, 8. [Google Scholar] [CrossRef]
  61. Li, J.; Zhu, L.; Guan, F.; Yan, Z.; Liu, D.; Han, W.; Chen, T. Relationship between schizophrenia and changes in the expression of the long non-coding RNAs Meg3, Miat, Neat1 and Neat2. J. Psychiatr. Res. 2018, 106, 22–30. [Google Scholar] [CrossRef]
  62. Pandey, A.; Sarkar, S.; Yadav, S.K.; Yadav, S.S.; Srikrishna, S.; Siddiqui, M.H.; Parmar, D.; Yadav, S. Studies on Regulation of Global Protein Profile and Cellular Bioenergetics of Differentiating SH-SY5Y Cells. Mol. Neurobiol. 2022, 59, 1799–1818. [Google Scholar] [CrossRef]
  63. Mahlab-Aviv, S.; Linial, N.; Linial, M. miRNA Combinatorics and its Role in Cell State Control-A Probabilistic Approach. Front Mol. Biosci. 2021, 8, 772852. [Google Scholar] [CrossRef] [PubMed]
  64. Tamang, S.; Acharya, V.; Roy, D.; Sharma, R.; Aryaa, A.; Sharma, U.; Khandelwal, A.; Prakash, H.; Vasquez, K.M.; Jain, A. SNHG12: An LncRNA as a Potential Therapeutic Target and Biomarker for Human Cancer. Front Oncol. 2019, 9, 901. [Google Scholar] [CrossRef] [PubMed]
  65. Yan, L.; Li, L.; Lei, J. Long noncoding RNA small nucleolar RNA host gene 12/microRNA-138-5p/nuclear factor I/B regulates neuronal apoptosis, inflammatory response, and oxidative stress in Parkinson’s disease. Bioengineered 2021, 12, 12867–12879. [Google Scholar] [CrossRef] [PubMed]
  66. Lan, T.; Ma, W.; Hong, Z.; Wu, L.; Chen, X.; Yuan, Y. Long non-coding RNA small nucleolar RNA host gene 12 (SNHG12) promotes tumorigenesis and metastasis by targeting miR-199a/b-5p in hepatocellular carcinoma. J. Exp. Clin. Cancer Res. 2017, 36, 11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Gu, R.; Liu, R.; Wang, L.; Tang, M.; Li, S.-R.; Hu, X. LncRNA RPPH1 attenuates Aβ25-35-induced endoplasmic reticulum stress and apoptosis in SH-SY5Y cells via miR-326/PKM2. Int. J. Neurosci. 2021, 131, 425–432. [Google Scholar] [CrossRef]
  68. Zhao, J.; Geng, L.; Chen, Y.; Wu, C. SNHG1 promotes MPP(+)-induced cytotoxicity by regulating PTEN/AKT/mTOR signaling pathway in SH-SY5Y cells via sponging miR-153-3p. Biol. Res. 2020, 53, 1. [Google Scholar] [CrossRef]
  69. Zhou, Y.K.; Shen, Z.A.; Yu, H.; Luo, T.; Gao, Y.; Du, P.F. Predicting lncRNA-Protein Interactions With miRNAs as Mediators in a Heterogeneous Network Model. Front Genet 2019, 10, 1341. [Google Scholar] [CrossRef]
  70. Imig, J.; Brunschweiger, A.; Brummer, A.; Guennewig, B.; Mittal, N.; Kishore, S.; Tsikrika, P.; Gerber, A.P.; Zavolan, M.; Hall, J. miR-CLIP capture of a miRNA targetome uncovers a lincRNA H19-miR-106a interaction. Nat. Chem. Biol. 2015, 11, 107–114. [Google Scholar] [CrossRef]
  71. Sanctuary, M.R.; Huang, R.H.; Jones, A.A.; Luck, M.E.; Aherne, C.M.; Jedlicka, P.; de Zoeten, E.F.; Collins, C.B. miR-106a deficiency attenuates inflammation in murine IBD models. Mucosal. Immunol. 2019, 12, 200–211. [Google Scholar] [CrossRef] [Green Version]
  72. Warburton, A.; Breen, G.; Rujescu, D.; Bubb, V.J.; Quinn, J.P. Characterization of a REST-regulated internal promoter in the schizophrenia genome-wide associated gene MIR137. Schizophr. Bull. 2015, 41, 698–707. [Google Scholar] [CrossRef] [PubMed]
  73. Fries, G.R.; Carvalho, A.F.; Quevedo, J. The miRNome of bipolar disorder. J. Affect. Disord. 2018, 233, 110–116. [Google Scholar] [CrossRef] [PubMed]
  74. Lu, N.; Yin, Y.; Yao, Y.; Zhang, P. SNHG3/miR-2682-5p/HOXB8 promotes cell proliferation and migration in oral squamous cell carcinoma. Oral. Dis. 2021, 27, 1161–1170. [Google Scholar] [CrossRef] [PubMed]
  75. Bhattacharjee, S.; Li, J.; Dashwood, R.H. Emerging crosstalk between long non-coding RNAs and Nrf2 signaling. Cancer Lett. 2020, 490, 154–164. [Google Scholar] [CrossRef]
  76. Padmavathi, G.; Ramkumar, K.M. MicroRNA mediated regulation of the major redox homeostasis switch, Nrf2, and its impact on oxidative stress-induced ischemic/reperfusion injury. Arch. Biochem. Biophys. 2021, 698, 108725. [Google Scholar] [CrossRef]
  77. Jayasuriya, R.; Ramkumar, K.M. Role of long non-coding RNAs on the regulation of Nrf2 in chronic diseases. Life Sci. 2021, 270, 119025. [Google Scholar] [CrossRef]
  78. Shipley, M.M.; Mangold, C.A.; Szpara, M.L. Differentiation of the SH-SY5Y Human Neuroblastoma Cell Line. J. Vis. Exp. 2016, 108, 53193. [Google Scholar] [CrossRef]
  79. Cheung, Y.-T.; Lau, W.K.-W.; Yu, M.-S.; Lai, C.S.-W.; Yeung, S.-C.; So, K.-F.; Chang, R.C.-C. Effects of all-trans-retinoic acid on human SH-SY5Y neuroblastoma as in vitro model in neurotoxicity research. Neurotoxicology 2009, 30, 127–135. [Google Scholar] [CrossRef]
  80. Schneider, L.; Giordano, S.; Zelickson, B.R.; Johnson, M.S.; Benavides, G.A.; Ouyang, X.; Fineberg, N.; Darley-Usmar, V.M.; Zhang, J. Differentiation of SH-SY5Y cells to a neuronal phenotype changes cellular bioenergetics and the response to oxidative stress. Free Radic. Biol. Med. 2011, 51, 2007–2017. [Google Scholar] [CrossRef] [Green Version]
  81. Lopes, F.M.; da Motta, L.L.; De Bastiani, M.A.; Pfaffenseller, B.; Aguiar, B.W.; de Souza, L.F.; Zanatta, G.; Vargas, D.M.; Schonhofen, P.; Londero, G.F.; et al. RA Differentiation Enhances Dopaminergic Features, Changes Redox Parameters, and Increases Dopamine Transporter Dependency in 6-Hydroxydopamine-Induced Neurotoxicity in SH-SY5Y Cells. Neurotox Res. 2017, 31, 545–559. [Google Scholar] [CrossRef] [Green Version]
  82. Meerloo, J.V.; Kaspers, G.J.; Cloos, J. Cell sensitivity assays: The MTT assay. In Cancer Cell Culture; Springer: Berlin/Heidelberg, Germany, 2011; pp. 237–245. [Google Scholar]
  83. Crowley, L.C.; Marfell, B.J.; Scott, A.P.; Waterhouse, N.J. Quantitation of Apoptosis and Necrosis by Annexin V Binding, Propidium Iodide Uptake, and Flow Cytometry. Cold Spring Harb. Protoc. 2016, 2016. [Google Scholar] [CrossRef] [PubMed]
  84. Sawai, H.; Domae, N. Discrimination between primary necrosis and apoptosis by necrostatin-1 in Annexin V-positive/propidium iodide-negative cells. Biochem. Biophys. Res. Commun. 2011, 411, 569–573. [Google Scholar] [CrossRef] [PubMed]
  85. Brazma, A.; Parkinson, H.; Sarkans, U.; Shojatalab, M.; Vilo, J.; Abeygunawardena, N.; Holloway, E.; Kapushesky, M.; Kemmeren, P.; Lara, G.G. ArrayExpress—A public repository for microarray gene expression data at the EBI. Nucleic Acids Res. 2003, 31, 68–71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Brown, J.; Pirrung, M.; McCue, L.A. FQC Dashboard: Integrates FastQC results into a web-based, interactive, and extensible FASTQ quality control tool. Bioinformatics 2017, 33, 3137–3139. [Google Scholar] [CrossRef] [Green Version]
  87. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  88. Zhang, P.; Hung, L.-H.; Lloyd, W.; Yeung, K.Y. Hot-starting software containers for STAR aligner. GigaScience 2018, 7, giy092. [Google Scholar] [CrossRef]
  89. Robinson, M.; McCarthy, D.; Chen, Y.; Smyth, G.K. edgeR: Differential expression analysis of digital gene expression data. User’s Guide 2013. Available online: http://rileylab.org/wp-content/uploads/2016/08/edgeRUsersGuide.pdf (accessed on 24 September 2022).
Figure 1. Flow cytometry analysis in cells under oxidative stress conditions in response to treatment with ladostigil. (A) Results from FACS analysis with 50k cells per experiment. The fraction of live and dead cells upon 40 μM H2O2 and ladostigil treatment was tested at 5 h following stimulus. Asterisks indicate the statistical significance of the results. The attenuation in the fraction of PI-positive cells by ladostigil (5.4 μM) was substantial (p-value = 0.02). (B) Cell viability upon increasing the concentration of Sin1 with and without the presence of ladostigil (5.4 μM, with mean and standard error; s.e.). Each data point is based on 8 replicates. (C) Results from 3 repeated experiments with mean and standard deviation (s.d.) from the flow cytometry run of SY-SY5Y cells treated with Sin1 for 5 h. All measured data were from FACS that was repeated in triplicate. The total FACS cells was considered as 100%, and the percentage of them stained with PI is indicated as a fraction. (D) Results from triplicate experiments with the mean and standard deviation (s.d.) from the flow cytometry run for Annexin-V of SY-SY5Y cells treated with Sin1 for 5 h. NT indicates cells non-treated with Sin1. Asterisks indicate the t-test statistical results with p-values of <0.005 (***), <0.01 (**), and <0.05 (*).
Figure 1. Flow cytometry analysis in cells under oxidative stress conditions in response to treatment with ladostigil. (A) Results from FACS analysis with 50k cells per experiment. The fraction of live and dead cells upon 40 μM H2O2 and ladostigil treatment was tested at 5 h following stimulus. Asterisks indicate the statistical significance of the results. The attenuation in the fraction of PI-positive cells by ladostigil (5.4 μM) was substantial (p-value = 0.02). (B) Cell viability upon increasing the concentration of Sin1 with and without the presence of ladostigil (5.4 μM, with mean and standard error; s.e.). Each data point is based on 8 replicates. (C) Results from 3 repeated experiments with mean and standard deviation (s.d.) from the flow cytometry run of SY-SY5Y cells treated with Sin1 for 5 h. All measured data were from FACS that was repeated in triplicate. The total FACS cells was considered as 100%, and the percentage of them stained with PI is indicated as a fraction. (D) Results from triplicate experiments with the mean and standard deviation (s.d.) from the flow cytometry run for Annexin-V of SY-SY5Y cells treated with Sin1 for 5 h. NT indicates cells non-treated with Sin1. Asterisks indicate the t-test statistical results with p-values of <0.005 (***), <0.01 (**), and <0.05 (*).
Ncrna 08 00072 g001
Figure 2. Partition of ncRNA biotypes in SH-SY5Y cells. (A) Partition of the RNA-seq results of untreated SH-SY5Y cells according to the major molecular biotypes. Altogether, there were 19,475 identified transcripts. Blue, coding genes and purple, ncRNAs portioned by biotype classes. Mito, mitochondrial transcript; TEC, to be experimentally confirmed; and PG, pseudogene. Combined related biotypes (e.g., transcribed unprocessed and unitary pseudogenes) are marked as multi. (B) Partition of 780 expressed ncRNA transcripts partitioned by major ncRNA classes. The ncRNAs were based on transcripts with an average TMM of >4. There were 11,588 aligned transcripts that met this threshold. (C) Top expressing ncRNAs (TMM > 100) colored by broad-sense functional annotations (inset). MT, mitochondria; snRNA, small nucleus RNA; and XCI, X-chromosome inactivation.
Figure 2. Partition of ncRNA biotypes in SH-SY5Y cells. (A) Partition of the RNA-seq results of untreated SH-SY5Y cells according to the major molecular biotypes. Altogether, there were 19,475 identified transcripts. Blue, coding genes and purple, ncRNAs portioned by biotype classes. Mito, mitochondrial transcript; TEC, to be experimentally confirmed; and PG, pseudogene. Combined related biotypes (e.g., transcribed unprocessed and unitary pseudogenes) are marked as multi. (B) Partition of 780 expressed ncRNA transcripts partitioned by major ncRNA classes. The ncRNAs were based on transcripts with an average TMM of >4. There were 11,588 aligned transcripts that met this threshold. (C) Top expressing ncRNAs (TMM > 100) colored by broad-sense functional annotations (inset). MT, mitochondria; snRNA, small nucleus RNA; and XCI, X-chromosome inactivation.
Ncrna 08 00072 g002
Figure 3. Differentially expressed (DE) ncRNAs following 24 h of exposure to Sin1 relative to untreated cells. (A) Transcripts associated with upregulation (Up) and downregulation (Down) relative to naïve, untreated cells. The number of genes associated with each expression trend is shown in parentheses. The rest of the genes (9719 genes) were unchanged. See the Materials and Methods for the expression trend definition. (B) Histogram showing the DE ncRNAs ratio for downregulated and upregulated transcripts. The dashed line marks the expression ratio (1.25 and 0.8) vs. naïve cells for Up and Down, respectively. The pie diagram applies to 94 genes that satisfy the threshold.
Figure 3. Differentially expressed (DE) ncRNAs following 24 h of exposure to Sin1 relative to untreated cells. (A) Transcripts associated with upregulation (Up) and downregulation (Down) relative to naïve, untreated cells. The number of genes associated with each expression trend is shown in parentheses. The rest of the genes (9719 genes) were unchanged. See the Materials and Methods for the expression trend definition. (B) Histogram showing the DE ncRNAs ratio for downregulated and upregulated transcripts. The dashed line marks the expression ratio (1.25 and 0.8) vs. naïve cells for Up and Down, respectively. The pie diagram applies to 94 genes that satisfy the threshold.
Ncrna 08 00072 g003
Figure 4. Expression of members of the small nucleolar RNA host genes (SNHGs). (A) Mean expressions (by log10TMM) of the 17 SNHG members identified by DE lncRNA analysis. (B) Fold changes (FCs) following treatment with 50 μM Sin1 are shown for each of the SNHGs. FCs below and above zero show downregulated and upregulation of the indicated genes. Note that the values are logFC. Despite the small effect of the differential expression genes, the filled bar shows statistically significant genes, and the striped bar indicates a gene belonging to the SNHG group, but failed to reach FDR (p-value < 0.05).
Figure 4. Expression of members of the small nucleolar RNA host genes (SNHGs). (A) Mean expressions (by log10TMM) of the 17 SNHG members identified by DE lncRNA analysis. (B) Fold changes (FCs) following treatment with 50 μM Sin1 are shown for each of the SNHGs. FCs below and above zero show downregulated and upregulation of the indicated genes. Note that the values are logFC. Despite the small effect of the differential expression genes, the filled bar shows statistically significant genes, and the striped bar indicates a gene belonging to the SNHG group, but failed to reach FDR (p-value < 0.05).
Ncrna 08 00072 g004
Figure 5. Sample of the lncRNAs by the measured expression levels in naïve untreated cells (NT), cells exposed to Sin1 (24 h) with or without ladostigil. (A) The transcripts shown are those with maximal induction by Sin 1 relative to NT. The average TMM and the standard deviations from biological triplicates are shown for each experimental group. (B) The transcripts shown are those with a maximal effect by ladostigil in the presence of Sin1. The average TMM and the standard deviations from biological triplicates are shown for each experimental group. The statistical significance with the p-value FDR ranging from 0.05–0.005 is depicted by * and <0.005 by **. (C) The genomic segments of AC105339.2 and RPARP-AS1 cover 10 and 12k nucleotides, respectively. AC105339.2 and RPARP-AS1 show the overlap between the lncRNA transcript and the gene in the vicinity. AC105339.2 acts as antisense to FSD2 (704 aa) and RPARP-AS1 is the antisense for C10orf95 (213 aa). Numerous transcripts of RPARP-AS1 are shown above the ENCODE annotations for enhancers (red) and proximal promotors (yellow).
Figure 5. Sample of the lncRNAs by the measured expression levels in naïve untreated cells (NT), cells exposed to Sin1 (24 h) with or without ladostigil. (A) The transcripts shown are those with maximal induction by Sin 1 relative to NT. The average TMM and the standard deviations from biological triplicates are shown for each experimental group. (B) The transcripts shown are those with a maximal effect by ladostigil in the presence of Sin1. The average TMM and the standard deviations from biological triplicates are shown for each experimental group. The statistical significance with the p-value FDR ranging from 0.05–0.005 is depicted by * and <0.005 by **. (C) The genomic segments of AC105339.2 and RPARP-AS1 cover 10 and 12k nucleotides, respectively. AC105339.2 and RPARP-AS1 show the overlap between the lncRNA transcript and the gene in the vicinity. AC105339.2 acts as antisense to FSD2 (704 aa) and RPARP-AS1 is the antisense for C10orf95 (213 aa). Numerous transcripts of RPARP-AS1 are shown above the ENCODE annotations for enhancers (red) and proximal promotors (yellow).
Ncrna 08 00072 g005
Figure 6. Differentially expressed (DE) ncRNAs following 24 h of exposure to Sin1 relative to untreated cells. (A) RT-PCR results on Agarose gel. The RT-PCR product of β-actin was used as an internal control. Each condition is shown for two biological duplicates. The indicated ncRNA results are for the untreated culture (N.T.) and cells in the presence of Sin1, and Sin1 and ladostigil (5.4 μM). (B) Quantification of the results in (A) by a heatmap visualization. Orange, comparison of the baseline of the untreated cells. Blue, changes in cells by Sin1, and Sin1 and ladostigil (5.4 μM).
Figure 6. Differentially expressed (DE) ncRNAs following 24 h of exposure to Sin1 relative to untreated cells. (A) RT-PCR results on Agarose gel. The RT-PCR product of β-actin was used as an internal control. Each condition is shown for two biological duplicates. The indicated ncRNA results are for the untreated culture (N.T.) and cells in the presence of Sin1, and Sin1 and ladostigil (5.4 μM). (B) Quantification of the results in (A) by a heatmap visualization. Orange, comparison of the baseline of the untreated cells. Blue, changes in cells by Sin1, and Sin1 and ladostigil (5.4 μM).
Ncrna 08 00072 g006
Table 1. Differentially expressed lncRNA genes in cells exposed to Sin1 (24 h).
Table 1. Differentially expressed lncRNA genes in cells exposed to Sin1 (24 h).
LncRNALncRNA—Full NameFDRDE (Fold)TMM a
MALAT1Metastasis-associated lung adenocarcinoma transcript 12.5 × 10−41.241963.7
HAND2-AS1HAND2 antisense RNA 16.7 × 10−251.33671.8
GABPB1-AS1GABPB1 antisense RNA 11.7 × 10−181.38280.8
MIATMyocardial infarction associated transcript6.3 × 10−211.35176.9
NEAT1Nuclear paraspeckle assembly transcript 13.5 × 10−51.49154.0
CASC15Cancer susceptibility 151.5 × 10−161.2785.8
AC093849.2Novel transcript2.9 × 10−331.4685.0
CCDC144NL-AS1CCDC144NL antisense RNA 11.6 × 10−121.3267.5
MIR137HGMIR137 host gene4.5 × 10−231.4758.7
GABPB1-IT1GABPB1 intronic transcript4.6 × 10−151.3545.6
LINC02268Long intergenic non-protein coding RNA 22681.2 × 10−161.4544.6
TTN-AS1TTN antisense RNA 11.2 × 10−61.2643.0
LINC01833Long intergenic non-protein coding RNA 18336.8 × 10−141.3242.0
LINC00632Long intergenic non-protein coding RNA 6325.0 × 10−81.3140.2
AL392083.1Novel transcript, sense intronic to ADARB21.8 × 10−331.7837.6
GATA3-AS1GATA3 antisense RNA 12.3 × 10−81.3236.5
AC012354.1Novel transcript2.1 × 10−81.2832.7
AC015813.1Novel transcript4.3 × 10−71.3531.8
AC136621.1Novel transcript (TEC)2.6 × 10−270.5637.7
AL645608.2Novel transcript4.7 × 10−80.7135.9
AC010931.1Novel transcript1.5 × 10−190.6834.9
LINC01128Long intergenic non-protein coding RNA 11281.2 × 10−160.7132.2
EP400P1EP400 pseudogene 14.3 × 10−60.7832.1
a TMM, trimmed mean of M-values.
Table 2. ncRNA upregulated by ladostigil treatment under oxidative stress conditions.
Table 2. ncRNA upregulated by ladostigil treatment under oxidative stress conditions.
ncRNAGH ScoreTotal ScoreElite GH aTSS
(kb)
# of DB Evidence TFsCoding TargetsncRNA TargetsSelected Target
AC010680.10.3250.7*64.01/51.223FKBP7
AC022400.61.9250.7**486.45/54.61013SEC24C
AC105339.22.111.3*24.15/59.875WHAMM
RPARP-AS12.1257.3**546.85/56.6135PPRC1
a GeneHancer (GH) definition of elite (*) gene based on the enhancer score and/or high confidence for the gene targets. High confidence scores are shown in bold. Genes that meet high confidence GH and Total scores are marked (**). Selected gene, a representative gene from the set of coding targets.
Table 3. Forward and reverse oligonucleotide primers used for gel-based RT-PCR assays.
Table 3. Forward and reverse oligonucleotide primers used for gel-based RT-PCR assays.
SymbolTranscriptForward PrimerReverse PrimerAmplicon (nt)
MALAT1NR_144568.1GCTCTGTGGTGTGGGATTGACTCGGGCGAGGCGTATTTAT386
NEAT1NR_028272.1GGGACAACATTGACCAACGCACCACGGTCCATGAAGCATT356
MIATNR_033321.2TCCCATTCCCGGAAGCTAGAGAGGCATGAAATCACCCCCA274
SNHG12NR_146383.1CCTTCTCTCGCTTCGGACTGATCTGCTTAAGTACGCCGGG167
ACTBNM_001101.5ACAGAGCCTCGCCTTTGCCGACATGCCCACCATCACGCCCTGG196
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zohar, K.; Giladi, E.; Eliyahu, T.; Linial, M. Oxidative Stress and Its Modulation by Ladostigil Alter the Expression of Abundant Long Non-Coding RNAs in SH-SY5Y Cells. Non-Coding RNA 2022, 8, 72. https://doi.org/10.3390/ncrna8060072

AMA Style

Zohar K, Giladi E, Eliyahu T, Linial M. Oxidative Stress and Its Modulation by Ladostigil Alter the Expression of Abundant Long Non-Coding RNAs in SH-SY5Y Cells. Non-Coding RNA. 2022; 8(6):72. https://doi.org/10.3390/ncrna8060072

Chicago/Turabian Style

Zohar, Keren, Eliran Giladi, Tsiona Eliyahu, and Michal Linial. 2022. "Oxidative Stress and Its Modulation by Ladostigil Alter the Expression of Abundant Long Non-Coding RNAs in SH-SY5Y Cells" Non-Coding RNA 8, no. 6: 72. https://doi.org/10.3390/ncrna8060072

APA Style

Zohar, K., Giladi, E., Eliyahu, T., & Linial, M. (2022). Oxidative Stress and Its Modulation by Ladostigil Alter the Expression of Abundant Long Non-Coding RNAs in SH-SY5Y Cells. Non-Coding RNA, 8(6), 72. https://doi.org/10.3390/ncrna8060072

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