Jiayue Zhou1,2†
Sifan Zheng
Xingchao Shentu- 1The Eye Center, Second Affiliated Hospital of School of Medicine, Zhejiang University, Hangzhou, China
- 2Zhejiang Provincial Key Lab of Ophthalmology, Second Affiliated Hospital of School of Medicine, Zhejiang University, Hangzhou, China
- 3GKT School of Medical Education, King’s College London, London, United Kingdom
- 4Institute of Translational Medicine, Zhejiang University School of Medicine, Hangzhou, China
Cold cataract is the reversible opacification of the lens when the temperature decreases. However, we observed that when temperature of the rats’ lens was maintained at a lower temperature for a prolonged time, the opacification of lens was only partly reversible. To review the potential molecular mechanism of the irreversible part of opacification under cold stimulation, we applied comparative transcriptomic and proteomic analysis to systematically investigate the molecular changes that occurred in the lens capsules of rats under low temperature treatments. The RNA sequencing based transcriptomic analysis showed a significant up-regulation of genes related to the lens structure and development in the Hypothermia Group. Hub genes were small heat shock proteins (sHSPs). Besides the same findings as the transcriptomic results, the liquid chromatography-tandem mass spectrometry based proteomic analysis also revealed the up-regulation of the apoptotic process. To further analyze the regulatory mechanism in this process, we subsequently performed integrated analysis and identified the down-regulation of Notch3/Hes1 and PI3K/Akt/Xiap signaling axis. Our research revealed the activation of the apoptotic process in rats’ lens under cold stimulation, and the sHSP related heat shock response as a potential protective factor through our transcriptomic and proteomic data.
Introduction
Cold cataract is a phenomenon described as a reversible opacification of the lens when the temperature decreases (Zigman and Lerman, 1964), and has been observed in the lenses of young rats, mice, fishes, cows, and humans (Zigman and Lerman, 1964; Benedek et al., 1979; Loewenstein and Bettelheim, 1979). Early studies revealed that the reversible part of the opacification was mainly caused by the temporary aggregation of γ crystalline (Zigman and Lerman, 1964; Siezen et al., 1985) or liquid-liquid phase separation (Tanaka et al., 1977). The severity of cold cataract was positively correlated with time and negatively correlated with temperature (Bermudez et al., 2011). However, we observed that when temperature of the rats’ lens was maintained at a lower temperature (4°C) for a prolonged time (8 h), the transparency of lens did not fully recover, and the tension of lens capsules increased. There lacks mechanistic research regarding the formation of these irreversible opacifications caused by hypothermal stimulation.
The lens is one of the most important refractive mediums in the eye. Keeping the lens transparent is important for vision and for accurate observation of the inside of the eye. The lens has no blood vessels or nerves, and the lens fiber cells in the core of the lens have degraded their organelles, including the nuclei (Bassnett et al., 2011). Lens epithelial cells (LECs) located under the anterior lens capsule are the most metabolically active part of the lens. LECs act as a barrier between the aqueous humor and send signals to the inside of the lens. LECs play a major role in maintaining lenses’ life activities and transparency. Abnormal changes and apoptosis of LECs are often observed across different types of cataract.
To identify the potential molecular mechanism of irreversible opacifications under cold treatment, we used RNA-seq and TMT based LC-MS/MS to systematically investigate the transcriptomic and proteomic changes in lens capsules of rats under low temperature treatments. To further analyze the regulatory mechanism in this process, we subsequently performed integrated analysis of the transcriptomic and proteomic changes.
Materials and Methods
Animal and Tissue Collection
All animal experiments complied with the Association for Research in Vision and Ophthalmology Statement for the Use of Animals in Ophthalmic and Vision Research. All animals were treated in accordance with the Zhejiang University Administration on Laboratory Animal Care. Male Sprague-Dawley rats (6 weeks, 200 g) were obtained from Shanghai SLAC Laboratory Animal Co., Ltd.
Rats were divided randomly into the Normal Group and the Hypothermia Group and subsequently euthanized. The lens capsules of rats in the Normal Group were collected directly after euthanasia. The lens of rats in Hypothermia Group were isolated and incubated in Hibernate-A medium (Thermo Fisher Scientific, MA, United States) at 4°C for 8 h. The lens capsules were separated from the lens cortex after being washed in 4°C phosphate-buffered solution (PBS). For RNA-seq and TMT, 5 rats were used from each group. Two lens capsules from a rat were combined as a sample (Figure 1A). The collected lens capsules were rapidly frozen in liquid nitrogen and stored at −80°C for further transcriptomic and proteomic analysis. For RT-qPCR, 6 rats were used from each group. 4 lens capsules from 2 rats were combined as a sample. The samples were put into a TRIzol reagent (Invitrogen) immediately after separation and stored at 4°C.
Figure 1. (A) Grouping of the rats and treatments of their lens. Strategy for transcriptomic and proteomic analysis of the lens capsules in each group. (B) The morphological changes induced by hypothermia and rewarming treatments.
RNA Sequencing
Sample Preparation
Total RNA of the Hypothermia Group (n = 5) and the Nomal Control Group (n = 5) were extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, United States). Agarose gel electrophoresis, NanoPhotometer spectrophotometer (IMPLEN, CA, United States) and RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, United States) were used for total RNA qualification. RNA samples with RIN ≥ 5.8 were select for sequencing.
Library Preparation
A total amount of 0.8 μg RNA per sample was used as input material for the library preparation. Sequencing libraries were generated using NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, United States) following manufacturer’s recommendations.
Sequencing
The library preparations were sequenced using the Novaseq 6000 platform (Illumina, United States) and 150 bp paired-end reads were generated.
Data Analysis
Clean data were obtained by removing reads containing adapter, reads containing ploy-N and low-quality reads from the raw data of fastq format. Q20, Q30 and GC content of the clean data were calculated at the same time. Paired-end clean reads were aligned to the rat reference genome Rnor_6.0 using Hisat2 (v2.0.5), featurecounts (v1.5.0-p3) was used to count the number of reads mapped to each gene. And then FPKM (expected number of Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced) of each gene was calculated based on the length of the gene and reads count mapped to this gene. Differential expression analysis of two groups was performed using DESeq2 v1.16.1. P-values were adjusted using Benjamini and Hochberg’s approach. Genes with | log2FoldChange | > 1 compared with Normal Control group and adjusted P < 0.05 were considered as DEGs for transcriptomic analysis.
Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS) Proteomic
Sample preparation: Samples were lysed and digested using pressure cycling technology (PCT), as described previously with modifications (Zhu and Guo, 2018). The lysis buffer contained 6 M urea (Sigma) and 2 M thiourea (Sigma) in 100 mM triethylammonium bicarbonate (TEAB, Thermo Fisher Scientific). The peptide concentration was measured using a ScanDrop spectrophotometer (Endress + Hauser Group, Lorrach, Germany).
Tandem Mass Tag (TMT) Labeling and Off-Line High pH Reversed-Phase Fractionation: Peptides were labeled with TMT 10-plex label reagents (Thermo Fisher Scientific, San Jose, United States) as described previously (Gao et al., 2020; Shen et al., 2020). Peptides were separated into 120 fractions, which were consolidated into 30 fractions. The fractions were subsequently dried and redissolved in 2% ACN/0.1% formic acid (FA).
DDA-MS Acquisition: The redissolved peptides were analyzed by LC-MS/MS with the Ultimate 3000 nanoLC-MS/MS system (Dionex LC-Packings, CA) coupled with a Q Exactive HF-X mass spectrometer (QE-HFX, Thermo Fisher Scientific) in data dependent acquisition (DDA) mode as described previously (Shen et al., 2020).
Data Analysis: The resultant MS data were analyzed with Proteome Discoverer (Version 2.4.1.15, Thermo Fisher Scientific) using a protein database composed of the Rat fasta database downloaded from UniProtKB on 12 May 2020, containing 9949 reviewed protein sequences. Enzyme, static modifications, precursor ion mass toleranceas were set as described previously (Shen et al., 2020). The peptide-spectrum-match allowed 1% target false discovery rate (FDR) (strict) and 5% target FDR (relaxed). Normalization was performed against the total peptide amount. Default setup was used for all other parameters. We used two tailed two-sided t-tests to identify the differential expression proteins between two groups. Differential expression proteins were defined by the thresholds based on adjusted P < 0.05, | log2FoldChange | > 0.25.
Integrated Analysis of Transcriptomic and Proteomic Results
Overlapped genes in the lists of DEGs and DEPs were taken to perform correlation analysis according to their log2FoldChanges using Prism v8.2.1 and PPI analysis using String v11.0. List of rat translation factors (TFs) and cofactors (co-TFs) was download from AnimalTFDB3.0 to identify the TFs or co-TFs in the overlapped genes. Then, KEGG analysis of these TFs and co-TFs was performed.
Real-Time Quantitative PCR
Total RNAs of the lens capsule samples of Hypothermia Group (n = 3) and the Nomal Control Group (n = 3) were extracted using Trizol Reagent (Thermo Fisher), four corneas served as one biological replicate. RNA concentrations and OD260/OD280 ratios (1.9–2.0) were determined using a Nanodrop 1,000 UV spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States). Approximately 1,500 ng RNA was used per reverse transcription reaction (PrimeScriptRT Reagent Kit; TaKaRa, Shiga, Japan). The TB Green Premix Ex Taq Kit (TaKaRa) was used for quantitative real-time PCR on a 7500 Fast Real-Time PCR System (ABI, Loma Linda, CA, United States), with a 2-step method. The expression levels of the target genes were normalized to that of β-actin and were calculated using the 2–ΔΔCT method and indicated as means ± SE, statistical analysis was conducted using t-test. The data are presented as fold-changes. The primers (TSINGKE, Beijing, China) are shown in Supplementary Table 1.
Data Analysis and Visualization
Gene Ontology (GO) and KEGG pathway enrichment analysis of differentially expressed genes (DEGs) implemented by The Database for Annotation, Visualization and Integrated Discovery (DAVID v6.8). The cut off value was FDR < 0.01 and results were visualized by R software (v3.6.1). Protein-protein interaction (PPI) analysis was implemented by String (v11.0), and visualized by Cytoscape (v3.7.2). The node degree were calculated using cytoHubba plugin in Cytoscape(v3.7.2)., and the gene with node degree ≥ 10 was set as the hub gene in PPI network.
Results
Irreversible Opacification in the Lens of Sprague-Dawley Rats After 8 h of Cold Treatment
The newly extracted lens of rats at normal body temperature (37°C) was completely transparent. After 8 h of low temperature treatment at 4°C, the lens showed severe opacification. Following rewarming to 37°C for 10 min, the cold cataract disappeared. However, compared with the lens before cold treatment, the whole lens still had irreversible opacity (Figure 1B).
Transcriptomic Analysis Reveals Up-Regulation of Lens-Specific Gene Expression and Small Heat Shock Proteins (sHSPs) in Cold Treated Lens Capsules
In order to identify the differential expressed mRNA under cold treatment, we compared the transcriptomic changes of rats’ lens capsules between the Hypothermia Group and the Normal Group. Approximately 45 million paired-end clean reads per sample on average were obtained, and mapped to a total number of 17,005 expressed genes.
DEGs analysis identified 627 significant changed genes between the two groups (adjusted P < 0.05, | log2FoldChange | > 1), of which 420 genes were up-regulated and 207 genes were down-regulated (Figures 2A,B).
Figure 2. Analysis of transcriptomic changes in cold treated lens capsules. (A) Volcano plots depict the log2FoldChange and statistical significance (|log2FoldChange| > 1), with the adjusted P-value expressed as −log10(adjusted P-value), the P-values were adjusted using the Benjamini and Hochberg method. Red points, green points, and yellow points represent up-regulated genes, down-regulated genes, and non-differentially expressed genes, respectively. (B) Heatmap of differentially expressed genes (DEGs) in the 10 samples. (C) Gene ontology (GO) analysis bar plot of DEGs. The top five GO terms with the lowest FDR values within each subcategory (BP, CC, MF) are listed. The number of DEGs annotated to the GO terms are shown on the bars. (D) Protein-protein interaction (PPI) network of DEGs annotated to the GO terms related to eyes. (E) Comparison of gene expression level of the hub genes in the PPI network between the two groups, adjusted RNA-seq means FPKM.
To find out the main functions of the DEGs and the biological processes they are involved in, we performed a GO enrichment analysis. The GO enrichment results suggested that many significantly enriched terms of the DEGs were related to the eye and the lens, such as lens development in camera type eye (GO:0002088, FDR = 2.13E-07), visual perception (GO:0007601, FDR = 1.91E-05), structural constituent of eye lens (GO:0005212, FDR = 1.48E-18) (Figure 2C). We subsequently performed a PPI analysis of DEGs enriched in those GO terms related to the eye (Figure 2D). We found that the gene cluster consisted of genes specifically expressed in the lens, including α, β, γ crystallin families, membrane protein Mip and Lim2, cytoskeleton protein Bfsp1, and Bfsp2. All these genes were up-regulated (Figure 2E). Surprisingly, among the hub genes, Cryaa and Cryab are sHSPs and Hsf4 is a heat shock transcriptome factor.
GO enrichment analysis of up-regulated and down-regulated DEGs separately shows interesting results. Up-regulated DEGs were enriched in the eye-related and cytoskeleton-related categories (Supplementary Table 2). Down-regulated DEGs were enriched in more comprehensive categories including cell adhesion, extracellular matrix, ion transmembrane transport, and response to stimuli such as drugs, starvation, wounds, and TGF beta (Supplementary Table 3).
Proteomic Changes Correspond With Transcriptomic Changes and Show the Up-Regulation of the Apoptotic Process
A total of 4,149 proteins were identified in the 10 samples. A total of 1,641 of these proteins were significant differentially expressed proteins (DEPs) (adjusted P < 0.05, | log2FoldChange | > 0.25). Among these DEPs, 306 were up-regulated and 1,335 were down-regulated (Figures 3A,B).
Figure 3. Analysis of proteomic changes in cold treated lens capsules. (A) Volcano plots depict the log2FoldChange and statistical significance, with the adjusted P-value expressed as −log10(adjusted P-value). Red points, green points, and yellow points represent up-regulated genes, down-regulated genes, and non-differentially expressed genes, respectively. (B) Heatmap of differentially expressed proteins (DEPs) in the 10 samples. (C) Biological process in Gene ontology (GO) analysis for the up- and down-regulated DEPs. Terms in orange bars were enriched from up-regulated proteins and terms in blue bars were enriched from down-regulated proteins.
Using the cut-off value of the adjusted P-value (Bonferroni) < 0.05, up-regulated DEPs were enriched in 3 significant biological processes Gene Ontology (GO) terms, while down-regulated DEPs were enriched in 25 terms. The results suggested the up-regulated DEPs were enriched in “lens development in camera-type eye” (GO:0002088, Padj = 0.039) in the category of “Biological Process.” The associated DEPs are Cryaa, Cryab, Cryba1, Crygb, Crygc, Crygs, and Mip, which are the same as the results of RNA-seq. “Structural constituent of eye lens” (GO:0005212, Padj = 7.21E-16) in the category of “Molecular Function” contains even more lens specific DEPs including Bfsp1 and Bfsp2 (Supplementary Table 3). Conversely, GO terms “negative regulation of apoptotic process” (GO:0043066, Padj = 1.77E-04), “cell division” (GO:0051301, Padj = 0.014), “cell cycle” (GO:0007049, Padj = 0.020) in which down-regulated DEPs are enriched, indicates the up-regulation of apoptosis and the stagnation of cell cycle and division (Figure 3C).
Integrated Analysis of Transcriptomic and Proteomic Results Reveals the Up-Regulation of Apoptosis
Protein expression levels are usually determined by the corresponding mRNA expression levels. To verify this regulatory relationship, we performed an integrated analysis of our transcriptomic and proteomic results. There were 421 overlapped genes belonging to both DEGs and DEPs in our study (adjusted P < 0.05, | log2FoldChange | > 0.25). Correlation analysis shows that the expression changing of protein and mRNA were positively correlated in general (r = 0.627, p < 0.0001; Figure 4A). Therefore, transcriptional regulation plays a major role when the lens is stimulated by hypothermia. We also compared the genes regulated in same direction in RNA and proteins with a database of translation factors (TFs) and cofactors (co-TFs) (Hu et al., 2019). 33 TFs or co-TFs were identified, among which 11 were up-regulated, and 22 were down-regulated (Figure 4B). To identify the key TFs and pathways, we performed a KEGG pathway enrichment analysis of these 33 genes The Notch signaling pathway (rno04330, including gene Notch3 and Hes1) and the Apoptosis pathway (rno04210, including gene Pik3r1 and Xiap) were identified (Figure 4C). We found that Notch3, Hes1, Pik3r1 and Xiap were all down-regulated DEGs and DEPs (Figure 4D). PPI of these four genes and their first neighbors in the PPI network of DEGs and DEPs regulated in the same direction shows a tight interaction (Figure 4E).
Figure 4. Integrated analysis of transcriptomic and proteomic changes in cold treated lens capsules of the Hypothermia Group. (A) Correlationship between transcriptomic and proteomic results, the horizontal coordinate indicates the log2FoldChange of DEGs, while the vertical coordinate indicatesthe log2FoldChange of DEPs. (B) Overlap of differentially expressed genes (DEGs, differentially expressed proteins (DEPs, and translation factors (TFs) and cofactors (co-TFs) database, up-regulated and down-regulated part were separated. (C) KEGG pathway enrichment analysis of differentially expressed TFs&co-TFs. (D) RNA and protein expression level of the four TFs&co-TFs in the Notch signaling pathway and Apoptosis KEGG pathways. (E) PPI network of Notch3, Hes1, Pik3r1, and Xiap and their first neighbors in the PPI network of DEGs and DEPs regulated in the same direction.
Confirmation of Gene Expression by RT-qPCR
We used RT-qPCR to confirm the results above. sHSPs Cryaa, Cryab, heat shock transcription factor Hsf4, transcription factor Hes1, transcription co-factors Notch3, Pik3r1, Xiap were selected. The RT-qPCR confirmed the expression change of these genes (Figure 5).
Figure 5. RT-qPCR confirmation of the expression level of Cryaa, Cryab, Hsf4, Notch3, Hes1, Pik3r1, and Xiap.
Discussion
We studied the transcriptomic and proteomic changes in cold-treated lens capsules of rats (Hypothermia Group) compared with the normal body temperature group (Normal Group. The transcriptomic and proteomic results showed high consistency as genes related to lens structure and development (α, β, γ crystallins, cytoskeletal protein BFSP1, and BFSP2, membrane proteins such as Mip and Lim2, and heat shock transcription factor Hsf4) were all significantly up-regulated. Based on the positive correlation between the transcriptomic and proteomic results, we hypothesized that transcriptional regulation played an important role and identified two apoptotic-related signaling pathways centered on four key regulatory factors. RT-qPCR was used to confirmed the gene expression level.
Previous studies have reported that mutations of these up-regulated genes in the lens were related to the formation of cataracts (Shiels and Hejtmancik, 2017). In addition to the membrane proteins Mip and Lim2 causing cataract formation by affecting cell permeability and adhesion (Louis et al., 1989; Ponnam et al., 2008; Kumari et al., 2013; Shentu et al., 2015), Cryaa, Cryab and Hsf4 are surprisingly all involved in the heat shock response. The α-crystallin subunits Cryaa and Cryab, also known as Hspb4 and Hspb5, are major protein components of the eye lens, and also belong to the sHSP (also called heat shock protein family B) family (Horwitz, 2003; Clark et al., 2012). sHSPs are a family of ATP-independent chaperones that can delay the onset of irreversible protein aggregation in response to cellular stressors (Janowska et al., 2019; Rosenzweig et al., 2019). α-crystallin appeared to have therapeutic effects in certain diseases. Injection of α-crystallin B prior to the intraocular pressure elevation markedly mitigated the neurodegenerative processes in rat glaucoma models (Anders et al., 2017). α-crystallin also protected axons from optic nerve degeneration in rats (Ying et al., 2008). In addition, we further screened our data and found that another two sHSPs, Hspb1 and Hspb2, were significantly up-regulated. Hspb1 could inhibit reactive oxygen species (ROS) and has a protection effect against atherosclerosis and coronary heart disease (Zhang et al., 2019). Hspb2 suppresses caspases-8 and 10 activation, thereby blocking downstream apoptotic events and inhibits apoptosis (Oshita et al., 2010). Upstream of these sHSPs is the heat shock transcription factor Hsf4, which is an important transcription factor in lens development and lens cell differentiation. Its target genes include the subunit of α, β, and γ crystallin, Bfsp1, Bfsp2 as well (Fujimoto et al., 2004; Shi et al., 2009). Hsf4 could protected human lens epithelial cells from oxidant stress by increasing the HMOX-1 expression in human lymphatic endothelial cells, which would then up-regulate antioxidant enzyme activity, and thus inhibit caspase family dependent apoptosis (Liao et al., 2018). In response to stress, the heat shock transcription factors can up regulate the cellular concentrations of some sHSPs (de Thonel et al., 2012; Zhong et al., 2016). sHSP Hspb1 and Cryab are target genes of Hsf4 (Fujimoto et al., 2004), which could be activated in response to cellular damage. Therefore, we postulate that the energy-saving heat shock response of sHSPs plays a major protective role in the homeostasis regulation of cold induced irreversible cataract formation.
In the integrated analysis of transcriptomic and proteomic results, we aimed to identify regulatory factors in the lens under cold stimulation, considering the positive correlation between the changing expression level of protein and mRNA. We discovered two important regulatory pathways, the anti-apoptotic pathway centered on Pik3r1 and Xiap, and the Notch signaling pathway centered on Notch3 and Hes1. Pik3r1 encodes the 85 kD regulatory subunit of Phosphoinositide-3-kinase (PI3K), which inhibits cell apoptosis through PI3K/Akt signaling pathway by activating NF-κB and inhibiting FOXO1 and BAD (Vivanco and Sawyers, 2002). The activation of PI3K/Akt signaling pathway was proven to inhibit oxidative stress and epithelial cell apoptosis in the lens of ultraviolet radiation-induced cataract rats (Yao and Yan, 2020). X-chromosome-linked inhibitor of apoptosis protein (Xiap) is the downstream target of Akt, which could be activated when Akt was phosphorylated (Castells et al., 2013). Xiap able to bind to and inhibit the caspases that mediate apoptotic cell death (Jost and Vucic, 2020). Akt1 and Akt2 were both detected as down-regulated in our proteomic data. Notch3 is a transmembrane, developmental signaling receptor, playing many crucial roles in developmental patterning, cell fate decisions, regulation of cell survival and proliferation (Hosseini-Alghaderi and Baron, 2020). Hes1 is one of the most studied targets of the Notch signaling pathway (Rani et al., 2016). Loss of Hes1 sensitizes mammalian cells to ER stress. Severe ER stress triggers cell death through expression of various pro-apoptotic genes (Lee et al., 2018). In a study of asthma pathogenesis, inhibition of Notch/Hes1 promotes PTEN expression (Li et al., 2020), which is an antagonist of PI3K/Akt pathway. On the contrary, promotion of the Notch3 production activated the downstream PI3K/Akt pathway and inhibited cardiomyocyte apoptosis, which alleviated cardiac ischemia/reperfusion injury (Zhang et al., 2015). Our results are consistent with those of past research, cold-induced apoptosis was observed in the myocardium (Kong et al., 2020) and rat hepatocytes (Vairetti et al., 2001). Apoptosis of lens epithelial cells was known to be an early event in the development of many types of cataracts (Li and Spector, 1996; Shiels and Hejtmancik, 2019). Previous studies suggested that the dominant-negative mutation of Pik3r1 could lead to cataract in human and mice (Solheim et al., 2017), and the activation of PI3K/AKT signaling pathway by circRNA HIPK3 or melatonin could inhibit apoptosis of HLECs (Bai et al., 2013; Cui et al., 2020). Notch2 was proved to be inhibited in the formation of age-related cataract (Fan et al., 2017). Hes1 was detected down-regulated in the cell model of glucocorticoid cataract (Gupta et al., 2005). Thus, the cold-induced irreversible cataract may be the result of lens epithelial cells apoptosis caused by the down regulation of Notch3/Hes1 and PI3K/Akt/Xiap signaling axis. Our study was a further proof of the importance of these two pathways in the formation of cataract. The down-regulated regulatory molecules Pik3r1, Xiap, Notch3, and Hes1 we’ve found out from our dataset need further research for their potential value.
In conclusion, the down-regulation of key molecules in the anti-apoptotic pathway and the Notch signaling pathway promotes the apoptosis of lens epithelial cells and promotes the occurrence of opacity. On the other hand, lens epithelial cells respond to cold stimulation through heat shock response. The expression of heat shock transcription factor Hsf4 was up-regulated, thereby enhancing the expression of sHSP family members.
We hypothesize that the balance was broken under cold stimulation. The up-regulated heat shock response was not strong enough to inhibit apoptosis successfully in physiological state, causing the lens to form irreversible opacity. There were some limitations of our study. First, the Normal Control Group was not treated with hibernation medium, but considering the RNA degradation of in vitro tissues at body temperature, we chose to collect the sample directly after euthanasia. Second, for some of the antibodies were difficult to obtain, we only verified the results by RT-qPCR. If possible, we would collect and sequence cold-treated lens capsule samples at more time points and different temperature in the future to obtain and analyze a series of constant data. Further research on the potential protective effects of the sHSPs and the anti-apoptotic Notch3/Hes1 and PI3K/Akt/Xiap signaling pathways in lens epithelial cells should also be considered.
Data Availability Statement
The transcriptomic data presented in the study are deposited in the GEO repository, accession number GSE173481. The proteomic data presented in the study are deposited in the iProX repository, accession number IPX0002974000.
Ethics Statement
The animal study was reviewed and approved by the Zhejiang University Administration on Laboratory Animal Care.
Author Contributions
XS designed the study and revised the manuscript. JZ and JW designed the study, collected the samples, and drafted the manuscript. JZ and DZ analyzed the data. DZ, SZ, and XC revised the manuscript. All authors read and approved the final manuscript.
Funding
This study was supported by Grant 81970781 to XS from the National Natural Science Foundation of China; Grant LD21H120001 to XS from the Natural Science Foundation of Zhejiang Province.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Xiao Yi and Tiannan Guo from Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University for technical support of mass spectrometric analysis. We thank Novogene for assistance in RNA-seq experiments.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.683056/full#supplementary-material
References
Anders, F., Liu, A., Mann, C., Teister, J., Lauzi, J., Thanos, S., et al. (2017). The small heat shock protein α-Crystallin B shows neuroprotective properties in a glaucoma animal model. Int. J. Mol. Sci. 18:2418. doi: 10.3390/ijms18112418
Bai, J., Dong, L., Song, Z., Ge, H., Cai, X., Wang, G., et al. (2013). The role of melatonin as an antioxidant in human lens epithelial cells. Free Radical Res. 47, 635–642. doi: 10.3109/10715762.2013.808743
Bassnett, S., Shi, Y., and Vrensen, G. F. J. M. (2011). Biological glass: structural determinants of eye lens transparency. Philos. Trans. R. Soc. London. Ser.B Biol. Sci. 366, 1250–1264. doi: 10.1098/rstb.2010.0302
Benedek, G. B., Clark, J. I., Serrallach, E. N., Young, C. Y., Mengel, L., Sauke, T., et al. (1979). Light scattering and reversible cataracts in the calf and human lenS. Philos. Trans. R. Soc. A Math. Phys. Engin. Sci. 293:329. doi: 10.1098/rsta.1979.0100
Bermudez, M. A., Vicente, A. F., Romero, M. C., Arcos, M. D., Abalo, J. M., and Gonzalez, F. (2011). Time course of cold cataract development in anesthetized mice. Curr. Eye Res. 36, 278–284. doi: 10.3109/02713683.2010.542868
Castells, M., Milhas, D., Gandy, C., Thibault, B., Rafii, A., Delord, J. P., et al. (2013). Microenvironment mesenchymal cells protect ovarian cancer cell lines from apoptosis by inhibiting XIAP inactivation. Cell Death Dis. 4:e887. doi: 10.1038/cddis.2013.384
Clark, A. R., Lubsen, N. H., and Slingsby, C. (2012). sHSP in the eye lens: crystallin mutations, cataract and proteostasis. Int. J. Biochem. Cell Biol. 44, 1687–1697. doi: 10.1016/j.biocel.2012.02.015
Cui, G., Wang, L., and Huang, W. (2020). Circular RNA HIPK3 regulates human lens epithelial cell dysfunction by targeting the miR-221-3p/PI3K/AKT pathway in age-related cataract. Exp. Eye Res. 198:108128. doi: 10.1016/j.exer.2020.108128
de Thonel, A., Le Mouël, A., and Mezger, V. (2012). Transcriptional regulation of small HSP-HSF1 and beyond. Int. J. Biochem. Cell Biol. 44, 1593–1612. doi: 10.1016/j.biocel.2012.06.012
Fan, F., Zhuang, J., Zhou, P., Liu, X., and Luo, Y. (2017). MicroRNA-34a promotes mitochondrial dysfunction-induced apoptosis in human lens epithelial cells by targeting Notch2. Oncotarget 8, 110209–110220. doi: 10.18632/oncotarget.22597
Fujimoto, M., Izu, H., Seki, K., Fukuda, K., Nishida, T., Yamada, S.-I., et al. (2004). HSF4 is required for normal cell growth and differentiation during mouse lens development. EMBO J. 23, 4297–4306. doi: 10.1038/sj.emboj.7600435
Gao, H., Zhang, F., Liang, S., Zhang, Q., Lyu, M., Qian, L., et al. (2020). Accelerated lysis and proteolytic digestion of biopsy-level fresh-frozen and FFPE tissue samples using pressure cycling technology. J. Proteome Res. 19, 1982–1990. doi: 10.1021/acs.jproteome.9b00790
Gupta, V., Galante, A., Soteropoulos, P., Guo, S., and Wagner, B. J. (2005). Global gene profiling reveals novel glucocorticoid induced changes in gene expression of human lens epithelial cells. Mol. Vis. 11, 1018–1040.
Hosseini-Alghaderi, S., and Baron, M. (2020). Notch3 in development, health and disease. Biomolecules 10:485. doi: 10.3390/biom10030485
Hu, H., Miao, Y.-R., Jia, L.-H., Yu, Q.-Y., Zhang, Q., and Guo, A.-Y. (2019). AnimalTFDB 3.0: a comprehensive resource for annotation and prediction of animal transcription factors. Nucleic Acids Res. 47, D33–D38. doi: 10.1093/nar/gky822
Janowska, M. K., Baughman, H. E. R., Woods, C. N., and Klevit, R. E. (2019). Mechanisms of small heat shock proteins. Cold Spring Harbor Perspec. Biol. 11:a034025. doi: 10.1101/cshperspect.a034025
Jost, P. J., and Vucic, D. (2020). Regulation of cell death and immunity by XIAP. Cold Spring Harbor Perspec. Biol. 12:a036426. doi: 10.1101/cshperspect.a036426
Kong, X., Liu, H., He, X., Sun, Y., and Ge, W. (2020). Unraveling the mystery of cold stress-induced myocardial injury. Front. Physiol. 11:580811. doi: 10.3389/fphys.2020.580811
Kumari, S. S., Gandhi, J., Mustehsan, M. H., Eren, S., and Varadaraj, K. (2013). Functional characterization of an AQP0 missense mutation, R33C, that causes dominant congenital lens cataract, reveals impaired cell-to-cell adhesion. Exp. Eye Res. 116, 371–385. doi: 10.1016/j.exer.2013.09.019
Lee, J. E., Morrison, W., and Hollien, J. (2018). Hairy and enhancer of split 1 (HES1) protects cells from endoplasmic reticulum stress-induced apoptosis through repression of. J. Biol. Chem. 293, 5947–5955. doi: 10.1074/jbc.RA118.002124
Li, W. C., and Spector, A. (1996). Lens epithelial cell apoptosis is an early event in the development of UVB-induced cataract. Free Radical Biol. Med. 20, 301–311. doi: 10.1016/0891-5849(96)02050-3
Li, X., Zou, F., Lu, Y., Fan, X., Wu, Y., Feng, X., et al. (2020). Notch1 contributes to TNF-α-induced proliferation and migration of airway smooth muscle cells through regulation of the Hes1/PTEN axis. Int. Immunopharmacol. 88:106911. doi: 10.1016/j.intimp.2020.106911
Liao, S., Qu, Z., Li, L., Zhou, B., Gao, M., Huang, M., et al. (2018). HSF4 transcriptional regulates HMOX-1 expression in HLECs. Gene 655, 30–34. doi: 10.1016/j.gene.2018.02.033
Loewenstein, M. A., and Bettelheim, F. A. (1979). Cold cataract formation in fish lenses. Exp. Eye Res. 28, 651–663. doi: 10.1016/0014-4835(79)90066-6
Louis, C. F., Hur, K. C., Galvan, A. C., TenBroek, E. M., Jarvis, L. J., Eccleston, E. D., et al. (1989). Identification of an 18,000-dalton protein in mammalian lens fiber cell membranes. J. Biol. Chem. 264, 19967–19973. doi: 10.1016/s0021-9258(19)47205-0
Oshita, S. E., Chen, F., Kwan, T., Yehiely, F., and Cryns, V. L. (2010). The small heat shock protein HspB2 is a novel anti-apoptotic protein that inhibits apical caspase activation in the extrinsic apoptotic pathway. Breast Cancer Res. Treat. 124, 307–315. doi: 10.1007/s10549-010-0735-0
Ponnam, S. P. G., Ramesha, K., Tejwani, S., Matalia, J., and Kannabiran, C. (2008). A missense mutation in LIM2 causes autosomal recessive congenital cataract. Mol. Vis. 14, 1204–1208.
Rani, A., Greenlaw, R., Smith, R. A., and Galustian, C. (2016). HES1 in immunity and cancer. Cytokine Growth Fact. Rev. 30, 113–117. doi: 10.1016/j.cytogfr.2016.03.010
Rosenzweig, R., Nillegoda, N. B., Mayer, M. P., and Bukau, B. (2019). The Hsp70 chaperone network. Nat. Rev. Mol. Cell Biol. 20, 665–680. doi: 10.1038/s41580-019-0133-3
Shen, B., Yi, X., Sun, Y., Bi, X., Du, J., Zhang, C., et al. (2020). Proteomic and metabolomic characterization of COVID-19 patient sera. Cell 182, 59–72. doi: 10.1016/j.cell.2020.05.032
Shentu, X., Miao, Q., Tang, X., Yin, H., and Zhao, Y. (2015). Identification and functional analysis of a novel MIP gene mutation associated with congenital cataract in a chinese family. PloS One 10:e0126679. doi: 10.1371/journal.pone.0126679
Shi, X., Cui, B., Wang, Z., Weng, L., Xu, Z., Ma, J., et al. (2009). Removal of Hsf4 leads to cataract development in mice through down-regulation of gamma S-crystallin and Bfsp expression. BMC Mol. Biol. 10:10. doi: 10.1186/1471-2199-10-10
Shiels, A., and Hejtmancik, J. F. (2017). Mutations and mechanisms in congenital and age-related cataracts. Exp. Eye Res. 156, 95–102. doi: 10.1016/j.exer.2016.06.011
Shiels, A., and Hejtmancik, J. F. (2019). Biology of inherited cataracts and opportunities for treatment. Ann. Rev. Vis. Sci. 5, 123–149. doi: 10.1146/annurev-vision-091517-034346
Siezen, R. J., Fisch, M. R., Slingsby, C., and Benedek, G. B. (1985). Opacification of gamma-crystallin solutions from calf lens in relation to cold cataract formation. Proc. Natl. Acad. Sci. U S A. 82, 1701–1705. doi: 10.1073/pnas.82.6.1701
Solheim, M. H., Clermont, A. C., Winnay, J. N., Hallstensen, E., Molven, A., Njølstad, P. R., et al. (2017). Iris malformation and anterior segment dysgenesis in mice and humans with a mutation in PI 3-Kinase. Invest. Ophthalmol. Vis. Sci. 58, 3100–3106. doi: 10.1167/iovs.16-21347
Tanaka, T., Ishimoto, C., and Chylack, L. T. (1977). Phase separation of a protein-water mixture in cold cataract in the young rat lens. Science 197, 1010–1012. doi: 10.1126/science.887936
Vairetti, M., Griffini, P., Pietrocola, G., Richelmi, P., and Freitas, I. (2001). Cold-induced apoptosis in isolated rat hepatocytes: protective role of glutathione. Free Radical Biol. Med. 31, 954–961. doi: 10.1016/s0891-5849(01)00670-0
Vivanco, I., and Sawyers, C. L. (2002). The phosphatidylinositol 3-Kinase AKT pathway in human cancer. Nat. Rev. Cancer 2, 489–501. doi: 10.1038/nrc839
Yao, L., and Yan, H. (2020). MiR-182 inhibits oxidative stress and epithelial cell apoptosis in lens of cataract rats through PI3K/Akt signaling pathway. Eur. Rev. Med. Pharmacol. Sci. 24, 12001–12008. doi: 10.26355/eurrev_202012_23988
Ying, X., Zhang, J., Wang, Y., Wu, N., Wang, Y., and Yew, D. T. (2008). Alpha-crystallin protected axons from optic nerve degeneration after crushing in rats. J. Mol. Neurosci. 35, 253–258. doi: 10.1007/s12031-007-9010-1
Zhang, H.-L., Jia, K.-Y., Sun, D., and Yang, M. (2019). Protective effect of HSP27 in atherosclerosis and coronary heart disease by inhibiting reactive oxygen species. J. Cell. Biochem. 120, 2859–2868. doi: 10.1002/jcb.26575
Zhang, M., Wang, C., Hu, J., Lin, J., Zhao, Z., Shen, M., et al. (2015). Notch3/Akt signaling contributes to OSM-induced protection against cardiac ischemia/reperfusion injury. Apoptosis 20, 1150–1163. doi: 10.1007/s10495-015-1148-7
Zhong, Y.-H., Cheng, H.-Z., Peng, H., Tang, S.-C., and Wang, P. (2016). Heat shock factor 2 is associated with the occurrence of lung cancer by enhancing the expression of heat shock proteins. Oncol. Lett. 12, 5106–5112. doi: 10.3892/ol.2016.5368
Zhu, Y., and Guo, T. (2018). High-throughput proteomic analysis of fresh-frozen biopsy tissue samples using pressure cycling technology coupled with SWATH mass spectrometry. Methods Mol. Biol. 1788, 279–287. doi: 10.1007/7651_2017_87
Keywords: hypothermia, lens, transcriptomics, proteomics, small heat shock proteins, apoptosis
Citation: Zhou J, Wu J, Zheng S, Chen X, Zhou D and Shentu X (2021) Integrated Transcriptomic and Proteomic Analysis Reveals Up-Regulation of Apoptosis and Small Heat Shock Proteins in Lens of Rats Under Low Temperature. Front. Physiol. 12:683056. doi: 10.3389/fphys.2021.683056
Received: 22 March 2021; Accepted: 25 May 2021;
Published: 17 June 2021.
Edited by:
Nathaniel J. Szewczyk, Ohio University, United StatesReviewed by:
Colleen Siobhan Deane, University of Exeter, United KingdomJoseph J. Bass, University of Nottingham, United Kingdom
Copyright © 2021 Zhou, Wu, Zheng, Chen, Zhou and Shentu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Daizhan Zhou, 3416252@zju.edu.cn; Xingchao Shentu, stxc@zju.edu.cn
†These authors have contributed equally to this work and share first authorship