1.Predicting intraoperative blood transfusion risk in hip fracture patients using explainable machine learning models
Fengting LU ; Xiaoming LI ; Dekui LI ; Xianyuan XIE ; Jiazhong WANG ; Qing YU ; Gan HUANG ; Jun SHEN
Chinese Journal of Blood Transfusion 2026;39(2):196-202
Objective: To investigate the factors influencing intraoperative blood transfusion in patients with hip fractures and to develop a machine learning (ML) model for predicting this risk. Methods: A total of 424 patients with hip fractures who underwent surgical treatment between November 2022 and March 2025 in our hospital were selected. Key feature variables of intraoperative blood transfusion risk were identified using the Boruta algorithm. Four different ML algorithms—support vector machine (SVM), linear discriminant analysis (LDA), mixed discriminant analysis (MDA), and extreme gradient boosting (XGBoost)—were used to develop predictive models for intraoperative blood transfusion risk. The predictive performance of the four ML models were evaluated using accuracy, precision, receiver operating characteristic (ROC) curves, precision-recall curves (PRC), precision-recall gain curves (PRGC), and F1 scores. Shapley additive interpretation (SHAP) was used to interpret the final model. Results: Among the 424 patients, 77(18.2%) received intraoperative blood transfusion. The Boruta algorithm identified albumin (ALB), activated partial thromboplastin time (APTT), types of anesthesia, types of fracture, and hemoglobin (Hb) as key feature variables for predicting intraoperative blood transfusion risk. In model evaluation, the SVM model outperforms the other three models across multiple metrics, including the area under the receiver operating characteristic curve (AUC), recall, recall gain, accuracy, precision, F1 score, and the area under the precision-recall curve (PRC-AUC). The SVM model, interpreted and visualized based on SHAP values, effectively predicted intraoperative blood transfusion risk in patients with hip fracture. A visual online application was developed based on the SVM model (https://pbo-nomogram.shinyapps.io/blood/). Conclusion: Preoperative low ALB and Hb levels, prolonged APTT, general anesthesia, and intertrochanteric fractures are risk factors for intraoperative blood transfusion in hip fracture patients. The risk prediction model for intraoperative blood transfusion constructed based on the SVM algorithm has optimal performance, which provides new ideas and methods for the clinical early identification of hip fracture patients with high transfusion risk and the implementation of targeted interventions.
2.Expert Consensus on the Ethical Requirements for Generative AI-Assisted Academic Writing
You-Quan BU ; Yong-Fu CAO ; Zeng-Yi CHANG ; Hong-Yu CHEN ; Xiao-Wei CHEN ; Yuan-Yuan CHEN ; Zhu-Cheng CHEN ; Rui DENG ; Jie DING ; Zhong-Kai FAN ; Guo-Quan GAO ; Xu GAO ; Lan HU ; Xiao-Qing HU ; Hong-Ti JIA ; Ying KONG ; En-Min LI ; Ling LI ; Yu-Hua LI ; Jun-Rong LIU ; Zhi-Qiang LIU ; Ya-Ping LUO ; Xue-Mei LV ; Yan-Xi PEI ; Xiao-Zhong PENG ; Qi-Qun TANG ; You WAN ; Yong WANG ; Ming-Xu WANG ; Xian WANG ; Guang-Kuan XIE ; Jun XIE ; Xiao-Hua YAN ; Mei YIN ; Zhong-Shan YU ; Chun-Yan ZHOU ; Rui-Fang ZHU
Chinese Journal of Biochemistry and Molecular Biology 2025;41(6):826-832
With the rapid development of generative artificial intelligence(GAI)technologies,their widespread application in academic research and writing is continuously expanding the boundaries of sci-entific inquiry.However,this trend has also raised a series of ethical and regulatory challenges,inclu-ding issues related to authorship,content authenticity,citation accuracy,and accountability.In light of the growing involvement of AI in generating academic content,establishing an open,controllable,and trustworthy ethical governance framework has become a key task for safeguarding research integrity and maintaining trust within the academic community.This expert consensus outlines ethical requirements across key stages of AI-assisted academic writing-including topic selection,data management,citation practices,and authorship attribution.It aims to clarify the boundaries and ethical obligations surrounding AI use in academic writing,ensuring that technological tools enhance efficiency without compromising in-tegrity.The goal is to provide guidance and institutional support for building a responsible and sustainable research ecosystem.
3.Construction and validation of machine learning-based dynamic early warning model for mortality risk in trauma-induced hypothermia patients
Yi-jing FU ; Jing YUAN ; Guan-jun LIU ; Qing-yan XIE ; Jia-meng XU ; Wei CHEN ; Guang ZHANG
Chinese Medical Equipment Journal 2025;46(3):9-14
Objective To propose a dynamic early warning model based on machine learning methods and validate its predi-ctive efficacy so as to achieve precise assessment and early warning of mortality risk in patients with traumatic hypothermia.Methods Firstly,a total of 480 patients who met inclusion criteria were retrospectively selected from the eICU database and randomly divided into training and test sets at an 8∶2 ratio.Secondly,physiological parameters were extracted from these patients,and five machine learning algorithms including XGBoost,AdaBoost,LightGBM,logistic regression(LR)and random forest(RF)were employed respectively to develop dynamic mortality risk warning models for traumatic hypothermia patients,utilizing a 1-hour observation window.Thirdly,receiver operating characteristic curves(ROC)were plotted using the test set data and the effects of different warning windows on the model performance were analyzed by calculating the AUC.Finally,the interpretability of the models was analyzed using the SHapley Additive exPlanations(SHAP)algorithm to elucidate the contribution of each feature to predictive performance.Results The optimal warning window for the dynamic warning model constructed using the eICU database was 12 hours,and in case of 12-hour warning window the logistic regression model achieved the highest AUC of 0.935 and showed optimal predictive performance.The results of the interpretability analysis by the SHAP algorithm showed that body temperature was the feature that had the greatest impact on the model results,and its reduction was positively correlated with the increased risk of death.Conclusion The machine learning-based dynamic warning model for mortality risk in traumatic hypothermia patients enables real-time dynamic risk assessment,providing robust support for clinicians to identify the patient's condition changes at an early stage and references for the adjustment of clinical treatment programs.[Chinese Medical Equipment Journal,2025,46(3):9-14]
4.Expert consensus on management of instrument separation in root canal therapy.
Yi FAN ; Yuan GAO ; Xiangzhu WANG ; Bing FAN ; Zhi CHEN ; Qing YU ; Ming XUE ; Xiaoyan WANG ; Zhengwei HUANG ; Deqin YANG ; Zhengmei LIN ; Yihuai PAN ; Jin ZHAO ; Jinhua YU ; Zhuo CHEN ; Sijing XIE ; He YUAN ; Kehua QUE ; Shuang PAN ; Xiaojing HUANG ; Jun LUO ; Xiuping MENG ; Jin ZHANG ; Yi DU ; Lei ZHANG ; Hong LI ; Wenxia CHEN ; Jiayuan WU ; Xin XU ; Jing ZOU ; Jiyao LI ; Dingming HUANG ; Lei CHENG ; Tiemei WANG ; Benxiang HOU ; Xuedong ZHOU
International Journal of Oral Science 2025;17(1):46-46
Instrument separation is a critical complication during root canal therapy, impacting treatment success and long-term tooth preservation. The etiology of instrument separation is multifactorial, involving the intricate anatomy of the root canal system, instrument-related factors, and instrumentation techniques. Instrument separation can hinder thorough cleaning, shaping, and obturation of the root canal, posing challenges to successful treatment outcomes. Although retrieval of separated instrument is often feasible, it carries risks including perforation, excessive removal of tooth structure and root fractures. Effective management of separated instruments requires a comprehensive understanding of the contributing factors, meticulous preoperative assessment, and precise evaluation of the retrieval difficulty. The application of appropriate retrieval techniques is essential to minimize complications and optimize clinical outcomes. The current manuscript provides a framework for understanding the causes, risk factors, and clinical management principles of instrument separation. By integrating effective strategies, endodontists can enhance decision-making, improve endodontic treatment success and ensure the preservation of natural dentition.
Humans
;
Root Canal Therapy/adverse effects*
;
Consensus
;
Root Canal Preparation/adverse effects*
5.Correlation between dynamic contrast-enhanced MRI imaging and clinical pathological features of invasive breast cancer and lymphovascular invasion
Shi-Qi GUO ; Yu-Jiao XIE ; Qing-Yang LI ; Si-Yi CHEN ; Jia-Hong SUN ; Zhao-Feng GAO ; Jun-Qing LIANG ; Yu-Hui CHEN ; Bao-Shi BAO ; Li ZHU ; Jian-Dong WANG
Medical Journal of Chinese People's Liberation Army 2025;50(7):847-854
Objective To explore the relationship between dynamic contrast-enhanced magnetic resonance imaging(DCE-MRI)and clinical pathological features of invasive breast cancer and lymphovascular invasion(LVI).Methods Imaging and clinical pathological data were retrospectively collected from 508 patients with invasive breast cancer who underwent breast DCE-MRI at the First Medical Center of Chinese PLA General Hospital from January 2019 to August 2021.Patients were divided into the LVI-positive(LVI+)group(n=79)and LVI-negative(LVI-)group(n=429)based on postoperative pathological results.Univariate and multivariate logistic regression analyses were used to identify risk factors for LVI.Results Compared with LVI-group,LVI+group had a higher proportion of patients aged<45 years(44.3%vs.27.0%,P=0.002),non-mass-like enhancement(NME)(31.7%vs.17.7%,P=0.004),Ki-67 expression rate(40.0%vs.30.0%,P<0.001),high Ki-67 expression(94.9%vs.78.1%,P=0.001),Luminal B subtype(76.0%vs.60.1%,P=0.008),and positive axillary lymph nodes rate(72.2%vs.31.5%,P<0.001),while the proportion of Luminal A subtype was lower(2.5%vs.21.5%,P<0.001).Univariate and multivariate logistic regression analyses showed that age≥45 years(OR=0.468,95%CI 0.280-0.783,P=0.004)was an independent protective factor for LVI,while NME(OR=1.987,95%CI 1.126-3.444,P=0.016)was an independent risk factor.Compared with Luminal A subtype,patients with Luminal B subtype(OR=10.482,95%CI 3.164-64.923,P=0.001),HER-2 overexpression subtype(OR=11.571,95%CI 2.755-79.341,P=0.003)and triple-negative subtypes(OR=8.433,95%CI 1.985-57.908,P=0.009)had a higher risk of LVI.Conclusions Age≥45 years is an independent protective factor for LVI,while NME is an independent risk factor.Among molecular subtypes,patients with Luminal B,HER-2 overexpression and triple-negative subtypes have a higher risk of LVI compared with the Luminal A subtype.
6.Percutaneous coronary intervention vs . medical therapy in patients on dialysis with coronary artery disease in China.
Enmin XIE ; Yaxin WU ; Zixiang YE ; Yong HE ; Hesong ZENG ; Jianfang LUO ; Mulei CHEN ; Wenyue PANG ; Yanmin XU ; Chuanyu GAO ; Xiaogang GUO ; Lin CAI ; Qingwei JI ; Yining YANG ; Di WU ; Yiqiang YUAN ; Jing WAN ; Yuliang MA ; Jun ZHANG ; Zhimin DU ; Qing YANG ; Jinsong CHENG ; Chunhua DING ; Xiang MA ; Chunlin YIN ; Zeyuan FAN ; Qiang TANG ; Yue LI ; Lihua SUN ; Chengzhi LU ; Jufang CHI ; Zhuhua YAO ; Yanxiang GAO ; Changan YU ; Jingyi REN ; Jingang ZHENG
Chinese Medical Journal 2025;138(3):301-310
BACKGROUND:
The available evidence regarding the benefits of percutaneous coronary intervention (PCI) on patients receiving dialysis with coronary artery disease (CAD) is limited and inconsistent. This study aimed to evaluate the association between PCI and clinical outcomes as compared with medical therapy alone in patients undergoing dialysis with CAD in China.
METHODS:
This multicenter, retrospective study was conducted in 30 tertiary medical centers across 12 provinces in China from January 2015 to June 2021 to include patients on dialysis with CAD. The primary outcome was major adverse cardiovascular events (MACE), defined as a composite of cardiovascular death, non-fatal myocardial infarction, and non-fatal stroke. Secondary outcomes included all-cause death, the individual components of MACE, and Bleeding Academic Research Consortium criteria types 2, 3, or 5 bleeding. Multivariable Cox proportional hazard models were used to assess the association between PCI and outcomes. Inverse probability of treatment weighting (IPTW) and propensity score matching (PSM) were performed to account for potential between-group differences.
RESULTS:
Of the 1146 patients on dialysis with significant CAD, 821 (71.6%) underwent PCI. After a median follow-up of 23.0 months, PCI was associated with a 43.0% significantly lower risk for MACE (33.9% [ n = 278] vs . 43.7% [ n = 142]; adjusted hazards ratio 0.57, 95% confidence interval 0.45-0.71), along with a slightly increased risk for bleeding outcomes that did not reach statistical significance (11.1% vs . 8.3%; adjusted hazards ratio 1.31, 95% confidence interval, 0.82-2.11). Furthermore, PCI was associated with a significant reduction in all-cause and cardiovascular mortalities. Subgroup analysis did not modify the association of PCI with patient outcomes. These primary findings were consistent across IPTW, PSM, and competing risk analyses.
CONCLUSION
This study indicated that PCI in patients on dialysis with CAD was significantly associated with lower MACE and mortality when comparing with those with medical therapy alone, albeit with a slightly increased risk for bleeding events that did not reach statistical significance.
Humans
;
Percutaneous Coronary Intervention/methods*
;
Male
;
Female
;
Coronary Artery Disease/drug therapy*
;
Retrospective Studies
;
Renal Dialysis/methods*
;
Middle Aged
;
Aged
;
China
;
Proportional Hazards Models
;
Treatment Outcome
7.Carvedilol to prevent hepatic decompensation of cirrhosis in patients with clinically significant portal hypertension stratified by new non-invasive model (CHESS2306)
Chuan LIU ; Hong YOU ; Qing-Lei ZENG ; Yu Jun WONG ; Bingqiong WANG ; Ivica GRGUREVIC ; Chenghai LIU ; Hyung Joon YIM ; Wei GOU ; Bingtian DONG ; Shenghong JU ; Yanan GUO ; Qian YU ; Masashi HIROOKA ; Hirayuki ENOMOTO ; Amr Shaaban HANAFY ; Zhujun CAO ; Xiemin DONG ; Jing LV ; Tae Hyung KIM ; Yohei KOIZUMI ; Yoichi HIASA ; Takashi NISHIMURA ; Hiroko IIJIMA ; Chuanjun XU ; Erhei DAI ; Xiaoling LAN ; Changxiang LAI ; Shirong LIU ; Fang WANG ; Ying GUO ; Jiaojian LV ; Liting ZHANG ; Yuqing WANG ; Qing XIE ; Chuxiao SHAO ; Zhensheng LIU ; Federico RAVAIOLI ; Antonio COLECCHIA ; Jie LI ; Gao-Jun TENG ; Xiaolong QI
Clinical and Molecular Hepatology 2025;31(1):105-118
Background:
s/Aims: Non-invasive models stratifying clinically significant portal hypertension (CSPH) are limited. Herein, we developed a new non-invasive model for predicting CSPH in patients with compensated cirrhosis and investigated whether carvedilol can prevent hepatic decompensation in patients with high-risk CSPH stratified using the new model.
Methods:
Non-invasive risk factors of CSPH were identified via systematic review and meta-analysis of studies involving patients with hepatic venous pressure gradient (HVPG). A new non-invasive model was validated for various performance aspects in three cohorts, i.e., a multicenter HVPG cohort, a follow-up cohort, and a carvediloltreating cohort.
Results:
In the meta-analysis with six studies (n=819), liver stiffness measurement and platelet count were identified as independent risk factors for CSPH and were used to develop the new “CSPH risk” model. In the HVPG cohort (n=151), the new model accurately predicted CSPH with cutoff values of 0 and –0.68 for ruling in and out CSPH, respectively. In the follow-up cohort (n=1,102), the cumulative incidences of decompensation events significantly differed using the cutoff values of <–0.68 (low-risk), –0.68 to 0 (medium-risk), and >0 (high-risk). In the carvediloltreated cohort, patients with high-risk CSPH treated with carvedilol (n=81) had lower rates of decompensation events than non-selective beta-blockers untreated patients with high-risk CSPH (n=613 before propensity score matching [PSM], n=162 after PSM).
Conclusions
Treatment with carvedilol significantly reduces the risk of hepatic decompensation in patients with high-risk CSPH stratified by the new model.
8.Carvedilol to prevent hepatic decompensation of cirrhosis in patients with clinically significant portal hypertension stratified by new non-invasive model (CHESS2306)
Chuan LIU ; Hong YOU ; Qing-Lei ZENG ; Yu Jun WONG ; Bingqiong WANG ; Ivica GRGUREVIC ; Chenghai LIU ; Hyung Joon YIM ; Wei GOU ; Bingtian DONG ; Shenghong JU ; Yanan GUO ; Qian YU ; Masashi HIROOKA ; Hirayuki ENOMOTO ; Amr Shaaban HANAFY ; Zhujun CAO ; Xiemin DONG ; Jing LV ; Tae Hyung KIM ; Yohei KOIZUMI ; Yoichi HIASA ; Takashi NISHIMURA ; Hiroko IIJIMA ; Chuanjun XU ; Erhei DAI ; Xiaoling LAN ; Changxiang LAI ; Shirong LIU ; Fang WANG ; Ying GUO ; Jiaojian LV ; Liting ZHANG ; Yuqing WANG ; Qing XIE ; Chuxiao SHAO ; Zhensheng LIU ; Federico RAVAIOLI ; Antonio COLECCHIA ; Jie LI ; Gao-Jun TENG ; Xiaolong QI
Clinical and Molecular Hepatology 2025;31(1):105-118
Background:
s/Aims: Non-invasive models stratifying clinically significant portal hypertension (CSPH) are limited. Herein, we developed a new non-invasive model for predicting CSPH in patients with compensated cirrhosis and investigated whether carvedilol can prevent hepatic decompensation in patients with high-risk CSPH stratified using the new model.
Methods:
Non-invasive risk factors of CSPH were identified via systematic review and meta-analysis of studies involving patients with hepatic venous pressure gradient (HVPG). A new non-invasive model was validated for various performance aspects in three cohorts, i.e., a multicenter HVPG cohort, a follow-up cohort, and a carvediloltreating cohort.
Results:
In the meta-analysis with six studies (n=819), liver stiffness measurement and platelet count were identified as independent risk factors for CSPH and were used to develop the new “CSPH risk” model. In the HVPG cohort (n=151), the new model accurately predicted CSPH with cutoff values of 0 and –0.68 for ruling in and out CSPH, respectively. In the follow-up cohort (n=1,102), the cumulative incidences of decompensation events significantly differed using the cutoff values of <–0.68 (low-risk), –0.68 to 0 (medium-risk), and >0 (high-risk). In the carvediloltreated cohort, patients with high-risk CSPH treated with carvedilol (n=81) had lower rates of decompensation events than non-selective beta-blockers untreated patients with high-risk CSPH (n=613 before propensity score matching [PSM], n=162 after PSM).
Conclusions
Treatment with carvedilol significantly reduces the risk of hepatic decompensation in patients with high-risk CSPH stratified by the new model.
9.Carvedilol to prevent hepatic decompensation of cirrhosis in patients with clinically significant portal hypertension stratified by new non-invasive model (CHESS2306)
Chuan LIU ; Hong YOU ; Qing-Lei ZENG ; Yu Jun WONG ; Bingqiong WANG ; Ivica GRGUREVIC ; Chenghai LIU ; Hyung Joon YIM ; Wei GOU ; Bingtian DONG ; Shenghong JU ; Yanan GUO ; Qian YU ; Masashi HIROOKA ; Hirayuki ENOMOTO ; Amr Shaaban HANAFY ; Zhujun CAO ; Xiemin DONG ; Jing LV ; Tae Hyung KIM ; Yohei KOIZUMI ; Yoichi HIASA ; Takashi NISHIMURA ; Hiroko IIJIMA ; Chuanjun XU ; Erhei DAI ; Xiaoling LAN ; Changxiang LAI ; Shirong LIU ; Fang WANG ; Ying GUO ; Jiaojian LV ; Liting ZHANG ; Yuqing WANG ; Qing XIE ; Chuxiao SHAO ; Zhensheng LIU ; Federico RAVAIOLI ; Antonio COLECCHIA ; Jie LI ; Gao-Jun TENG ; Xiaolong QI
Clinical and Molecular Hepatology 2025;31(1):105-118
Background:
s/Aims: Non-invasive models stratifying clinically significant portal hypertension (CSPH) are limited. Herein, we developed a new non-invasive model for predicting CSPH in patients with compensated cirrhosis and investigated whether carvedilol can prevent hepatic decompensation in patients with high-risk CSPH stratified using the new model.
Methods:
Non-invasive risk factors of CSPH were identified via systematic review and meta-analysis of studies involving patients with hepatic venous pressure gradient (HVPG). A new non-invasive model was validated for various performance aspects in three cohorts, i.e., a multicenter HVPG cohort, a follow-up cohort, and a carvediloltreating cohort.
Results:
In the meta-analysis with six studies (n=819), liver stiffness measurement and platelet count were identified as independent risk factors for CSPH and were used to develop the new “CSPH risk” model. In the HVPG cohort (n=151), the new model accurately predicted CSPH with cutoff values of 0 and –0.68 for ruling in and out CSPH, respectively. In the follow-up cohort (n=1,102), the cumulative incidences of decompensation events significantly differed using the cutoff values of <–0.68 (low-risk), –0.68 to 0 (medium-risk), and >0 (high-risk). In the carvediloltreated cohort, patients with high-risk CSPH treated with carvedilol (n=81) had lower rates of decompensation events than non-selective beta-blockers untreated patients with high-risk CSPH (n=613 before propensity score matching [PSM], n=162 after PSM).
Conclusions
Treatment with carvedilol significantly reduces the risk of hepatic decompensation in patients with high-risk CSPH stratified by the new model.
10.International clinical practice guideline on the use of traditional Chinese medicine for functional dyspepsia (2025).
Sheng-Sheng ZHANG ; Lu-Qing ZHAO ; Xiao-Hua HOU ; Zhao-Xiang BIAN ; Jian-Hua ZHENG ; Hai-He TIAN ; Guan-Hu YANG ; Won-Sook HONG ; Yu-Ying HE ; Li LIU ; Hong SHEN ; Yan-Ping LI ; Sheng XIE ; Jin SHU ; Bin-Fang ZENG ; Jun-Xiang LI ; Zhen LIU ; Zheng-Hua XIAO ; Jing-Dong XIAO ; Pei-Yong ZHENG ; Shao-Gang HUANG ; Sheng-Liang CHEN ; Gui-Jun FEI
Journal of Integrative Medicine 2025;23(5):502-518
Functional dyspepsia (FD), characterized by persistent or recurrent dyspeptic symptoms without identifiable organic, systemic or metabolic causes, is an increasingly recognized global health issue. The objective of this guideline is to equip clinicians and nursing professionals with evidence-based strategies for the management and treatment of adult patients with FD using traditional Chinese medicine (TCM). The Guideline Development Group consulted existing TCM consensus documents on FD and convened a panel of 35 clinicians to generate initial clinical queries. To address these queries, a systematic literature search was conducted across PubMed, EMBASE, the Cochrane Library, China National Knowledge Infrastructure (CNKI), VIP Database, China Biology Medicine (SinoMed) Database, Wanfang Database, Traditional Medicine Research Data Expanded (TMRDE), and the Traditional Chinese Medical Literature Analysis and Retrieval System (TCMLARS). The evidence from the literature was critically appraised using the Grading of Recommendations Assessment, Development, and Evaluation (GRADE) approach. The strength of the recommendations was ascertained through a consensus-building process involving TCM and allopathic medicine experts, methodologists, pharmacologists, nursing specialists, and health economists, leveraging their collective expertise and empirical knowledge. The guideline comprises a total of 43 evidence-informed recommendations that span a range of clinical aspects, including the pathogenesis according to TCM, diagnostic approaches, therapeutic interventions, efficacy assessments, and prognostic considerations. Please cite this article as: Zhang SS, Zhao LQ, Hou XH, Bian ZX, Zheng JH, Tian HH, Yang GH, Hong WS, He YY, Liu L, Shen H, Li YP, Xie S, Shu J, Zeng BF, Li JX, Liu Z, Xiao ZH, Xiao JD, Zheng PY, Huang SG, Chen SL, Fei GJ. International clinical practice guideline on the use of traditional Chinese medicine for functional dyspepsia (2025). J Integr Med. 2025; 23(5):502-518.
Dyspepsia/drug therapy*
;
Humans
;
Medicine, Chinese Traditional/methods*
;
Practice Guidelines as Topic
;
Drugs, Chinese Herbal/therapeutic use*

Result Analysis
Print
Save
E-mail