1.Explainable Machine Learning Model for Predicting Prognosis in Patients with Malignant Tumors Complicated by Acute Respiratory Failure: Based on the eICU Collaborative Research Database in the United States
Zihan NAN ; Linan HAN ; Suwei LI ; Ziyi ZHU ; Qinqin ZHU ; Yan DUAN ; Xiaoting WANG ; Lixia LIU
Medical Journal of Peking Union Medical College Hospital 2026;17(1):98-108
To develop and validate a model for predicting intensive care unit (ICU) mortality risk in patients with malignant tumors complicated by acute respiratory failure (ARF) based on an explainable machine learning framework. Clinical data of patients with malignant tumors and ARF were extracted from the eICU Collaborative Research Database in the United States, including demographic characteristics, comorbidities, vital signs, laboratory test indicators, and major interventions within the first 24 hours after ICU admission.The study outcome was ICU death.Enrolled patients were randomly divided into a training set and a validation set at a ratio of 7:3.Predictor variables were selected using least absolute shrinkage and selection operator (LASSO) regression.Five machine learning algorithms-extreme gradient boosting (XGBoost), support vector machine (SVM), Logistic regression, multilayer perceptron (MLP), and C5.0 Decision Tree-were employed to construct predictive models.Model performance was evaluated based on the area under the receiver operating characteristic curve (AUC), accuracy, sensitivity, and other metrics.The optimal model was further interpreted using the Shapley additive explanations (SHAP) algorithm. A total of 3196 patients with malignant tumors complicated by ARF were included.The training set comprised 2, 261 patients and the validation set 935 patients; 683 patients died during ICU stay, while 2513 survived.LASSO regression ultimately selected 12 variables closely associated with patient ICU outcomes, including sepsis comorbidity, use of vasoactive drugs, and within the first 24 hours after ICU admission: minimum mean arterial pressure, maximum heart rate, maximum respiratory rate, minimum oxygen saturation, minimum serum bicarbonate, minimum blood urea nitrogen, maximum white blood cell count, maximum mean corpuscular volume, maximum serum potassium, and maximum blood glucose.After model evaluation, the XGBoost model demonstrated the best performance.The AUCs for predicting ICU mortality risk in the training and validation sets were 0.940 and 0.763, respectively; accuracy was 88.3% and 81.2%;sensitivity was 98.5% and 95.9%.Its predictive performance also remained optimal in sensitivity analyses.SHAP analysis indicated that the top five variables contributing to the model's predictions were minimum oxygen saturation, minimum serum bicarbonate, minimum mean arterial pressure, use of vasoactive drugs, and maximum white blood cell count. This study successfully developed a mortality risk prediction model for ICU patients with malignant tumors complicated by ARF based on a large-scale dataset and performed explainability analysis.The model aids clinicians in early identification of high-risk patients and implementing individualized interventions.
2.Targeting GYS1: From Metabolic Regulatory Mechanisms to Precision Therapeutic Strategies
Jia-Nan ZHAO ; Yu-Xuan LI ; Jie ZHU ; Hong LI ; Xiao-Feng JIN
Progress in Biochemistry and Biophysics 2026;53(7):1807-1825
Glycogen synthase 1 (GYS1) is the rate-limiting enzyme responsible for glycogen synthesis in skeletal muscle, heart, brain, and other extrahepatic tissues, playing a central role in systemic energy homeostasis. The human GYS1 gene maps to chromosome 19q13.33, comprises 16 exons, and encodes a 737-amino-acid polypeptide that is highly conserved across mammals. GYS1 activity is subject to multilayered and precisely coordinated regulation. At the transcriptional level, the GYS1 promoter contains a hypoxia response element (HRE) that mediates HIF-1α-dependent induction under low-oxygen conditions, as well as a muscle-specific enhancer harboring MEF2 and MyoD binding sites that confers tissue-restricted expression. At the post-translational level, a hierarchical phosphorylation cascade serves as the primary activity switch: glycogen synthase kinase 3β (GSK3β) sequentially phosphorylates four C-terminal serine residues following casein kinase II priming, while protein kinase A (PKA) and AMP-activated protein kinase (AMPK) provide parallel inhibitory inputs at both N- and C-terminal sites. Dephosphorylation and reactivation are mediated by protein phosphatase 1 (PP1) through tissue-specific glycogen-targeting regulatory subunits such as PPP1R3A and PPP1R3B, which anchor PP1 to glycogen particles and direct its activity toward GYS1. The allosteric activator glucose-6-phosphate (G6P) binds at the dimer interface, simultaneously enhancing catalytic efficiency and promoting dephosphorylation susceptibility, thereby establishing a feed-forward activation loop that couples substrate availability to glycogen synthesis. Beyond phosphorylation, GYS1 is regulated by ubiquitination (mediated by the E3 ligase PJA1), acetylation, O-linked β-N-acetylglucosamine (O-GlcNAc) modification, and SUMOylation, which collectively modulate protein stability, subcellular localization, and protein-protein interactions. Epigenetic mechanisms, including CpG island methylation and histone acetylation dynamics, govern chromatin accessibility at the GYS1 locus, while muscle-specific microRNAs such as miR-1 and miR-206 fine-tune GYS1 expression at the post-transcriptional level. Dysregulation of GYS1 has been identified as a central pathogenic driver in a spectrum of human diseases. In inherited glycogen storage disorders—including Lafora disease, adult polyglucosan body disease (APBD), and Pompe disease—loss of upstream regulatory control leads to GYS1 hyperactivation and the accumulation of structurally abnormal or excessive glycogen, resulting in progressive neurodegeneration, myopathy, and multiorgan dysfunction. In type 2 diabetes mellitus (T2DM), impaired insulin signaling through the PI3K-AKT-GSK3β axis maintains GYS1 in a hyperphosphorylated inactive state in skeletal muscle, compromising postprandial glucose disposal and exacerbating hyperglycemia. In oncology, GYS1 exhibits context-dependent roles across multiple cancer types. In hepatocellular carcinoma, FMO2+ cancer-associated fibroblasts stabilize GYS1 by competitively inhibiting PJA1-mediated ubiquitination, and stabilized GYS1 subsequently activates NF‑κB/CCL19 signaling to promote tertiary lymphoid structure formation and enhance anti-PD-1 immunotherapy responsiveness. In clear cell renal cell carcinoma, GYS1 promotes tumor progression through non-canonical NF‑κB pathway activation via the scaffold protein RPS27A. In triple-negative breast cancer, GYS1 has been identified as a trigger of disulfidptosis and an activator of NF-κB signaling through non-enzymatic facilitation of IκBα degradation. In colorectal cancer, mitochondrial fission deficiency drives AMPK-dependent GYS1 upregulation and glycogen accumulation as a compensatory survival mechanism, while in cervical cancer, GYS1-maintained glycogen reserves fuel the pentose phosphate pathway to generate NADPH for ROS clearance, thereby conferring cisplatin resistance in cancer stem cells. Therapeutic strategies targeting GYS1 have gained substantial momentum across these disease contexts. For glycogen storage disorders, antisense oligonucleotides, small interfering RNAs (e.g., ABX1100), and small-molecule inhibitors (e.g., MZ-101) have demonstrated preclinical and early clinical efficacy in reducing pathological glycogen accumulation. For T2DM, pharmacological activation of GYS1 through GSK3β inhibition or enhancement of PP1-mediated dephosphorylation is being explored to restore insulin-stimulated glycogen synthesis. In cancer, GYS1-directed interventions—including targeted silencing to sensitize tumors to chemotherapy and immune microenvironment modulation to enhance immunotherapy—represent emerging precision oncology approaches. This review provides a comprehensive and integrated account of GYS1 gene structure, tissue-specific distribution, regulatory networks, and pathogenic roles in metabolic disorders and malignancies, with the aim of establishing a theoretical framework for the development of GYS1-targeted precision therapies.
3.Recommendations for the clinical use of anti-amyloid-β monoclonal antibody for Alzheimer's disease(2025)
Nan ZHI ; Jinwen XIAO ; Rujing REN ; Binyin LI ; Jintao WANG ; Jieli GENG ; Wenwei CAO ; Yaying SONG ; Hualong WANG ; Shuguang CHU ; Guoping PENG ; Jun LIU ; Xiaoyun LIU ; Fang YUAN ; Wen WANG ; Ronghua DOU ; Xia LI ; Ling YUE ; Wenshi WEI ; Xiaoling PAN ; Xiangyang ZHU ; Dian HE ; Weinü FAN ; Jingping SHI ; Nan ZHANG ; Hui ZHAO ; Qin CHEN ; Cuibai WEI ; Xiaochun CHEN ; Gang WANG
Journal of Chongqing Medical University 2025;50(9):1133-1140
In recent years,significant breakthroughs have been achieved in the immunotherapy for Alzheimer's disease.In line with global advancements,two anti-amyloid-β monoclonal antibodies have been approved and successfully launched in China for clinical use.Lecanemab and Donanemab were officially used in June 2024 and April 2025 in China,respectively.In order to standardize the rational and safe application of anti-amyloid-β monoclonal antibodies for Alzheimer's disease in China,this article integrates recom-mendations from the clinical trials and real-world experience from the author's team and domestic peers to further update the recom-mendations for the clinical use of anti-amyloid-β monoclonal antibody based on the 2024 version.It includes indications for therapy,pre-treatment evaluation and preparation,administration protocols and safety measures during treatment,and post-treatment monitor-ing strategies.
4.Simultaneous Determination of Perfluorooctanoic Acid and Perfluorooctane Sulfonate Isomers in Seawater by Online Solid Phase Extraction Coupled with Liquid Chromatography-Tandem Mass Spectrometry
Jun-Hui CHEN ; Nan SHEN ; Tong-Zhu HAN ; Xiu-Ping HE ; Xian-Guo LI
Chinese Journal of Analytical Chemistry 2025;53(7):1146-1157
A new method was developed for simultaneous and efficient determination of linear perfluorooctanoic acid(n-PFOA)and linear perfluorooctane sulfonate(n-PFOS),and their typical branched isomers in seawater by online solid phase extraction-liquid chromatography-tandem mass spectrometry(Online SPE-LC-MS/MS).Only centrifugation of the seawater sample was required to remove the particulate matter,and then the seawater sample was directly injected and analyzed by online SPE-LC-MS/MS.An Eclipse Plus-C18 guard column was selected as SPE column for online enrichment of linear and branched isomers,and a F5 PFP column(150 mm×2.1 mm,2.7 μm)was used as the analytical column.Under the optimized experimental conditions,the separation and detection of all PFOA and PFOS linear and branched isomers could be completed within 20 min.The spiked recoveries of various target compounds ranged from 82.9%to 107.7%with detection limits and limits of quantification of 0.10-1.05 ng/L and 0.30-2.11 ng/L,respectively.The method was characterized by good precision(RSD≤9.10%)and linearity(R2≥0.990).Subsequently,linear and branched isomers of PFOA and PFOS in surface and bottom seawater samples collected from the Laizhou Bay of China were determined.The results showed that the detection rate of all the four branched PFOA isomers were 100%,with the highest average concentration of 25.85 ng/L found for 6m-PFOA,which accounted for 11.79%of the∑PFOA.For the five branched isomers of PFOS,the highest detection rate of 90.84%was found for 5m-PFOS.The highest average concentration of 0.64 ng/L was observed for 3m-PFOS,accounting for 19.88%of ∑PFOS.The proposed method provided an effective detection tool for qualitative and quantitative detection of PFOA and PFOS isomers in the marine aquatic environment.
5.Prevalence and influencing factors of chronic somatic comorbidities in patients with schizophrenia in Shanghai communities
Wei-bo ZHANG ; Jiang-nan LI ; Yan-li LIU ; Yi-zhou JIANG ; Yi ZHU ; Na WANG ; Jun CAI
Fudan University Journal of Medical Sciences 2025;52(4):484-491,499
Objective To investigate the prevalence of chronic somatic comorbidities in schizophrenic patients in Shanghai communities,and to explore the factors influencing comorbidities.Methods Based on Shanghai community-based severe mental disorders cohort,5 422 patients with schizophrenia(SCZ)were included in the study.12 common chronic somatic diseases,defined by patients'self-report,were selected to analyze the prevalence of comorbidity,and Logistic regression model was used to analyze the factors influencing the number of somatic comorbidities.Results The total prevalence of somatic comorbidity was 37.0%in 5 422 patients with SCZ,with the highest prevalence of hypertension(22.6%)and diabetes mellitus(13.1%)among 12 somatic diseases selected.Older age was the main factor associated with chronic somatic comorbidities in community schizophrenic patients.The risk of 1-2 comorbidities in patients aged≥60 years was 3.34(95%CI:2.74-4.07)times higher than those aged<45 years,while the risk of≥3 comorbidities was 3.27(95%CI:2.11-5.09)times higher,correspondingly.Female gender,marriage,smoking,and longer duration of illness were also risk factors for comorbidity.Women after menopause had higher risk of comorbidity than perimenopausal women.Conclusion Cardiovascular and metabolic diseases were common somatic comorbidities among schizophrenic patients in Shanghai communities.Older age,female gender,marriage,smoking,and longer duration of illness were risk factors for increasing number of comorbidities.
6.Clinical and genetic characteristics of SCN2A gene related developmental delay
Jialu GU ; Shaofang SHANGGUAN ; Jianhong WANG ; Jiayi LI ; Hua XIE ; Xia QU ; Nan PENG ; Xi WANG ; Qi XU ; Yike ZHU ; Xinghui LI ; Xuefeng SUN ; Xiaoli CHEN ; Lin WANG
Chinese Journal of Preventive Medicine 2025;59(5):667-676
Objective:To explore the genotype and the clinical phenotype of SCN2A-related developmental delay in children. Methods:A case series study was adopted. Collect clinical data from 10 cases of children with SCN2A gene variants diagnosed with global developmental delay/intellectual disability who were admitted to the Children′s Hospital between July 2019 and March 2023. Summarize the clinical phenotype and genotype based on clinical data such as general information, clinical manifestations, imaging examinations, laboratory tests, genetic testing results, and comprehensive pediatric neuropsychological development assessment. Results:A total of 10 patients were recruited, including 7 males and 3 females, with an age range of 27 days to 5 years and 9 months. 9 patients underwent children′s neuropsychological and behavioral assessments, and the results were consistent with global developmental delay, including 2 mild cases, 4 moderate cases, and 3 severe cases. 3 cases had autism spectrum disorder, and 2 cases had epilepsy. 6 patients underwent complete head MRI examination, and 4 of them showed abnormalities, including delayed myelination, widening of the local extra brain space in the frontal lobe, and abnormal frontal lobe morphology. All 10 cases had point variants. Among them, 9 cases are de novo and 1 case is maternal inheritance. Out of 10 cases, there were 5 cases with copy number variations, but all of them were of unknown significance. Among the 10 variants, 8 have been reported and 2 have not been reported, namely c.4145A>T(p.N1382I) and c.4937T>A(p.I1646N). In this study, 4 out of 10 patients with SCN2A variants had variation sites located in the S4 segment of domain which constitute Nav1.2, the sodium ion channel encoded by SCN2A. The developmental quotient level was lower when the variation sites were located in the S4 segment of domain, and the difference was statistically significant ( t=-3.101, P=0.017), indicating that the severity of developmental delay may be related to the localization of amino acids corresponding to variant sites within the protein domain. Conclusion:SCN2A mutations are strongly associated with diverse neurodevelopmental disorders. In this study, the phenotypic spectrum of SCN2A variants encompassed epilepsy, global developmental delay, and autism spectrum disorder. Affected individuals exhibited early-onset developmental delays, predominantly moderate to severe in severity. Voltage-sensing domain dysfunction in sodium channels may constitute a critical pathomechanism underlying neurodevelopmental impairments. Further electrophysiological characterization and molecular mechanistic studies are warranted todelineate the genotype-phenotype correlations between specific variant loci and clinical severity.
7.Real-world efficacy and safety of azvudine in hospitalized older patients with COVID-19 during the omicron wave in China: A retrospective cohort study.
Yuanchao ZHU ; Fei ZHAO ; Yubing ZHU ; Xingang LI ; Deshi DONG ; Bolin ZHU ; Jianchun LI ; Xin HU ; Zinan ZHAO ; Wenfeng XU ; Yang JV ; Dandan WANG ; Yingming ZHENG ; Yiwen DONG ; Lu LI ; Shilei YANG ; Zhiyuan TENG ; Ling LU ; Jingwei ZHU ; Linzhe DU ; Yunxin LIU ; Lechuan JIA ; Qiujv ZHANG ; Hui MA ; Ana ZHAO ; Hongliu JIANG ; Xin XU ; Jinli WANG ; Xuping QIAN ; Wei ZHANG ; Tingting ZHENG ; Chunxia YANG ; Xuguang CHEN ; Kun LIU ; Huanhuan JIANG ; Dongxiang QU ; Jia SONG ; Hua CHENG ; Wenfang SUN ; Hanqiu ZHAN ; Xiao LI ; Yafeng WANG ; Aixia WANG ; Li LIU ; Lihua YANG ; Nan ZHANG ; Shumin CHEN ; Jingjing MA ; Wei LIU ; Xiaoxiang DU ; Meiqin ZHENG ; Liyan WAN ; Guangqing DU ; Hangmei LIU ; Pengfei JIN
Acta Pharmaceutica Sinica B 2025;15(1):123-132
Debates persist regarding the efficacy and safety of azvudine, particularly its real-world outcomes. This study involved patients aged ≥60 years who were admitted to 25 hospitals in mainland China with confirmed SARS-CoV-2 infection between December 1, 2022, and February 28, 2023. Efficacy outcomes were all-cause mortality during hospitalization, the proportion of patients discharged with recovery, time to nucleic acid-negative conversion (T NANC), time to symptom improvement (T SI), and time of hospital stay (T HS). Safety was also assessed. Among the 5884 participants identified, 1999 received azvudine, and 1999 matched controls were included after exclusion and propensity score matching. Azvudine recipients exhibited lower all-cause mortality compared with controls in the overall population (13.3% vs. 17.1%, RR, 0.78; 95% CI, 0.67-0.90; P = 0.001) and in the severe subgroup (25.7% vs. 33.7%; RR, 0.76; 95% CI, 0.66-0.88; P < 0.001). A higher proportion of patients discharged with recovery, and a shorter T NANC were associated with azvudine recipients, especially in the severe subgroup. The incidence of adverse events in azvudine recipients was comparable to that in the control group (2.3% vs. 1.7%, P = 0.170). In conclusion, azvudine showed efficacy and safety in older patients hospitalized with COVID-19 during the SARS-CoV-2 omicron wave in China.
8.DNA methylation age prediction model based on ovarian granulosa cells
Peng LIU ; Bowen ZHU ; Yuheng LI ; Liang WANG ; Nan LIU ; Ningxia SUN
Chinese Journal of Reproduction and Contraception 2025;45(5):442-447
Objective:To develop a novel ovarian biological age prediction model based on DNA methylation for ovarian aging assessment.Methods:A prospective cohort study method was used. From March 2024 to January 2025, we collected 96 ovarian granulosa cell samples of infertility patients due to fallopian tube factors or male factors from the Reproductive Medicine Center of the Second Affiliated Hospital of Naval Medical University. Then we analyzed DNA methylation patterns across five age-associated gene regions ( ELOVL2, miR29B2C, TRIM59, KLF14 and FHL2) in a discovery cohort comprising 63 human ovarian granulosa cell samples. Targeted bisulfite sequencing was performed through PCR amplification followed by next-generation DNA sequencing. Leveraging elastic net regression analysis, we developed a predictive model incorporating 29 methylation sites that demonstrated strong age correlation. The model was subsequently validated using an independent cohort comprising 33 human ovarian granulosa cell samples. Results:The DNA methylation age prediction model based on ovarian granulosa cells showed the following results in the discovery cohort as follows: median absolute error (MAE) was 2.534 ( R=0.742, P<0.001). In the independent validation cohort, MAE was 3.019 ( R=0.729, P<0.001). Conclusion:In this study, we utilized human ovarian granulosa cells as experimental samples to develop a novel DNA methylation-based model for predicting ovarian biological age. By integrating multiple methylation sites across five age-related gene regions, this model serves as a robust indicator of ovarian aging status.
9.Analysis of the viral molecular characteristics in a dengue fever outbreak in Jiangsu province in 2023
Yuanfang QIN ; Nan ZHANG ; Qian ZHEN ; Zhifeng LI ; Hao JU ; Liguo ZHU ; Jianguang FU
Chinese Journal of Experimental and Clinical Virology 2025;39(1):81-85
Objective:To analyze the molecular characteristics of the virus in a local outbreak of dengue fever in Jiangsu province in 2023, and to provide a basis for the prevention and control of the outbreak.Methods:Serum samples were collected from suspected dengue patients in the acute phase of the outbreak for virus detection and serotyping by real-time fluorescence quantitative RT-PCR (RT-qPCR). Positive specimens were amplified with full-length genomic fragments and subjected to second-generation sequencing and related evolutionary analyses.Results:Four confirmed cases of dengue were found in Changzhou city, Jiangsu province, from October 18 to 21, 2023, with epidemiological association between the cases, which was recognized as a dengue outbreak. The serum RT-qPCR result of the four cases were all dengue type 1, and the whole genome sequences of three of the cases were obtained. The evolutionary tree of the E gene and the whole genome showed that the three sequences were located in the 3rd branch of the 1-I genotype, which is similar to the genotype 1-I. The genome-wide sequences of the E gene and the genome-wide evolution tree showed that the three sequences were located in the 3rd branch of the 1-I genotype, which is similar to the genome-wide genotype 1-Ⅰ. The E gene and the genome-wide evolutionary tree showed that all three sequences were located on branch 3 of genotype 1-Ⅰ, with high sequence similarity to the dengue virus epidemic strains in Guangdong and Yunnan provinces in 2023. Amino acid variant site analysis showed that there were 16 branch-specific amino acid site changes in the sequences of the three cases, among which the structural proteins, C protein and prM protein, had one variant site each, E protein had two, and the non-structural proteins had the largest number of NS5 variant sites (9).Conclusions:The local outbreak in Jiangsu was caused by dengue fever type 1 virus, with high nucleotide sequence similarity to strains from other regions of China, and amino acid site alterations.
10.Clinical and genetic characteristics of SCN2A gene related developmental delay
Jialu GU ; Shaofang SHANGGUAN ; Jianhong WANG ; Jiayi LI ; Hua XIE ; Xia QU ; Nan PENG ; Xi WANG ; Qi XU ; Yike ZHU ; Xinghui LI ; Xuefeng SUN ; Xiaoli CHEN ; Lin WANG
Chinese Journal of Preventive Medicine 2025;59(5):667-676
Objective:To explore the genotype and the clinical phenotype of SCN2A-related developmental delay in children. Methods:A case series study was adopted. Collect clinical data from 10 cases of children with SCN2A gene variants diagnosed with global developmental delay/intellectual disability who were admitted to the Children′s Hospital between July 2019 and March 2023. Summarize the clinical phenotype and genotype based on clinical data such as general information, clinical manifestations, imaging examinations, laboratory tests, genetic testing results, and comprehensive pediatric neuropsychological development assessment. Results:A total of 10 patients were recruited, including 7 males and 3 females, with an age range of 27 days to 5 years and 9 months. 9 patients underwent children′s neuropsychological and behavioral assessments, and the results were consistent with global developmental delay, including 2 mild cases, 4 moderate cases, and 3 severe cases. 3 cases had autism spectrum disorder, and 2 cases had epilepsy. 6 patients underwent complete head MRI examination, and 4 of them showed abnormalities, including delayed myelination, widening of the local extra brain space in the frontal lobe, and abnormal frontal lobe morphology. All 10 cases had point variants. Among them, 9 cases are de novo and 1 case is maternal inheritance. Out of 10 cases, there were 5 cases with copy number variations, but all of them were of unknown significance. Among the 10 variants, 8 have been reported and 2 have not been reported, namely c.4145A>T(p.N1382I) and c.4937T>A(p.I1646N). In this study, 4 out of 10 patients with SCN2A variants had variation sites located in the S4 segment of domain which constitute Nav1.2, the sodium ion channel encoded by SCN2A. The developmental quotient level was lower when the variation sites were located in the S4 segment of domain, and the difference was statistically significant ( t=-3.101, P=0.017), indicating that the severity of developmental delay may be related to the localization of amino acids corresponding to variant sites within the protein domain. Conclusion:SCN2A mutations are strongly associated with diverse neurodevelopmental disorders. In this study, the phenotypic spectrum of SCN2A variants encompassed epilepsy, global developmental delay, and autism spectrum disorder. Affected individuals exhibited early-onset developmental delays, predominantly moderate to severe in severity. Voltage-sensing domain dysfunction in sodium channels may constitute a critical pathomechanism underlying neurodevelopmental impairments. Further electrophysiological characterization and molecular mechanistic studies are warranted todelineate the genotype-phenotype correlations between specific variant loci and clinical severity.

Result Analysis
Print
Save
E-mail