ObjectiveTo analyze the current status and existing challenges in talent cultivation for sports rehabilitation in China based on World Health Organization rehabilitation competency framework (RCF). MethodsData were collected from 104 higher education institutions nationwide that offer undergraduate programs in sports rehabilitation, including enrollment scale, regional distribution, degree conferral and training curriculum. Descriptive statistics and content analysis were used to examine the status and characteristics from three dimensions: institutional distribution, competency development and resource allocation， in terms of institutional type, regional distribution density, curriculum structure differences, competency-oriented training approaches and resource allocation patterns. ResultsTalent cultivation in sports rehabilitation in China currently faced a three-dimensional dilemma involving competency, resources, and public cognition. In terms of competency structure, a disconnect existed between medical fundamentals and exercise practice: physical education institutions provided insufficient medical training, whereas medical institutions lacked systematic instruction in exercise techniques. Regarding resource allocation, significant regional disparities were observed, with institutions heavily concentrated in Southwest (24.0%), East China (19.2%), and North China (15.4%), while Northwest China accounted for only 3.8%, forming a pronounced east-west gap. At the cognitive level, the public generally perceived sports rehabilitation as a service exclusive to athletes, and within the discipline, divergences persisted among the medical-oriented, sports-oriented, and integrated schools of thought. These challenges collectively hindered the quality of talent supply and regional balance. ConclusionTalent cultivation in sports rehabilitation in China is at a pivotal stage of transition from rapid expansion to quality enhancement. It is necessary to realign the knowledge structure and competency system of training programs with the five core competency domains of RCF.

Objective To investigate the expression profile, prognostic value, gene co-expression network, and immunomodulatory role of BRF1 in a pan-cancer context, and to explore its biological functions and molecular regulatory mechanisms in esophageal squamous cell carcinoma (ESCC). Methods The pan-cancer dataset from The Cancer Genome Atlas (TCGA) was utilized to analyze the differential expression of BRF1 in tumor versus normal tissues, its association with patient survival, pathway enrichment for co-expressed genes, and immune features (including immune checkpoints, cytokines, and immune cell infiltration). The expression profile of BRF1 in ESCC was validated using the Gene Expression Omnibus (GEO) database. In vitro, BRF1 was knocked down in ESCC cells using siRNA. Cell proliferation and migration were assessed by MTT and Transwell assays, respectively. The expression levels of proliferation- and migration-related proteins were detected by Western blotting. The correlation between BRF1 and ferroptosis was analyzed using TCGA data. Results BRF1 was significantly upregulated in over 20 types of cancer, and its high expression was associated with poor prognosis in patients with adrenocortical carcinoma and prostate adenocarcinoma. BRF1 was found to positively regulate the T-cell-mediated cell death pathway in esophageal adenocarcinoma and was associated with the circadian rhythm regulation pathway in pancreatic adenocarcinoma. The correlation of BRF1 with immune checkpoints, cytokine networks, and immune cell infiltration was found to be cancer type-specific. In vitro experiments demonstrated that knocking down BRF1 significantly inhibited the proliferation of ESCC cells, accompanied by the downregulation of the proliferation marker PCNA. Cell migration was also significantly impaired, with decreased expression of Vimentin and MMPs and increased expression of E-cadherin. Furthermore, the expression of BRF1 was positively correlated with that of ferroptosis-antagonizing genes, such as GPX4, HSPA5, and SLC7A11. Conclusion BRF1 plays complex roles in pan-cancer, participating in the regulation of tumorigenesis, progression, and immune infiltration. BRF1 promotes the proliferation and migration of ESCC cells, a mechanism potentially associated with the regulation of ferroptosis resistance. These findings suggest that BRF1 could be a potential therapeutic target for ESCC.

Esophageal cancer is a prevalent malignant tumor of the digestive tract in China, and radical surgery remains the cornerstone of its comprehensive treatment. However, multifactorial challenges such as postoperative gastrointestinal tract reconstruction, traumatic stress, and tumor-related metabolic disturbances render esophageal cancer patients highly susceptible to malnutrition. Perioperative nutritional support therapy plays a crucial role in enhancing surgical safety, improving clinical outcomes, and elevating patients' quality of life by regulating metabolic homeostasis, preserving organ function, and optimizing the immune microenvironment. This article reviews the mechanisms underlying malnutrition in esophageal cancer, methods for nutritional status assessment, and precision intervention pathways based on multi-omics evaluations. The aim is to strengthen clinicians' awareness of standardized perioperative nutritional management for esophageal cancer patients and promote its clinical implementation, thereby facilitating postoperative recovery and improving long-term quality of life.

OBJECTIVE: To investigate the current status of clinical genetics specialization development and the diagnostic and therapeutic capabilities for hereditary diseases across medical institutions in Shanghai, and to assess the necessity and feasibility of establishing training bases for clinical genetics specialists. METHODS: By employing a cross-sectional survey design, the Clinical Genetics Committee of Shanghai Medical Association has conducted questionnaire surveys from March to April 2025 across 54 healthcare institutions in Shanghai (including 33 tertiary hospitals and 21 secondary hospitals). The survey involved administrative departments and medical personnel from 15 clinical specialties. The survey has covered current genetic disease diagnosis and treatment practices, relevant and specialised disease types, genetic department establishment, testing capabilities, personnel teams, and training requirements. RESULTS: The results revealed that 78.0% of clinical departments surveyed had treated patients with hereditary disorders. Shanghai possesses diagnostic and therapeutic expertise for over 95% of hereditary diseases listed in its rare disease catalogue, reflecting both the practical clinical demand for such conditions and the city's overall diagnostic and therapeutic strengths in this field. Nevertheless, significant disparities exist in the development of genetics departments across different tiers of healthcare institutions. Resources for genetic testing capabilities (including molecular, cellular, and biochemical testing) are also unevenly distributed across different tiers of hospitals. The survey further revealed that only 26.0% of departments believe that their current physician structure fully meets the diagnostic and treatment demands. Over 90% of departments consider standard training for clinical genetic specialists necessary, with 74.0% expressing willingness to participate in establishing training bases. Based on above findings and thorough deliberation, the Clinical Genetics Committee of the Shanghai Medical Association proposes advancing specialist training and discipline development through establishing a standard training system. The committee has drafted a three-year training protocol featuring a "joint training"-centered model, recommending a pilot-first, dynamically optimized strategy for steadily advancing training base development. CONCLUSION Shanghai faces substantial demand for genetic disease diagnosis and treatment, yet exhibits shortcomings in clinical genetics specialization development, resource allocation, and talent pipeline cultivation. To establish a standard training system holds significant practical importance and is underpinned by a broad demand.
Humans ; China ; Surveys and Questionnaires ; Genetic Diseases, Inborn/genetics* ; Cross-Sectional Studies ; Genetics, Medical/education* ; Genetic Testing

Neurocritical care involves complex pathophysiological mechanisms, and its incidence is higher, injuries are more severe, and treatment is more challenging in high-altitude environments. This consensus, based on the latest domestic and international evidence-based medical data, establishes a standardized, goal-oriented framework for neurocritical care management applicable in high-altitude regions and nationwide. The consensus was developed following international standards for evidence quality assessment and underwent two rounds of Delphi expert consultation, resulting in 32 recommendation statements covering three parts: management systems, monitoring and assessment, and core strategies. Key updates include: advocating for the establishment of independent neurocritical care units and implementing precise tiered diagnosis and treatment based on the "Five Differences in Critical Care" concept; constructing a "trinity" multimodal brain monitoring system centered on cerebral blood flow, cerebral oxygenation, and brain function, emphasizing routine bedside transcranial Doppler ultrasound, cerebral oximetry, and continuous electroencephalography monitoring; shifting management strategies from mild hypothermia therapy to targeted temperature management, and defining the "446" target management pathway for the supercritical stage; emphasizing the assessment of static and dynamic cerebrovascular autoregulation functions through multimodal methods to achieve individualized optimal mean arterial pressure management; elevating cerebrospinal fluid management goals to the level of "glymphatic system" function maintenance; implementing a multidisciplinary collaborative, whole-process management model focusing on patients' long-term neurological functional outcomes; de-escalation criteria include multidimensional indicators such as recovery of brain structure, restoration of cerebrovascular autoregulation, improvement in cerebrospinal fluid dynamics, and reduction in biomarker levels; and integrating cutting-edge technologies like artificial intelligence into post-critical care management and rehabilitation planning. This consensus systematically integrates the entire process of neurocritical care management, reflecting the modern connotation of goal-oriented, dynamic, and multimodal integration in neurocritical care medicine. It aims to adapt to new trends such as deepening understanding of pathophysiological mechanisms, the integration of medicine and engineering, and the empowerment of artificial intelligence, thereby further advancing the discipline of critical care medicine.

ObjectiveSpinal cord injury (SCI) directly impairs the regulatory function of the autonomic nervous system, induces intestinal dysfunction, and significantly reduces patients’ quality of life. Preclinical studies have shown that electroacupuncture (EA) therapy can regulate the brain-gut axis and is used to treat central nervous system diseases such as major depressive disorder, Alzheimer’s disease and Parkinson’s disease. Recent research has established that fecal microbiota transplantation (FMT) from EA-treated SCI rats restored intestinal motility and colonic morphology. However, it remains unclear whether the regulation of gut microbiota by EA therapy directly contributes to neural repair after SCI. This study aims to explore whether gut microbiota mediates the neuroprotective effect of EA in the treatment of SCI and its possible mechanism. MethodsThe study employed RNA transcriptome analysis of spinal cord tissue to characterize gene expression profiles and to identify key signaling pathways following EA treatment for SCI. Hematoxylin-Eosin (HE) staining and Nissl staining were used to observe the morphological changes in spinal cord tissue. Western blot (WB) and enzyme-linked immunosorbent assay (ELISA) were applied to detect the effects of EA on the expression of proteins related to nucleotide-binding domain leucine-rich repeat and pyrin domain-containing receptor 3 (NLRP3) -dependent pyroptosis. Using 16S rDNA sequencing, the study observed alterations in gut microbiota diversity and community composition in SCI rats. Prior to establishing SCI models, rats were pretreated with an antibiotic cocktail to induce gut dysbiosis, and the effects on intestinal function and spinal cord neural repair were evaluated. FMT was performed to investigate the regulatory effects of post-EA FMT on motor function, general status, liver and spleen indices, and NLRP3-mediated pyroptosis in SCI rats. ResultsEA improved motor function and reduced regulated neuronal cell death in SCI rats. Transcriptomic analysis demonstrated the activation of immune- and inflammation-related pathways post-SCI, including NOD-like receptors, nuclear factor-kappa B(NF-κB), and Toll-like receptor (TLR) pathways. EA primarily influenced intestinal inflammation and autoimmune functions. 16S rDNA sequencing illustrated that EA did not alter the diversity of gut microbiota. However, EA altered the gut microbiota composition in SCI rats, increasing Lactobacillus and Akkermansia genera while rebalancing the Firmicutes/Bacteroidetes ratio. Furthermore, depletion of gut microbiota by antibiotics disrupted the intestinal barrier, reduced the expression of intestinal barrier proteins Zonula Occludens-1 (ZO-1) and Occludin, elevated serum lipopolysaccharide-binding protein (LBP) levels, exacerbated spinal cord tissue damage, and hindered motor function recovery in SCI rats. FMT from donors treated with EA reduced LBP levels in the intestine, blood, and spinal cord of rats, inhibited the TLR4 myeloid differentiation primary response protein 88 (MyD88)-NF‑κB pathway and NLRP3-dependent pyroptosis, and improved motor function. On the other hand, FMT treatment resulted in decreased body weight and food intake, whereas FMT using EA-treated donors effectively alleviated these alterations. ConclusionEA effectively alleviated neuroinflammatory responses in rats with SCI, primarily through regulating the gut microbiota and suppressing the NLRP3-dependent pyroptosis signaling pathway.

ObjectiveChronic atrophic gastritis (CAG) is usually diagnosed by gastroscopy and histopathological biopsy. These procedures remain the reference standard, but their invasive nature and resource requirements may limit their use in large-scale population screening and repeated follow-up. A convenient and reproducible method for noninvasive auxiliary screening may help identify individuals who require further endoscopic assessment. Fingertip photoplethysmography (PPG) provides a noninvasive recording of peripheral pulse waves and allows harmonic features to be extracted from the signal. In this study, the so-called meridian-related variables were defined as PPG-derived harmonic parameters labelled according to meridian nomenclature, rather than as direct measurements of meridian physiology. This study aimed to compare these harmonic parameters between patients with CAG and non-CAG controls, identify parameters that remained different after age adjustment, and develop a multivariable model for noninvasive auxiliary screening and pre-endoscopic risk stratification of CAG. MethodsA total of 343 participants were included, comprising 171 patients with CAG and 172 non-CAG controls. CAG diagnosis was established using gastroscopy and histopathology as the reference standard. Fingertip PPG signals were collected using a PPG-based pulse acquisition device. Eight PPG-derived harmonic parameters labelled according to meridian nomenclature were extracted for analysis. Between-group differences were first assessed using nonparametric tests. Age-adjusted analyses were then performed to reduce potential confounding by age. The false discovery rate (FDR) method was applied for multiple-comparison correction. A multivariable logistic regression model integrating age and multiple harmonic parameters was constructed. Model performance was evaluated using receiver operating characteristic (ROC) analysis and the area under the curve (AUC). Internal validation performance was assessed using stratified five-fold cross-validation and bootstrap optimism correction. Threshold performance was examined using both a high-specificity strategy and a Youden index-based cutoff. Decision curve analysis was used to evaluate the model’s net clinical benefit across a range of threshold probabilities. ResultsAll eight harmonic parameters were non-normally distributed. In the univariate analysis, the stomach-labelled harmonic parameter (ST), bladder-labelled harmonic parameter (BL), and liver-labelled harmonic parameter (LR) differed between the CAG and non-CAG groups. After age adjustment and FDR correction, only ST and BL remained statistically significant. Compared with non-CAG controls, patients with CAG showed higher ST values and lower BL values. This finding indicates an associated differential harmonic pattern that was not fully explained by age distribution. However, the discriminative ability of a single harmonic parameter was limited. The best-performing single indicator was ST, with an AUC of 0.652 (95% CI: 0.595-0.707). The multivariable model integrating age and multiple harmonic parameters achieved an AUC of 0.791 (95% CI: 0.743-0.835), representing an improvement of 0.139 over ST alone. In internal validation, stratified five-fold cross-validation yielded a mean AUC of 0.753 (95% CI: 0.715-0.781), and the bootstrap optimism-corrected AUC was 0.748. These results suggest that the model retained moderate discriminative performance after internal validation.At a specificity of at least 95%, the model achieved a sensitivity of only 40.4% (95% CI: 25.7%-49.7%). This high-specificity cutoff may be suboptimal as the preferred threshold for an initial screening setting because of the potential risk of missed CAG cases. The Youden index-based optimal cutoff was 0.419, corresponding to a sensitivity of 80.7% and a specificity of 62.8%. This threshold may better match the practical aim of noninvasive auxiliary screening, where sensitivity is usually prioritized to reduce missed cases. Decision curve analysis showed that, within a threshold probability range of 10%-55%, the model provided higher net clinical benefit than the reference strategies of recommending gastroscopy for all participants or for none. ConclusionPatients with CAG showed associated harmonic differences in fingertip PPG-derived features, mainly characterized by higher ST and lower BL values after age adjustment and FDR correction. Compared with a single harmonic parameter, the multivariable model showed better overall discrimination and retained moderate internal validation performance. These findings suggest that PPG-derived harmonic parameters labelled according to meridian nomenclature may provide auxiliary information for noninvasive auxiliary screening and front-line triage before gastroscopic confirmation in CAG. The present results support further validation rather than immediate clinical implementation. External validation in independent, multicenter, and preferably prospective screening cohorts is needed to assess the model’s generalizability, screening performance, and potential clinical utility.

ObjectiveSpinal cord injury (SCI) directly impairs the regulatory function of the autonomic nervous system, induces intestinal dysfunction, and significantly reduces patients’ quality of life. Preclinical studies have shown that electroacupuncture (EA) therapy can regulate the brain-gut axis and is used to treat central nervous system diseases such as major depressive disorder, Alzheimer’s disease and Parkinson’s disease. Recent research has established that fecal microbiota transplantation (FMT) from EA-treated SCI rats restored intestinal motility and colonic morphology. However, it remains unclear whether the regulation of gut microbiota by EA therapy directly contributes to neural repair after SCI. This study aims to explore whether gut microbiota mediates the neuroprotective effect of EA in the treatment of SCI and its possible mechanism. MethodsThe study employed RNA transcriptome analysis of spinal cord tissue to characterize gene expression profiles and to identify key signaling pathways following EA treatment for SCI. Hematoxylin-Eosin (HE) staining and Nissl staining were used to observe the morphological changes in spinal cord tissue. Western blot (WB) and enzyme-linked immunosorbent assay (ELISA) were applied to detect the effects of EA on the expression of proteins related to nucleotide-binding domain leucine-rich repeat and pyrin domain-containing receptor 3 (NLRP3) -dependent pyroptosis. Using 16S rDNA sequencing, the study observed alterations in gut microbiota diversity and community composition in SCI rats. Prior to establishing SCI models, rats were pretreated with an antibiotic cocktail to induce gut dysbiosis, and the effects on intestinal function and spinal cord neural repair were evaluated. FMT was performed to investigate the regulatory effects of post-EA FMT on motor function, general status, liver and spleen indices, and NLRP3-mediated pyroptosis in SCI rats. ResultsEA improved motor function and reduced regulated neuronal cell death in SCI rats. Transcriptomic analysis demonstrated the activation of immune- and inflammation-related pathways post-SCI, including NOD-like receptors, nuclear factor-kappa B(NF-κB), and Toll-like receptor (TLR) pathways. EA primarily influenced intestinal inflammation and autoimmune functions. 16S rDNA sequencing illustrated that EA did not alter the diversity of gut microbiota. However, EA altered the gut microbiota composition in SCI rats, increasing Lactobacillus and Akkermansia genera while rebalancing the Firmicutes/Bacteroidetes ratio. Furthermore, depletion of gut microbiota by antibiotics disrupted the intestinal barrier, reduced the expression of intestinal barrier proteins Zonula Occludens-1 (ZO-1) and Occludin, elevated serum lipopolysaccharide-binding protein (LBP) levels, exacerbated spinal cord tissue damage, and hindered motor function recovery in SCI rats. FMT from donors treated with EA reduced LBP levels in the intestine, blood, and spinal cord of rats, inhibited the TLR4 myeloid differentiation primary response protein 88 (MyD88)-NF‑κB pathway and NLRP3-dependent pyroptosis, and improved motor function. On the other hand, FMT treatment resulted in decreased body weight and food intake, whereas FMT using EA-treated donors effectively alleviated these alterations. ConclusionEA effectively alleviated neuroinflammatory responses in rats with SCI, primarily through regulating the gut microbiota and suppressing the NLRP3-dependent pyroptosis signaling pathway.

ObjectiveChronic atrophic gastritis (CAG) is usually diagnosed by gastroscopy and histopathological biopsy. These procedures remain the reference standard, but their invasive nature and resource requirements may limit their use in large-scale population screening and repeated follow-up. A convenient and reproducible method for noninvasive auxiliary screening may help identify individuals who require further endoscopic assessment. Fingertip photoplethysmography (PPG) provides a noninvasive recording of peripheral pulse waves and allows harmonic features to be extracted from the signal. In this study, the so-called meridian-related variables were defined as PPG-derived harmonic parameters labelled according to meridian nomenclature, rather than as direct measurements of meridian physiology. This study aimed to compare these harmonic parameters between patients with CAG and non-CAG controls, identify parameters that remained different after age adjustment, and develop a multivariable model for noninvasive auxiliary screening and pre-endoscopic risk stratification of CAG. MethodsA total of 343 participants were included, comprising 171 patients with CAG and 172 non-CAG controls. CAG diagnosis was established using gastroscopy and histopathology as the reference standard. Fingertip PPG signals were collected using a PPG-based pulse acquisition device. Eight PPG-derived harmonic parameters labelled according to meridian nomenclature were extracted for analysis. Between-group differences were first assessed using nonparametric tests. Age-adjusted analyses were then performed to reduce potential confounding by age. The false discovery rate (FDR) method was applied for multiple-comparison correction. A multivariable logistic regression model integrating age and multiple harmonic parameters was constructed. Model performance was evaluated using receiver operating characteristic (ROC) analysis and the area under the curve (AUC). Internal validation performance was assessed using stratified five-fold cross-validation and bootstrap optimism correction. Threshold performance was examined using both a high-specificity strategy and a Youden index-based cutoff. Decision curve analysis was used to evaluate the model’s net clinical benefit across a range of threshold probabilities. ResultsAll eight harmonic parameters were non-normally distributed. In the univariate analysis, the stomach-labelled harmonic parameter (ST), bladder-labelled harmonic parameter (BL), and liver-labelled harmonic parameter (LR) differed between the CAG and non-CAG groups. After age adjustment and FDR correction, only ST and BL remained statistically significant. Compared with non-CAG controls, patients with CAG showed higher ST values and lower BL values. This finding indicates an associated differential harmonic pattern that was not fully explained by age distribution. However, the discriminative ability of a single harmonic parameter was limited. The best-performing single indicator was ST, with an AUC of 0.652 (95% CI: 0.595-0.707). The multivariable model integrating age and multiple harmonic parameters achieved an AUC of 0.791 (95% CI: 0.743-0.835), representing an improvement of 0.139 over ST alone. In internal validation, stratified five-fold cross-validation yielded a mean AUC of 0.753 (95% CI: 0.715-0.781), and the bootstrap optimism-corrected AUC was 0.748. These results suggest that the model retained moderate discriminative performance after internal validation.At a specificity of at least 95%, the model achieved a sensitivity of only 40.4% (95% CI: 25.7%-49.7%). This high-specificity cutoff may be suboptimal as the preferred threshold for an initial screening setting because of the potential risk of missed CAG cases. The Youden index-based optimal cutoff was 0.419, corresponding to a sensitivity of 80.7% and a specificity of 62.8%. This threshold may better match the practical aim of noninvasive auxiliary screening, where sensitivity is usually prioritized to reduce missed cases. Decision curve analysis showed that, within a threshold probability range of 10%-55%, the model provided higher net clinical benefit than the reference strategies of recommending gastroscopy for all participants or for none. ConclusionPatients with CAG showed associated harmonic differences in fingertip PPG-derived features, mainly characterized by higher ST and lower BL values after age adjustment and FDR correction. Compared with a single harmonic parameter, the multivariable model showed better overall discrimination and retained moderate internal validation performance. These findings suggest that PPG-derived harmonic parameters labelled according to meridian nomenclature may provide auxiliary information for noninvasive auxiliary screening and front-line triage before gastroscopic confirmation in CAG. The present results support further validation rather than immediate clinical implementation. External validation in independent, multicenter, and preferably prospective screening cohorts is needed to assess the model’s generalizability, screening performance, and potential clinical utility.

Background The contamination of public baths with Legionella pneumophila contamination has become a growing public health concern in recent years. However, research on its association with bacterial community characteristics in water samples remains limited. The integration of 16S rRNA sequencing and random forest modeling provides a new approach to elucidate the bacterial community characteristics of public bath water and their association with Legionella pneumophila contamination. Objective To investigate the bacterial community structure and diversity of public bath water in Shanghai, explore the association between Legionella pneumophila contamination and bacterial community characteristics, and identify key bacterial genera associated with contamination, thereby providing a scientific basis for formulating hygiene management regulations for public bath water. Methods From February to March 2023, water samples were collected from ten public baths in Shanghai which were selected based on business scale, regional distribution, and functional differences. Water quality parameters were evaluated, and the samples were categorized into Legionella-positive and Legionella-negative groups based on the detection results of Legionella pneumophila. The bacterial community structure, α-diversity, and β-diversity were analyzed using 16S rRNA sequencing. Redundancy analysis (RDA) was employed to examine the relationship between physicochemical factors and bacterial community diversity. A random forest model was employed to identify key bacterial genera distinguishing the two groups, with the importance of genera being evaluated based on the mean decrease accuracy (MDA). Results The oxygen consumption in the Legionella-positive group was significantly lower than that in the Legionella-negative group (mean values: 1.85 mg·L⁻¹ vs. 6.81 mg·L⁻¹, P< 0.05), while no significant differences were observed in other physicochemical indicators. The sequencing results revealed a total of 27 bacterial phyla and 454 bacterial genera, with Proteobacteria (63.00%) being the dominant phylum. The dominant genera included Pelomonas (8.50%), Acidovorax (8.13%), Mycobacterium (7.93%), and Acinetobacter (6.59%). The α-diversity analysis indicated that bacterial community richness (Chao1 and ACE indices) was significantly higher in the Legionella-positive group than in the Legionella-negative group (P<0.01). The β-diversity analysis showed no significant difference in the bacterial community structure between the two groups (P>0.05). The RDA analysis demonstrated that the bacterial community diversity was positively correlated with pH and negatively correlated with oxygen consumption and free residual chlorine. The RDA1 and RDA2 explained 23.92% and 21.30% of the bacterial community diversity, respectively. The random forest model identified 20 key genera significantly influencing the microbial community distribution between the two groups, including unclassified_Bradyrhizobiaceae (MDA=2.42), Meiothermus (MDA=2.37), and Flavihumibacter (MDA=2.26). Conclusion The diversity of bacterial communities in public bath water is influenced by pH, oxygen consumption, and free residual chlorine. Samples contaminated with Legionella pneumophila exhibit greater microbial richness and contain characteristic key bacterial genera that contribute to community differences. Machine learning random forest technology helps identify these distinctive key bacterial genera. The findings provide a basis for carrying out risk early warning strategies in such settings.