1.Development of an Analytical Software for Forensic Proteomic SAP Typing
Feng HU ; Meng-Jiao WANG ; Jia-Lei WU ; Dong-Sheng DING ; Zhi-Yuan YANG ; An-Quan JI ; Lei FENG ; Jian YE
Progress in Biochemistry and Biophysics 2025;52(9):2406-2416
ObjectiveThe proteome of biological evidence contains rich genetic information, namely single amino acid polymorphisms (SAPs) in protein sequences. However, due to the lack of efficient and convenient analysis tools, the application of SAP in public security still faces many challenges. This paper aims to meet the application requirements of SAP analysis for forensic biological evidence’s proteome data. MethodsThe software is divided into three modules. First, based on a built-in database of common non-synonymous single nucleotide polymorphisms (nsSNPs) and SAPs in East Asian populations, the software integrates and annotates newly identified exonic nsSNPs as SAPs, thereby constructing a customized SAP protein sequence database. It then utilizes a pre-installed search engine—either pFind or MaxQuant—to perform analysis and output SAP typing results, identifying both reference and variant types, along with their corresponding imputed nsSNPs. Finally, SAPTyper compares the proteome-based typing results with the individual’s exome-derived nsSNP profile and outputs the comparison report. ResultsSAPTyper accepts proteomic DDA mass spectrometry raw data (DDA acquisition mode) and exome sequencing results of nsSNPs as input and outputs the report of SAPs result. The pFind and Maxquant search engines were used to test the proteome data of 2 hair shafts of2 individuals, and both obtained SAP results. It was found that the results of the Maxquant search engine were slightly less than those of pFind. This result shows that SAPTyper can achieve SAP fingding function. Moreover, the pFind search engine was used to test the proteome data of 3 hair shafts from 1 European person and 1 African person in the literature. Among the sites fully matched by the literature method, sites detected by SAPTyper are also included; for semi-matching sites, that is, nsSNPs are heterozygous, both literature method and SAPTyper method had the risk of missing detection for one type of the allele. Comparing the analysis results of SAPTyper with the SAP test results reported in the literature, it was found that some imputed nsSNP sites identified by the literature method but not detected by SAPTyper had a MAF of less than 0.1% in East Asian populations, and therefore they were not included in the common nsSNP database of East Asian populations constructed by this software. Since the database construction of this software is based on the genetic variation information of East Asian populations, it is currently unable to effectively identify representative unique common variation sites in European or African populations, but it can still identify SAP sites shared by these populations and East Asian populations. ConclusionAn automated SAP analysis algorithm was developed for East Asian populations, and the software named SAPTyper was developed. This software provides a convenient and efficient analysis tool for the research and application of forensic proteomic SAP and has important application prospects in individual identification and phenotypic inference based on SAP.
2.Identification of Brucella strains in Ordos City, Inner Mongolia based on whole genome sequencing
WANG Nana ; JIA Yuan ; WANG Sheng
China Tropical Medicine 2025;25(3):333-
Objective To explore the whole genome structure of Brucella in Ordos city, conduct bioinformatics analysis, and explore a variety of Brucella species type identification methods, so as to provide a theoretical basis for the study of Brucella in Ordos City. Methods The whole genome was sequenced using illumina Miseq gene sequencer. The Microobench microbial gene analysis workstation was used to analyze the genome components, genome annotation, 16Sr DNA, SNP sites, insertions and deletions, MLST sites, and phylogenetic tree of 15 strains of Brucella. Results The average sequencing depth of 15 Brucella strains was 170.04×, the average genome size was 3 290 725 bp, the genome coverage was 99.93%, and the average GC content of assembly files was 54.69%. The results of 16Sr DNA analysis showed that the homology of 15 Brucella strains and Brucella Maltese was 100%. The analysis results of SNP sites and insert-deletion variation showed that the average number of SNPs of the 15 strains of Brucella was 2 412, the average number of polynucleotide mutations was 39, the average number of insertional variants was 131, the average number of missing variants was 153, the average number of complex variants was 43, and the average number of total variants was 2 778. MLST locus analysis showed that all 15 strains of Brucella were ST-8. Conclusion Based on the whole genome sequencing of 15 strains of Brucella in Ordos City, this study carried out a variety of type identification methods. Different methods have different focuses, and multiple methods are needed to confirm each other, so as to obtain more accurate experimental results and provide more research basis for the study of Brucella, to provide more data and theoretical support for the control and prevention of brucellosis in Ordos City.
3.Effects of Huayu Tongluo moxibustion on cognitive function and insulin resistance in patients with type 2 diabetes mellitus and cognitive decline: a randomized controlled trial.
Min YE ; Aihong YUAN ; Lele ZHANG ; Qiqi YANG ; Hongyu XIE ; Xia GE ; Wenjing KAN ; Sheng LI ; Jun YANG
Chinese Acupuncture & Moxibustion 2025;45(11):1541-1548
OBJECTIVE:
To investigate the effects of Huayu Tongluo (transforming stasis and unblocking collaterals) moxibustion on cognitive function and insulin resistance in patients with type 2 diabetes mellitus (T2DM) and cognitive decline.
METHODS:
Ninety patients with T2DM and cognitive decline were randomly divided into a moxibustion group (n=45, 3 cases dropped out, 2 cases were eliminated) and a waiting moxibustion group (n=45, 2 cases dropped out). Both groups received routine hypoglycemic treatment for 12 weeks. The moxibustion group additionally received Huayu Tongluo moxibustion at Baihui (GV20), Shenting (GV24), and Dazhui (GV14). Pressing moxibustion was applied to Baihui (GV20) for 20 min, while suspended moxibustion was applied to Shenting (GV24) and Dazhui (GV14) for 20 min each. Treatments of moxibustion were administered every other day (three times per week) for 12 weeks. All patients were followed up for 12 weeks, during which their original hypoglycemic medication regimen was maintained. Before treatment, after 12 weeks of treatment, and at the 12-week follow-up, the scores of Montreal cognitive assessment (MoCA), mini-mental state examination (MMSE), Addenbrooke's cognitive examination Ⅲ (ACE-Ⅲ), symbol digit modalities test (SDMT), and Athens insomnia scale (AIS) and the insulin resistance index (HOMA-IR) were observed in the two groups.
RESULTS:
Compared with before treatment, the MoCA scores, MMSE scores, ACE-Ⅲ subscale scores (attention, memory, language fluency, language, visuospatial ability) and total scores, and SDMT scores were increased (P<0.01), while the AIS scores were decreased (P<0.05) in the moxibustion group after treatment and at follow-up. Compared with before treatment, the MMSE score, ACE-Ⅲ subscale scores (memory, attention) and total score after treatment, as well as the ACE-Ⅲ subscale scores (language, memory, attention) and total score, and SDMT score at follow-up were increased (P<0.05, P<0.01) in the waiting moxibustion group. Compared with before treatment, HOMA-IR was decreased in both groups after treatment and at follow-up (P<0.01). At follow-up, ACE-Ⅲ subscale scores (attention, memory), and the total score in the moxibustion group were lower than those after treatment (P<0.05, P<0.01), and the ACE-Ⅲ language subscale score, total ACE-Ⅲ score, and SDMT score in the waiting moxibustion group were higher than those after treatment (P<0.01, P<0.05). After treatment and at follow-up, compared with the waiting moxibustion group, the moxibustion group had higher MoCA scores, MMSE scores, SDMT scores, ACE-Ⅲ subscale scores (attention, memory, language fluency) and total scores (P<0.05, P<0.01), and lower HOMA-IR (P<0.05).
CONCLUSION
Huayu Tongluo moxibustion can effectively improve cognitive function in patients with T2DM and cognitive decline. This improvement may be associated with the reduction in insulin resistance.
Humans
;
Insulin Resistance
;
Diabetes Mellitus, Type 2/complications*
;
Male
;
Female
;
Moxibustion
;
Middle Aged
;
Aged
;
Cognition
;
Acupuncture Points
;
Adult
;
Cognitive Dysfunction/therapy*
4.Chemical composition and efficacy of warming lung and resolving fluid retention of Asarum forbesii grown under different shading conditions.
Lu LIAO ; Li-Xian LU ; Hong-Zhuan SHI ; Qiao-Sheng GUO ; Cheng-Hao FEI ; Kun ZHAO ; Yuan-Yuan XING ; Yong SU ; Chang LIU ; Xin-Yue YUAN
China Journal of Chinese Materia Medica 2025;50(2):384-394
Asarum forbesii is a perennial herb born in a shaded and humid environment, which is warm in nature. With the efficacy of warming lung, resolving fluid retention, and relieving coughs, it can be used to treat the syndrome of cold fluid accumulating in lung. To investigate the effects of different shading conditions on the composition and efficacy of A. forbesii, this study planted A. forbesii under 20% natural light(NL20), 40% natural light(NL40), 60% natural light(NL60), and 80% natural light(NL80) and utilized ultra performance liquid chromatography(UPLC) and micro broth 2-fold dilution method to detect the volatile chemical compounds and the minimum inhibitory concentration. At the same time, the study investigated the effects of A. forbesii grown under different shading conditions on the signs, pathological changes of lung tissues, serum cytokine levels, activities of mitochondrial respiratory chain complexes Ⅰ-Ⅴ in lung tissues, and relative expression of related genes of mice with syndrome of cold fluid accumulating in lung. The results indicated that with the increase of shading, the content of kakuol, methyl eugenol, and asarinin in A. forbesii and the antibacterial effect showed a tendency of increasing first and then decreasing, and the NL40 group was significantly better than the other groups. Under the conditions of NL20 and NL40, A. forbesii significantly alleviated the pathological damage to lung tissues, restored the homeostasis of the lung, and enhanced the energy metabolism level of mice with syndrome of cold fluid accumulating in lung. In addition, A. forbesii planted under the two conditions reduced the content of interleukin-8(IL-8), interleukin-13(IL-13), tumor necrosis factor-α(TNF-α), and mucin 5AC(MUC5AC), increased the levels of interleukin-10(IL-10) and aquaporin 1(AQP1), lowered the expression of MMP9, VEGF, TGF-β, and MAPK3. In conclusion, the therapeutic effect of A. forbesii on the syndrome of cold fluid accumulating in lung was positively correlated with the degree of shading, and the chemical composition and efficacy of warming lung and resolving fluid retention were optimal under the conditions of NL20-NL40. This study can provide reference for the pharmacological research and cultivation of A. forbesii.
Animals
;
Mice
;
Lung/pathology*
;
Drugs, Chinese Herbal/administration & dosage*
;
Male
;
Light
;
Cytokines/genetics*
;
Humans
5.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*
6.Research progress in effect of traditional Chinese medicine on aerobic glycolysis in colorectal cancer.
Xu MA ; Sheng-Long LI ; Guang-Rong ZHENG ; Da-Cheng TIAN ; Gang-Gang LU ; Jie GAO ; Yu-Qi AN ; Li-Yuan CAO ; Liang LI ; Xiao-Yong TANG
China Journal of Chinese Materia Medica 2025;50(6):1496-1506
Colorectal cancer(CRC) is a common malignant tumor worldwide. Due to the treatment intolerance and side effects, CRC rank the top among various cancers regarding the incidence and mortality rates. Therefore, exploring new therapies is of great significance for the treatment of CRC. Aerobic glycolysis(AEG) plays an important role in the microenvironment formation, proliferation, metastasis, and recurrence of CRC and other tumor cells. It has been confirmed that intervening in the AEG pathway can effectively curb CRC. The active ingredients and compound prescriptions of traditional Chinese medicine(TCM) can effectively inhibit the proliferation, metastasis, and drug resistance and regulate the apoptosis of tumor cells by modulating AEG-associated transport proteins [eg, glucose transporters(GLUT)], key enzymes [hexokinase(HK) and phosphofructokinase(PFK)], key genes [hypoxia-inducible factor 1(HIF-1) and oncogene(c-Myc)], and signaling pathways(MET/PI3K/Akt/mTOR). Accordingly, they can treat CRC, reduce the recurrence, and improve the prognosis of CRC. Although AEG plays a key role in the development and progression of CRC, the specific mechanisms are not yet fully understood. Therefore, this article delves into the intrinsic connection of the targets and mechanisms of the AEG pathway with CRC from the perspective of tumor cell glycolysis and explores how active ingredients(oxymatrine, kaempferol, and dioscin) and compound prescriptions(Quxie Capsules, Jiedu Sangen Decoction, and Xianlian Jiedu Prescription) of TCM treat CRC by intervening in the AEG pathway. Additionally, this article explores the shortcomings in the current research, aiming to provide reliable targets and a theoretical basis for treating CRC with TCM.
Humans
;
Colorectal Neoplasms/genetics*
;
Drugs, Chinese Herbal/therapeutic use*
;
Glycolysis/drug effects*
;
Animals
;
Medicine, Chinese Traditional
;
Signal Transduction/drug effects*
7.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
8.Effects of total extract of Anthriscus sylvestris on immune inflammation and thrombosis in rats with pulmonary arterial hypertension based on TGF-β1/Smad3 signaling pathway.
Ya-Juan ZHENG ; Pei-Pei YUAN ; Zhen-Kai ZHANG ; Yan-Ling LIU ; Sai-Fei LI ; Yuan RUAN ; Yi CHEN ; Yang FU ; Wei-Sheng FENG ; Xiao-Ke ZHENG
China Journal of Chinese Materia Medica 2025;50(9):2472-2483
This study aimed to explore the effects and mechanisms of total extracts from Anthriscus sylvestris on pulmonary hypertension in rats. Sixty male SD rats were divided into normal(NC) group, model(M) group, positive drug sildenafil(Y) group, low-dose A. sylvestris(ES-L) group, medium-dose A. sylvestris(ES-M) group, and high-dose A. sylvestris(ES-H) group. On day 1, rats were intraperitoneally injected with monocrotaline(60 mg·kg~(-1)) to induce pulmonary hypertension, and the rat model was established on day 28. From days 15 to 28, intragastric administration of the respective treatments was performed. After modeling and treatment, small animal echocardiography was used to detect the right heart function of the rats. Arterial blood gas was measured using a blood gas analyzer. Hematoxylin and eosin(HE) staining and Masson staining were performed to observe cardiopulmonary pathological damage. Flow cytometry was used to detect apoptosis in the lung and myocardial tissues and reactive oxygen species(ROS) levels. Western blot was applied to detect the expression levels of transforming growth factor-β1(TGF-β1), phosphorylated mothers against decapentaplegic homolog 3(p-Smad3), Smad3, tissue plasminogen activator(t-PA), and plasminogen activator inhibitor-1(PAI-1) in lung tissue. A blood routine analyzer was used to measure inflammatory immune cell levels in the blood. Enzyme-linked immunosorbent assay(ELISA) was used to detect the expression levels of P-selectin and thromboxane A2(TXA2) in plasma. The results showed that, compared with the NC group, right heart hypertrophy index, right ventricular free wall thickness, right heart internal diameter, partial carbon dioxide pressure(PaCO_2), apoptosis in cardiopulmonary tissue, and ROS levels were significantly increased in the M group. In contrast, the ratio of pulmonary blood flow acceleration time(PAT)/ejection time(PET), right cardiac output, change rate of right ventricular systolic area, systolic displacement of the tricuspid ring, oxygen partial pressure(PaO_2), and blood oxygen saturation(SaO_2) were significantly decreased in the M group. After administration of the total extract of A. sylvestris, right heart function and blood gas levels were significantly improved, while apoptosis in cardiopulmonary tissue and ROS levels significantly decreased. Further testing revealed that the total extract of A. sylvestris significantly decreased the levels of interleukin-1β(IL-1β), interleukin-6(IL-6), and PAI-1 proteins in lung tissue, while increasing the expression of t-PA. Additionally, the extract reduced the levels of inflammatory cells such as leukocytes, lymphocytes, granulocytes, and monocytes in the blood, as well as the levels of P-selectin and TXA2 in plasma. Metabolomics results showed that the total extract of A. sylvestris significantly affected metabolic pathways, including arginine biosynthesis, tyrosine metabolism, and taurine and hypotaurine metabolism. In conclusion, the total extract of A. sylvestris may exert an anti-pulmonary hypertension effect by inhibiting the TGF-β1/Smad3 signaling pathway, thereby alleviating immune-inflammatory responses and thrombosis.
Animals
;
Male
;
Smad3 Protein/metabolism*
;
Transforming Growth Factor beta1/metabolism*
;
Rats, Sprague-Dawley
;
Rats
;
Signal Transduction/drug effects*
;
Hypertension, Pulmonary/genetics*
;
Thrombosis/immunology*
;
Drugs, Chinese Herbal/administration & dosage*
;
Humans
;
Apoptosis/drug effects*
9.Identification of characteristics, supply channels, and imperial court processing of Arecae Semen in the Qing court.
Feng-Yuan LI ; Hua-Sheng PENG ; Xue-Ling GUAN ; Yan JIN ; Ting YAO ; Yuan YUAN ; Lu-Qi HUANG
China Journal of Chinese Materia Medica 2025;50(11):2924-2930
Qing court records show that Arecae Semen was extensively applied. The royal medical records of the Qing Dynasty document nine types of Arecae Semen, with the Palace Museum preserving seven kinds, totaling twelve cultural relics. Historical documents and physical artifacts corroborate each other, providing evidence for the study of the supply channels and court processing of Arecae Semen in the Qing court. According to relevant Qing court archival records, the sources of Arecae Semen used in the imperial court were diverse, including tributes from foreign countries such as Vietnam and Gurkha, annual tributes from local governments in Guangdong, gifts from close aides, and commodities purchased by the Imperial Household Department from civilian shops. The imperial physicians of the Qing court placed great emphasis on the specifications of Arecae Semen slices and were extremely meticulous about their processing. The variety of Arecae Semen slices used in the Qing palace exceeded those recorded in the botanical texts of the era. Compared with the commonly used processing methods for Arecae Semen in the Qing Dynasty, the imperial physicians adjusted the properties and efficacy of the herbs through different processing techniques, based on the patient's condition, constitution, and other factors, in order to meet the clinical treatment needs of the court. The slicing of Arecae Semen in the Qing court required strict control of thickness, with an average thickness of 0.44 mm, which is significantly thinner than the Arecae Semen slices found in today's markets. The texture was softer, making them easier to chew and absorb. Both the Qing court Arecae Semen slices and the Muxiang Binglang Pills focused on the use of authentic medicinal materials, ensuring the quality of the medicine and enhancing the efficacy of Arecae Semen through meticulous selection and preparation.
China
;
Drugs, Chinese Herbal/history*
;
Humans
;
Medicine, Chinese Traditional/history*
;
History, 19th Century
;
History, Ancient
;
History, 17th Century
;
History, 18th Century
10.Tracing origin of "Qinggong Maidong" production area based on analysis of literature and historical materials and identification of characteristics of Qinggong medicinal materials and cultural relics.
Ao-Yu REN ; Ting YAO ; Feng-Yuan LI ; Hua-Sheng PENG
China Journal of Chinese Materia Medica 2025;50(11):2931-2937
Maidong products are categorized into "Hang Maidong" and "Chuan Maidong". Since the Ming and Qing Dynasties, "Hang Maidong" has been regarded as having superior quality, but currently, it remains in name only in the market. This article reviewed historical materia medica and local chronicles from the Ming and Qing Dynasties and analyzed the historical evolution of Maidong production areas. The history of Maidong production in Zhejiang can be traced back to the Song Dynasty, and cultivation had already developed by at least the Ming Dynasty. During the Ming and Qing Dynasties, it was consistently used as a tribute. Ming Dynasty chronicles record "Chuan Maidong", which had already been cultivated on a large scale by the Qing Dynasty. "Hang Maidong" and "Chuan Maidong" share the same origin, with the former identifiable by the "gourd waist" shape of its tuberous root. Based on this, it can be inferred that the "Maimendong" herb illustrated in the Origins of Materia Medica(Ben Cao Yuan Shi) and the Maidong stored in the Qing Palace Imperial Pharmacy were both "Hang Maidong". The protection and development of the authentic "Hang Maidong" medicinal herb are urgently needed.
China
;
Drugs, Chinese Herbal/history*
;
History, 17th Century
;
History, Ancient
;
Medicine, Chinese Traditional/history*
;
History, Medieval
;
History, 16th Century
;
History, 18th Century
;
History, 15th Century
;
Plants, Medicinal/chemistry*
;
History, 19th Century
;
History, 20th Century
;
Humans
;
Materia Medica/history*
;
History, 21st Century

Result Analysis
Print
Save
E-mail