- Research article
- Open access
- Published:
MiRNAs regulate oxidative stress related genes via binding to the 3′ UTR and TATA-box regions: a new hypothesis for cataract pathogenesis
BMC Ophthalmology volume 17, Article number: 142 (2017)
Abstract
Background
Age-related cataracts are related to oxidative stress. However, the genome-wide screening of cataract related oxidative stress related genes are not thoroughly investigated. Our study aims to identify cataract regulated miRNA target genes that are related to oxidative stress and to propose a new possible mechanism for cataract formation.
Methods
Microarrays were used to determine the mRNA expression profiles of both transparent and cataractous lenses. The results were analyzed by significance analyses performed by the microarray software, and bioinformatics analysis was further conducted using Molecular Annotation System. The Eukaryotic Promoter Database (EPD) was used to retrieve promoter sequences and identify TATA-box motifs. Online resource miRWalk was exploited to screen for validated miRNAs targeting mRNAs related to oxidative stress. RNAhybrid online tool was applied to predict the binding between significantly regulated miRNAs in cataract lenses and target mRNAs.
Results
Oxidative stress pathway was significantly regulated in cataractous lens samples. Pro-oxidative genes were half up-regulated (11/20), with a small number of genes down-regulated (4/20) and the rest of them with no significant change (5/20). Anti-oxidative genes were partly up-regulated (17/69) and partly down-regulated (17/69). Four down-regulated miRNAs (has-miR-1207-5p, has-miR-124-3p, has-miR-204-3p, has-miR-204-5p) were found to target 3′ UTR of pro-oxidative genes and could also bind to the TATA-box regions of anti-oxidative genes (with the exception of has-miR-204-3p), whilst two up-regulated miRNAs (has-miR-222-3p, has-miR-378a-3p) were found to target 3′ UTR of anti-oxidative genes and could simultaneously bind to the TATA-box regions of pro-oxidative genes.
Conclusions
We propose for the first time a hypothesis that cataract regulated miRNAs could contribute to cataract formation not only by targeting 3′ UTR but also by targeting TATA-box region of oxidative stress related genes. This results in the subsequent elevation of pro-oxidative genes and inhibition of anti-oxidative genes. This miRNA-TATA-box/3′ UTR-gene-regulation network may contribute to cataract pathogenesis.
Background
Age-related cataract is till now still the dominant cause of blindness worldwide. Although cataract surgery is a satisfying solution for cataract caused vision impairment, this procedure is not easily available in developing countries [1]. There have been studies on the mechanism of age-related cataract aiming to discover new targets for cataract treatment. However, the specific molecular pathway of this disease still merits further investigation. It has been suggested that oxidative stress plays an important role in cataract formation. Lens proteins undergo non-enzymatic, post-translational modification, accumulate fluorescent chromophores and are more susceptible to oxidation along with aging [2], therefore oxidative stress is partially responsible for cataract formation. It has also been reported that oxidation plays a key role in nuclear cataract, but its causal roles in cortical and posterior subcapsular cataracts are substantially less important [3]. Other groups have published gene profiling results from either cataract tissue samples or cell lines [4, 5]. There are also studies on the function of specific genes in cataract formation, particularly oxidative stress related genes, such as SOD1, TXNIP, etc. [6,7,8]. Thereupon oxidative stress related gene regulation is worth investigating.
MicroRNAs (miRNAs) are a group of small non-coding RNAs that modulate many pathways related to various diseases. Our previous work has shown that multiple miRNAs are regulated in human cataractous lenses [9]. However, the targets of these miRNAs still warrant further investigation. In this study, we used gene profiling to screen for differentially regulated mRNAs in human cataractous lenses compared with post mortem clear lenses. Then bioinformatics analysis was applied to examine the gene-pathway distribution and the regulation of oxidative stress related genes in cataract lenses. We compared our previously published cataract regulated miRNAs and these differentially regulated mRNAs, screening for miRNA 3′ UTR binding targets that could contribute to cataract pathogenesis. Our previous work indicated miRNAs could bind to the TATA-box region and promote gene transcription as well as 3′ UTR binding-mediated gene silencing [10, 11], so we also screened for TATA-box targets of the miRNAs reported in our previous publication [9]. Our study establishes a miRNA-TATA-box/3′ UTR-gene-regulation network that may contribute to cataract formation.
Methods
Human lens epithelium-capsule sample collection
Lens samples were obtained from Department of Ophthalmology, the First Affiliated Hospital of Medical School of Xi’an Jiaotong University. Epithelium-capsule samples were obtained after capsulorrhexis during surgery (15 age-related cataract patients, age ranges from 60 to 68, no other ocular diseases) or obtained from postmortem transparent lenses within 24 h after death (15 donors, age ranges from 58 to 65, no history of ocular diseases). The degrees of lenticular opacification of patients were determined by the Lens Opacities Classification System III (LOCSIII) [12, 13]. The lenticular opacification of postmortem lenses were grade 1 ~ 2 and cataract lenses were grade 4 ~ 6 (Additional file 1: Table S1). There was no significant difference between transparent and cataractous human lens samples respect to the age of donors in two groups (Additional file 1: Table S1, unpaired t test, p > 0.05).
Total RNA preparation
Five samples were randomly pooled into one lens sample to insure the quality and quantity of extracted total RNA. Total RNA of six pooled samples (three cataractous and three transparent samples) was extracted with Trizol reagent (Invitrogen, Shanghai, China) according to the manufacturer’s instructions. The quality of RNA samples was confirmed by measuring OD260/280 ratio using NanoDrop spectrophotometer (Thermo Scientific, Shanghai, China) and the integrity of RNA samples was verified by agarose electrophoresis.
Microarray analysis
Gene expression profiling was performed for each pooled RNA sample using a Gene Chip-Human Genome Array (HG-U133 Plus 2.0, Affymetrix, Santa Clara, CA) at CapitalBio Corporation (Beijing, China). Gene Chip microarray service of CapitalBio Corporation is certificated by Affymetrix. Microarray processing was performed according to Gene Chip Expression Analysis Technical Manual provided by Affymetrix. In brief, 1 μg of total RNA was used to synthesize double-stranded cDNA. Biotin-tagged complementary RNA (cRNA) was produced using the MessageAmp II aRNA Amplification Kit (Ambion, Austin, TX). The resulting cRNA was fragmented and hybridized to the microarray. After hybridization the GeneChip arrays were washed, stained and then scanned on a GeneChip Scanner 3000 (Affymetrix).
Data processing and analysis
Microarray raw data were processed at CapitalBio Corporation. Briefly, the hybridization data were analyzed using GeneChip Operating software (GCOS 1.4, Affymetrix) and normalized using a DNA-chip analyzer. Significant Analysis of Microarray software (SAM) was used to identify mRNAs that exhibited significant differences in expression between transparent and cataractous samples (average fold change >2 or <0.5). Results were screened for mRNAs with FDR (false discovery rate)-corrected p values < 0.05.
Bioinformatics analysis
Bioinformatics analysis was conducted via the online Molecular Annotation System (MAS 3.0) provided by CapitalBio Corporation. Gene symbols of mRNAs with average fold change >2 of <0.5 were uploaded in the MAS 3.0 system for Gene Ontology (GO) and gene-pathway network analysis. Heatmaps were either provided by CapitalBio Corporation (Fig. 1) or generated by using Heatmap illustrator 1.0 according to the user’s manual (Fig. 4) [14]. The Eukaryotic Promoter Database (EPD) [15, 16] was used to retrieve promoter sequences of selected oxidative stress related mRNAs and to identify TATA-box motifs as described in our previous publications [10, 11]. Online resource miRWalk [17, 18] was used to screen for validated miRNAs targeting mRNAs related to oxidative stress. RNAhybrid online tool [19] was applied to predict the binding between significantly regulated miRNAs in cataract lenses in our previous findings [9] and target mRNAs.
Statistical analysis
Differences between the two groups in microarray data were assessed using FDR (false discovery rate)-corrected p values (p < 0.05) provided by SAM. Age differences between the two groups were estimated using unpaired t test (p < 0.05 was considered statistically significant).
Results
Gene profiling in microarrays shows oxidative stress pathway is regulated in cataractous lenses
Gene profiling was determined by six separate microarray assays, which identified 2223 up-regulated genes and 1829 down-regulated genes in cataractous lenses compared with transparent lenses (average fold change >2 or <0.5) (Fig. 1). Gene symbols of the top 100 up-regulated genes and top 100 down-regulated genes were uploaded to MAS 3.0 system to generate gene-pathway network graph using GenMAPP database, as it was too complex to generate graph by Graphviz if more than 200 gene symbols were included (Fig. 2). It is shown that oxidative stress related pathway was involved in cataract related gene regulation (Fig. 2). The overall distribution of cataract-regulated genes within different categories of the Gene Ontology (GO) classification system (the “biological process” and “pathway in Kegg database” categories) was examined (Fig. 3). Several classes of cataract-regulated genes in the biological process such as signal transduction, transcription, ion transport, etc. were consistent with previous work [4]. Interestingly, oxidation reduction was one of the top eight classes of GO enriched functional categories (Fig. 3a, Additional file 1: Table S2). We also found several pathways that were involved in cataract-mediated gene regulation, such as calcium signaling pathway, MAPK signaling pathway, TGF-beta signaling pathway, etc. (Fig. 3b).
Expression of oxidative stress pathway related genes in cataractous lenses
To further investigate the regulation of oxidative stress related genes in cataractous lenses, we generated heatmaps consisting only genes related to oxidative stress (Fig. 4). As it is shown there were more significantly up-regulated pro-oxidative genes (11 of 20) than down-regulated genes (4 of 20) in cataractous lenses (Fig. 4a, Additional file 1: Table S3). The rest 5 genes were not significantly regulated. So the oxidative stress related genes were generally up-regulated in cataractous samples. However, there were equal numbers of up-regulated and down-regulated anti-oxidative genes (17 of 69 for both) (Fig. 4b, Additional file 1: Table S3). The up-regulation of anti-oxidative genes was probably due to the response to elevated oxidative stress in cataractous lenses. The down-regulated anti-oxidative genes and up-regulated pro-oxidative genes could lead to the progression of oxidative stress, thus may contribute to cataract formation.
MiRNAs bind to the TATA-box/3′ UTR of oxidative stress related genes
Since miRNA has been suggested to play a role in cataract formation by our previous work as well as publications of other colleagues [9, 20, 21], we used miRWalk database [17, 18] to screen for miRNA-target oxidative stress related genes. Four down-regulated miRNAs (has-miR-1207-5p, has-miR-124-3p, has-miR-204-3p, has-miR-204-5p) were found to target 3′ UTR of pro-oxidative genes whilst two up-regulated miRNAs (has-miR-222-3p, has-miR-378a-3p) were found to target 3′ UTR of anti-oxidative genes (Tables 1 and 2). To further investigate the part played by these miRNAs, we retrieved the promoter sequences of oxidative stress related genes and predicted the binding between the aforementioned miRNAs and the TATA-box regions of promoters. We found three out of four down-regulated miRNAs could specifically bind to the TATA-box regions of anti-oxidative genes, whilst the two up-regulated miRNAs could target the TATA-box regions of pro-oxidative genes (Table 1). For instance, miR-204-5p could specifically bind to the 3′ UTR of TXNIP, which is a regulator of the bioavailability of thioredoxin in the lens, and may promote oxidative stress-induced apoptosis [7]. MiR-204-5p could also bind to the TATA-box region of ALDH1A3, which is a member of gene family protecting the eye against cataract formation via detoxification function [8]. Since miR-204-5p was down-regulated in cataractous lenses, it is hypothesized miR-204-5p up-regulates TXNIP expression through the reduction of post-transcriptional gene silencing, while down-regulates ALDH1A3 expression via reduced promoter activation-mediated transcription (Fig. 5a). On the other hand, miR-378a-3p could bind to the 3′ UTR of SOD1, a gene preventing hydrogen peroxide -induced oxidative damage to the lens [6]. MiR-378a-3p could also specifically bind to the TATA-box region of TXNIP, a pro-oxidative gene [7]. Because miR-378a-3p was up-regulated in cataractous lenses, we propose that miR-378a-3p down-regulates SOD1 expression via post-transcriptional gene silencing and up-regulates TXNIP expression through promoter activation-mediated transcription (Fig. 5b). Our results suggest up-regulated miRNAs down-regulate anti-oxidative genes via 3′ UTR binding, meanwhile up-regulate pro-oxidative genes via TATA-box binding-mediated transcription activation, and it is the opposite for down-regulated miRNAs. This results in the elevation of pro-oxidative genes and inhibition of anti-oxidative genes, which may lead to cataract (Fig. 6).
Discussion
Age-related cataract is believed to be the result of post-translational modification, the accumulation of fluorescent chromophores, increasing susceptibility to oxidation, etc. [2]. Oxidative stress plays an important role in nuclear cataract formation, therefore in this study we exploited genome microarray to determine the genes regulated in cataractous lenses, especially genes related with oxidative stress. Pro-oxidative genes were nearly half up-regulated (11/20), with a small number of genes down-regulated (4/20) and the rest of them with no significant change (5/20). Interestingly, anti-oxidative genes were partly up-regulated (17/69) and partly down-regulated (17/69). The possible explanation of this discrepancy is oxidative stress is not the only cause of nuclear cataract and other mechanism might contribute to this disease as well. Furthermore, there may be a compensation mechanism towards oxidative damage, and the up-regulation of some anti-oxidative genes and the down-regulation of some pro-oxidative genes may act as a response to oxidative stress. Finally, there are limitations in our study such as the limited number of genes selected related to oxidative stress, lack of experimental validations, etc.
MiRNAs have been linked to cataract pathogenesis in studies from us and others as well [9, 20,21,22,23,24,25]. It would be worth exploring the relationship between these miRNAs and oxidative stress related genes. The classic pathway of miRNA regulated gene expression is via 3′ UTR binding-mediated post-transcription gene silencing. Our previous work indicated miRNAs could also bind to the TATA-box region and act as a gene transcription activator [10, 11]. Therefore, we took advantage of bioinformatics web tools and online resources to screen for targets of cataract regulated miRNAs. In this study, we found cataract regulated miRNAs could contribute to cataract formation not only by targeting 3′ UTR, but also by binding to the TATA-box region of oxidative stress related genes. This resulted in the subsequent elevation of pro-oxidative genes and inhibition of anti-oxidative genes. The elevated level of oxidative stress may lead to cataract formation.
Conclusions
In conclusion, we propose a hypothesis that this miRNA-TATA-box/3′ UTR-gene-regulation network may contribute to cataract pathogenesis. Our next step would be to validate these aforementioned miRNAs and mRNAs, to certify the target relationship between these miRNAs and their corresponding mRNAs, and to elucidate the specific mechanism of miRNA-TATA-box/3′ UTR-gene-regulation network during cataract formation.
Abbreviations
- cRNA:
-
Biotin-tagged complementary RNA
- EPD:
-
The Eukaryotic Promoter Database
- FDR:
-
False discovery rate
- GO:
-
Gene ontology
- LOCSIII:
-
Lens Opacities Classification System III
- MAS 3.0:
-
Molecular Annotation System
- miRNAs:
-
MicroRNAs
- SAM:
-
Significant Analysis of Microarray software
References
Vrensen GF. Early cortical lens opacities: a short overview. Acta Ophthalmol. 2009;87(6):602–10.
Michael R, Bron AJ. The ageing lens and cataract: a model of normal and pathological ageing. Philos Trans R Soc Lond Ser B Biol Sci. 2011;366(1568):1278–92.
Beebe DC, Holekamp NM, Shui YB. Oxidative damage and the prevention of age-related cataracts. Ophthalmic Res. 2010;44(3):155–65.
Ruotolo R, Grassi F, Percudani R, Rivetti C, Martorana D, Maraini G, Ottonello S. Gene expression profiling in human age-related nuclear cataract. Mol Vis. 2003;9:538–48.
Zhao W, Zhao J, Wang D, Li J. Screening of potential target genes for cataract by analyzing mRNA expression profile of mouse Hsf4-null lens. BMC Ophthalmol. 2015;15:76.
Lin D, Barnett M, Grauer L, Robben J, Jewell A, Takemoto L, Takemoto DJ. Expression of superoxide dismutase in whole lens prevents cataract formation. Mol Vis. 2005;11:853–8.
Liyanage NP, Fernando MR, Lou MF. Regulation of the bioavailability of thioredoxin in the lens by a specific thioredoxin-binding protein (TBP-2). Exp Eye Res. 2007;85(2):270–9.
Lassen N, Bateman JB, Estey T, Kuszak JR, Nees DW, Piatigorsky J, Duester G, Day BJ, Huang J, Hines LM, et al. Multiple and additive functions of ALDH3A1 and ALDH1A1: cataract phenotype and ocular oxidative damage in Aldh3a1(−/−)/Aldh1a1(−/−) knock-out mice. J Biol Chem. 2007;282(35):25668–76.
Wu C, Lin H, Wang Q, Chen W, Luo H, Zhang H. Discrepant expression of microRNAs in transparent and cataractous human lenses. Invest Ophthalmol Vis Sci. 2012;53(7):3906–12.
Zhang Y, Yin Y, Zhang S, Luo H, Zhang H. HIV-1 infection-induced suppression of the let-7i/IL-2 Axis contributes to CD4(+) T cell death. Sci Rep. 2016;6:25341.
Zhang Y, Fan M, Zhang X, Huang F, Wu K, Zhang J, Liu J, Huang Z, Luo H, Tao L, et al. Cellular microRNAs up-regulate transcription via interaction with promoter TATA-box motifs. RNA. 2014;20(12):1878–89.
Chylack LT Jr, Wolfe JK, Singer DM, Leske MC, Bullimore MA, Bailey IL, Friend J, McCarthy D, Wu SY. The lens opacities classification system III. The longitudinal study of cataract study group. Arch Ophthalmol. 1993;111(6):831–6.
Grewal DS, Brar GS, Grewal SP. Correlation of nuclear cataract lens density using Scheimpflug images with lens opacities classification system III and visual function. Ophthalmology. 2009;116(8):1436–43.
Deng W, Wang Y, Liu Z, Cheng H, Xue Y. HemI: a toolkit for illustrating heatmaps. PLoS One. 2014;9(11):e111988.
Dreos R, Ambrosini G, Perier RC, Bucher P. The eukaryotic promoter database: expansion of EPDnew and new promoter analysis tools. Nucleic Acids Res. 2015;43(Database issue):D92–6.
Dreos R, Ambrosini G, Cavin Perier R, Bucher P. EPD and EPDnew, high-quality promoter resources in the next-generation sequencing era. Nucleic Acids Res. 2013;41(Database issue):D157–64.
Dweep H, Sticht C, Pandey P, Gretz N. miRWalk--database: prediction of possible miRNA binding sites by “walking” the genes of three genomes. J Biomed Inform. 2011;44(5):839–47.
Dweep H, Gretz N. miRWalk2.0: a comprehensive atlas of microRNA-target interactions. Nat Methods. 2015;12(8):697.
Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R. Fast and effective prediction of microRNA/target duplexes. RNA. 2004;10(10):1507–17.
Qin Y, Zhao J, Min X, Wang M, Luo W, Wu D, Yan Q, Li J, Wu X, Zhang J. MicroRNA-125b inhibits lens epithelial cell apoptosis by targeting p53 in age-related cataract. Biochim Biophys Acta. 2014;1842(12 Pt A):2439–47.
Dong Y, Zheng Y, Xiao J, Zhu C, Zhao M. MicroRNA let-7b induces lens epithelial cell apoptosis by targeting leucine-rich repeat containing G protein-coupled receptor 4 (Lgr4) in age-related cataract. Exp Eye Res. 2016;147:98–104.
Li Y, Liu S, Zhang F, Jiang P, Wu X, Liang Y. Expression of the microRNAs hsa-miR-15a and hsa-miR-16-1 in lens epithelial cells of patients with age-related cataract. Int J Clin Exp Med. 2015;8(2):2405–10.
Dong N, Tang X, Xu B. miRNA-181a inhibits the proliferation, migration, and epithelial-mesenchymal transition of lens epithelial cells. Invest Ophthalmol Vis Sci. 2015;56(2):993–1001.
Dong N, Xu B, Benya SR, Tang X. MiRNA-26b inhibits the proliferation, migration, and epithelial-mesenchymal transition of lens epithelial cells. Mol Cell Biochem. 2014;396(1–2):229–38.
Wang Y, Li W, Zang X, Chen N, Liu T, Tsonis PA, Huang Y. MicroRNA-204-5p regulates epithelial-to-mesenchymal transition during human posterior capsule opacification by targeting SMAD4. Invest Ophthalmol Vis Sci. 2013;54(1):323–32.
Acknowledgements
Not applicable.
Funding
The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results. This study was supported by the Fundamental Research Funds for the Basic Research Operating Expenses Program of Central College sponsored by Xi′an Jiaotong University (xjj2013067), Youth Foundation of the First Affiliated Hospital, Medical College, Xi′an Jiaotong University (2014YK7) and the National Natural Science Foundation of China (No. 81602698, No. 81470614).
Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files.
Author information
Authors and Affiliations
Contributions
YY and CW conceived and designed the experiments; CW performed the experiments; CW and YY analyzed the data; ZL, LM, CP, LQ, NG and JL performed capsulorrhexis during cataract surgery and provided clinical samples for this study; CW and YY wrote the paper. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Authors’ information
Not applicable.
Ethics approval and consent to participate
This research was approved by the Ethics Review Board of Xi’an Jiaotong University and performed in accordance with the tenets of the Declaration of Helsinki.
Consent for publication
Not applicable.
Competing interests
The authors declare they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file
Additional file 1: Table S1.
Degrees of Lenticular Opacification Determined by Lens Opacities Classification System III (LOCSIII). Table S2. Top 8 Populated Gene Ontology Terms of the Biological Process Category. Table S3. Selected Oxidative Stress Related Genes (DOCX 23 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Wu, C., Liu, Z., Ma, L. et al. MiRNAs regulate oxidative stress related genes via binding to the 3′ UTR and TATA-box regions: a new hypothesis for cataract pathogenesis. BMC Ophthalmol 17, 142 (2017). https://doi.org/10.1186/s12886-017-0537-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12886-017-0537-9