- 1Eye Institute and Department of Ophthalmology, Eye & ENT Hospital, Fudan University, Shanghai, China
- 2Key laboratory of Myopia and Related Eye Diseases, NHC, Key laboratory of Myopia and Related Eye Diseases, Chinese Academy of Medical Sciences, Shanghai, China
- 3Shanghai Key Laboratory of Visual Impairment and Restoration, Shanghai, China
- 4Department of Neurology and Geriatrics, Fujian Institute of Geriatrics, Fujian Medical University Union Hospital, Fuzhou, Fujian, China
Objective: This study revealed a core regulator and common upstream mechanisms for the multifaceted pathological processes of age-related macular degeneration (AMD) and provided proof-of-concept for this new therapeutic target.
Methods: Comprehensive gene expression analysis was performed using RNA sequencing of eye cup from old mice as well as laser-induced choroidal neovascularization (CNV) mouse model. Through integrative analysis and protein-protein interaction (PPI) analysis, common pathways and key transcription factor was identified simultaneously engaged in age-related retinal degeneration and CNV, the two typical pathological process of AMD. Subsequently, the expression changes of Spi1, the key regulator, as well as the alternation of the downstream mechanisms were validated in both models through qRT-PCR, Elisa, flow cytometry and immunofluorescence. Further, we assessed the impact of Spi1 knockdown in vitro and in vivo using gene intervention vectors carried by adeno-associated virus or lentivirus to test its potential as a therapeutic target.
Results: Compared to corresponding controls, we found 1,939 and 1,319 genes differentially expressed in eye cups of old and CNV mice respectively. The integrative analysis identified a total of 275 overlapping DEGs, of which 150 genes were co-upregulated. PPI analysis verified a central transcription factor, SPI1. The significant upregulation of Spi1 expression was then validated in both models, accompanied by macrophage polarization towards the M1 phenotype. Finally, SPI1 suppression significantly inhibited M1 polarization of BMDMs and attenuated neovascularization in CNV mice.
Conclusion: This study demonstrates that SPI1 exerts a pivotal role in AMD by regulation of macrophage polarization and innate immune response, offering promise as an innovative target for treating AMD.
Introduction
AMD as a degenerative eye condition that mainly impacts the elderly, is a major contributor to worldwide vision impairment (1–3). This condition arises when the macula gets damaged due to the process of aging. AMD is expected to influence approximately 288 million individuals worldwide by 2040 (4). The precise reason behind AMD remains incompletely comprehended. It is generally recognized to arise from a combination of genetic and environmental elements, with advancing age being the determinant. Risk factors identified so far include smoking (5), obesity, hypertension, hyperlipidaemia (6), and certain genetic variations (7, 8). Drusen and basal linear deposits (BLinD) are the defining features of early-stage AMD. As for late-stage AMD, typically, there are two primary classifications: geographic atrophy (GA) and neovascular AMD (nAMD) (9).
Photoreceptors, retinal pigment epithelium (RPE), and choriocapillaris displayed combined degenerative alterations in GA (9). Over time, the GA area progressively expands, resulting in progressive central vision decline (10). Neovascularization may occur in the subretinal/sub-RPE space, or choroid at any stage of AMD, leading to CNV, which is a hallmark lesion of nAMD. When these newly formed blood vessels leak or rupture, nAMD enters the exudative phase. Accumulated fluid and/or hemorrhage would lead to distorted vision (11). If left untreated, CNV area typically undergo extensive fibrosis, which further deteriorates central vision (3). Over 80% of vision loss in AMD patients is attributed to nAMD (12). VEGF serves as the primary initiator for angiogenesis in CNV (13). And the introduction of anti-VEGF biologics marks a significant advancement in AMD therapy. Since the approval of anti-VEGF therapy in 2006, the number of visually impaired patients aged 80 and above with nAMD has significantly decreased. However, due to individual variability, not all patients accepting anti-VEGF treatment can maintain good vision over time stably. GA patients, on the other hand, did not respond to this therapy due to different pathogenesis. Therefore, exploring the upstream mechanisms of AMD and searching for available therapeutic targets remain urgent.
In our study, transcriptome sequencing of RPE-choroid-sclera complex (eye cup) from old mice as well as laser-induced CNV mice suggested that activated innate inflammatory response is a common mechanism of neovascularization and systemic aging, which was validated in subsequent experiments. Integrative bioinformatics analysis identified PU.1/SPI1 as a core transcription factor (TF) crucial for the development of both CNV and age-related retinal degeneration. The significant rescue in CNV model by Spi1 interference has also highlighted its potential as an innovative target for treating AMD.
Methods
Animals and ethics statement
We obtained 2-month-old male C57BL/6J mice from GemPharmatech and housed them in the animal facility at Fudan University (Shanghai, China). The old C57BL/6J mice (≥18-month-old) were all provided by the research group of Professor Shi-qing Cai from the Center for Excellence in Brain Science and Intelligence Technology, Institute of Neuroscience, Chinese Academy of Sciences. The mice were provided with free access to water and solid diets, and they were maintained under standard housing conditions. Specific experiments indicated the ages and quantities of animals used, detailed either in the figure legends or the methods section. The Animal Care and Use Committee of Fudan University reviewed and approved all animal protocols, ensuring compliance with relevant regulatory standards.
Electroretinography
The retinal function of mice was examined through electroretinography (ERG). ERG assessments were conducted following an overnight dark adaptation. Mice were situated on a heating pad throughout the experiment to ensure their body temperature. All procedures were conducted under red light. ERG was recorded in both eyes by placing electrodes on the corneas, a reference electrode was positioned on the cheek, and a neutral electrode was inserted in the tail. Mice were positioned beneath the ERG dome, and scotopic ERGs were recorded at various light intensities (0.01, 0.1, 1, 3, 10 cd.s/m2). The amplitudes of the a-wave and b-wave were analyzed using Diagnosys version 6.63 software.
OCT
OCT examinations were conducted utilizing ISOCT (Optoprobe Science Ltd). The horizontal scan image was analyzed with the optic nerve head as the center. We used the built-in software to measure total retinal thickness at distances of 0.2, 0.4, 0.6, 0.8, and 1.0 mm from the optic nerve head (nasally and temporally).
RNA-Seq and bioinformatics analysis
We conducted RNA-seq analysis at Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China). TRIzol was utilized to extract total RNA from the eye cup. The cDNA library was prepared by 1μg of total RNA. Differentially expressed genes (DEGs) between two distinct groups were identified based on the transcripts per million reads (TPM) method. The datasets PRJNA1102502 and PRJNA1102537 in this study are available on NCBI. DEG selection criteria included |log2(FC)| > 1.0 and FDR < 0.05. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were conducted utilizing Metascape database (https://metascape.org/). The heatmap was generated on https://www.bioinformatics.com.cn. In addition, Venn’s diagram was made on https://bioinfogp.cnb.csic.es/tools/venny/index.html. We performed PPI network analysis among proteins corresponding to DEGs via STRING database (http://string-db.org). We determined key proteins in the interaction network based on betweenness centrality, and the PPI network was established through Cytoscape software (v3.8.2).
Laser-induced CNV
Mice were given anesthesia through intraperitoneal injection of 1.25% tribromoethanol (0.3 mL/20 g) and their pupils were dilated using tropicamide. Approximately 2 minutes later, both eyes were topically treated with carbomer gel, and a small circular glass slide was placed at the center of the cornea. Using a laser photocoagulator with a wavelength of 532nm (VITRA, France), 4 or 6 laser burns were induced at 2 PD equal distances from the optic disc (4 burns for immunofluorescence experiments and 6 burns for sequencing and molecular experiments). The laser parameters were set at a power of 110mW and duration of 140ms. Successful disruption of the Bruch’s membrane was indicated by the generation of white bubbles upon photocoagulation. The laser burns were positioned to avoid retinal blood vessels. If the laser did not disrupt the Bruch’s membrane or caused retinal hemorrhage, it was excluded from analysis.
Intravitreal injection
Intravitreal injections of AAV NC and shSpi1 (Obio Technology, Shanghai) were performed in the left and right eyes of 2-month-old C57BL/6J mice following pupil dilation respectively. Three weeks later, CNV model was conducted.
Immunofluorescence
Following euthanasia of the mice, the eye globes were removed and immersed in paraformaldehyde for 30 minutes. The eye cup was dissected to remove the neuroretina and incubated in PBS containing 5% BSA and 0.5% Triton X-100 for 1 hour. Primary antibodies were applied for overnight incubation thereafter. FITC coagulated Isolectin B4 (1:100, Sigma) was used to examine blood vessels, while Iba1 (1:500, ab178846, Abcam) was utilized to identify mononuclear phagocytes. Following multiple rinses with PBS, the eye cup was then exposed to a secondary Cy3-conjugated goat anti-rabbit antibody for a duration of 1 hour. After that, the tissue was placed on a glass slide and photographed with Olympus FV3000 confocal microscope. The Z-stacks were combined into a single image for examination.
BMDM culture
Macrophages derived from bone marrow of mice of various ages were isolated using established methods, then cultured in DMEM including 10% FBS, 1% penicillin/streptomycin, and 10ng/ml M-CSF (PeproTech). Enriched bone marrow–derived macrophages (BMDMs) can be collected for experiments after 7 days of cell culture. For Spi knockdown, BMDMs were infected with NC or shSpi1 lentivirus (Obio Technology, Shanghai) for 48 hours and harvested after 24-hour LPS stimulation.
ELISA
We utilized ELISA kits from Multisciences Biotech to assess the protein concentrations of IL-6, IL-1β, and TNF-α in both cell culture supernatant and eye cup samples, according to the guidelines provided by the manufacturers.
qRT-PCR
Tissues and cells were processed to extract RNA with EZBioscience kits, followed by cDNA transcription with the Color Reverse Transcription Kit from EZBioscience, according to the guidelines provided by the manufacturers. Roche LightCycler 480 II system (Roche, Switzerland) and SYBR Green (EZBioscience) were utilized to perform quantitative real-time PCR (qRT-PCR). For PCR, 2μl of cDNA was included in the 10μl overall volume as the template. The temperature cycling parameters consisted of heating to 95°C for 5 minutes, then undergoing 40 cycles of heating to 95°C for 10 seconds followed by cooling to 60°C for 30 seconds. The primer efficiency was confirmed. Following the reactions, the raw data were exported and analyzed by the 2−△△CT method. The primers used are listed as follows. Il-6: 5’-TAGTCCTTCCTACCCCAATTTCC-3’ and 5’-TTGGTCCTTAGCCACTCCTTC-3’; Il-1β: 5’-GCAACTGTTCCTGAACTCAACT-3’ and 5’-ATCTTTTGGGGTCCGTCAACT-3’; Tnf-α: 5’-CCCTCACACTCAGATCATCTTCT-3’ and 5’-GCTACGACGTGGGCTACAG-3’; Spi1: 5’-ATGTTACAGGCGTGCAAAATGG-3’ and 5’-TGATCGCTATGGCTTTCTCCA-3’; Actb: 5’-GGCTGTATTCCCCTCCATCG-3’ and 5’-CCAGTTGGTAACAATGCCATGT-3’.
Flow cytometry
After washing single-cell suspensions of BMDM or eye cups, they were incubated on ice for 30 minutes with specific fluorescent antibodies (PE-anti-F4/80, Biolegend, 123110; FITC-anti-CD11b, Biolegend, 101206; and APC-anti-CD11c, MULTI Sciences, AM011C05–250; or Pacific Blue-anti-CD86, Biolegend, 105123) in the dark. Samples were examined using a FACS Vantage flow cytometer from BD Biosciences, and FlowJo software was utilized to process the data.
Statistical analysis
Analysis and figure generation were performed using GraphPad Prism 9 software. Mean ± SEM. values are depicted on the graphs. Statistical comparisons between two groups were conducted using a Student’s t-test with unequal variance assumed, while ERG and OCT data were analyzed using two-way ANOVA. P < 0.05 was considered statistically significant.
Results
Degenerative phenotypes and transcriptome changes in eyes of old mice
As mice age, the retina experiences degenerative alterations, including a reduction in total retinal thickness (Figures 1A, B) and decrease of both the a and b wave amplitude in electroretinogram (Figures 1C, D). To investigate the underlying causes of these changes, we conducted transcriptome sequencing on the eye cups of 30-month-old elderly mice (n=5) and 3-month-old young controls (n=5). Principal component analysis (PCA) indicated that samples separated based on their experimental category (Figure 2A). We found 1,939 genes that were differentially expressed in old eye cup, including 1,332 upregulated and 607 downregulated genes (Figure 2B). A heatmap was generated for the top 20 genes that showed increased and decreased expression levels (Figure 2C). GO and KEGG analysis was conducted for the DEGs to illustrate functional enrichment. Top 20 upregulated GO enrichments are shown in Figures 2D, including positive regulation of immune/inflammation response, cytokine production, leukocyte activation, and immune effector process. These findings align with previous research suggesting the significant involvement of the immune response in eye aging. Among the top 20 downregulated GO enrichment (Figure 2E), extracellular matrix (ECM), neuron projection development, regulation of synapse structure or activity, vasculature development and cell junction organization were enriched. Figure 2F demonstrated that the KEGG analysis of increased DEGs in the old eye cup was predominantly enriched in cytokine-cytokine receptor interaction, complement and coagulation cascades, phagosome, NK cell mediated cytotoxicity, Th1 and Th2 cell differentiation, IL17, JAK-STAT and NF-κB signaling pathway. The main pathways involving down-regulated DEGs included ECM–receptor interaction, calcium signal pathway, proteoglycans in cancer, as well as protein digestion and absorption (Figure 2F).
Figure 1 Degenerative changes in the retina of old mice. (A) Representative OCT images of old (18-month) and young (3-month) mice. (B) Total retinal thickness at varying distances (–1.0 mm to +1.0 mm) from the optic disk center (0 μm) in increments of 0.2 mm for both old and young mice, n =7. (C, D) Amplitudes of the a-waves (C) and b-waves (D) at different intensities in old (18-month) and young (3-month) mice, n=8. *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001. All data were demonstrated as mean ± SEM.
Figure 2 Transcriptomic analysis of the eye cups from old mice. (A) Principal component analysis of old (30-month) and young (3-month) eye cup, n=5. Old eye cups are marked by red and young eye cup are marked by blue. (B) Plot displaying DEGs in a volcano format. (C) Heatmap illustrating the top 20 DEGs identified between old and young eye cup. (D, E) Top 20 enriched GO terms among upregulated (D) and downregulated (E) DEGs in old eye cup. (F) Top 15 upregulated and 10 downregulated KEGG pathways in old eye cup.
Transcripts and biological pathways altered in CNV eye cup
To investigate the pathogenesis of CNV, the other subtype of AMD, we used laser photocoagulation to construct a mouse model of CNV and performed transcriptome sequencing on eye cups from CNV mice (n=4) and age-matched control mice (n=4). We found that 1,319 genes were differentially expressed in eye cups with CNV, of which 747 were upregulated and 572 were downregulated (Figures 3A, B). A heatmap format displayed the top 20 genes that were upregulated and downregulated (Figure 3C). Figures 3D, E displayed top 20 and top 15 GO enrichments of DEGs. We found that extracellular matrix organization, leukocyte activation, inflammation/immune response and positive regulation of cytokine and chemotaxis production were enriched among the upregulated genes. Conversely, the downregulated DEGs were found to be predominantly associated with carboxylic acid transport, sensory perception of light stimulus and vitamin transport. As shown in Figure 3F, the KEGG indicated that upregulated DEGs were predominantly concentrated in complement and coagulation cascades, cytokine-cytokine receptor interaction, NF-κB and toll-like receptor signaling pathway. The main pathways involving downregulated genes in CNV group included the vitamin digestion and absorption, carbon metabolism, tight junction, and retinol metabolism. Notably, ECM receptor interactions were enriched among both upregulated and downregulated genes, suggesting that genes within the same term could also participate in CNV development in diverse ways.
Figure 3 Transcriptomic analysis of the CNV eye cup. (A) Principal component analysis of CNV and control eye cup, n=4. Control eye cups are marked by green and CNV eye cups are marked by red. (B) Plot displaying DEGs in a volcano format. (C) Heatmap depicting the top 20 DEGs between CNV and control eye cup. (D, E) The top 20 enriched GO terms among upregulated (D) and downregulated (E) DEGs in CNV group. (F) Top 15 upregulated and 10 downregulated KEGG pathways in CNV group.
Common regulators involved in pathology of both aging and CNV
Considering that CNV and age-related retinal degeneration are two representative aspects of pathologic changes of AMD, we conducted integrative analysis on the 2 sets of transcriptomic data described above to reveal their common pathogenesis. We identified a total of 275 overlapped DEGs (Figure 4A). Out of the total, 150 genes showed simultaneous upregulation, while 13 genes exhibited simultaneous downregulation in both old and CNV mice. We conducted a PPI network analysis on the proteins corresponding to the 150 co-upregulated DEGs in both old eyes and CNV eyes to determine the main pathological factors contributing to AMD at the molecular level (Figure 4B). The calculation of betweenness centrality for each node identified the key proteins that act as communication bottlenecks in the network. Node size reflects the degree of connections while node color represents the betweenness centrality. Red markers indicate the core proteins with the highest values in the figure. Generally, high betweenness centrality implies important biological roles and could be a hint of potential therapeutic targets. Top 10 genes with highest betweenness centrality are Tyrobp, Spi1, Ctss, Rac2, Cxcl10, Fcgr3, Fcer1g, Cd68, Timp1 and Fcgr1. The betweenness centrality and fold change of these genes are listed in Figure 4C. Subsequently, the elevated mRNA expression level of a TF, Spi1, was verified in eyes from old mice and CNV mice compared with their respective controls (Figures 4D, E).
Figure 4 Integrative analysis of the old and CNV transcriptomics. (A) Venn diagram for overlapped DEGs in the old and CNV transcriptomes. (B) PPI network analysis of the proteins corresponding to the 150 co-upregulated DEGs in old and CNV eyes. Red: high betweenness centrality, yellow: moderate betweenness centrality, purple: low betweenness centrality. (C) Top 10 DEGs with highest betweenness centrality and their fold change in old and CNV transcriptomes. (D) mRNA expression of Spi1 in old (30-month) and young (3-month) eye cups by qRT-PCR, n=5. (E) mRNA expression of Spi1 in CNV and control eye cups by qRT-PCR, n=5. *P<0.05, **P<0.01. All data were demonstrated as mean ± SEM.
Systemic pro-inflammatory status in old mice
According to the results of scRNA-seq in the database, Spi1 is highly expressed in macrophages, positively correlated to innate immune responses. We performed co-staining of IB4 and Iba1 on choroid flat mounts of young and old mice to label activated mononuclear phagocytes (MNPs) dislocated in the subretinal space and quantitatively analyze them (Figures 5A, B). Unsurprisingly, young mice had almost no or very few activated MNPs in the subretinal space, whereas old mice exhibited an elevation in the number. qRT-PCR analysis of pro-inflammatory factors of eye cup, including Il-6, Il-1β, and Tnf-α revealed a significant elevation in old group (Figure 5C). BMDMs were cultured from both young and old mice for in vitro experiments. And we conducted additional investigations on the expression differences in Spi1 and pro-inflammatory factors at cellular level. We confirmed that Spi1 and pro-inflammatory factors were elevated in BMDMs from old mice (Figure 6A). Additionally, the concentrations of inflammatory proteins in the cell supernatant were notably higher in the older group (Figure 6B). Flow cytometry analysis was conducted on BMDMs utilizing F4/80 and CD11b markers to distinguish macrophages, as well as the CD11c marker to identify M1 macrophages (Figures 6C, D). As shown in Figure 6E, the study found that the percentage of M1 macrophages in the elderly group was 1.5 times greater than the young mice,.
Figure 5 Eye cup from old mice exhibited a pro-inflammatory state. (A) Representative confocal images of IB4 and Iba1 stained activated MNPs of choroid flat mounts in old (24-month) and young (3-month) mice. (B) Numbers of IB4+Iba1+ MNPs in subretinal space of young and old mice, n=5. (C) Pro-inflammatory factors in old (30-month) and young (3-month) eye cup quantified by qRT-PCR, n=5. *P<0.05, ***P<0.001. All data were demonstrated as mean ± SEM.
Figure 6 BMDMs from old mice presented a pro-inflammatory state. (A, B) Spi1 and pro-inflammatory factors in old (21-month) and young (3-month) BMDMs quantified by qRT-PCR (A) and Elisa (B), n=5. (C) Flow cytometry analysis strategy of BMDMs from old and young mice. (D) Representative flow plots of M1 macrophage in BMDMs of old and young mice. (E) Quantification of the proportion of M1 macrophages in BMDMs from old and young mice, n=5. *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001. All data were demonstrated as mean ± SEM.
Pro-inflammatory alternation and M1 polarization were shared in CNV eye cup
To investigate whether SPI1 also affects the inflammatory status of the eye cup in CNV mice, qRT-PCR and ELISA analysis were performed, revealing an elevation in pro-inflammatory factors (Figures 7A, B). Flow cytometry analysis of eye cup was conducted on control and CNV mice, showing an increased trend in the proportion of F4/80+CD11b+ macrophages (Figures 7C, D) as well as a higher proportion of M1 macrophages (Figures 7E, F) in eye cup from CNV mice.
Figure 7 Eye cup from CNV mice displayed a pro-inflammatory state. (A, B) Pro-inflammatory factors in CNV and control eye cup quantified by qRT-PCR (A) and Elisa (B), n=5. (C) Representative flow plots of macrophages in eye cup of CNV and control mice. (D) The proportion of macrophages in CNV and control group, n=5. (E) Representative flow plots of M1 macrophages in eye cup of CNV and control mice. (F) Quantification of the proportion of M1 macrophages in CNV and control mice, n=5. *P<0.05, **P<0.01, ****P<0.0001, ns, not significant. All data were demonstrated as mean ± SEM.
Knockdown of Spi1 inhibited M1 polarization of macrophages and alleviated CNV
To further validate the role of SPI1 in inflammatory response and AMD, we assessed the neovascular process in CNV mice with Spi1 knockdown. 2-month-old mice were intravitreally injected with AAV carrying either NC or Spi1 knockdown constructs in their left and right eyes, respectively. After three weeks, CNV modeling was performed, and choroid flat mounts were obtained one week later for IB4 staining to label the area of neovascularization. It was observed that Spi1 knockdown significantly reduced laser-induced neovascularization area (Figures 8A, B). In vitro cell experiments were also performed on extracted BMDMs with Spi1-knockdown lentivirus. After 48 hours, the cells were exposed to LPS for an additional 24 hours before being harvested to assess polarization. It was found that Spi1 knockdown markedly inhibited macrophage polarization towards the M1 phenotype (Figures 8C, D), leading to a notable decrease in pro-inflammatory factors (Figures 8E, F).
Figure 8 Spi1 knockdown alleviated CNV by inhibiting M1 polarization. (A) Representative confocal images of IB4-stained neovascularization area in NC and shSpi1 eye cups of CNV mice. (B) Quantification of the IB4-labeled CNV area in NC and shSpi1 groups (n=8). (C) Representative flow plots of M1 macrophages in NC and shSpi1 BMDMs. (D) The proportion of M1 macrophages in NC and shSpi1 BMDMs, (n=3). (E, F) Pro-inflammatory factors in NC and shSpi1 BMDMs quantified by qRT-PCR (E) and Elisa (F), (n=3). *P<0.05, **P<0.01, ***P<0.001. All data were demonstrated as mean ± SEM.
Discussion
Inflammation activation is similar in age-related retinal disorders and CNV, which are two major pathological changes of AMD. We performed integrative analysis of old and CNV transcriptomic features to confirm the broad proinflammatory state and identify a core transcription factor, SPI1. Our findings revealed a significant upregulation of SPI1 expression in both old mice and CNV models, suggesting that it is closely related to aging and neovascularization. Through gene knockdown experiments targeting Spi1 in BMDM and CNV mice, we observed a significant inhibition of macrophage M1 polarization and a reduction in neovascularization, suggesting it as a potential therapeutic target for AMD.
Our study of the old mice provided insights into both the phenotypic and transcriptomic changes associated with aging in the eye. The process of aging is intricate, impacting different facets of the structure and function of the retina, ultimately raising the likelihood of developing age-related retinal conditions. The phenotypic findings we obtained aligned with prior research, revealing significant alterations in the aging mouse retina such as reduced thickness of retinal layers and diminished a and b wave amplitudes in electroretinogram recordings (14, 15). Furthermore, these alterations are accordance with either functional decline or structural abnormalities often seen in age-related retinal disorders, such as AMD (15), providing a basis for understanding the development of retinal pathologies with natural aging. Complementing the phenotypic changes, our transcriptomic analysis revealed significant alterations in gene expression profiles in the eye cup from old mice. It was discovered that osteocytes could promote fat formation and constrain bone formation of bone marrow mesenchymal stem cells through the secretion of neuropeptide Y (NPY). and the expression of NPY rises with aging and in cases of osteoporosis (16). We found that NPY increased to the highest degree in old mice, which may lead to the occurrence of degenerative diseases due to the disturbance of retinal metabolic regulation. Previous studies employed senescence-accelerated mouse prone 8 (SAMP8) as an experimental model to investigate neurodegenerative diseases. It was found that MOBP expression is downregulated in SAMP8 mice, potentially playing a role in age-related brain diseases (17). However, our sequencing revealed that MOBP is significantly upregulated in the eye cup of old mice. This could be due to tissue differences as well as differences in normal aging and aging induced by genetic editing. Among the remaining significantly upregulated genes, some are functionally related to aging immune regulation, such as Trgc, Pdcd1, S100a8, and others. Administering PD-1 antibody to aging mice or non-alcoholic steatohepatitis (NASH) mice lead to the reduction on the number of p16+ cells as well as the activated CD8+ T cell-dependent PD-L1+ population in vivo, ameliorating various aging-related phenotypes (18). S100A8-positive microglia are upregulated in tau disease mouse models and SAMP8 mice, and S100A8+ microglia were found in brain biopsies from individuals with AD and elderly patients with no reported pathology, which may be associated with cognitive dysfunction (19). GO and KEGG reflected a wide range of biological processes and molecular pathways implicated in ocular aging and AMD. Aligning with previous studies, inflammation response and ECM remodeling turned to be the most dysregulated pathways, highlighting their irreplaceable role (20, 21). Under typical circumstances, the blood-retinal barrier (BRB) collaborates with immune effectors to uphold retinal balance. With aging, both the retinal vascular endothelial cells and RPE cells experience age-related functional impairments, including the weakening of physical barriers such as the inner and outer BRBs (22, 23). As for the immune effectors of the retina, including microglia and the complement system, they also experience notable alterations with age (24, 25). On one side, senescent microglia show reduced capacity in immune surveillance and synaptic maintenance, bringing about the deposition of defective neurons and fragments. Conversely, microglia and the complement system exhibit an increased pro-inflammatory response to damaged neurons and debris, contributing to immune pathology. These aging-related alterations would violate intraocular microenvironment and make it more prone to conditions like AMD, diabetic retinopathy, and glaucoma (26, 27).
In addition to the degenerative changes occurring with aging, CNV, which is identified by the abnormal expansion of blood vessels from choroid into retina, is another major subtype of AMD that leads to over 80% of vision impairment in affected individuals. We conducted transcriptome sequencing of eye cups from a mouse model with laser induced CNV to discover important genes and pathways related to CNV. Significant changes in gene expression pattern were revealed, encompassing multiple signaling pathways mainly associated with inflammation, tight junction, and ECM remodeling (11, 26, 28). Among them, changes in various cytokines were the most striking, including IL-1β, IL-2, and IL-6. It has been reported that they not only participate in neovascularization and extensive fibrosis but are also closely associated with dry AMD. The effectiveness of IL-1Ra in reducing CNV is attributed to blocking IL-1α and IL-1β (29–31). In our study, we found a significant upregulation of Il1rn in CNV mice. This may be due to a protective response mounted by the organism in response to the disease state. However, as the disease progresses and exceeds the body’s self-regulatory capacity, this response may become insufficient. As reported in the literature, in the early stages of CNV, local administration of IL-1ra significantly inhibits CNV development. However, treatment initiated one week after CNV induction only transiently suppresses the vascular response (32). IL-2 is pivotal in fibrosis of macular degeneration via participation in the migration of RPE cells, synthesis of ECM, and expression of TGF-β2 (33). IL-6 holds importance in promoting subretinal fibrosis and is involved in the repair of damaged organs (34). Serum IL-6 levels are associated with the course of GA (35). Additionally, studies have shown that the activation of STAT3 through IL-6 receptor contributes to the development of CNV, with higher levels of IL-6 in the aqueous humor of AMD patients correlating with increased CNV size and activity (34–36). In the DEGs of CNV mice, a significant downregulation of several genes related to neuronal function is observed. OMG, LRTM2, GRIK1, and CDH9 are predominantly expressed in the nervous system, participating in neuronal and synaptic functions, has a significant impact on the occurrence of various neurological disorders. LRAT and RDH10 play crucial roles in the visual process, participating in the vitamin A metabolic pathway, which is essential for normal visual system function. Mutations in LRAT and RDH10 genes are essential for the occurrence of retinal degenerative diseases.
The identification of shared KEGG pathways involved in inflammation and ECM remodeling altered in eye cups from both old mice and the CNV model suggests common pathogenic mechanisms underlying. Considering that aging-related retinal degeneration and CNV progression are two classic branches of AMD lesions, we performed integrative analysis of these two set of transcriptomics in search of potential upstream therapeutic targets for AMD. It was revealed that 150 genes were significantly co-upregulated. Subsequently, PPI analysis was conducted and the 150 DEGs were sorted according to their betweenness centrality. While the top 10 core DEGs have distinct functions, they share some common characteristics. Overall, these genes collectively contribute to immune responses, inflammation, and immune cell functions, coinciding the common pathology shared by two subtypes of AMD. Their coordinated actions are crucial for maintaining immune homeostasis and defending against infectious agents as well as other challenges to the immune system. It is noteworthy that Timp1 is drastically upregulated in the CNV group to a much greater extent than the other nine genes. TIMP1 exerts its function by forming a complex with metalloproteinases, inhibiting their activity, and thereby maintaining the stability of the ECM and the normal structure of tissues. It has significant effects on cell migration, tumor development, and angiogenesis. Epoxomicin may lead to an imbalance between MMP-2 and TIMP-1 in ARPE-19 cells, resulting in the occurrence of retinal fibrosis (37). Since in most cases, transcription factors (TFs) are always the key regulators of gene expression and cellular functions, and their dysregulation has been implicated in various pathological conditions including neovascular diseases (38, 39), we focused on the only TF among the top 10 pivotal DEGs, SPI1, for further study.
SPI1, also referred to as PU.1, is mainly found in hematopoietic cells and lymphocytes, and is recognized as a key controller of myeloid cell development and macrophage terminal differentiation (40). Prior investigations indicated a positive association between SPI1 expression and the progression of AD as well as Aβ deposition. To be specific, SPI1 binding motifs were found to be enriched in the promoters of genes associated with disease, and SPI1 was observed to promote disease advancement by transcriptionally regulating these downstream genes (41). Additionally, SPI1 can modulate liver metabolic function by influencing macrophage polarization, making it a promising therapeutic target for type 2 diabetes (T2D) and NAFLD/NASH (42). To date, the majority of research concerning SPI1 in the retina has predominantly centered on ocular inflammation, including conditions such as uveitis, optic nerve injury, and photoreceptor malnutrition (43–45). Significant inhibition of IFN-γ and IL-2 production stimulated by IRBP was observed with SPI1 knockdown, alleviating autoimmune uveoretinitis (44). What’s more, SPI1 can synergistically regulate the activation of microglia with IRF8 and is involved in neurodegeneration caused by traumatic injury (45). Sun et al. also discovered that blocking SPI1 could reverse the damage to RPE caused by Aβ by decreasing the production of ROS and alleviating mitochondrial dysfunction (46). Given the active involvement of macrophages in the progress of CNV and aging which was confirmed in our study, we were curious whether SPI1 also contributed significantly to the onset of AMD through regulating mononuclear macrophage system, which was finally proved in subsequent experiments. As reported in the literature and observed in our experiments, macrophages are essential in the development of both CNV and age-related retinal degenerative alterations, particularly the M1 phenotype being associated with pro-inflammatory responses and angiogenesis promotion (47, 48). This is to be expected, as prolonged intraocular damage would lead to the overactivation of microglia, specific macrophages resident in the eye, secretion of multiple pro-inflammatory and cytotoxic factors as well. They create a pro-inflammatory environment and in turn further recruit more retinal microglia as well as exogenous monocytes. Sustained inflammatory damage is central to the initiation of GA, providing a breeding ground for angiogenesis (49, 50). Our findings further highlighted the involvement of SPI1 in macrophage polarization and its impact on disease progression. The suppression of Spi1 expression led to significant reduction in CNV as well as alleviation of M1 polarization in BMDMs. This suggests that by attenuating the M1 polarization of macrophages, SPI1 inhibition may help modulate the inflammatory microenvironment and reduce the neovascularization process, potentially offering a new treatment approach for AMD.
Although our study focused on SPI1 and its impact on CNV, it is important to acknowledge that CNV pathogenesis is a complex process involving multiple factors and cellular interactions. Other transcription factors, signaling pathways, and cellular components may also contribute to CNV development and warrant further investigation. In addition, we only performed Spi1 knockdown in BMDMs to examine its regulation of the production of pro-inflammatory factors in vitro. Phenotypic analysis of old mice with SPI1 interference in vivo is helpful to further understand its role in immunosenescence.
In summary, using transcriptomic sequencing in eye cups from old and CNV mice, we demonstrated the hyperactivation of innate immune response as a shared alternation of biological process, and identified a core TF SPI1. Subsequently, Spi1 was validated to be upregulated not only in the eye cup from two groups but also senescent BMDMs, accompanied by polarization of macrophages towards the M1 phenotype. Knocking down Spi1 resulted in the inhibition of M1 polarization of macrophages, relieved the pro-inflammatory microenvironment and alleviated CNV. That’s to say, SPI1 holds great potential as a new therapeutic target for AMD because of its regulatory role in the general pathologic mechanism upstream.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The animal study was approved by The Animal Care and Use Committee of Fudan University. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
SQ: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. YZ: Conceptualization, Formal analysis, Investigation, Methodology, Project administration, Software, Supervision, Writing – original draft, Writing – review & editing. LK: Data curation, Formal analysis, Investigation, Methodology, Software, Supervision, Visualization, Writing – review & editing. DB: Investigation, Methodology, Validation, Visualization, Writing – review & editing. HK: Writing – review & editing, Investigation, Methodology, Software. SZ: Supervision, Validation, Writing – review & editing, Visualization. CZ: Funding acquisition, Project administration, Resources, Supervision, Validation, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by National Natural Science Foundation of China (82020108006 and 81730025); the Open Research Funds of the State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University.
Acknowledgments
We thank Professor Shiqing Cai for providing the old C57BL/6J (≥18-month-old) mice and experimental platform. We also thank online website (https://www.bioinformatics.com.cn and https://bioinfogp.cnb.csic.es/tools/venny/index.html) for providing data analysis and visualization support.
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Jonas JB, Cheung CMG, Panda-Jonas S. Updates on the epidemiology of age-related macular degeneration. Asia Pac J Ophthalmol (Phila). (2017) 6:493–7. doi: 10.22608/apo.2017251
2. Coleman HR, Chan CC, Ferris FL 3rd, Chew EY. Age-related macular degeneration. Lancet. (2008) 372:1835–45. doi: 10.1016/S0140-6736(08)61759-6
3. Lim LS, Mitchell P, Seddon JM, Holz FG, Wong TY. Age-related macular degeneration. Lancet. (2012) 379:1728–38. doi: 10.1016/S0140-6736(12)60282-7
4. Wong WL, Su X, Li X, Cheung CM, Klein R, Cheng CY, et al. Global prevalence of age-related macular degeneration and disease burden projection for 2020 and 2040: a systematic review and meta-analysis. Lancet Glob Health. (2014) 2:e106–16. doi: 10.1016/S2214-109X(13)70145-1
5. Seddon JM, Willett WC, Speizer FE, Hankinson SE. A prospective study of cigarette smoking and age-related macular degeneration in women. Jama. (1996) 276:1141–6. doi: 10.1001/jama.1996.03540140029022
6. Cheung CM, Wong TY. Is age-related macular degeneration a manifestation of systemic disease? New prospects for early intervention and treatment. J Intern Med. (2014) 276:140–53. doi: 10.1111/joim.12227
7. Seddon JM, Francis PJ, George S, Schultz DW, Rosner B, Klein ML. Association of CFH Y402H and LOC387715 A69S with progression of age-related macular degeneration. Jama. (2007) 297:1793–800. doi: 10.1001/jama.297.16.1793
8. Fritsche LG, Igl W, Bailey JN, Grassmann F, Sengupta S, Bragg-Gresham JL, et al. A large genome-wide association study of age-related macular degeneration highlights contributions of rare and common variants. Nat Genet. (2016) 48:134–43. doi: 10.1038/ng.3448
9. Ferris FL 3rd, Wilkinson CP, Bird A, Chakravarthy U, Chew E, Csaky K, et al. Clinical classification of age-related macular degeneration. Ophthalmology. (2013) 120:844–51. doi: 10.1016/j.ophtha.2012.10.036
10. Fleckenstein M, Mitchell P, Freund KB, Sadda S, Holz FG, Brittain C, et al. The progression of geographic atrophy secondary to age-related macular degeneration. Ophthalmology. (2018) 125:369–90. doi: 10.1016/j.ophtha.2017.08.038
11. Fleckenstein M, Keenan TDL, Guymer RH, Chakravarthy U, Schmitz-Valckenberg S, Klaver CC, et al. Age-related macular degeneration. Nat Rev Dis Primers. (2021) 7:31. doi: 10.1038/s41572-021-00265-2
12. Klein R, Peto T, Bird A, Vannewkirk MR. The epidemiology of age-related macular degeneration. Am J Ophthalmol. (2004) 137:486–95. doi: 10.1016/j.ajo.2003.11.069
13. Grossniklaus HE, Ling JX, Wallace TM, Dithmar S, Lawson DH, Cohen C, et al. Macrophage and retinal pigment epithelium expression of angiogenic cytokines in choroidal neovascularization. Mol Vis. (2002) 8:119–26.
14. Ferdous S, Liao KL, Gefke ID, Summers VR, Wu W, Donaldson KJ, et al. Age-related retinal changes in wild-type C57BL/6J mice between 2 and 32 months. Invest Ophthalmol Vis Sci. (2021) 62:9. doi: 10.1167/iovs.62.7.9
15. Ryoo NK, Ahn SJ, Park KH, Ahn J, Seo J, Han JW, et al. Thickness of retina and choroid in the elderly population and its association with Complement Factor H polymorphism: KLoSHA Eye study. PloS One. (2018) 13:e0209276. doi: 10.1371/journal.pone.0209276
16. Zhang Y, Chen CY, Liu YW, Rao SS, Tan YJ, Qian YX, et al. Neuronal induction of bone-fat imbalance through osteocyte neuropeptide Y. Adv Sci (Weinh). (2021) 8:e2100808. doi: 10.1002/advs.202100808
17. Zhang S, Li H, Zheng L, Li H, Feng C, Zhang W. Identification of functional tRNA-derived fragments in senescence-accelerated mouse prone 8 brain. Aging (Albany NY). (2019) 11:10485–98. doi: 10.18632/aging.102471
18. Wang TW, Johmura Y, Suzuki N, Omori S, Migita T, Yamaguchi K, et al. Blocking PD-L1-PD-1 improves senescence surveillance and ageing phenotypes. Nature. (2022) 611:358–64. doi: 10.1038/s41586-022-05388-4
19. Gruel R, Bijnens B, Van Den Daele J, Thys S, Willems R, Wuyts D, et al. S100A8-enriched microglia populate the brain of tau-seeded and accelerated aging mice. Aging Cell. (2024):e14120. doi: 10.1111/acel.14120
20. Yi W, Lu Y, Zhong S, Zhang M, Sun L, Dong H, et al. A single-cell transcriptome atlas of the aging human and macaque retina. Natl Sci Rev. (2021) 8:nwaa179. doi: 10.1093/nsr/nwaa179
21. Chen M, Luo C, Zhao J, Devarajan G, Xu H. Immune regulation in the aging retina. Prog Retin Eye Res. (2019) 69:159–72. doi: 10.1016/j.preteyeres.2018.10.003
22. López-Luppo M, Catita J, Ramos D, Navarro M, Carretero A, Mendes-Jorge L, et al. Cellular senescence is associated with human retinal microaneurysm formation during aging. Invest Ophthalmol Vis Sci. (2017) 58:2832–42. doi: 10.1167/iovs.16-20312
23. Chen M, Rajapakse D, Fraczek M, Luo C, Forrester JV, Xu H. Retinal pigment epithelial cell multinucleation in the aging eye - a mechanism to repair damage and maintain homoeostasis. Aging Cell. (2016) 15:436–45. doi: 10.1111/acel.12447
24. Chen M, Xu H. Parainflammation, chronic inflammation, and age-related macular degeneration. J Leukoc Biol. (2015) 98:713–25. doi: 10.1189/jlb.3RI0615-239R
25. Xu H, Chen M. Targeting the complement system for the management of retinal inflammatory and degenerative diseases. Eur J Pharmacol. (2016) 787:94–104. doi: 10.1016/j.ejphar.2016.03.001
26. Ardeljan D, Chan CC. Aging is not a disease: distinguishing age-related macular degeneration from aging. Prog Retin Eye Res. (2013) 37:68–89. doi: 10.1016/j.preteyeres.2013.07.003
27. Luu J, Palczewski K. Human aging and disease: Lessons from age-related macular degeneration. Proc Natl Acad Sci USA. (2018) 115:2866–72. doi: 10.1073/pnas.1721033115
28. de Jong S, Tang J, Clark SJ. Age-related macular degeneration: A disease of extracellular complement amplification. Immunol Rev. (2023) 313:279–97. doi: 10.1111/imr.13145
29. Lavalette S, Raoul W, Houssier M, Camelo S, Levy O, Calippe B, et al. Interleukin-1β inhibition prevents choroidal neovascularization and does not exacerbate photoreceptor degeneration. Am J Pathol. (2011) 178:2416–23. doi: 10.1016/j.ajpath.2011.01.013
30. Olson JL, Courtney RJ, Rouhani B, Mandava N, Dinarello CA. Intravitreal anakinra inhibits choroidal neovascular membrane growth in a rat model. Ocul Immunol Inflamm. (2009) 17:195–200. doi: 10.1080/09273940802710705
31. Celkova L, Doyle SL, Campbell M. NLRP3 inflammasome and pathobiology in AMD. J Clin Med. (2015) 4:172–92. doi: 10.3390/jcm4010172
32. Dana MR, Zhu SN, Yamada J. Topical modulation of interleukin-1 activity in corneal neovascularization. Cornea. (1998) 17:403–9. doi: 10.1097/00003226-199807000-00011
33. Jing R, Qi T, Wen C, Yue J, Wang G, Pei C, et al. Interleukin-2 induces extracellular matrix synthesis and TGF-β2 expression in retinal pigment epithelial cells. Dev Growth Differ. (2019) 61:410–8. doi: 10.1111/dgd.12630
34. Sato K, Takeda A, Hasegawa E, Jo YJ, Arima M, Oshima Y, et al. Interleukin-6 plays a crucial role in the development of subretinal fibrosis in a mouse model. Immunol Med. (2018) 41:23–9. doi: 10.1080/09114300.2018.1451609
35. Qin T, Gao S. Inhibition of proteasome activity upregulates IL-6 expression in RPE cells through the activation of P38 MAPKs. J Ophthalmol. (2018) 2018:5392432. doi: 10.1155/2018/5392432
36. Izumi-Nagai K, Nagai N, Ozawa Y, Mihara M, Ohsugi Y, Kurihara T, et al. Interleukin-6 receptor-mediated activation of signal transducer and activator of transcription-3 (STAT3) promotes choroidal neovascularization. Am J Pathol. (2007) 170:2149–58. doi: 10.2353/ajpath.2007.061018
37. Ramos de Carvalho JE, Verwoert MT, Vogels IMC, Reits EA, Van Noorden CJF, Klaassen I, et al. Involvement of the ubiquitin-proteasome system in the expression of extracellular matrix genes in retinal pigment epithelial cells. Biochem Biophys Rep. (2018) 13:83–92. doi: 10.1016/j.bbrep.2018.01.005
38. Hata M, Andriessen E, Hata M, Diaz-Marin R, Fournier F, Crespo-Garcia S, et al. Past history of obesity triggers persistent epigenetic changes in innate immunity and exacerbates neuroinflammation. Science. (2023) 379:45–62. doi: 10.1126/science.abj8894
39. Liu TT, Yang H, Zhuo FF, Yang Z, Zhao MM, Guo Q, et al. Atypical E3 ligase ZFP91 promotes small-molecule-induced E2F2 transcription factor degradation for cancer therapy. EBioMedicine. (2022) 86:104353. doi: 10.1016/j.ebiom.2022.104353
40. Hromas R, Orazi A, Neiman RS, Maki R, Van Beveran C, Moore J, et al. Hematopoietic lineage- and stage-restricted expression of the ETS oncogene family member PU.1. Blood. (1993) 82:2998–3004. doi: 10.1182/blood.V82.10.2998.2998
41. Huang KL, Marcora E, Pimenova AA, Di Narzo AF, Kapoor M, Jin SC, et al. A common haplotype lowers PU.1 expression in myeloid cells and delays onset of Alzheimer's disease. Nat Neurosci. (2017) 20:1052–61. doi: 10.1038/nn.4587
42. Liu Q, Yu J, Wang L, Tang Y, Zhou Q, Ji S, et al. Inhibition of PU.1 ameliorates metabolic dysfunction and non-alcoholic steatohepatitis. J Hepatol. (2020) 73:361–70. doi: 10.1016/j.jhep.2020.02.025
43. Karlstetter M, Ebert S, Langmann T. Microglia in the healthy and degenerating retina: insights from novel mouse models. Immunobiology. (2010) 215:685–91. doi: 10.1016/j.imbio.2010.05.010
44. Umazume A, Kezuka T, Matsuda R, Usui Y, Takahashi H, Yamakawa N, et al. Role of PU.1 expression as an inflammatory marker in experimental autoimmune uveoretinitis. Ocul Immunol Inflamm. (2018) 26:951–63. doi: 10.1080/09273948.2017.1299867
45. Zhou N, Liu K, Sun Y, Cao Y, Yang J. Transcriptional mechanism of IRF8 and PU.1 governs microglial activation in neurodegenerative condition. Protein Cell. (2019) 10:87–103. doi: 10.1007/s13238-018-0599-3
46. Sun J, Chen J, Li T, Huang P, Li J, Shen M, et al. ROS production and mitochondrial dysfunction driven by PU.1-regulated NOX4-p22(phox) activation in Aβ-induced retinal pigment epithelial cell injury. Theranostics. (2020) 10:11637–55. doi: 10.7150/thno.48064
47. Xu Y, Cui K, Li J, Tang X, Lin J, Lu X, et al. Melatonin attenuates choroidal neovascularization by regulating macrophage/microglia polarization via inhibition of RhoA/ROCK signaling pathway. J Pineal Res. (2020) 69:e12660. doi: 10.1111/jpi.12660
48. Zhu Y, Tan W, Demetriades AM, Cai Y, Gao Y, Sui A, et al. Interleukin-17A neutralization alleviated ocular neovascularization by promoting M2 and mitigating M1 macrophage polarization. Immunology. (2016) 147:414–28. doi: 10.1111/imm.12571
49. Wooff Y, Man SM, Aggio-Bruce R, Natoli R, Fernando N. IL-1 family members mediate cell death, inflammation and angiogenesis in retinal degenerative diseases. Front Immunol. (2019) 10:1618. doi: 10.3389/fimmu.2019.01618
Keywords: SPI1, transcription factor, AMD, CNV, macrophage polarization, innate immunity
Citation: Qi S, Zhang Y, Kong L, Bi D, Kong H, Zhang S and Zhao C (2024) SPI1-mediated macrophage polarization aggravates age-related macular degeneration. Front. Immunol. 15:1421012. doi: 10.3389/fimmu.2024.1421012
Received: 21 April 2024; Accepted: 28 May 2024;
Published: 24 June 2024.
Edited by:
Chaofeng Han, Second Military Medical University, ChinaReviewed by:
Juan Gao, Shanghai University, ChinaHui Liu, Case Western Reserve University, United States
Copyright © 2024 Qi, Zhang, Kong, Bi, Kong, Zhang and Zhao. 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: Chen Zhao, ZHJfemhhb2NoZW5AZnVkYW4uZWR1LmNu; Shujie Zhang, c2h1amllemhhbmdAZnVkYW4uZWR1LmNu
†These authors have contributed equally to this work