Objective:To construct an in situ pancreatic tumor model in mice and explore the mechanism of neuromedin U (NMU) remodeling the metabolic microenvironment of pancreatic tumors via the Hedgehog signaling pathway.Methods:C57BL/6 NMU -/- heterozygous mice were bred in one cage with a female-to-male ratio of 2:1. The tail tissues of offspring rats were taken, and knockout primers were designed according to the sequence structure of NMU gene. The C57BL/6 NMU -/- mouse genotypes were identified by polymerase chain reaction (PCR) and agarose gel electrophoresis. Five C57BL/6 NMU -/- homozygotes and five wild type mice aged 8 to 9 weeks were selected, and Panc02 cell suspension of pancreatic cancer in logarithmic growth phase was injected into pancreatic tail tissue to construct tumor bearing mice model of pancreatic tumor in situ. After 3 weeks of tumor loading, flow cytometry was used to detect the abundance changes of CD8 + T cells in the spleen tissues of two groups of mice. GSE102238 (data published in August 2017) and GSE62452 (data published in July 2016) datasets were download from the Gene Expression Omnibus (GEO) database, including 119 samples of pancreatic cancer tissue and 111 samples of normal pancreatic tissue adjacent to cancer, and the differential expression factors of the two groups of samples were screened. R language limma package was used to analyze the differential expression of NMU mRNA in pancreatic cancer tissues and adjacent normal tissues; the relative expression level of NMU mRNA in pancreatic cancer patients of different ages and tumor grades was analyzed by ggpubr package. The relative expression level data of NMU mRNA in pancreatic tumor patients were extracted, and the patients were divided into high expression group (>2.48) and low expression group (<2.48) based on the median value (2.48). The CIBERSORT algorithm was used to calculate the difference in the content of infiltrating immune cells in pancreatic tumors between the two groups of patients. The overall survival of patients with different relative expression levels of NMU mRNA in the OncoLnc database (59 cases in the high expression group and 59 cases in the low expression group) was compared. The data of 178 patients with pancreatic cancer in TCGA-PAAD of The Cancer Genome Atlas (TCGA) database (updated in March 2024) were download. According to the median value of relative expression of NMU mRNA (2.740), the patients were divided into the high expression group (>2.740) and the low expression group (<2.740), with 89 cases in each group. GSEA software was used to analyze the biological functions and pathways in which the genes enriched significantly. The molecules and key targets for molecular docking were explored using the large molecule protein interaction tool PDBePISA. Spearman correlation analysis was performed using R language data packets. Results:The results of C57BL/6 NMU -/- mouse genotypes analyzed by using the gel electrophoresis pattern of PCR amplification products of tail tissues showed that the wild type was one 380 bp electrophoretic band, the homozygous type was one 450 bp electrophoretic band, and the heterozygous type was 450 bp and 380 bp electrophoretic bands. After identifying and screening pure strains mouse, reproduction was carried out, and C57BL/6 NMU -/- homozygous mice were successfully obtained. The results of flow cytometry analysis showed that the proportion of CD8 + T cells in the spleen tissues of C57BL/6 NMU -/- wild-type and homozygous pancreatic tumor bearing mice was (30.38±0.37)% and (37.00±0.48)%, with a statistically significant difference between the two groups ( t = 9.79, P < 0.001). The absolute values of CD8 + T cells were (8.36±0.27)×10 4 and (10.20±0.28)×10 4, respectively, and the difference was statistically significant ( t = 3.71, P = 0.013). In the GEO database GSE102238 and GSE62452 datasets, there were differentially expressed genes in pancreatic cancer tissues compared with adjacent tissues, including 219 up-regulated genes and 325 down-regulated genes in pancreatic cancer tissues; among them, NMU was an upregulated differentially expressed gene. The relative expression of NMU mRNA in patients aged > 65 year or grade G3-G4 was higher than that in patients aged ≤ 65 years or grade G1-G2, and the differences were statistically significant (both P < 0.05). The analysis results of CIBERSORT algorithm showed that the expression abundance of CD8 + T cells, initial CD4 + T cells, activated memory CD4 + T cells, and γδ T cells in the high expression group of NMU was lower than that in the low expression group (all P < 0.05). The overall survival of the high expression group of NMU was worse than that of the low expression group in the OncoLnc database ( P = 0.010). GSEA enrichment analysis and Spearman correlation analysis revealed significant changes in the Hedgehog signaling pathway, biosynthesis of unsaturated fatty acids and pyrimidine nucleotide metabolism, which affected tumor metabolic remodeling. Conclusions:NMU may remodel the tumor microenvironment of pancreatic cancer by regulating Hedgehog signaling pathway.

Objective:To explore the spousal correlations of total cholesterol(TC),total triglyceride(TG),low-density lipoprotein cholesterol(LDL-C),and high-density lipoprotein cholesterol(HDL-C),and to investigate the reasons behind these spousal correlations.Methods:Participants and data were from the baseline survey of family-based cohort studies in Fangshan,Beijing and Tulou,Fujian.The ori-gin of spousal correlations were explored from perspectives of convergence,assortative mating,social ho-mogamy.Pearson's correlation and generalized linear models(GLM)were used to estimate the spousal correlation.Convergence was assessed by Pearson's correlation between the phenotypic differences be-tween couples and the duration of marriage,with GLM used for further validation.Pearson's correlation of genetic risk scores(GRS)and couple-specific Mendelian randomization(MR)were calculated to assess the genetic correlation and possible causal relationships between spouses.Two-independent-sample t-tests were used to compare GRS consistency across subgroups divided by education attainment,couple-specific MR and Q statistics used to test assortative mating in subgroups and intergroup differences.Results:In the study,342 couples(287 couples from Fangshan and 55 couples from Fujian)were included,with the average age of(64.91±8.76)years.Spousal correlations of TC,TG,HDL-C,and LDL-C showed statistically significant associations both before and after adjusting for covariates,with effect sizes of 0.229(95％CI:0.125-0.327),0.257(95％CI:0.155-0.354),0.179(95％CI:0.074-0.280),and 0.181(95％CI:0.076-0.282).For convergence,for each additional year of marriage,ΔTC increased by 0.016 mmol/L(95％CI:0.001-0.033 mmol/L),and ΔLDL-C increased by 0.017 mmol/L(95％CI:0.002-0.031 mmol/L).For assortative mating,GRS correlations and results of couple specific MR didn't show any statistical significance.For social homogamy,no differences in GRS or assortative mating were found between subgroups stratified by education attainment.Conclusion:The blood lipid in participants exhibit spousal phenotypic correlations,however,no effects of convergence,assortative mating or social homogamy were observed.More independent studies with larger sample sizes are warranted to further validate these findings in the future.

AIM:To explore the impact of chronic stress on postoperative cognitive dysfunction in rats and to elucidate the mechanistic link to histone deacetylase 2(HDAC2).METHODS:A repeated social defeat stress model and a prolonged isoflurane anesthesia model were established in mice.The rats were randomly assigned to four groups:control(Ctrl)group,isoflurane anesthesia(Iso)group,chronic social defeat stress(RSDS)group,and chronic social de-feat stress combined with isoflurane anesthesia(RSDS+Iso)group.Anxiety-like behaviors were evaluated using the social avoidance test and the novelty-suppressed feeding test.Cognitive function was assessed through the novel object recogni-tion test and the Morris water maze.Plasma corticosterone levels were measured via enzyme-linked immunosorbent assay(ELISA).Primary hippocampal neurons were isolated from fetal mouse hippocampi and classified into four groups:con-trol group,chronic stress combined with prolonged isoflurane anesthesia(Cort+Iso)group,CAY-10683 intervention(CAY),and CAY-10683 treatment(CAY+Cort+Iso)group.Cell viability was determined using CCK-8 assay,and pro-tein expression levels of brain-derived neurotrophic factor(BDNF)and HDAC2 were analyzed by Western blot.RE-SULTS:The RSDS mouse model was successfully established,with ELISA results indicating a significant increase in plasma corticosterone levels in mice subjected to chronic stress combined with prolonged isoflurane anesthesia.Behavioral assessments and Western blot analyses revealed that mice exposed to prolonged isoflurane anesthesia following chronic stress showed marked declines in cognitive function and hippocampal BDNF protein expression levels.Additionally,chronic stress significantly elevated HDAC2 protein expression in the hippocampi of mice undergoing prolonged isoflurane anesthesia.Treatment with an HDAC2 inhibitor reduced HDAC2 protein expression in primary hippocampal neurons sub-jected to chronic stress combined with prolonged isoflurane anesthesia,concurrently increasing BDNF protein expression levels.CONCLUSION:Chronic stress significantly worsens postoperative cognitive dysfunction induced by prolonged isoflurane anesthesia,associated with increased HDAC2 protein expression in the hippocampus.Inhibition of HDAC2 ef-fectively counteracts the reduction in BDNF,a protein crucial for cognitive function,caused by the combination of chronic stress and prolonged isoflurane anesthesia.

Objective:To explore the spousal correlations of total cholesterol(TC),total triglyceride(TG),low-density lipoprotein cholesterol(LDL-C),and high-density lipoprotein cholesterol(HDL-C),and to investigate the reasons behind these spousal correlations.Methods:Participants and data were from the baseline survey of family-based cohort studies in Fangshan,Beijing and Tulou,Fujian.The ori-gin of spousal correlations were explored from perspectives of convergence,assortative mating,social ho-mogamy.Pearson's correlation and generalized linear models(GLM)were used to estimate the spousal correlation.Convergence was assessed by Pearson's correlation between the phenotypic differences be-tween couples and the duration of marriage,with GLM used for further validation.Pearson's correlation of genetic risk scores(GRS)and couple-specific Mendelian randomization(MR)were calculated to assess the genetic correlation and possible causal relationships between spouses.Two-independent-sample t-tests were used to compare GRS consistency across subgroups divided by education attainment,couple-specific MR and Q statistics used to test assortative mating in subgroups and intergroup differences.Results:In the study,342 couples(287 couples from Fangshan and 55 couples from Fujian)were included,with the average age of(64.91±8.76)years.Spousal correlations of TC,TG,HDL-C,and LDL-C showed statistically significant associations both before and after adjusting for covariates,with effect sizes of 0.229(95％CI:0.125-0.327),0.257(95％CI:0.155-0.354),0.179(95％CI:0.074-0.280),and 0.181(95％CI:0.076-0.282).For convergence,for each additional year of marriage,ΔTC increased by 0.016 mmol/L(95％CI:0.001-0.033 mmol/L),and ΔLDL-C increased by 0.017 mmol/L(95％CI:0.002-0.031 mmol/L).For assortative mating,GRS correlations and results of couple specific MR didn't show any statistical significance.For social homogamy,no differences in GRS or assortative mating were found between subgroups stratified by education attainment.Conclusion:The blood lipid in participants exhibit spousal phenotypic correlations,however,no effects of convergence,assortative mating or social homogamy were observed.More independent studies with larger sample sizes are warranted to further validate these findings in the future.

AIM:To explore the impact of chronic stress on postoperative cognitive dysfunction in rats and to elucidate the mechanistic link to histone deacetylase 2(HDAC2).METHODS:A repeated social defeat stress model and a prolonged isoflurane anesthesia model were established in mice.The rats were randomly assigned to four groups:control(Ctrl)group,isoflurane anesthesia(Iso)group,chronic social defeat stress(RSDS)group,and chronic social de-feat stress combined with isoflurane anesthesia(RSDS+Iso)group.Anxiety-like behaviors were evaluated using the social avoidance test and the novelty-suppressed feeding test.Cognitive function was assessed through the novel object recogni-tion test and the Morris water maze.Plasma corticosterone levels were measured via enzyme-linked immunosorbent assay(ELISA).Primary hippocampal neurons were isolated from fetal mouse hippocampi and classified into four groups:con-trol group,chronic stress combined with prolonged isoflurane anesthesia(Cort+Iso)group,CAY-10683 intervention(CAY),and CAY-10683 treatment(CAY+Cort+Iso)group.Cell viability was determined using CCK-8 assay,and pro-tein expression levels of brain-derived neurotrophic factor(BDNF)and HDAC2 were analyzed by Western blot.RE-SULTS:The RSDS mouse model was successfully established,with ELISA results indicating a significant increase in plasma corticosterone levels in mice subjected to chronic stress combined with prolonged isoflurane anesthesia.Behavioral assessments and Western blot analyses revealed that mice exposed to prolonged isoflurane anesthesia following chronic stress showed marked declines in cognitive function and hippocampal BDNF protein expression levels.Additionally,chronic stress significantly elevated HDAC2 protein expression in the hippocampi of mice undergoing prolonged isoflurane anesthesia.Treatment with an HDAC2 inhibitor reduced HDAC2 protein expression in primary hippocampal neurons sub-jected to chronic stress combined with prolonged isoflurane anesthesia,concurrently increasing BDNF protein expression levels.CONCLUSION:Chronic stress significantly worsens postoperative cognitive dysfunction induced by prolonged isoflurane anesthesia,associated with increased HDAC2 protein expression in the hippocampus.Inhibition of HDAC2 ef-fectively counteracts the reduction in BDNF,a protein crucial for cognitive function,caused by the combination of chronic stress and prolonged isoflurane anesthesia.

Objective:To construct an in situ pancreatic tumor model in mice and explore the mechanism of neuromedin U (NMU) remodeling the metabolic microenvironment of pancreatic tumors via the Hedgehog signaling pathway.Methods:C57BL/6 NMU -/- heterozygous mice were bred in one cage with a female-to-male ratio of 2:1. The tail tissues of offspring rats were taken, and knockout primers were designed according to the sequence structure of NMU gene. The C57BL/6 NMU -/- mouse genotypes were identified by polymerase chain reaction (PCR) and agarose gel electrophoresis. Five C57BL/6 NMU -/- homozygotes and five wild type mice aged 8 to 9 weeks were selected, and Panc02 cell suspension of pancreatic cancer in logarithmic growth phase was injected into pancreatic tail tissue to construct tumor bearing mice model of pancreatic tumor in situ. After 3 weeks of tumor loading, flow cytometry was used to detect the abundance changes of CD8 + T cells in the spleen tissues of two groups of mice. GSE102238 (data published in August 2017) and GSE62452 (data published in July 2016) datasets were download from the Gene Expression Omnibus (GEO) database, including 119 samples of pancreatic cancer tissue and 111 samples of normal pancreatic tissue adjacent to cancer, and the differential expression factors of the two groups of samples were screened. R language limma package was used to analyze the differential expression of NMU mRNA in pancreatic cancer tissues and adjacent normal tissues; the relative expression level of NMU mRNA in pancreatic cancer patients of different ages and tumor grades was analyzed by ggpubr package. The relative expression level data of NMU mRNA in pancreatic tumor patients were extracted, and the patients were divided into high expression group (>2.48) and low expression group (<2.48) based on the median value (2.48). The CIBERSORT algorithm was used to calculate the difference in the content of infiltrating immune cells in pancreatic tumors between the two groups of patients. The overall survival of patients with different relative expression levels of NMU mRNA in the OncoLnc database (59 cases in the high expression group and 59 cases in the low expression group) was compared. The data of 178 patients with pancreatic cancer in TCGA-PAAD of The Cancer Genome Atlas (TCGA) database (updated in March 2024) were download. According to the median value of relative expression of NMU mRNA (2.740), the patients were divided into the high expression group (>2.740) and the low expression group (<2.740), with 89 cases in each group. GSEA software was used to analyze the biological functions and pathways in which the genes enriched significantly. The molecules and key targets for molecular docking were explored using the large molecule protein interaction tool PDBePISA. Spearman correlation analysis was performed using R language data packets. Results:The results of C57BL/6 NMU -/- mouse genotypes analyzed by using the gel electrophoresis pattern of PCR amplification products of tail tissues showed that the wild type was one 380 bp electrophoretic band, the homozygous type was one 450 bp electrophoretic band, and the heterozygous type was 450 bp and 380 bp electrophoretic bands. After identifying and screening pure strains mouse, reproduction was carried out, and C57BL/6 NMU -/- homozygous mice were successfully obtained. The results of flow cytometry analysis showed that the proportion of CD8 + T cells in the spleen tissues of C57BL/6 NMU -/- wild-type and homozygous pancreatic tumor bearing mice was (30.38±0.37)% and (37.00±0.48)%, with a statistically significant difference between the two groups ( t = 9.79, P < 0.001). The absolute values of CD8 + T cells were (8.36±0.27)×10 4 and (10.20±0.28)×10 4, respectively, and the difference was statistically significant ( t = 3.71, P = 0.013). In the GEO database GSE102238 and GSE62452 datasets, there were differentially expressed genes in pancreatic cancer tissues compared with adjacent tissues, including 219 up-regulated genes and 325 down-regulated genes in pancreatic cancer tissues; among them, NMU was an upregulated differentially expressed gene. The relative expression of NMU mRNA in patients aged > 65 year or grade G3-G4 was higher than that in patients aged ≤ 65 years or grade G1-G2, and the differences were statistically significant (both P < 0.05). The analysis results of CIBERSORT algorithm showed that the expression abundance of CD8 + T cells, initial CD4 + T cells, activated memory CD4 + T cells, and γδ T cells in the high expression group of NMU was lower than that in the low expression group (all P < 0.05). The overall survival of the high expression group of NMU was worse than that of the low expression group in the OncoLnc database ( P = 0.010). GSEA enrichment analysis and Spearman correlation analysis revealed significant changes in the Hedgehog signaling pathway, biosynthesis of unsaturated fatty acids and pyrimidine nucleotide metabolism, which affected tumor metabolic remodeling. Conclusions:NMU may remodel the tumor microenvironment of pancreatic cancer by regulating Hedgehog signaling pathway.

The complex chemical composition and limited research ideas of traditional Chinese medicine （TCM） have led to the unclear material basis and mechanism of the medicinal effects， which is a common problem hindering the modernization of TCM in China. The introduction of computer virtual technology has provided a new perspective for TCM research. In this study， we established the research method of structure-activity omics to study the relationships between the structures and effects of different compounds in TCM based on the chemical structures of TCM components and to analyze and predict the material basis and multitarget synergistic mechanism of TCM. Furthermore， a structure-activity omics study was carried out with the anti-inflammatory and analgesic effects of Qizhi Weitong granules as an example. This study provides support for screening the pharmacodynamic components and analyzing the active ingredients of TCM and gives insights into the research on the material basis and mechanism of compound efficacy and the development of lead compounds of TCM， thus promoting the modern research and the innovative development of TCM.

ObjectiveTo explain the anti-inflammatory and analgesic effects of Corydalis Rhizoma by the means of structure-activity omics. MethodOn the basis of the previous in vitro screening study， we studied the in vivo efficacy of the alkaloids in Corydalis Rhizoma. With the targets as a bridge， the structures of chemical components in Corydalis Rhizoma were connected with the efficacy. The molecular docking of the alkaloids in Corydalis Rhizoma with the targets of inflammation and pain was carried out. According to the docking scores and the differences in the structural nucleus of Corydalis Rhizoma alkaloids， a study of structure-activity omics was carried out to summarize the rules of their connection. ResultThe alkaloids in Corydalis Rhizoma had good anti-inflammatory and analgesic effects in vivo， involving 53 chemical components and 73 targets. There were 3 074 targets associated with inflammation and pain， and 42 targets of direct action were shared by the chemical components and the disease. The protein-protein interaction （PPI） and molecular docking analysis predicted that the main active components of Corydalis Rhizoma were tetrahydropalmatine and palmatine， and the core targets were prostaglandin endoperoxide synthase 2 （PTGS2）， glutamate receptor metabotropic 5 （GRM5）， estrogen receptor 1 （ESR1）， solute carrier family 6 member 4 （SLC6A4）， and fusion oncoproteins （FOS）. According to the differences of mother nucleus， the 53 alkaloid components of Corydalis Rhizoma were classified into 8 categories， including protoberberine， berberine， and aporphine， which had high binding affinities with PTGS2， GRM5 and other targets. The relationship between the structures of Corydalis Rhizoma alkaloids and docking scores in each group showed the same law. In protoberberine， appropriate substituents with hydroxyl， alkoxy or methyl groups on the A and D rings of the parent ring were conducive to enhancing the binding activities with the two targets. In berberine， the structure containing a methyl group on position 13 had strong binding affinities with the two targets. It is hypothesized that the methyl fragment changes the binding mode between the component structure and amino acid residues， which greatly improves the binding affinity. ConclusionThis study employs the method of structure-activity omics to analyze the material basis for the anti-inflammatory and analgesic effects of alkaloids in Corydalis Rhizoma， and the structure-activity omics provides new ideas for revealing the pharmacodynamic substances of traditional Chinese medicine.

Objective:To evaluate oncological and reproductive outcomes of women ≤40 years undergoing fertility-sparing surgery (FSS) for stage Ⅱ or Ⅲ borderline ovarian tumor (BOT).Methods:The patients with BOT and ≤40 years old with stage Ⅱ-Ⅲ BOT who underwent FSS enrolled from the First Affiliated Hospital of Zhengzhou University between January 2011 and March 2023 were analyzed retrospectively. The clinical data and follow-up results were obtained and analyzed. The univariate and multivariate Cox proportional hazard regression analysis were used to explore high-risk factors associated with prognosis. Additionally, pregnancy outcomes were also analyzed.Results:(1) A total of 79 patients with stage Ⅱ-Ⅲ BOT who have been treated with FSS were conducted, with an average age of (27.5±6.7) years old. The median tumor maximum diameter were 10.4 cm (range: 4.8-90.0 cm). The International Federation of Gynecology and Obstetrics (FIGO) stage was stage Ⅱ in 45 cases and stage Ⅲ in 34 cases. According to the pathological types, there were 48 cases of serous tumor, 21 cases of mucinous tumor, 1 case of endometrioid tumor, and 9 cases of mixed types. There were 41 cases of unilateral ovarian involvement, 38 cases of bilateral ovarian involvement. There were 5 cases of microinvasion, 17 cases of micropapillary subtype. Extra-ovarian invasive implants were found in 5 cases, and there were 31 cases of merged ascites. (2) Tumor outcomes: the median follow-up time from primary cytoreduction were 58 months (range: 8-146 months). At the end of the observation period, 24 cases (30%, 24/79) recurred, among them 5 cases had two recurrences and 2 cases had three recurrences. There were 2 cases (3%, 2/79) of death and 1 case (1%, 1/79) of survival with tumor. (3) Analysis of prognostic risk factors: the results of univariate analysis showed that mucinous tumor, tumor maximum diameter >13.15 cm, FIGO stage Ⅲ, merged ascites, micropapillary subtype, invasive implantation, and bilateral ovarian involvement were independent risk factors (all P<0.05) for disease-free survival (DFS). FIGO stage Ⅲ ( HR=4.555, 95% CI: 1.525-13.607; P=0.007) and micropapillary subtype ( HR=2.396, 95% CI: 1.003-5.725; P=0.049) were found to be related to DFS through the multivariable Cox proportional hazard regression analysis. (4) Pregnancy outcomes: among the patients with fertility intentions 36 cases (46%,36/79), 29 cases (81%, 29/36) had successful pregnancies, and 27 cases (75%, 27/36) had successful births. Conclusions:Patients with stage Ⅱ-Ⅲ BOT underwent FSS have favorable survival and pregnancy rates. Micropapillary subtypes and FIGO staging (stage Ⅲ) are the significant risk factors of DFS.

Objective:To explore the association between polymorphisms of transforming growth factor-β(TGF-β)signaling pathway and non-syndromic cleft lip with or without cleft palate(NSCL/P)among Asian populations,while considering gene-gene interaction and gene-environment interaction.Methods:A total of 1 038 Asian NSCL/P case-parent trios were ascertained from an international consortium,which conducted a genome-wide association study using a case-parent trio design to investigate the genes affec-ting risk to NSCL/P.After stringent quality control measures,343 single nucleotide polymorphism(SNP)spanning across 10 pivotal genes in the TGF-β signaling pathway were selected from the original genome-wide association study(GWAS)dataset for further analysis.The transmission disequilibrium test(TDT)was used to test for SNP effects.The conditional Logistic regression models were used to test for gene-gene interaction and gene-environment interaction.Environmental factors collected for the study in-cluded smoking during pregnancy,passive smoking during pregnancy,alcohol intake during pregnancy,and vitamin use during pregnancy.Due to the low rates of exposure to smoking during pregnancy and al-cohol consumption during pregnancy(＜3％),only the interaction between maternal smoking during pregnancy and multivitamin supplementation during pregnancy was analyzed.The threshold for statistical significance was rigorously set at P=1.46 × 10-4,applying Bonferroni correction to account for multiple testing.Results:A total of 23 SNPs in 4 genes yielded nominal association with NSCL/P(P＜0.05),but none of these associations was statistically significant after Bonferroni's multiple test correction.How-ever,there were 6 pairs of SNPs rs4939874(SMAD2)and rs1864615(TGFBR2),rs2796813(TGFB2)and rs2132298(TGFBR2),rs4147358(SMAD3)and rs1346907(TGFBR2),rs4939874(SMAD2)and rs1019855(TGFBR2),rs4939874(SMAD2)and rs12490466(TGFBR2),rs2009112(TGFB2)and rs4075748(TGFBR2)showed statistically significant SNP-SNP interaction(P＜1.46 × 10-4).In contrast,the analysis of gene-environment interactions did not yield any significant results after being cor-rected by multiple testing.Conclusion:The comprehensive evaluation of SNP associations and interac-tions within the TGF-β signaling pathway did not yield any direct associations with NSCL/P risk in Asian populations.However,the significant gene-gene interactions identified suggest that the genetic architec-ture influencing NSCL/P risk may involve interactions between genes within the TGF-β signaling path-way.These findings underscore the necessity for further investigations to unravel these results and further explore the underlying biological mechanisms.