1.Assessing High-density Y-SNP Panels for Paternal Haplogroup Assignment in Forensic Practice
De-Qin ZHANG ; Chun-Nian WANG ; Lin-Lin LOU ; Meng NI ; Jing GAO ; Jiang HUANG ; Li JIANG
Progress in Biochemistry and Biophysics 2026;53(2):458-469
ObjectiveThe accuracy of Y-chromosome haplogroup assignment is crucial for tracing paternal lineage in male samples. With the advancement of high-throughput sequencing technologies, high-density Y-SNP genotyping from whole-genome or array-based data has become a standard method for determiningY-chromosome haplogroups. This study systematically evaluated the performance of 4 commonly used high-density SNP genotyping systems—namely, the Global Screening Array (GSA), Chinese Genotyping Array (CGA), Affymetrix array, and the 1240K capture panel—for haplogroup assignment. This work provides a reference for data comparison across different systems. MethodsWe extracted genotype data for the 4 Y-SNP panels from 30× whole-genome sequencing (WGS) data of 1 590 male samples from the 1000 Genomes Project. Additionally, GSA array genotype data from 384 relative pairs (spanning 1st- to 12th-degree relationships) from 109 Chinese Han families were collected. Haplogroup assignment was performed using Y-LineageTracker v1.3.0 software. We assessed the concordance and resolution of haplogroup assignments between the four Y-SNP panels and the WGS data. The consistency and resolution of haplogroup assignments were also evaluated for both the 1000 Genomes Project samples and the 109 family samples collected in this study. Furthermore, the impact of varying numbers of Y-SNPs on haplogroup assignment was examined. ResultsThe GSA and CGA panels demonstrated superior resolution and discrimination of haplogroup subclades compared with the other two panels. The haplogroup assignments from the GSA, CGA, and 1240K panels showed high concordance with WGS data, with consistency rates exceeding 88.70%, whereas the Affymetrix platform exhibited a significantly lower consistency rate of 61.89%. Specifically, the GSA and CGA panels consistently demonstrated superior performance compared with the other two panels in the assignment of haplogroups O-M175 and H-L901, achieving complete concordance (100%) for both haplogroups. In contrast, the Affymetrix panel erroneously assigned all individuals belonging to haplogroup O-M175 to haplogroup K2-M526. Furthermore, its accuracy for haplogroup H-L901 was exceedingly low, at merely 1.41%. This poor performance was characterized by the misassignment of 98.59% of H-L901 samples—specifically, 1.41% to J-M304 and a predominant 97.18% to F-M89. For haplogroup R-M207, all four panels exhibited uniformly high levels of consistency, with concordance values exceeding 94.00%. Notably, for haplogroup E-M96, the 1240K and Affymetrix panels outperformed the GSA and CGA panels in terms of concordance, representing the first instance in which these two panels surpassed the latter. Conversely, for haplogroups J-M304, Q-M242, and I-M170, all 4 panels showed relatively elevated misclassification rates, with the Affymetrix array demonstrating the poorest overall performance. None of the four panels showed any discordant haplogroup assignments among the familial relative pairs analyzed. A positive correlation was observed between the number of Y-SNPs (ranging from 1 000 to 10 000) and classification consistency; however, classification consistency plateaued when the number of Y-SNPs exceeded 10 000. Furthermore, a random sampling analysis conducted on the GSA and CGA panels demonstrated that the haplogroup misclassification rate exhibited negligible fluctuation across the Y-SNP range of 500 to 1 000. Conversely, a marked enhancement in classification consistency was observed as the number of markers increased from 1 000 to 5 000, ultimately reaching a plateau within the interval of 5 000 to 8 000 markers. ConclusionThese findings indicate that the GSA and CGA panels provide high resolution and concordance, delivering reliable Y-haplogroup assignment for forensic investigations.
2.Assessing High-density Y-SNP Panels for Paternal Haplogroup Assignment in Forensic Practice
De-Qin ZHANG ; Chun-Nian WANG ; Lin-Lin LOU ; Meng NI ; Jing GAO ; Jiang HUANG ; Li JIANG
Progress in Biochemistry and Biophysics 2026;53(2):458-469
ObjectiveThe accuracy of Y-chromosome haplogroup assignment is crucial for tracing paternal lineage in male samples. With the advancement of high-throughput sequencing technologies, high-density Y-SNP genotyping from whole-genome or array-based data has become a standard method for determiningY-chromosome haplogroups. This study systematically evaluated the performance of 4 commonly used high-density SNP genotyping systems—namely, the Global Screening Array (GSA), Chinese Genotyping Array (CGA), Affymetrix array, and the 1240K capture panel—for haplogroup assignment. This work provides a reference for data comparison across different systems. MethodsWe extracted genotype data for the 4 Y-SNP panels from 30× whole-genome sequencing (WGS) data of 1 590 male samples from the 1000 Genomes Project. Additionally, GSA array genotype data from 384 relative pairs (spanning 1st- to 12th-degree relationships) from 109 Chinese Han families were collected. Haplogroup assignment was performed using Y-LineageTracker v1.3.0 software. We assessed the concordance and resolution of haplogroup assignments between the four Y-SNP panels and the WGS data. The consistency and resolution of haplogroup assignments were also evaluated for both the 1000 Genomes Project samples and the 109 family samples collected in this study. Furthermore, the impact of varying numbers of Y-SNPs on haplogroup assignment was examined. ResultsThe GSA and CGA panels demonstrated superior resolution and discrimination of haplogroup subclades compared with the other two panels. The haplogroup assignments from the GSA, CGA, and 1240K panels showed high concordance with WGS data, with consistency rates exceeding 88.70%, whereas the Affymetrix platform exhibited a significantly lower consistency rate of 61.89%. Specifically, the GSA and CGA panels consistently demonstrated superior performance compared with the other two panels in the assignment of haplogroups O-M175 and H-L901, achieving complete concordance (100%) for both haplogroups. In contrast, the Affymetrix panel erroneously assigned all individuals belonging to haplogroup O-M175 to haplogroup K2-M526. Furthermore, its accuracy for haplogroup H-L901 was exceedingly low, at merely 1.41%. This poor performance was characterized by the misassignment of 98.59% of H-L901 samples—specifically, 1.41% to J-M304 and a predominant 97.18% to F-M89. For haplogroup R-M207, all four panels exhibited uniformly high levels of consistency, with concordance values exceeding 94.00%. Notably, for haplogroup E-M96, the 1240K and Affymetrix panels outperformed the GSA and CGA panels in terms of concordance, representing the first instance in which these two panels surpassed the latter. Conversely, for haplogroups J-M304, Q-M242, and I-M170, all 4 panels showed relatively elevated misclassification rates, with the Affymetrix array demonstrating the poorest overall performance. None of the four panels showed any discordant haplogroup assignments among the familial relative pairs analyzed. A positive correlation was observed between the number of Y-SNPs (ranging from 1 000 to 10 000) and classification consistency; however, classification consistency plateaued when the number of Y-SNPs exceeded 10 000. Furthermore, a random sampling analysis conducted on the GSA and CGA panels demonstrated that the haplogroup misclassification rate exhibited negligible fluctuation across the Y-SNP range of 500 to 1 000. Conversely, a marked enhancement in classification consistency was observed as the number of markers increased from 1 000 to 5 000, ultimately reaching a plateau within the interval of 5 000 to 8 000 markers. ConclusionThese findings indicate that the GSA and CGA panels provide high resolution and concordance, delivering reliable Y-haplogroup assignment for forensic investigations.
3.The SMAD-Pathway Mediates HMGB1-Induced Proliferation and Metastatic Progression in Cutaneous Squamous Cell Carcinoma Cells
De-De LIAN ; Xue Mei LI ; Yu-Xi JIA ; Ming-Wei ZHOU ; Xiang-Ru CHEN ; Yang-Yang TIAN ; Min LI ; Ming-Hui SUN ; Ye ZHAO ; Hong-Jun LI ; Qing-Ling ZHANG
Annals of Dermatology 2026;38(1):51-58
Background:
High-mobility group box protein 1 (HMGB1) is a chromatin-binding protein involved in arthritis, ischemia, sepsis, atherosclerosis, neurodegenerative disorders, meningitis, and cancer. HMGB1 exhibits dual roles in cancer, acting as either a tumor suppressor or oncoprotein depending on context.
Objective:
This research aimed to elucidate HMGB1’s functional significance in cutaneous squamous cell carcinoma (cSCC).
Methods:
We overexpressed HMGB1 in cSCC cell lines using recombinant adenovirus and examined its effects on cell proliferation, colony formation, and cell migration.
Results:
Immunohistochemical analysis revealed elevated HMGB1 expression levels in cSCC tissue relative to normal epidermis. To assess the influence of HMGB1, we employed recombinant adenoviruses expressing HMGB1 to transduce SCC cell lines (SCC12 and SCC13). Enhanced HMGB1 expression significantly promoted cellular proliferation and colony formation capacity.Notably, HMGB1 overexpression elevated the levels of proliferation regulators, including P63, SOX2, CDK4 and CDK6. Furthermore, HMGB1 overexpression substantially enhanced tumor invasiveness, accompanied by upregulation of epithelial-mesenchymal transition (EMT) biomarkers. Mechanistically, overexpression of HMGB1 enhanced transforming growth factor-β signaling by increasing phosphorylation of SMAD2/3, the key mediators of EMT.
Conclusion
These data imply that HMGB1 acts as a tumor-promoting factor in cSCC.
4.Identification and infection rate of densovirus in Culex pipiens pallens in Beijing in 2023
Xiu-yan XU ; Ting YAN ; Si-jie ZHU ; Jing LI ; Mei-de LIU ; Hong-jiang ZHANG ; Ting LIU ; Qiu-hong LI ; Xiao-jie ZHOU ; Ying TONG ; Yong ZHANG
Acta Parasitologica et Medica Entomologica Sinica 2026;33(1):25-30
Objective This study conducted molecular biological identification of the viruses carried by Culex pipiens pallens specimens collected in Shunyi District, Beijing in 2023, and observed the changes in the infection rate of the viruses carried by Cx. pipiens pallens at different collection times. Methods Cx. pipiens pallens were collected using carbon dioxide mosquito traps. The mosquito samples were ground in batches and analyzed by molecular biology technologies. The virus infection rate at different collection times was analyzed statistically. Results 17 strains of Culex pipiens pallens densovirus(CppDNV)were identified from Cx. pipiens pallens samples collected in Shunyi District, Beijing, in 2023. The nucleotide sequence analysis of the virus genome coding region showed that CppDNV was a single-stranded DNA virus with a total length of 3 335 nt, encoding 2 non-structural proteins(NS1, NS2)and 1 capsid protein(VP). The nucleotide(amino acid)sequence lengths of the three proteins were 2 376 nt(791 aa),1 092 nt(363 aa)and 1 071 nt(356 aa), respectively. Phylogenetic analysis showed that CppDNV was located in genus Brevihamaparvovirus. Statistical analysis showed no significant difference in infection rates across collection times(χ2=4.429, P=0.194). Conclusions CppDNV was identified in Cx. pipiens pallens in Beijing, and it was stably maintained in this natural population.
5.Research on Discrimination of Degradation Levels in Shipwreck Archaeological Wood Based on Microscale Attenuated Total Reflection Fourier Transform Infrared Spectroscopy
Ren LI ; Man-Li SUN ; Li-Chao JIAO ; Ya-Fang YIN ; Zhi-Guo ZHANG ; Fu-De TIE
Chinese Journal of Analytical Chemistry 2025;53(6):967-975
After the wooden shipwreck was recovered from the marine underwater environment,the wooden components undergo varying degrees of degradation,therefore,accurately determining the extent of degradation is a fundamental scientific issue for implementing effective preservation strategies.In this work,the wooden remains of Pinus massoniana excavated from the"Nanhai No.1"shipwreck(Southern Song Dynasty)were investigated and compared with the modern wood to discriminate the degradation levels of archaeological wood using attenuated total reflection Fourier transform infrared(ATR-FTIR)spectroscopy.The residual sugar content within wood cell walls was determined using a non-invasive automated microscale ATR-FTIR method to extract chemical information from the wood tangential section.Microstructural characterization of wood samples was conducted by super depth of field microscopy and scanning electron microscopy.FTIR spectral analysis was performed to evaluate the degradation state and elucidate changes in cellulose crystallinity.Finally,the combination of FTIR spectroscopy with the sparse partial least squares discriminant analysis(sPLS-DA)model facilitated the rapid discrimination of degradation levels in shipwreck archaeological wood,and the performance of the model was evaluated using receiver operating characteristic(ROC)curves and area under the curve(AUC).The results showed that the higher the degree of wood degradation,the lower the residual sugar content in the wood cell wall,and the residual glucose content of highly degraded wood was only 4.7%.Significant differences were observed in both the tangential section microstructure and FTIR characteristic absorption patterns across degradation levels,and as the degradation advanced,progressive cell wall loosening occurred alongside selective removal of polysaccharide components,and the relative lignin content was increased,resulting in an elevated A1509/A1370 ratio in FTIR spectra.The sPLS-DA model achieved excellent discrimination performance with AUC values exceeding 0.9,confirming that the combination of FTIR spectroscopy with sPLS-DA enabled accurate assessment of degradation levels in shipwreck archaeological wood.This study developed a rapid and accurate methodology for assessing degradation levels in shipwreck archaeological wood based on microscale ATR-FTIR spectroscopy,which would help to promote the accurate assessment of the preservation state of waterlogged wooden artifacts.
6.Clinical features and immunotherapy for children with loss-of-function/gain-of-function mutations in the STAT gene: an analysis of 10 cases.
Hong-Wei LI ; Yan-Hong WANG ; Shang-Zhi WU ; Bi-Yun ZHANG ; Shi-Hui XU ; Jia-Xing XU ; Zhan-Hang HUANG ; Cheng-Yu LU ; De-Hui CHEN
Chinese Journal of Contemporary Pediatrics 2025;27(8):951-958
OBJECTIVES:
To investigate the clinical features of children with STAT gene mutations, and to explore corresponding immunotherapy strategies.
METHODS:
A retrospective analysis was performed for the clinical data of 10 children with STAT gene mutations who were admitted to the Department of Pediatrics of the First Affiliated Hospital of Guangzhou Medical University, from October 2015 to October 2024. Exploratory immunotherapy was implemented in some refractory cases, and the changes in symptoms, imaging manifestations, and cytokine levels were assessed after treatment.
RESULTS:
For the 10 children, the main clinical manifestations were recurrent rash since birth (7/10), cough (8/10), wheezing (5/10), expectoration (4/10), and purulent nasal discharge (4/10). Genotyping results showed that there was one child with heterozygous loss-of-function (LOF) mutation in the STAT1 gene, four children with heterozygous LOF mutation in the STAT3 gene, and five children with heterozygous gain-of-function (GOF) mutation in the STAT3 gene. Two children with LOF mutation in the STAT3 gene showed decreased interleukin-6 levels and improved clinical symptoms and imaging findings after omalizumab treatment. Three children with GOF mutation in the STAT3 gene achieved effective disease control after treatment with methylprednisolone (0.5 mg/kg per day). Two children with GOF mutation in the STAT3 gene received treatment with JAK inhibitor and then showed some improvement in symptoms.
CONCLUSIONS
STAT gene mutation screening should be considered for children with recurrent rash and purulent respiratory tract infections. Targeted immunotherapy may improve prognosis in patients with no response to conventional treatment.
Humans
;
Male
;
Immunotherapy
;
Female
;
Child, Preschool
;
Child
;
Gain of Function Mutation
;
Retrospective Studies
;
Infant
;
Loss of Function Mutation
;
STAT Transcription Factors/genetics*
7.Exploring urban versus rural disparities in atrial fibrillation: prevalence and management trends among elderly Chinese in a screening study.
Wei ZHANG ; Yi CHEN ; Lei-Xiao HU ; Jia-Hui XIA ; Xiao-Fei YE ; Wen-Yuan-Yue WANG ; Xin-Yu WANG ; Quan-Yong XIANG ; Qin TAN ; Xiao-Long WANG ; Xiao-Min YANG ; De-Chao ZHAO ; Xin CHEN ; Yan LI ; Ji-Guang WANG ; FOR THE IMPRESSION INVESTIGATORS AND COORDINATORS
Journal of Geriatric Cardiology 2025;22(2):246-254
BACKGROUND:
Atrial fibrillation (AF) is a common cardiac arrhythmia in the elderly. This study aimed to evaluate urban-rural disparities in its prevalence and management in elderly Chinese.
METHODS:
Consecutive participants aged ≥ 65 years attending outpatient clinics were enrolled for AF screening using handheld single-lead electrocardiogram (ECG) from April 2017 to December 2022. Each ECG rhythm strip was reviewed from the research team. AF or uninterpretable single-lead ECGs were referred for 12-lead ECG. Primary study outcome comparison was between rural and urban areas for the prevalence of AF. The Student's t-test was used to compare mean values of clinical characteristics between rural and urban participants, while the Pearson's chi-square test was used to compare between-group proportions. Multivariate stepwise logistic regression analysis was performed to estimate the association between AF and various patient characteristics.
RESULTS:
The 29,166 study participants included 13,253 men (45.4%) and had a mean age of 72.2 years. The 7073 rural participants differed significantly (P ≤ 0.02) from the 22,093 urban participants in several major characteristics, such as older age, greater body mass index, and so on. The overall prevalence of AF was 4.6% (n = 1347). AF was more prevalent in 7073 rural participants than 22,093 urban participants (5.6% vs. 4.3%, P < 0.01), before and after adjustment for age, body mass index, blood pressure, pulse rate, cigarette smoking, alcohol consumption and prior medical history. Multivariate logistic regression analysis identified overweight/obesity (OR = 1.35, 95% CI: 1.17-1.54) in urban areas and cigarette smoking (OR = 1.62, 95% CI: 1.20-2.17) and alcohol consumption (OR = 1.42, 95% CI: 1.04-1.93) in rural areas as specific risk factors for prevalent AF. In patients with known AF in urban areas (n = 781) and rural areas (n = 338), 60.6% and 45.9%, respectively, received AF treatment (P < 0.01), and only 22.4% and 17.2%, respectively, received anticoagulation therapy (P = 0.05).
CONCLUSIONS
In China, there are urban-rural disparities in AF in the elderly, with a higher prevalence and worse management in rural areas than urban areas. Our study findings provide insight for health policymakers to consider urban-rural disparity in the prevention and treatment of AF.
8.The Valvular Heart Disease-specific Age-adjusted Comorbidity Index (VHD-ACI) score in patients with moderate or severe valvular heart disease.
Mu-Rong XIE ; Bin ZHANG ; Yun-Qing YE ; Zhe LI ; Qing-Rong LIU ; Zhen-Yan ZHAO ; Jun-Xing LV ; De-Jing FENG ; Qing-Hao ZHAO ; Hai-Tong ZHANG ; Zhen-Ya DUAN ; Bin-Cheng WANG ; Shuai GUO ; Yan-Yan ZHAO ; Run-Lin GAO ; Hai-Yan XU ; Yong-Jian WU
Journal of Geriatric Cardiology 2025;22(9):759-774
BACKGROUND:
Based on the China-VHD database, this study sought to develop and validate a Valvular Heart Disease- specific Age-adjusted Comorbidity Index (VHD-ACI) for predicting mortality risk in patients with VHD.
METHODS & RESULTS:
The China-VHD study was a nationwide, multi-centre multi-centre cohort study enrolling 13,917 patients with moderate or severe VHD across 46 medical centres in China between April-June 2018. After excluding cases with missing key variables, 11,459 patients were retained for final analysis. The primary endpoint was 2-year all-cause mortality, with 941 deaths (10.0%) observed during follow-up. The VHD-ACI was derived after identifying 13 independent mortality predictors: cardiomyopathy, myocardial infarction, chronic obstructive pulmonary disease, pulmonary artery hypertension, low body weight, anaemia, hypoalbuminaemia, renal insufficiency, moderate/severe hepatic dysfunction, heart failure, cancer, NYHA functional class and age. The index exhibited good discrimination (AUC, 0.79) and calibration (Brier score, 0.062) in the total cohort, outperforming both EuroSCORE II and ACCI (P < 0.001 for comparison). Internal validation through 100 bootstrap iterations yielded a C statistic of 0.694 (95% CI: 0.665-0.723) for 2-year mortality prediction. VHD-ACI scores, as a continuous variable (VHD-ACI score: adjusted HR (95% CI): 1.263 (1.245-1.282), P < 0.001) or categorized using thresholds determined by the Yoden index (VHD-ACI ≥ 9 vs. < 9, adjusted HR (95% CI): 6.216 (5.378-7.184), P < 0.001), were independently associated with mortality. The prognostic performance remained consistent across all VHD subtypes (aortic stenosis, aortic regurgitation, mitral stenosis, mitral regurgitation, tricuspid valve disease, mixed aortic/mitral valve disease and multiple VHD), and clinical subgroups stratified by therapeutic strategy, LVEF status (preserved vs. reduced), disease severity and etiology.
CONCLUSION
The VHD-ACI is a simple 13-comorbidity algorithm for the prediction of mortality in VHD patients and providing a simple and rapid tool for risk stratification.
9.Performance assessment of computed tomographic angiography fractional flow reserve using deep learning: SMART trial summary.
Wei ZHANG ; You-Bing YIN ; Zhi-Qiang WANG ; Ying-Xin ZHAO ; Dong-Mei SHI ; Yong-He GUO ; Zhi-Ming ZHOU ; Zhi-Jian WANG ; Shi-Wei YANG ; De-An JIA ; Li-Xia YANG ; Yu-Jie ZHOU
Journal of Geriatric Cardiology 2025;22(9):793-801
BACKGROUND:
Non-invasive computed tomography angiography (CTA)-based fractional flow reserve (CT-FFR) could become a gatekeeper to invasive coronary angiography. Deep learning (DL)-based CT-FFR has shown promise when compared to invasive FFR. To evaluate the performance of a DL-based CT-FFR technique, DeepVessel FFR (DVFFR).
METHODS:
This retrospective study was designed for iScheMia Assessment based on a Retrospective, single-center Trial of CT-FFR (SMART). Patients suspected of stable coronary artery disease (CAD) and undergoing both CTA and invasive FFR examinations were consecutively selected from the Beijing Anzhen Hospital between January 1, 2016 to December 30, 2018. FFR obtained during invasive coronary angiography was used as the reference standard. DVFFR was calculated blindly using a DL-based CT-FFR approach that utilized the complete tree structure of the coronary arteries.
RESULTS:
Three hundred and thirty nine patients (60.5 ±10.0 years and 209 men) and 414 vessels with direct invasive FFR were included in the analysis. At per-vessel level, sensitivity, specificity, accuracy, positive predictive value (PPV) and negative predictive value (NPV) of DVFFR were 94.7%, 88.6%, 90.8%, 82.7%, and 96.7%, respectively. The area under the receiver operating characteristics curve (AUC) was 0.95 for DVFFR and 0.56 for CTA-based assessment with a significant difference (P < 0.0001). At patient level, sensitivity, specificity, accuracy, PPV and NPV of DVFFR were 93.8%, 88.0%, 90.3%, 83.0%, and 95.8%, respectively. The computation for DVFFR was fast with the average time of 22.5 ± 1.9 s.
CONCLUSIONS
The results demonstrate that DVFFR was able to evaluate lesion hemodynamic significance accurately and effectively with improved diagnostic performance over CTA alone. Coronary artery disease (CAD) is a critical disease in which coronary artery luminal narrowing may result in myocardial ischemia. Early and effective assessment of myocardial ischemia is essential for optimal treatment planning so as to improve the quality of life and reduce medical costs.
10.Exploring Mechanism of Electroacupuncture in Modulating Neuroinflammation Based on Intestinal Flora and Its Metabolites.
Hai-Min YE ; Zhuo-Yan LI ; Peng ZHANG ; Zhen KANG ; De-Sheng ZHOU
Chinese journal of integrative medicine 2025;31(2):183-192
Neuroinflammatory responses play an important role in the pathogenesis of various diseases, particularly those affecting the central nervous system. Inhibition of neuroinflammation is a crucial therapeutic strategy for the management of central nervous system disorders. The intestinal microbial-gut-brain axis serves as a key regulatory pathway that modulates neuroinflammatory processes. Intestinal flora metabolites such as short-chain fatty acids, indoles and their derivatives, lipopolysaccharides, trimethylamine oxide, and secondary bile acids exert direct or indirect effects on neuroinflammation. Studies have shown that electroacupuncture (EA) modulates the composition of the intestinal microbiota and its metabolites, while also suppressing neuroinflammation by targeting the TLR4/NF- κ B, NLRP3/caspase-1, and microglial cell M2-type transformation pathways. This review discusses the mechanisms by which EA regulates neuroinflammation via intestinal microbiota and its metabolites, providing information and a foundation for further investigation of the precise therapeutic mechanisms of EA in neurological disorders.
Humans
;
Gastrointestinal Microbiome
;
Electroacupuncture
;
Neuroinflammatory Diseases/metabolism*
;
Animals


Result Analysis
Print
Save
E-mail