1.Comparison of sleep EEG power spectral density between depressive episode patients and schizophrenia patients with suicidal behavior
Jingwen LIU ; Yunfei ZHOU ; Jingchu HU ; Jiaoyan ZHOU ; Junwei YANG ; Jie LIANG ; Hong XU ; Yu CANG ; Shimeng MA
Sichuan Mental Health 2026;39(1):50-57
BackgroundPatients with depressive episode and schizophrenia have a high risk of suicide. The sleep electroencephalogram power spectral density characteristics of patients with depressive episode accompanied by suicidal behavior and those with schizophrenia may be different, but there is currently a lack of direct comparative studies on these two groups of patients. ObjectiveTo compare the sleep electroencephalogram power spectral density between depressive episode and schizophrenic patients with suicidal behavior, in order to provide references for exploring predictive indicators of suicidal behavior. MethodsFrom June 2018 to December 2020, 20 patients with depressive episode and 20 patients with schizophrenia who had committed suicide within the past month and were treated at the outpatient department of Shenzhen Kangning Hospital were selected. All of them met the diagnostic criteria for depressive episode or schizophrenia as defined in the International Classification of Diseases, tenth edition (ICD-10). Using a random sampling method, 20 volunteers with matching gender and age to the patient groups were selected from the Cuiping community in Shenzhen as the control group. The subjective sleep of the patients was evaluated using the Insomnia Severity Index (ISI), the Dysfunctional Belief and Attitude about Sleep (DBAS), the Disturbing Dreams and Nightmare Severity Index (DDNSI), and the Epworth Somnolence Scale (ESS). The objective sleep of the patients was assessed using polysomnography. The sleep electroencephalogram was filtered and the power spectral density of the brain wave was analyzed and processed for all the subjects. The subjective and objective sleep conditions of the two patient groups were compared, and the sleep electroencephalogram power spectral density of the patient groups and the control group were also compared. ResultsA comparison of subjective and objective sleep conditions between patients with depressive episode accompanied by suicidal behavior and patients with schizophrenia accompanied by suicidal behavior showed no statistically significant differences (P>0.05). Comparisons of sleep electroencephalogram power spectral density in the W stage (average power of α wave, total power of δ wave, average power of δ wave, average power of θ wave), N1 stage (average power of β wave, total power of α wave, total power of δ wave), N2 stage (total power of α wave, average power of α wave, total power of δ wave, average power of δ wave), N3 stage (average power of α wave, average power of δ wave), and R stage (total power of α wave, average power of α wave, total power of δ wave, average power of δ wave) between patients with depressive episode accompanied by suicidal behavior, patients with schizophrenia accompanied by suicidal behavior, and the control group showed statistically significant differences (P<0.05 or 0.01). The total power of δ wave in the W stage and the average power of β wave and δ wave in the N1 stage were higher in two patient groups were higher than those of the control group. The total power of α wave and the average power of α wave in the N2 stage were lower than those of the control group, while the average power of δ wave was higher than that of the control group. The average power of α wave in the N3 stage of both patient groups were lower than that of the control group, while the average power of δ wave was higher than that of the control group. The total power and average power of α wave in the R stage were lower than those of the control group, while the total power and average power of δ wave were higher than those of the control group. All the differences were statistically significant. Patients with depressive episode accompanied by suicidal behavior had higher average powers of α wave, δ wave, and θ wave in the W stage compared with the control group, while the total power of α wave in the N1 stage was lower in the former group. All these differences were statistically significant (P<0.05). ConclusionThe depressive episode patients accompanied by suicidal behavior have highly overlapping sleep electroencephalogram abnormal patterns with those of schizophrenia patients, mainly manifested as a general decrease in α wave power (N2, N3, R stage) and a general increase in δ wave power (W, N1, N2, N3, R stage) as well as β wave power in N1 stage. At the same time, patients with depressive episode accompanied by suicidal behavior also show specific changes, including an increase in the average power of α and θ waves during the wakefulness period (W stage), and a decrease in the total power of α wave in N1 stage. [Funded by Guangdong Province High-level Clinical Key Specialty (with supporting funds from Shenzhen City) (number, SZGSP013); Shenzhen Key Medical Discipline (number, SZXK041); Shenzhen Clinical Medicine Research Center Project (number, 20210617155253001)]
2.Plasma miRNA testing in the differential diagnosis of very early-stage hepatocellular carcinoma: a multicenter real-world study
Jie HU ; Ying XU ; Ao HUANG ; Lei YU ; Zheng WANG ; Xiaoying WANG ; Xinrong YANG ; Zhenbin DING ; Qinghai YE ; Yinghong SHI ; Shuangjian QIU ; Huichuan SUN ; Qiang GAO ; Jia FAN ; Jian ZHOU
Chinese Journal of Clinical Medicine 2025;32(3):350-354
Objective To explore the application of plasma 7 microRNA (miR7) testing in the differential diagnosis of very early-stage hepatocellular carcinoma (HCC). Methods This study is a multicenter real-world study. Patients with single hepatic lesion (maximum diameter≤2 cm) who underwent plasma miR7 testing at Zhongshan Hospital, Fudan University, Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Anhui Provincial Hospital, and Peking University People’s Hospital between January 2019 and December 2024 were retrospectively enrolled. Patients were divided into very early-stage HCC group and non-HCC group, and the clinical pathological characteristics of the two groups were compared. The value of plasma miR7 levels, alpha-fetoprotein (AFP), and des-gamma-carboxy prothrombin (DCP) in the differential diagnosis of very early-stage HCC was evaluated using receiver operating characteristic (ROC) curves and area under the curve (AUC). In patients with both negative AFP and DCP (AFP<20 ng/mL, DCP<40 mAU/mL), the diagnostic value of plasma miR7 for very early-stage HCC was analyzed. Results A total of 64 528 patients from 4 hospitals underwent miR7 testing, and 1 682 were finally included, of which 1 073 were diagnosed with very early-stage HCC and 609 were diagnosed with non-HCC. The positive rate of miR7 in HCC patients was significantly higher than that in non-HCC patients (67.9% vs 24.3%, P<0.001). ROC curves showed that the AUCs for miR7, AFP, and DCP in distinguishing HCC patients from the non-HCC individuals were 0.718, 0.682, and 0.642, respectively. The sensitivities were 67.85%, 43.71%, and 44.45%, and the specificities were 75.70%, 92.78%, and 83.91%, respectively. The pairwise comparison of AUCs showed that the diagnostic efficacy of plasma miR7 detection was significantly better than that of AFP or DCP (P<0.05). Although its specificity was slightly lower than AFP and DCP, the sensitivity was significantly higher. Among patients negative for both AFP and DCP, miR7 maintained an AUC of 0.728 for diagnosing very early-stage HCC, with 67.82% sensitivity and 77.73% specificity. Conclusions Plasma miR7 testing is a potential molecular marker with high sensitivity and specificity for the differential diagnosis of small hepatic nodules. In patients with very early-stage HCC lacking effective molecular markers (negative for both AFP and DCP), miR7 can serve as a novel and effective molecular marker to assist diagnosis.
3.Prediction of Protein Thermodynamic Stability Based on Artificial Intelligence
Lin-Jie TAO ; Fan-Ding XU ; Yu GUO ; Jian-Gang LONG ; Zhuo-Yang LU
Progress in Biochemistry and Biophysics 2025;52(8):1972-1985
In recent years, the application of artificial intelligence (AI) in the field of biology has witnessed remarkable advancements. Among these, the most notable achievements have emerged in the domain of protein structure prediction and design, with AlphaFold and related innovations earning the 2024 Nobel Prize in Chemistry. These breakthroughs have transformed our ability to understand protein folding and molecular interactions, marking a pivotal milestone in computational biology. Looking ahead, it is foreseeable that the accurate prediction of various physicochemical properties of proteins—beyond static structure—will become the next critical frontier in this rapidly evolving field. One of the most important protein properties is thermodynamic stability, which refers to a protein’s ability to maintain its native conformation under physiological or stress conditions. Accurate prediction of protein stability, especially upon single-point mutations, plays a vital role in numerous scientific and industrial domains. These include understanding the molecular basis of disease, rational drug design, development of therapeutic proteins, design of more robust industrial enzymes, and engineering of biosensors. Consequently, the ability to reliably forecast the stability changes caused by mutations has broad and transformative implications across biomedical and biotechnological applications. Historically, protein stability was assessed via experimental methods such as differential scanning calorimetry (DSC) and circular dichroism (CD), which, while precise, are time-consuming and resource-intensive. This prompted the development of computational approaches, including empirical energy functions and physics-based simulations. However, these traditional models often fall short in capturing the complex, high-dimensional nature of protein conformational landscapes and mutational effects. Recent advances in machine learning (ML) have significantly improved predictive performance in this area. Early ML models used handcrafted features derived from sequence and structure, whereas modern deep learning models leverage massive datasets and learn representations directly from data. Deep neural networks (DNNs), graph neural networks (GNNs), and attention-based architectures such as transformers have shown particular promise. GNNs, in particular, excel at modeling spatial and topological relationships in molecular structures, making them well-suited for protein modeling tasks. Furthermore, attention mechanisms enable models to dynamically weigh the contribution of specific residues or regions, capturing long-range interactions and allosteric effects. Nevertheless, several key challenges remain. These include the imbalance and scarcity of high-quality experimental datasets, particularly for rare or functionally significant mutations, which can lead to biased or overfitted models. Additionally, the inherently dynamic nature of proteins—their conformational flexibility and context-dependent behavior—is difficult to encode in static structural representations. Current models often rely on a single structure or average conformation, which may overlook important aspects of stability modulation. Efforts are ongoing to incorporate multi-conformational ensembles, molecular dynamics simulations, and physics-informed learning frameworks into predictive models. This paper presents a comprehensive review of the evolution of protein thermodynamic stability prediction techniques, with emphasis on the recent progress enabled by machine learning. It highlights representative datasets, modeling strategies, evaluation benchmarks, and the integration of structural and biochemical features. The aim is to provide researchers with a structured and up-to-date reference, guiding the development of more robust, generalizable, and interpretable models for predicting protein stability changes upon mutation. As the field moves forward, the synergy between data-driven AI methods and domain-specific biological knowledge will be key to unlocking deeper understanding and broader applications of protein engineering.
4.Analysis of risk prevention behaviors and influencing factors of HIV infection among young students with MSM
Chinese Journal of School Health 2025;46(9):1323-1327
Objective:
To explore the potential category patterns of risk prevention and control behaviors of HIV infection among young students who have sex with men (MSM) and their impact on HIV infection and late detection, aiming to optimize intervention strategies.
Methods:
From September 2017 to December 2024, a total of 1 637 MSM young students in Tianjin were recruited through both online and offline channels. Latent class analysis was applied to classify 11 HIV risk prevention and control behaviors [condom use during the most recent anal sex in the past 6 months, consistent condom use, use of water based lubricants, abstinence from recreational drugs, regular on site professional testing, fixed sexual partners, partner testing, awareness of partner s HIV testing results, testing before sexual activity, nucleic acid testing, and use of pre exposure prophylaxis (PrEP) or post exposure prophylaxis (PEP)]. Multivariate Logistic regression analyzed associations between demographic characteristics/intervention services factors and latent classes. Differences in HIV infection and late detection across behavior patterns were compared.
Results:
HIV risk prevention and control behaviors among MSM students were classified into three latent classes:condom dependent group (38.42%), low prevention group (27.73%), and comprehensive prevention group (33.85%). Students who received condom promotion/testing services were more likely to belong to the comprehensive prevention group ( OR =5.58), while those who received peer education were less likely to the comprehensive prevention group ( OR =0.43) (both P <0.01). Among the MSM student population, the HIV infection rate was 4.83%, with 2.26% of cases detected late. The HIV infection rate (1.45%) and late detection proportion (0.82%) in the comprehensive prevention group were lower than those in the low prevention group (7.89% and 3.83%, respectively) ( χ 2=16.20, 7.31, both P <0.01).
Conclusions
HIV risk prevention and control behaviors among MSM young students exhibit significant heterogeneity. Comprehensive prevention strategies can effectively reduce HIV infection and late detection risks. It is necessary to optimize peer education content and improve the accessibility of diversified prevention measures such as PrEP/PEP to enhance HIV prevention and control.
5.Application and data analysis of the cardio-cerebrovascular events monitoring system in Yichang
Zhengchao FANG ; Jiajuan YANG ; Chi HU ; Chan WU ; Yaling DENG ; Zhiying YU ; Jie ZHU ; Ling ZHANG
Journal of Public Health and Preventive Medicine 2025;36(6):95-98
Objective To analyze the monitoring data of cardio-cerebrovascular diseases prevention and control system in Yichang in 2022, and to provide data support and experience for the precise prevention and treatment of cardio-cerebrovascular diseases. Methods Acute cardiovascular and cerebrovascular event data were collected from the Yichang Cardio-cerebrovascular Events Monitoring System from January 1, 2022 to December 31, 2022. Descriptive analysis was conducted for the data collected. Statistical analysis was performed using SPSS 20.0 software, and a chi-square test was used to analyze the count data. Results A total of 37,217 cases of cardio-cerebrovascular events were monitored in Yichang in 2022. The crude incidence and the standardized incidence were 983.84/100,000 and 541.55/100,000, respectively. The incidence in males was higher than females (554.93/100,000 vs 428.91/100,000,χ2 =464.52,P<0.05). The top three diseases were cerebral infarction, acute myocardial infarction, and cerebral hemorrhage. The incidence of events increased with age, and 79.80% of the cases were over 60 years old. The main onset time was from May to August. Conclusion The use of the cardio-cerebrovascular events monitoring system in Yichang and the implementation of “mandatory reporting card” monitoring can timely obtain the epidemic characteristics of the diseases, provide support for the precise formulation of prevention and control strategies and measures, reduce underreporting rates, and improve the monitoring system, which is worthy of reference and promotion.
6.Staged treatment for 32 cases of ankylosing spondylitis with different acupuncture methods based on jingjin theory.
Ganggang LIU ; Ke YANG ; Jun GENG ; Cuiping LUO ; Ya ZHAO ; Jie ZHAO ; Yiyan ZHOU ; Yu SUN
Chinese Acupuncture & Moxibustion 2025;45(2):156-158
OBJECTIVE:
To observe the clinical efficacy of different acupuncture methods in the treatment of ankylosing spondylitis (AS) at different stages.
METHODS:
Thirty-two patients with AS were treated with acupuncture at ashi points (tendon focus), triple acupuncture method was used in the acute phase, accompanied needling method was used in the remission phase, once a day, 5 times a week for 4 weeks. The Bath ankylosing spondylitis disease activity index (BASDAI) score, Bath ankylosing spondylitis functional index (BASFI) score and TCM symptom grade quantitative standard score before and after treatment were observed, and the efficacy was evaluated.
RESULTS:
After treatment, the BASDAI score, BASFI score and TCM symptom grade quantitative standard score were reduced compared with those before treatment (P<0.01), the effective rate was 96.9% (31/32).
CONCLUSION
Different acupuncture methods in the treatment of AS at different stages could improve signs and symptoms, such as pain, activity limitation and morning stiffness.
Humans
;
Spondylitis, Ankylosing/physiopathology*
;
Acupuncture Therapy/methods*
;
Male
;
Female
;
Adult
;
Middle Aged
;
Young Adult
;
Acupuncture Points
;
Treatment Outcome
7.Effect of miR-129-3p mimetic on bone loss in tail-suspended mice
Yi WU ; Zi-dong AN ; Yong-jie PANG ; Li-qiang WANG ; Xin-yang WANG ; Yu-hai GAO ; Xue-yan LI ; Ke-ming CHEN
Chinese Pharmacological Bulletin 2025;41(4):703-709
Aim To study whether intravenous injec-tion of miR-129-3p mimetic(agomir)can resist bone loss caused by hind limb disuse,and to provide new i-deas for preventing bone loss in microgravity environ-ment.Methods Forty-eight C57BL/6J male mice were randomly divided into the control group(CON),tail suspension model group(TS),tail suspension+miR-129-3p agomir administration group(miRNA)and tail suspension+miR-129-3p negative sequence agomir control group(NC).The miRNA group was given 4 mg·kg-1 miR-129-3p agomir by intravenous injection into the medial canthus twice a week.The NC agomir group were consistent with those in the miR-129-3p agomir group,and the CON and TS groups were given only equal volumes of normal saline.After four weeks,all mice were sacrificed and samples were collected.Micro-CT scan of femur,three-point femur bending test,serum bone metabolism index detection,oxidative stress index detection and osteogenesis-related protein expression analysis in bone tissue were per-formed.Results After four weeks,the number of tra-becular bone in the TS group was significantly re-duced,and Tb.BMD,Tb.Th,Tb.N,Tb.BS/TV and Tb.BV/TV were significantly lower than those in the CON group(P<0.01).While Tb.Sp TS group was significantly higher than the CON group(P<0.05),the maximum load and flexural strength of the femur significantly decreased(P<0.01),the content of ser-um bone formation index PINP was significantly lower than that of the CON group(P<0.01),and the con-tent of bone resorption index CTX-I was significantly higher than that of the CON group(P<0.01),the content of serum oxidative damage indexes 8-iso-PGF2α and 8-OHdG significantly increased(P<0.01),and the expression of osteogenesis-related pro-teins in bone tissue markedly decreased(P<0.01).However,the increase or decrease of all indexes in miRNA group was significantly lower than that in TS group.Conclusions miR-129-3p mimetic can signifi-cantly reduce bone loss caused by hind limb disuse.This experiment provides a new idea and method for preventing bone loss in microgravity environment.
8.Effects of Rutong Ruanjian Tablets on angiogenesis in a rat model of preneoplastic breast cancer of Liver-Qi Stagnation and Blood Stasis Pattern via DLL4/Notch1/Hes1 pathway
Hua YANG ; Jun-yao LONG ; Jie GONG ; Bing-bing LU ; Xi ZOU ; Yu-rong WU ; Li-fang LIU ; Hui LIU ; Qi-hua CHEN
Chinese Traditional Patent Medicine 2025;47(3):774-781
AIM To investigate the effects of Rutong Ruanjian Tablets on angiogenesis in cancer tissues of rats with preneoplastic breast cancer(PBC).METHODS 60 female SD rats were randomly divided into a blank group of 10 rats and a model group of 50 rats for the establishment of the PBC models of Liver-Qi Stagnation and Blood Stasis Pattern with 9 weeks of oral administration of 7,12-dimethylbenz[a]anthracene(DMBA)and cervical ligation.After successful modeling,the rats were randomly divided into the model group,the tamoxifen group(3.2 mg/kg),the Rutong Ruanjian Tablets group(128 mg/kg),the 3,5-difluorobenzoyl group(DAPT,5 mg/kg),and the Rutong Ruanjian Tablets(128 mg/kg via gavage)+DAPT(5 mg/kg intraperitoneal injection)group,for 1 month corresponding drug administration,with 10 rats in each group.Then the rats had their cancer progression and syndrome scores observed;their angiogenesis evaluated by assessment of microvascular density(MVD);their vascular endothelial growth factor(VEGF)expression assessed by immunohistochemistry;and their mRNA and protein expressions of proteins related to the DLL4/Notch1/Hes1 pathway measured using RT-qPCR,immunohistochemistry and Western blot.RESULTS During carcinogenesis of rats induced by DMBA,there was gradual disappearance of E-cadherin expression and consistency of HE staining result with the PBC progression confirming the success of the modeling.Compared with the blank group,the model group showed increased MVD values,mRNA expression of Notch1 and Hes1,and protein expressions of VEGF,DLL4,Notch1 and Hes1(P<0.05,P<0.01).Compared with the model group,the Rutong Ruanjian Tablets group exhibited reduced MVD values,mRNA expression of Notch1 and Hes1,and protein expressions of VEGF,DLL4,Notch1 and Hes1(P<0.05,P<0.01).The Rutong Ruanjian Tablets+DAPT group showed reduced mRNA expression of Notch1 and Hes1,and protein expressions of DLL4,Notch1 and Hes1 compared to the Rutong Ruanjian Tablets group(P<0.05,P<0.01).CONCLUSION Rutong Ruanjian Tablets can inhibit angiogenesis and attenuate cancer progression in PBC rats of Liver-Qi Stagnation and Blood Stasis Pattern,and the mechanism may lie in the downregulation of DLL4/Notch1/Hes1 signaling pathway related proteins.
9.The value of total volume response and total mass response in the therapeutic evaluation of lung metastasis of hepatocarcinoma
Jun-cheng WAN ; Cai-hong YU ; Chang-yu LI ; Yong-jie ZHOU ; Wei ZHANG ; Jian-hua WANG ; Zhi-ping YAN ; Guo-wei YANG ; Zhuo-yang FAN ; Xu-dong QU
Fudan University Journal of Medical Sciences 2025;52(2):201-208,231
Objective To analyze the correlation between lesion volume,lesion mass,and maximum lesion diameter in the assessment of advanced hepatocarcinoma with lung metastasis,and to evaluate the application value of total volume response and total mass response of lung metastatic lesions in efficacy assessment.Methods A retrospective analysis was conducted on the CT imaging data of 20 patients clinically confirmed with hepatocarcinoma and lung metastases,followed by subsequent follow-up to monitor their survival outcomes.Volume measurement software was used to measure the volume of lesions before and after treatment.We recored lesion diameter,volume measurements and CT values,calculated the mass of the lesions.The correlation between lesion volume,mass and diameter was analyzed,as well as the correlation between the change rates of volume,mass and lesion diameter.Additionally,the total volume and total mass of all lesions were calculated.The correlation between the change rates of total volume/total mass and the change rate of pulmonary lesion diameter under the RECIST 1.1 criteria,as well as the correlation with changes in patients'tumor markers,were analyzed.Furthermore,the overall volume response and overall mass response of lesions were evaluated based on changes in total volume and total mass,and their consistencies with the RECIST 1.1 criteria for efficacy evaluation were analyzed.Finally,univariate Cox regression analysis was performed to explore the association between these variables and patient survival outcomes.Results There was strong correlation between lesion volume,mass and tumor diameter(r=0.771,0.775),between the rate of change in mass and the rate of change in lesion diameter(r=0.846),and between the rates of change in total volume/total mass and the rate of change in pulmonary lesion diameter under the RECIST 1.1 criteria(r=0.800,0.896).The correlation between the rates of change in total volume/total mass and patients'tumor markers was not statistically significant.There was moderate correlation between the rate of change in volume and the rate of change in lesion diameter(r=0.692).The evaluation results of total volume response and total mass response for pulmonary lesions in advanced hepatocarcinoma with lung metastasis were generally consistent with the RECIST 1.1 criteria(Kappa=0.486,0.426).Univariate Cox regression analysis revealed that total lesion volume(P=0.047)and total lesion mass(P=0.049)were independent prognostic factors for survival outcomes.Conclusion Lesion volume,mass,and diameter,as well as their respective change rates,were found to be interrelated.Furthermore,total lesion volume and total lesion mass were identified as independent prognostic factors for survival outcomes.The total volume response and total mass response are promising evaluation methods in evaluating the efficacy of lung metastasis of hepatocarcinoma,which are different from the RECIST 1.1 evaluation criteria.
10.Effect of TPCK trypsin on proliferation of porcine epidemic diarrhea virus in Vero cells
Damei ZHANG ; Liu YANG ; Guangliang GAO ; Shaomei LI ; Jie LUO ; Lizhi FU ; Yuandi YU ; Guoyang XU
Chinese Journal of Veterinary Science 2025;45(5):919-925
The purpose of this study is to investigate the effect of TPCK trypsin on the proliferation pattern of porcine epidemic diarrhea virus in Vero cells.The TPCK trypsin and conventional tryp-sin were added for virus proliferation,and RT-qPCR technology was used to analyze the changes in virus adsorption and invasion in Vero cells.The replication ability of porcine epidemic diarrhea vi-rus in Vero cells was explored through growth curve drawing,IFA identification,and cell activity detection.The results showed that the optimal concentrations of TPCK trypsin and conventional trypsin were 1 mg/L and 6 mg/L,respectively.The virus showed a decreasing trend with the pro-longation of TPCK trypsin and conventional trypsin pretreatment time.Adding different pancreatic enzymes during the virus proliferation process did not promote the virus invasion in Vero cells.Af-ter 4 h of invasion,the virus particles of each group gradually increased.By plotting the growth curve,it was found that the virus content of the TPCK trypsin group reached its highest level at 24 h(lgTCID50=(6.30±0.14)/0.1 mL),followed by a decreasing trend at 36 h,and the fluorescence intensity produced at 24 h was higher than that of conventional trypsin.In summary,TPCK trypsin has a better promoting effect on the proliferation of porcine epidemic diarrhea virus in Vero cells.It provided theoretical basis for further research on the mechanism of TPCK trypsin affecting porcine epidemic diarrhea virus proliferation,and also provided data support for the isola-tion and purification of porcine epidemic diarrhea virus epidemic strains.


Result Analysis
Print
Save
E-mail