1.Pre-operative risk assessment of hepatocellular carcinoma recurrence in liver transplant recipients by non-invasive detection of pre-existing genetic lesions
Suqin YANG ; Sunbin LING ; Jianhua LI ; Yan WANG ; Jiapei WANG ; Qiwei HUANG ; Fanming LIU ; Yiqi ZHUANG ; Yingyu ZHENG ; Rui WANG ; Zhe YANG ; Xiaoping ZHENG ; Kai WANG ; Zhikun LIU ; Jun CHEN ; Jianguo WANG ; Haiyang XIE ; Lin ZHOU ; Leiming CHEN ; Guoqiang CAO ; Dandan CHEN ; Junfang JI ; Bin ZHAO ; Chao JIANG ; Di LU ; Xuyong WEI ; Hangjin JIANG ; Qiaonan SHAN ; Hengbo SHI ; Yong-Zhen XU ; Shusen ZHENG ; Zhengxin WANG ; Shengda LIN ; Xiao XU
Clinical and Molecular Hepatology 2026;32(2):884-903
Background/Aims:
Liver transplantation (LT) following total hepatectomy is a life-saving treatment for hepatocellular carcinoma (HCC). The HCC recurrence after LT hinders the effectiveness of the procedure. The objective of this study is to develop a pre-operative risk stratification model based on a liquid biopsy.
Methods:
We conducted a comprehensive multi-omics study of 260 HCC patients from three centers, including clinical data, low-coverage whole-genome sequencing of cell-free DNA (cfDNA) from plasma, as well as whole-exome, single-nucleus RNA, and spatial transcriptomics from matched tumor and non-tumor tissues.
Results:
We identified cfDNA-derived copy number alteration (CNA) signatures associated with post-transplant recurrence. By integrating cfDNA-derived CNA profiles with single-cell transcriptomic data, we traced recurrence-associated cfDNA to a distinct subpopulation of malignant cells within the primary tumor. These cells were embedded in a pro-metastatic microenvironment of specialized endothelial subtypes and cancer-associated fibroblasts. Notably, most recurrence-associated lesions were detectable in cfDNA prior to liver transplantation (LT). Building on these insights, we developed the ZJU Criteria based on CNA fragments and tumor markers, a pre-LT risk prediction tool that integrates conventional clinical factors with cfDNA-derived CNA signatures, and validated it using internal and independent external cohorts.
Conclusion
Our findings suggest that post-transplant recurrence commonly originates from advanced subclones that emerge late during tumor evolution. The ZJU Criteria provides an accurate, non-invasive strategy that significantly improves pre-LT risk stratification and clinical decision-making for patients with HCC.
2.Erratum: Author correction to "Generation of αGal-enhanced bifunctional tumor vaccine" Acta Pharm Sin B 12 (2022) 3177-3186.
Jian HE ; Yu HUO ; Zhikun ZHANG ; Yiqun LUO ; Xiuli LIU ; Qiaoying CHEN ; Pan WU ; Wei SHI ; Tao WU ; Chao TANG ; Huixue WANG ; Lan LI ; Xiyu LIU ; Yong HUANG ; Yongxiang ZHAO ; Lu GAN ; Bing WANG ; Liping ZHONG
Acta Pharmaceutica Sinica B 2025;15(2):1207-1207
[This corrects the article DOI: 10.1016/j.apsb.2022.03.002.].
3.Evaluation value of C-reactive protein/albumin ratio combined with platelet count and Glasgow coma scale for prognosis of patients with heat stroke.
Shanshan SHI ; Zhengzhen WU ; Yong HUANG ; Xianglei FU
Chinese Critical Care Medicine 2025;37(2):160-164
OBJECTIVE:
To explore the prognostic value of C-reactive protein (CRP)/albumin (Alb) ratio combined with platelet count (PLT) and Glasgow coma score (GCS) in patients with heat stroke (HS).
METHODS:
A retrospective analysis was conducted on the clinical data of HS patients admitted to the department of intensive care unit (ICU) of Nanchong Central Hospital from May 1, 2020 to October 31, 2023. This included general information, admission GCS, laboratory indicators and 28-day prognosis. The differences in the above indicators were compared between two groups of patients with different prognoses. Statistically significant indicators from univariate analysis were included in multivariate Logistic regression analysis to screen for factors influencing 28-day mortality in HS patients. The predictive value of various influencing factors on the 28 days prognosis of HS patients were analyzed by receiver operator characteristic curve (ROC curve).
RESULTS:
A total of 73 HS patients were included, of whom 41 survived for 28-day and 32 died. There were no statistically significant differences in gender and age between the two groups of HS patients with different prognoses. The white blood cell count (WBC), neutrophil count (NEU), aspartate aminotransferase (AST), alanine aminotransferase (ALT), CRP, and CRP/Alb ratio in the death group were significantly higher than those of the survival group, and the admission GCS score, platelet count (PLT), total bilirubin (TBil) and Alb were significantly lower than the survival group [WBC (×109/L): 14.80 (11.44, 17.15) vs. 11.96 (9.47, 14.82), NEU (×109/L): 13.05 (8.56, 15.67) vs. 9.50 (6.68, 12.09), AST (U/L): 108.00 (52.70, 291.50) vs. 64.50 (38.25, 110.50), ALT (U/L): 62.00 (19.50, 159.00) vs. 34.50 (20.75, 70.75), CRP (mg/L): 22.49 (3.42, 58.93) vs. 3.68 (1.01, 11.46), CRP/Alb ratio: 0.53 (0.08, 1.77) vs. 0.08 (0.02, 0.44), GCS score: 7.0 (5.0, 8.0) vs. 8.5 (7.0, 11.0), PLT (×109/L): 107.00 (73.50, 126.00) vs. 131.50 (107.50, 176.25), TBil (mmol/L): 15.60 (10.00, 25.30) vs. 21.40 (14.80, 30.05), Alb (g/L): 32.65 (32.53, 49.30) vs. 38.70 (36.20, 40.40), all P < 0.05]. Binary Logistic regression analysis showed that the GCS score [odds ratio (OR) = 0.686, 95% confidence interval (95%CI) was 0.491-0.959, P = 0.028], PLT (OR = 0.973, 95%CI was 0.954-0.992, P = 0.005), NEU (OR = 1.312, 95%CI was 1.072-1.606, P = 0.009) and CRP/Alb ratio (OR = 7.652, 95%CI was 1.632-35.881, P = 0.010) were independent influencing factors for 28-day mortality in HS patients. ROC curve analysis showed that the area under the curve (AUC) of GCS score, PLT, and CRP/Alb ratio for single prediction of 28-day prognosis in HS patients was 0.705, 0.752, and 0.729, and the combination of all three predicted the highest AUC of 28-day prognosis in HS patients (0.917), with a sensitivity and specificity of 86.2% and 81.2%, respectively.
CONCLUSION
CRP/Alb ratio, PLT, and GCS score are independent influencing factors affecting the prognosis of HS patients, and all of them have a certain predictive value for the prognosis of HS patients, in which the combination of the three has a higher predictive value for the prognosis of HS patients.
Humans
;
C-Reactive Protein/analysis*
;
Prognosis
;
Glasgow Coma Scale
;
Retrospective Studies
;
Heat Stroke/diagnosis*
;
Platelet Count
;
Male
;
Female
;
Serum Albumin/analysis*
;
Middle Aged
;
Aged
;
Adult
;
ROC Curve
4.The Impacts of Climate Change on the Environment and Human Health in China: A Call for more Ambitious Action.
Shi Lu TONG ; Yu WANG ; Yong Long LU ; Cun de XIAO ; Qi Yong LIU ; Qi ZHAO ; Cun Rui HUANG ; Jia Yu XU ; Ning KANG ; Tong ZHU ; Dahe QIN ; Ying XU ; Buda SU ; Xiao Ming SHI
Biomedical and Environmental Sciences 2025;38(2):127-143
As global greenhouse gases continue rising, the urgency of more ambitious action is clearer than ever before. China is the world's biggest emitter of greenhouse gases and one of the countries affected most by climate change. The evidence about the impacts of climate change on the environment and human health may encourage China to take more decisive action to mitigate greenhouse gas emissions and adapt to climate impacts. This article aimed to review the evidence of environmental damages and health risks posed by climate change and to provide a new science-based perspective for the delivery of sustainable development goals. Over recent decades, China has experienced a strong warming pattern with a growing frequency of extreme weather events, and the impacts of climate change on China's environment and human health have been consistently observed, with increasing O 3 air pollution, decreases in water resources and availability, land degradation, and increased risks for both communicable and non-communicable diseases. Therefore, China's climate policy should target the key factors driving climate change and scale up strategic measures to curb carbon emissions and adapt to inevitable increasing climate impacts. It provides new insights for not only China but also other countries, particularly developing and emerging economies, to ensure climate and environmental sustainability whilst pursuing economic growth.
Climate Change
;
China
;
Humans
;
Greenhouse Gases
;
Air Pollution
;
Sustainable Development
;
Environment
5.The multi-center mid-term clinical outcomes of combined complete preservation of chordal structure mitral valve replacement with total anatomic complete arterial myocardial revascularization for coronary patients with moderate-to-severe or severe ischemic mitral regurgitation
Ke GUO ; Xujun CHEN ; Baoshi ZHENG ; Chao SHI ; Keli HUANG ; Yong CAO ; Chengquan LIAO ; Jingwei CHEN ; Yu LIN ; Chengxin LIU ; Quansheng CAO ; Lin SHEN ; Zhendong WANG
Chinese Journal of Surgery 2025;63(1):58-67
Objective:To evaluate the clinical outcomes of combined complete preservation of chordal structure mitral valve replacement (C-MVR) with total anatomical arterial myocardial revascularization (TACR) in coronary patients with moderate-to-severe or severe ischemic mitral regurgitation (IMR).Methods:This is a retrospective multi-center case series study. Data were retrospectively collected from 127 patients with coronary artery disease with moderate to severe or severe IMR who received TACR with C-MVR from July 2015 to April 2024 in 13 hospitals in China. There were 90 males and 37 females, aged (56.5±10.7) years (range: 33 to 74 years). Perioperative data and follow-up data including left ventricular ejection fraction, left ventricular end-diastolic diameter, and patency rate of arterial grafts of patients were collected. Comparisons were made using paired sample t-test or χ2 test. Results:In this cohort of 127 patients, 67 underwent concurrent tricuspid valve repair. During surgery, 113 grafts of the left internal mammary artery (LIMA), 127 grafts of the left radial artery, 80 grafts of the right radial artery, and 110 grafts of the right internal mammary artery (RIMA) were harvested. The number of the distal anastomosis was 4.2±0.4 (range: 3 to 5). The aortic cross-clamp time and cardiopulmonary bypass time were (97.5±23.4) minutes (range: 90 to 161 minutes) and (145.4±19.2) minutes (range: 101 to 210 minutes), respectively. There was one operative death. Intraoperative placement of an intra-aortic balloon pump was performed in 21 patients to improve the left ventricular ejection. No sternal ischemic occurred. All patients completed follow-up, with a mean follow-up period of (64.3±7.5) months (range: 4 to 110 months). No major cerebrovascular events occurred during the follow-up period, and all patients survived. Left ventricular ejection fraction improved postoperatively (55.0%±5.3% vs. 41.0%±15.3%, t=17.23, P<0.01). The proportion of patients with New York Heart Association functional class ≤2 increased postoperatively (23.6% (30/127) vs. 87.3% (110/126), χ2=103.77, P<0.01). The proportion of patients with Canadian Cardiovascular Society Angina Classification ≤3 decreased postoperatively (4.8% (6/126) vs. 78.7% (100/127), χ2=142.19, P<0.01). The left ventricular end-diastolic diameter decreased postoperatively ((5.70±4.50) cm vs. (6.10±0.23) cm, t=12.15, P<0.01). Coronary multi-detector computed tomography angiography (MDCTA) follow-up was conducted for (60.5±11.7) months (range: 6 to 109 months) postoperatively. MDCTA confirmed the patency rates of the grafts: 96.4% (108/112) for the LIMA grafts, 88.9% (112/126) for the left radial artery grafts, 93.7% (74/79) for the right radial artery grafts, and 90.9% (100/110) for the free RIMA grafts. No significant differences in graft patency rates were observed between the arterial grafts ( χ2=5.24, P=0.155). Conclusion:The results of this multi-centre study demonstrate satisfactory mid-term results of C-MVR with TACR for the treatment of coronary artery disease with moderate to severe or severe IMR.
6.The multi-center mid-term clinical outcomes of combined complete preservation of chordal structure mitral valve replacement with total anatomic complete arterial myocardial revascularization for coronary patients with moderate-to-severe or severe ischemic mitral regurgitation
Ke GUO ; Xujun CHEN ; Baoshi ZHENG ; Chao SHI ; Keli HUANG ; Yong CAO ; Chengquan LIAO ; Jingwei CHEN ; Yu LIN ; Chengxin LIU ; Quansheng CAO ; Lin SHEN ; Zhendong WANG
Chinese Journal of Surgery 2025;63(1):58-67
Objective:To evaluate the clinical outcomes of combined complete preservation of chordal structure mitral valve replacement (C-MVR) with total anatomical arterial myocardial revascularization (TACR) in coronary patients with moderate-to-severe or severe ischemic mitral regurgitation (IMR).Methods:This is a retrospective multi-center case series study. Data were retrospectively collected from 127 patients with coronary artery disease with moderate to severe or severe IMR who received TACR with C-MVR from July 2015 to April 2024 in 13 hospitals in China. There were 90 males and 37 females, aged (56.5±10.7) years (range: 33 to 74 years). Perioperative data and follow-up data including left ventricular ejection fraction, left ventricular end-diastolic diameter, and patency rate of arterial grafts of patients were collected. Comparisons were made using paired sample t-test or χ2 test. Results:In this cohort of 127 patients, 67 underwent concurrent tricuspid valve repair. During surgery, 113 grafts of the left internal mammary artery (LIMA), 127 grafts of the left radial artery, 80 grafts of the right radial artery, and 110 grafts of the right internal mammary artery (RIMA) were harvested. The number of the distal anastomosis was 4.2±0.4 (range: 3 to 5). The aortic cross-clamp time and cardiopulmonary bypass time were (97.5±23.4) minutes (range: 90 to 161 minutes) and (145.4±19.2) minutes (range: 101 to 210 minutes), respectively. There was one operative death. Intraoperative placement of an intra-aortic balloon pump was performed in 21 patients to improve the left ventricular ejection. No sternal ischemic occurred. All patients completed follow-up, with a mean follow-up period of (64.3±7.5) months (range: 4 to 110 months). No major cerebrovascular events occurred during the follow-up period, and all patients survived. Left ventricular ejection fraction improved postoperatively (55.0%±5.3% vs. 41.0%±15.3%, t=17.23, P<0.01). The proportion of patients with New York Heart Association functional class ≤2 increased postoperatively (23.6% (30/127) vs. 87.3% (110/126), χ2=103.77, P<0.01). The proportion of patients with Canadian Cardiovascular Society Angina Classification ≤3 decreased postoperatively (4.8% (6/126) vs. 78.7% (100/127), χ2=142.19, P<0.01). The left ventricular end-diastolic diameter decreased postoperatively ((5.70±4.50) cm vs. (6.10±0.23) cm, t=12.15, P<0.01). Coronary multi-detector computed tomography angiography (MDCTA) follow-up was conducted for (60.5±11.7) months (range: 6 to 109 months) postoperatively. MDCTA confirmed the patency rates of the grafts: 96.4% (108/112) for the LIMA grafts, 88.9% (112/126) for the left radial artery grafts, 93.7% (74/79) for the right radial artery grafts, and 90.9% (100/110) for the free RIMA grafts. No significant differences in graft patency rates were observed between the arterial grafts ( χ2=5.24, P=0.155). Conclusion:The results of this multi-centre study demonstrate satisfactory mid-term results of C-MVR with TACR for the treatment of coronary artery disease with moderate to severe or severe IMR.
7.Simultaneous Determination of Four Off-flavors in Freshwater Fish by Gas Chromatography-Mass Spectrometry Combined with Silica Solid Phase Extraction
Liang-Liang TIAN ; Dong-Mei HUANG ; Yuan WANG ; Xuan-Yun HUANG ; Yong-Fu SHI ; Hong-Li YE
Chinese Journal of Analytical Chemistry 2025;53(7):1158-1166
An effective method for simultaneously detecting four semivolatile earthy-musty odors in freshwater fish by gas chromatography-mass spectrometry(GC-MS)was developed.The concurrent extraction of geosmin(GSM),2-methylisoborneol(MIB),2-isopropyl-3-methoxypyrazine(IPMP),and 2-isobutyl-3-methoxypyrazine(IBMP)in fish tissue was conducted with n-hexane.The optimized QuEChERS material was implemented,and it was found that C18,primary secondary amine(PSA)and MgSO4 could adsorb the target analytes in n-hexane.So only the graphitized carbon black(GCB)could be used to purify the extraction.The adsorption rates of different materials for the four kinds of odors materials were explored in n-hexane and ethyl acetate.The experimental results revealed that the adsorption rates of silica for the four targets were 99.5%-100%in n-hexane and 0.7%-5.0%in ethyl acetate respectively.Then the silica solid phase extraction(SPE)method was utilized to eluent the compounds using 1.0 mL n-hexane/ethyl acetate in different proportions.The results of the comparative analysis demonstrated that n-hexane/ethyl acetate(4∶1,V/V)was the optimized eluent.Based on the obtained results,n-hexane extraction and GCB purification combined with silica SPE were used to isolate GSM,MIB,IPMP and IBMP from fish and the method was validated.It was found that the method showed good linearity in the range of 0.5-200 ng/mL,and with detection limits of 0.6 μg/kg for GSM and MIB,0.2 μg/kg for IPMP and IBMP.The limits of quantitation(LOQ)were 1.0 μg/kg for GSM and MIB,0.6 μg/kg for IPMP and IBMP.Good recoveries(77.5%-112.0%)and relative standard deviations(1.56%-9.42%)were also obtained.The use of silica SPE greatly mitigated the issue that the off-flavor compounds were easily lost in the gas blowing concentration process.There was no cross contamination in this method because the sample pretreatments were conducted separately,which was different with the most commonly used HS-SPME method for detecting semi-volatile substances.The sensitivity of this method was high enough to produce good quantitative results below the odor thresholds of the examined off-flavor compounds.
8.Development and validation of a recognition and classification system for portal hypertensive gastropathy based on deep learning
Haowen GU ; Jie YANG ; Yong XIAO ; Xinyue WAN ; Wei HU ; Xianmu XIE ; Dingpeng HUANG ; Chengming YAO ; Xinliang SHI ; Shiqian LIU ; Li HUANG ; Chi ZHANG ; Biqing ZHENG ; Mingkai CHEN
Chinese Journal of Digestive Endoscopy 2025;42(10):789-795
Objective:To develop a deep learning-based system for real-time recognition and classification of portal hypertensive gastropathy (PHG) and evaluate its ability to assist junior endoscopists.Methods:A total of 2 848 gastroscopy images from 832 patients with liver cirrhosis were selected from Digestive Endoscopy Center databases of Renmin Hospital of Wuhan University, Wuhan Hospital of Traditional Chinese and Western Medicine, and the Second Hospital of Jingzhou from January 2015 to October 2023. This system referred to 3 endoscopic features of Baveno Ⅱ scoring system. Three models were developed respectively for gastric antral vascular ectasia (GAVE), mosaic-like pattern (MLP), and red marks (RM). The specific classification references were as follows: (1) GAVE model: 0 no, 1 yes; (2) MLP model: 0 no, 1 mild, 2 severe; (3) RM model: 0 no, 1 isolated, 2 fused. The classification results for endoscopic characteristics of PHG of 3 endoscopy experts were taken as the gold standard. The yolov8-m model was used for training. The training dataset, validation dataset, and test dataset were allocated at a ratio of 8∶1∶1. The test dataset was used to evaluate the performance of models and their auxiliary effects on endoscopists. The accuracy, recall, precision, specificity and Kappa coefficient were calculated. Results:The accuracy, recall, specificity of GAVE model were 96.0% (48/50), 87.5% (7/8) and 97.6% (41/42). There was no significant difference between its accuracy and the gold standard ( χ2=316.226, P=1.000). The precision of GAVE1 and GAVE0 were 87.5% (7/8) and 97.6% (41/42) respectively. The accuracy of MLP model was 84.1% (132/157), and there was no significant difference compared with the gold standard ( χ2=3.286, P=0.193). The precision and recall of MLP2 were 88.2% (15/17) and 75.0% (15/20). The precision and recall of MLP1 were 77.9% (60/77) and 88.2% (60/68). The precision and recall of MLP0 were 90.5% (57/63) and 82.6% (57/69). The accuracy of RM model was 87.9% (123/140), and there was no significant difference compared with the gold standard ( χ2=2.891, P=0.409). The precision and recall of RM2 were 94.7% (18/19) and 78.3% (18/23). The precision and recall of RM1 were 72.2% (26/36) and 81.3% (26/32). The precision and recall of RM0 were 92.9% (79/85) and 92.9% (79/85). The mean accuracy of the three junior endoscopists, with and without the assistance of the GAVE model, MLP model, and RM model, respectively increased from 95.3% to 99.3%, from 83.9% to 91.9%, and from 81.9% to 83.1%. The overall consistency analysis of the 3 junior endoscopists with the gold standard indicated that the consistency of the GAVE model before and after assistance was extremely strong (both an overall Kappa of 1.000); the consistency before assistance of the MLP model was moderate (with an overall Kappa of 0.601), which increased to extremely strong after assistance (with an overall Kappa of 0.964); and the consistency of the RM model before and after assistance was also relatively strong (with an overall Kappa of 0.792 before and 0.798 after). Conclusion:The deep learning system accurately identifies and classifies PHG features and significantly enhances diagnostic performance of junior endoscopists.
9.Multimodal Data-Driven Prediction of Gynecological Surgery Duration
Yong HUANG ; Zhilin YONG ; Banghua WU ; Xueying ZHOU ; Xiaoling LANG ; Yuming LI ; Miye WANG ; Qingke SHI ; Li RAO
Journal of Sichuan University (Medical Sciences) 2025;56(5):1392-1398
Objective Focusing on gynecological surgery,we constructed a prediction model for surgical duration by extracting features from unstructured surgical planning texts and integrating multimodal data via artificial intelligence technology.Methods The clinical data of 34 614 patients who underwent gynecologic surgeries at West China Second University Hospital,Sichuan University between January 2022 and October 2024 were collected.An embedding-transformer model was constructed to convert surgical planning texts into a one-dimensional numerical feature,referred to as the step feature.The predictive value of the step feature was assessed by comparing the performance improvements of linear regression,random forest,eXtreme Gradient Boosting(XGBoost),support vector regression,K-nearest neighbor regression,and artificial neural network algorithms in two scenarios—with and without the step feature as an input.The out-of-sample prediction accuracy of the models was assessed using mean absolute error(MAE),root mean squared error(RMSE),and R-squared(R2).Furthermore,the model interpretability was examined using SHapley Additive exPlanations(SHAP)values.Results SHAP results showed that the step feature had the highest predictive contribution.Temporal factors in surgical scheduling also influenced gynecological surgery duration.The XGBoost model demonstrated optimal performance on the test set,significantly improving prediction accuracy with a 40.43%increase in R2,while reducing MAE and RMSE by 21.27%and 20.13%,respectively,compared to the baseline model without the step feature.Conclusion The embedding-transformer model developed in this study effectively extracts features from surgical planning texts and enhances the predictive performance of machine learning models.The XGBoost prediction model can assist hospital administrators in implementing more refined management of gynecological surgeries and improving the utilization efficiency of surgical resources.
10.Research on microscopic histochemical localization of Astragalus membranaceus and its application
Yong HUANG ; Shengchao SHI ; Jian ZHANG ; Chaofeng ZHANG
Journal of China Pharmaceutical University 2025;56(5):566-571
This study compared the chemical composition differences among various root tissues of Astragalus membranaceus through optimization of histochemical localization methods, thereby providing some reference for targeted extraction and enrichment of bioactive components. The chloral hydrate permeabilization technique was enhanced with sodium nitrite-aluminum nitrate-sodium hydroxide reagent for flavonoid staining. Comparative analyses of total saponins, total flavonoids, and astragaloside IV content in xylem, phloem, and periderm tissues were performed using ultraviolet-visible spectroscopy (UV-Vis) and high-performance liquid chromatography (HPLC). By optimizing staining protocols, interference from natural pigmentation was successfully eliminated, thereby enhancing the visualization of tissue-specific chemical distribution. Quantitative analysis revealed that the phloem contained the highest total saponins (6.92%), while the periderm exhibited peak total flavonoids (0.876%) and exclusive enrichment of astragaloside IV (0.850%). The refined histochemical localization method enables precise characterization of phytochemical distribution across root tissues, which can offer guidance for selective extraction and enrichment strategies, and comprehensively promote the utilization and development of Astragalus membranaceus resources.

Result Analysis
Print
Save
E-mail