Environmental pollution is a significant public health challenge worldwide, and investigating the causal relationship between environmental exposure and population health outcomes is a key objective of environmental epidemiology research. In recent years, the complexity of environmental exposures has increasingly come to the forefront, making it challenging for observational studies that dominate environmental epidemiology to accurately estimate causal effects. Causal inference methods are particularly advantageous in controlling for confounding factors, thus holding great potential in environmental epidemiology research. Researchers can use appropriate causal inference methods to simulate the process of randomization, providing strong support for revealing the causal relationship between environmental exposure and health outcomes. However, there is a lack of reviews on the application of causal inference methods in environmental epidemiology studies in China. Therefore, this study introduced the basic principles of common causal inference statistical methods in environmental epidemiology, summarized the applicable conditions, advantages and disadvantages of various methods, and provided R software implementation codes for these methods, aiming to offer guidance for optimizing research design and practicing causal inference statistical methods.

Background Ozone (O₃) is a major air pollutant. The existing monitoring system has uneven distribution of sites, insufficient coverage in underdeveloped areas, and low temporal resolution, making it difficult to obtain hourly data. This limits the dynamic identification of pollution and the formulation of prevention and control strategies. Objective To construct an hourly O₃ concentration estimation model based on ensemble machine learning, aiming to improve the accuracy of pollution exposure assessment and explore O₃ health impacts. Methods This study integrated land use regression modeling with modern machine learning techniques, employing random forest and XGBoost algorithms to construct base models, and stacking integration using non-negative least squares. The ensemble model was trained and validated across China using high-resolution, multi-source geographic data (e.g., meteorologicaldata, population density, land cover types, and aerosol optical thickness). It was tested in Taiyuan City, combined with a distributed lag non-linear model to analyze the association between O₃ and emergency admissions. Results The constructed ensemble model performed well in predicting O₃ concentration, with a higher coefficient of determination (R²) and a lower root-mean-square deviation (RMSE) compared to the single models. The R² improved from 0.90 to 0.92, and the RMSE decreased from 11.41 to 10.62, enhancing both prediction accuracy and generalization ability. In the application to Taiyuan City, the model successfully imputed the hourly-level data for the entire year. The distributed lag non-linear model analysis revealed that the relative risk (RR) values for the 6th to 8th days following O₃ exposure were 1.14 (95%CI: 1.01, 1.29), 1.16 (95%CI: 1.02, 1.31), and 1.14 (95%CI: 1.01, 1.29), respectively, which were significantly higher than 1, indicating a significant lagged association (lagged 6-8 d) between O₃ and the number of emergency room visits. Conclusion A high-precision, hourly-level O₃ concentration estimation model is successfully constructed by combining the land use regression model with an ensemble machine learning approach to provide a scientific basis for environmental policy formulation and public health intervention. The application of the model verifies its generalization ability and practical application value, which can provide a new technical framework for subsequent environmental health research.

Objective: To understand the correlation between autistic traits in children with autism spectrum disorder(ASD) and their parental behaviors, so as to provide evidence for improving social behavior in children with ASD. Methods: A total of 62 children aged 4-7 years diagnosed with ASD at Guangdong Provincial Work Injury Rehabilitation Hospital and Guangzhou Yuexiu District Children s Hospital from October 2021 to May 2025 were selected as the ASD group. A random number table method was used to select 62 healthy children of the same age from 3 ordinary kindergartens in Guangzhou as the control (typically developing,TD) group. The Parental Behavior Inventory and Social Responsiveness Scale were used to evaluate parental behaviors and autistic traits in both groups. The t-test and multiple linear regression were used to analyze the relationship between autistic traits in children with ASD and parental behaviors. Results: Compared with the TD group, the ASD group had significantly higher scores in social awareness, social cognition, social communication, social motivation, autistic behavior patterns, and the total score of social responsiveness, with statistically significant differences ( t =35.83, 46.17, 64.36, 41.45, 49.46,101.15, all P <0.01). The score of parental support/participation in the ASD group (24.61±4.30) was lower than that in the TD group (31.27±4.58), while the score of parental hostility/coercion in the ASD group (25.51±2.72) was higher than that in the TD group (12.75±2.72), with statistically significant differences ( t = -8.34, 26.14, both P <0.01). The total score of the Social Responsiveness Scale was negatively correlated with the dimension score of parental support/participation behaviors, and positively correlated with the dimension score of parental hostility/coercion behaviors ( r =-0.60, 0.91,both P <0.05). After adjusting for covariates such as children s age, gender, and primary caregiver, ASD children with exclusive breastfeeding ( β =-8.79) and parental support/participation behaviors ( β = -0.79 ) had a lower risk of autistic traits, while children with parental hostility/coercion behaviors ( β =4.62) had a higher risk of autistic traits (all P <0.05). Conclusion Autistic traits in children with ASD may be related to their parental behaviors and early feeding mode.

Objective:To compare the efficacy of anrikefon and tegileridine for analgesia in patients with moderate-to-severe pain after abdominal surgery with general anesthesia.Methods:In this multicenter, randomized, double-blind, active-controlled clinical trial, 101 patients with moderate to severe pain (numeric pain rating scale [NRS] score ≥4 within 4 h after operation) after abdominal surgery with general anesthesia between February 24 and April 1, 2025, aged 18-70 yr, with a body mass index of 18-40 kg/m 2, were assigned to anrikefon group ( n=50) and tegileridine group ( n=51) in a 1∶1 ratio using stratified blocked randomization. Double-dummy design was employed to maintain blinding. Each group received an initial intravenous injection of anrikefon 1 μg/kg or tegileridine 1 mg, followed by connection to a patient-controlled intravenous analgesia (PCIA) pump (the PCIA solution contained normal saline in anrikefon group; the PCIA solution contained tegileridine 5 mg in tegileridine pump) within 10 min. If the patient′s NRS score ≥4 at 8 and 16 h after the initial injection, anrikefon 1 μg/kg was intravenously injected in anrikefon group, and tegileridine group received the equal volume of normal saline. The primary efficacy endpoint was the sum of pain intensity difference (SPID) over the first 24 h after the initial dose (SPID 0-24h). The secondary efficacy endpoints included the incidence and severity of vomiting and nausea, incidence of postoperative nausea and vomiting(PONV), the proportion of patients who received antiemetic treatment, and total consumption of antiemetics within 0-24 h after the initial dose, NRS score at rest ≤ 1 at 24 h after the initial dose, and NRS score at rest ≤ 3 over the first 24 h after the initial dose. Safety indicators included adverse events, vital signs, physical examination findings, 12-lead ECG and laboratory test indicators, and adverse events of special interest. Results:Compared with tegileridine group, no significant change was found in the SPID 0-24h ( P>0.05), and the incidence of vomiting, PONV, proportion of patients requiring antiemetic medication, and total consumption of antiemetics were significantly decreased within the first 24 h after the initial dose in tegileridine group ( P<0.05). One treatment-emergent adverse event of Common Terminology Criteria for Adverse Events grade 3 or higher occurred in tegileridine group, while no treatment-emergent adverse events of Common Terminology Criteria for Adverse Events grade 3 or higher were found in anrikefon group. Among the adverse events of special interest, one case of respiratory depression and one case of cough occurred in tegileridine group, while one case of cough occurred in anrikefon group, with no respiratory depression. Conclusions:Anrikefon and tegileridine provide comparable analgesic efficacy for moderate-to-severe pain after abdominal surgery with general anesthesia. However, anrikefon exhibits an advantage in reducing the risk of PONV, with a superior safety profile.

OBJECTIVE: To assess the value of combined detection strategies using multiple technologies for the genetic testing and prenatal diagnosis for pedigrees affected with Duchenne muscular dystrophy (DMD) for optimizing genetic counseling and reproductive guidance. METHODS: This study has involved 142 subjects from 65 suspected DMD families who had visited the First Affiliated Hospital of Chongqing Medical University from January 2018 to December 2023. A combination of multiple ligation-dependent probe amplification (MLPA), quantitative fluorescence PCR, and next-generation sequencing (NGS) was used. After confirming the genetic diagnosis of the probands, prenatal diagnosis was provided for carrier mothers. This study was approved by the Medical Ethics Committee of the hospital (Ethics No.: 2021-264). RESULTS: Among the 142 subjects tested, 73 cases of large deletions/duplications and 15 cases of small variants of the DMD gene were detected. The hotspot regions for the variants were exons 45 to 55. A total of 41 variant types were identified, of which 3 were previously unreported. In 19 families with suspected patients, 7 exonic deletions, 2 exonic duplications, and 3 small variants were identified. Prenatal diagnosis was performed on 48 fetuses from 46 families, revealing 16 affected male fetuses (including 12 with deletion variants, 2 with duplication variants, and 2 with small variants). Seven carrier females were identified among the 16 female fetuses (including 6 with deletions and 1 with duplication). Among the couples with an affected fetus, 16 had opted to terminate the pregnancy, while the parents of 32 fetuses had chosen to continue with the pregnancy. In families undergoing prenatal diagnosis, 53 (79.1%) pregnant women and their family members were found to carry mutations of the DMD gene. CONCLUSION The combined detection strategy of MLPA, qPCR, and NGS can encompass large deletions/duplications and small variants of the DMD gene, providing timely and accurate prenatal diagnosis for families affected by DMD. In conjunction with genetic counseling, this can effectively reduce the risk of producing affected offspring, which is crucial for the prevention and control of this disease.
Humans ; Muscular Dystrophy, Duchenne/diagnosis* ; Prenatal Diagnosis/methods* ; Female ; Male ; Pregnancy ; Pedigree ; Genetic Testing/methods* ; Dystrophin/genetics* ; Adult ; Genetic Counseling ; High-Throughput Nucleotide Sequencing/methods* ; Exons

ObjectiveTo investigate the ameliorative effect of paeonol on acute alcohol-induced hepatic inflammation in mice via the regulation of the short-chain fatty acids （SCFAs）-specific receptor GPR43/mitogen-activated protein kinase （MAPK） signaling pathway. MethodsC57BL/6 mice were randomly divided into five groups： blank control group， model group， low-dose paeonol group （120 mg·kg^-1）， high-dose paeonol group （480 mg·kg^-1）， and silybin group （36.8 mg·kg^-1）. A mouse model of alcohol-induced liver disease （ALD） was established by ad libitum administration of a Lieber-DeCarli alcohol liquid diet. Serum lipid levels， liver function， inflammatory cytokines， and oxidative stress markers were measured. Liver hematoxylin-eosin （HE） staining and Oil Red O staining were performed to validate successful modeling. Western blot analysis was used to assess the expression levels of zonula occludens-1 （ZO-1）， Claudin-1， and proteins related to the GPR43/MAPK signaling pathway in the colonic tissue. Immunohistochemistry was employed to detect the protein expression of GPR43， ZO-1， and Claudin-1 in the colon. Then 16S rDNA sequencing was performed to analyze differences in intestinal flora between the model group and the high-dose paeonol group. Additionally， fecal microbiota transplantation （FMT） experiments were conducted to validate the regulatory effect of paeonol on ALD via modulation of intestinal flora. ResultsCompared with the blank control group， the model group showed significantly elevated serum lipid levels， oxidative stress， and inflammatory cytokine expression （P<0.01）. Liver histology revealed increased inflammatory infiltration and lipid droplet accumulation. Colonic mucosal injury and impaired intestinal barrier function were observed. Levels of MAPK pathway-related proteins in the colonic tissue were upregulated （P<0.01）， while GPR43， ZO-1， and Claudin-1 protein expression levels were significantly decreased （P<0.01）. The composition and abundance of the intestinal flora were markedly altered， with a reduced Bacteroidetes-to-Firmicutes ratio and decreased relative abundances of Eubacterium， Parabacteroides， Erysipelothrix， and Adlercreutzia， alongside increased abundances of Clostridium butyricum， Enterococcus， and Helicobacter pylori in the model group. Compared with the model group， paeonol significantly reduced serum lipid levels， oxidative stress responses， and the expression of inflammatory cytokines in ALD mice （P<0.05， P<0.01）. It also attenuated hepatic lipid accumulation， restored intestinal barrier function， and repaired the structural integrity of liver and colonic tissues. The protein expression levels of ZO-1， Claudin-1， and GPR43 in the colonic tissue were significantly increased （P<0.05， P<0.01）， while those of MAPK pathway-related proteins were significantly decreased （P<0.05， P<0.01）. The intestinal flora dysbiosis was effectively alleviated， rendering its composition closer to that of normal mice. The efficacy of paeonol in modulating ALD was further confirmed by FMT experiments， supporting its mechanistic involvement in the SCFAs-GPR43/MAPK signaling pathway. ConclusionPaeonol exerts a protective effect against ALD in mice， which may be mediated through regulation of the SCFAs-GPR43/MAPK signaling pathway， thereby achieving anti-inflammatory effects and improving intestinal barrier function.

ObjectiveTo investigate the ameliorative effect of paeonol on acute alcohol-induced hepatic inflammation in mice via the regulation of the short-chain fatty acids （SCFAs）-specific receptor GPR43/mitogen-activated protein kinase （MAPK） signaling pathway. MethodsC57BL/6 mice were randomly divided into five groups： blank control group， model group， low-dose paeonol group （120 mg·kg^-1）， high-dose paeonol group （480 mg·kg^-1）， and silybin group （36.8 mg·kg^-1）. A mouse model of alcohol-induced liver disease （ALD） was established by ad libitum administration of a Lieber-DeCarli alcohol liquid diet. Serum lipid levels， liver function， inflammatory cytokines， and oxidative stress markers were measured. Liver hematoxylin-eosin （HE） staining and Oil Red O staining were performed to validate successful modeling. Western blot analysis was used to assess the expression levels of zonula occludens-1 （ZO-1）， Claudin-1， and proteins related to the GPR43/MAPK signaling pathway in the colonic tissue. Immunohistochemistry was employed to detect the protein expression of GPR43， ZO-1， and Claudin-1 in the colon. Then 16S rDNA sequencing was performed to analyze differences in intestinal flora between the model group and the high-dose paeonol group. Additionally， fecal microbiota transplantation （FMT） experiments were conducted to validate the regulatory effect of paeonol on ALD via modulation of intestinal flora. ResultsCompared with the blank control group， the model group showed significantly elevated serum lipid levels， oxidative stress， and inflammatory cytokine expression （P<0.01）. Liver histology revealed increased inflammatory infiltration and lipid droplet accumulation. Colonic mucosal injury and impaired intestinal barrier function were observed. Levels of MAPK pathway-related proteins in the colonic tissue were upregulated （P<0.01）， while GPR43， ZO-1， and Claudin-1 protein expression levels were significantly decreased （P<0.01）. The composition and abundance of the intestinal flora were markedly altered， with a reduced Bacteroidetes-to-Firmicutes ratio and decreased relative abundances of Eubacterium， Parabacteroides， Erysipelothrix， and Adlercreutzia， alongside increased abundances of Clostridium butyricum， Enterococcus， and Helicobacter pylori in the model group. Compared with the model group， paeonol significantly reduced serum lipid levels， oxidative stress responses， and the expression of inflammatory cytokines in ALD mice （P<0.05， P<0.01）. It also attenuated hepatic lipid accumulation， restored intestinal barrier function， and repaired the structural integrity of liver and colonic tissues. The protein expression levels of ZO-1， Claudin-1， and GPR43 in the colonic tissue were significantly increased （P<0.05， P<0.01）， while those of MAPK pathway-related proteins were significantly decreased （P<0.05， P<0.01）. The intestinal flora dysbiosis was effectively alleviated， rendering its composition closer to that of normal mice. The efficacy of paeonol in modulating ALD was further confirmed by FMT experiments， supporting its mechanistic involvement in the SCFAs-GPR43/MAPK signaling pathway. ConclusionPaeonol exerts a protective effect against ALD in mice， which may be mediated through regulation of the SCFAs-GPR43/MAPK signaling pathway， thereby achieving anti-inflammatory effects and improving intestinal barrier function.

Objective:To compare the efficacy of anrikefon and tegileridine for analgesia in patients with moderate-to-severe pain after abdominal surgery with general anesthesia.Methods:In this multicenter, randomized, double-blind, active-controlled clinical trial, 101 patients with moderate to severe pain (numeric pain rating scale [NRS] score ≥4 within 4 h after operation) after abdominal surgery with general anesthesia between February 24 and April 1, 2025, aged 18-70 yr, with a body mass index of 18-40 kg/m 2, were assigned to anrikefon group ( n=50) and tegileridine group ( n=51) in a 1∶1 ratio using stratified blocked randomization. Double-dummy design was employed to maintain blinding. Each group received an initial intravenous injection of anrikefon 1 μg/kg or tegileridine 1 mg, followed by connection to a patient-controlled intravenous analgesia (PCIA) pump (the PCIA solution contained normal saline in anrikefon group; the PCIA solution contained tegileridine 5 mg in tegileridine pump) within 10 min. If the patient′s NRS score ≥4 at 8 and 16 h after the initial injection, anrikefon 1 μg/kg was intravenously injected in anrikefon group, and tegileridine group received the equal volume of normal saline. The primary efficacy endpoint was the sum of pain intensity difference (SPID) over the first 24 h after the initial dose (SPID 0-24h). The secondary efficacy endpoints included the incidence and severity of vomiting and nausea, incidence of postoperative nausea and vomiting(PONV), the proportion of patients who received antiemetic treatment, and total consumption of antiemetics within 0-24 h after the initial dose, NRS score at rest ≤ 1 at 24 h after the initial dose, and NRS score at rest ≤ 3 over the first 24 h after the initial dose. Safety indicators included adverse events, vital signs, physical examination findings, 12-lead ECG and laboratory test indicators, and adverse events of special interest. Results:Compared with tegileridine group, no significant change was found in the SPID 0-24h ( P>0.05), and the incidence of vomiting, PONV, proportion of patients requiring antiemetic medication, and total consumption of antiemetics were significantly decreased within the first 24 h after the initial dose in tegileridine group ( P<0.05). One treatment-emergent adverse event of Common Terminology Criteria for Adverse Events grade 3 or higher occurred in tegileridine group, while no treatment-emergent adverse events of Common Terminology Criteria for Adverse Events grade 3 or higher were found in anrikefon group. Among the adverse events of special interest, one case of respiratory depression and one case of cough occurred in tegileridine group, while one case of cough occurred in anrikefon group, with no respiratory depression. Conclusions:Anrikefon and tegileridine provide comparable analgesic efficacy for moderate-to-severe pain after abdominal surgery with general anesthesia. However, anrikefon exhibits an advantage in reducing the risk of PONV, with a superior safety profile.

Objective Many environmental epidemiological studies have shown the associations between short-term exposure of air pollution and daily deaths.However,the generally available population data only contain a small number of measured confounding factors,which is faced with the problem that a large number of unobserved confounding factors are not included in the model,resulting in biased estimates.The instrumental variable method can solve the problem of estimating the effects caused by unobserved confounders.In this paper,We used the instrumental variable method to estimate the effects of PM2.5 on daily mortality.Methods We collected daily PM2.5 concentrations,meteorological data,and nonaccidental daily deaths in a Chinese city from 2016 to 2019.We used boundary layer height and wind speed as instrumental variables to estimate the effects of PM2.5 on nonaccidental daily mortality.Negative exposure control was used to test the hypothesis of instrumental variables.Meanwhile time series bootstrap method was used to estimate confidence interval.We compared the results of the generalized additive model and instrumental variable method.Results The instrumental variable method showed that PM2.5 was significantly related to daily deaths.For every 10 μg/m3 increase of PM2.5 concentrations,the nonaccidental daily deaths increased by 0.94%(95%CI:0.39%～1.55%).Negative exposure control results showed no correlation between negative exposure and nonaccidental deaths(P=0.19),indicated that the aforementioned instrumental variable model was not affected by unmeasured and uncontrolled confounders.The traditional generalized additive model estimated that for every 10 μg/m3 increase in PM2.5 concentrations,the nonaccidental deaths would increase by 0.24%(95%CI:0.01%～0.47%).Conclusion The instrumental variable method estimated that PM2.5 concentrations were significantly correlated with the nonaccidental daily deaths.Boundary layer height and wind speed can be used as instrumental variables to estimate the effects of PM2.5 concentrations on nonaccidental deaths.

Objective To analyze the clinical features and surgical treatment strategies of T4a thyroid cancer.Methods We retrospectively analyzed patients with thyroid cancer treated in the Department of Head and Neck Surgery of Sichuan Cancer Hospital from January 2004 to May 2021.A total of 303 cases were included and statistically analyzed for pathological type,invaded organs,surgical approach,survival time,and overall survival rate.The postoperative survival curves of the patients were analyzed using the Kaplan Meier method.Results Of the 303 patients enrolled,53 patients were lost to follow-up,and the 1-year,3-year and 5-year overall survival rates were 98.4％(246/250),97.0％(224/231)and 90.2％(92/102),respectively.Of the 94 patients with recurrent laryngeal nerve invasion only,13 were lost to follow-up,and the 1-year,3-year and 5-year overall survival rates were 100％(81/81),98.7％(77/78)and 97.4％(38/39),respectively.There were 151 patients with invasion of recurrent laryngeal nerve and tracheal/laryngeal/esophageal nerve,31 of them were lost to follow-up,and the 1-year,3-year and 5-year overall survival rates were 96.7％(116/120),95.3％(101/106)and 82.2％(37/45),respectively.According to the survival curve analysis,the group with recurrent laryngeal nerve invasion only had an advantage in overall survival time over the group with recurrent laryngeal nerve and tracheal/laryngeal/esophageal invasion.Conclusion Surgical resection is supposed to be preferred for T4a thyroid cancer if there is a chance of surgery.A reasonable surgical strategy,radical surgery while preserving the vital tissues and organs,and one-stage repair and reconstruction can bring patients a better quality of life and prognosis.