1.SIRT5 Potentiates Hepatocarcinogenesis by Modulating Protein Acylation in Mice
Yu ZHANG ; Feng-Rui REN ; Jia-Yun LI ; Xiang-Yu CHEN ; Zi-Yi WANG ; Qi SUN ; Jun-Cheng ZHAO ; Ye ZHANG ; Zhen HUANG ; Hao HU ; Tao-Tao WEI ; Min XIAO
Progress in Biochemistry and Biophysics 2026;53(6):1712-1722
ObjectiveHepatocellular carcinoma (HCC) represents 90% of all primary liver cancers. The main risk factors associated with HCC include viral hepatitis (B and/or C), alcohol abuse, and metabolic dysfunction-associated steatotic liver disease (MASLD), which progressively advance to liver fibrosis, cirrhosis, and ultimately evolve into HCC. Surgical resection represents the most effective treatment for HCC, while recent advances in immunotherapy, including immune checkpoint inhibitors and adoptive cell therapies, have provided improved treatment prospects for patients with unresectable HCC. However, the complex metabolic heterogeneity of HCC limits the therapeutic efficacy. Metabolic intermediates acyl-CoA not only provide energy and substrates for numerous biochemical reactions but also serve as donors for protein lysine acylation, a major class of post-translational modification (PTM). Therefore, a deeper understanding of the molecular mechanisms underlying protein lysine acylation and hepatocarcinogenesis is urgently needed. MethodsThe levels of protein lysine acylation and silence information regulator 5 (SIRT5) expression levels in clinical HCC samples were analyzed by Western blot. Quantitative malonylome and succinylome of HCC samples were analyzed by antibody-based affinity enrichment coupled with tandem mass spectrometry. The proliferation of HCC cells was analyzed with Cell Counting Kit-8 (CCK-8) assays, the apoptosis was quantified by Annexin V-FITC/propidium iodide (PI) staining coupled with flow cytometry, and the ability of cells to migrate was assayed by Transwell assays. The enzymatic activity of glutathione S-transferase Mu 1 (GSTM1) was quantified. Transgenic mice with hepatic overexpression of SIRT5 were constructed using CRISPR-Cas9, and primary hepatocarcinogenesis was induced by administration of diethylnitrosamine. ResultsWestern blot analysis indicated that the expression level of SIRT5 was elevated in clinical samples from HCC patients, and the levels of lysine malonylation, glutarylation, and succinylation were significantly reduced in HCC tissues. Knockout of SIRT5 in MHCC-97H and MHCC-97L hepatoma cells suppressed cell proliferation, and increased the percentage of apoptotic cells significantly. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the differentially malonylome and succinylome of HCC samples revealed significant enrichment in two major classes of biological processes: core energy metabolism (e.g., glycolysis/gluconeogenesis, tricarboxylic acid metabolic process, fatty acid beta oxidation) and detoxification and oxidative stress response (e.g., response to toxic substance, chemical carcinogenesis, reactive oxygen species (ROS)). SIRT5 removes malonylation from lysine residues in GSTM1 and restores its detoxification activity, which is crucial for the survival of hepatocytes under stressed conditions. More importantly, in vivo experiment indicated that hepatic-specific overexpression of SIRT5 in mice accelerated diethylnitrosamine-induced liver fibrosis and hepatocarcinogenesis, indicating the critical role of SIRT5 in HCC progression. ConclusionThis study highlights the previously unrecognized SIRT5-GSTM1 axis as a key regulator in hepatocarcinogenesis, and suggests a potential target for the treatment of patients with HCC.
2.Research on Electrical Impedance and Microwave Dual-modality Tomography Algorithm Based on Conditional Diffusion Models
Jin-Zhen LIU ; Xiang-Qian MENG ; Hui XIONG ; Li-Min ZHOU ; Chun-Chan LI
Progress in Biochemistry and Biophysics 2026;53(6):1780-1792
ObjectiveStroke poses a heavy burden due to its high mortality and morbidity rates. Accurate and real-time detection of lesions is pivotal for prompt clinical intervention and favorable prognosis. Electrical impedance tomography (EIT) and microwave tomography (MWT) have emerged as compelling alternatives for stroke screening, owing to their non-ionizing, non-invasive and portable nature. EIT provides information on tissue conductivity, and MWT offers high sensitivity to changes in dielectric properties. However, single-modality imaging is inherently limited, EIT suffers from low sensitivity to deep-seated tissues and severe ill-posedness of inverse problems, whereas MWT is challenged by strong nonlinearity in inverse scattering and susceptibility to modeling errors. Consequently, the clinical utility of standalone EIT or MWT for stroke diagnosis remains constrained by poor spatial resolution and imaging artifacts. To improve the accuracy and robustness of stroke imaging, a dual-modality fusion conditional denoising diffusion probabilistic model (DM-DDPM) was proposed for high-precision dual-modality image reconstruction. MethodsA dual-encoder network with a symmetric architecture and independently trained parameters was constructed to extract heterogeneous features separately from EIT boundary voltage measurements and MWT scattered field signals. Attentional feature fusion (AFF) is employed to integrate complementary information from the two modalities adaptively, generating robust fused priors that suppress redundant noise while preserving key physical characteristics. Subsequently, the fused priors are embedded into a Transformer-based diffusion model via a cross attention mechanism to guide the reverse denoising process. This approach effectively reduces artifacts and enhances the stability of conductivity distribution reconstruction. Time step embedding is introduced to enable the network to perceive the diffusion stage and further improve the accuracy of noise prediction. ResultsSimulated experiments demonstrated that DM-DDPM significantly outperforms single-modality and multi-modality networks under various noise levels. A head model simulation dataset was constructed based on COMSOL Multiphysics, and tests were carried out under 50 dB, 40 dB and 30 dB signal-to-noise ratio levels. At 30 dB, the average relative error (RE) was below 0.20, while the structural similarity index measure (SSIM) and correlation coefficient (CC) remained above 0.90 and 0.89, respectively. Compared with single-modality and multi-modality networks, artifacts were significantly reduced, lesion edges were clearer, and localization was more accurate. The model maintains high reconstruction quality and strong robustness for single, double, and triple lesions simultaneously. Furthermore, physical experiments were conducted using a 16-electrode EIT system and a 16-antenna MWT system with asynchronous data acquisition. These experiments confirmed the feasibility of the method in real-world scenarios and demonstrated that it can robustly reconstruct simulated lesions despite environmental interference and measurement noise, validating its reliability for practical clinical applications. ConclusionThe proposed method effectively combines complementary dual-modality information with a conditional diffusion model. Low accuracy and poor noise resistance in single-modality imaging were effectively addressed, while the noise amplification issue caused by direct multimodal data fusion was avoided. The proposed algorithm exhibits strong anti-noise interference ability and high imaging stability in both simulation and physical experiments. Precise localization of stroke lesions with different quantities was achieved, providing a high-precision, and practical technical support for clinical stroke detection.
3.Resistance of Aedes albopictus to common insecticides in Zhejiang Province, 2024
Qin-mei LIU ; Jin-na WANG ; Tian-qi LI ; Ming-yu LUO ; Zhou GUAN ; Zhen-yu GONG ; Ji-min SUN
Acta Parasitologica et Medica Entomologica Sinica 2026;33(2):96-100
Objective In order to understand the current resistance status of Aedes albopictus(Ae. albopictus) adult population to commonly used insecticides and the resistance dynamics development in Zhejiang Province and provide a theoretical basis for the scientific and rational use of hygienic insecticides. Methods Resistance determination was carried out using the diagnostic dose method(GB/T 26347-2010)recommended by WHO in the adult mosquito contact cylinder method to monitor the resistance of Ae. albopictus field populations in 11 cities and districts in Zhejiang Province in 2024. Results The monitoring result showed that adult Ae. albopictus in 11 areas were resistant or suspected resistant populations to both cis-cypermethrin and permethrin; suspected and resistant populations to the organophosphate insecticide malathion were as high as 81.82%; while suspected resistance to the carbamate insecticide residual carbofuran accounted for 45.45%, and the rest of the monitoring sites were sensitive to it. Conclusions The result suggest that the generalized resistance of Ae. albopictus populations to pyrethroids and organophosphorus insecticides in the field in Zhejiang Province needs to be highly emphasized. In future mosquito vector and mosquito-borne disease control, insecticides should be selected based on the resistance monitoring information, while following scientific guidelines for insecticide use.
4.Diagnosis and treatment of colorectal liver metastases: Chinese expert consensus-based multidisciplinary team (2024 edition).
Wen ZHANG ; Xinyu BI ; Yongkun SUN ; Yuan TANG ; Haizhen LU ; Jun JIANG ; Haitao ZHOU ; Yue HAN ; Min YANG ; Xiao CHEN ; Zhen HUANG ; Weihua LI ; Zhiyu LI ; Yufei LU ; Kun WANG ; Xiaobo YANG ; Jianguo ZHOU ; Wenyu ZHANG ; Muxing LI ; Yefan ZHANG ; Jianjun ZHAO ; Aiping ZHOU ; Jianqiang CAI
Chinese Medical Journal 2025;138(15):1765-1768
5.The role of microglia activated by the deletion of immune checkpoint receptor CD200R1 gene in a mouse model of Parkinson's disease.
Jia-Li GUO ; Tao-Ying HUANG ; Zhen ZHANG ; Kun NIU ; Xarbat GONGBIKAI ; Xiao-Li GONG ; Xiao-Min WANG ; Ting ZHANG
Acta Physiologica Sinica 2025;77(1):13-24
The study aimed to investigate the effect of the CD200R1 gene deletion on microglia activation and nigrostriatal dopamine neuron loss in the Parkinson's disease (PD) process. The CRISPR-Cas9 technology was applied to construct the CD200R1-/- mice. The primary microglia cells of wild-type and CD200R1-/- mice were cultured and treated with bacterial lipopolysaccharide (LPS). Microglia phagocytosis level was assessed by a fluorescent microsphere phagocytosis assay. PD mouse model was prepared by nigral stereotaxic injection of recombinant adeno-associated virus vector carrying human α-synuclein (α-syn). The changes in the motor behavior of the mice with both genotypes were evaluated by cylinder test, open field test, and rotarod test. Immunohistochemical staining was used to assess the loss of dopamine neurons in substantia nigra. Immunofluorescence staining was used to detect the expression level of CD68 (a key molecule involved in phagocytosis) in microglia. The results showed that CD200R1 deletion markedly enhanced LPS-induced phagocytosis in vitro by the microglial cells. In the mouse model of PD, CD200R1 deletion exacerbated motor behavior impairment and dopamine neuron loss in substantia nigra. Fluorescence intensity analysis results revealed a significant increase in CD68 expression in microglia located in the substantia nigra of CD200R1-/- mice. The above results suggest that CD200R1 deletion may further activates microglia by promoting microglial phagocytosis, leading to increased loss of the nigrostriatal dopamine neurons in the PD model mice. Therefore, targeting CD200R1 could potentially serve as a novel therapeutic target for the treatment of early-stage PD.
Animals
;
Microglia/physiology*
;
Mice
;
Phagocytosis
;
Parkinson Disease/genetics*
;
Disease Models, Animal
;
Receptors, Cell Surface/physiology*
;
Dopaminergic Neurons/pathology*
;
Antigens, CD/metabolism*
;
Gene Deletion
;
Substantia Nigra
;
Mice, Inbred C57BL
;
Mice, Knockout
;
Cells, Cultured
;
Male
;
alpha-Synuclein
;
CD68 Molecule
;
Orexin Receptors
6.Transcriptomic analysis of key genes involved in sex differences in intellectual development.
Jia-Wei ZHANG ; Xiao-Li ZHENG ; Hai-Qian ZHOU ; Zhen ZHU ; Wei HAN ; Dong-Min YIN
Acta Physiologica Sinica 2025;77(2):211-221
Intelligence encompasses various abilities, including logical reasoning, comprehension, self-awareness, learning, planning, creativity, and problem-solving. Extensive research and practical experience suggest that there are sex differences in intellectual development, with females typically maturing earlier than males. However, the key genes and molecular network mechanisms underlying these sex differences in intellectual development remain unclear. To date, Genome-Wide Association Studies (GWAS) have identified 507 genes that are significantly associated with intelligence. This study first analyzed RNA sequencing data from different stages of brain development (from BrainSpan), revealing that during the late embryonic stage, the average expression levels of intelligence-related genes are higher in males than in females, while the opposite is observed during puberty. This study further constructed interaction networks of intelligence-related genes with sex-differential expression in the brain, including the prenatal male network (HELP-M: intelligence genes with higher expression levels in prenatal males) and the pubertal female network (HELP-F: intelligence genes with higher expression levels in pubertal females). The findings indicate that the key genes in both networks are Ep300 and Ctnnb1. Specifically, Ep300 regulates the transcription of 53 genes in both HELP-M and HELP-F, while Ctnnb1 regulates the transcription of 45 genes. Ctnnb1 plays a more prominent role in HELP-M, while Ep300 is more crucial in HELP-F. Finally, this study conducted sequencing validation on rats at different developmental stages, and the results indicated that in the prefrontal cortex of female rats during adolescence, the expression levels of the intelligence genes in HELP-F, as well as key genes Ep300 and Ctnnb1, were higher than those in male rats. These genes were also involved in neurodevelopment-related biological processes. The findings reveal a sex-differentiated intelligence gene network and its key genes, which exhibit varying expression levels during the neurodevelopmental process.
Female
;
Intelligence/physiology*
;
Male
;
Sex Characteristics
;
Animals
;
Brain/growth & development*
;
E1A-Associated p300 Protein/physiology*
;
beta Catenin/physiology*
;
Transcriptome
;
Rats
;
Gene Expression Profiling
;
Genome-Wide Association Study
7.c-Met-targeted chimeric antigen receptor T cells inhibit human serous ovarian cancer cell SKOV-3 in vitro.
Na-Na DU ; Yan-Jun ZHANG ; Yan-Qiu LI ; Lu ZHANG ; Ran AN ; Xiang-Cheng ZHEN ; Jing-Ting MIN ; Zheng-Hong LI
Acta Physiologica Sinica 2025;77(2):241-254
The study aimed to construct the second and third generation chimeric antigen receptor T cells (CAR-T) targeting the c-mesenchymal-epithelial transition factor (c-Met) protein, and observe their killing effect on human serous ovarian cancer cell SKOV-3. The expression of MET gene in ovarian serous cystadenocarcinoma, the correlation between MET gene expression and the abundance of immune cell infiltration, and the effect of MET gene expression on the tissue function of ovarian serous cystadenocarcinoma were analyzed by bioinformatics. The expression of c-Met in ovarian cancer tissues and adjacent tissues was detected by immunohistochemical staining. The second and third generation c-Met CAR-T cells, namely c-Met CAR-T(2G/3G), were prepared by lentivirus infection, and the cell subsets and infection efficiency were detected by flow cytometry. Using CD19 CAR-T and activated T cells as control groups and A2780 cells with c-Met negative expression as Non target groups, the kill efficiency on SKOV-3 cells with c-Met positive expression, cytokine release and cell proliferation of c-Met CAR-T(2G/3G) were explored by lactate dehydrogenase (LDH) release, ELISA and CCK-8 respectively. The results showed that MET gene expression was significantly up-regulated in ovarian cancer tissues compared with normal tissues, which was consistent with the immunohistochemistry results. However, in all pathological stages, there was no obvious difference in MET expression and no correlation between MET gene expression and the race and age of ovarian cancer patients. The second generation and third generation c-Met CAR-T cells were successfully constructed. After lentivirus infection, the proportion of CD8+ T cells in c-Met CAR-T(2G) was upregulated, while there was no significant change in the cell subsets of c-Met CAR-T(3G). The LDH release experiment showed that the kill efficiency of c-Met CAR-T(2G/3G) on SKOV-3 increased with the increase of effect-target ratio. When the effect-target ratio was 20:1, the kill efficiency of c-Met CAR-T(2G) reached (42.02 ± 5.17)% (P < 0.05), and the kill efficiency of c-Met CAR-T(3G) reached (51.40 ± 2.71)% (P < 0.05). ELISA results showed that c-Met CAR-T released more cytokine compared to CD19 CAR-T and activated T cells (P < 0.05). Moreover, the cytokine release of c-Met CAR-T(3G) was higher than c-Met CAR-T(2G) (P < 0.01). The CCK-8 results showed that after 48 h, the cell number of c-Met CAR-T(2G) was higher than that of c-Met CAR-T(3G) (P < 0.01). In conclusion, both the second and third generation c-Met CAR-T can target and kill c-Met-positive SKOV-3 cells, with no significant difference. c-Met CAR-T(2G) has stronger proliferative ability, and c-Met CAR-T(3G) releases more cytokines.
Humans
;
Female
;
Ovarian Neoplasms/immunology*
;
Proto-Oncogene Proteins c-met/metabolism*
;
Receptors, Chimeric Antigen/immunology*
;
Cell Line, Tumor
;
Cystadenocarcinoma, Serous/immunology*
;
T-Lymphocytes/immunology*
8.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
9.Biomarkers of hepatotoxicity in rats induced by aqueous extract of Dictamni Cortex based on urine metabolomics.
Hui-Juan SUN ; Rui GAO ; Meng-Meng ZHANG ; Ge-Yu DENG ; Lin HUANG ; Zhen-Dong ZHANG ; Yu WANG ; Fang LU ; Shu-Min LIU
China Journal of Chinese Materia Medica 2025;50(9):2526-2538
This paper aimed to use non-targeted urine metabolomics to reveal the potential biomarkers of toxicity in rats with hepatic injury induced by aqueous extracts of Dictamni Cortex(ADC). Forty-eight SD rats were randomly assigned to a blank group and high-dose, medium-dose, and low-dose ADC groups, with 12 rats in each group(half male and half female), and they were administered orally for four weeks. The hepatic injury in SD rats was assessed by body weight, liver weight/index, biochemical index, L-glutathione(GSH), malondialdehyde(MDA), and pathological alterations. The qPCR was utilized to determine the expression of metabolic enzymes in the liver and inflammatory factors. Differential metabolites were screened using principal component analysis(PCA) and partial least squares-discriminant analysis(PLS-DA), followed by a metabolic pathway analysis. The Mantel test was performed to assess differential metabolites and abnormally expressed biochemical indexes, obtaining potential biomarkers. The high-dose ADC group showed a decrease in body weight and an increase in liver weight and index, resulting in hepatic inflammatory cell infiltration and hepatic steatosis. In addition, this group showed elevated levels of MDA, cytochrome P450(CYP) 3A1, interleukin-1β(IL-1β), and tumor necrosis factor-α(TNF-α), as well as lower levels of alanine transaminase(ALT) and GSH. A total of 76 differential metabolites were screened from the blank and high-dose ADC groups, which were mainly involved in the pentose phosphate pathway, tryptophan metabolism, purine metabolism, pentose and glucuronic acid interconversion, galactose metabolism, glutathione metabolism, and other pathways. The Mantel test identified biomarkers of hepatotoxicity induced by ADC in SD rats, including glycineamideribotide, dIDP, and galactosylglycerol. In summary, ADC induced hepatotoxicity by disrupting glucose metabolism, ferroptosis, purine metabolism, and other pathways in rats, and glycineamideribotide, dIDP, and galactosylglycerol could be employed as the biomarkers of its toxicity.
Animals
;
Male
;
Rats, Sprague-Dawley
;
Rats
;
Metabolomics
;
Biomarkers/metabolism*
;
Liver/metabolism*
;
Drugs, Chinese Herbal/adverse effects*
;
Female
;
Chemical and Drug Induced Liver Injury/metabolism*
;
Glutathione/metabolism*
;
Humans
10.Original plants, production areas, and spread of Sanqi: based on historical materials of the Ming and Qing Dynasties.
Mei WANG ; Min-Zhen YIN ; Hua-Sheng PENG
China Journal of Chinese Materia Medica 2025;50(11):2938-2944
Sanqi is first recorded in the Compendium of Materia Medica(Ben Cao Gang Mu) in the Ming Dynasty. During the Ming and Qing Dynasties, Sanqi, as a precious Dao-di herb, was successively spread and introduced for cultivation. This study verified the germplasm resources, production areas, and spread of Sanqi in the Ming and Qing Dynasties by systematically reviewing the historical materials, such as materia medica works and local chronicles, and the modern distribution of production areas. In the Ming and Qing Dynasties, the original plants of Sanqi included Panax notoginseng, P. japonicus, P. bipinnatifidus, P. zingiberensis, P. stipuleanatus, and Gynura japonica. Among them, the production area of P. notoginseng has changed. From 1578 to 1593, the main production areas of P. notoginseng were Nandan county, Hechi city in Guangxi Zhuang autonomous region and Guangnan county and Funing county, Wenshan prefecture in Yunnan province. From 1683 to 1755, the production areas of P. notoginseng additionally included Yizhou district, Tian'e county, and Huanjiang county in Hechi city, and Tianyang district and Tiandong county in Baise city, Xincheng county and Gongcheng county in Guangxi Zhuang autonomous region. From 1765 to 1892, the production areas additionally included Youjiang district, Debao county, Napo county, and Jingxi city in Baise city, and Tiandeng county in Guangxi Zhuang autonomous region, and Wenshan city, Malipo county, Yanshan county, Xichou county, and Maguan county in Wenshan prefecture, and Baoshan city, Dali prefecture, Lincang city, Honghe prefecture, Mangshi city, and Lushui city in Yunnan province. During the Wanli period of the Ming Dynasty, Sanqi was introduced to Zhejiang province. During the Qianlong period of the Qing Dynasty, it was introduced to Fujian province. During the Daoguang period of the Qing Dynasty, it was introduced to Hunan province. By comprehensively reviewing the materia medica works, local chronicles, and novel historical materials, this study restores the development history of the Sanqi industry in the Ming and Qing Dynasties. Historical data show that the introduction of Dao-di herbs should consider the biological characteristics of medicinal plants and avoid blind introduction.
China
;
Drugs, Chinese Herbal/history*
;
History, 17th Century
;
History, 16th Century
;
Plants, Medicinal/chemistry*
;
Medicine, Chinese Traditional/history*
;
History, 18th Century


Result Analysis
Print
Save
E-mail