Objective To evaluate the effectiveness of schistosomiasis surveillance in Jiangsu Province during the stage moving from transmission control to transmission interruption, and to analyze the current risk and challenges, so as to provide the evidence for achieving the target of schistosomiasis elimination. Methods Schistosomiasis surveillance data were collected from Jiangsu Province from 2012 to 2024, and the endemic areas, Schistosoma japonicum infections in humans and livestock, Oncomelania hupensis snail distribution and implementation of integrated interventions were descriptively analyzed. In addition, the trends in areas with snails, seroprevalence of human S. japonicum infections and numbers of advanced schistosomiasis cases were assessed using a Joinpoint regression model. Results The endemic areas of schistosomiasis continued to shrink in Jiangsu Province from 2012 to 2024, with the number of schistosomiasis-eliminated counties (cities, districts) increasing from 53 (75.71%) to 63 (96.92%), and interruption of schistosomiasis transmission was achieved across the province. A total of 4 600 300 person-times were tested for serum antibodies against S. japonicum, with 28 719 person-times positive detected; and 616 500 person-times were tested S. japonicum infections among local residents in Jiangsu Province from 2012 to 2024, with only 3 egg-positives detected, and no egg-positives found since 2017. A total of 187 600 herd-times were tested for schistosomiasis in livestock, and no S. japonicum infections were found. O. hupensis snail survey was performed covering 1 018 408.97 hm², and a total of 35 556.35 hm² was found with snail-infested habitats, including 174.40 hm2 of emerging snail-infested habitats. A total of 1 102 800 O. hupensis snails were identified for S. japonicum infections, and no infections were found. The areas of snail-infested habitats appeared a tendency towards a rise in Jiangsu Province from 2019 to 2023 (APC = 23.67%, P < 0.05), and the actual areas of snail-infested habitats appeared a tendency towards a decline from 2012 to 2015 (APC = −22.77%, P < 0.05), and towards a rise from 2015 to 2023 (APC = 9.76%, P < 0.01). The seroprevalence of anti-S. japonicum antibodies appeared a tendency towards a decline among residents in Jiangsu Province from 2017 to 2023 (APC = −14.92%, P < 0.01). In addition, the number of newly diagnosed advanced schistosomiasis cases appeared a tendency towards a decline from 2012 to 2024 (APC = −12.02%, P < 0.01), and the numbers of advanced schistosomiasis patients requiring treatment showed a tendency towards a decline from 2012 to 2021 (APC = −10.56%, P < 0.01) and from 2021 to 2023 (APC = −20.06%, P < 0.01). Conclusions Great progresses had been achieved in schistosomiasis control in Jiangsu Province following transmission control, and transmission interruption had been achieved; however, there are still snail-infested habitats. High-intensity surveillance and integrated control are required to be maintained to advance the achievement of the target of schistosomiasis elimination in Jiangsu Province.

Objective To investigate the species of sandflies and the prevalence of Leishmania infections in sandflies from selected areas of northern and northwestern China, so as to provide insights into identification of leishmaniasis vectors and assessment of epidemiological trends of leishmaniasis in China. Methods Sandfly samples were collected from Mentougou District of Beijing Municipality, Xiangning County in Linfen City of Shanxi Province, Ejin Banner in Alxa League of Inner Mongolia Autonomous Region, and Payzawat County of Karamay District of Karamay City, Gaochang District of Turpan City in Xinjiang Uygur Autonomous Region from July 2023 to July 2024. Approximately 100 intact female sandfly samples were randomly selected from each site and the species of sandflies was identified according to morphological characteristics and molecular assays. Female sandflies originating from the same habitat were grouped into pools of 10 individuals. Leishmania infection was detected using polymerase chain reaction (PCR) assay targeting the internal transcribed spacer 1 (ITS-1) gene, and the prevalence of Leishmania infection was calculated in sandflies from different sampling sites using the minimum infection rate (MIR) method. In addition, positive amplicons were sequenced and subjected to phylogenetic analysis. Results A total of 6 155 sandflies were collected from different environments at sampling sites across the six aforementioned regions from July 2023 to July 2024. Phlebotomus chinensis (96.00%) was the dominant sandfly species in Mentougou District, Beijing Municipality, with a small proportion of Ph. sergenti (4.00%), and only Ph. chinensis was found in Xiangning County, Linfen City, Shanxi Province. Ph. wui was the only sandfly species detected in Ejin Banner, Alxa League, Inner Mongolia Autonomous Region, and Payzawat County, Kashgar City, Xinjiang Uygur Autonomous Region, and Ph. caucasicus (97.70%) was the dominant sandfly species in Karamay District, Karamay City, Xinjiang Uygur Autonomous Region, with a small proportion of Ph. wui (2.30%), while Ph. alexandri was the only species in Gaochang District, Turpan City, Xinjiang Uygur Autonomous Region. A total of 40, 60, 34, 18, 18, and 22 pools of sandfly samples were tested from Mentougou District in Beijing Municipality, Xiangning County in Linfen City of Shanxi Province, Ejin Banner in Alxa League of Inner Mongolia Autonomous Region, Payzawat County in Kashgar City, Karamay District in Karamay City, and Gaochang District in Turpan City of Xinjiang Uygur Autonomous Region, respectively. L. infantum was detected in Ph. chinensis samples from Mentougou District in Beijing Municipality, and Xiangning County of Linfen City in Shanxi Province, with MIR of 0.25% to 1.00%, and L. donovani was detected in Ph. wui from Ejin Banner in Alxa League of Inner Mongolia Autonomous Region, and Payzawat County in Kashgar City of Xinjiang Uygur Autonomous Region, with MIR of 0.56% to 0.88%; however, no Leishmania infection was detected in Ph. caucasicus from Karamay District in Karamay City or Ph. alexandri from Gaochang District in Turpan City of Xinjiang Uygur Autonomous Region. Phylogenetic analysis showed that the Leishmania ITS-1 gene sequences obtained from Mentougou District in Beijing Municipality and Xiangning County in Linfen City of Shanxi Province were clustered into the same clade with the reference sequences of L. infantum ITS-1 gene, while the Leishmania ITS-1 gene sequences obtained from Ejin Banner in Alxa League of Inner Mongolia Autonomous Region and Payzawat County in Kashgar City of Xinjiang Uygur Autonomous Region were clustered into the same clade with the reference sequences of L. donovani ITS-1 gene. Conclusions There are variations in sandfly species in selected areas of northern and northwestern China, and variations in the species of Leishmania infecting sandflies. Improved surveillance of sandfly vectors and targeted control strategies with adaptations to geographical features and leishmaniasis vectors are recommended.

Objective To examine the impact of extreme temperature and drought stress on the survival of Bulinus globosus, so as to provide the theoretical evidence for the genomic research of Bulinus in absence of reference genes. Methods B. globosus snail samples were collected from Kiwani Shehia in Pemba Island, Zanzibar, Tanzania, and offspring snails were obtained through laboratory breeding and reproduction. A total of 120 10-week-old B. globosus snails from the same generation were selected and randomly assigned into four groups, including the high-temperature drought (HD) group, normal temperature drought (D) group, low-temperature drought (LD) group, and the control (C) group, of 30 snails in each group. Snails in HD, D, and LD groups were placed in beakers containing dry soil at the bottom and subsequently housed in climate chambers at 35, 26 ℃, and 10 ℃, respectively, while snails in Group C were maintained in 500 mL petri dishes containing dechlorinated tap water at 26 ℃. Following 3 days of breeding, living snails in each group were collected, and soft tissues were dissected and isolated. Total RNA was extracted from snail soft tissues for library construction, followed by high-throughput sequencing on the Illumina HiSeq 4000 sequencing system. De novo transcriptome assembly was performed using the Trinity software, and the longest transcripts were selected as unigenes. Gene functional annotations of unigenes were conducted using the Diamond software against Gene Ontology (GO) knowledgebase, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database, NCBI non-redundant (NR) protein sequences database, Protein Family (Pfam) database, and UniProtKB/Swiss-Prot (Swiss-Prot) knowledgebase. GO and KEGG enrichment analyses of differentially expressed genes (DEGs) were performed using the topGO and clusterProfiler software, respectively. In addition, four relevant genes were selected for validation using a real-time quantitative PCR (qRT-PCR) assay to verify the reliability of transcriptome sequencing results. Results Following 3 days of breeding, there were 7, 20, 28, and 30 survival B. globosus snails in HD, LD, D, and C groups, with corresponding survival rates of 23.33% (7/30), 66.67% (20/30), 93.33% (28/30), and 100.00% (30/30), respectively (χ² = 52.72, P < 0.001). De novo transcriptome assembly generated 176 942 unigenes, with annotation rates of 0.98%, 13.49%, 26.46%, 12.48%, and 14.39% against GO knowledgebase, KEGG pathway database, NR protein sequences database, Pfam database, and Swiss-Prot knowledgebase, respectively. There were 33 up-regulated and 72 down-regulated genes in Group D, 483 up-regulated and 815 down-regulated genes in Group HD, and 245 up-regulated and 172 down-regulated genes in Group LD relative to in Group C. Following removal of overlapping genes across groups and unmatched genes, 11 candidate genes were identified. GO and KEGG analyses revealed 3 heat shock protein (HSP)-related DEGs in these 11 candidate genes, which were annotated as HSP12.2, HSP70, and HSP20 genes and were all significantly up-regulated in each treatment group. Three immune and nervous system-related DEGs were identified, and were all significantly down-regulated in each treatment group, which were involved in the neural cell adhesion molecule L1-like protein pathway, fibrinogen binding protein pathway, and leukocyte elastase inhibitor-like protein pathway. qRT-PCR assay quantified that the expression trends of four genes related to temperature and drought stress across different treatment groups were highly consistent with transcriptome sequencing data. Conclusion The survival rate of B. globosus significantly reduces under combined stresses of extreme temperature and drought, possibly due to an imbalance in its cellular homeostasis regulatory system.

Objective Objective To analyze the potential spatial factors affecting the genetic differentiation of Oncomelania hupensis robertsoni in Yunnan Province. Methods A total of 13 administrative villages were selected from schistosomiasis-endemic areas of Yunnan Province as O. hupensis snail sampling sites. At least 200 snails were collected in each site, and the spatial variable data of each site were recorded, including longitude, latitude and altitude. Thirty active and Schistosoma japonicum uninfected O. hupensis snails were selected from each sampling site by means of the crawling method and the cercarial shedding method. Genomic DNA was extracted from O. hupensis snails. Following PCR amplification, purification of PCR amplification products and sequencing, the gene sequences of O. hupensis snail samples were spliced and edited using the DNAstar software and the NCBI database to yield the complete mitochondrial sequences of O. hupensis snails at each sampling site, and the mitochondrial genetic distance matrix of O. hupensis robertsoni was calculated at each sampling site. The geographical coordinates of each sampling site were marked using the software ArcGIS 10.2, and the straight-line geographical distance between each sampling site was calculated. The altitude difference, longitude difference and latitude difference between each sampling site were calculated using the Excel software, and the correlation between the mitochondrial genetic distance matrix of O. hupensis robertsoni and each spatial variable matrix was examined by using the Mantel test at 13 sampling sites in Yunnan Province. Results Among the 13 O. hupensis snail sampling sites in Yunnan Province, the largest mitochondrial genetic distance of O. hupensis robertsoni snail populations was seen between Anding Village, Nanjian Yi Autonomous County and Caizhuang Village, Midu County (26.244 2), and the largest geographical distance was seen between Dongyuan Village, Gucheng District and Cangling Village, Chuxiong County (272.64 km). The highest altitude difference was seen between Anding Village, Nanjian Yi Autonomous County and Dongyuan Village, Gucheng District (1 086.10 m), and the largest longitude difference was found between Qiandian Village, Eryuan County and Cangling Village, Chuxiong County (1.86°), while the largest latitude difference was measured between Leqiu Village, Nanjian Yi Autonomous County and Dongyuan Village, Gucheng District (1.81°). In addition, the mitochondrial genetic distance of O. hupensis robertsoni snail populations was positively correlated with altitude at 13 snail sampling sites in Yunnan Province (r = 0.542 8, P < 0.001), and showed no significant correlations with geographical distance (r = 0.093 4, P > 0.05), longitude (r = −0.199 5, P > 0.05) or latitude (r = 0.205 7, P > 0.05). Conclusion Altitude may be a potential spatial factor affecting the genetic differentiation of O. hupensis robertsoni in Yunnan Province.

Background Meteorological factors are among the key extrinsic triggers for the onset and exacerbation of cardiovascular and cerebrovascular diseases (CVD). Against the backdrop of sustained global warming, elucidating the impact of ambient temperature and atmospheric pressure on CVD, especially on pre-hospital CVD emergent events, has become imperative for evidence-based prevention and emergency preparedness. Objective To quantify the temporal trends of daily mean temperature and atmospheric pressure and their associations with pre-hospital CVD emergent events in Yantai, and to explore effect modification by demographic subgroups and geographic areas, thereby providing an empirical basis for the rational allocation of emergency medical resources. Methods Pre-hospital CVD emergency data from January 1, 2016 to December 31, 2022 were selected from the Yantai 120 Emergency Medical Command System. Synchronous meteorological factors and environmental pollutant data were obtained from the websites of the National Oceanic and Atmospheric Administration and the National Centers for Environmental Information of the United States. Time-series analysis combined with distributed lag non-linear model was used to analyze the association between daily temperature, atmospheric pressure, and pre-hospital CVD emergencies. Average annual percentage changes (AAPC) were calculated using Joinpoint (version 5.2.0.0) to reflect temporal trends. Spearman correlation analysis was employed to screen variables with low collinearity for inclusion in the multi-pollutant adjusted models. Results From 2016 to 2022, a total of 268275 pre-hospital CVD emergency cases were recorded in Yantai City, showing an increasing annual trend (AAPC=7.16%, P=0.002). Emergency volume, daily temperature, atmospheric pressure, and PM₁₀ displayed marked seasonal patterns, with summer bottoms and winter peaks for emergencies, daily temperature higher in summer and lower in winter, and the inverse for atmospheric pressure. Correlation analysis showed that both daily temperature and atmospheric pressure were positively correlated with O₃, with correlation coefficients of 0.745 and 0.723, respectively; negatively correlated with PM_2.5, with correlation coefficients of −0.246 and −0.235, respectively; and negatively correlated with PM₁₀, with correlation coefficients of −0.241 and −0.229, respectively (P < 0.05). A U-shaped correlation was revealed between daily temperature and pre-hospital CVD emergency risk. Using 24 °C as the reference, risk was elevated at both low and high daily temperatures, with the strongest cold effect at –11 ℃ (RR=1.53, 95%CI: 1.35, 1.73) and strongest heat effect at 30 ℃ (RR=1.06, 95%CI: 1.01, 1.11); cold effects were prolonged (lag0–14d) compared to heat effects. Atmospheric pressure exhibited a positive association with emergencies; using a median of 1012.2 hPa as the reference, the maximum risk occurred at 1037 hPa (RR=1.25, 95%CI: 1.03, 1.52), with a 0–5 d lag. In the subgroup analysis, the impact of low temperature was more pronounced among females (RR_{−9 ℃}=1.66), people over 60 years (RR_{−9 ℃}=1.91), and rural residents (RR_{−9 ℃}=1.76). High atmospheric pressure significantly increased the risk of pre-hospital CVD emergency among females (RR_{1035 hPa}=1.54), people under 18 years (RR_{1035 hPa}=3.62), and rural residents (RR_{1035 hPa}=1.46). The sensitivity analysis results indicated that, after excluding the impact of pandemic and the degrees of freedom for time variations, the effects of daily average temperature and atmospheric pressure on pre-hospital CVD emergency volume were consistent with the primary analysis findings. Conclusion Both non-optimal temperature and atmospheric pressure significantly increase the pre-hospital CVD emergency volume; therefore, emergency medical services should develop targeted response protocols, optimize resource allocation, and intensify public health education to mitigate the pre-hospital burden attributable to temperature and pressure fluctuations.

ObjectiveTo delineate imaging findings using an imaging platform and investigate the correlation between MRI characteristics of spinal cord atrophy and clinical diagnosis in children with spinal cord injury (SCI). MethodsImaging data of 150 children with SCI admitted to Beijing Bo'ai Hospital, China Rehabilitation Research Center, from January, 2002 to March, 2024 were collected and imported into the imaging platform. The anteroposterior and transverse diameters of the middle part of the spinal cord at the cross-section with the most severe atrophy were measured, and the relevant indicators of the previous normal spinal cord segment were measured as controls; the radiomic features were extracted. Clinical data of the children including gender, age, cause of injury, sensory level, motor level, spinal cord injury level, injury severity and disease course were collected. ResultsSpinal cord atrophy was identified in 81 cases (54%), among which 78 cases (96%) were American Spinal Injury Association Impairment Scale (AIS) grade A and 3 cases (4%) were AIS grade C. The upper boundary of the spinal cord atrophy site strongly correlated with the injury level, motor level and sensory level (r > 0.8, P < 0.001). ConclusionMore than half of children with SCI may develop secondary spinal cord atrophy, the vast majority of whom suffer from complete spinal cord injury; the upper boundary of spinal cord atrophy is correlated with the injury level.

ObjectiveCleft palate (CP) is a common congenital deformity often associated with velopharyngeal insufficiency (VPI), which disrupts the physiological coupling between respiration and speech. Conventional clinical assessments, such as nasometry and spirometry, provide limited static data and fail to visualize the dynamic spatiotemporal distribution of lung ventilation during phonation. This study introduces spatiotemporal electrical impedance tomography (ST-EIT) to evaluate speech-respiratory functional features in CP patients compared to normal controls (NC). The aim is to characterize multi-domain respiratory patterns and to validate an interpretable machine learning framework for providing objective, quantitative evidence for clinical assessment. MethodsSeventy-five participants were enrolled in this study, comprising 37 patients with surgically repaired CP and 38 healthy volunteers matched for age, gender, and body mass index (BMI). All subjects performed standardized sustained phonation tasks while undergoing synchronous monitoring with a 16-electrode EIT system and a pneumotachograph. A comprehensive feature engineering pipeline was developed to extract physiological parameters across 3 complementary domains. (1) Temporal domain: including inspiratory/expiratory phase duration (tPhase), time constants (Tau), and inspiratory-to-expiratory time ratios (TI/TE); (2) airflow domain: comprising mean flow, peak flow, and instantaneous flow at 25%, 50%, and 75% of tidal volume; and (3) spatial domain: quantifying global and regional tidal impedance variation (TIV), global inhomogeneity (GI), and center of ventilation (CoV). Extreme Gradient Boosting (XGBoost) classifiers were trained using 5 distinct data sources (Spirometry, Nasometry, Inspiratory-EIT, Expiratory-EIT, and fused ST-EIT). Model performance was rigorously evaluated via stratified 5-fold cross-validation, and Shapley additive explanations (SHAP) were employed to quantify global and local feature contributions. ResultsThe CP group exhibited a distinct respiratory phenotype compared to controls. In the temporal domain, CP patients showed significantly shorter inspiratory (1.60 s vs.1.85 s, P<0.001) and expiratory phase durations (2.45 s vs. 3.95 s, P<0.001), indicating a rapid, shallow breathing rhythm. In the airflow domain, while inspiratory flows were comparable, the CP group demonstrated significantly elevated mean and peak flows during the expiratory phase (P<0.001), reflecting compensatory respiratory effort. Spatially, CP patients presented significant ventilation redistribution, characterized by higher regional TIV in the right-anterior (ROI1) and left-posterior (ROI4) quadrants, but lower TIV in the left-anterior (ROI2) quadrant. In terms of diagnostic accuracy, the multi-modal ST-EIT model achieved the highest performance (AUC: 0.915±0.012, Accuracy: 0.843±0.019, F1-score: 0.872±0.017), substantially outperforming models based on spirometry (AUC: 0.721) or nasometry (AUC: 0.625) alone. Interpretability analysis revealed that spatial domain features were the most critical, contributing 53.4% to the model’s decision-making, followed by temporal (25.0%) and airflow (21.6%) features. ConclusionST-EIT successfully captures the temporal, airflow, and spatial deviations in CP speech respiration that are undetectable by conventional methods—specifically, rapid phase transitions, hyperdynamic expiratory airflow, and regional ventilation heterogeneity. This study validates ST-EIT as a robust, non-invasive, and radiation-free tool for characterizing speech-respiratory dysfunction, offering high clinical value for bedside screening, rehabilitation planning, and longitudinal monitoring of patients with cleft palate.

ObjectiveCleft palate (CP) is a common congenital deformity often associated with velopharyngeal insufficiency (VPI), which disrupts the physiological coupling between respiration and speech. Conventional clinical assessments, such as nasometry and spirometry, provide limited static data and fail to visualize the dynamic spatiotemporal distribution of lung ventilation during phonation. This study introduces spatiotemporal electrical impedance tomography (ST-EIT) to evaluate speech-respiratory functional features in CP patients compared to normal controls (NC). The aim is to characterize multi-domain respiratory patterns and to validate an interpretable machine learning framework for providing objective, quantitative evidence for clinical assessment. MethodsSeventy-five participants were enrolled in this study, comprising 37 patients with surgically repaired CP and 38 healthy volunteers matched for age, gender, and body mass index (BMI). All subjects performed standardized sustained phonation tasks while undergoing synchronous monitoring with a 16-electrode EIT system and a pneumotachograph. A comprehensive feature engineering pipeline was developed to extract physiological parameters across 3 complementary domains. (1) Temporal domain: including inspiratory/expiratory phase duration (tPhase), time constants (Tau), and inspiratory-to-expiratory time ratios (TI/TE); (2) airflow domain: comprising mean flow, peak flow, and instantaneous flow at 25%, 50%, and 75% of tidal volume; and (3) spatial domain: quantifying global and regional tidal impedance variation (TIV), global inhomogeneity (GI), and center of ventilation (CoV). Extreme Gradient Boosting (XGBoost) classifiers were trained using 5 distinct data sources (Spirometry, Nasometry, Inspiratory-EIT, Expiratory-EIT, and fused ST-EIT). Model performance was rigorously evaluated via stratified 5-fold cross-validation, and Shapley additive explanations (SHAP) were employed to quantify global and local feature contributions. ResultsThe CP group exhibited a distinct respiratory phenotype compared to controls. In the temporal domain, CP patients showed significantly shorter inspiratory (1.60 s vs.1.85 s, P<0.001) and expiratory phase durations (2.45 s vs. 3.95 s, P<0.001), indicating a rapid, shallow breathing rhythm. In the airflow domain, while inspiratory flows were comparable, the CP group demonstrated significantly elevated mean and peak flows during the expiratory phase (P<0.001), reflecting compensatory respiratory effort. Spatially, CP patients presented significant ventilation redistribution, characterized by higher regional TIV in the right-anterior (ROI1) and left-posterior (ROI4) quadrants, but lower TIV in the left-anterior (ROI2) quadrant. In terms of diagnostic accuracy, the multi-modal ST-EIT model achieved the highest performance (AUC: 0.915±0.012, Accuracy: 0.843±0.019, F1-score: 0.872±0.017), substantially outperforming models based on spirometry (AUC: 0.721) or nasometry (AUC: 0.625) alone. Interpretability analysis revealed that spatial domain features were the most critical, contributing 53.4% to the model’s decision-making, followed by temporal (25.0%) and airflow (21.6%) features. ConclusionST-EIT successfully captures the temporal, airflow, and spatial deviations in CP speech respiration that are undetectable by conventional methods—specifically, rapid phase transitions, hyperdynamic expiratory airflow, and regional ventilation heterogeneity. This study validates ST-EIT as a robust, non-invasive, and radiation-free tool for characterizing speech-respiratory dysfunction, offering high clinical value for bedside screening, rehabilitation planning, and longitudinal monitoring of patients with cleft palate.

ObjectiveProtein-protein interactions (PPIs) are fundamental to the execution of biological functions within living cells. However, traditional biochemical methods, such as co-immunoprecipitation (Co-IP), often fail to capture transient, weak, or membrane-associated interactions due to the stringent detergent requirements for cell lysis. Proximity labeling (PL) has emerged in recent years as a transformative technology for mapping the proteomes of specific subcellular compartments and identifying dynamic interactomes in situ. Golgi protein 73 (GP73, also known as GOLPH2), a resident type II Golgi transmembrane protein, is a well-recognized clinical biomarker for liver diseases, including hepatocellular carcinoma (HCC). Despite its clinical significance, the comprehensive physiological and pathological functions of GP73 remain partially understood. This study aims to establish an APEX2-mediated proximity labeling system specifically targeting GP73 to map its interactome in a living cellular environment, thereby providing new insights into its molecular roles and regulatory mechanisms. MethodsTo achieve spatial specificity, we first constructed a stable cell line expressing a fusion protein consisting of GP73 and the engineered soybean peroxidase APEX2. The localization of the GP73-APEX2 fusion protein was validated to ensure it correctly targeted the Golgi apparatus. The proximity labeling reaction was initiated by incubating the cells with biotin-phenol (BP) for 30 min, followed by a brief (1 min) treatment with1 mmol/L hydrogen peroxide (H₂O₂). This catalytic reaction converts BP into highly reactive, short-lived biotin-phenoxyl radicals that covalently attach to endogenous proteins within a small labeling radius of the GP73-APEX2 enzyme. Subsequently, the cells were quenched, and biotinylated proteins were enriched using high-affinity streptavidin-coated magnetic beads. The captured “neighbor” proteins were subjected to on-bead digestion and analyzed via liquid chromatography-tandem mass spectrometry (LC-MS/MS) for high-throughput identification. Rigorous bioinformatics analysis, including Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, and protein-protein interaction network mapping, was performed to interpret the biological significance of the identified candidates. ResultsOur results demonstrate the successful establishment of a robust and sensitive APEX2-based proximity labeling system for GP73. We identified a total of 95 high-confidence interacting proteins that were significantly enriched in the GP73 proximity proteome compared to control groups. Bioinformatics analysis revealed that these interactors were predominantly associated with biological processes such as vesicular transport, protein localization, and, most notably, molecular functions related to “ribosome binding” and “translation regulation”. This suggested an unexpected role for the Golgi-resident GP73 in the cellular translation machinery. To validate these findings, we performed targeted biochemical assays which confirmed a direct interaction between GP73 and the subunits of the eukaryotic translation initiation factor 3 (eIF3) complex, specifically EIF3G and EIF3I. Furthermore, functional validation using the surface sensing of translation (SUnSET) assay—a non-radioactive method to monitor protein synthesis—revealed that the overexpression of GP73 significantly promoted global protein translation levels in the cell, whereas its depletion or inhibition resulted in reduced translation efficiency. ConclusionThis study successfully utilized APEX2-mediated proximity labeling to provide the first systematic map of GP73 interactome in living cells. Our findings uncover a novel, unconventional function of GP73 as a regulator of cellular protein translation, likely mediated through its interaction with the eIF3 complex. This discovery significantly broadens our understanding of the biological roles of GP73 beyond its traditional function in the Golgi apparatus and suggests that it may act as a bridge between Golgi-related trafficking and the protein synthesis machinery. Furthermore, the technical framework established in this study provides a valuable template for investigating other complex organelle-associated protein networks and resolving transient macromolecular interactions in various physiological and pathological contexts.

ObjectivePhotoacoustic tomography (PAT) holds significant potential for high-resolution deep-tissue imaging. In preclinical research, custom-designed concave arc-shaped ultrasound transducer arrays are often used to maximize the detection aperture. However, manufacturing limitations and assembly tolerances frequently cause the actual physical positions of array elements to deviate from their theoretical design. Additionally, concave arrays are typically covered with an acoustic lens, which introduces a mismatch in the speed of sound between the coupling medium and the lens material. The combination of these geometric and acoustic-phase errors leads to severe image artifacts, reduced contrast, and degraded resolution. This study proposes a systematic two-step calibration strategy to address these issues and substantially improve image quality. MethodsFirst, a high-intensity isotropic photoacoustic point source was constructed using a multi-mode optical fiber coated with carbon nanotubes (CNTs) to acquire high signal-to-noise ratio calibration data. The Akaike information criterion (AIC) was employed to accurately determine the time of arrival (ToA) of photoacoustic signals. Subsequently, a geometric calibration algorithm based on nonlinear least-squares (NLS) estimation was developed. This algorithm iteratively solves for the true spatial coordinates of each array element by minimizing the residual between theoretical and measured acoustic path lengths. To further address sound-speed inhomogeneity caused by the acoustic lens, a phase compensation algorithm based on bilinear interpolation was proposed. This algorithm computes a pixel-specific phase delay map across the imaging region and performs point-by-point signal correction during delay-and-sum (DAS) reconstruction. The proposed methods were validated using a custom 96-channel concave arc-shaped array (center frequency: 12  MHz) through both phantom imaging and in vivo mouse tumor models. ResultsPhantom experiments showed that at an imaging depth of14 mm, the reconstruction position deviation of the point source in the uncalibrated system reached up to 1  mm. After applying the combined calibration, the lateral resolution (full width at half maximum, FWHM) at the focal point of the arc array reached 95  μm—representing a 85% reduction compared to the uncalibrated state and a 79% reduction compared to geometric calibration alone without phase compensation. In vivo experiments demonstrated that the calibrated system clearly resolved the microvascular network of subcutaneous tumors in mice. Photoacoustic signals were strictly confined within tumor boundaries delineated by ultrasound imaging (USI), eliminating the vascular spillover artifacts commonly observed in uncalibrated images. Furthermore, after intravenous injection of indocyanine green (ICG), the system successfully detected weak photoacoustic signals at a depth of 5 mm, performing significantly better than the uncalibrated system. ConclusionThe proposed calibration method, which integrates nonlinear least-squares estimation with phase compensation, significantly improves image fidelity and spatial resolution consistency across a wide field of view by correcting systemic geometric errors and acoustic phase aberrations. This approach demonstrates high robustness and provides a reliable technical foundation for the clinical translation of photoacoustic probes with non-standard geometries.