Skip to main content

RNA sequencing and bioinformatics analysis of human lens epithelial cells in age-related cataract



Age-related cataract (ARC) is the main cause of blindness in older individuals but its specific pathogenic mechanism is unclear. This study aimed to identify differentially expressed genes (DEGs) associated with ARC and to improve our understanding of the disease mechanism.


Anterior capsule samples of the human lens were collected from ARC patients and healthy controls and used for RNA sequencing to detect DEGs. Identified DEGs underwent bioinformatics analyses, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. Subsequently, reverse transcription quantitative RT-qPCR was used to validate the different expression levels of selected genes.


A total of 698 up-regulated DEGs and 414 down-regulated DEGs were identified in ARC patients compared with controls by transcriptome analysis. Through GO and KEGG bioinformatics analysis, the functions of significantly DEGs and their possible molecular mechanisms were determined. Sequencing results were verified by RT-qPCR as being accurate and reliable.


This study identified several genes associated with ARC, which improves our knowledge of the disease mechanism.

Peer Review reports


Age-related cataract (ARC) is a common visual disorder caused by lens opacity, and the main cause of blindness in individuals aged over 50 years [1]. With an increasingly aging society, the incidence of ARC is rising sharply [2]. Currently, surgery is the only effective treatment for ARC, but this is associated with an economic burden for patients and may cause complications [3]. Therefore, alternative effective methods are necessary to prevent and treat the occurrence and development of ARC.

Risk factors such as genetics, oxidative stress, ultraviolet radiation, and drug use were previously shown to be closely related to the pathogenesis of ARC [4]. Of these, changes in the expression of multiple genes under oxidative stress damage seem to be key events [5]. Lens epithelial cells are the only nucleated cells in the lens, and as such contain most of the metabolic, protective, permeable, and other regulatory systems [6]. They communicate with the underlying lens fiber cells and respond to cataract-related damage by altering gene expression [7]. Although several genes associated with the pathogenesis of ARC have been identified, many have not been explored. Therefore, the identification of candidate genes for ARC progression in lens epithelial cells will help further understand the disease mechanism and provide an important basis for reducing disease occurrence and for targeted therapy.

In the present study, we performed RNA sequencing on human lens anterior capsule tissue and compared transcriptome data between ARC patients and healthy controls. Differentially expressed genes (DEGs) were analyzed using bioinformatics such as Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. The accuracy of RNA sequencing was verified using reverse transcription quantitative RTqPCR for 10 selected DEGs. Several genes were identified as being associated with ARC, which furthers our knowledge of the disease mechanism.


Patient and tissue sample collection

Twelve cases of age-related cataract anterior capsule tissue samples were collected from Hongqi Hospital affiliated to Mudanjiang Medical College from January 2019 to March 2019. All the patients were cortical cataracts, including 6 males and 6 females, aged 55 to 64 years. The diagnostic criteria for ARC are non-congenital cataract, non-posterior cataract, and non-metabolic cataract, no hypertension, no diabetes, no glaucoma, no fundus lesions, no uveitis, no eye trauma, and no intraocular surgery history. Twelve cases of normal anterior capsule tissue samples were collected from corneal donors and patients with lens trauma. The information of the samples was displayed in Table 1. Human anterior lens capsule tissue samples were obtained by an ophthalmologist through a continuous annular capsulorhexis of the lens and randomly divided 4 samples into a group. Tissue samples were immediately frozen in liquid nitrogen and stored at − 80 °C for RNA sequencing. This study followed the principles of the Declaration of Helsinki and all participants signed informed consent.

Table 1 The information of our patients and normal samples

RNA extraction

Since the miRNA isolation kit can isolate total RNA containing small RNAs from cultured cells or tissues for better extraction, we selected the mirVana™ miRNA ISOlation Kit (Ambion-1561) to isolate total RNA from tissue samples according to the manufacturer’s protocol. The RNA quality was determined by OD 260/280 absorbance assessment using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, UK). RNA integrity was evaluated using agarose gel electrophoresis stained with ethidium bromide. All RNA samples were stored at − 80 °C for further application.

RNA library preparation and transcriptome sequencing

The 12 healthy controls samples were randomly divided into 3 groups, with 4 samples in each group. Similarly, the 12 samples from the ARC patients were also randomly divided into 3 groups of 4 samples per group. Since the sample tissue of the anterior lens capsule we used was very small, and less RNA was extracted from a single sample, we combined the samples for RNA sequencing analysis. The total RNA extracted from each group was about 0.8 μg, and the RIN value was above 7. Therefore, we used a low start volume library construction method (NEBNext® Ultra™ II Directional RNA Library Prep Kit for Illumina) for library construction, which includes mRNA enrichment analysis. Briefly, RNA Library Preparation mainly has the following steps: RNA purification and fragmentation, Synthesize First Strand cDNA, Synthesize Second Strand cDNA, Adenylate 3′Ends, Ligate Adapters, Enrich DNA Fragments. The Agilent bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA) was used to check the size and purity of the library. Sequencing was performed on Illumina HiSeq X Ten to generate 125 bp/150 bp paired-end reads.

Sequence data analysis

The raw data in the fastq format (raw reads) was first processed using Trimmomatic. Then, clean data (clean reads) were obtained by removing low-quality reads, reads containing adapter, and poly-N from raw data. Q30 content represents the percentage of bases with a Qphred value greater than 30 to the total bases. All subsequent analyses were performed using clean data with high quality. HISAT2 was used to compare the sequence of Clean Reads with the specified reference genome [8]. Fragments Per kb Per Million Reads (FPKM) were calculated to quantify gene expression levels. Differentially expressed genes (DEGs) analysis was performed using the DESeq R package, which based DEGs on the negative binomial distribution. The screening conditions for DEGs were P < 0.05 with fold change (FC) > 2.

GO and KEGG enrichment analysis of DEGs

Gene Ontology (GO) analysis ( was used to categorize the function of the DEGs. And the functional description was mainly divided into biological process, cellular component, and molecular function. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis [9] ( was used to predict the signaling pathways in which these DEGs may be involved, KEGG software from the Kanehisa laboratory was used with prior permiossion. The significantly enriched GO and KEGG terms met the criterion of corrected P < 0.05.

Quantitative real-time PCR validation

The remaining RNA after transcriptome was further used for validation studies by qPCR. Total RNA was extracted from the anterior capsule of the normal lens and cataract patients for RNA transcriptome sequencing. Detected the concentration of RNA with a NanoDrop spectrophotometer. Agarose gel electrophoresis detected the integrity of RNA. Quantification was carried out through a two-step reaction process of reverse transcription (RT) and PCR. Each RT reaction consisted of a total of 10 μl of 0.5 μg RNA, 2 μl of 5x Transx all-in-one SuperMix for qPCR, and 0.5 μl of gDNA remover. The reaction was performed on GeneAmp® PCR System 9700 (Applied Biosystems, USA) at 42 °C, 5 s, and 85 °C for 15 min. Next, 10 μl of the RT reaction mixture was diluted 10-fold with nuclease-free water and kept at − 20 °C. Use LightCycler®480II Real-Time PCR instrument (Roche, Switzerland) with 1 μl cDNA, 5 μl 2xPerfectStartTM Green qPCR SuperMix, 0.2 μl forward primer, 0.2 μl PCR reverse primer, and 3.6 μl nuclease-free water. The reaction was incubated in a 384-well optical plate (Roche, Switzerland) at 94 °C for 30 s, and then at 94 °C for 5 s and 60 °C for 30 s for 45 cycles. The primer sequence is synthesized according to the mRNA sequence obtained from the NCBI database (Table 2). Each sample was analyzed in triplicate. The 2-ΔΔCt method was used for data quantification [10].

Table 2 Primer sequences of the selected DEGs

Statistical analysis

The statistical significance of the experimental data was conducted by a Student’s t-test. The acceptable significance was set as *p < 0.05; **p < 0.01.


Descriptive analysis of sequencing data quality assessment

Anterior capsule tissue samples were collected from ARC patients (Disease Group, JBZ) and healthy controls (Control group, ZCR). Figure 1 shows an example of a cortical cataract visible as a cloudy lens cortex in a patient with ARC.

Fig. 1

Representative slit-lamp image of the lens with cortical opacity of cataract. From the representative patient, the slit lamp biomicroscopy examination showed the cortical cataract with age. Because the normal subjects with transparent lens were enrolled from corneal transplant donors and ocular trauma patients with lens detachment, a picture of normal control cannot be obtained

We used the Illumina platform to obtain paired-end sequencing data. Because data error can have a major impact on the results, we employed Trimmomatic software to perform quality preprocessing of the original data, and statistically summarize the number of reads throughout the quality control process (Table 3). A total of 55.59 million raw reads were generated from the JBZ1 group, 48.28 million from JBZ2, 48.25 million from JBZ3, 58.89 million from ZCR1, 47.17 million from ZCR2, and 47.97 million from ZCR3. A total of 52.73, 44.64, and 44.59 million clean reads were obtained for JBZ1–3, respectively, and 55.68, 42.94, and 42.95 million for ZCR1–3, respectively. Clean reads ratios were more than 89% for both patients and controls, and over 77% of clean reads could be mapped to human genome sequences. The highest proportion of sequences with only one alignment position on the reference sequence was 89.10% in the JBZ group and 90.47% in the ZCR group. Both JBZ and ZCR groups had approximately 3% multiple alignment positions, and both groups exceeded 90% for Q30 content, representing the percentage of bases with a Qphred value greater than 30. These data indicate that high-quality sequencing was achieved in this study, so subsequent transcriptional analysis could continue.

Table 3 Statistics of the sequencing data quality results from all samples

Analysis of differentially expressed genes (DEGs)

According to the criteria of P < 0.05 and FC > 2, a total of 1112 DEGs were identified between patients and controls. Of these 698 were up-regulated and 414 were down-regulated in the JBZ group compared with the ZCR group (Fig. 2a). Volcano plots, MA plots, and heatmaps of hierarchical clustering showed that the gene expression levels were distinguished (Fig. 2b–d). The top 20 up-regulated genes and top 20 down-regulated genes are listed in Table 4.

Fig. 2

Differentially expressed genes (DEGs) in anterior capsule tissues between the two groups. a The number of differentially expressed genes (DEGs) in the JBZ group and ZCR group (P < 0.05, Fold change (FC) > 2). The red histogram represented upregulated gens and a blue histogram represented downregulated gens. b The MA plots visualizing the DEGs between the two groups. The X-axis represented the average expression of all samples for comparison after normalization; the Y-axis represented log2FoldChange. The red dots above the red line represented the up-regulated DEGs, and the red dots below the red line represented the down-regulated DEGs. c The volcano plot of the DEGs between two groups. Red dots represented significantly up-regulated differential genes, green dots represented significantly down-regulated differential genes, and gray dots represented non-significantly differential genes. d Heat-maps of hierarchical clustering analysis of DEGs between two groups. The horizontal ordinate represents the clustering of samples from the JBZ and ZCR groups. The longitudinal coordinates represent the DEGs and the clustering of genes. The red color indicates up-regulated genes, while the blue color indicates down-regulated genes

Table 4 Top 20 upregulated and downregulated differentially expressed genes (DEGs) in the JBZ group compared to the ZCR group

Analysis of GO enrichment

GO is used to annotate gene function and standardize the description of gene products according to biological process, cellular component, and molecular function categories. We performed GO enrichment analysis by screening GO terms with more than two DEGs in the three gene product classifications, then sorting based on P < 0.05 from the largest to the smallest –log10p value, and selecting the top 10 in each category. A total of 848 DEGs (76.3% of all identified DEGs) were mapped to 3646 GO terms, and the top 30 GO terms are shown in Fig. 3. Of these, the most significantly enriched in biological process and molecular function categories were trigeminal nerve structural organization (GO: 0021637, P = 9.01E-06) and the structural constituent of the eye lens (GO: 0005212, P = 5.51E-05), respectively, while the most enriched in the cellular component category were the plasma membrane (GO: 0005886, P = 4.65E-07) and the integral component of the plasma membrane (GO: 0005887, P = 3.96E-06). These findings confirmed the accuracy of our sequencing results.

Fig. 3

GO enrichment analysis of differentially expressed genes (DEGs). The GO enrichment analysis of the top 30 GO terms. The X-axis represented the different GO terms, the Y-axis represented -log10Pvalue

Analysis of the KEGG pathway

KEGG is a database for the pathway analysis of DEGs to understand the biological functions of identified genes. In this study, KEGG was divided into the following six classifications: Cellular processing (four categories), Environmental information processing (three categories), Genetic information processing (four categories), Human diseases (eight categories), Metabolism (11 categories), and Organismal systems (10 categories). A total of 409 DEGs (36.8% of all identified DEGs) were assigned to 286 KEGG pathways. The top 20 KEGG enrichment pathways are shown in Fig. 4, and include transforming growth factor (TGF)-β, p53, mitogen-activated protein kinase (MAPK), tumor necrosis factor, and other common signaling pathways. A cancer-related pathway had the highest number of DEGs, suggesting that it might be associated with cell proliferation and apoptosis. KEGG analysis indicated that ARC involves multiple processes and complex signaling pathways, which require further exploration.

Fig. 4

KEGG enrichment analysis of differentially expressed genes (DEGs). As per the citation guidelines (, KEGG analysis was used to predict the signaling pathways in which these DEGs may be involved [9], KEGG software from the Kanehisa laboratory was used with prior permiossion. The KEGG enrichment analysis of the top 20 was shown in the figure. The number of genes was represented by the size of the circle, and the P-value was represented by the color

RT-qPCR validation

To confirm RNA sequencing results, 10 DEG transcripts were selected for validation using RT-qPCR. DEG expression levels in the ARC group compared with the control group were consistent with RNA sequencing data, indicating the accuracy of our findings (Fig. 5).

Fig. 5

qRT-PCR validation of ten differentially expressed genes. The relative expression levels of genes in ARC patients and normal samples. (*p < 0.05; **P < 0.01)


ARC is the leading cause of blindness worldwide, and its incidence is increasing annually. Till date, ARC pathogenesis has mainly been studied at the molecular and cellular levels, with particular focus on the epigenetic changes occurring for multiple genes [11]. Although many genes associated with ARC have been identified, several are undiscovered. In the present study, we conducted in-depth research into ARC-related genes using the bioinformatics technology of transcriptome sequencing.

Here, our study showed that the clean ratio in both groups exceeds 89%, and the highest mapping ratio is 93.5%. In all samples, we identified a total of 17,348 genes, including 1112 genes whose expression level was 2 times or more different than that of the transparent lens. A comparative transcriptome analysis between normal and ARC samples showed that 698 DEGs in ARC were up-regulated and 414 DEGs were down-regulated. These data can provide necessary tissue-specific insights for the ARC transcriptome, such as the number, the expression level, and the distribution of DEGs. It can be used to examine the expression of target genes in ARC, which helps explain genetic-related research. Besides, it can also be used to select candidate genes that are more functionally related to ARC development.

As the only transparent, avascular, and nerveless organ in the human body, the lens is composed of three parts: the lens capsule, lens epithelial cells, and lens fibers [12]. Lens epithelial cells are the only cells that contain a nucleus capable of aerobic metabolism and are closely connected to the anterior lens capsule. Because the stability of lens epithelial cells is important for maintaining lens transparency and metabolic homeostasis [13], we used these cells for sequencing analysis. We obtained a clean reads ratio in both ARC patients and controls of more than 89%, with a highest mapping ratio of 93.5%. From all samples, we identified a total of 17,348 genes, including 1112 whose expression differed by more than 2-fold than that of the transparent lens. Comparative transcriptome analysis between ARC and control samples showed that 698 DEGs were up-regulated and 414 were down-regulated in ARC. These data provide important tissue-specific insights for the ARC transcriptome, and can be used to examine the expression of target genes that are more functionally related to ARC development.

GO term analysis showed that structural constituents of the eye lens were highly enriched in the present study, which undoubtedly verifies the reliability of our results. DEGs in this GO term included BFSP1, LIM2, CRYBA4, CRYBA2, MIP, TMOD1, GJA1, LIM2, HSF4, and WNT2B, several of which are closely associated with cataract development. Additionally, some genes that we identified were previously reported to be related to the pathogenesis of ARC, such as CRYA2 and HSF4 [14, 15], and we have also verified these genes used RT-qPCR that are consistent with the results. Our transcriptome results also allowed us to identify and verify some genes rarely associated with ARC. Yan et al. found that the increased expression of laminin α4 in the anterior lens capsule is the main mechanism leading to the occurrence of ARC [16], which is consistent with our transcriptome sequencing results. Moreover, the inhibition of NLRP3 was shown to reduce oxidative stress-induced damage to human lens epithelial cells, while NLRP3 protein expression was significantly increased in lens epithelial cells under oxidative stress [17], which was also consistent with our verification results. However, more than 80% of the genes identified by transcriptome sequencing in this study have not previously been associated with ARC, so should be explored in future in-depth studies.

Some genes of interest identified here could be studied as future candidate genes of ARC. As a member of the cell adhesion molecule (CAM) family, the NrCAM gene plays a key role in the development of the nervous system [18, 19]. NrCAM expression was reported to be essential for maintaining contact between lens cells, and a lack of NrCAM caused lens fiber disorder which ultimately led to the formation of mouse cataracts [20]. However, NrCAM expression changes have not been investigated in ARC. We showed by transcriptome sequencing that NrCAM was down-regulated in ARC patients compared with controls, but RT-qPCR found it to be elevated in patients, which we speculate might reflect its low expression level. With this in mind, we aim to explore its function in ARC in future studies. We also identified potential ARC candidate genes previously shown to be involved in the occurrence and development of a variety of eye diseases, but not ARC, such as ANGPTL4. Aberrant ANGPTL4 expression levels are associated with branch retinal vein occlusion, macular edema, and age-related macular degeneration diseases, and are considered to be a biomarker and therapeutic target for retinopathy caused by ischemia [21, 22]. However, no studies have directly linked ANGPTL4 with ARC, so we will investigate this in our future research.

The apoptosis of lens epithelial cells is a common cellular basis for the occurrence and development of ARC [23]. A number of common apoptosis pathways such as p53, MAPK, and TGF-β signaling pathways are closely related to the occurrence of ARC, and several appear in our top 20 KEGG pathways [24,25,26]. Some genes of interest in these pathways include DUSP1, ID1, and CAV2, which all function in multiple signaling pathways [27,28,29]. We also verified their expression by RT-qPCR and obtained results consistent with transcriptome sequencing analysis.


The present study aimed to explore the transcription of genes expressed in the human anterior lens capsule of ARC patients to aid future research into the pathogenesis of ARC and novel therapies. We identified 1112 DEGs between ARC patients and controls using RNA sequencing technology of which many were novel as reported from this study. Although the functions of most identified DEGs have not been explored in ARC, we selected several potential ARC candidate genes for future studies.

Availability of data and materials

All data files are available from the NCBI database (accession number PRJNA678236),



Age-related cataract


Gene Ontology


Kyoto Encyclopedia of Genes and Genomes


Differentially expressed genes


Reverse transcription-quantitative polymerase chain reaction


Cell adhesion molecules


Fold change


Branch retinal vein occlusion


Macular edema


Age-related macular degeneration


  1. 1.

    Yao L, Yang L, Song H, et al. MicroRNA miR-29c-3p modulates FOS expression to repress EMT and cell proliferation while induces apoptosis in TGF-beta2-treated lens epithelial cells regulated by lncRNA KCNQ1OT1. Biomed Pharmacother. 2020;129:110290.

    CAS  Article  Google Scholar 

  2. 2.

    Knopfel EB, Vilches C, Camargo SMR, et al. Dysfunctional LAT2 amino acid transporter is associated with cataract in mouse and humans. Front Physiol. 2019;10:688.

    Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Ono K, Hiratsuka Y, Murakami A. Global inequality in eye health: country-level analysis from the global burden of disease study. Am J Public Health. 2010;100(9):1784–7.

    Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Wang Z, Su D, Sun Z, Liu S, sun L, Li Q, et al. MDM2 phosphorylation mediates H2O2-induced lens epithelial cells apoptosis and age-related cataract. Biochem Biophys Res Commun. 2020;528(1):112–9.

  5. 5.

    Asbell PA, Dualan I, Mindel J, et al. Age-related cataract. Lancet. 2005;365(9459):599–609.

    Article  PubMed  Google Scholar 

  6. 6.

    Hawse JR, Hejtmancik JF, Huang QL, et al. Identification and functional clustering of global gene expression differences between human age-related cataract and clear lenses. Mol Vis. 2003;9:515–37.

  7. 7.

    Hejtmancik JF, Kantorow M. Molecular genetics of age-related cataract. Exp Eye Res. 2004;79(1):3–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Kantorow M, Kays T, Horwitz J, et al. Differential display detects altered gene expression between cataractous and normal human lenses. Invest Ophthalmol Vis Sci. 1998;39(12):2344–54.

    CAS  PubMed  Google Scholar 

  12. 12.

    Vendra VPR, Khan I, Chandani S, et al. Gamma crystallins of the human eye lens. Biochimica et Biophysica Acta. 2016;1860(1 Pt B):333–43.

    CAS  Article  Google Scholar 

  13. 13.

    Ghosh KS, Chauhan P. Crystallins and their complexes. Subcell Biochem. 2019;93:439–60.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Puk O, Ahmad N, Wagner S, Hrabé de Angelis M, Graw J. First mutation in the betaA2-crystallin encoding gene is associated with small lenses and age-related cataracts. Invest Ophthalmol Vis Sci. 2011;52(5):2571–6.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Shi Y, Shi X, Jin Y, et al. Mutation screening of HSF4 in 150 age-related cataract patients. Mol Vis. 2008;14:1850–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Yan Y, Yu HY, Sun LY, Liu H, Wang C, Wei X, et al. Laminin α4 overexpression in the anterior lens capsule may contribute to the senescence of human lens epithelial cells in age related cataract. Aging. 2019;11(9):2699–723.

  17. 17.

    Zou YY, Cui BJ, Liang P, Tian X, Ma Y, Zhao S. Inhibition of NLRP3 protects human Lens epithelial cells against oxidative stress-induced apoptosis by NF-kappaB signaling. Ophthalmic Res. 2020;63(2):174–81.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Grurnet M, Mauro V, Burgoon MP, et al. Structure of a new nervous system glycoprotein, nr-CAM, and its relationship to subgroups of neural cell adhesion molecules. J Cell Biol. 1991;113(6):1399–412.

    Article  Google Scholar 

  19. 19.

    Kayyem JF, Roman JM, de la Rosa Ei J, et al. Bravo/nr-CAM is closely related to the cell adhesion molecules L1 and ng-CAM and has a similar heterodimer structure. The Journal of CeU Biolog. 1992;118(5):1259–70.

    CAS  Article  Google Scholar 

  20. 20.

    More MI, Kirsch FP, Rathjen FG. Targeted ablation of NrCAM or ankyrin-B results in disorganized lens fibers leading to cataract formation. J Cell Biol. 2001;154(1):187–96.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Kim JH, Shin JP, Kim IT, et al. Aqueous Angiopoietin-like 4 levels correlate with nonperfusion area and macular edema in branch retinal vein occlusion. Invest Ophthalmol Vis Sci. 2016;57(1):6–11.

    CAS  PubMed  Google Scholar 

  22. 22.

    Kim JH, Shin JP, Kim IT, Park DH. Angiopoietin-like 4 correlates with response to intravitreal ranibizumab injections in neovascular age-related macular degeneration. Retina. 2018;38(3):523–30.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Li WC, Kuszak JR, Dunn K, Wang RR, Ma W, Wang GM, et al. Lens epithelial cell apoptosis appears to be a common cellular basis for non-congenital cataract development in humans and animals. J Cell Biol. 1995;130(1):169–81.

  24. 24.

    Qin Y, Zhao JY, Min XJ, et al. MicroRNA-125b inhibits lens epithelial cell apoptosis by targeting p53 in age-related cataract. Biochimica et Biophysica Acta. 2014;1842(12 Pt A):2439–47.

    CAS  Article  Google Scholar 

  25. 25.

    Du SS, Shao JZ, Xie DD, et al. Decorin inhibits glucose-induced lens epithelial cell apoptosis via suppressing p22phox-p38 MAPK signaling pathway. PLoS One. 2020;15(4):e0224251.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Hales AM, Chamberlain CG, McAvoy JW. Susceptibility to tgfb2-induced cataract increases with aging in the rat. Invest Ophthalmol Vis Sci. 2000;41(11):3544–51.

    CAS  PubMed  Google Scholar 

  27. 27.

    Yang JX, Sun LG, Han J, Zheng W, Peng W. DUSP1/MKP-1 regulates proliferation and apoptosis in keratinocytes through the ERK/Elk-1/Egr-1 signaling pathway. Life Sci. 2019;223:47–53.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Huang YH, Hu J, Chen F, Lecomte N, Basnet H, David CJ, et al. ID1 mediates escape from TGFbeta tumor suppression in pancreatic Cancer. Cancer Discov. 2020;10(1):142–57.

  29. 29.

    Liu F, Shangli Z, Hu ZL. CAV2 promotes the growth of renal cell carcinoma through the EGFR/PI3K/Akt pathway. Onco Targets Ther. 2018;11:6209–16.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


The authors thank all patients and their family members for their participation.


This study was supported by the National Natural Science Foundation of China (31601116), the Postdoctoral Science Foundation of China (2016 M601447), the Basic research projects of provincial institutions of higher learning in Heilongjiang Province, China (No. 2019-KYYWFMY-0002) and the Mudanjiang Medical College Graduate Innovation Research Project (No. 2019YJSCX-10MY).

Author information




SH, ZS, and XM designed this study. ZW wrote the manuscript. ZW, DS, SL, GQ Z, GB Z, and TC collected the samples and performed data analyses. ZW and SL performed the experiments. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Xu Ma or Zhaoyi Sun or Shanshan Hu.

Ethics declarations

Ethics approval and consent to participate

The present study complied with the Declaration of Helsinki for research involving human subjects and was approved by the research ethics committee of Heilongjiang Provincial People’s Hospital. Written informed consent was obtained from each patient.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, Z., Su, D., Liu, S. et al. RNA sequencing and bioinformatics analysis of human lens epithelial cells in age-related cataract. BMC Ophthalmol 21, 152 (2021).

Download citation


  • Age-related cataract
  • Differentially expressed genes
  • Bioinformatics analysis
  • Pathogenic genes