1.An Attention-weighted Tri-modal Ultrasound Network (TUS-Net) for Screening of Atypical Hepatocellular Carcinoma From LR-M Liver Nodules
He-Chong ZHANG ; Liang-Hui HUANG ; Xue-Hua WANG ; Shang-Lin JIANG ; Ying-Ying CHEN ; Ya-Guang ZENG ; Wei ZHENG
Progress in Biochemistry and Biophysics 2026;53(5):1485-1498
ObjectiveDiscriminating atypical hepatocellular carcinoma (HCC) from other malignancies in liver nodules classified as Liver Imaging Reporting and Data System category M (LR-M) remains a significant diagnostic challenge on conventional ultrasound examination. The LR-M category, originally intended to capture non-HCC malignancies, paradoxically contains up to 63% of atypical HCCs that deviate from classic enhancement patterns, leading to potential misdiagnosis and suboptimal treatment planning. While deep learning has shown promise in HCC diagnosis, most existing models rely exclusively on single-modality ultrasound, overlooking the diagnostic benefits of integrating complementary information from multiple imaging sources. To address this gap, we propose a novel attention-weighted tri-modal ultrasound network (TUS-Net) that integrates contrast-enhanced ultrasound (CEUS), B-mode ultrasound (BUS), and time-intensity curves (TICs) to improve diagnostic accuracy for these clinically challenging lesions. MethodsOur framework incorporates a three-dimensional convolutional neural network (C3D) backbone to extract spatiotemporal features from CEUS videos, capturing dynamic vascular patterns critical for lesion characterization. To effectively fuse complementary modalities, we introduce a dual-channel feature fusion module (DCFFM) that adaptively combines features from CEUS and BUS through channel-wise attention mechanisms, allowing the model to dynamically weigh the contribution of each modality based on diagnostic relevance. Additionally, we propose a temporal intensity feature fusion module (TIFFM) that leverages quantitative hemodynamic information from TICs to guide the model’s attention toward diagnostically critical temporal phases, such as arterial wash-in and portal venous washout. The model is further enhanced by automated lesion localization using YOLOX and class activation mapping for interpretability, ensuring that predictions align with clinically meaningful imaging features. ResultsEvaluated on a tri-modal ultrasound dataset comprising 161 patients with pathologically confirmed LR-M nodules (131 atypical HCC and 30 non-HCC malignancies), our model achieved an accuracy of 86.83%, a sensitivity of 92.50%, a specificity of 75.50%, and an AUC of 89.32% in screening atypical HCC. Compared to single-modality baselines, TUS-Net demonstrated superior specificity, a clinically critical metric given the higher risk associated with misclassifying non-HCC malignancies. Ablation studies confirmed the contribution of each module, with the full model outperforming both standard C3D and 3D ResNet backbones integrated with attention mechanisms. A reader study involving junior and senior radiologists further validated the clinical utility of AI assistance, showing consistent improvements in specificity and inter-reader consistency, particularly for less experienced clinicians. ConclusionThese results surpass existing benchmark models and demonstrate the potential of our approach to enhance diagnostic precision in clinically specific cases. By intelligently fusing multi-modal ultrasound data with attention-guided mechanisms, TUS-Net offers a reliable and interpretable tool that holds promise for improving the non-invasive diagnosis of atypical HCC in challenging LR-M liver nodules.
2.An Attention-weighted Tri-modal Ultrasound Network (TUS-Net) for Screening of Atypical Hepatocellular Carcinoma From LR-M Liver Nodules
He-Chong ZHANG ; Liang-Hui HUANG ; Xue-Hua WANG ; Shang-Lin JIANG ; Ying-Ying CHEN ; Ya-Guang ZENG ; Wei ZHENG
Progress in Biochemistry and Biophysics 2026;53(5):1485-1498
ObjectiveDiscriminating atypical hepatocellular carcinoma (HCC) from other malignancies in liver nodules classified as Liver Imaging Reporting and Data System category M (LR-M) remains a significant diagnostic challenge on conventional ultrasound examination. The LR-M category, originally intended to capture non-HCC malignancies, paradoxically contains up to 63% of atypical HCCs that deviate from classic enhancement patterns, leading to potential misdiagnosis and suboptimal treatment planning. While deep learning has shown promise in HCC diagnosis, most existing models rely exclusively on single-modality ultrasound, overlooking the diagnostic benefits of integrating complementary information from multiple imaging sources. To address this gap, we propose a novel attention-weighted tri-modal ultrasound network (TUS-Net) that integrates contrast-enhanced ultrasound (CEUS), B-mode ultrasound (BUS), and time-intensity curves (TICs) to improve diagnostic accuracy for these clinically challenging lesions. MethodsOur framework incorporates a three-dimensional convolutional neural network (C3D) backbone to extract spatiotemporal features from CEUS videos, capturing dynamic vascular patterns critical for lesion characterization. To effectively fuse complementary modalities, we introduce a dual-channel feature fusion module (DCFFM) that adaptively combines features from CEUS and BUS through channel-wise attention mechanisms, allowing the model to dynamically weigh the contribution of each modality based on diagnostic relevance. Additionally, we propose a temporal intensity feature fusion module (TIFFM) that leverages quantitative hemodynamic information from TICs to guide the model’s attention toward diagnostically critical temporal phases, such as arterial wash-in and portal venous washout. The model is further enhanced by automated lesion localization using YOLOX and class activation mapping for interpretability, ensuring that predictions align with clinically meaningful imaging features. ResultsEvaluated on a tri-modal ultrasound dataset comprising 161 patients with pathologically confirmed LR-M nodules (131 atypical HCC and 30 non-HCC malignancies), our model achieved an accuracy of 86.83%, a sensitivity of 92.50%, a specificity of 75.50%, and an AUC of 89.32% in screening atypical HCC. Compared to single-modality baselines, TUS-Net demonstrated superior specificity, a clinically critical metric given the higher risk associated with misclassifying non-HCC malignancies. Ablation studies confirmed the contribution of each module, with the full model outperforming both standard C3D and 3D ResNet backbones integrated with attention mechanisms. A reader study involving junior and senior radiologists further validated the clinical utility of AI assistance, showing consistent improvements in specificity and inter-reader consistency, particularly for less experienced clinicians. ConclusionThese results surpass existing benchmark models and demonstrate the potential of our approach to enhance diagnostic precision in clinically specific cases. By intelligently fusing multi-modal ultrasound data with attention-guided mechanisms, TUS-Net offers a reliable and interpretable tool that holds promise for improving the non-invasive diagnosis of atypical HCC in challenging LR-M liver nodules.
3.Identification of GSK3 family and regulatory effects of brassinolide on growth and development of Nardostachys jatamansi.
Yu-Yan LEI ; Zheng MA ; Jing WEI ; Wen-Bing LI ; Ying LI ; Zheng-Ming YANG ; Shao-Shan ZHANG ; Jing-Qiu FENG ; Hua-Chun SHENG ; Yuan LIU
China Journal of Chinese Materia Medica 2025;50(2):395-403
This study identified 8 members including NjBIN2 of the GSK3 family in Nardostachys jatamansi by bioinformatics analysis. Moreover, the phylogenetic tree revealed that the GKS3 family members of N. jatamansi had a close relationship with those of Arabidopsis. RT-qPCR results showed that NjBIN2 presented a tissue-specific expression pattern with the highest expression in roots, suggesting that NjBIN2 played a role in root growth and development. In addition, the application of epibrassinolide or the brassinosteroid(BR) synthesis inhibitor(brassinazole) altered the expression pattern of NjBIN2 and influenced the photomorphogenesis(cotyledon opening) and root development of N. jatamansi, which provided direct evidence about the functions of NjBIN2. In conclusion, this study highlights the roles of BIN2 in regulating the growth and development of N. jatamansi by analyzing the expression pattern and biological function of NjBIN2. It not only enriches the understanding about the regulatory mechanism of the growth and development of N. jatamansi but also provides a theoretical basis and potential gene targets for molecular breeding of N. jatamansi with improved quality in the future.
Brassinosteroids/metabolism*
;
Steroids, Heterocyclic/metabolism*
;
Gene Expression Regulation, Plant/drug effects*
;
Plant Proteins/metabolism*
;
Phylogeny
;
Nardostachys/metabolism*
;
Plant Growth Regulators/pharmacology*
;
Plant Roots/drug effects*
4.Exploration of pharmacodynamic substances and potential mechanisms of Huazhuo Sanjie Chubi Decoction in treatment of gouty arthritis based on UPLC-Q-Exactive Orbitrap-MS technology and network pharmacology.
Yan XIAO ; Ting ZHANG ; Ying-Jie ZHANG ; Bin HUANG ; Peng CHEN ; Xiao-Hua CHEN ; Ming-Qing HUANG ; Xue-Ting CHEN ; You-Xin SU ; Jie-Mei GUO
China Journal of Chinese Materia Medica 2025;50(2):444-488
Based on ultra-high performance liquid chromatography-quadrupole-Exactive Orbitrap mass spectrometry(UPLC-Q-Exactive Orbitrap-MS) technology and network pharmacology, this study explored the pharmacodynamic substances and potential mechanisms of Huazhuo Sanjie Chubi Decoction in the treatment of gouty arthritis(GA). UPLC-Q-Exactive Orbitrap-MS technology was used to identify the components in Huazhuo Sanjie Chubi Decoction, and the qualitative analysis of its active ingredients was carried out, with a total of 184 active ingredients identified. A total of 897 active ingredient targets were screened through the PharmMapper database, and 491 GA-related disease targets were obtained from the OMIM, GeneCards, CTD databases. After Venn analysis, 60 intersecting targets were obtained. The component target-GA target network was constructed through the Cytoscape platform, and the STRING database was used to construct a protein-protein interaction network, with 16 core targets screened. The core targets were subjected to Gene Ontology(GO) and Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway enrichment analyses, and the component-target-pathway network was constructed. It was found that the main active ingredients of the formula for the treatment of GA were phenols, flavonoids, alkaloids, and terpenoids, and the key targets were SRC, MMP3, MMP9, REN, ALB, IGF1R, PPARG, MAPK1, HPRT1, and CASP1. Through GO analysis, it was found that the treatment of GA mainly involved biological processes such as lipid response, bacterial response, and biostimulus response. KEGG analysis showed that the pathways related to the treatment of GA included lipids and atherosclerosis, neutrophil extracellular traps(NETs), IL-17, and so on. In summary, phenols, flavonoids, alkaloids, and terpenoids may be the core pharmacodynamic substances of Huazhuo Sanjie Chubi Decoction in the treatment of GA, and the pharmacodynamic mechanism may be related to SRC, MMP3, MMP9, and other targets, as well as lipids and atherosclerosis, NETs, IL-17, and other pathways.
Drugs, Chinese Herbal/therapeutic use*
;
Network Pharmacology
;
Arthritis, Gouty/metabolism*
;
Chromatography, High Pressure Liquid/methods*
;
Humans
;
Mass Spectrometry/methods*
;
Protein Interaction Maps/drug effects*
5."Component-effect" correlations in traditional Chinese medicine from holistic view: taking discovery of gintonin from ginseng as an example.
Xin-Ming YU ; Chen-Yu YU ; Hua-Ying WANG ; Wei-Sheng YUE ; Zhu-Bin ZHANG ; Wei WU ; Xiao-Bin JIA ; Bing YANG ; Liang FENG
China Journal of Chinese Materia Medica 2025;50(7):2001-2012
The holistic view is the key in the study of traditional Chinese medicine(TCM). The component structure theory is based on the holistic view to investigate the correlation between material basis and efficiency, which enriches the holistic "component-effect" research of TCM. Gintonin is a newly isolated non-saponin component of ginseng. Compared to ginsenosides, gintonin has many different pharmacological activities, and it provides new knowledge for the holistic research of ginseng. Thus, taking the discovery of gintonin from ginseng as an example, this paper explored the linkage between ginsenosides and gintonin from the perspective of "component-effect" correlations and systematically sorted out the similarities and differences between them in terms of structural characteristics, modes of action, and pharmacological activities. Starting from the collaborative interaction of TCM compounds, the study discussed the application and value of the holistic view in TCM "component-effect" research in the light of the component structure theory to provide new thoughts for the development of modern TCM research.
Panax/chemistry*
;
Drugs, Chinese Herbal/pharmacology*
;
Medicine, Chinese Traditional
;
Humans
;
Ginsenosides/pharmacology*
;
Animals
6.Randomized, double-blind, parallel-controlled, multicenter, equivalence clinical trial of Jiuwei Xifeng Granules(Os Draconis replaced by Ostreae Concha) for treating tic disorder in children.
Qiu-Han CAI ; Cheng-Liang ZHONG ; Si-Yuan HU ; Xin-Min LI ; Zhi-Chun XU ; Hui CHEN ; Ying HUA ; Jun-Hong WANG ; Ji-Hong TANG ; Bing-Xiang MA ; Xiu-Xia WANG ; Ai-Zhen WANG ; Meng-Qing WANG ; Wei ZHANG ; Chun WANG ; Yi-Qun TENG ; Yi-Hui SHAN ; Sheng-Xuan GUO
China Journal of Chinese Materia Medica 2025;50(6):1699-1705
Jiuwei Xifeng Granules have become a Chinese patent medicine in the market. Because the formula contains Os Draconis, a top-level protected fossil of ancient organisms, the formula was to be improved by replacing Os Draconis with Ostreae Concha. To evaluate whether the improved formula has the same effectiveness and safety as the original formula, a randomized, double-blind, parallel-controlled, equivalence clinical trial was conducted. This study enrolled 288 tic disorder(TD) of children and assigned them into two groups in 1∶1. The treatment group and control group took the modified formula and original formula, respectively. The treatment lasted for 6 weeks, and follow-up visits were conducted at weeks 2, 4, and 6. The primary efficacy endpoint was the difference in Yale global tic severity scale(YGTSS)-total tic severity(TTS) score from baseline after 6 weeks of treatment. The results showed that after 6 weeks of treatment, the declines in YGTSS-TSS score showed no statistically significant difference between the two groups. The difference in YGTSS-TSS score(treatment group-control group) and the 95%CI of the full analysis set(FAS) were-0.17[-1.42, 1.08] and those of per-protocol set(PPS) were 0.29[-0.97, 1.56], which were within the equivalence boundary [-3, 3]. The equivalence test was therefore concluded. The two groups showed no significant differences in the secondary efficacy endpoints of effective rate for TD, total score and factor scores of YGTSS, clinical global impressions-severity(CGI-S) score, traditional Chinese medicine(TCM) response rate, or symptom disappearance rate, and thus a complete evidence chain with the primary outcome was formed. A total of 6 adverse reactions were reported, including 4(2.82%) cases in the treatment group and 2(1.41%) cases in the control group, which showed no statistically significant difference between the two groups. No serious suspected unexpected adverse reactions were reported, and no laboratory test results indicated serious clinically significant abnormalities. The results support the replacement of Os Draconis by Ostreae Concha in the original formula, and the efficacy and safety of the modified formula are consistent with those of the original formula.
Adolescent
;
Child
;
Child, Preschool
;
Female
;
Humans
;
Male
;
Double-Blind Method
;
Drugs, Chinese Herbal/therapeutic use*
;
Tic Disorders/drug therapy*
;
Treatment Outcome
7.Evaluation of nutritional value of three kinds of medicinal snakes based on content of 15 amino acids.
Xi WANG ; Ye-Yuan LIN ; Wen-Ting ZHONG ; Zhi-Guo MA ; Meng-Hua WU ; Hui CAO ; Ying ZHANG
China Journal of Chinese Materia Medica 2025;50(9):2411-2421
A high-performance liquid chromatography method using pre-column derivatization with 6-aminoquinolyl-N-hydroxysuccinimidyl carbamate was developed to determine the content of 15 amino acids in the medicinal snakes Bungarus Parvus, Agkistrodon, and Zaocys. The results showed that the total amino acid(TAA) content ranged from 277.13 to 515.05 mg·g~(-1), with the top four amino acids in all three species being glutamic acid(Glu), glycine(Gly), aspartic acid(Asp), and lysine(Lys). The essential amino acid(EAA) content ranged from 74.56 to 203.94 mg·g~(-1), with Agkistrodon exhibiting the highest content. The non-essential amino acid(NEAA), semi-essential amino acid(semi-EAA), and medicinal amino acid(MAA) content ranged from 189.06 to 318.23, 12.89 to 33.53, and 179.83 to 342.33 mg·g~(-1), respectively, with Zaocys having the highest content in these categories. Amino acid nutritional value was evaluated using the amino acid ratio(RAA), amino acid ratio coefficient(RCAA), and amino acid ratio coefficient score(SRCAA), and the results indicated that all three medicinal snakes possessed good nutritional value. The amino acid composition was similar across the species, though significant differences in content were observed. Based on these differences, an orthogonal partial least squares-discriminant analysis(OPLS-DA) model was established, which could clearly distinguish between the three medicinal snake species. The key differences in amino acid content included Gly, tyrosine(Tyr), Glu, and serine(Ser), which may be related to the observed clinical application differences among the species. Further research into the mechanisms of these differential amino acids is expected to provide more insights into the clinical application disparities of these three medicinal snake species.
Amino Acids/chemistry*
;
Animals
;
Nutritive Value
;
Chromatography, High Pressure Liquid
;
Snakes/classification*
;
Bungarus
8.Analysis of impact of host plants on quality of Taxilli Herba based on widely targeted metabolomics.
Dong-Lan ZHOU ; Zi-Shu CHAI ; Mei RU ; Fei-Ying HUANG ; Xie-Jun ZHANG ; Min GUO ; Yong-Hua LI
China Journal of Chinese Materia Medica 2025;50(12):3281-3290
This study aims to explore the impact of host plants on the quality of Taxilli Herba and provide a theoretical basis for the quality control of Taxilli Herba. The components of Taxilli Herba from three different host plants(Morus alba, Salix babylonica, and Cinnamomum cassia) and its 3 hosts(mulberry branch, willow branch, and cinnamon branch) were detected by widely targeted metabolomics based on ultra-high performance liquid chromatography-tandem mass spectrometry(UPLC-MS/MS). Principal component analysis(PCA), orthogonal partial least squares discriminant analysis(OPLS-DA), and Venn diagram were employed for analysis. A total of 717 metabolites were detected in Taxilli Herba from the three host plants and the branches of these host plants by UPLC-MS/MS. The results of PCA and OPLS-DA of Taxilli Herba from the three different host plants showed an obvious separation trend due to the different effects of host plants. The Venn diagram showed that there were 32, 8, and 26 characteristic metabolites in samples of Taxilli Herba from M. alba host, S. babylonica host, and C. cassia host, respectively. It was found by comparing the characteristic metabolites of Taxilli Herba and its hosts that each host transmits its characteristic components to Taxilli Herba, so that the Taxilli Herba contains the characteristic components of the host. The Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway analysis showed that the differential metabolites of Taxilli Herba from the three hosts were mainly enriched in flavonoid biosynthesis, arginine and proline metabolism, and glycolysis/gluconeogenesis pathways. Furthermore, the differential metabolites enriching pathways of Taxilli Herba from the three hosts were different depending on the host. In a word, host plants have a significant impact on the metabolites of Taxilli Herba, and it may be an important factor for the quality of Taxilli Herba.
Metabolomics/methods*
;
Drugs, Chinese Herbal/chemistry*
;
Chromatography, High Pressure Liquid
;
Tandem Mass Spectrometry
;
Quality Control
;
Salix/chemistry*
;
Cinnamomum aromaticum/metabolism*
;
Principal Component Analysis
9.Mechanism related to bile acids metabolism of liver injury induced by long-term administration of emodin.
Jing-Zhuo TIAN ; Lian-Mei WANG ; Yan YI ; Zhong XIAN ; Nuo DENG ; Yong ZHAO ; Chun-Ying LI ; Yu-Shi ZHANG ; Su-Yan LIU ; Jia-Yin HAN ; Chen PAN ; Chen-Yue LIU ; Jing MENG ; Ai-Hua LIANG
China Journal of Chinese Materia Medica 2025;50(11):3079-3087
Emodin is a hydroxyanthraquinone compound that is widely distributed and has multiple pharmacological activities, including anti-diarrheal, anti-inflammatory, and liver-protective effects. Research indicates that emodin may be one of the main components responsible for inducing hepatotoxicity. However, studies on the mechanisms of liver injury are relatively limited, particularly those related to bile acids(BAs) metabolism. This study aims to systematically investigate the effects of different dosages of emodin on BAs metabolism, providing a basis for the safe clinical use of traditional Chinese medicine(TCM)containing emodin. First, this study evaluated the safety of repeated administration of different dosages of emodin over a 5-week period, with a particular focus on its impact on the liver. Next, the composition and content of BAs in serum and liver were analyzed. Subsequently, qRT-PCR was used to detect the mRNA expression of nuclear receptors and transporters related to BAs metabolism. The results showed that 1 g·kg~(-1) emodin induced hepatic damage, with bile duct hyperplasia as the primary pathological manifestation. It significantly increased the levels of various BAs in the serum and primary BAs(including taurine-conjugated and free BAs) in the liver. Additionally, it downregulated the mRNA expression of farnesoid X receptor(FXR), retinoid X receptor(RXR), and sodium taurocholate cotransporting polypeptide(NTCP), and upregulated the mRNA expression of cholesterol 7α-hydroxylase(CYP7A1) in the liver. Although 0.01 g·kg~(-1) and 0.03 g·kg~(-1) emodin did not induce obvious liver injury, they significantly increased the level of taurine-conjugated BAs in the liver, suggesting a potential interference with BAs homeostasis. In conclusion, 1 g·kg~(-1) emodin may promote the production of primary BAs in the liver by affecting the FXR-RXR-CYP7A1 pathway, inhibit NTCP expression, and reduce BA reabsorption in the liver, resulting in BA accumulation in the peripheral blood. This disruption of BA homeostasis leads to liver injury. Even doses of emodin close to the clinical dose can also have a certain effect on the homeostasis of BAs. Therefore, when using traditional Chinese medicine or formulas containing emodin in clinical practice, it is necessary to regularly monitor liver function indicators and closely monitor the risk of drug-induced liver injury.
Emodin/administration & dosage*
;
Bile Acids and Salts/metabolism*
;
Animals
;
Male
;
Liver/injuries*
;
Chemical and Drug Induced Liver Injury/genetics*
;
Drugs, Chinese Herbal/adverse effects*
;
Humans
;
Rats, Sprague-Dawley
;
Mice
;
Rats
10.Textual study of Baihuasheshecao (Hedyotis diffusa).
Dong-Min JIANG ; Chu-Chu ZHONG ; Pang-Chui SHAW ; Bik-San LAU ; Tai-Wai LAU ; Guang-Hao XU ; Ying ZHANG ; Zhi-Guo MA ; Hui CAO ; Meng-Hua WU
China Journal of Chinese Materia Medica 2025;50(15):4386-4396
Baihuasheshecao(Hedyotis diffusa) is a commonly used traditional Chinese medicine derived from the whole herb of H. diffusa and has been widely utilized in folk medicine. It possesses anti-tumor, antibacterial, and anti-inflammatory properties, making it one of the frequently used herbs in TCM clinical practice. However, Shuixiancao(H. corymbosa) and Xianhuaercao(H. tenelliflora), species of the same genus, are often used as substitutes for Baihuasheshecao. To substantiate the medicinal basis of Baihuasheshecao, this study systematically reviewed classical herbal texts and modern literature, examining its nomenclature, botanical origin, harvesting, processing, properties, meridian tropism, pharmacological effects, and clinical applications. The results indicate that Baihuasheshecao was initially recorded as "Shuixiancao" in Preface to the Indexes to the Great Chinese Botany(Zhi Wu Ming Shi Tu Kao). Based on its morphological characteristics and habitat description, it was identified as H. diffusa in the Rubiaceae family. Subsequent records predominantly refer to it as Baihuasheshecao as its official name. In most regions, Baihuasheshecao is recognized as the authentic medicinal material, distinct from Shuixiancao and Xianhuaercao. Baihuasheshecao is harvested in late summer and early autumn, and the dried whole plant, including its roots, is used medicinally. The standard processing method involves cutting. It is known for its effects in clearing heat, removing toxins, reducing swelling and pain, and promoting diuresis to resolve abscesses. Initially, it was mainly used for treating appendicitis, intestinal abscesses, and venomous snake bites, and later, it became a treatment for cancer. The excavation of its clinical value followed a process in which overseas Chinese introduced the herb from Chinese folk medicine to other countries. After its unique anti-cancer effects were recognized abroad, it was reintroduced to China and gradually became a crucial TCM for cancer treatment. The findings of this study help clarify the historical and contemporary uses of Baihuasheshecao, providing literature support and a scientific basis for its rational development and precise clinical application.
Humans
;
China
;
Drugs, Chinese Herbal/chemistry*
;
Hedyotis/classification*
;
Medicine, Chinese Traditional/history*

Result Analysis
Print
Save
E-mail