Background Air pollution is a major environmental factor threatening respiratory health. Different pollutants exhibit varying degrees of lag effects on respiratory diseases, and synergistic effects may exist among multiple pollutants. There is an urgent need to identify the key air pollutants influencing respiratory diseases and their interactive effects at specific lags. Objective To identify key pollutants affecting hospital admissions for respiratory diseases, to analyze their lag effect characteristics, and to quantify the impact of multi-pollutant synergistic effects on respiratory disease admissions. Methods Daily air pollution data, meteorological data, and respiratory disease hospitalization records were collected from multiple national monitoring stations in Hefei City from 2019 to 2024. A two-stage analytical framework was employed. First, a distributed lag model (DLM) was used to construct pollutant lag matrices, followed by least absolute shrinkage and selection operator (LASSO) regression to select key variables among fine particulate matter (PM_2.5), sulfur dioxide (SO₂), nitrogen dioxide (NO₂), carbon monoxide (CO), and ozone (O₃). Second, a generalized additive model (GAM) was established, incorporating product interaction terms and excess relative risk (ERI) to quantitatively assess synergistic effects among the selected pollutants. Results Through LASSO regression, 24 pollutant lag terms with non-zero coefficients were identified, among which NO₂, PM_2.5, and SO₂ accounted for 66.7% of the total positive effects and exhibited distinct lag patterns. Exposure to NO₂ showed acute risk, with a relative risk of 1.040 (95%CI: 1.023, 1.057) at lag0. Conversely, PM_2.5 and SO₂ exhibited delayed effects, with peak impacts observed at lag7 (RR=1.012, 95%CI: 1.0002, 1.024) and lag3 (RR=1.015, 95%CI: 1.0004, 1.03), respectively. Short-term NO₂ exposure (lag0–1) exhibited a significant positive interaction with medium-term PM_2.5 exposure (lag4–7). When both pollutants coexisted at high exposure levels, the combined health risk exceeded the sum of their independent effects by 4.40%. Stratified analysis showed that the interaction effects of pollutants were higher in females (ERI=5.00%), children (<18 years, ERI=5.50%), and older adults (≥65 years, ERI=3.90%)Conclusion There are significant lag heterogeneity and synergistic mechanisms in the health effects of air pollution in Hefei. It is recommended to implement real-time early warning responses for NO₂, prioritize medium-to-long-term management for PM_2.5 and SO₂, and adopt coordinated emission reduction measures when NO₂ and PM_2.5 concentrations rise simultaneously. Concurrently, enhanced health protection should be provided for sensitive populations such as women, children, and the elderly.

Glycogen synthase 1 (GYS1) is the rate-limiting enzyme responsible for glycogen synthesis in skeletal muscle, heart, brain, and other extrahepatic tissues, playing a central role in systemic energy homeostasis. The human GYS1 gene maps to chromosome 19q13.33, comprises 16 exons, and encodes a 737-amino-acid polypeptide that is highly conserved across mammals. GYS1 activity is subject to multilayered and precisely coordinated regulation. At the transcriptional level, the GYS1 promoter contains a hypoxia response element (HRE) that mediates HIF-1α-dependent induction under low-oxygen conditions, as well as a muscle-specific enhancer harboring MEF2 and MyoD binding sites that confers tissue-restricted expression. At the post-translational level, a hierarchical phosphorylation cascade serves as the primary activity switch: glycogen synthase kinase 3β (GSK3β) sequentially phosphorylates four C-terminal serine residues following casein kinase II priming, while protein kinase A (PKA) and AMP-activated protein kinase (AMPK) provide parallel inhibitory inputs at both N- and C-terminal sites. Dephosphorylation and reactivation are mediated by protein phosphatase 1 (PP1) through tissue-specific glycogen-targeting regulatory subunits such as PPP1R3A and PPP1R3B, which anchor PP1 to glycogen particles and direct its activity toward GYS1. The allosteric activator glucose-6-phosphate (G6P) binds at the dimer interface, simultaneously enhancing catalytic efficiency and promoting dephosphorylation susceptibility, thereby establishing a feed-forward activation loop that couples substrate availability to glycogen synthesis. Beyond phosphorylation, GYS1 is regulated by ubiquitination (mediated by the E3 ligase PJA1), acetylation, O-linked β-N-acetylglucosamine (O-GlcNAc) modification, and SUMOylation, which collectively modulate protein stability, subcellular localization, and protein-protein interactions. Epigenetic mechanisms, including CpG island methylation and histone acetylation dynamics, govern chromatin accessibility at the GYS1 locus, while muscle-specific microRNAs such as miR-1 and miR-206 fine-tune GYS1 expression at the post-transcriptional level. Dysregulation of GYS1 has been identified as a central pathogenic driver in a spectrum of human diseases. In inherited glycogen storage disorders—including Lafora disease, adult polyglucosan body disease (APBD), and Pompe disease—loss of upstream regulatory control leads to GYS1 hyperactivation and the accumulation of structurally abnormal or excessive glycogen, resulting in progressive neurodegeneration, myopathy, and multiorgan dysfunction. In type 2 diabetes mellitus (T2DM), impaired insulin signaling through the PI3K-AKT-GSK3β axis maintains GYS1 in a hyperphosphorylated inactive state in skeletal muscle, compromising postprandial glucose disposal and exacerbating hyperglycemia. In oncology, GYS1 exhibits context-dependent roles across multiple cancer types. In hepatocellular carcinoma, FMO2⁺ cancer-associated fibroblasts stabilize GYS1 by competitively inhibiting PJA1-mediated ubiquitination, and stabilized GYS1 subsequently activates NF‑κB/CCL19 signaling to promote tertiary lymphoid structure formation and enhance anti-PD-1 immunotherapy responsiveness. In clear cell renal cell carcinoma, GYS1 promotes tumor progression through non-canonical NF‑κB pathway activation via the scaffold protein RPS27A. In triple-negative breast cancer, GYS1 has been identified as a trigger of disulfidptosis and an activator of NF-κB signaling through non-enzymatic facilitation of IκBα degradation. In colorectal cancer, mitochondrial fission deficiency drives AMPK-dependent GYS1 upregulation and glycogen accumulation as a compensatory survival mechanism, while in cervical cancer, GYS1-maintained glycogen reserves fuel the pentose phosphate pathway to generate NADPH for ROS clearance, thereby conferring cisplatin resistance in cancer stem cells. Therapeutic strategies targeting GYS1 have gained substantial momentum across these disease contexts. For glycogen storage disorders, antisense oligonucleotides, small interfering RNAs (e.g., ABX1100), and small-molecule inhibitors (e.g., MZ-101) have demonstrated preclinical and early clinical efficacy in reducing pathological glycogen accumulation. For T2DM, pharmacological activation of GYS1 through GSK3β inhibition or enhancement of PP1-mediated dephosphorylation is being explored to restore insulin-stimulated glycogen synthesis. In cancer, GYS1-directed interventions—including targeted silencing to sensitize tumors to chemotherapy and immune microenvironment modulation to enhance immunotherapy—represent emerging precision oncology approaches. This review provides a comprehensive and integrated account of GYS1 gene structure, tissue-specific distribution, regulatory networks, and pathogenic roles in metabolic disorders and malignancies, with the aim of establishing a theoretical framework for the development of GYS1-targeted precision therapies.

Diabetic retinopathy (DR) is one of the most prevalent and vision-threatening microvascular complications of diabetes mellitus, yet its pathogenesis extends far beyond vascular injury alone. As the retina is among the most energy-demanding tissues in the body, its neurons, glial cells, pigment epithelial cells, pericytes, and endothelial cells are highly dependent on mitochondrial oxidative phosphorylation to maintain visual signal transduction, ionic homeostasis, and neurovascular integrity. This review summarizes current evidence indicating that mitochondrial dysfunction is not merely a downstream consequence of chronic hyperglycemia, but a central pathogenic hub that initiates, amplifies, and perpetuates retinal neurovascular degeneration in DR. Persistent hyperglycemia activates multiple abnormal metabolic pathways, including the polyol pathway, hexosamine pathway, protein kinase C signaling, advanced glycation end-product formation, and angiotensin II-related responses. Although these pathways differ mechanistically, they converge on excessive reactive oxygen species (ROS) generation, antioxidant depletion, and mitochondrial injury. Under diabetic stress, electron transport chain overload promotes mitochondrial ROS leakage, damages mitochondrial DNA, disrupts membrane potential, and impairs the transcription of key respiratory chain components. In parallel, mitochondrial quality-control systems become progressively compromised. The balance between fusion and fission shifts toward pathological fragmentation through reduced MFN1/2 and OPA1 activity and enhanced DRP1-mediated fission. Mitochondrial biogenesis is suppressed through inhibition of the AMPK/SIRT1/PGC-1α/NRF1/TFAM axis, while mitophagy changes from an early compensatory response to a later state of autophagic flux blockade and accumulation of dysfunctional mitochondria. Importantly, damaged mitochondria serve as signal amplifiers linking metabolic stress to inflammation and programmed cell death. Mitochondrial ROS, oxidized mitochondrial DNA, calcium overload, cardiolipin exposure, and membrane permeabilization activate interrelated death pathways, including intrinsic apoptosis, ferroptosis, and pyroptosis. Cytochrome C and apoptosis-inducing factor promote caspase-dependent and caspase-independent apoptosis; iron dyshomeostasis, glutathione depletion, GPX4 dysfunction, and lipid peroxidation drive ferroptosis; and mitochondrial danger signals activate the NLRP3 inflammasome and gasdermin-dependent pyroptosis. These pathways jointly damage the retinal neurovascular unit and contribute to pericyte loss, endothelial barrier breakdown, Müller cell dysfunction, retinal ganglion cell apoptosis, retinal pigment epithelial injury, and photoreceptor degeneration. This review also emphasizes the role of epigenetic regulation in stabilizing mitochondrial pathology. DNA methylation, histone modifications, and non-coding RNAs interact to silence mitochondrial protective genes, alter antioxidant responses, and maintain the “metabolic memory” of DR even after glycemic normalization. Therefore, mitochondrial dysfunction should be understood as a dynamic, multidimensional network rather than a single pathological event. Current clinical approaches, such as laser photocoagulation, intravitreal anti-VEGF therapy, and vitrectomy, mainly target advanced vascular lesions and are limited by invasiveness, incomplete responsiveness, recurrence, and potential adverse effects. Therapeutically, strategies targeting mitochondrial ROS, restoring mitochondrial dynamics, enhancing biogenesis, regulating mitophagy, inhibiting inflammasome activation, correcting epigenetic abnormalities, and improving targeted delivery systems show promising potential. However, major translational barriers remain, including retinal cell heterogeneity, stage-specific mitochondrial responses, insufficient organelle-specific drug delivery, and long-term safety concerns. A deeper understanding of mitochondrial regulatory networks may support earlier, more precise, and multi-target interventions for preventing or slowing DR progression.

Owing to its inherent advantages—such as being non-destructive, rapid, highly molecule-specific, and minimally interfered with by moisture—Raman spectroscopy has been widely adopted in the fields of skin barrier function assessment, monitoring the transdermal penetration of active cosmetic ingredients, and the identification and quality control of cosmetic products. Despite these strengths, the practical application of this technique faces a significant bottleneck: the strong fluorescence background generated by endogenous skin components and exogenous cosmetic additives. Endogenous skin substances, such as structural proteins (e.g., collagen and elastin), metabolic coenzymes (e.g., nicotinamide adenine dinucleotide), and pigments (e.g., melanin), together with exogenous cosmetic constituents like organic colorants, chemical sunscreens, and fragrances, often possess strong absorption and emission characteristics. When excited by lasers, these components produce a fluorescence background that can be 10⁶ to 10⁸ times stronger than the Raman scattering signals, effectively masking the inherently weak vibrational fingerprint information. In recent years, driven by the rapid development of optoelectronic hardware and artificial intelligence algorithms, fluorescence suppression strategies have evolved from isolated, single-method approaches into comprehensive, multi-level synergistic systems. These systems are categorized into three distinct tiers: sample preparation, signal acquisition, and data processing. At the sample preparation level, techniques such as photobleaching and surface-enhanced Raman spectroscopy (SERS) are employed to eliminate or bypass the generation of fluorescence at the source. At the signal acquisition level, instrumental improvements—including the use of long-wavelength near-infrared excitation (typically 785 nm or 1 064 nm), confocal spatial filtering, and shifted excitation Raman difference spectroscopy (SERDS)— are utilized to physically isolate Raman signals from the fluorescence background. Furthermore, at the data processing level, numerical baseline correction methods such as polynomial fitting, penalized least squares (e.g., airPLS, arPLS), wavelet transform, and derivative algorithms are increasingly integrated into the analytical pipeline to extract Raman spectral features from mixed signals without increasing hardware costs or acquisition time. This review provides a systematic categorization and critical evaluation of these fluorescence suppression methods, detailing their underlying principles, technical advantages, and inherent limitations in diverse experimental setups. By focusing on critical application scenarios—including skin barrier assessment, percutaneous absorption monitoring, the routine quality control of cosmetics, and the emerging field of portable on-site detection—this paper explores the current state of technique selection and optimization. Finally, the article discusses future development trends, emphasizing the necessity of constructing adaptive, tiered suppression strategies, developing intelligent and automated data processing algorithms, and promoting the integration of portable, multi-modal diagnostic devices. The objective of this review is to provide a comprehensive technical reference to facilitate the transition of Raman spectroscopy from a specialized laboratory tool into a routine, robust analytical platform for advancements in skin science and cosmetic research.

ObjectiveTo investigate the relationship between periodontitis severity and brain functional indices-amplitude of low-frequency fluctuations （ALFF） and fractional ALFF （fALFF） in cognitively normal middle-aged and older adults， and to examine their associations with cognitive performance. Methods42 cognitively normal middle-aged and older participants were enrolled. Periodontal data—including the number of remaining teeth （NTP）， bleeding on probing positivity rate （BOP%）， probing depth （PD）， and clinical attachment loss （CAL） as well as resting-state functional magnetic resonance imaging data were collected. Participants were divided into a mild periodontitis group and a moderate-to-severe periodontitis group according to periodontitis severity. Group differences in ALFF and fALFF were compared， and correlation analyses were performed. ResultsCompared with the mild periodontitis group， the moderate-to-severe group showed decreased ALFF values in the right lingual gyrus， superior occipital gyrus， paracentral lobule， and in the left calcarine cortex as well as the precentral and postcentral gyri. Decreased fALFF values were observed in the right lingual gyrus， cuneus， and left calcarine cortex. Indicators related to periodontitis progression ［BOP%， PD=4-5 mm （%）， PD≥6 mm （%）， CAL=3-4 mm （%） and CAL≥5 mm （%）］ were negatively correlated with ALFF， fALFF and cognitive performance， whereas indicators reflecting better periodontal function ［PD=1-3 mm （%） and CAL=0-2 mm （%）］ were positively correlated with these measures. ConclusionIn cognitively normal middle-aged and older adults， periodontitis severity is associated with alterations in brain functional activity and cognitive performance.

Objective: To investigate the incidence and influencing factors of new onset suicidal behaviors among junior and senior high school students, so as to provide evidence for early warning and intervention. Methods: Using a cluster sampling method, baseline data were collected from September to November 2023 among 2 379 adolescents selected from 2 junior high schools, 2 general senior high schools and 2 vocational senior high schools in Wuhan, Hubei Province, with follow up carried out from September to November 2024. Data were collected using the School Adaptation Scale, the Essential Resilience Scale, the Personal Social Support Scale, the Adverse Childhood Experiences Questionnaire, the Adolescent Life Events Scale, the Pittsburgh Sleep Quality Index, the Center for Epidemiologic Studies Depression Scale, and the Columbia-Suicide Severity Rating Scale. Multivariate Logistic regression analysis was used to identify factors influencing new onset suicidal behavior among junior and senior high school students. Results: Among junior and senior high school students without suicide risk at baseline, 53.55% remained stable risk free during follow up, while 13.28% developed new suicidal behavior during follow up, of which suicidal ideation, suicide plan and suicide attempts accounted for 11.69%, 8.83% and 6.52%, respectively. Multivariate Logistic regression revealed that increased depressive symptoms at follow up elevated the risk of overall new onset suicidal behaviors ( OR =4.14), ideation ( OR = 4.59), plans ( OR =4.24), and attempts ( OR =4.09); poor sleep quality also increased the risk of overall new onset suicidal behaviors ( OR =2.23), plans ( OR =2.03), and attempts ( OR =2.65) (all P <0.05). Non boarding junior and senior high school students had lower risks of overall new onset suicidal behaviors and ideation ( OR=0.70 and 0.49, respectively); secondary school students with academic performance ranking in the top 20%-30% of their class had a lower risk of new onset ideation ( OR =0.27); and those whose place of origin was outside the city had a lower risk of new onset plans ( OR =0.49) (all P <0.05). Conclusions New onset suicidal behaviors among junior and senior high school students exhibit dynamic and stage specific differences. Depressive symptoms serve as a common core risk factor across all stages, whereas sleep quality, academic performance, boarding status, and place of origin exert differential influences at distinct stages.

ObjectiveStroke remains the leading cause of acquired disability in adults worldwide, imposing a substantial socioeconomic burden on patients, families, and healthcare systems. Among the various neurological deficits caused by stroke, motor dysfunction——particularly contralateral hemiplegia——is the most common and debilitating consequence, and it predominantly results from damage to the cortical motor conduction pathways. When motor cortex neurons or their projecting fibers are damaged, the distal axons and their associated myelin sheaths undergo a series of progressive, programmed degeneration, fragmentation, and eventual phagocytic clearance. This pathological process is termed Wallerian degeneration (WD). However, the temporal progression of WD following ischemic stroke and the associated alterations in the local glial microenvironment remain incompletely understood. The present study was designed to systematically elucidate the time course of corticospinal tract (CST) axonal WD after ischemic stroke in mice, and to simultaneously characterize the dynamic pathological changes of astrocytes, microglia, oligodendrocyte lineage cells, and myelinating oligodendrocytes during the WD process. MethodsEight-week-old male C57BL/6J mice were randomly assigned to control and ischemic stroke groups. A focal cortical ischemic stroke model was established in the right sensorimotor cortex using the photothrombotic method. This method induces localized intravascular thrombosis in the target cortical region through precisely controlled interaction between the photosensitizer and laser light of specific wavelength. At day 1 post‑surgery, 2,3,5-triphenyltetrazolium chloride (TTC) staining was performed to evaluate infarct volume, hematoxylin and eosin (HE) staining for histopathological examination, and Nissl staining for assessment of neuronal injury. For axonal degeneration time-course analysis, rAAV-CAG-mCherry anterograde viral tracing was employed to specifically label CST neurons. Cervical spinal cord sections at the CST region were collected at day 3, 5, 7, 14, and 60 post-injury. Morphological changes of mCherry-positive CST fibers were observed under fluorescence microscopy, and the fluorescence signal area fraction within the defined region was quantitatively analyzed to evaluate the dynamic process of axonal degeneration. For glial response assessment, immunofluorescence staining against Iba1 (microglia), GFAP (astrocytes), Olig2 (oligodendrocyte lineage cells), and APC (myelinating oligodendrocytes) was performed at day 7, 14, and 60 post‑surgery. The numbers of positive cells were counted to quantitatively analyze the dynamic changes in each glial population at different time points after ischemia, thereby revealing the evolution of the glial microenvironment during WD. ResultsThe photothrombotic method stably and reliably established an ischemic brain injury model confined to the sensorimotor cortex. Viral tracing revealed that the proportion of mCherry‑positive CST fibers was significantly reduced at day 7 after ischemic injury to the neuronal soma. Immunofluorescence staining demonstrated that during WD, the numbers of microglia and astrocytes in the CST region were significantly increased, whereas the numbers of oligodendrocyte lineage cells and myelinating oligodendrocytes showed no significant changes. ConclusionWD of CST axons is observed at day 7 post-ischemic stroke in a delayed manner. This degenerative process selectively activates microglia and astrocytes, while leaving the number of oligodendrocyte lineage cells largely unaffected.

ObjectiveTo investigate the impact of hepatocellular carcinoma （HCC） on the prognosis of patients with liver cirrhosis undergoing emergency endoscopic therapy for esophagogastric variceal bleeding， as well as independent influencing factors for the prognosis of liver cirrhosis patients without HCC after emergency endoscopic therapy for esophagogastric variceal bleeding. MethodsA total of 117 liver cirrhosis patients without HCC and 119 liver cirrhosis patients with HCC who underwent emergency endoscopic therapy for esophagogastric variceal bleeding in Renmin Hospital of Wuhan University from January 2017 to July 2023 were enrolled. Basic information including age and sex was collected from all patients， as well as the presence or absence of chronic diseases such as hypertension， diabetes， and coronary heart disease， the time of emergency endoscopy after admission， and liver function parameters including international normalized ratio， albumin， creatinine， sodium， total bilirubin， alanine aminotransferase， and aspartate aminotransferase （AST）. The independent-samples t test was used for comparison of normally distributed continuous variables between two groups， and the Wilcoxon rank-sum test was used for comparison of non-normally distributed continuous variables between two groups； the chi-square test was used for comparison of categorical variables between groups. The covariance analysis and the multivariate logistic regression analysis were used for comparison of outcome variables after control of baseline variables， and the Kaplan-Meier survival curve was plotted for each group. The univariate and multivariate Cox regression analyses were performed for survival time in the non-HCC group to investigate the independent influencing factors for survival time， and then the Kaplan-Meier curve analysis and the log-rank test were performed to validate such independent influencing factors and analyze the independent influencing factors for secondary outcomes. ResultsCompared with the non-HCC group， the HCC group had significantly higher red blood cell transfusion units （6.00［2.00～9.00］ vs 4.00［1.75～7.00］， Z=-2.050， P=0.040， F=4.869， adjusted P=0.028）， a significantly shorter survival time （29.77±16.01 days vs 38.07±11.43 days， t=4.574， P<0.001， F=17.294， adjusted P<0.001）， and a significantly higher 5-day rebleeding rate （22.69% vs 6.84%， χ²=11.736， P<0.001， adjusted P=0.021）. The Kaplan-Meier curve analysis showed that the risk of 42-day mortality in the HCC group was 3.897 （95% confidence interval ［CI］： 2.338 ‍— 6.495， P<0.001） times that in the non-HCC group. The multivariate Cox regression analysis of the non-HCC group showed that the total length of hospital stay （hazard ratio ［HR］=0.793， 95%CI： 0.644 ‍— ‍0.976， P=0.029） was an independent protective factor for 42-day survival. The Kaplan-Meier curve analysis showed that a length of hospital stay of >9 days was beneficial for the prognosis of patients （HR=4.302， 95%CI： 1.439 — 12.870， P=0.037）. Blood sodium level （odds ratio ［OR］=0.523， 95%CI： 0.289 ‍— ‍0.945， P=0.032） and MELD-Na score （OR=0.495， 95%CI： 0.257 ‍— ‍0.954， P=0.036） were independent protective factors against 5-day rebleeding， while AST level was an independent risk factor for 5-day rebleeding （OR=1.023， 95%CI： 1.002 ‍— ‍1.043， P=0.028） and in-hospital death （OR=1.036， 95%CI： 1.001‍— ‍1.073， P=0.045）. ConclusionLiver cirrhosis patients with variceal bleeding and HCC tend to have a worse prognosis， and for the non-HCC group， in-hospital mortality rate increases with the increase in AST level. The total length of hospital stay is an independent protective factor for survival time in the non-HCC group， and it is recommended to appropriately prolong the length of hospital stay for such patients.

Objective To summarize the clinical and genetic characteristics of end-stage renal disease caused by PLCE1 gene mutations.Methods A retrospective analysis of the clinical and genetic features of three children from a family with PLCE1 gene mutations was conducted,along with a literature review of hereditary kidney disease cases caused by PLCE1 gene mutations.Results The proband was an 8-year-old male presenting with nephrotic syndrome stage 4 chronic kidney disease.Renal biopsy showed focal segmental glomerulosclerosis.Two years and five months after kidney transplantation,the patient had persistent negative proteinuria and normal renal function.Whole-exome sequencing identified two pathogenic heterozygous variants:c.961C>T and c.3255_3256delinsT,with c.3255_3256delinsT being a novel mutation.Family screening revealed no renal involvement in the parents,but among five siblings,one brother died at age of 4 years from end-stage renal disease.A 7-year-old sister presented with proteinuria and bilateral medullary sponge kidney,with proteinuria resolving after one year of follow-up.A 3-year-old brother died after kidney transplantation due to severe pneumonia.The literature review included 45 patients with hereditary kidney disease caused by PLCE1 gene mutations.The main clinical phenotype was nephrotic syndrome(87%,39/45),and renal pathology predominantly showed focal segmental glomerulosclerosis(57%,16/28).No mutation hotspots were identified.Conclusions Compound heterozygous mutations in the PLCE1 gene can lead to rapid progression of the disease to end-stage renal disease,with favorable outcomes following kidney transplantation.Family screening is crucial for early diagnosis,and medullary sponge kidney may be a novel phenotype associated with these gene mutations.Citaion:[Chinese Journal of Contemporary Pediatrics,2025,27(5):580-587]

Combined with elastic network model(ENM),the perturbation response scanning(PRS)has emerged as a robust technique for pinpointing allosteric interactions within proteins.Here,we proposed the PRS analysis of drug-target networks(DTNs),which could provide a promising avenue in network medicine.We demonstrated the utility of the method by introducing a deep learning and network perturbation-based framework,for drug repurposing of multiple sclerosis(MS).First,the MS comorbidity network was constructed by performing a random walk with restart algorithm based on shared genes between MS and other diseases as seed nodes.Then,based on topological analysis and functional annotation,the neurotransmission module was identified as the"therapeutic module"of MS.Further,perturbation scores of drugs on the module were calculated by constructing the DTN and introducing the PRS analysis,giving a list of repurposable drugs for MS.Mechanism of action analysis both at pathway and structural levels screened dihydroergocristine as a candidate drug of MS by targeting a serotonin receptor of se-rotonin 2B receptor(HTR2B).Finally,we established a cuprizone-induced chronic mouse model to evaluate the alteration of HTR2B in mouse brain regions and observed that HTR2B was significantly reduced in the cuprizone-induced mouse cortex.These findings proved that the network perturbation modeling is a promising avenue for drug repurposing of MS.As a useful systematic method,our approach can also be used to discover the new molecular mechanism and provide effective candidate drugs for other complex diseases.