1.DYRK2:a novel therapeutic target for rheumatoid arthritis combined with osteoporosis based on East Asian and European populations
Zhilin WU ; Qin HE ; Pingxi WANG ; Xian SHI ; Song YUAN ; Jun ZHANG ; Hao WANG
Chinese Journal of Tissue Engineering Research 2026;30(6):1569-1579
BACKGROUND:Studies have shown that rheumatoid arthritis and osteoporosis are positively correlated,but the causal relationship and related mechanisms have not yet been confirmed.With the cross-fertilization of computer science and life sciences,Mendelian randomization and bioinformatics analyses based on genome-wide association study(GWAS)and transcriptome sequencing data can assess the causal relationship between two diseases,explore the related mechanisms,and mine the therapeutic targets,which will be beneficial to the precision treatment of rheumatoid arthritis combined with osteoporosis.OBJECTIVE:To explore the causal relationship between rheumatoid arthritis and osteoporosis using two-sample Mendelian randomization and to mine potential co-morbid targets and potential targeted drugs through summary-data-based Mendelian randomization and bioinformatics analyses,aiming to provide theoretical basis for mechanism exploration and precision treatment in the field of rheumatoid arthritis combined with osteoporosis.METHODS:(1)Firstly,GWAS data of rheumatoid arthritis,osteoporosis,and cis-expression quantitative trait locus(cis-eQTL)in Asian and European populations were downloaded from the GWAS Catalog,IEU Open GWAS,FinnGen,and eQTLGen databases,and were used for two-sample Mendelian randomization analysis and summary-data-based Mendelian randomization analysis.(2)Transcriptome sequencing data of rheumatoid arthritis(GSE93272 and GSE15573)were downloaded from the GEO database for bioinformatics analysis.(3)Subsequently,forward and inverse Mendelian randomization analyses between rheumatoid arthritis and osteoporosis were performed,and inverse variance weighted was used as the main metric for the analyses,and the results were corroborated with MR Egger,simple mode,weighted median and weighted mode.(4)Then,the genes closely related to rheumatoid arthritis and osteoporosis were identified based on the summary-data-based Mendelian randomization analysis,and the co-disease targets of rheumatoid arthritis and osteoporosis were mined based on cross-analysis.Meanwhile,the biological functions of the co-morbid targets were verified based on bioinformatics analysis and cellular experiments.(5)In addition,a rheumatoid arthritis risk prediction nomogram was constructed based on DYRK2,and its prediction performance was verified by receiver operating characteristic curve,correction curve and decision curve.Finally,the target potential drugs were mined based on Enrichr database and molecular docking was performed.RESULTS AND CONCLUSION:(1)Forward Mendelian randomization analysis of rheumatoid arthritis and osteoporosis showed statistically significant results except for GCST90044540 and GCST90086118,and all other results indicated a significant causal relationship and positive correlation between rheumatoid arthritis and osteoporosis.(2)Inverse Mendelian randomization analysis suggested that no significant causal relationship was seen between osteoporosis and rheumatoid arthritis.(3)Summary-data-based Mendelian randomization analysis identified a total of 412 and 344 genes positively associated with rheumatoid arthritis and osteoporosis,and 421 and 347 genes negatively associated.Based on the cross-analysis,26 co-morbid genes were subsequently obtained.Among them,DYRK2 was a potential therapeutic target,and subsequent bioinformatics analysis and cellular experiments confirmed its important role in the progression of rheumatoid arthritis and osteoporosis.(4)Furthermore,the constructed nomogram has excellent predictive performance.Finally,four potential DYRK2-targeting drugs(undecanoic acid,metyrapone,JNJ-38877605,and ACA)were discovered and molecular docking also demonstrated reliable targeting ability.(5)In conclusion,based on GWAS data from Asian and European populations,we successfully demonstrated that rheumatoid arthritis and osteoporosis are causally related at the genetic level,DYRK2 is a potential therapeutic target,and four small molecules are potential target drugs.
2.DYRK2:a novel therapeutic target for rheumatoid arthritis combined with osteoporosis based on East Asian and European populations
Zhilin WU ; Qin HE ; Pingxi WANG ; Xian SHI ; Song YUAN ; Jun ZHANG ; Hao WANG
Chinese Journal of Tissue Engineering Research 2026;30(6):1569-1579
BACKGROUND:Studies have shown that rheumatoid arthritis and osteoporosis are positively correlated,but the causal relationship and related mechanisms have not yet been confirmed.With the cross-fertilization of computer science and life sciences,Mendelian randomization and bioinformatics analyses based on genome-wide association study(GWAS)and transcriptome sequencing data can assess the causal relationship between two diseases,explore the related mechanisms,and mine the therapeutic targets,which will be beneficial to the precision treatment of rheumatoid arthritis combined with osteoporosis.OBJECTIVE:To explore the causal relationship between rheumatoid arthritis and osteoporosis using two-sample Mendelian randomization and to mine potential co-morbid targets and potential targeted drugs through summary-data-based Mendelian randomization and bioinformatics analyses,aiming to provide theoretical basis for mechanism exploration and precision treatment in the field of rheumatoid arthritis combined with osteoporosis.METHODS:(1)Firstly,GWAS data of rheumatoid arthritis,osteoporosis,and cis-expression quantitative trait locus(cis-eQTL)in Asian and European populations were downloaded from the GWAS Catalog,IEU Open GWAS,FinnGen,and eQTLGen databases,and were used for two-sample Mendelian randomization analysis and summary-data-based Mendelian randomization analysis.(2)Transcriptome sequencing data of rheumatoid arthritis(GSE93272 and GSE15573)were downloaded from the GEO database for bioinformatics analysis.(3)Subsequently,forward and inverse Mendelian randomization analyses between rheumatoid arthritis and osteoporosis were performed,and inverse variance weighted was used as the main metric for the analyses,and the results were corroborated with MR Egger,simple mode,weighted median and weighted mode.(4)Then,the genes closely related to rheumatoid arthritis and osteoporosis were identified based on the summary-data-based Mendelian randomization analysis,and the co-disease targets of rheumatoid arthritis and osteoporosis were mined based on cross-analysis.Meanwhile,the biological functions of the co-morbid targets were verified based on bioinformatics analysis and cellular experiments.(5)In addition,a rheumatoid arthritis risk prediction nomogram was constructed based on DYRK2,and its prediction performance was verified by receiver operating characteristic curve,correction curve and decision curve.Finally,the target potential drugs were mined based on Enrichr database and molecular docking was performed.RESULTS AND CONCLUSION:(1)Forward Mendelian randomization analysis of rheumatoid arthritis and osteoporosis showed statistically significant results except for GCST90044540 and GCST90086118,and all other results indicated a significant causal relationship and positive correlation between rheumatoid arthritis and osteoporosis.(2)Inverse Mendelian randomization analysis suggested that no significant causal relationship was seen between osteoporosis and rheumatoid arthritis.(3)Summary-data-based Mendelian randomization analysis identified a total of 412 and 344 genes positively associated with rheumatoid arthritis and osteoporosis,and 421 and 347 genes negatively associated.Based on the cross-analysis,26 co-morbid genes were subsequently obtained.Among them,DYRK2 was a potential therapeutic target,and subsequent bioinformatics analysis and cellular experiments confirmed its important role in the progression of rheumatoid arthritis and osteoporosis.(4)Furthermore,the constructed nomogram has excellent predictive performance.Finally,four potential DYRK2-targeting drugs(undecanoic acid,metyrapone,JNJ-38877605,and ACA)were discovered and molecular docking also demonstrated reliable targeting ability.(5)In conclusion,based on GWAS data from Asian and European populations,we successfully demonstrated that rheumatoid arthritis and osteoporosis are causally related at the genetic level,DYRK2 is a potential therapeutic target,and four small molecules are potential target drugs.
3.Correlation between liver fibrosis degree and carotid plaque in patients with lean metabolic dysfunction-associated fatty liver disease
Shuai ZHANG ; Shoulu JIN ; Wanqing LI ; Xijing SHI ; Hao LIANG ; Hao DONG ; Dailong LU ; Ying ZHU ; Xiaoxing XIANG ; Jun LIU
Journal of Clinical Hepatology 2026;42(2):319-325
ObjectiveTo investigate the association between noninvasive liver fibrosis markers and carotid plaque (CP) in patients with lean metabolic dysfunction-associated fatty liver disease (MAFLD), and to provide a basis for screening high-risk populations. MethodsA total of 957 patients with lean MAFLD who underwent physical examination in Subei People’s Hospital from January 2021 to June 2023 was enrolled as the observation cohort, with the presence or absence of CP as the outcome, and fibrosis-4 (FIB-4) index and nonalcoholic fatty liver disease fibrosis score (NFS) were used to assess liver fibrosis degree. The independent-samples t test was used for comparison of normally distributed continuous data between two groups, and the Mann-Whitney U test was used for comparison of non-normally distributed continuous data between two groups; the chi-square test was used for comparison of categorical data between two groups. The multivariate logistic regression analysis, the restricted cubic spline analysis, the receiver operating characteristic curve, and the mediation effect analysis were used to investigate the association between liver fibrosis degree and CP. ResultsThe prevalence rate of CP was 36.6% in the lean MAFLD population. Compared with the non-CP group(n=607), the CP group (n=350) had a significantly higher proportion of male patients, a significantly higher proportion of patients with smoking/diabetes/hypertension, and significantly higher levels of age, creatinine, blood urea nitrogen, triglycerides, fasting blood glucose, aspartate aminotransferase, aspartate aminotransferase/alanine aminotransferase ratio, NFS, and FIB-4 index, as well as significantly lower levels of platelet count and albumin (all P<0.05). The multivariate logistic regression analysis showed that after adjustment for confounding factors, FIB-4 index (odds ratio[OR]=2.979, 95% confidence interval[CI]:2.141 — 4.219, P<0.001) and NFS (OR=1.747, 95%CI: 1.499 — 2.046, P<0.001) were positively correlated with CP. Both FIB-4 index and NFS had a good value in predicting CP. Hypertension had a significant indirect effect on the prevalence rate of CP through its impact on liver fibrosis markers, and its mediating effect accounted for 39.5% — 40.8% of the total effect (P<0.001). ConclusionIn patients with lean MAFLD, NFS and FIB-4 index are significantly positively correlated with the prevalence rate of CP, and they can be used as potential epidemiological predictive indicators. Liver fibrosis markers may play a mediating role in the association between hypertension and CP. Interventions targeting hypertension and liver fibrosis markers may help to prevent and delay the progression of CP.
4.Association of liver fibrosis markers and inflammation markers with the risk of gallstones in patients with metabolic dysfunction-associated fatty liver disease
Shuai ZHANG ; Shoulu JIN ; Wanqing LI ; Xijing SHI ; Hao LIANG ; Hao DONG ; Dailong LU ; Ying ZHU ; Xiaoxing XIANG ; Jun LIU
Journal of Clinical Hepatology 2026;42(3):579-585
ObjectiveTo investigate the association of liver fibrosis scores and inflammation markers with gallstones in patients with metabolic dysfunction-associated fatty liver disease (MAFLD), as well as the mediating role of liver fibrosis scores in the relationship between inflammation markers and gallstones. MethodsA total of 14 567 patients who received physical examination and were diagnosed with MAFLD in Subei People’s Hospital from January 2014 to June 2023 were enrolled in this study, and according to the results of abdominal color Doppler ultrasound, they were divided into gallstone group with 1 724 patients and non-gallstone group with 12 843 patients. Related clinical data were collected from all patients, including demographic data, medical history, family history, physical examination, Color Doppler ultrasound, and biochemical parameters. The biomarkers associated with metabolic disorders and insulin resistance included triglyceride-glucose index (TyG), TyG-body mass index (BMI) index, atherogenic index of plasma (AIP), and non-high-density lipoprotein cholesterol-to-high-density lipoprotein cholesterol ratio (NHHR); the biomarkers associated with inflammation and nutritional status included neutrophil-to-lymphocyte ratio (NLR), neutrophil percentage-to-albumin ratio (NPAR), and monocyte-to-lymphocyte ratio (MLR); the biomarkers for assessing liver fibrosis degree and liver function included albumin-bilirubin (ALBI) score, NAFLD fibrosis score (NFS), fibrosis-4 (FIB-4) index, and aspartate aminotransferase-to-platelet ratio index (APRI). The independent-samples t test was used for comparison of normally distributed continuous data between two groups, while the Mann-Whitney U test was used for comparison of non-normally distributed continuous data between two groups; the chi-square test was used for comparison of categorical data between two groups. Multivariate Logistic regression analysis, restricted cubic spline analysis, and mediating effect analysis were used to assess the association of liver fibrosis markers and inflammation markers with the risk of gallstones. ResultsThe prevalence rate of gallstones was 11.8% among the MAFLD patients. There were significant differences between the gallstone group and the non-gallstone group in sex, age, smoking history, diabetes, hypertension, lymphocytes, platelets, glucose, albumin, serum uric acid, alanine aminotransferase, aspartate aminotransferase, red blood cell, NLR, NPAR, MLR, NFS, FIB-4 index, and ALBI score (all P<0.05). The multivariate Logistic regression analysis showed that NLR (odds ratio [OR]=1.091, 95% confidence interval [CI]: 1.028 — 1.160, P<0.05), NPAR (OR=1.073, 95%CI: 1.042 — 1.105, P<0.05), MLR (OR=1.142, 95%CI: 1.057 — 1.232, P<0.05), NFS (OR=1.239, 95%CI: 1.190 — 1.291, P<0.05), and FIB-4 index (OR=1.326, 95%CI: 1.241 — 1.417, P<0.05) were influencing factors for the prevalence rate of gallstones. The restricted cubic spline analysis showed a significant non-linear association between NFS/FIB-4 index and the risk of gallstone (non-linear P<0.05). The mediating effect analysis further showed that the association of NLR, MLR, and NPAR with gallstones was partially mediated by NFS or FIB-4 index, with a mediating effect accounting for 36.79%、28.09%、29.67% and 18.31%、17.70、11.57%, respectively. ConclusionNFS and FIB-4 index have a non-linear association with the prevalence rate of gallstones in MAFLD patients, and they also mediate the association of NLR, NPAR, and MLR with the risk of gallstone.
5.Construction of a risk prediction model for cardiovascular events in community hypertensive patients based on remote ambulatory blood pressure parameters
Guiqiu ZHU ; Yihong WU ; Hao ZHANG ; Jun SUN ; Yajuan ZHANG ; Xiaohong WANG ; Zongquan ZHAO
Journal of Public Health and Preventive Medicine 2026;37(3):85-89
Objective To explore the risk prediction model of major adverse cardiovascular events (MACE) in community patients with hypertension based on remote ambulatory blood pressure parameters. Methods From November 2023 to October 2024, 486 community patients with hypertension who received standardized management in Nanjing Medical University Affiliated to Suzhou Hospital were retrospectively selected. All patients wore remote ambulatory blood pressure monitor to obtain 24-hour ambulatory blood pressure data. Clinical data were collected and remote ambulatory blood pressure parameters [24-hour systolic blood pressure variability (SBPV), 24-hour diastolic blood pressure variability (DBPV), nighttime SBPV, nighttime DBPV, daytime SBPV, daytime DBPV] were extracted. The patients were followed up for 12 months, and were classified into MACE group (n=42) and non-MACE group (n=444) according to whether MACE occurred during follow-up. Multivariate Logistic regression analysis was adopted to screen the influencing factors for MACE. Based on the above factors, a risk prediction model was constructed and verified by receiver operating characteristic (ROC) curve. Results MACE occurred in 42 cases among 486 patients, with an incidence rate of 8.64%. Multivariate Logistic regression analysis suggested that nighttime DBPV (OR=1.119, 95%CI: 1.030-1.214), 24h-SBPV (OR=1.115, 95%CI: 1.007-1.235), nighttime SBPV (OR=1.116, 95%CI: 1.016-1.226) and diabetes mellitus (OR=2.762, 95%CI: 1.059-7.203) were independent factors for MACE (P<0.05). The model validation results revealed that the area under the ROC curve was 0.905 (95%CI: 0.854-0.956 ), and the model had a good discrimination degree. Conclusion Nighttime DBPV, 24h-SBPV, nighttime SBPV and diabetes mellitus are independent risk factors for MACE in community patients with hypertension. The clinical prediction model based on these variables exhibits certain predictive value on MACE risk.
6.Investigation of the prevalence of cognitive frailty in patients with chronic obstructive pulmonary disease and analysis of associated factors
Yihong WU ; Hao ZHANG ; Jun SUN ; Yajuan ZHANG ; Jufan TIAN ; Guiqiu ZHU ; Zongquan ZHAO
Journal of Public Health and Preventive Medicine 2026;37(3):99-103
Objective To investigate the prevalence of cognitive frailty and its associated factors among middle-aged and elderly patients with chronic obstructive pulmonary disease (COPD) in Gusu District, Suzhou City, to provide a basis for intervention measures. Methods COPD patients managed at Pingjiang New Town Community Health Service Center in Gusu District from January to December 2024 were selected. Data were collected using the Fried Frailty Phenotype, Mini-Mental State Examination, and a general information questionnaire. Chi-square tests, one-way analysis of variance, and multivariate Logistic regression were used for data analysis. Results A total of 1 388 patients were included, with 916 males (65.99%) and 472 females (34.01%). Cognitive frailty was observed in 554 cases (prevalence: 39.91%). Univariate analysis showed significant differences in gender, education, visual impairment, hearing impairment, sleep disorders, and disease duration between the cognitive frailty and non-cognitive frailty groups (all P<0.05). Multivariate Logistic regression indicated that education was negatively associated with cognitive frailty (illiteracy OR=6.969, P<0.001; primary school OR=1.499, P=0.005), while visual impairment (OR=2.347, P<0.001), hearing impairment (OR=1.705, P<0.001), sleep disorders (OR=2.488, P<0.001), and disease duration (≤5 years OR=0.329, P<0.001; 6-10 years OR=0.487, P<0.001) were positively associated with cognitive frailty. Conclusion The prevalence of cognitive frailty among middle-aged and elderly COPD patients in Gusu District, Suzhou City, is high (39.91%). Epositively associated with cognitive frailtyducation, sensory impairments, sleep disorders, and disease duration are associated factors. Clinical attention to high-risk groups and implementation of early screening and targeted interventions are recommended.
7.Analysis of the frequency of X-ray diagnostic examinations and CT radiation doses in public hospitals of a district in Ningbo City, China
Shuxia HAO ; Mengxue LI ; Yong WANG ; Shengnan FAN ; Jingguo ZHANG ; Xueying WANG ; Jun DENG ; Quanfu SUN
Chinese Journal of Radiological Health 2025;34(3):324-330
Objective To systematically analyze the medical radiation exposure levels in a district of Ningbo City and to provide a scientific basis for the reasonable and effective control of medical radiation exposure. Methods Based on the radiological diagnosis frequency and dose information system, basic medical radiation exposure data were collected, such as radiation doses received by patients in various X-ray diagnostic examinations, from all 13 public medical institutions in a district of Ningbo City from January 1 to December 31, 2020. The data were analyzed for the percentage and collective effective dose of various diagnostic examinations, the distribution of examinations by sex and age, and the number of patients undergoing two or more examinations and their cumulative doses within multiple time intervals. Results Among medical X-ray diagnostic examinations in the district, the percentages of CT examination and routine photography examination were 50.88% and 47.93%, respectively, and the collective effective dose of CT examination accounted for 97.75%. By age and sex, the frequency of examination was the highest in the age group of 45-54 years, and the frequency of examination in the male was higher than that in the female before age 55. The annual effective dose for two patients exceeded 100 mSv. Conclusion In this study, CT examination accounted for up to 50.88% of all medical X-ray diagnostic examinations, and contributed 97.75% of the collective effective dose, highlighting the need for particular attention to the justification of medical radiation exposure from CT.
8.Status quo of cognitive frailty in community elderly patients with chronic obstructive pulmonary disease and its association with sleep quality, anxiety and depression
Zongquan ZHAO ; Hao ZHANG ; Jun SUN ; Yajuan ZHANG ; Xiangfan TIAN ; Xiaohong WANG ; Zhenyuan TANG ; Zhiying ZHANG
Journal of Public Health and Preventive Medicine 2025;36(5):71-75
Objective To analyze the status quo of cognitive frailty (CF) in community elderly patients with chronic obstructive pulmonary disease (COPD) and its correlation with sleep quality, anxiety and depression. Methods Elderly patients with COPD receiving health management in the center were selected from July 2023 to June 2024. The general data of patients were collected and Mini-Mental State Examination (MMSE), Fried Frailty Phenotype (FP), Pittsburgh Sleep Quality Index (PSQI) and Hospital Anxiety and Depression Scale (HADS) were used for investigation, and the above scores were analyzed. All patients were divided into CF group (n=129) and non-CF group (n=319) according to MMSE score and FP score. Univariate and multivariate logistic analyses were used to analyze the influencing factors of CF in elderly COPD patients. Results Pearson correlation analysis showed that MMSE score was significantly negatively correlated with PSQI score and HADS score in elderly COPD patients (P<0.05), and FP score was significantly positively correlated with PSQI score and HADS score (P<0.05). After logistic regression analysis, it was found that education level, marital status and sleep time were protective factors of CF in elderly COPD patients (P<0.05), and PSQI score and HADS score were risk factors of CF in elderly patients with COPD (P<0.05). Conclusion CF in community elderly COPD patients is related to sleep quality, sleep duration and anxiety and depression. It is necessary to take clinical measures to improve the sleep quality and psychological status, so as to avoid or slow down the occurrence of CF.
9.Treatment of erectile dysfunction based on the "brain-heart-kidney-essence chamber" axis and the meridian-zangfu relationship.
Dicheng LUO ; Jun GUO ; Hao WANG ; Dongyue MA ; Ziwei ZHAO ; Yang LIU ; Hongyuan CHANG ; Jiwei ZHANG ; Wenxiao YU
Chinese Acupuncture & Moxibustion 2025;45(5):609-613
Based on the pathogenesis of erectile dysfunction (ED) from the meridian-zangfu relationship and the "brain-heart-kidney-essence chamber" axis, it proposes that dysfunction of the "brain-heart-kidney-essence chamber" axis is closely related to the occurrence of ED. Among these, brain-heart disharmony is the key pathogenic factor, kidney deficiency and essence depletion constitute an important basis, and essence chamber stasis is a critical mechanism. The treatment approach emphasizes harmonizing the brain and heart, regulating the mind, tonifying the kidney and replenishing qi, unblocking qi and blood to harmonize the essence chamber. The primary acupoints include Baihui (GV20)-Neiguan (PC6)-Shenmen (HT7), Taixi (KI3)-Guanyuan (CV4)-Sanyinjiao (SP6), and Zhongji (CV3)-Dahe (KI12)-Gongsun (SP4), with additional acupoints selected based on syndrome differentiation. This approach aims to restore the clarity of the brain and heart, replenish kidney qi, and unblock the essence chamber, thereby facilitating the restoration of normal functions of the brain, heart, kidney, and essence chamber, and alleviating ED symptoms and improving overall clinical efficacy.
Humans
;
Male
;
Meridians
;
Erectile Dysfunction/physiopathology*
;
Kidney/physiopathology*
;
Brain/physiopathology*
;
Acupuncture Therapy
;
Acupuncture Points
;
Heart/physiopathology*
10.Discussion on the Treatment of Coronary Heart Disease with Erectile Dysfunction Based on the"Brain-Heart-Kidney-Seminal Chamber"Axis Theory
Dicheng LUO ; Liping PAN ; Hao WANG ; Yang LIU ; Yunzhi LI ; Jiwei ZHANG ; Shengjing LIU ; Jun GUO
Chinese Journal of Information on Traditional Chinese Medicine 2025;32(7):167-171
Based on the"brain-heart-kidney-seminal chamber"axis theory to understand the pathogenesis of coronary heart disease(CHD)complicated with erectile dysfunction(ED),it is considered that the dysfunction of the"brain-heart-kidney-seminal chamber"axis is closely related to the occurrence of CHD complicated with ED,in which kidney deficiency is the initial factor;heart brain disorder,loss of consciousness is the key to the pathogenesis;the loss of seminal chamber is an important link in the pathogenesis.Under the guidance of the concept of"brain-heart-kidney-seminal chamber"axis,Cistanches Herba-Rehmannize Radix et Praeparata-Cuscutae Semen were used to invigorate the kidney and regulate the essence of the brain,heart,kidney and seminal chamber;Salviea Miltiorrhizae Radix et Rhizoma-Chuanxiong Rhizoma-Schisandrae Chinensis Fructus can reach the heart and brain,calm the heart and calm the mind;Cyathulae Radix-Scolopendra-Spatholobi Caulis can adjust and restore the seminal chamber,moistening the tendons and vibrating the flaccidity.In clinical practice,the whole medication should be used to make the brain,heart and kidney in harmony with the seminal chamber at the same time,and the drugs should be flexibly added and subtracted according to the specific symptoms of patients to enhance the curative effect of CHD patients with ED.


Result Analysis
Print
Save
E-mail