OBJECTIVE To predict and validate the mechanisms of Chaiqi yigan granules （CQYG） improving hepatocellular carcinoma （HCC）. METHODS The signaling pathways of CQYG intervention in HCC were predicted using network pharmacology. A mice model of transplanted hepatocellular carcinoma was established by injecting H22 hepatoma cells into the axilla. Successfully modeled mice were randomly divided into model group （normal saline）， sorafenib group （positive control， 50 mg/kg）， and CQYG low-， medium- and high-dose groups （24.83， 49.66， 99.32 g/kg）， with 10 mice in each group. Mice in each group were administered the corresponding drug solution or normal saline intragastrically， once a day， for 14 consecutive days. After last administration， pathological morphological changes in the tumor tissues of mice were observed in each group. Immunohistochemical staining was performed to detect the expression of the nuclear proliferation antigen Ki-67 in tumor tissues of mice. Western blot assay was used to measure the expression of proteins related to epithelial-mesenchymal transition （EMT） [N-cadherin， E-cadherin， Vimentin， matrix metalloproteinase 7 （MMP7）] and the mitogen-activated protein kinase （MAPK） signaling pathway [p38 MAPK， phosphorylated p38 MAPK， c-Jun N-terminal kinase （JNK）， phosphorylated JNK， extracellular regulated protein kinase 1/2 （ERK1/2）， phosphorylated ERK1/2] in tumor tissue of mice. RESULTS Network pharmacology analysis revealed that metabolic pathways， pathways in cancer， and the MAPK signaling pathway were key signaling pathways through which CQYG exert their anti-hepatocellular carcinoma effects. In animal experiments， the tumor tissues of mice in the model group exhibited dense tumor cells and vigorous growth. Compared with model group， CQYG high-dose group showed a decreased density of tumor cells in the tumor tissues of mice. Moreover， the expression levels of Ki-67， N-cadherin， MMP7 and Vimentin proteins， along with the phosphorylation levels of ERK1/2 and JNK proteins， were all significantly reduced （ P ＜0.05）. The expression level of E-cadherin protein was significantly increased （ P ＜0.05）， the phosphorylation level of p38 MAPK protein was increased， the difference was not statistically significant （ P ＞0.05）. CONCLUSIONS CQYG can inhibit EMT by regulating the MAPK signaling pathway， thereby suppressing tumor cell invasion and metastasis and ultimately exerting a therapeutic effect in improving HCC.

ObjectiveThis study aimed to provide a scientific basis for the safety management of drinking water by assessing the contamination levels of 21 per- and polyfluoroalkyl substances (PFASs) in water samples collected from the Yangtze River Delta region. MethodsSolid-phase extraction coupled with ultra-high-performance liquid chromatography-tandem mass spectrometry (UHPLC-MS/MS) was employed to determine 21 PFASs in 6 source water samples, 20 treated water samples, and 60 tap water samples collected from 5 administrative districts in the Yangtze River Delta region across four seasons in 2023. The Kruskal-Wallis H test was used to compare differences in PFAS concentrations among the three types of water samples, followed by post hoc pairwise comparisons with Bonferroni correction. Removal efficiencies of target compounds were also calculated. ResultsAmong the 86 samples, 14 PFASs were detected, with 7 compounds detected in 100% of samples. The total concentrations of PFASs in source water, treated water, and tap water ranged from 63.85 to 190.54 ng·L^-¹, 55.75 to 774.36 ng·L^-¹, and 40.72 to 689.80 ng·L^-¹, respectively. The compounds with the highest median concentrations were perfluorobutanoic acid (PFBA) (40.19, 15.20, and 14.05 ng·L^-¹), perfluorohexanoic acid (PFHxA) (29.69, 32.62, and 34.24 ng·L^-¹), perfluorooctanoic acid (PFOA) (17.65, 39.45, and 36.83 ng·L^-¹), and potassium perfluorobutanesulfonic acid (PFBS) (16.99, 24.38, and 23.89 ng·L^-¹). Notably, concentrations of PFOA (7.71‒71.63 ng·L^-¹) and PFOS (-¹) in tap water were below the regulatory limits specified in GB 5749‒2022 Standards for Drinking Water Quality (PFOA: 0.000 08 mg·L^-1, PFOS: 0.000 04 mg·L^-1). Significant seasonal variations in PFAS concentrations were observed, with total PFAS levels in March being higher than those in September and November (adjusted P<0.001, and P=0.006), and levels in May being higher than those in September (adjusted P=0.037). The removal efficiencies of PFASs in the water treatment plants ranged from -14% to 34%. ConclusionPFAS contamination is widespread in source water, treated water, and tap water in the Yangtze River Delta region, and concentrations of PFOA and PFOS in all tested water samples comply with the Chinese national standards for drinking water quality (GB 5749‒2022). The use of PFASs should be continuously regulated to ensure water safety.

Background With the progression of industrialization, an increasing number of emerging contaminants are entering aquatic environments, posing significant threats to the safety of drinking water. Therefore, establishing a system for identifying unknown hazardous factors and implementing safety warning mechanisms for drinking water is of paramount importance. Among these efforts, non-target screening plays a critical role, but its effectiveness is largely constrained by the scope of coverage of sample pre-treatment methods. Objective To integrate modern chromatography/mass spectrometry techniques with advanced data mining methods to develop a non-discriminatory sample pre-treatment method for comprehensive enrichment of unknown contaminants in drinking water, laying a technical foundation for the discovery and identification of unknown organic hazardous factors in drinking water. Methods A non-discriminatory pre-treatment method based on supramolecular and solid-phase extraction was developed. The final target compounds including 333 pesticides, 194 pharmaceuticals and personal care products (PPCPs), and 59 per- and polyfluoroalkyl substances (PFASs) were used for optimizing the pre-treatment method, confirming its coverage. The impacts of different eluents on the absolute recovery rates of target compounds were compared to select the conditions with the highest recovery for sample pre-treatment. The effects of different supramolecular solvents and salt concentrations on target compound recovery were also evaluated to determine the most suitable solvent and salt concentration. Results The solid-phase extraction elution solvents, supramolecular extraction solvents, and salt concentrations were optimized based on the target compound recovery rates. The optimal recovery conditions were achieved using 2 mL methanol, 2 mL methanol (containing 1% formic acid), 2 mL ethyl acetate, 2 mL dichloromethane, hexanediol supramolecular solvent, and 426 mg salt. The detection method developed based on these conditions showed a good linear relationship for all target compounds in the range of 0.1-100.0 ng·mL⁻¹, with R² > 0.99. The method’s limit of detection ranged from 0.01 ng⁻¹ to 0.95 ng⁻¹, and 95% of target compounds were recovered in the range of 20%-120%, with relative standard deviation (RSD) less than 30%, indicating good precision. Conclusion The combined pre-treatment method of solid-phase extraction and supramolecular solvent extraction can effectively enrich contaminants in drinking water across low, medium, and high polarities, enabling broad-spectrum enrichment of diverse trace contaminants in drinking water. It provides technical support for broad-spectrum, high-throughput screening and identification of organic pollutants in drinking water, and also serves as a reference for establishing urban drinking water public safety warning systems.

This study aims to describe the population and regional distribution characteristics of smoking behavior among adult twins in the China Twin Registry (CNTR), as well as the concordance rates for smoking behavior in monozygotic and dizygotic twins, and estimate the heritability. The study population included adult twins in CNTR who had smoking questionnaire data. A random-effects regression model was used to describe the distribution of smoking behavior among different subgroups based on various characteristics. The concordance of smoking behavior between different zygosity groups was calculated, and heritability was estimated. A total of 28 444 twin pairs were included in this study, with an average age of (36.6±12.0) years. Among male twins, 41.2% were current smokers, while only 1.2% of females smoked. Higher smoking rates were observed among male smokers in the 50-59 age group ( z=23.0, P<0.001), northern regions ( z=2.9, P<0.01), rural areas ( z=-5.2, P<0.001), those who were divorced/widowed ( z=3.8, P<0.001), and first-born twins ( z=-4.3, P<0.001), while lower smoking rates were found in those with higher education ( z=-16.1, P<0.001) and unmarried individuals ( z=-16.0, P<0.001). The smoking concordance rate for male monozygotic twins was 69.6%, significantly higher than the 57.3% concordance rate for dizygotic twins ( χ 2=105.0, P<0.05). The heritability of smoking behavior in male twins was estimated at 28.9% (95% CI: 24.3%-33.4%). Stratified analyses showed differences in heritability across regions and age groups: the heritability in northern regions was 32.6% (95% CI: 27.3%-38.0%), higher than the 21.0% (95% CI: 12.4%-29.5%) observed in southern regions; the highest heritability of 35.1% (95% CI: 26.3%-43.9%) was found in the 18-29 age group, with heritability decreasing with age. In conclusion, the smoking rate and influencing factors in the twin population are similar to those in the general population, with unique characteristics, such as higher smoking rates in first-born twins. Genetic factors have a significant impact on smoking behavior.

Objective:To describe the distribution characteristics of alcohol consumption in adult twins in the Chinese National Twin Registry (CNTR), and further explore the influence of genetic factors on alcohol consumption in adult twins.Methods:The subjects of the study were twins registered by CNTR in 11 project areas across China from 2010 to 2018. A total of 56 966 twins (28 483 pairs) aged 18 years and above who answered questions about drinking behavior were included, and the random effect model was used to describe the population and regional distribution characteristics of alcohol consumption. Intra-pair analysis was performed to calculate the concordance rate and heritability of their alcohol consumption.Results:The age of all subjects was (36.6±12.0) years, and current drinkers accounted for 16.6% (9 461/56 966) of all subjects. In men, those aged 50-59 years, those in northern China, those living in rural area, those with low education level and those with high BMI, the proportions of current drinkers were higher. After excluding 468 pairs of twins who had stopped alcohol use and 21 764 pairs of twins who had no drink or had small amount drink, an intra-pair analysis was conducted in 4 929 pairs of same-sex twins, and found that the concordance rate of alcohol consumption was 64.0% (2 059/3 215) in monozygotic twins, and 52.6% (902/1 714) in dizygotic twins, the difference was significant ( P<0.001), and the heritability of alcohol consumption was 24.1% (95% CI: 18.9%- 29.3%). The further stratified analysis found that in southern men, the heritability was highest in those aged 40-49 years (36.1%, 95% CI: 21.6%-50.7%), while in northern men, the heritability was highest in those aged 50-59 years (34.2%, 95% CI: 18.1%-50.3%). Conclusions:In adult twins in China, there were population and regional differences in the distribution of alcohol consumption behavior, and alcohol consumption was influenced by genetic factors, and gender, age and region had potential modifying effects.

Objective:To explore the relationship between obesity indicators and DNA methylation clocks acceleration,and to analyze their temporal sequence.Methods:Data were obtained from two sur-veys conducted in 2013 and 2017-2018 by the Chinese National Twin Registry.Peripheral blood DNA methylation data were measured using the Illumina Infinium Human Methylation 450K BeadChip and EPIC BeadChip.DNA methylation clocks/acceleration metrics(GrimAA,PCGrimAA and Dunedin-PACE)were calculated using the DNA methylation online tool(https://dnamage.genetics.ucla.edu/)or R code provided by researchers.Obesity indicators included weight,body mass index(BMI),waist circumference,waist-hip ratio,and waist-height ratio.A total of 1 070 twin individuals were included in the cross-sectional analysis,comprising 378 monozygotic(MZ)twin pairs and 155 dizygotic(DZ)twin pairs for within-pair analysis.Mixed-effects models were used to examine the associations between obesity indicators and DNA methylation clocks,as well as their acceleration measures.The longitudinal analysis included 314 twin individuals,comprising 95 MZ twin pairs and 62 DZ twin pairs for within-pair analy-sis.Cross-lagged panel models were applied to further explore the temporal relationships between obesity and DNA methylation clock indicators.All analyses were conducted both in the full twin sample and separately within MZ and DZ twin pairs.Results:In the cross-sectional analysis population,monozygotic twins accounted for 71.0％,males for 68.0％,and the mean chronological age was(49.9±12.1)years.In the longitudinal analysis population,monozygotic twins accounted for 60.5％,males for 60.8％,with a mean baseline chronological age of(50.4±10.2)years and a mean follow-up duration of(4.6±0.6)years.Except for the waist-to-hip ratio,which was significantly higher at follow-up com-pared with baseline,no statistically significant differences were observed in the means of other obesity in-dicators between baseline and follow-up.Correlation analysis revealed that weight,BMI,waist circumfe-rence,waist-hip ratio(WHR),and waist-height ratio(WHtR)were positively correlated with Dunedin-PACE in all the twins,with WHtR showing the strongest association(β=0.21,95％CI:0.11 to 0.31).Weight and BMI were negatively associated with GrimAA(β=-0.03,95％CI:-0.05 to-0.01;β=-0.07,95％CI:-0.12 to-0.02),while weight was negatively associated with PCGrim-AA(β=-0.02,95％CI:-0.03 to 0.00).However,within-twin-pair analyses showed no statistically significant correlations.Cross-lagged panel model analysis indicated that higher baseline weight might lead to increased GrimAA at follow-up,while elevated baseline weight,BMI,and waist circumference might increase PCGrimAA.Higher baseline WHR was associated with increased DunedinPACE at follow-up.Conclusion:Obesity indicators correlate with DNA methylation clock acceleration metrics.Baseline obesity may influence changes in certain DNA methylation clock indicators over time,suggesting that obesity could exert long-term health effects by accelerating DNA methylation aging.However,these associations may be confounded by shared genetic or environmental factors among the twins.

This study aims to describe the population and regional distribution characteristics of smoking behavior among adult twins in the China Twin Registry (CNTR), as well as the concordance rates for smoking behavior in monozygotic and dizygotic twins, and estimate the heritability. The study population included adult twins in CNTR who had smoking questionnaire data. A random-effects regression model was used to describe the distribution of smoking behavior among different subgroups based on various characteristics. The concordance of smoking behavior between different zygosity groups was calculated, and heritability was estimated. A total of 28 444 twin pairs were included in this study, with an average age of (36.6±12.0) years. Among male twins, 41.2% were current smokers, while only 1.2% of females smoked. Higher smoking rates were observed among male smokers in the 50-59 age group ( z=23.0, P<0.001), northern regions ( z=2.9, P<0.01), rural areas ( z=-5.2, P<0.001), those who were divorced/widowed ( z=3.8, P<0.001), and first-born twins ( z=-4.3, P<0.001), while lower smoking rates were found in those with higher education ( z=-16.1, P<0.001) and unmarried individuals ( z=-16.0, P<0.001). The smoking concordance rate for male monozygotic twins was 69.6%, significantly higher than the 57.3% concordance rate for dizygotic twins ( χ 2=105.0, P<0.05). The heritability of smoking behavior in male twins was estimated at 28.9% (95% CI: 24.3%-33.4%). Stratified analyses showed differences in heritability across regions and age groups: the heritability in northern regions was 32.6% (95% CI: 27.3%-38.0%), higher than the 21.0% (95% CI: 12.4%-29.5%) observed in southern regions; the highest heritability of 35.1% (95% CI: 26.3%-43.9%) was found in the 18-29 age group, with heritability decreasing with age. In conclusion, the smoking rate and influencing factors in the twin population are similar to those in the general population, with unique characteristics, such as higher smoking rates in first-born twins. Genetic factors have a significant impact on smoking behavior.

Objective:To explore the relationship between obesity indicators and DNA methylation clocks acceleration,and to analyze their temporal sequence.Methods:Data were obtained from two sur-veys conducted in 2013 and 2017-2018 by the Chinese National Twin Registry.Peripheral blood DNA methylation data were measured using the Illumina Infinium Human Methylation 450K BeadChip and EPIC BeadChip.DNA methylation clocks/acceleration metrics(GrimAA,PCGrimAA and Dunedin-PACE)were calculated using the DNA methylation online tool(https://dnamage.genetics.ucla.edu/)or R code provided by researchers.Obesity indicators included weight,body mass index(BMI),waist circumference,waist-hip ratio,and waist-height ratio.A total of 1 070 twin individuals were included in the cross-sectional analysis,comprising 378 monozygotic(MZ)twin pairs and 155 dizygotic(DZ)twin pairs for within-pair analysis.Mixed-effects models were used to examine the associations between obesity indicators and DNA methylation clocks,as well as their acceleration measures.The longitudinal analysis included 314 twin individuals,comprising 95 MZ twin pairs and 62 DZ twin pairs for within-pair analy-sis.Cross-lagged panel models were applied to further explore the temporal relationships between obesity and DNA methylation clock indicators.All analyses were conducted both in the full twin sample and separately within MZ and DZ twin pairs.Results:In the cross-sectional analysis population,monozygotic twins accounted for 71.0％,males for 68.0％,and the mean chronological age was(49.9±12.1)years.In the longitudinal analysis population,monozygotic twins accounted for 60.5％,males for 60.8％,with a mean baseline chronological age of(50.4±10.2)years and a mean follow-up duration of(4.6±0.6)years.Except for the waist-to-hip ratio,which was significantly higher at follow-up com-pared with baseline,no statistically significant differences were observed in the means of other obesity in-dicators between baseline and follow-up.Correlation analysis revealed that weight,BMI,waist circumfe-rence,waist-hip ratio(WHR),and waist-height ratio(WHtR)were positively correlated with Dunedin-PACE in all the twins,with WHtR showing the strongest association(β=0.21,95％CI:0.11 to 0.31).Weight and BMI were negatively associated with GrimAA(β=-0.03,95％CI:-0.05 to-0.01;β=-0.07,95％CI:-0.12 to-0.02),while weight was negatively associated with PCGrim-AA(β=-0.02,95％CI:-0.03 to 0.00).However,within-twin-pair analyses showed no statistically significant correlations.Cross-lagged panel model analysis indicated that higher baseline weight might lead to increased GrimAA at follow-up,while elevated baseline weight,BMI,and waist circumference might increase PCGrimAA.Higher baseline WHR was associated with increased DunedinPACE at follow-up.Conclusion:Obesity indicators correlate with DNA methylation clock acceleration metrics.Baseline obesity may influence changes in certain DNA methylation clock indicators over time,suggesting that obesity could exert long-term health effects by accelerating DNA methylation aging.However,these associations may be confounded by shared genetic or environmental factors among the twins.

Objective:To describe the distribution characteristics of alcohol consumption in adult twins in the Chinese National Twin Registry (CNTR), and further explore the influence of genetic factors on alcohol consumption in adult twins.Methods:The subjects of the study were twins registered by CNTR in 11 project areas across China from 2010 to 2018. A total of 56 966 twins (28 483 pairs) aged 18 years and above who answered questions about drinking behavior were included, and the random effect model was used to describe the population and regional distribution characteristics of alcohol consumption. Intra-pair analysis was performed to calculate the concordance rate and heritability of their alcohol consumption.Results:The age of all subjects was (36.6±12.0) years, and current drinkers accounted for 16.6% (9 461/56 966) of all subjects. In men, those aged 50-59 years, those in northern China, those living in rural area, those with low education level and those with high BMI, the proportions of current drinkers were higher. After excluding 468 pairs of twins who had stopped alcohol use and 21 764 pairs of twins who had no drink or had small amount drink, an intra-pair analysis was conducted in 4 929 pairs of same-sex twins, and found that the concordance rate of alcohol consumption was 64.0% (2 059/3 215) in monozygotic twins, and 52.6% (902/1 714) in dizygotic twins, the difference was significant ( P<0.001), and the heritability of alcohol consumption was 24.1% (95% CI: 18.9%- 29.3%). The further stratified analysis found that in southern men, the heritability was highest in those aged 40-49 years (36.1%, 95% CI: 21.6%-50.7%), while in northern men, the heritability was highest in those aged 50-59 years (34.2%, 95% CI: 18.1%-50.3%). Conclusions:In adult twins in China, there were population and regional differences in the distribution of alcohol consumption behavior, and alcohol consumption was influenced by genetic factors, and gender, age and region had potential modifying effects.

Background Per- and polyfluoroalkyl substances (PFAS) are a class of persistent organic pollutants that possess potential toxicity to the human body. The production and utilization of diverse emerging PFAS have resulted in widespread human exposure. Therefore, it is imperative to establish a quantitative methodology encompassing a wide range of PFAS for a comprehensive assessment of human exposure to these compounds. Objective To establish a high-throughput quantitative method for the simultaneous determination of 53 PFAS in human serum based on ultra-high-performance liquid chromatography-Q Exactive high resolution mass spectrometry (UPLC-Q Exactive HRMS). Methods The extraction recoveries of hydrophilic-lipophilic balance (HLB) column, weak anionexchange (WAX) column, and 96-well WAX μElution plate were compared to select the SPE column with the highest recovery. The retention time and peak shape of the target compounds were compared between ACQUITY UPLC BEH C₁₈ column and Accucore aQ column, and the more cost-effective column was chosen. The effects of adding different levels of ammonium formate (0, 2, 5 and 10 mmol·L⁻¹) in mobile phase on peak shape and target response were compared to determine the optimal buffer salt concentration. The optimal spray voltage was obtained by comparing −2 kV and −4 kV. The proposed method was validated from the aspects of selectivity, standard curve, limits of detection, precision, accuracy, and matrix effect. The method was applied to 142 umbilical serum samples. Results The best recovery rate (64%-118%) was achieved by using 96-well WAX μElution plate. The optimal separation and peak shape were obtained by utilizing Accucore aQ column with H₂O-methanol (containing 5 mmol·L⁻¹ ammonium formate) as the mobile phase. Less in-source collision and better target response were observed when the spray voltage was set to −2 kV. All target analytes had a good linearity, with R² > 0.99. The limits of detection ranged from 0.01 to 0.50 μg·L⁻¹, and the recovery ranged from 69% to 127% with the precision less than 26%. A total of 31 PFAS were detected in the 142 actual samples, among which 14 PFAS had a detection frequency over 50%. Perfluorooctanoic acid showed the highest median concentration of 4.16 μg·L⁻¹, followed by 6:2 chlorinated polyfluorinated ether sulfonate and perfluorooctane sulfonates (3.50 μg·L⁻¹ and 1.59 μg·L⁻¹, respectively). Conclusion In this study, we establish a UPLC-Q Excative HRMS method for simutanious determination of 53 PFAS concentrations in serum. This method has the advantages of wide coverage of PFAS, good selectivity, and easy operation, and is suitable for biological detection with a large sample size.