1.Bioinformatics-based analysis of autophagy-related genes and prediction of potential Chinese medicines in diabetic kidney disease
Yufeng XING ; Zining PENG ; Chaoyang YE
Digital Chinese Medicine 2025;8(1):90-99
[Objective] :
To predict the autophagy-related pathogenesis and key diagnostic genes of diabetic kidney disease (DKD) through bioinformatics analysis, and to identify related Chinese medicines.
[Methods] :
Data from sequencing microarrays GSE30528, GSE30529, and GSE1009 in the Gene Expression Omnibus (GEO) were employed. Differentially expressed genes (DEGs) with adjusted P < 0.05 from GSE30528 and GSE30529 were identified. Combining these DEGs with the human autophagy gene database, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, and protein-protein interaction (PPI) network analysis were conducted on the obtained DKD autophagy-related genes. Subsequently, the least absolute shrinkage and selection operator (LASSO) regression and support vector machine-recursive feature elimination (SVM-RFE) algorithms were adopted to select autophagy-related genes. The diagnostic capability of these genes was assessed through analysis with the external validation set from microarray GSE1009, and relevant Chinese medicines were inversely predicted using the SymMap database.
[Results] :
A total of 2 014 DEGs were selected from GSE30528 and GSE30529, leading to the
identification of 37 DKD autophagy-related genes. GO analysis indicated 681 biological
mechanisms, including autophagy regulation and plasma membrane microdomain activity.
KEGG enrichment analysis identified 112 related signaling pathways. PPI network analysis
showed a marked enrichment of autophagy-related genes in DKD. Through LASSO regression and SVM-RFE, four core diagnostic genes for autophagy in DKD were identified: protein
phosphatase 1 regulatory subunit 15A (PPP1R15A), hypoxia inducible factor 1 alpha subunit
(HIF1α), deleted in liver cancer 1 (DLC1), and ceroid lipofuscinosis neuronal 3 (CLN3). The
external validation set demonstrated high diagnostic efficiency for these genes. Finally, 146
kinds of potential Chinese medicines were predicted using the SymMap database, with heatclearing and detoxifying medicine and blood-activating and stasis-eliminating medicine accounting for the largest proportion (25/146 and 13/146, respectively).
[Conclusion]
This study analyzed and validated bioinformatics sequencing databases to elucidate the potential molecular mechanisms of DKD autophagy and predicted key diagnostic genes, potential therapeutic targets, and related Chinese medicines, laying a solid foundation for clinical research and application.
2.Study on the Pathological Mechanism-Syndrome-Treatment Patterns of Approved Chinese Patent Medicines Targeting Collateral Disorders
Pengli SU ; Peng XU ; Yanhong WANG ; Yaqi ZU ; Run YUAN ; Kun LI ; Yufeng ZHAO
Journal of Traditional Chinese Medicine 2025;66(16):1711-1718
ObjectiveTo explore the pathological mechanism-syndrome-treatment patterns of approved Chinese patent medicines (CPMs) that treat collateral disorders, providing a reference for the principle of "treating different diseases with the same therapy" in collateral pathology. MethodsCPMs that apply treatment strategies based on collateral disorders were identified from the Pharmacodia database by extracting information from the "efficacy" or "indications" sections of drug package inserts. A database was established to extract the names and compositions of the CPMs, as well as their indications, related traditional Chinese medicine (TCM) symptoms, disease locations (affected areas), and pathological factors. Frequency statistics were performed. Using the Apriori algorithm, an association rule analysis was conducted on CPMs and disease-location combinations related to the top three most frequent pathological factor combinations. Core formulas for these combinations were identified and analyzed through drug network analysis and MCODE module clustering. ResultsA total of 660 CPMs targeting collateral disorders were retrieved, involving 299 indications, 323 TCM symptoms, 21 disease locations, 19 pathological factors, and 124 pathological factor combinations. The most frequent pathological factor combinations were blood stasis (involved in 109 CPMs, 16.52%), exogenous wind (外风) -cold-dampness (involved in 43 CPMs, 6.52%), and qi deficiency-blood stasis (involved in 42 CPMs, 6.36%). Analysis of the core formulas for these combinations revealed common ingredients such as Honghua (Carthami Flos), Chuanxiong (Chuanxiong Rhizoma), Danggui (Angelicae Sinensis Radix), and Dilong (Pheretima). ConclusionCollateral disorders involve a wide range of pathogenesis and represent a fundamental mechanism in the onset and development of various diseases, characterized by obstruction and stagnation. The primary therapeutic principle is unblocking of the collaterals. Blood stasis obstructing the collaterals is the core pathological basis, and the strategy of activating blood circulation and resolving stasis to unblock the collaterals should be central to the treatment. The core medication pattern involves combining blood-activating and stasis-resolving herbs with insect-derived medicinals that unblock collaterals. Exogenous wind is often the initiating patholo-gical factor in colla-teral disorders, and the appropriate addition of wind-dispelling herbs can enrich the treatment strategies for such conditions.
3.Effect of the impaction of posterior wall on the prognosis following open reduction and internal fixation for fractures of acetabular posterior wall
Yufeng GE ; Feng GAO ; Chao TU ; Gang LIU ; Minghui YANG ; Xu SUN ; Zhelun TAN ; Yimin CHEN ; Weidong PENG ; Xinbao WU
Chinese Journal of Orthopaedic Trauma 2024;26(3):188-193
Objective:To evaluate the effect of the impaction of posterior wall on the prognosis following open reduction and internal fixation for fractures of acetabular posterior wall.Methods:A retrospective study was conducted to analyze the data from the 83 patients with fracture of acetabular posterior wall who had been consecutively treated by open reduction and internal fixation at Department of Orthopaedics and Traumatology, Beijing Jishuitan Hospital from January 2017 to December 2020. The patients were divided into 2 groups based on involvement of posterior wall impaction. In the impaction group of 33 cases, there were 26 males and 7 females with an age of (47.4±11.6) years; in the non-impaction group of 50 cases, there were 43 males and 7 females with an age of (41.3±12.0) years. The quality of postoperative fracture reduction, the function of the affected hip at the last follow-up, and the complication rate during follow-up were compared between the 2 groups. Multifactorial binary logistic regression and age subgroups were used to analyze the effects of posterior wall impaction on functional outcomes.Results:The age, rate of associated injuries in other body parts, and rate of posterior wall comminution in the impaction group were significantly higher than those in the non-impaction group ( P<0.05), but there was no statistically significant difference in other general data of patients between the 2 groups ( P>0.05). All patients were followed up for (44.5±13.3) months after surgery. The rate of anatomical reduction in the non-impaction group (96.0%, 48/50) was significantly higher than that in the impaction group (57.6%, 19/33) ( P<0.05), and the good and excellent rate by the modified Merle d'Aubigné & Postel scale at the last follow-up in the non-impaction group (84.0%, 42/50) was significantly higher than that in the impaction group (51.5%, 17/33) ( P<0.05). There was no significant difference in the incidence of complications between the 2 groups ( P>0.05). After adjusting for age and gender, the difference in hip function was still significantly different between the 2 groups ( OR=0.23, 95% CI: 0.06 to 0.79, P=0.020). The effect of posterior wall impaction on functional outcomes was statistically significant in patients aged ≥50 years ( P=0.008), whereas the difference was not statistically significant in patients aged <50 years ( P=0.194). Conclusions:Compared with non-impaction ones, acetabular fractures of posterior wall impaction tend to lead to poorer quality of reduction, which in turn affects the postoperative recovery of hip joint function. The impact of impaction fractures on functional recovery is more significant in patients aged 50 years and above.
4.Analysis of Medication Regularity of Prescriptions Containing Cangzhu (Atractylodes Lancea) - Huangbai (Phellodendron Amurense) Based on Data Mining
Pengli SU ; Peng XU ; Yanhong WANG ; Yaqi ZU ; Kun LI ; Yufeng ZHAO
Journal of Traditional Chinese Medicine 2024;65(13):1392-1398
ObjectiveTo analyze medication regularity of presciptions containing Cangzhu (Atractylodes Lancea) - Huangbai (Phellodendron Amurense) based on data mining. MethodsAll the prescriptions containing Cangzhu-Huangbai in Dictionary of Chinese Medical Formulas (《中医方剂大辞典》) and the fifth edition of Chinese Pharmacopoeia (《中华医典》) were retrieved, setting up a database, counting the number of prescriptions involved, the composition of formulas in the prescriptions, the properties of the medicinal, the number of types of diseases and syndromes treated by the prescriptions, and the frequency of their appearances; drawing a network diagram of the medicinal combinations; analysing the medicinal combinations based on the association rules of the Apriori algorithm, and at the same time, analysing the association rules for the medicinal involved in the high-frequency diseases, and the core prescriptions containing Cangzhu-Huangbai pairs for specific diseases were extracted, and the medicinal in the core prescriptions were classified into four categories of raising yang and eliminating dampness, fortifying spleen and eliminating dampness, clearing dampness-heat, and specialised disease-use. ResultsA total of 323 prescriptions were identified, involving 318 Chinese herbal medicines, 52 syndromes, 200 diseases, among which arthralgia, beriberi, atrophy-flaccidity, leukorrhea, hernia-type diseases, headache, and haemorrhoids appeared more than 10 times. Through analyzing the core prescriptions for high-frequency diseases, we had four findings. Firstly, medicinal of spleen-fortifying and dampness-removing were the most frequent used in core prescriptions of Cangzhu-Huangbai pair (353 times), and the medication mode was combination of qi-blood supplement medicinal and qi rectifying medicinal, usually using Huangqi (Astragalus mongholicus), Renshen (Panax ginseng), Processed Gancao (Glycyrrhiza glabra), Baizhu (Atractylodes macrocephala), Danggui (Levisticum officinale), Chenpi (Citrus reticulata), Xiangfu (Cyperus rotundus). Then followed by the combination of heat-clearing and dampness-draining medicinal (213 times), usually using Fangji (Stephania tetrandra), Niuxi (Achyranthes bidentata), and Banxia (Pinellia ternata), and this type of heat-clearing and dampness-draining medicinal is most frequently used for haemorrhoids, and beriberi is the most frequent disease. The combination of yang-raising and dampness-removing medicinal (212 times), often using Qianghuo (Hansenia weberbaueriana), Chaihu (Elephantopus scaber), Fangfeng (Carum carvi), and this type of medicinal is most frequently used for headache, and atrophy-flaccidity is the most frequent disease. The combination of medicinal for specialised diseases was used the least frequently (39 times). ConclusionCangzhu-Huangbai pair were often combined with medicinal of yang-raising and dampness-removing, spleen-fortifying and dampness-removing, and heat-clearing and dampness-draining to expel dampness pathogen, at the same time combined with disease-specific medicinal to treat disease and syndrome together.
5.Tim4 deficiency reduces CD301b+macrophage and aggravates periodontitis bone loss
Wang ZIMING ; Zeng HAO ; Wang CAN ; Wang JIAOLONG ; Zhang JING ; Qu SHUYUAN ; Han YUE ; Yang LIU ; Ni YUEQI ; Peng WENAN ; Liu HUAN ; Tang HUA ; Zhao QIN ; Zhang YUFENG
International Journal of Oral Science 2024;16(2):280-292
Periodontitis is a common chronic inflammatory disease that causes the periodontal bone destruction and may ultimately result in tooth loss.With the progression of periodontitis,the osteoimmunology microenvironment in periodontitis is damaged and leads to the formation of pathological alveolar bone resorption.CD301b+macrophages are specific to the osteoimmunology microenvironment,and are emerging as vital booster for conducting bone regeneration.However,the key upstream targets of CD301b+macrophages and their potential mechanism in periodontitis remain elusive.In this study,we concentrated on the role of Tim4,a latent upstream regulator of CD301b+macrophages.We first demonstrated that the transcription level of Timd4(gene name of Tim4)in CD301b+macrophages was significantly upregulated compared to CD301b-macrophages via high-throughput RNA sequencing.Moreover,several Tim4-related functions such as apoptotic cell clearance,phagocytosis and engulfment were positively regulated by CD301b+macrophages.The single-cell RNA sequencing analysis subsequently discovered that Cd301b and Timd4 were specifically co-expressed in macrophages.The following flow cytometric analysis indicated that Tim4 positive expression rates in total macrophages shared highly synchronized dynamic changes with the proportions of CD301b+macrophages as periodontitis progressed.Furthermore,the deficiency of Tim4 in mice decreased CD301b+macrophages and eventually magnified alveolar bone resorption in periodontitis.Additionally,Tim4 controlled the p38 MAPK signaling pathway to ultimately mediate CD301b+macrophages phenotype.In a word,Tim4 might regulate CD301b+macrophages through p38 MAPK signaling pathway in periodontitis,which provided new insights into periodontitis immunoregulation as well as help to develop innovative therapeutic targets and treatment strategies for periodontitis.
6.The osteoclastic activity in apical distal region of molar mesial roots affects orthodontic tooth movement and root resorption in rats
Zheng WENHAO ; Lu XIAOFENG ; Chen GUANGJIN ; Shen YUFENG ; Huang XIAOFEI ; Peng JINFENG ; Wang JIAJIA ; Yin YING ; Song WENCHENG ; Xie MENGRU ; Yu SHAOLING ; Chen LILI
International Journal of Oral Science 2024;16(2):322-332
The utilization of optimal orthodontic force is crucial to prevent undesirable side effects and ensure efficient tooth movement during orthodontic treatment.However,the sensitivity of existing detection techniques is not sufficient,and the criteria for evaluating optimal force have not been yet established.Here,by employing 3D finite element analysis methodology,we found that the apical distal region(A-D region)of mesial roots is particularly sensitive to orthodontic force in rats.Tartrate-resistant acidic phosphatase(TRAP)-positive osteoclasts began accumulating in the A-D region under the force of 40 grams(g),leading to alveolar bone resorption and tooth movement.When the force reached 80 g,TRAP-positive osteoclasts started appearing on the root surface in the A-D region.Additionally,micro-computed tomography revealed a significant root resorption at 80 g.Notably,the A-D region was identified as a major contributor to whole root resorption.It was determined that 40 g is the minimum effective force for tooth movement with minimal side effects according to the analysis of tooth movement,inclination,and hyalinization.These findings suggest that the A-D region with its changes on the root surface is an important consideration and sensitive indicator when evaluating orthodontic forces for a rat model.Collectively,our investigations into this region would aid in offering valuable implications for preventing and minimizing root resorption during patients'orthodontic treatment.
7.Expression and clinical significance of long non-coding RNA differentially expressed genes, micro RNA-181a-5p and autophagy-related protein 5 in primary gouty arthritis
Jianwei GUO ; Tianyi LEI ; Peng WANG ; Zeng ZHANG ; Guilin JIAN ; Quanbo ZHANG ; Yufeng QING
Chinese Journal of Rheumatology 2024;28(5):303-311
Objective:To explore the expression and clinical significance of long non-coding RNA colorectal neoplasia differentially expressed (lncRNA CRNDE), microRNA-181a-5p (miR-181a-5p) and autophagy related 5 (ATG5) in the peripheral blood of patients with gouty arthritis (GA) patients.Methods:The clinical data, laboratory parameters and peripheral blood samples were collected from 40 patients with acute gout (AG), 40 patients with intermittent gout (IG) and 50 healthy subjects (HC). The expression levels of lncRNA CRNDE, miR-181a-5p and ATG5 mRNA were detected by real-time fluorescence quantification (RT-qPCR) and the expression level of ATG5 protein was detected by Western-blot. The expression levels of lncRNA CRNDE, miR-181a-5p, ATG5 mRNA were compared among the three groups and correlated with clinical indices, and a subject operating characteristic curve (ROC) was constructed to assess the value of lncRNA CRNDE, miR-181a-5p, ATG5 mRNA in the diagnosis of gout. Measurements conforming to normal distribution were analyzed using t test or ANOVA, data with non-normal distribution was analyzed using Mann-Whitney U test or Kruskal-Wallis H test, correlation analysis between variables was analyzed using Spearman's analysis, and the diagnostic value of each indicator was analyzed using ROC curve. Results:① The differences in the expression of lncRNA CRNDE, miR-181a-5p, and ATG5 mRNA between the three groups were statistically significant ( H=32.12, 57.73, 68.32, all P<0.001). Among them, lncRNA CRNDE expression level in the AG group was significantly higher than that in the IG group and healthy control group [61.95(11.39, 108.30)×10 -3, 25.71(15.40, 38.40)×10 -3, 13.80(3.97, 23.99)×10 -3; Z=-3.24, P=0.001; Z=-5.03, P<0.001], and the expression level of IG group was higher than that of healthy control group( Z=-3.56, P<0.001); miR-181a-5p and ATG5 mRNA expression levels in AG group were significantly lower than those in IG group and healthy control group [miR-181a-5p: 39.81(31.22, 69.38)×10 -3, 60.74(44.19, 90.35)×10 -3, 121.30(101.50, 316.90)×10 -3; Z=-3.01, P=0.030; Z=-6.93, P<0.001. ATG5 mRNA: 4.52(2.31, 26.63)×10 -3, 43.63(13.72, 102.70)×10 -3, 153.90(66.62, 365.80)×10 -3; Z=-5.47, -7.36, all P<0.001)], which were expressed at lower levels in the IG group than in the healthy controls ( Z=-5.25, -4.47, all P<0.001). The difference of ATG5 protein expression level among the three groups expressed was statistically significant ( F=6.24, P=0.030), and the AG group was higher than the healthy control group, and the difference was statistically significant [(0.96±0.13) vs.(0.61±0.04), t=4.25, P=0.013], but the difference between the IG group (0.78±0.15) and the AG group and the HC group was not statistically significant ( t=1.51, P=0.206; t=1.85, P=0138). ② Spearman correlation analysis showed that lncRNA CRNDE was negatively correlated with the expression levels of miR-181a-5p and ATG5 mRNA in gout patients ( r=-0.49, P<0.001; r=-0.35, P=0.002); miR-181a-5p was positively correlated with ATG5 mRNA expression levels ( r=0.64, P<0.001); lncRNA CRNDE expression level was positively correlated with ESR and WBC ( r=0.49, P<0.001; r=0.43, P=0.001); miR-181a-5p expression level was negatively correlated with ESR and WBC ( r=-0.29, P=0.009; r=-0.35, P=0.002), and ATG5 mRNA expression levels were negatively correlated with ESR, WBC, and GR ( r=-0.26, P=0.021; r=-0.26, P=0.024; r=-0.27, P=0.021). In the AG group lncRNA CRNDE was positively correlated with ESR and WBC ( r=0.36, P=0.022; r=0.36, P=0.026) and miR-181a-5p was negatively correlated with WBC ( r=-0.34, P=0.038) ③ ROC curve showed that the areas under ROC curve of lncRNA CRNDE, miR-181a-5p and ATG5 mRNA expression levels to predict gout were 0.764, 0.875 and 0.864, respectively. The area under ROC curve of gout predicted by the three combined was 0.928. Conclusion:lncRNA CRNDE, miR-181a-5p, and ATG5 may be involved in the pathoge-nesis of primary gouty arthritis, and are potential biological parameters for studying the pathogenesis of gout.
8.Prevalence of parasitic infections in human stool samples from a hospital in Chenzhou City of Hunan Province
Yufeng PENG ; Xiaomei LIAO ; Lepan ZHU ; Yangnan ZHANG
Chinese Journal of Schistosomiasis Control 2023;35(3):291-293
Objective To investigate the prevalence of parasitic infections in human stool samples from a hospital in Chenzhou City, Hunan Province, so as to provide insights into the management of intestinal parasitic diseases. Methods Stool samples were collected from patients admitted to a hospital in Chenzhou City from September 2020 to March 2021, subjected to physiological saline smearing and microscopy for detection of intestinal parasites. The prevalence of parasitic infections and the species of parasites were descriptively analyzed. Results The overall prevalence of intestinal parasitic infections was 1.61% in the 10 728 stool samples, and there were 3 samples with mixed infections of two parasite species. A total of seven parasite species were identified, including Blastocystis hominis (162 cases, 1.55%), Giardia lamblia (5 cases, 0.05%), Dientamoeba fragilis (5 cases, 0.05%), Endolimax nana (one case, 0.01%), Iodamoeba bütschlii (one case, 0.01%), Strongyloides stercoralis (one case, 0.01%) and Trichomonas hominis (one case, 0.01%). The prevalence of intestinal parasitic infection was significantly higher among women than in men (2.14% vs. 1.25%; χ2 = 13.01, P < 0.01), and a high prevalence rate was seen among patients at ages of 20 to 30 years (2.99%) and 80 years and older (2.86%); however, no age-specific prevalence of intestinal parasitic infection was detected (χ2 = 12.45, P > 0.05). Conclusions The overall prevalence of intestinal parasitic infection was low among patients admitted to a hospital in Chenzhou City, and gender-specific prevalence was found. Food-borne and opportunistic parasites were predominant intestinal parasites, including B. hominis, G. lamblia and D. fragilis.
9.Value of number of negative lymph nodes in predicting the prognosis of patients with esophageal cancer after neoadjuvant therapy and the construction of nomogram prodiction model
Yueyang YANG ; Peng TANG ; Zhentao YU ; Haitong WANG ; Hongdian ZHANG ; Mingquan MA ; Yufeng QIAO ; Peng REN ; Xiangming LIU ; Lei GONG
Chinese Journal of Digestive Surgery 2023;22(3):371-382
Objective:To investigate the value of number of negative lymph nodes (NLNs) in predicting the prognosis of patients with esophageal cancer after neoadjuvant therapy and the construction of nomogram prodiction model.Methods:The retrospective cohort study was conducted. The clinicopathological data of 1 924 patients with esophageal cancer after neoadjuvant therapy uploaded to the Surveillance, Epidemiology, and End Results Database of the National Cancer Institute from 2004 to 2015 were collected. There were 1 624 males and 300 females, aged 63 (range, 23?85)years. All 1 924 patients were randomly divided into the training dataset of 1 348 cases and the validation dataset of 576 cases with a ratio of 7:3 based on random number method in the R software (3.6.2 version). The training dataset was used to constructed the nomogram predic-tion model, and the validation dataset was used to validate the performance of the nomogrram prediction model. The optimal cutoff values of number of NLNs and number of examined lymph nodes (ELNs) were 8, 14 and 10, 14, respectively, determined by the X-tile software (3.6.1 version), and then data of NLNs and ELNs were converted into classification variables. Observation indicators: (1) clinicopathological characteristics of patients in the training dataset and the validation dataset; (2) survival of patients in the training dataset and the validation dataset; (3) prognostic factors analysis of patients in the training dataset; (4) survival of patients in subgroup of the training dataset; (5) prognostic factors analysis in subgroup of the training dataset; (6) construction of nomogram prediction model and calibration curve. Measurement data with normal distribution were represented as Mean± SD, and comparison between groups was conducted using the t test. Measurement data with skewed distribution were represented as M(range), and comparison between groups was conducted using the Mann-Whitney U test. Count data were described as absolute numbers, and comparison between groups was conducted using the chi-square test. The Kaplan-Meier method was used to draw survival curve and Log-Rank test was used for survival analysis. The COX proportional hazard model was used for univariate and multivariate analyses. Based on the results of multivariate analysis, the nomogram prediction model was constructed. The prediction efficacy of nomogram prediction model was evaluated using the area under curve (AUC) of the receiver operating characteristic curve and the Harrell′s c index. Errors of the nomogram prediction model in predicting survival of patients for the training dataset and the validation dataset were evaluated using the calibration curve. Results:(1) Clinicopathological characteristics of patients in the training dataset and the validation dataset. There was no significant difference in clinicopatholo-gical characteristics between the 1 348 patients of the training dataset and the 576 patients of the validation dataset ( P>0.05). (2) Survival of patients in the training dataset and the validation dataset. All 1 924 patients were followed up for 50(range, 3?140)months, with 3-year and 5-year cumulative survival rate as 59.4% and 49.5%, respectively. The 3-year cumulative survival rate of patients with number of NLNs as <8, 8?14 and >14 in the training dataset was 46.7%, 62.0% and 66.0%, respectively, and the 5-year cumulative survival rate was 38.1%, 52.1% and 59.7%, respectively. There was a significant difference in the survival of these patients in the training dataset ( χ2=33.70, P<0.05). The 3-year cumulative survival rate of patients with number of NLNs as <8, 8?14 and >14 in the validation dataset was 51.1%, 54.9% and 71.2%, respectively, and the 5-year cumulative survival rate was 39.3%, 42.5% and 55.7%, respectively. There was a significant difference in the survival of these patients in the validation dataset ( χ2=14.49, P<0.05). The 3-year cumulative survival rate of patients with number of ELNs as <10, 10?14 and >14 in the training dataset was 53.9%, 60.0% and 62.7%, respectively, and the 5-year cumulative survival rate was 44.7%, 49.1% and 56.9%, respectively. There was a significant difference in the survival of these patients in the training dataset ( χ2=9.88, P<0.05). The 3-year cumulative survival rate of patients with number of ELNs as <10, 10?14 and >14 in the validation dataset was 56.2%, 47.9% and 69.3%, respectively, and the 5-year cumula-tive survival rate was 44.9%, 38.4% and 51.9%, respectively. There was a significant difference in the survival of these patients in the validation dataset ( χ2=9.30, P<0.05). (3) Prognostic factors analysis of patients in the training dataset. Results of multivariate analysis showed that gender, neoadjuvant pathological (yp) T staging, ypN staging (stage N1, stage N2, stage N3) and number of NLNs (8?14, >14) were independent influencing factors for the prognosis of patients with esophageal cancer after neoadjuvant therapy ( hazard ratio=0.65, 1.44, 1.96, 2.41, 4.12, 0.69, 0.56, 95% confidence interval as 0.49?0.87, 1.17?1.78, 1.59?2.42, 1.84?3.14, 2.89?5.88, 0.56?0.86, 0.45?0.70, P<0.05). (4) Survival of patients in subgroup of the training dataset. Of the patients with NLNs in the training dataset, the 3-year cumulative survival rate of patients with number of NLNs as <8, 8?14 and >14 was 61.1%, 71.6% and 76.8%, respectively, and the 5-year cumulative survival rate was 50.7%, 59.9% and 70.1%, respectively. There was a significant difference in the survival of these patients in the training dataset ( χ2=12.66, P<0.05). Of the patients with positive lymph nodes in the training dataset, the 3-year cumulative survival rate of patients with number of NLNs as <8, 8?14 and >14 was 26.1%, 42.9% and 44.7%, respectively, and the 5-year cumulative survival rate was 20.0%, 36.5% and 39.3%, respectively. There was a significant difference in the survival of these patients in the training dataset ( χ2=20.39, P<0.05). (5) Prognostic factors analysis in subgroup of the training dataset. Results of multivariate analysis in patients with NLNs in the training dataset showed that gender, ypT staging and number of NLNs (>14) were independent influencing factors for the prognosis of patients with esophageal cancer after neoadju-vant therapy ( hazard ratio=0.67, 1.44, 0.56, 95% confidence interval as 0.47?0.96, 1.09?1.90, 0.41?0.77, P<0.05). Results of multi-variate analysis in patients with positive lymph nodes in the training dataset showed that race as others, histological grade as G2, ypN staging as stage N3 and number of NLNs (8?14, >14) were independent influencing factors for the prognosis of patients with esophageal cancer after neoadjuvant therapy ( hazard ratio=2.73, 0.70, 2.08, 0.63, 0.59, 95% confidence interval as 1.43?5.21, 0.54?0.91, 1.44?3.02, 0.46?0.87, 0.44?0.78, P<0.05). (6) Construction of nomogram prediction model and calibration curve. Based on the multivariate analysis of prognosis in patients of the training dataset ,the nomogram prediction model for the prognosis of patients with esophageal cancer after neoadju-vant treatment was constructed based on the indicators of gender, ypT staging, ypN staging and number of NLNs. The AUC of nomogram prediction model in predicting the 3-, 5-year cumulative survival rate of patients in the training dataset and the validation dataset was 0.70, 0. 70 and 0.71, 0.71, respectively. The Harrell′s c index of nomogram prediction model of patients in the training dataset and the validation dataset was 0.66 and 0.63, respectively. Results of calibration curve showed that the predicted value of the nomogram prediction model of patients in the training dataset and the validation dataset was in good agreement with the actual observed value. Conclusion:The number of NLNs is an independent influencing factor for the prognosis of esophageal cancer patients after neoadjuvant therapy, and the nomogram prediction model based on number of NLNs can predict the prognosis of esophageal cancer patients after neoadjuvant therapy.
10.Differential expression profiles analysis of DNA methylation between “disease” and “syndrome” in coronary heart disease-induced unstable angina patientswith Qi deficiency and blood stasis syndrome
WU Huaying ; HU Hongchun ; LIU Yufeng ; LI Liang ; LI Jing ; HAN Yuming ; XIAO Changjiang ; PENG Qinghua
Digital Chinese Medicine 2023;6(4):451-466
Objective:
To explore the differential expression profiles of DNA methylation sites/regions and potential molecular mechanisms in the peripheral blood of coronary heart disease (CHD)-induced unstable angina pectoris patients with or without Qi deficiency and blood stasis syndrome, and to provide scientific evidence for the conbination of disease and syndrome.
Methods:
According to the pre-determined inclusion and exclusion criteria, the study subjects were enrolled and divided into two groups namely CHD-induced unstable angina group (G group) and healthy control group (J group) to conduct “disease” analysis, while G group was further divided into Qi deficiency and blood stasis syndrome group (case group) and non-Qi deficiency blood stasis syndrome group (control group) to perform “syndrome” analysis. The general data and clinical information of the study subjects were collected. The peripheral venous blood was extracted on an empty stomach, and the Illumina Infinium MethylationEPIC BeadChip (850K methylation chip) was used to detect the differential expressionprofiles of DNA methylation in each group, ChAMP software (V 2.14.0) was used for the differential methylation data analysis, with a threshold of the adjusted P value (adj.P.val) < 0.01. Gene Ontology (GO) and Kyoto Encyclopedia of Genomes (KEGG) were employed for the functional and pathway enrichment analyses of related mapped genes.
Results:
A total of 263 differentially methylated CpG positions (DMPs) were screened out between G and J groups, including 191 hypermethylated positions such as cg05845204 and cg08906898, and 72 hypomethylated positions such as cg26919182 and cg13149459. These positions were mainly mapped to 148 genes encompassing RNA binding motif protein 39 (RBM39), acetyl-CoA acyltransferase 2 (ACAA2), protein phosphatase 1 regulatory subunit 12B (PPP1R12B), and the dual-specificity tyrosine phosphorylation-regulated kinase 2 (DYRK2). GO functional enrichment analysis revealed that the genes of the DMPs were primarily enriched in protein localization to chromosomes, regulation of cell morphogenesis, negative regulation of calcium-mediated signals, etc. KEGG pathway analysis suggested that the genes were mainly enriched in fatty acid metabolism and endocytosis pathways. In addition, a total of 23 differential methylation regions (DMRs) were identified, with overlapping genes such as transmembrane protein 232 (TMEM232), ribosomal protein large P1 (RPLP1), peroxisomal biogenesis factor 10 (PEX10), and forkhead box N3 (FOXN3) recognized. It was found that GO functions were mainly enriched in the negative regulation of Ras protein signal transduction, small GTPase-mediated signal transduction, negative regulation, etc. A total of 1 703 differential methylation sites were screened out between case and control groups, including 444 increased methylation positions such as cg05573767 and 1 259 decreased methylationpositions such as cg19938535, and cg03893872. These positions were mapped to 1 108 genes such as ribosomal protein S6 kinase A2 (RPS6KA2), leucine rich repeat containing 16A (LRRC16A), and hedgehog acyltransferase (HHAT). According to the GO functional enrichment analysis, the genes relating to the DMPs were mainly enriched in biological functions such as transmembrane receptor protein serine/threonine kinase signaling pathway and axonogenesis. The KEGG pathway enrichment analysis suggested the involvement of Rap1 signaling pathway, adenosine 5’-monophosphate-activated protein kinase (AMPK) signaling pathway, etc. A total of 21 DMRs were identified, including 22 overlapping genes such as mucin 4 (MUC4), three prime repair exonuclease 1 (TREX1), and LIM homeobox 6 (LHX6). GO analysis demonstrated that the genes primarily participated in molecular functions such as positive regulation of transmembrane transport, regulation of fatty acid metabolism, and copper ion binding.
Conclusion
This study reveals the methylation patterns of DMPs and DMRs in patients with Qi deficiency and blood stasis syndrome caused by CHD-induced unstable angina pectoris. Potential epigenetic regulation of fatty acid metabolism, Rap1 signaling, and other molecular functions are involved in the development of CHD between the "disease" and "syndrome".

Result Analysis
Print
Save
E-mail