1.Compact Fundus Imaging System Using Shack-Hartmann Wavefront Sensing for High-speed Auto-focus
Zhe-Kai LIN ; Long CHEN ; Geng-Yong ZHENG ; Jin-Tian HUANG ; Jia-Xin DONG ; Shang-Pan YANG ; Wen-Zheng DING ; Ding-An HAN ; Xue-Hua WANG ; Ya-Guang ZENG
Progress in Biochemistry and Biophysics 2026;53(4):1076-1086
ObjectiveThe widespread adoption of portable fundus cameras for primary care and community screening is hindered by limitations in current autofocus(AF) technologies. Image-based methods relying on sharpness evaluation require iterative searches, resulting in slow convergence, while projection-based techniques are susceptible to optical artifacts and calibration errors. To address these challenges, this study introduces a novel AF system based on direct wavefront sensing, designed to deliver simultaneous high speed, high precision, and operational robustness within the compact form factor essential for portable ophthalmic devices. MethodsOur approach fundamentally reimagines the AF process by directly measuring the ocular wavefront aberration. We developed a custom portable fundus camera integrating a miniaturized Shack-Hartmann wavefront sensor (SHWS) into the optical path. An 850 nm laser diode projects a point source onto the retina via oblique illumination to minimize corneal reflections. Light scattered from this spot carries the eye’s refractive error through the imaging optics and is directed to the SHWS, positioned at a plane optically conjugate to the primary color CMOS imaging sensor. A microlens array within the SHWS samples the incident wavefront, generating a pattern of focal spots on a CCD. Real-time centroid analysis of these spots provides a map of local wavefront slopes. These measurements are processed through a singular value decomposition (SVD) algorithm to fit a Zernike polynomial basis set, enabling real-time reconstruction of the wavefront phase. The defocus component (S) is extracted from the second-order Zernike coefficients, providing a direct, quantitative measure of the refractive error in diopters. This value serves as a precise error signal in a closed-loop control system, which commands a voice-coil actuated focusing lens to its null position in a single, deterministic step, eliminating the need for iterative search algorithms. ResultsComprehensive evaluation demonstrated the system’s high performance. Testing on a calibrated model eye (OEMI-7) established a highly linear relationship between the computed defocus S and the focusing lens position across a ±20 Diopter (D) compensation range, achievable within a 5 mm mechanical travel. The system achieved a focusing precision of 0.08 D, corresponding to an 18-fold improvement over a conventional projection spot-size method tested under identical conditions. The total focus acquisition time, encompassing wavefront measurement, computation, and lens actuation, averaged under 0.5 s. Clinical validation with 25 human volunteers (50 eyes, refractive range -15 D to +10 D) confirmed practical efficacy. The wavefront-sensing AF succeeded in 92% of attempts with a mean time of 0.5 s, substantially outperforming a projection-based benchmark which achieved only a 32% success rate with an average time of 4.25 s. The system provided instantaneous directional guidance and maintained stability during minor ocular movements. Objective assessment of image quality, via amplitude contrast of retinal vasculature, showed consistent and significant enhancement following AF correction across the entire tested diopter range. ConclusionThis work successfully implements and validates a direct wavefront-sensing autofocus paradigm for portable fundus cameras. By directly quantifying and compensating for the optical defocus aberration, this method bypasses the fundamental limitations of image-processing and projection-based techniques, enabling rapid, precise, and deterministic diopter compensation. The developed system delivers an exceptional combination of a wide operational range (±20 D), high accuracy (0.08 D), fast convergence (0.5 s), and a compact physical footprint. This technology provides a practical and high-performance focusing solution capable of enhancing the reliability, throughput, and diagnostic utility of portable retinal imaging in large-scale screening applications. Future efforts will be directed towards system cost optimization and performance adaptation for diverse ocular conditions.
2.An Attention-weighted Tri-modal Ultrasound Network (TUS-Net) for Screening of Atypical Hepatocellular Carcinoma From LR-M Liver Nodules
He-Chong ZHANG ; Liang-Hui HUANG ; Xue-Hua WANG ; Shang-Lin JIANG ; Ying-Ying CHEN ; Ya-Guang ZENG ; Wei ZHENG
Progress in Biochemistry and Biophysics 2026;53(5):1485-1498
ObjectiveDiscriminating atypical hepatocellular carcinoma (HCC) from other malignancies in liver nodules classified as Liver Imaging Reporting and Data System category M (LR-M) remains a significant diagnostic challenge on conventional ultrasound examination. The LR-M category, originally intended to capture non-HCC malignancies, paradoxically contains up to 63% of atypical HCCs that deviate from classic enhancement patterns, leading to potential misdiagnosis and suboptimal treatment planning. While deep learning has shown promise in HCC diagnosis, most existing models rely exclusively on single-modality ultrasound, overlooking the diagnostic benefits of integrating complementary information from multiple imaging sources. To address this gap, we propose a novel attention-weighted tri-modal ultrasound network (TUS-Net) that integrates contrast-enhanced ultrasound (CEUS), B-mode ultrasound (BUS), and time-intensity curves (TICs) to improve diagnostic accuracy for these clinically challenging lesions. MethodsOur framework incorporates a three-dimensional convolutional neural network (C3D) backbone to extract spatiotemporal features from CEUS videos, capturing dynamic vascular patterns critical for lesion characterization. To effectively fuse complementary modalities, we introduce a dual-channel feature fusion module (DCFFM) that adaptively combines features from CEUS and BUS through channel-wise attention mechanisms, allowing the model to dynamically weigh the contribution of each modality based on diagnostic relevance. Additionally, we propose a temporal intensity feature fusion module (TIFFM) that leverages quantitative hemodynamic information from TICs to guide the model’s attention toward diagnostically critical temporal phases, such as arterial wash-in and portal venous washout. The model is further enhanced by automated lesion localization using YOLOX and class activation mapping for interpretability, ensuring that predictions align with clinically meaningful imaging features. ResultsEvaluated on a tri-modal ultrasound dataset comprising 161 patients with pathologically confirmed LR-M nodules (131 atypical HCC and 30 non-HCC malignancies), our model achieved an accuracy of 86.83%, a sensitivity of 92.50%, a specificity of 75.50%, and an AUC of 89.32% in screening atypical HCC. Compared to single-modality baselines, TUS-Net demonstrated superior specificity, a clinically critical metric given the higher risk associated with misclassifying non-HCC malignancies. Ablation studies confirmed the contribution of each module, with the full model outperforming both standard C3D and 3D ResNet backbones integrated with attention mechanisms. A reader study involving junior and senior radiologists further validated the clinical utility of AI assistance, showing consistent improvements in specificity and inter-reader consistency, particularly for less experienced clinicians. ConclusionThese results surpass existing benchmark models and demonstrate the potential of our approach to enhance diagnostic precision in clinically specific cases. By intelligently fusing multi-modal ultrasound data with attention-guided mechanisms, TUS-Net offers a reliable and interpretable tool that holds promise for improving the non-invasive diagnosis of atypical HCC in challenging LR-M liver nodules.
3.An Attention-weighted Tri-modal Ultrasound Network (TUS-Net) for Screening of Atypical Hepatocellular Carcinoma From LR-M Liver Nodules
He-Chong ZHANG ; Liang-Hui HUANG ; Xue-Hua WANG ; Shang-Lin JIANG ; Ying-Ying CHEN ; Ya-Guang ZENG ; Wei ZHENG
Progress in Biochemistry and Biophysics 2026;53(5):1485-1498
ObjectiveDiscriminating atypical hepatocellular carcinoma (HCC) from other malignancies in liver nodules classified as Liver Imaging Reporting and Data System category M (LR-M) remains a significant diagnostic challenge on conventional ultrasound examination. The LR-M category, originally intended to capture non-HCC malignancies, paradoxically contains up to 63% of atypical HCCs that deviate from classic enhancement patterns, leading to potential misdiagnosis and suboptimal treatment planning. While deep learning has shown promise in HCC diagnosis, most existing models rely exclusively on single-modality ultrasound, overlooking the diagnostic benefits of integrating complementary information from multiple imaging sources. To address this gap, we propose a novel attention-weighted tri-modal ultrasound network (TUS-Net) that integrates contrast-enhanced ultrasound (CEUS), B-mode ultrasound (BUS), and time-intensity curves (TICs) to improve diagnostic accuracy for these clinically challenging lesions. MethodsOur framework incorporates a three-dimensional convolutional neural network (C3D) backbone to extract spatiotemporal features from CEUS videos, capturing dynamic vascular patterns critical for lesion characterization. To effectively fuse complementary modalities, we introduce a dual-channel feature fusion module (DCFFM) that adaptively combines features from CEUS and BUS through channel-wise attention mechanisms, allowing the model to dynamically weigh the contribution of each modality based on diagnostic relevance. Additionally, we propose a temporal intensity feature fusion module (TIFFM) that leverages quantitative hemodynamic information from TICs to guide the model’s attention toward diagnostically critical temporal phases, such as arterial wash-in and portal venous washout. The model is further enhanced by automated lesion localization using YOLOX and class activation mapping for interpretability, ensuring that predictions align with clinically meaningful imaging features. ResultsEvaluated on a tri-modal ultrasound dataset comprising 161 patients with pathologically confirmed LR-M nodules (131 atypical HCC and 30 non-HCC malignancies), our model achieved an accuracy of 86.83%, a sensitivity of 92.50%, a specificity of 75.50%, and an AUC of 89.32% in screening atypical HCC. Compared to single-modality baselines, TUS-Net demonstrated superior specificity, a clinically critical metric given the higher risk associated with misclassifying non-HCC malignancies. Ablation studies confirmed the contribution of each module, with the full model outperforming both standard C3D and 3D ResNet backbones integrated with attention mechanisms. A reader study involving junior and senior radiologists further validated the clinical utility of AI assistance, showing consistent improvements in specificity and inter-reader consistency, particularly for less experienced clinicians. ConclusionThese results surpass existing benchmark models and demonstrate the potential of our approach to enhance diagnostic precision in clinically specific cases. By intelligently fusing multi-modal ultrasound data with attention-guided mechanisms, TUS-Net offers a reliable and interpretable tool that holds promise for improving the non-invasive diagnosis of atypical HCC in challenging LR-M liver nodules.
4.Clinical characteristics and survival analysis of pediatric Hodgkin lymphoma: a multicenter study.
Ying LIN ; Li-Li PAN ; Shao-Hua LE ; Jian LI ; Bi-Yun GUO ; Yu ZHU ; Kai-Zhi WENG ; Jin-Hong LUO ; Gao-Yuan SUN ; Yong-Zhi ZHENG
Chinese Journal of Contemporary Pediatrics 2025;27(6):668-674
OBJECTIVES:
To investigate the clinicopathological characteristics and prognostic factors of pediatric Hodgkin lymphoma (HL).
METHODS:
A retrospective analysis was conducted on the clinical data of children with newly diagnosed HL from January 2011 to December 2023 at four hospitals: Fujian Medical University Union Hospital, Fujian Medical University Zhangzhou Hospital, First Affiliated Hospital of Xiamen University, and Fujian Children's Hospital. Patients were categorized into low-risk (R1), intermediate-risk (R2), and high-risk (R3) groups based on HL staging and pre-treatment risk factors. The patients received ABVD regimen or Chinese Pediatric HL-2013 regimen chemotherapy. Early treatment response and long-term efficacy were assessed, and prognostic factors were analyzed using the Cox proportional hazards regression model.
RESULTS:
The overall complete response (CR) rates after 2 and 4 cycles of chemotherapy were 42% and 68%, respectively. Compared with the ABVD regimen group, patients treated with the HL-2013 regimen in the R1 group showed significantly higher CR rates after both 2 and 4 cycles (P<0.05). However, no statistically significant differences in CR rates were observed between the two regimens in the R2 and R3 groups (P>0.05). The 5-year event-free survival (EFS) rate, overall survival rate, and freedom from treatment failure rate were 83%±4%, 97%±2%, and 88%±4%, respectively. Cox analysis indicated that the presence of a large tumor mass at diagnosis and failure to achieve CR after 4 cycles of chemotherapy were independent risk factors for lower EFS rates (P<0.05).
CONCLUSIONS
Pediatric HL generally has a favorable prognosis. The presence of a large tumor mass at diagnosis and failure to achieve CR after 4 cycles of chemotherapy indicate poor prognosis.
Humans
;
Hodgkin Disease/pathology*
;
Male
;
Child
;
Female
;
Adolescent
;
Retrospective Studies
;
Child, Preschool
;
Antineoplastic Combined Chemotherapy Protocols/therapeutic use*
;
Prognosis
;
Proportional Hazards Models
;
Survival Analysis
;
Infant
5.RNA Sequencing Reveals Molecular Alternations of Splenocytes Associated with Anti-FⅧ Immune Response in Hemophilia A Murine Model.
Chen-Chen WANG ; Ya-Li WANG ; Yuan-Hua CAI ; Qiao-Yun ZHENG ; Zhen-Xing LIN ; Ying-Yu CHEN
Journal of Experimental Hematology 2025;33(5):1476-1485
OBJECTIVE:
To investigate the molecular alterations of splenocytes associated with anti-factor Ⅷ (FⅧ) immune response and the underlying mechanisms based on hemophilia A (HA) murine model via RNA sequencing (RNA-seq) technology.
METHODS:
Severe HA mice were immunized with recombinant human factor Ⅷ (rhF8) weekly for 4 weeks to establish an FⅧ inhibitor model. High quality raw data were obtained by using bulk RNA-seq and CASAVA base identification technology, and the differentially expressed genes (DEGs) were identified. The DEGs were statistically classified by gene ontology (GO) annotation to obtain information on the major signaling pathways and biological processes involved in anti-FⅧ immune response in HA mouse splenocytes. The cell clusters, genes, and signaling pathway datasets were comprehensively analyzed by GO, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and single cell RNA-seq (ScRNA-seq) analysis, respectively. Flow cytometry analysis was used to verify the changes in T follicular helper cells (Tfh) and regulatory T cells (Treg).
RESULTS:
A total of 3731 DEGs was identified, including 2275 genes with up-regulated expression and 1456 genes with down-regulated expression. The DEGs were enriched in helper T cell differentiation, cytokine receptor, T cell receptor signaling pathway, ferroptosis, etc. Uniform Manifold Approximation and Project (UMAP) downscaling and visualization analysis yielded a total number of 11 T/NK cell subsets, visualizing the overall expression distribution of C-X-C chemokine-specific receptor gene cxcr5 among these T/NK cell subsets. Higher expression of cxcr5 was found in activated Tfh from FⅧ inhibitor mice, in comparison to the control group. The visualization using Upset plot R language showed a close interaction between Tfh and Treg. Moreover, the increased frequencies of Tfh and the decreased frequencies of Treg in inhibitor mouse splenocytes were further verified by flow cytometry analysis.
CONCLUSION
Multiple immune cell subsets, signaling pathways, and characteristic genes may be involved in the process of anti-FⅧ immune response in HA mouse splenocytes. The molecules involved in the regulation of Tfh/Treg may play key roles, which provide potential biological targets and therapeutic strategies for HA patients with inhibitors in the future.
Animals
;
Hemophilia A/genetics*
;
Mice
;
Sequence Analysis, RNA
;
Disease Models, Animal
;
Spleen/cytology*
;
T-Lymphocytes, Regulatory/immunology*
;
Humans
;
Signal Transduction
;
Factor VIII/immunology*
;
T-Lymphocytes, Helper-Inducer/immunology*
6.Prognostic value of quantitative flow ratio measured immediately after percutaneous coronary intervention for chronic total occlusion.
Zheng QIAO ; Zhang-Yu LIN ; Qian-Qian LIU ; Rui ZHANG ; Chang-Dong GUAN ; Sheng YUAN ; Tong-Qiang ZOU ; Xiao-Hui BIAN ; Li-Hua XIE ; Cheng-Gang ZHU ; Hao-Yu WANG ; Guo-Feng GAO ; Ke-Fei DOU
Journal of Geriatric Cardiology 2025;22(4):433-442
BACKGROUND:
The clinical impact of post-percutaneous coronary intervention (PCI) quantitative flow ratio (QFR) in patients treated with PCI for chronic total occlusion (CTO) was still undetermined.
METHODS:
All CTO vessels treated with successful anatomical PCI in patients from PANDA III trial were retrospectively measured for post-PCI QFR. The primary outcome was 2-year vessel-oriented composite endpoints (VOCEs, composite of target vessel-related cardiac death, target vessel-related myocardial infarction, and ischemia-driven target vessel revascularization). Receiver operator characteristic curve analysis was conducted to identify optimal cutoff value of post-PCI QFR for predicting the 2-year VOCEs, and all vessels were stratified by this optimal cutoff value. Cox proportional hazards models were employed to calculate the hazard ratio (HR) with 95% CI.
RESULTS:
Among 428 CTO vessels treated with PCI, 353 vessels (82.5%) were analyzable for post-PCI QFR. 31 VOCEs (8.7%) occurred at 2 years. Mean value of post-PCI QFR was 0.92 ± 0.13. Receiver operator characteristic curve analysis shown the optimal cutoff value of post-PCI QFR for predicting 2-year VOCEs was 0.91. The incidence of 2-year VOCEs in the vessel with post-PCI QFR < 0.91 (n = 91) was significantly higher compared with the vessels with post-PCI QFR ≥ 0.91 (n = 262) (22.0% vs. 4.2%, HR = 4.98, 95% CI: 2.32-10.70).
CONCLUSIONS
Higher post-PCI QFR values were associated with improved prognosis in the PCI practice for coronary CTO. Achieving functionally optimal PCI results (post-PCI QFR value ≥ 0.91) tends to get better prognosis for patients with CTO lesions.
7.Qingda Granule Attenuates Hypertension-Induced Cardiac Damage via Regulating Renin-Angiotensin System Pathway.
Lin-Zi LONG ; Ling TAN ; Feng-Qin XU ; Wen-Wen YANG ; Hong-Zheng LI ; Jian-Gang LIU ; Ke WANG ; Zhi-Ru ZHAO ; Yue-Qi WANG ; Chao-Ju WANG ; Yi-Chao WEN ; Ming-Yan HUANG ; Hua QU ; Chang-Geng FU ; Ke-Ji CHEN
Chinese journal of integrative medicine 2025;31(5):402-411
OBJECTIVE:
To assess the efficacy of Qingda Granule (QDG) in ameliorating hypertension-induced cardiac damage and investigate the underlying mechanisms involved.
METHODS:
Twenty spontaneously hypertensive rats (SHRs) were used to develope a hypertension-induced cardiac damage model. Another 10 Wistar Kyoto (WKY) rats were used as normotension group. Rats were administrated intragastrically QDG [0.9 g/(kg•d)] or an equivalent volume of pure water for 8 weeks. Blood pressure, histopathological changes, cardiac function, levels of oxidative stress and inflammatory response markers were measured. Furthermore, to gain insights into the potential mechanisms underlying the protective effects of QDG against hypertension-induced cardiac injury, a network pharmacology study was conducted. Predicted results were validated by Western blot, radioimmunoassay immunohistochemistry and quantitative polymerase chain reaction, respectively.
RESULTS:
The administration of QDG resulted in a significant decrease in blood pressure levels in SHRs (P<0.01). Histological examinations, including hematoxylin-eosin staining and Masson trichrome staining revealed that QDG effectively attenuated hypertension-induced cardiac damage. Furthermore, echocardiography demonstrated that QDG improved hypertension-associated cardiac dysfunction. Enzyme-linked immunosorbent assay and colorimetric method indicated that QDG significantly reduced oxidative stress and inflammatory response levels in both myocardial tissue and serum (P<0.01).
CONCLUSIONS
Both network pharmacology and experimental investigations confirmed that QDG exerted its beneficial effects in decreasing hypertension-induced cardiac damage by regulating the angiotensin converting enzyme (ACE)/angiotensin II (Ang II)/Ang II receptor type 1 axis and ACE/Ang II/Ang II receptor type 2 axis.
Animals
;
Drugs, Chinese Herbal/therapeutic use*
;
Hypertension/pathology*
;
Renin-Angiotensin System/drug effects*
;
Rats, Inbred SHR
;
Oxidative Stress/drug effects*
;
Male
;
Rats, Inbred WKY
;
Blood Pressure/drug effects*
;
Myocardium/pathology*
;
Rats
;
Inflammation/pathology*
8.CHAF1B promotes the progression of lung squamous-cell carcinoma by inhibiting SETD7 expression.
Zhuo ZHENG ; Yongfang LIN ; Hua GUO ; Zheng LIU ; Xiaoliang JIE ; Guizhen WANG ; Guangbiao ZHOU
Frontiers of Medicine 2025;19(2):318-328
The p60 subunit of the chromatin assembly factor-1 complex, that is, chromatin assembly factor-1 subunit B (CHAF1B), is a histone H3/H4 chaperone crucial for the transcriptional regulation of cell differentiation and self-renewal. CHAF1B is overexpressed in several cancers and may represent a potential target for cancer therapy. However, its expression and clinical significance in lung squamous-cell carcinoma (LUSC) remain unclear. In this study, we performed weighted gene correlation network analysis to analyze the Gene Expression Omnibus GSE68793 LUSC dataset and identified CHAF1B as one of the most important driver gene candidates. Immunohistochemical analysis of 126 LUSC tumor samples and 80 adjacent normal lung tissues showed the marked upregulation of CHAF1B in tumor tissues and the negative association of its expression level with patient survival outcomes. Silencing of CHAF1B suppressed LUSC proliferation in vitro and LUSC tumor growth in vivo. Furthermore, bulk RNA sequencing of CHAF1B knockdown cells indicated SET domain containing 7 (SETD7) as a significant CHAF1B target gene. In addition, CHAF1B competitively binds to the SETD7 promoter region and represses its transcription. Altogether, these results imply that CHAF1B plays a vital role in LUSC tumorigenesis and may represent a potential molecular target for this deadly disease.
Humans
;
Lung Neoplasms/metabolism*
;
Histone-Lysine N-Methyltransferase/metabolism*
;
Carcinoma, Squamous Cell/metabolism*
;
Gene Expression Regulation, Neoplastic
;
Disease Progression
;
Cell Proliferation/genetics*
;
Cell Line, Tumor
;
Chromatin Assembly Factor-1/metabolism*
;
Animals
;
Mice
;
Male
;
Female
9.W 18O 49 Crystal and ICG Labeled Macrophage: An Efficient Targeting Vector for Fluorescence Imaging-guided Photothermal Therapy.
Yang BAI ; Guo Qing FENG ; Muskan Saif KHAN ; Qing Bin YANG ; Ting Ting HUA ; Hao Lin GUO ; Yuan LIU ; Bo Wen LI ; Yi Wen WU ; Bin ZHENG ; Nian Song QIAN ; Qing YUAN
Biomedical and Environmental Sciences 2025;38(1):100-105
10.Generalized Functional Linear Models: Efficient Modeling for High-dimensional Correlated Mixture Exposures.
Bing Song ZHANG ; Hai Bin YU ; Xin PENG ; Hai Yi YAN ; Si Ran LI ; Shutong LUO ; Hui Zi WEIREN ; Zhu Jiang ZHOU ; Ya Lin KUANG ; Yi Huan ZHENG ; Chu Lan OU ; Lin Hua LIU ; Yuehua HU ; Jin Dong NI
Biomedical and Environmental Sciences 2025;38(8):961-976
OBJECTIVE:
Humans are exposed to complex mixtures of environmental chemicals and other factors that can affect their health. Analysis of these mixture exposures presents several key challenges for environmental epidemiology and risk assessment, including high dimensionality, correlated exposure, and subtle individual effects.
METHODS:
We proposed a novel statistical approach, the generalized functional linear model (GFLM), to analyze the health effects of exposure mixtures. GFLM treats the effect of mixture exposures as a smooth function by reordering exposures based on specific mechanisms and capturing internal correlations to provide a meaningful estimation and interpretation. The robustness and efficiency was evaluated under various scenarios through extensive simulation studies.
RESULTS:
We applied the GFLM to two datasets from the National Health and Nutrition Examination Survey (NHANES). In the first application, we examined the effects of 37 nutrients on BMI (2011-2016 cycles). The GFLM identified a significant mixture effect, with fiber and fat emerging as the nutrients with the greatest negative and positive effects on BMI, respectively. For the second application, we investigated the association between four pre- and perfluoroalkyl substances (PFAS) and gout risk (2007-2018 cycles). Unlike traditional methods, the GFLM indicated no significant association, demonstrating its robustness to multicollinearity.
CONCLUSION
GFLM framework is a powerful tool for mixture exposure analysis, offering improved handling of correlated exposures and interpretable results. It demonstrates robust performance across various scenarios and real-world applications, advancing our understanding of complex environmental exposures and their health impacts on environmental epidemiology and toxicology.
Humans
;
Environmental Exposure/analysis*
;
Linear Models
;
Nutrition Surveys
;
Environmental Pollutants
;
Body Mass Index

Result Analysis
Print
Save
E-mail