1.Study on the epidemiological characteristics and influencing factors of long COVID among previously infected individuals in two communities in Shanghai
Junhong YUE ; Chen CHEN ; Qingqing JIA ; Xiaoxia LIU ; Huiting WANG ; Fei WU ; Yanlu YIN ; Jiajie ZANG ; Yanfei GUO ; Fan WU
Shanghai Journal of Preventive Medicine 2025;37(7):597-605
ObjectiveTo analyze the epidemiological characteristics of long COVID and to investigate its main influencing factors by examining individuals infected with SARS-CoV-2 between March and June 2022 in two communities in Shanghai, to lay the foundation for further research on the mechanism and clinical treatment of long COVID, and to provide the basis for the development of inexpensive, convenient, and feasible prevention and intervention strategies. MethodsA cross-sectional study was conducted, enrolling 6 410 individuals infected with SARS-CoV-2. Data were collected through a questionnaire survey. The incidence and common symptoms of long COVID were analyzed, along with their associations with demographic characteristics, medical history, and behavioral factors. A logistic regression model was used to identify the major factors associated with the development of long COVID symptoms. ResultsThe overall incidence rate of long COVID among the study population was 13.9%. The most commonly reported symptoms included fatigue (65.1%), attention disorders (23.1%), and cough (16.9%). The analysis showed that having underlying chronic diseases (OR=2.580, 95%CI: 2.165‒3.074), a history of allergies (OR=1.418, 95%CI: 1.003‒1.971), current smoking (OR=1.461, 95%CI: 1.013‒2.079), ever smoking (OR=2.462, 95%CI: 1.687‒3.551), a greater number of symptoms during the acute phase [1 symptom (OR=1.778, 95%CI: 1.459‒2.162), 2 symptoms (OR=2.749, 95%CI: 2.209‒3.409), ≥3 symptoms (OR=7.792, 95%CI: 6.333‒9.593)] and aggravated symptoms during the acute phase (OR=1.082, 95%CI: 1.070‒1.094) were factors associated with a higher risk of developing long COVID symptoms. Additionally, individuals who had consumed alcohol in the past year (OR=1.914, 95%CI: 1.344‒2.684) were more prone to objective long COVID symptoms. Among individuals under 50 years of age, females (OR=1.427, 95%CI: 1.052‒1.943) were more likely to develop objective long COVID symptoms. ConclusionThis study has identified the diversity of long COVID symptoms, which involve multiple organs and systems, including fatigue, attention disorders, cough, and joint pain. It has also revealed associations between long COVID and various demographic factors (e.g., age, gender), personal medical history (e.g., underlying chronic diseases, history of allergies), acute-phase characteristics (e.g., number and severity of symptoms), and behavioral factors (e.g., smoking, alcohol consumption). These findings highlight the need for further research and ongoing surveillance of long COVID and may inform the development of more targeted health management strategies for specific populations.
2.Characteristics of malignant tumor incidence and mortality in cancer registration areas of Zhejiang Province in 2021
DU Lingbin ; QIU Yu ; LI Huizhang ; LI Runhua ; ZHU Chen ; WANG Le ; QIU Yanfei
Journal of Preventive Medicine 2025;37(10):973-978
Objective:
To investigate the characteristics of malignant tumor incidence and mortality in cancer registration areas of Zhejiang Province in 2021.
Methods:
Based on the 2021 cancer registration data from 22 national cancer registries in Zhejiang Province, the crude incidence, crude mortality, and cumulative rate for 0 to 74 years were calculated. Age standardized was performed using the age composition of the standard population from the Fifth National Population Census in 2000 and Segi's world standard population. The incidence and mortality characteristics of malignant tumor in different genders, urban/rural areas and ages were described. The order of crude incidence and mortality of malignant tumor were analyzed.
Results:
In 2021, there were 116 639 new malignant tumor cases in Zhejiang Province. The crude, Chinese population standardized, and world population-standardized incidences were 530.93/100 000, 304.83/100 000, and 288.20/100 000, respectively. The cumulative incidence for 0 to 74 years was 31.92%. There were 40 475 death cases. The crude, Chinese population-standardized, and world population-standardized mortalities were 184.24/100 000, 79.40/100 000, and 78.97/100 000, respectively. The cumulative mortality for 0 to 74 years was 8.37%. The Chinese population-standardized incidence for males and females were 286.34/100 000 and 323.45/100 000, respectively, and the Chinese population-standardized mortality were 106.25/100 000 and 54.17/100 000, respectively. The Chinese population-standardized incidence for urban and rural were 316.85/100 000 and 285.11/100 000, respectively, and the Chinese population-standardized mortality were 75.59/100 000 and 85.48/100 000, respectively. The crude incidence and crude mortality of malignant tumor both increased with age, peaking in the groups aged 80-<80 and ≥85 years at 1 845.06/100 000 and 1 656.88/100 000, respectively. The top ten malignant tumors with the highest crude incidence were, in descending order: lung cancer, thyroid cancer, breast cancer, prostate cancer, colorectal cancer, gastric cancer, liver cancer, cervical cancer, lymphoma, and brain tumors. These accounted for 80.34% of all new malignant tumor cases. The top ten malignant tumors with the highest crude mortality were, in descending order: lung cancer, liver cancer, colorectal cancer, gastric cancer, pancreatic cancer, prostate cancer, breast cancer, esophageal cancer, lymphoma, and gallbladder cancer. These accounted for 82.26% of all malignant tumor deaths.
Conclusions
In 2021, the cancer registration areas of Zhejiang Province were characterized by a relatively high malignant tumor incidence and a comparatively low mortality compared to national data. Males and the elderly emerged as key populations for targeted prevention and control. It is recommended to enhance screening, early diagnosis, and early treatment for lung cancer, thyroid cancer, prostate cancer, breast cancer, and gastric cancer malignancies.
3.Impact of palliative care on medication use and medical utilization in patients with advanced cancer.
Dingyi CHEN ; Haoxin DU ; Yichen ZHANG ; Yanfei WANG ; Wei LIU ; Yuanyuan JIAO ; Luwen SHI ; Xiaodong GUAN ; Xinpu LU
Journal of Peking University(Health Sciences) 2025;57(5):996-1001
OBJECTIVE:
To evaluate the effect of palliative care on drug use, medical service utilization and medical expenditure of patients with advanced cancer.
METHODS:
A cohort of patients including pal-liative care and standard care was constructed using the medical records of the patients in Peking University Cancer Hospital from 2018 to 2020, and coarsened exact matching was used to match the two groups of patients. The average monthly opioid consumption, hospitalization rate, intensive care unit (ICU) rate and operation rate, and the average monthly total cost were selected to evaluate drug use, medical service utilization and medical expenditure. Chi-square test and Wilcoxon signed rank test were used to compare the differences between the two groups before and after exposure and the change in the palliative care group. The net impact of palliative care on the patients was calculated using the difference-in-differences analysis.
RESULTS:
In this study, 180 patients in the palliative care group and 3 101 patients in the stan-dard care group were finally included in the matching, and the matching effect of the two groups was good (L1 < 0.1). Before and after exposure, the average monthly opioid consumption in the palliative care group was significantly higher than that in the standard care group (Before exposure: 0.3 DDD/person-month vs. 0.1 DDD/person-month, P < 0.01; After exposure: 0.7 DDD/person-month vs. 0.1 DDD/person-month, P < 0.01; DDD refers to defined daily dose), palliative care significantly increased the average monthly opioid consumption in the patients (0.3 DDD/person-month, P < 0.01). The hospitalization rate (48.9% vs. 74.3%, P < 0.01) and operation rate (3.9% vs. 8.8%, P < 0.01) of the patients in palliative care group were significantly lower than those in standard care group, and the ICU rate became similar between the two groups (1.1% vs. 1.6%, P=0.634). Palliative care significantly reduced the patients ' hospitalization rate (-25.6%, P < 0.01), ICU rate (-4.9%, P < 0.01) and operation rate (-14.5%, P < 0.01). Before and after exposure, the average monthly total costs of pal-liative care group were slightly higher than those of standard care group (Before exposure: 20 092.3 yuan vs. 19 132.8 yuan, P=0.725; After exposure: 9 719.8 yuan vs. 8 818.8 yuan, P=0.165). Palliative care increased the average monthly total cost by 2 208.8 yuan, but it was not statistically significant (P=0.316).
CONCLUSION
Palliative care can increase the opioid consumption in advanced cancer patients, reduce the rates of hospitalization, ICU and surgery, but has no significant effect on medical expenditure.
Humans
;
Palliative Care/economics*
;
Neoplasms/drug therapy*
;
Analgesics, Opioid/economics*
;
Male
;
Female
;
Middle Aged
;
Aged
;
Hospitalization/economics*
;
Intensive Care Units/statistics & numerical data*
;
Health Expenditures/statistics & numerical data*
;
Adult
;
Drug Utilization/statistics & numerical data*
;
Patient Acceptance of Health Care/statistics & numerical data*
4.AI-integrated IQPD framework of quality prediction and diagnostics in small-sample multi-unit pharmaceutical manufacturing: Advancing from experience-driven to data-driven manufacturing.
Kaiyi WANG ; Xinhai CHEN ; Nan LI ; Huimin FENG ; Xiaoyi LIU ; Yifei WANG ; Yanfei WU ; Yufeng GUO ; Shuoshuo XU ; Lu YAO ; Zhaohua ZHANG ; Jun JIA ; Zhishu TANG ; Zhisheng WU
Acta Pharmaceutica Sinica B 2025;15(8):4193-4209
The pharmaceutical industry faces challenges in quality digitization for complex multi-stage processes, especially in small-sample systems. Here, an intelligent quality prediction and diagnostic (IQPD) framework was developed and applied to Tong Ren Tang's Niuhuang Qingxin Pills, utilizing four years of data collected from four production units, covering the entire process from raw materials to finished products. In this framework, a novel path-enhanced double ensemble quality prediction model (PeDGAT) is proposed, which combines a graph attention network and path information to encode inter-unit long-range and sequential dependencies. Additionally, the double ensemble strategy enhances model stability in small samples. Compared to global traditional models, PeDGAT achieves state-of-the-art results, with an average improvement of 13.18% and 87.67% in prediction accuracy and stability on three indicators. Additionally, a more in-depth diagnostic model leveraging grey correlation analysis and expert knowledge reduces reliance on large samples, offering a panoramic view of attribute relationships across units and improving process transparency. Finally, the IQPD framework integrates into a Human-Cyber-Physical system, enabling faster decision-making and real-time quality adjustments for Tong Ren Tang's Niuhuang Qingxin Pills, a product with annual sales exceeding 100 million CNY. This facilitates the transition from experience-driven to data-driven manufacturing.
5.Role of O-sialoglycoprotein endopeptidase in hepatic ischemia-reperfusion injury in mice: relationship with oxidative stress
Tengjuan ZHANG ; Wanqing ZHOU ; Cheng CHEN ; Qian ZHANG ; Yanfei ZHAO ; Dehao HE ; Zhi YE ; Pingping XIA
Chinese Journal of Anesthesiology 2024;44(1):85-90
Objective:To evaluate the role of O-sialoglycoprotein endopeptidase (OSGEP) in hepatic ischemia-reperfusion injury (HIRI) and the relationship with oxidative stress in mice.Methods:Experiment Ⅰ Twenty-four SPF healthy male C57BL/6 mice, 12 wild-type and 12 OSGEP knockdown, aged 6-8 weeks, weighing 18-22 g, were divided into 4 groups ( n=6 each) by the random number table method: wild-type shamoperation group (Sham group), wild-type HIRI group (HIRI group), OSGEP knockdown+ sham operation group (Sham+ KD group) and OSGEP knockdown+ HIRI group (HIRI+ KD group). Ischemia-reperfusion model was prepared by blocking the hepatic artery and portal vein for 60 min followed by reperfusion in anesthetized animals, the blood vessels were only exposed without occlusion in Sham group and Sham+ KD group, and the blood vessels were clamped for 60 min followed by reperfusion in HIRI group and HIRI+ KD group. The mice were sacrificed after 6-h reperfusion to extract liver tissue samples for microscopic examination of histopathological changes (with an optical microscope after HE staining) which were evaluated using Suzuki score and for determination of the serum concentrations of alanine aminotransferase (ALT) and aspartate aminotransferase (AST), level of reactive oxygen species (ROS) (using the DCFH-DA fluorescent probe method), contents of malondialdehyde (MDA) and glutathione(GSH) in liver tissues (using a colorimetric method) and expression of OSGEP (using Western blot). Experiment Ⅱ The well-growing AML12 cells were divided into 4 groups ( n=30 each) using a random number table method: control group (C group), oxygen-glucose deprivation/restoration (OGD/R) group, OGD/R+ OSGEP knockdown group (OGD/R+ KD group), and OGD/R+ OSGEP knockdown negative control group (OGD/R+ NC group). Group C was cultured under normal conditions. Group OGD/R was subjected to O 2-glucose deprivation for 6 h followed by restoration of O 2-glucose supply for 24 h in OGD/R group. In OGD/R+ KD group, stable transfection of AML12 cells with OSGEP knockdown was performed prior to the experiment, and the other procedures were the same as those previously described. The cell survival rate was measured by the CCK-8 assay, the release of lactate dehydrogenase (LDH) was measured, the DCFH-DA method was used to detect the levels of ROS, and the contents of MDA and GSH were determined using a colorimetric method. Results:Experiment Ⅰ Compared with Sham group, the expression of OSGEP was significantly down-regulated, the serum concentrations of AST and ALT, Suzuki score, levels of ROS and content of MDA were increased, and the GSH content was decreased in HIRI group ( P<0.05), and no significant change was found in each parameter in Sham+ KD group ( P>0.05). Compared with HIRI group, the serum concentrations of AST and ALT, Suzuki score, levels of ROS and content of MDA were significantly increased, and the GSH content was decreased in HIRI+ KD group ( P<0.05). Experiment Ⅱ Compared with group C, the expression of OSGEP was significantly down-regulated, the cell survival rate and GSH content were decreased, and the release of LDH, levels of ROS and content of MDA were increased in group OGD/R ( P<0.05). Compared with OGD/R group, the cell survival rate and GSH content were significantly decreased, and the release of LDH, levels of ROS and content of MDA were increased in OGD/R+ KD group ( P<0.05), and no significant change was found in each parameter in OGD/R+ NC group ( P>0.05). Conclusions:OSGEP plays an endogenous protective role in HIRI by inhibiting oxidative stress in mice.
6.Establishment of a Method for Galvanic Vestibular Stimulation-vestibular Evoked Myogenic Potentials in Healthy Children
Zichen CHEN ; Juan HU ; Feiyun CHEN ; Hui YANG ; Yanfei CHEN ; Tingting XUE ; Fangyuan YANG ; Yuzhong ZHANG ; Qiong WU ; Yulian JIN ; Xiaoyong REN ; Qing ZHANG
Journal of Audiology and Speech Pathology 2024;32(2):100-106
Objective To establish the methods of galvanic vestibular stimulation-vestibular evoked myogenic potentials(GVS-VEMPs)in healthy children and to obtain the normal value of GVS-cVEMP and GVS-oVEMP in these children in China.Methods Twenty(3~14 years)healthy children and 24 healthy adults(18~30 years)were enrolled for conventional examinations of GVS-cVEMP and GVS-oVEMP.Using the galvanic stimulation in-tensity under 3 mA/1 ms for children and 5 mA/1 ms for adults.The characteristics of elicitation and parameter re-sults of GVS-cVEMP and GVS-oVEMP in children and adults,as well as the pain scores and the elicitation of differ-ent stimulus intensities in the two age groups were recorded.Results The elicitation of GVS-cVEMP and GVS-oVEMP were both 100.0%in children and adult groups.The p1 latency,n1 latency and p1-n1 interval latency of GVS-cVEMP were 10.46±1.84 ms,16.98±2.12 ms and 6.52±1.42 ms respectively in children group,the n1 la-tency and p1-n1 interval latency were significantly shorter than the adult group(P<0.05).The n1 latency,p1 la-tency and p1-n1 interval latency of GVS-oVEMP were 8.87±1.40 ms,12.25±1.80 ms and 3.39±1.07 ms re-spectively in children group with no significant difference between the two groups.The thresholds of GVS-cVEMP and GVS-oVEMP in children group were significantly lower than adult group(P<0.01),but no differences were found in adult group regarding on the amplitude and interaural amplitude asymmetry ratio.In addition,with the in-crease of the intensity of galvanic stimulation,the correlation between pain scores and the elicitation rates of GVS-cVEMP and GVS-oVEMP also increased.Conclusion Using appropriate stimulus intensity and recording methods,GVS-cVEMP and GVS-oVEMP could be successfully assessed and detected in healthy children over 3 years old and adolescents.The latency of GVS-cVEMP in children is slightly shorter than that in adults,therefore we recommend selecting the matched age group for assessment in the children group.
7.Willingness to preventive treatments and related factors among college freshmen with latent tuberculosis infection in Changzhou
Chinese Journal of School Health 2024;45(12):1802-1806
Objective:
To investigate the willingness to accept preventive treatments and its related factors among college freshmen with latent tuberculosis infection (LTBI), so as to provide the evidence for preventive treatment intervention measures for students with LTBI.
Methods:
Cluster sampling method was used to select 368 LTBI freshmen from 8 colleges and universities in Changzhou in September 2023, who conducted a questionnaire survey on the willingness to receive preventive treatment. General demographic data were collected and relevant data were collected using tuberculosis knowledge scale, General Anxiety Disorder-7 (GAD-7), Patient Health Questionnaire-9 (PHQ-9), Adaptation, Partnership, Growth, Affection and Resolve (APGAR), and a self developed Stigma Scale. A binary Logistic regression model was constructed with the willingness to accept preventive treatment as the dependent variable to analyze the willingness to accept preventive treatment and the influencing factors.
Results:
A total of 253 LTBI college freshmen were willing to take preventive treatment, the acceptance rate was 68.75%. The rate of willingness to accept preventive treatment for LTBI was higher among students whose fathers had an education level of high school, compared to those whose fathers had an education level of junior high school or below ( OR =2.16, P <0.05). LTBI students whose per capita family income was >5 000-10 000 yuan and >10 000 yuan were more willing to accept LTBI preventive treatment than those whose per capita family income was <3 000 yuan ( OR =2.72, 4.46, P <0.05). LTBI students who engaged in physical exercise for more than 2 hours per week were more willing to accept than those who exercised less than 0.5 hours per week ( OR =1.91, P <0.05). LTBI students with high levels of tuberculosis knowledge and stigma were more likely to receive preventive treatment ( OR =1.18, 1.11, P < 0.05). LTBI students with high PHQ-9 ( OR =0.85) and GAD-7 ( OR =0.92) scores were more likely to refuse preventive treatment ( P <0.05).
Conclusion
The present study revealed a moderate level of willingness of LTBI students to preventive treatment in Changzhou City, and the acceptance is affected by family factors, healthy lifestyles, tuberculosis knowledge and psychological status.
8.The current status of palliative sedation for end-of-life care of elderly patients and research progress in China and the world
Tonghui FENG ; Mengting GU ; Qiaozhen XIANG ; Xinyu CHEN ; Yanfei XIA
Chinese Journal of Geriatrics 2024;43(1):98-102
With the rapid acceleration of aging in China, there is a huge need for elderly patients to have improved quality of life in the terminal stage.Palliative sedation is an integral part of palliative care and can alleviate painful refractory symptoms, and its use in patients in various terminal illnesses is being explored across the world.Attention is focused on its indications and implementation.In China, palliative sedation in clinical practice is in an early exploratory stage and relevant criteria and guidelines have yet to be established.A review of the current practice and research progress concerning palliative sedation for patients' end-of-life care in China and the rest of the world will offer insight and strategic considerations in the initial pursuit and accelerated acceptance in the future in China.
9.Construction and external validation of a risk prediction model for unplanned interruption during continuous renal replacement therapy
Hongyan XU ; Qi REN ; Lihong ZHU ; Juan LIN ; Shangzhong CHEN ; Caibao HU ; Yanfei SHEN ; Guolong CAI
Chinese Critical Care Medicine 2024;36(5):520-526
Objective:To identify the independent factors of unplanned interruption during continuous renal replacement therapy (CRRT) and construct a risk prediction model, and to verify the clinical application effectiveness of the model.Methods:A retrospective study was conducted on critically ill adult patients who received CRRT treatment in the intensive care unit (ICU) of Zhejiang Hospital from January 2021 to August 2022 for model construction. According to whether unplanned weaning occurred, the patients were divided into two groups. The potential influencing factors of unplanned CRRT weaning in the two groups were compared. The independent influencing factors of unplanned CRRT weaning were screened by binary Logistic regression and a risk prediction model was constructed. The goodness of fit of the model was verified by a Hosmer-Lemeshow test and its predictive validity was evaluated by receiver operator characteristic curve (ROC curve). Then embed the risk prediction model into the hospital's ICU multifunctional electronic medical record system for severe illness, critically ill patients with CRRT admitted to the ICU of Zhejiang Hospital from November 2022 to October 2023 were prospectively analyzed to verify the model's clinical application effect.Results:① Model construction and internal validation: a total of 331 critically ill patients with CRRT were included to be retrospectively analyzed. Among them, there were 238 patients in planned interruption group and 93 patients in unplanned interruption group. Compared with the planned interruption group, the unplanned interruption group was shown as a lower proportion of males (80.6% vs. 91.6%) and a higher proportion of chronic diseases (60.2% vs. 41.6%), poor blood purification catheter function (31.2% vs. 6.3%), as a higher platelet count (PLT) before CRRT initiation [×10 9/L: 137 (101, 187) vs. 109 (74, 160)], lower level of blood flow rate [mL/min: 120 (120, 150) vs. 150 (140, 180)], higher proportion of using pre-dilution (37.6% vs. 23.5%), higher filtration fraction [23.0% (17.5%, 32.9%) vs. 19.1% (15.7%, 22.6%)], and frequency of blood pump stops [times: 19 (14, 21) vs. 9 (6, 13)], the differences of the above 8 factors between the two groups were statistically significant (all P < 0.05). Binary Logistic regression analysis showed that chronic diseases [odds ratio ( OR) = 3.063, 95% confidence interval (95% CI) was 1.200-7.819], blood purification catheter function ( OR = 4.429, 95% CI was 1.270-15.451), blood flow rate ( OR = 0.928, 95% CI was 0.900-0.957), and frequency of blood pump stops ( OR = 1.339, 95% CI was 1.231-1.457) were the independent factors for the unplanned interruption of CRRT (all P < 0.05). These 4 factors were used to construct a risk prediction model, and ROC curve analysis showed that the area under the curve (AUC) predicted by the model was 0.952 (95% CI was 0.930-0.973, P = 0.003 0), with a sensitivity of 88.2%, a specificity of 89.9%, and a maximum value of 1.781 for the Youden index. ② External validation: prospective inclusion of 110 patients, including 63 planned interruption group and 47 unplanned interruption group. ROC curve analysis showed that the AUC of the risk prediction model was 0.919 (95% CI was 0.870-0.969, P = 0.004 3), with a sensitivity of 91.5%, a specificity of 79.4%, and a maximum value of the Youden index of 1.709. Conclusion:The risk prediction model for unplanned interruption during CRRT has a high predictive efficiency, allowing for rapid and real-time identification of the high risk patients, thus providing references for preventative nursing.
10.Angle of progression measured with transperineal ultrasound in the second stage of labor for predicting delivery mode:Meta-analysis
Sen LIU ; Zhenyu CHEN ; Yanfei XING
Chinese Journal of Medical Imaging Technology 2024;40(5):752-757
Objective To observe the value of angle of progression(AOP)measured with transperineal ultrasound in the second stage of labor for predicting delivery mode using meta-analysis.Methods Articles published in Chinese or English focused on the prediction of delivery mode using AOP measured with ultrasound during the second stage of labor were searched in databases from January 1,2018 to August 31,2023.The literature were screened by 2 researchers according to the inclusion and exclusion criteria,and information were extracted to assess the value of AOP measured in the second stage of labor with transperineal ultrasound for predicting delivery mode.The quality of the included literature were evaluated using Quality Evaluation Tool for Diagnostic Tests-2(QUADAS-2).Results Eleven studies including 2 315 parturient women were enrolled.Meta-analysis showed that the combined sensitivity(Sen),combined specificity(Spe),combined positive likelihood ratio(+LR),combined negative likelihood ratio(-LR),and combined diagnostic ratio(DOR)of AOP measured with transperineal ultrasound in the second stage of labor for predicting delivery mode was 0.80[95%CI(0.76,0.83)],0.83[95%CI(0.78,0.88)],4.8[95%CI(3.5,6.6)],0.24[95%CI(0.20,0.29)]and 20[95%CI(13,30)],respectively,while area under the summary receiver operating characteristic curve of 0.88.Conclusion AOP measured with transperineal ultrasound in the second stage of labor had fine efficiency for predicting delivery mode,which could be regarded as a useful predictor for probability of natural delivery.


Result Analysis
Print
Save
E-mail