Galectin-1-dependent ceRNA network in HRMECs revealed its association with retinal neovascularization
BMC Genomics volume 24, Article number: 327 (2023)
Retinal neovascularization (RNV) is a leading cause of blindness worldwide. Long non-coding RNA (lncRNA) and competing endogenous RNA (ceRNA) regulatory networks play vital roles in angiogenesis. The RNA-binding protein galectin-1 (Gal-1) participates in pathological RNV in oxygen-induced retinopathy mouse models. However, the molecular associations between Gal-1 and lncRNAs remain unclear. Herein, we aimed to explore the potential mechanism of action of Gal-1 as an RNA-binding protein.
A comprehensive network of Gal-1, ceRNAs, and neovascularization-related genes was constructed based on transcriptome chip data and bioinformatics analysis of human retinal microvascular endothelial cells (HRMECs). We also conducted functional enrichment and pathway enrichment analyses. Fourteen lncRNAs, twenty-nine miRNAs, and eleven differentially expressed angiogenic genes were included in the Gal-1/ceRNA network. Additionally, the expression of six lncRNAs and eleven differentially expressed angiogenic genes were validated by qPCR in HRMECs with or without siLGALS1. Several hub genes, such as NRIR, ZFPM2-AS1, LINC0121, apelin, claudin-5, and C-X-C motif chemokine ligand 10, were found to potentially interact with Gal-1 via the ceRNA axis. Furthermore, Gal-1 may be involved in regulating biological processes related to chemotaxis, chemokine-mediated signaling, the immune response, and the inflammatory response.
The Gal-1/ceRNA axis identified in this study may play a vital role in RNV. This study provides a foundation for the continued exploration of therapeutic targets and biomarkers associated with RNV.
Retinal neovascular diseases, such as proliferative diabetic retinopathy, retinopathy of prematurity, and retinal vein occlusion, are the major causes of blindness . Retinal angiogenesis is the process by which new blood vessels sprout from preexisting retinal vessels and invade the retina . During normal retinal vascular development, vascular endothelial cells proliferate and migrate through the extracellular matrix under the action of various cytokines, forming new vessels in a highly ordered manner. The unique vascular supply of the retina as well as the ability to visualize its vasculature in vivo and selectively express its genes make the retina an ideal model system for studying the molecular mechanisms of angiogenesis. Retinal neovascularization (RNV), also known as abnormal angiogenesis, is a key pathological process in ocular blinding diseases . These newly formed blood vessels are fragile and readily lead to retinal edema, leakage, hemorrhage, and detachment, resulting in catastrophic visual impairment and blindness . Angiogenesis is a complex process that depends on precisely regulated gene expression. Several important factors related to angiogenesis have been identified, including vascular endothelial growth factor (VEGF), platelet-derived growth factor, and hypoxia-inducible factors [5,6,7]. However, the biochemical events responsible for RNV development and progression are not fully understood. Although anti-VEGF therapy is considered a principal clinical treatment regime, the associated therapeutic effects can be insufficient in some patients and are accompanied by a high recurrence rate after a single treatment as well as rapid allergic reactions . Therefore, it is necessary to develop other therapeutic targets and biomarkers. However, to achieve this, studies on the upstream and downstream regulatory mechanisms of these neovascularization-related factors are needed to elucidate additional potential targets for the treatment of ischemic retinopathy.
Recently, an increasing number of studies have found that RNA-binding proteins (RBPs) and non-coding RNAs play important regulatory roles in major physiological and pathological changes . RBPs can regulate basal RNA transcription and post-transcriptional processes, including variable cleavage, modification, transport, translation, degradation, and metabolism by binding to RNA. Abnormal expression, or loss of function, of RBPs leads to the occurrence and development of various diseases . In particular, RBPs have important functions in eye development and ocular diseases, making them potential therapeutic targets for neovascularization .
MicroRNAs (miRNAs) are typically ~ 22 nt in length and negatively regulate target messenger RNAs (mRNAs) by promoting mRNA degradation or inhibiting translation. Long non-coding RNAs (lncRNAs) are > 200 nt in length and exhibit no, or limited, protein-coding potential. Accumulating evidence indicates that lncRNAs play crucial regulatory roles in a wide range of biological and pathological processes including angiogenesis . For example, lncRNA-Meg3 knockout mice exhibit increased expression of VEGF pathway genes and increased cortical microvessel density . In fact, many abnormally expressed lncRNAs are associated with eye disorders . For example, MALAT1 up-regulation is an important pathogenic mechanism of diabetes-induced microvascular dysfunction. Hence, inhibition of MALAT1 may be a potential target for anti-angiogenic therapy of diabetes-related microvascular complications . Moreover, the lncRNA TUG1 reduces the expression of VEGFA in OIR mice by competitively adsorbing miR-299, thus participating in the regulation of RNV . Additionally, numerous studies have revealed that competing endogenous RNA (ceRNA) networks participate in gene regulation at multiple levels . lncRNAs are predominantly located in the nucleus, which underscores their regulatory role in gene transcription. Studies have indicated that lncRNAs can be targeted for therapeutic purposes. In particular, the role of ceRNA regulatory networks among lncRNAs, miRNAs, and mRNAs in RNV progression has drawn increasing attention. For instance, knockdown of the lncRNA MIAT considerably improves diabetic retinal microvascular dysfunction in vivo and inhibits endothelial cell proliferation, migration, and tube formation in vitro. MIAT functions as a competitive endogenous RNA and forms a feedback loop with vascular endothelial growth factor and Mir-150-5p to regulate endothelial cell function . Meanwhile, MALAT1 knockdown has an important regulatory role in diabetic retinopathy by inhibiting proliferation, migration, tubule formation, and vascular permeability of human retinal microvascular endothelial cells (HRMECs) induced by mercury through up-regulation of Mir-125b .
Galectins are expressed in a wide range of species, suggesting their important role in basic cellular mechanisms. Galectin-1 (Gal-1), encoded by LGALS1, is a carbohydrate-binding protein that shares a conserved carbohydrate-recognition domain . It has been proposed to mediate cell adhesion and migration and is involved in modulating various aspects of cellular biology, including cell cycle, proliferation, migration, apoptosis, angiogenesis, and cell adhesion . A study indicated that Gal-1 played a central role in maintaining angiogenesis in endometriosis and may be a potential therapeutic target for this disease . We previously demonstrated that siLGALS1 inhibits pathological RNV and promotes retinal vascular normalization in an oxygen-induced retinopathy mouse model, suggesting its crucial role in triggering vascular signaling programs and mediating the anti-angiogenic treatment response . Recently, several independent experiments showed that Gal-1 could achieve its biological functions, including angiogenesis via binding to lncRNAs and mRNAs, indicating its potential RNA-binding activity [24,25,26]. In addition, Gal-1 was found to interact with angiogenesis-related mRNA, including VEGFA, EGR1, and LAMA5. ShLGALS1 inhibits capillary formation and alters the expression levels of several Gal-1 binding angiogenesis-related mRNAs within in vitro angiogenesis assays . However, the potential downstream mechanisms of Gal-1 in the RNV process have not been extensively studied.
In this study, we explored the potential RNA-binding activity of Gal-1. RBPs are associated with several cellular functions that interact with various RNAs. For instance, RBP FUS can regulate glioma angiogenesis through the FUS/circ_002136/Mir-138-5p/SOX13 feedback loop . Additionally, lncRNA-422 reportedly inhibits the proliferation and growth of colorectal cancer cells by interacting with RBP SFPQ . Meanwhile, cumulative evidence has demonstrated that lncRNAs are involved in nearly all cellular processes. Given the complex function of RBP-lncRNA pairs, aberrantly expressed lncRNAs may be partially responsible for Gal-1-mediated RNV . Continued studies on the function of abnormally expressed RBPs in retinal neovascular diseases and their possible upstream and downstream mechanisms of action are required to develop personalized treatment regimens. Indeed, downstream molecules of Gal-1 may serve as an effective indicator of disease states and might represent a superior therapeutic target for existing protein-coding genes; however, further in-depth research is required to examine the toxicity and pharmacokinetics of lncRNAs.
Herein, we construct a comprehensive regulatory network of lncRNAs, miRNAs, and mRNAs in siLGALS-knockdown HRMECs. This study describes a previously unrecognized relationship between Gal-1 and the ceRNA regulatory network and provides a foundation for the development of a novel treatment for neovascularization-related diseases.
LGALS1 regulates the expression of genes associated with angiogenesis in HRMECs
We successfully constructed siLGALS1 HRMECs models for further investigation (Fig. 1). Si-1, si-2, si-3, si-4, and si-5 groups showed considerable down-regulation in LGALS1 mRNA, and si-3, si-4, and si-5 groups showed considerable down-regulation in Gal-1 protein levels. The original, unprocessed electrophoretic gel images are shown in the Supplementary file (Fig. S1). Taking all the results into account, we finally chose si-4 group to continue the experiments. We analyzed differentially-expressed genes (DEGs) between the LGALS1-silenced group (siLGALS1) and the control group (Ctrl) in HRMECs (Fig. 2a). Compared with those in the Ctrl, the downregulated genes in the siLGALS1 group were nearly thrice the number of upregulated genes (Fig. 2b). Upregulated genes were primarily concentrated in G protein-coupled receptor signaling pathway, cellular adhesion, and transcriptional regulation by RNA polymerase II pathway according to GO analysis, whereas downregulated genes were mainly concentrated in inflammatory response and chemokine-mediated signaling pathway (Fig. S2).
Here, we focused on DEGs associated with angiogenesis. Through DEGs and gene overlap contained in the gene ontology (GO) term of angiogenesis (GO:0,001,525), we found a considerable change in 14 angiogenesis-related genes, among which 10 were upregulated and 4 were downregulated (Fig. 2c). The 14 angiogenesis-related genes included EMILIN1, CLEC14A, APLN, and CLDN5 (downregulated) and AMOT, GNA13, C-X-C motif chemokine ligand 10 (CXCL10), ADAMTS1, IL1B, IL1A, CCL8, C3AR1, GPNMB, and DLL1 (upregulated).
Analysis of differentially expressed lncRNAs between LGALS1-knockdown and control HRMECs
LGALS1 may influence angiogenesis by regulating lncRNA expression [27, 31]. However, the Gal-1-lncRNA interaction map in angiogenesis remains unclear. LGALS1-associated lncRNAs were identified and analyzed in the current study.
The known and newly predicted lncRNAs were combined and a Venn diagram of the detected known lncRNAs and novel lncRNAs in siLGALS1 was constructed (Fig. 3a). LncRNAs with fragments per kilobase of exon per million (FPKM) ≥ 0.2 in at least one sample were considered to be detected in this group (Fig. S3). In total, 65 lncRNAs were significantly altered in HRMECs with LGALS1 knockdown compared with those in the control HRMECs (fold change [FC] ≥ 1.5 or ≤ 2/3, P ≤ 0.05). Among these, 21 upregulated and 44 downregulated annotated lncRNAs were identified through differential analysis in siLGALS1 cells (Fig. 3b). The heatmap shows the expression profiles of the top 20 most significant differentially expressed (DE) lncRNAs (Fig. 3c). Based on the number of co-expressed DE lncRNAs and DE mRNAs between siLGALS1 and control samples (Fig. 3d), the top 10 most enriched GO terms (biological processes) revealed that DE mRNAs may be closely associated with chemotaxis, immune responses, and inflammatory responses (Fig. 3e). The co-expression network between DE lncRNAs and DE mRNAs, which are involved in the 10 GO terms shown in red (Fig. S3g), was constructed. LINC01905, AJ006998.2, RP11-39,312.4, LINC02362, ZFPM2-AS1, RP11-726G1.2, RP11-322E11.5, D21S2088E, XLOC_OO2601, and RP11-74E22.8 were included. The top 10 most enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways by DE mRNAs co-expressed with DE lncRNAs between siLGALS1 and the control samples are demonstrated (Fig. S3f). The top three pathways were cytokine-cytokine receptor interaction, viral protein interaction with cytokine and cytokine receptor, and toll-like receptor signaling pathway.
Regulatory network of lncRNA-miRNA-angiogenesis genes mediated by LGALS1 deregulation in HRMEC
We constructed co-expression network of DE lncRNAs and DE angiogenesis-related genes (Fig. S4a). Eight DE lncRNAs (LINC02840, AJ006998.2, RP11-74E22.8, RP11-96A15.1, LUCAT1, LINC02362, LINC00592, and LINC01140) involved in the regulatory network were further verified (Fig. S4b). LncRNAs act as sponges of related miRNAs. Together, lncRNAs and miRNAs form an lncRNA–miRNA–angiogenesis-related gene axis mediated by LGALS1 deregulation in HRMECs. A regulatory network was established, comprising 14 lncRNAs, 29 miRNAs, and 11 angiogenic DEGs, and a network diagram was drawn using Cytoscape (Fig. 4a). The lncRNA-miRNA-mRNA pathways in this regulatory network diagram were selected for validation and follow-up studies. Expression profile supported six lncRNAs (WT1-AS1, ZFPM2-AS1, NRIR, XLOC_039889, THAP7-AS1, and LINC01215) as DE lncRNAs (Fig. 4b). Some DE lncRNAs were situated in the center of co-expression networks which more likely have an impact on angiogenesis regulation in siLGALS1 HRMECs. Among these angiogenesis-related genes, expression of APLN, CLDN5, and CLEC14A were up-regulated after siLGALS1 and they may participate in G protein-coupled receptor signaling pathway, cell adhesion, and neuroactive ligand-receptor interaction (Fig. S2a, c). Meanwhile, expression of AMOT, GNA13, IL1B, and CXCL10 were down-regulated in siLGALS1 cells and they may impact inflammatory response, chemokine-mediated signaling pathway, and viral protein interaction with cytokine and cytokine receptor (Fig. S2b, d).
Validation of the expression of lncRNA and angiogenesis-related genes in HRMECs with or without knockdown LGALS1
We further verified 6 crucial lncRNAs (down-regulated expression: NRIR, ZFPM2-AS1, and LINC0121; no significant differences: WT1-AS, THAP7-AS1, XLOC_039889) and 11 angiogenesis-related mRNAs (down-regulated expression: IL1A, GNA13, IL1B, ADAMTS1, CXCL10, C3AR1, DLL1; up-regulated expression: APLN, CLEC14A, CLDN5; no significant differences: AMOT) by qPCR in HRMECs with or without knockdown LGALS1 (Fig. 5 and Fig. S5).
Gal-1 is an endogenous lectin that exists in peripheral lymphoid organs and inflammatory sites. Several studies have demonstrated an intricate interplay between Gal-1 and lncRNAs. To illustrate, the lncRNA ESCCAL-1 may promote cell cycle progression by interacting with, and stabilizing, Gal-1 in esophageal squamous cell carcinoma . Moreover, the lncRNA small nucleolar RNA host gene 22 (SNHG22) impacts cancer initiation and progression via the miR-2467/Gal-1 signaling pathway . Nevertheless, the mechanism of action of most ncRNAs, as well as their interaction with RBPs, remain to be elucidated; in particular, the role of Gal-1-ncRNA interaction in the regulation of retinal neovascularization remains to be explored.
In the current study, we systemically profiled lncRNA expression and analyzed DE lncRNAs and neovascularization-related genes between LGALS1 knockdown and control groups in HEMECs using next-generation sequencing. We found a marked change in 14 angiogenesis-related genes, among which 10 were upregulated and 4 were downregulated (Fig. 2c). Apelin (APLN), angiomotin (AMOT), CXCL10, and claudin 5 (CLDN5) were located in the central nodes of the network (Fig. 4a). APLN is a recently discovered regulator of neovascularization that binds to the APJ receptor to promote proliferation, differentiation, and migration of vascular endothelial cells, while also forming lumen-like structures. APLN also plays an important role in retinal vascularization , while also impacting the progression of osteoarthritis by regulating VEGF-dependent angiogenesis and mir-150-5p expression . Meanwhile, deletion of AMOT can inhibit the migration and expansion of physiological and pathological vascular networks . In addition, the lncRNA UCA1 may directly regulate AMOT to activate target genes in epithelial ovarian cancer, thus highlighting a potential interaction between AMOT and lncRNAs . Additionally, CXCL10 plays a critical role in vascular remodeling following spinal cord injury, whereas angiogenesis is enhanced following anti-CXCL10 treatment of spinal cord injuries . CLDN5, a tight junction protein, is vital for maintaining integrity of endothelial cells and vasostabilizing . It also participates in regulating permeability of the blood-tumor barrier via lncRNA NEAT1/miR-181d-5p-mediated expression changes . We picked several genes to conduct further vertification. Interestingly, the expression levels of APLN and CLDN5 were elevated, while the expression level of CXCL10 was reduced. However, the expression level of AMOT showed no significant tendency. This possibly indicates that increased sample numbers are warranted to eliminate result deviation that may occur due to insufficient sample size.
The CeRNA network plays an important role in the regulation of retinal angiogenesis . In fact, the ceRNA hypothesis suggests a new mechanism of interaction among RNAs, in which lncRNAs block the binding of miRNAs to mRNAs by competitively sponging miRNAs, thereby regulating mRNA expression. Understanding this system of RNA communication will provide new insights into gene regulatory networks and mechanisms of human disease . We examined lncRNA-miRNA interactions to reveal the possible downstream regulatory mechanism of Gal-1 and constructed an lncRNA–angiogenesis-related mRNA co-expression network. Our results show that Gal-1 may be involved in regulating biological processes related to chemotaxis, chemokine-mediated signaling pathway, immune response, inflammatory response, positive regulation of cytosolic calcium ion concentration, cell–cell signaling, neutrophil chemotaxis, G protein-coupled receptor signaling pathway, and regulation of cell population proliferation (Fig. 3e). This is consistent with the reported function of Gal-1, that is, participating in inflammatory regulation and affecting both cell proliferation as well as chemotaxis [42, 43]. These pathways are closely related to the formation of new blood vessels. Gal-1 has emerged as a multifaceted protein that exerts a wide spectrum of regulatory effects both intracellularly and extracellularly . Whether Gal-1 can regulate RNV by acting either extracellularly through specific binding to cell surface glycan structures or intracellularly through modulation of pathways remains largely unexplored and is worth further exploration. Moreover, whether some agents, such as lactose, that inhibit the binding of galectins to the cell surface, or recombinant Gal-1 would affect the expression of these differentially expressed genes, when added to the cells, is worth further exploration. In the present study, we hypothesized that Gal-1 could bind RNA directly and focused on its potential intracellular regulatory functions first.
We also predicted key lncRNAs and miRNAs that may regulate the expression of angiogenesis-related genes. These key lncRNAs and miRNAs were then screened as potential prognostic biomarkers and therapeutic targets for angiogenesis-related diseases. We further verified the most likely crucial lncRNAs (WT1-AS, NRIR, ZFPM2-AS1, LINC01215, THAP7-AS1, and XLOC_039889). NRIR, ZFPM2-AS1, and LINC01215 expression levels were down-regulated in siLGALS1 HEMECs. ZFPM2-AS1 has been shown to play crucial roles in the progression of retinoblastoma, an angiogenesis-associated disease . The interferon (IFN) family of cytokines serves as pivotal regulators of angiogenesis , while NRIR was found to be a negative Regulator of IFN Response [47, 48]. Tumor growth of Clear Cell Renal Cell Carcinoma was suppressed by silenced LINC01215 via miR-184 . However, whether it can regulate RNV by ceRNA network remains unknown. Our study revealed lncRNA-miRNA-mRNA co-expression regulatory networks in siLGALS1 HRMECs, indicating the possible downstream modulation mechanism of Gal-1 in RNV. Whether the regulatory role of the ceRNA network itself is really valid when Gal-1 is normally expressed needs further exploration. Collectively, these findings suggest that the Gal-1/ceRNA pathway is a promising target for the diagnosis and treatment of neovascularization-associated diseases. More specifically, regulating the expression of both Gal-1 and downstream lncRNAs may maximize the efficiency of neovascularization inhibition.
In conclusion, our research provides insights regarding the possible Gal-1/ceRNA regulatory network in HRMECs with siLGALS1. However, additional studies elucidating the ceRNA regulatory networks requires validation. This lncRNA-miRNA-mRNA network may provide clues to the regulatory pathways in the RNV for the cascade-amplifying synergistic effects of Gal-1. Verification of this hypothesis is expected to further improve the understanding of the molecular mechanism underlying retinal neovascularization at the gene and molecular levels, thereby providing a new target for the prevention and treatment of retinal neovascularization diseases and providing a theoretical basis for the clinical inhibition of neovascularization and promotion of ischemic retinal reperfusion.
Cell culture and transfection
HRMECs lines (Angio-Proteomie, Boston, Mass, USA) were cultured at 37 °C with 5% CO2 in ECM supplemented with 10% fetal bovine serum (FBS), 100 µg/mL streptomycin, and 100 U/mL penicillin. To silence the expression of LGALS1 in HRMECs, an siRNA-containing plasmid was constructed using a vector. Transfected cells were harvested after 48 h for RT-qPCR. All siRNA duplexes were purchased from Gemma (Suzhou, China). Non-targeting control siRNA (siNegative): 5′-UUCUCCGAACGUGUCACGUTT-3′ (sense). siRNA targeting LGALS1 (siLGALS1): CCAGCAACCUGAAUCUCAATT (sense). siRNA transfection of HRMECs was performed using Lipofectamine RNAiMAX Transfection Reagent (Invitrogen, Carlsbad, CA, USA), according to the manufacturer’s protocol.
Assessment of gene expression
GAPDH (glyceraldehyde-3-phosphate dehydrogenase) was used as a control gene to assess the effects of LGALS1 knockdown. cDNA synthesis was performed using standard procedures, and RT-qPCR was performed on a Bio-Rad S1000 with Hieff™ qPCR SYBR® Green Master Mix (Low Rox Plus; YEASEN, Shanghai, China). Information on the primers is presented in Table S1. The concentration of each transcript was normalized to the GAPDH mRNA level using the 2−ΔΔCT method .
HRMECs were lysed in ice-cold wash buffer (1 × PBS, 0.1% SDS, 0.5% NP-40, and 0.5% sodium deoxycholate) supplemented with protease inhibitor cocktail (Roche) and incubated on ice for 30 min. Samples were boiled for 10 min in boiling water with 1X SDS sample buffer and separated using 10% SDS-PAGE. Membranes were incubated with TBST buffer (20 mM Tris-buffered saline and 0.1% Tween-20) containing 5% non-fat milk powder for 1 h at room temperature. Membranes were incubated with primary antibody: LGALS1 antibody (A18040, 1:1,000, ABclonal, Boston, MA, USA) and actin (10,298–1-AP, 1:2000, ABclonal, Boston, MA, USA), followed by incubation with HRP-conjugated secondary antibody. Bound secondary antibodies (anti-rabbit 1:10,000, Abcam, Cambridge, England) were detected using an enhanced chemiluminescence (ECL) reagent (Bio-Rad, CA, USA).
Total RNA was extracted using TRIZOL (Ambion, Austin, Texas, USA) and purified using two phenol–chloroform treatments followed by treated with RQ1 DNase (Promega, Madison, WI, USA) to remove DNA. The quality and quantity of the purified RNA were determined by measuring the absorbance at 260 nm/280 nm (A260/A280) using a Nanodrop spectrophotometer (Thermo, Waltham, MA, USA). RNA integrity was verified by 1.0% agarose gel electrophoresis.
Reads alignment and DEG analysis
Clean reads were aligned to the human GRch38 genome using HISAT2 . Uniquely mapped reads were used to calculate the read number and FPKM mapped for each gene. The expression levels of these genes were evaluated by FPKM. DEseq2 models were employed to detect DEGs. Briefly, the original reads were used and a scale factor was applied to explain the difference in library depth. DEseq2 then estimated gene dispersion and reduced the estimates to generate more accurate dispersion estimates to model the read count. Finally, the model of the negative binomial distribution was fitted using DEseq2, and the hypothesis was tested using the Wald test or likelihood ratio test. DEseq2 can be used to analyze the differential expression between two or more samples, and the analysis results can be used to determine whether a gene is differentially expressed by FC, absolute ratio of expression change, and p-value. The criteria for significant differential expression were: FC ≥ 1.5 or ≤ 0.66, p-value ≤ 0.05.
LncRNA prediction and direction identification
To systematically analyze the lncRNA expression pattern, we used a pipeline for lncRNA identification similar to that previously reported , which was constructed based on the StringTie software . The predicted lncRNAs were compared to the noncoding RNA database (NONCODEv6, http://www.noncode.org/) to identify the lncRNAs that were already present.
Construction of regulatory lncRNA-miRNA-angiogenesis-related gene networks
First, co-expression analysis was performed for 80 differentially expressed lncRNAs and 14 differentially expressed angiogenesis-related genes, requiring an absolute value of Pearson correlation coefficient ≥ 0.95 and p value ≤ 0.01. LncRNAs and angiogenesis-related genes were identified in this study. Then, according to the TargetScan (http://www.targetscan.org/vert_80/) and miRDB (http://mirdb.org/) databases, miRNAs that can target the angiogenesis genes were obtained. Transcription sequences of lncRNAs were extracted, and miRanda (score ≥ 150) and RNAhybrid (p-value ≤ 0.05) were used to predict the target of miRNAs. Finally, the results of the two methods were compared, and partial miRNAs and targeted lncRNAs were identified. A regulatory network was established for lncRNAs, miRNAs, and angiogenesis-related genes using Cytoscape (https://cytoscape.org/).
Gene Ontology and KEGG pathway enrichment analyses
To sort assess functional categories of DEGs, GO terms and KEGG pathways were identified using KOBAS 2.0 . The hypergeometric test and Benjamini–Hochberg FDR control procedure were used to define the enrichment of each term.
Other statistical analysis
The pheatmap package (https://cran.r-project.org/web/packages/pheatmap/index.html) in R was used to perform clustering based on Euclidean distance. The Student’s t-test was used for comparison between the two groups.
Availability of data and materials
The datasets generated and analyzed during the current study are available in the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/) database (Accession Number: GSE221037) (The following secure token has been created to allow review of record GSE221037: kxwfoiqcbngvnaf); TargetScan (http://www.targetscan.org/vert_80/); miRDB (http://mirdb.org/); Gene Ontology(http://geneontology.org/); Kyoto Encyclopedia of Genes and Genomes (http://www.kegg.jp/); NONCODE (http://www.noncode.org/).
Competing endogenous RNA
C-X-C motif chemokine ligand 10
Differentially expressed gene
Fetal bovine serum
Fragments per kilobase of exon per million
Kyoto Encyclopedia of Genes and Genomes
Long non-coding RNA
Pearson’s correlation coefficient
Small nucleolar RNA host gene 22
Vascular endothelial growth factor
Wong TY, Cheung CM, Larsen M, Sharma S, Simó R. Diabetic retinopathy. Nat Rev Dis Primers. 2016;2:16012.
Carmeliet P, Jain RK. Molecular mechanisms and clinical applications of angiogenesis. Nature. 2011;473:298–307.
Dorrell M, Uusitalo-Jarvinen H, Aguilar E, Friedlander M. Ocular neovascularization: basic mechanisms and therapeutic advances. Surv Ophthalmol. 2007;52(Suppl 1):S3-19.
Gariano RF, Gardner TW. Retinal angiogenesis in development and disease. Nature. 2005;438:960–6.
Zeng M, Shen J, Liu Y, Lu LY, Ding K, Fortmann SD, et al. The HIF-1 antagonist acriflavine: visualization in retina and suppression of ocular neovascularization. J Mol Med (Berl). 2017;95:417–29.
Ding K, Eaton L, Bowley D, Rieser M, Chang Q, Harris MC, et al. Generation and characterization of ABBV642, a dual variable domain immunoglobulin molecule (DVD-Ig) that potently neutralizes VEGF and PDGF-BB and is designed for the treatment of exudative age-related macular degeneration. MAbs. 2017;9:269–84.
Ozaki H, Seo MS, Ozaki K, Yamada H, Yamada E, Okamoto N, et al. Blockade of vascular endothelial cell growth factor receptor signaling is sufficient to completely prevent retinal neovascularization. Am J Pathol. 2000;156:697–707.
Stitt AW, Curtis TM, Chen M, Medina RJ, McKay GJ, Jenkins A, et al. The progress in understanding and treatment of diabetic retinopathy. Prog Retin Eye Res. 2016;51:156–86.
Chang SH, Hla T. Gene regulation by RNA binding proteins and microRNAs in angiogenesis. Trends Mol Med. 2011;17:650–8.
Hentze MW, Castello A, Schwarzl T, Preiss T. A brave new world of RNA-binding proteins. Nat Rev Mol Cell Biol. 2018;19:327–41.
Platania CBM, Pittalà V, Pascale A, Marchesi N, Anfuso CD, Lupo G, et al. Novel indole derivatives targeting HuR-mRNA complex to counteract high glucose damage in retinal endothelial cells. Biochem Pharmacol. 2020;175: 113908.
Mercer TR, Dinger ME, Mattick JS. Long non-coding RNAs: insights into functions. Nat Rev Genet. 2009;10:155–9.
Kumar MM, Goyal R. LncRNA as a Therapeutic Target for Angiogenesis. Curr Top Med Chem. 2017;17:1750–7.
Zhang N, Cao W, He X, Xing Y, Yang N. Long Non-Coding RNAs in Retinal Ganglion Cell Apoptosis. Cell Mol Neurobiol. 2022;43:561–74.
Wang Y, Wang X, Wang YX, Ma Y, Di Y. Effect and mechanism of the long noncoding RNA MALAT1 on retinal neovascularization in retinopathy of prematurity. Life Sci. 2020;260: 118299.
Wang Y, Wang X, Wang YX, Ma Y, Di Y. The Long-Noncoding RNA TUG1 Regulates Oxygen-Induced Retinal Neovascularization in Mice via MiR-299. Invest Ophthalmol Vis Sci. 2022;63:37.
Bossi L, Figueroa-Bossi N. Competing endogenous RNAs: a target-centric view of small RNA regulation in bacteria. Nat Rev Microbiol. 2016;14:775–84.
Yan B, Yao J, Liu JY, Li XM, Wang XQ, Li YJ, et al. lncRNA-MIAT regulates microvascular dysfunction by functioning as a competing endogenous RNA. Circ Res. 2015;116:1143–56.
Liu P, Jia SB, Shi JM, Li WJ, Tang LS, Zhu XH, et al. LncRNA-MALAT1 promotes neovascularization in diabetic retinopathy through regulating miR-125b/VE-cadherin axis. Biosci Rep. 2019;39:BSR20181469.
Hirbe AC, Gutmann DH. Understanding a complicated Gal-1. Neuro Oncol. 2019;21:1341–3.
Mathieu V, de Lassalle EM, Toelen J, Mohr T, Bellahcene A, Van Goietsenoven G, et al. Galectin-1 in melanoma biology and related neo-angiogenesis processes. J Invest Dermatol. 2012;132:2245–54.
Bastón JI, Barañao RI, Ricci AG, Bilotas MA, Olivares CN, Singla JJ, et al. Targeting galectin-1-induced angiogenesis mitigates the severity of endometriosis. J Pathol. 2014;234:329–37.
Yang N, Zhang W, He T, Xing Y. Silencing of galectin-1 inhibits retinal neovascularization and ameliorates retinal hypoxia in a murine model of oxygen-induced ischemic retinopathy. Exp Eye Res. 2017;159:1–15.
Castello A, Fischer B, Eichelbaum K, Horos R, Beckmann BM, Strein C, et al. Insights into RNA biology from an atlas of mammalian mRNA-binding proteins. Cell. 2012;149:1393–406.
Cui Y, Yan M, Wu W, Lv P, Wang J, Huo Y, et al. ESCCAL-1 promotes cell-cycle progression by interacting with and stabilizing galectin-1 in esophageal squamous cell carcinoma. NPJ Precis Oncol. 2022;6:12.
Liu Y, Ma J, Cui D, Fei X, Lv Y, Lin J. LncRNA MEG3-210 regulates endometrial stromal cells migration, invasion and apoptosis through p38 MAPK and PKA/SERCA2 signalling via interaction with Galectin-1 in endometriosis. Mol Cell Endocrinol. 2020;513: 110870.
Wei J, Li DK, Hu X, Cheng C, Zhang Y. Galectin-1-RNA interaction map reveals potential regulatory roles in angiogenesis. FEBS Lett. 2021;595:623–36.
He Z, Ruan X, Liu X, Zheng J, Liu Y, Liu L, et al. FUS/circ_002136/miR-138-5p/SOX13 feedback loop regulates angiogenesis in Glioma. J Exp Clin Cancer Res. 2019;38:65.
Meng Y, Li S, Zhang Q, Ben S, Zhu Q, Du M, et al. LncRNA-422 suppresses the proliferation and growth of colorectal cancer cells by targeting SFPQ. Clin Transl Med. 2022;12: e664.
Ferrè F, Colantoni A, Helmer-Citterich M. Revealing protein-lncRNA interaction. Brief Bioinform. 2016;17:106–16.
Zhou W, Chen X, Hu Q, Chen X, Chen Y, Huang L. Galectin-3 activates TLR4/NF-κB signaling to promote lung adenocarcinoma cell proliferation through activating lncRNA-NEAT1 expression. BMC Cancer. 2018;18:580.
Zhang PF, Wu J, Luo JH, Li KS, Wang F, Huang W, et al. SNHG22 overexpression indicates poor prognosis and induces chemotherapy resistance via the miR-2467/Gal-1 signaling pathway in epithelial ovarian carcinoma. Aging (Albany NY). 2019;11:8204–16.
Kojima Y, Quertermous T. Apelin-APJ signaling in retinal angiogenesis. Arterioscler Thromb Vasc Biol. 2008;28:1687–8.
Wang YH, Kuo SJ, Liu SC, Wang SW, Tsai CH, Fong YC, et al. Apelin affects the progression of osteoarthritis by regulating VEGF-dependent angiogenesis and miR-150–5p expression in human synovial fibroblasts. Cells. 2020;9:594.
Zhang Y, Zhang Y, Kameishi S, Barutello G, Zheng Y, Tobin NP, et al. The Amot/integrin protein complex transmits mechanical forces required for vascular expansion. Cell Rep. 2021;36: 109616.
Lin X, Spindler TJ, de Souza Fonseca MA, Corona RI, Seo JH, Dezem FS, et al. Super-Enhancer-Associated LncRNA UCA1 Interacts Directly with AMOT to Activate YAP Target Genes in Epithelial Ovarian Cancer. iScience. 2019;17:242–55.
Glaser J, Gonzalez R, Perreau VM, Cotman CW, Keirstead HS. Neutralization of the chemokine CXCL10 enhances tissue sparing and angiogenesis following spinal cord injury. J Neurosci Res. 2004;77:701–8.
Aplin AC, Nicosia RF. Tissue oxygenation stabilizes neovessels and mitigates hemorrhages in human atherosclerosis-induced angiogenesis. Angiogenesis. 2023;26:63–76.
Guo J, Cai H, Zheng J, Liu X, Liu Y, Ma J, et al. Long non-coding RNA NEAT1 regulates permeability of the blood-tumor barrier via miR-181d-5p-mediated expression changes in ZO-1, occludin, and claudin-5. Biochim Biophys Acta Mol Basis Dis. 2017;1863:2240–54.
Xia F, Xu Y, Zhang X, Lyu J, Zhao P. Competing endogenous RNA network associated with oxygen-induced retinopathy: Expression of the network and identification of the MALAT1/miR-124-3p/EGR1 regulatory axis. Exp Cell Res. 2021;408: 112783.
Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014;505:344–52.
Russo AJ, Vasudevan SO, Méndez-Huergo SP, Kumari P, Menoret A, Duduskar S, et al. Intracellular immune sensing promotes inflammation via gasdermin D-driven release of a lectin alarmin. Nat Immunol. 2021;22:154–65.
Freitag N, Tirado-González I, Barrientos G, Herse F, Thijssen VL, Weedon-Fekjær SM, et al. Interfering with Gal-1-mediated angiogenesis contributes to the pathogenesis of preeclampsia. Proc Natl Acad Sci U S A. 2013;110:11451–6.
Dhirapong A, Lleo A, Leung P, Gershwin ME, Liu FT. The immunological potential of galectin-1 and -3. Autoimmun Rev. 2009;8:360–3.
Ni W, Li Z, Ai K. lncRNA ZFPM2-AS1 promotes retinoblastoma progression by targeting microRNA miR-511-3p/paired box protein 6 (PAX6) axis. Bioengineered. 2022;13:1637–49.
Ciccarese F, Grassi A, Pasqualini L, Rosano S, Noghero A, Montenegro F, et al. Genetic perturbation of IFN-α transcriptional modulators in human endothelial cells uncovers pivotal regulators of angiogenesis. Comput Struct Biotechnol J. 2020;18:3977–86.
Mariotti B, Servaas NH, Rossato M, Tamassia N, Cassatella MA, Cossu M, et al. The Long Non-coding RNA NRIR Drives IFN-Response in Monocytes: Implication for Systemic Sclerosis. Front Immunol. 2019;10:100.
Song L, Chen J, Lo CZ, Guo Q, Feng J, Zhao XM. Impaired type I interferon signaling activity implicated in the peripheral blood transcriptome of preclinical Alzheimer’s disease. EBioMedicine. 2022;82: 104175.
Xu H, Wang X, Wu J, Ji H, Chen Z, Guo H, et al. Long Non-coding RNA LINC01094 Promotes the Development of Clear Cell Renal Cell Carcinoma by Upregulating SLC2A3 via MicroRNA-184. Front Genet. 2020;11: 562967.
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:402–8.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.
Liu S, Wang Z, Chen D, Zhang B, Tian R, Wu J, et al. Annotation and cluster analysis of spatiotemporal- and sex-related lncRNA expression in rhesus macaque brain. Genome Res. 2017;27:1608–20.
Kovaka S, Zimin AV, Pertea GM, Razaghi R, Salzberg SL, Pertea M. Transcriptome assembly from long-read RNA-seq alignments with StringTie2. Genome Biol. 2019;20:278.
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39:316–22.
We thank Jiaojiao Wei for her constructive suggestions.
This research was funded by the Hubei Key Laboratories Opening Project (grant no. 2021KFY055), Provincial Natural Science Foundation of Hubei Province (grant no. 2020CFB240), and Fundamental Research Funds for the Central Universities (grant no. 2042020kf0065).
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Fig S1. LGALS1 expression in normal control and siLGALS1 groups quantified by Western blot. a, b The original, unprocessed electrophoretic gel images. Fig S2. LGALS1 regulates the expression of genes associated with angiogenesis in HRMECs. a Top 10 most enriched Gene Ontology (GO) terms associated with upregulated genes between the siLGALS1 and control groups. b Top 10 most enriched GO terms associated with downregulated genes between the siLGALS1 and control groups. c Top 10 most enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) terms associated with upregulated genes between the siLGALS1 and control groups. d Top 10 most enriched KEGG terms of downregulated genes between the siLGALS1 and control groups. HRMECs, human retinal microvascular endothelial cells. Fig S3. Analysis of differentially expressed lncRNAs in HRMECs between LGALS1-knockdown and control groups. a Detected known lncRNAs (left panel) and novel lncRNAs(right panel) in the siLGALS1 and control groups. LncRNAs with FPKM ≥ 0.2 in at least one sample were considered differentially expressed. b Distribution of exon counts of known lncRNAs, novel lncRNAs, and protein-coding RNAs. c Distribution of exon lengths of known lncRNAs, novel lncRNAs, and protein-coding RNAs. d Density of the length distribution of known lncRNAs, novel lncRNAs, and protein-coding RNA. The length density distribution was generated using the density function in R software. e Expression profile of DE lncRNAs. f Top 10 most enriched KEGG pathways associated with DE mRNAs co-expressed with DE lncRNAs in the siLGALS1 and control groups. g Co-expression network of DE lncRNAs and DE mRNAs associated with the top 10 GO terms, shown in red in f. LncRNAs are on the left, co-expressed mRNAs are in the center, and enriched GO terms associated with mRNAs are on the right. DE, differentially expressed; HRMECs, human retinal microvascular endothelial cells; FPKM, fragments per kilobase of exon per million fragments mapped; KEGG, Kyoto Encyclopedia of Genes and Genomes; GO, Gene Ontology. Fig S4. Regulatory network of lncRNA–miRNA–angiogenesis-related genes mediated by LGALS1 deregulation in HRMECs. a Co-expression network of DE lncRNAs and DE angiogenesis-related genes. Orange triangles indicate lncRNAs, and red rhombuses indicate angiogenesis-related genes. Cutoffs of p ≤ 0.01 and PCC ≥ 0.95 were applied to identify co-expression pairs. b Box plots showing the expression profiles of the eight DE lncRNAs involved in the regulatory network. HRMECs, human retinal microvascular endothelial cells; DE, differentially expressed; PCC, Pearson’s correlation coefficient. Fig S5. Validation of the expression of angiogenesis-related genes in HRMECs with or without knockdown LGALS1. a-g Expression of differentially expressed angiogenesis-related genes involved in the network quantified by qPCR. ****p < 0.0001. Table S1. Primer sequences. Table S2. siRNA sequences.
About this article
Cite this article
Yang, N., Zhang, N., Wang, Z. et al. Galectin-1-dependent ceRNA network in HRMECs revealed its association with retinal neovascularization. BMC Genomics 24, 327 (2023). https://doi.org/10.1186/s12864-023-09352-y