Molecular Vision 2019; 25:734-744
Received 10 December 2018 | Accepted 11 November 2019 | Published 14 November 2019
The first two authors contributed equally to this work.
1Department of Ophthalmology, Affiliated Hospital of Nantong University, Nantong, Jiangsu, China; 2Department of Ophthalmology, Lixiang Eye Hospital of Soochow University, Suzhou, Jiangsu, China.
Correspondence to: Huaijin Guan, Department of Ophthalmology, Affiliated Hospital of Nantong University, Nantong, Jiangsu, China. firstname.lastname@example.org.
Purpose: Age-related cataract (ARC) is the leading cause of visual impairment and blindness worldwide. The apoptosis of lens epithelial cells (LECs) induced by oxidative damage is a major contributing factor to ARC. Long noncoding RNAs (lncRNAs) play important roles in various biologic processes. We aimed to explore the role of glutathione peroxidase 3 (GPX3)-antisense (AS) in ARCs.
Methods: We extracted total RNAs from transparent and age-matched cataractous human lenses and detected lncRNA expression profiles using high-throughput RNA sequencing. The expression of GPX3-AS and GPX3 was detected by quantitative real-time PCR (qRT-PCR). Apoptotic proteins were detected by western blot and immunofluorescence. We treated SRA01/04 cells with H2O2 to mimic oxidative stress and induce cell apoptosis, which was analyzed by flow cytometry and TdT-mediated dUTP Nick-End Labeling (TUNEL) assay. The cell counting kit-8 (CCK-8) assay was used to detect the viability of SRA01/04 cells. The location of GPX3-AS was determined by fluorescence in situ hybridization (FISH) and cell nuclear and cytoplasmic RNA separation.
Results: The lncRNA GPX3-AS, which is located in the nuclei of LECs, was downregulated in cataractous human lenses compared with control lenses, and proapoptotic proteins were expressed at high levels in the anterior lens capsules of ARC tissues. An in vitro study suggested that GPX3-AS inhibited H2O2-induced SRA01/04 cell apoptosis. As GPX3-AS is transcribed from the AS strand of the GPX3 gene locus, we further revealed its regulatory role in GPX3 expression. GPX3-AS was positively correlated with GPX3 expression. In addition, GPX3-AS inhibited H2O2-induced SRA01/04 cell apoptosis by upregulating GPX3 expression.
Conclusions: In summary, our study revealed that GPX3-AS downregulated the apoptosis of LECs via promoting GPX3 expression, implying a novel therapeutic target for ARCs.
Age-related cataract (ARC) is a common cause of visual impairment and blindness among elderly individuals worldwide . Cataract surgery is currently the most effective therapeutic method for treating ARC . However, surgery is invasive, and it inevitably brings multiple complications, including posterior capsule tearing and vitreous prolapse . Furthermore, many people, especially those in developing countries, become blind from cataracts because of inadequate surgery availability or unaffordable surgical expenses. Therefore, ascertaining the determinant factor in cataract pathogenesis is crucial for reducing the cataract incidence and rate of blindness.
The pathogenesis of ARC is not completely understood. Substantial evidence shows that oxidative stress is a major predisposing factor in ARC [4-7]. H2O2 is a major reactive oxygen species (ROS) present in the aqueous humor and lens . Cataract patients may have deficient defense systems, such as oxidative stress at the onset of the disease . Such stress triggers the apoptosis of lens epithelial cells (LECs), thereby initiating cataract development [10,11]. Therefore, LEC apoptosis induced by oxidative stress appears to be a common cellular basis for the development of noncongenital cataracts .
Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with similar structures to protein-coding mRNAs, which have little or no protein-coding capability . LncRNAs play regulatory roles in diverse biological processes, including cell apoptosis [14,15]. LncRNA dysfunction is involved in multiple human diseases, such as cancer, neurological problems, and cardiovascular diseases. Recently, several studies have confirmed that lncRNAs have a close relationship with ARCs [16-18]. To date, several isoforms of glutathione peroxidase (GPX) proteins have been identified. Of these, only GPX3, which scavenges H2O2 and peroxidized organic molecules to reduce systemic oxidative stress, is secreted . Cortical neurons subjected to O2 deprivation and low glucose (ODLG) display a loss of mitochondrial respiration. However, despite ROS production, neither necrosis nor apoptosis occurs. The absence of cellular death is a consequence of increased antioxidants, such as superoxide dismutase-1 (SOD1) and GPX3 . In addition, GPX3 has been shown to regulate the antioxidative effects of retinoic acid and promote the viability of human muscle stem cells .
The relationship between GPX3-antisense (AS) and ARC is unclear. In the present study, we identified the novel ARC-associated lncRNA GPX3-AS. We aimed to reveal the roles of GPX3-AS in ARC and seek a potential lncRNA-based therapeutic target.
The transparent lens epithelium samples were collected from the patients with vitreoretinal diseases (no other ocular diseases or systemic diseases,5 for RNA sequence and 20 for qRT-PCR verification) and ARC patients (no other ocular diseases or systemic diseases,6 for RNA sequence and 60 for qRT-PCR verification). We divided ARC patients into three groups according to the location of the lens opacity: the age-related cortical cataract (ARCC), age-related nuclear cataract (ARNC), and age-related posterior subcapsular cataract (ARPC) groups, which had 20 patients each. The specific criteria can be referred to our previous studies . All of them had their lenses removed at the Affiliated Hospital of Nantong University (Nantong, China) from January 2017 to June 2018 and the anterior lens capsules that had torn off during surgery were collected immediately.
This study protocol was approved by the Affiliated Hospital of Nantong University. All experiments were performed in accordance with the Declaration of Helsinki. Informed consent was obtained from all patients prior to the study.
The collected samples were stored in liquid nitrogen. Total RNA was extracted using TRIzol reagent (Invitrogen; Thermo Fisher, Waltham, MA). The yield of the RNA was determined using a NanoDrop 2000 spectrophotometer (Thermo Fisher). LncRNA high-throughput sequencing analysis, which included labeling, hybridization, scanning, normalization, and data analysis, was performed by Ribobio Technology (Guangdong, China).
A quality check of the input total RNA was performed by running an aliquot on an Agilent Bioanalyzer (Agilent Technologies, Santa Clara, CA) to confirm its integrity, and the RNA concentration was measured using an ultramicrospectrophotometer. A 3-μg quantity of RNA was selected as the initial amount of each sample from which to construct an lncRNA library, following removal of the rRNA in the sample using Ribo-Zero™ Gold kits (Cat. No. GZG122; Epicenter, Madison, WI). Different index tags were selected according to the New England Biolabs (NEB) Next Ultra database direct operating protocol (NEB Next Ultra Directional RNA Library Prep kit for Illumina; Cat. No. E7760; NEB, Ipswich, MA). The specific steps of library construction were as follows: First, the Ribo-Zero™ Gold kit was used to remove rRNA. Then, fragmentation buffer was added to the reaction system to fragment the RNA into short segments, and the short fragments were used as templates for first-strand cDNA synthesis using random six-base primers (random hexamers). This was followed by second-strand cDNA synthesis using DNA polymerase I and RNaseH. The cDNA fragments were then subjected to 3′-end repair, the addition of a single “A” base, and the ligation of Illumina multiplexing adapters. The products were purified and enriched by PCR to create the final cDNA sequencing library. During PCR, the thermocycling conditions were as follows: step 1, 65 °C for 15 min; step 2, 30 °C for 10 min, 42–50 °C for 45 min, and 95 °C for 5 min.
The raw data were filtered to provide high-quality, clean reads, and a follow-up analysis was conducted. The data processing steps were as follows: i) removal of contaminated reads, ii) removal of low-quality reads, and iii) removal of reads containing >5% bases with nitrogen content .
Only lncRNAs were analyzed in the present study. According to the characteristics of the lncRNA, a series of strict filters were used. The filter conditions were a length of ≥200 bp, an exon number of ≥2, and minimum coverage of ≥3 transcription reads. Known mRNA transcripts, known noncoding RNA transcripts, homologous transcripts, and transcripts with protein-coding capability were removed.
In this study, we used TopHat software version 2.0.12 (Center for Computational Biology, Johns Hopkins, Baltimore, MD) with the default parameters to compare the reference sequence to hg19 .
The cataract and normal samples from each group were subjected to differential expression analysis using DESeq2. Differential expression analysis was performed following standardization and variance estimation. The standards used when screening the differences in mRNA and lncRNA expression were p<0.001, a false discovery rate (FDR) of <0.001, and a |log2 ratio| >2 .
The SRA01/04 (Jennio Biotech Company, Guangzhou, China) human lens epithelial cell (HLEC) line was grown in Dulbecco’s modified Eagle medium (DMEM, Gibco) with 10% fetal bovine serum (FBS, Gibco), penicillin (100 U/ml), and streptomycin (100 mg/ml). The cells were maintained at 37 °C in a humidiﬁed incubator in 5% CO2. At 80–90% confluence, the cells were exposed to H2O2 (200 µM) for 24 h to mimic oxidative stress.
Total RNA was extracted from SRA01/04 cells with TRIzol reagent (Invitrogen), and quantitative real-time PCR (qRT-PCR) was performed with the powerUP SYBR Green Master Mix (Cat. No. A25742, Applied Biosystems, Foster City, CA). Relative gene expression was calculated using the 2-△△CT method, and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) served as an internal control. The AS oligos (ASOs) and overexpressed plasmid were purchased from RiboBio (Guangzhou, China) and GenePharma (Suzhou, China; Table 1), respectively. Cell transfection was performed with Lipofectamine 2000 (Invitrogen).
Nineteen short tandem repeat (STR) loci plus the gender-determining locus, amelogenin, were amplified using the commercially available EX20 Kit from AGCU (Cat. No. 600,010, Suzhou, China). A sample from the cell line was processed using the ABI Prism® 3500 Genetic Analyzer. The data were analyzed using GeneMapper® ID-X v 1.2 software (Applied Biosystems, Suzhou, China). Appropriate positive and negative controls were run and confirmed for each sample submitted.
Cell viability was examined by Cell Counting Kit-8 (CCK-8) assay (Cat. No. CK04–500T, Dojindo Laboratories, Japan). Cells were seeded in 96-well plates at a density of 4 × 104 cells/well. After transfection or treatment, 10 µl of CCK-8 solution was added to each well, and the cells were incubated at 37 °C for 2 h. Finally, the absorbance at 450 nm was calculated.
After H2O2 treatment or transfection, SRA01/04 cells were washed twice with PBS (in 800ml distilled H2O. 8g of NaCl; 0.2g of KCl; 1.44g of Na2HPO4; 0.24g of KH2PO4) and collected for apoptosis assay. Apoptosis was examined by annexin V-fluorescein isothiocyanate (FITC)/propidium iodide (PI) labeling with an apoptosis detection kit (Cat. No. AD10–50T, Dojindo Laboratories). The samples were analyzed within 60 min by flow cytometry (BD Biosciences) after incubation for 15 min at room temperature (RT) in the dark.
Proteins from human lens tissue and SRA01/04 cell lysates were separated by 12% sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) and then transferred to a polyvinylidene fluoride (PVDF) membrane. Primary mouse antibodies against caspase-3 (Abcam, UK), BCL2 associated X(Bax), B-cell lymphoma-2(Bcl-2), and GAPDH (Proteintech), as well as a secondary antibody, horseradish peroxidase (HRP)-conjugated anti-mouse immunoglobulin G (IgG; Sigma), were used. Band intensities were quantified by Quantity One software. The protein levels were normalized to the GAPDH level.
To detect GPX3-AS expression in SRA01/04 cells, Cyanine 3 Dye (Cy3)-labeled probes for GPX3-AS, 18S, and U6 (RiboBio, Guangzhou, China) were used for fluorescence in situ hybridization (FISH). Briefly, SRA01/04 cells were fixed with 4% paraformaldehyde at RT for 15 min, permeabilized in 0.5% Triton X-100 for 10 min on ice, and then treated with prehybridization buffer for 30 min at 37 °C. Next, we hybridized the cells with Cy3-labeled probes in a moist, dark chamber at 37 °C. After 12 h, the cells were stained with 4′,6-diamidino-2-phenylindole (DAPI). Finally, the FISH results were obtained by a confocal microscope (SP8, Leica, Germany).
An In Situ Cell Death Detection Kit (TUNEL fluorescence FITC kit, Cat. No. 11,684,795,910, Roche) was used to detect cell apoptosis. SRA01/04 cells were placed in 24-well culture plates. After treatment, the cells were fixed with 4% paraformaldehyde and permeabilized with 0.1% Triton X-100 on ice for 5 min. Then, 50 μl of TUNEL reaction mixture was added to the samples and incubated at 37 °C for 60 min. The nuclei were stained with DAPI. TUNEL staining was assessed via fluorescence microscopy (Leica, Germany).
Nuclear and cytoplasmic fractions were segregated using a protein and RNA isolation system (PARIS) Kit (Cat. No. AM1556, Life Technologies). Cells were washed twice with PBS, resuspended in precooled cell fractionation buffer (300 μl), and incubated for 10 min on ice. Following centrifugation for 5 min at 500 ×g and 4 °C, the supernatant (cytoplasm) and pellet (nucleus) were segmented. RNA was isolated from both the cytoplasmic and nuclear fractions by a binding and washing procedure. The RNA was prepared for the following reverse transcription reaction, and GPX3-AS was verified using qRT-PCR. U6 RNA and GAPDH mRNA were selected as the nuclear and cytoplasmic control transcripts, respectively.
SRA01/04 cells were fixed with 4% paraformaldehyde for 20 min at RT after treatment and then permeabilized with 0.5% Triton X-100, blocked with bovine serum albumin (BSA), and incubated with primary antibody for GPX3 (Proteintech) at 4 °C overnight. After washing with PBS, the LECs were incubated with the corresponding DyLight 488-conjugated secondary antibodies (EarthOx, San Francisco, CA). Cellular nuclei were counterstained with DAPI. Immunofluorescence was observed with a fluorescence microscope (Leica).
GraphPad Prism 5.0 (GraphPad Software, San Diego, CA) was applied for statistical analysis. The Student t-test was used for comparisons between groups. All presented data were based on biologic triplicates. Data are shown as the mean ± standard error of the mean (SEM), and p<0.05 indicates statistical significance.
To identify ARC-associated lncRNAs, we analyzed RNA sequencing data from six anterior lens capsules from cataract patients and five matched normal capsular membranes. A total of 231 differentially expressed lncRNAs, including 49 downregulated and 182 upregulated lncRNAs, were detected (Figure 1A). GPX3-AS (NONCODE Transcript ID: NONHSAT204874.1; NONCODE Gene ID: NONHSAG091317.1) was decreased in the anterior lens capsules of ARCs compared with its expression in normal tissues, as shown by the sequencing results and verified to be consistent with the sequencing results of the expanded samples by qRT-PCR (Figure 1B). Subgroup analysis performed to identify which ARC subtype was mostly correlated with GPX3-AS suggested that GPX3-AS expression was most downregulated in the ARCC and ARNC groups (Figure 1C). We also measured apoptosis-related protein expression in the anterior lens capsules by immunofluorescence and western blot. The expression levels of active caspase-3 and Bax were upregulated, while Bcl-2 was decreased (Figure 1D-F). These data imply that GPX3-AS may be correlated with cell apoptosis during cataract formation.
We identified the subcellular localization of GPX3-AS. RNA FISH and nuclear separation experiments showed that GPX3-AS was expressed in HLECs (SRA01/04 cells) and mainly localized in the nucleus (Figure 2A,B).
SRA01/04 cells were exposed to H2O2 (200 μM) for 24 h to mimic oxidative stress, following which, GPX3-AS levels decreased (Figure 3A). In addition, the SRA01/04 cell viability was inhibited by H2O2 treatment (Figure 3B). Correspondingly, the TUNEL assay showed that H2O2 induced SRA01/04 cell apoptosis (Figure 3C), and flow cytometry revealed increased apoptosis after H2O2 treatment (Figure 3F,G). Furthermore, western blotting showed that active caspase-3 and Bax were upregulated, while Bcl-2 was decreased (Figure 3D,E). The data suggest that H2O2 treatment induced a decrease in GPX3-AS levels and increased SRA01/04 cell apoptosis, implying a relationship between GPX3-AS and SRA01/04 cell apoptosis.
To verify the role of GPX3-AS in H2O2-induced SRA01/04 cell apoptosis, we used transfection to modulate the GPX3-AS levels in SRA01/04 cells. Endogenous GPX3-AS was downregulated in SRA01/04 cells transfected with three ASOs for human GPX3-AS; among them, ASO-#2 was the most efficient and was used for the subsequent experiments (Figure 4A). GPX3-AS was significantly upregulated in SRA01/04 cells transfected with the pcDNA-GPX3-AS plasmid (Figure 4D). Then, the cells were exposed to H2O2 for 24 h, and the viability of GPX3-AS knockdown cells decreased. After GPX3-AS was overexpressed, the opposite results were observed (Figure 4B,E). TUNEL staining showed that the proportion of apoptotic SRA01/04 cells increased following GPX3-AS knockdown, while GPX3-AS-overexpressing SRA01/04 cells exhibited less apoptosis (Figure 4C,F). The results of flow cytometry were similar to those of TUNEL staining (Figure 4G,H). Taken together, the findings indicated that GPX3-AS inhibited H2O2-induced SRA01/04 cell apoptosis.
GPX3-AS is transcribed from the AS strand of the GPX3 gene locus 5q33.1 . We assessed whether GPX3-AS affected GPX3 expression. We first measured GPX3 expression in both the anterior lens capsules of ARCs and SRA01/04 cells treated with H2O2. Like the expression of GPX3-AS, GPX3 expression also decreased (Figure 5A,B). GPX3 mRNA and protein levels increased in cells overexpressing GPX3-AS and decreased in cells following GPX3-AS knockdown (Figure 5C-F). These results suggested that GPX3 expression was upregulated by GPX3-AS.
To further explore the role of GPX3 in the GPX3-AS-mediated antiapoptotic effect on SRA01/04 cells, the endogenous expression of GPX3 was remarkably decreased in SRA01/04 cells by transfection with GPX3 siRNA, and siRNA1 was selected for the subsequent study (Figure 6A). However, GPX3-AS expression did not change after GPX3 siRNA transfection, indicating that GPX3-AS is upstream of GPX3 (Figure 6B). Before exposure to H2O2, SRA01/04 cells were transfected with GPX3 siRNA in the absence or presence of pcDNA-GPX3-AS. Twenty-four hours after transfection, SRA01/04 cells were treated with H2O2 for 24 h, and cell apoptosis was then assessed by flow cytometry. Cotransfection with pcDNA-GPX3-AS and GPX3 siRNA significantly reversed GPX3 mRNA levels induced by pcDNA-GPX3-AS transfection (Figure 6C). Compared with the pcDNA and siRNA controls, there was little change in the apoptosis of SRA01/04 cells cotransfected with pcDNA-GPX3-AS and GPX3 siRNA (Figure 6D,E). Furthermore, GPX3 knockdown in SRA01/04 cells resulted in decreased cell viability and an increased cellular apoptosis ratio. These results suggest that GPX3-AS inhibits lens epithelial cell apoptosis by upregulating GPX3.
It is generally accepted that the apoptosis of lens epithelial cells is associated with the development of cataracts, and oxidative stress is regarded as a major cause of lens opacity . Therefore, inhibiting LEC apoptosis and enhancing LEC antioxidant capacity is a promising strategy for the treatment of ARCs. In our study, we used H2O2-stimulated SRA01/04 cells for mimicking oxidative stress in vitro; these cells exhibited increased cell apoptosis, which is consistent with previous studies [8,28].
LncRNAs are a novel heterogeneous class of ncRNAs of more than 200 nucleotides in length and without protein-coding activity. They are important regulators of gene expression. The roles of lncRNA dysregulation in ocular diseases have been well verified . In the present study, we used high-throughput RNA sequencing to obtain lncRNA expression profiles in cataract lenses and transparent lenses and found that GPX3-AS expression was significantly decreased in the cataract lenses compared with that in the transparent lenses. Our study demonstrated that GPX3-AS was mainly expressed in the nuclei of SRA01/04 cells (an HLEC line). By using SRA01/04 cells as an in vitro model to mimic oxidative stress in vivo, we determined that GPX3-AS levels were downregulated and associated with cell apoptosis. An in vitro study demonstrated that knockdown of the lncRNA GPX3-AS weakened cell viability and promoted apoptosis, while the overexpression of GPX3-AS enhanced cell viability and inhibited cell apoptosis. Moreover, the inhibition of cell apoptosis by GPX3-AS overexpression was reversed by GPX3 siRNA. Hence, GPX3-AS is a negative regulator of ARCs; the underlying mechanism of this regulation is the upregulation of GPX3 expression by GPX3-AS, which thereby reduces apoptosis in SRA01/04 cells.
According to their relative gene locations, lncRNAs are typically categorized into the five following classes: AS, intergenic, overlapping, intronic, and full lapping . Among these categories of lncRNAs, AS lncRNAs, which are reverse complements of their endogenous sense counterparts, frequently do not encode proteins and comprise a substantial proportion of the entire long noncoding transcriptome [30-32]. Previous studies have shown that these natural AS transcripts play critical roles in multiple physiological and pathological processes, ranging from transcriptional regulation by gene promoter activation to posttranscriptional regulation by controlling mRNA stability and translation, through regulating sense gene expression [33,34]. For example, the AS lncRNAs keratin 7 (KRT7)-AS, fork head box C2 (FOXC2)-AS1 and zinc finger E-box binding homeobox 1 (ZEB1)-AS1 can alter the stability and expression of KRT7, FOXC2, and ZEB1, respectively [32,35,36]. Our study revealed that GPX3-AS upregulated GPX3 expression at both the mRNA and protein levels. KRT7-AS and KRT7 mRNA form a duplex RNA–RNA complex that protects KRT7 mRNA from degradation; thus, KRT7-AS upregulates KRT7 expression . Therefore, we speculate that GPX3-AS could bind to GPX3 mRNA to promote GPX3 expression, although this is a question for future studies.
Several studies have demonstrated that the tumor-suppressive function of GPX3 in prostate, colon, and lung cancers is mediated by the induction of apoptosis [26,37,38]. However, high GPX3 expression is related to decreased overall patient survival, and GPX3 facilitates ovarian cancer cell viability by protecting cells from extracellular factors of oxidative stress . Our study also found that GPX3 inhibited SRA01/04 cell apoptosis via protecting cells from oxidative stress. These data indicate dual roles for GPX3 in cell apoptosis depending on the cell type, oxidative stress status, and cell microenvironment. The detailed mechanisms underlying the inhibition of SRA01/04 cell apoptosis by GPX3 require further investigation. Taken together, our results suggest that GPX3-AS overexpression inhibits apoptosis in HLECs, possibly by upregulating GPX3 expression. The present study may provide a better understanding of lncRNA function in ARC development, and GPX3-AS could be a target for ARC therapies.
The study was supported by the National Natural Science Foundation of China (no. 81770906).