The continuous negative external pressure (cNEP) applied on the submental surface is a method of non-surgical treatment for obstructive sleep apnea (OSA), which can effectively widen the airway in some OSA patients. However, it cannot effectively improve airway collapse in obese patients and its mechanism remains unclear. In this study, we aim to analyze the reasons for the ineffectiveness of cNEP treatment in OSA patients with obesity. Based on magnetic resonance imaging (MRI), three-dimensional models of the head and neck were constructed for a healthy subject, an OSA patient with enlarged tongue, and an OSA patient with the tongue adjusted to normal size. By performing the one step staggered fluid–structure interaction numerical simulations, we analyzed the collapse of the airway in these three models under the influence of cNEP. Restoring the tongue to normal size in the OSA patient significantly improves the airway critical closing pressure under cNEP treatment compared to the patient with enlarged tongue. The enlargement of the tongue in the OSA patient hindered the widening of the velopharyngeal airway under the action of cNEP. The numerical results reveal that cNEP treatment can effectively widen the laryngopharyngeal airway, thus providing a potential therapeutic option for OSA patients with laryngopharyngeal obstruction. Tongue enlargement in OSA patients is a critical factor influencing the efficacy of cNEP treatment. This study reveals the reasons for cNEP treatment failure in obese patients and the potential value of cNEP targeted therapy.
{"title":"The impact of tongue size on submental negative pressure treatment of airway obstruction revealed by fluid-structure interaction simulations","authors":"Yuhang Tian , Huahui Xiong , Hui Tong , Changjin Ji , Xiaoqing Huang , Yaqi Huang","doi":"10.1016/j.bbe.2025.08.004","DOIUrl":"10.1016/j.bbe.2025.08.004","url":null,"abstract":"<div><div>The continuous negative external pressure (cNEP) applied on the submental surface is a method of non-surgical treatment for obstructive sleep apnea (OSA), which can effectively widen the airway in some OSA patients. However, it cannot effectively improve airway collapse in obese patients and its mechanism remains unclear. In this study, we aim to analyze the reasons for the ineffectiveness of cNEP treatment in OSA patients with obesity. Based on magnetic resonance imaging (MRI), three-dimensional models of the head and neck were constructed for a healthy subject, an OSA patient with enlarged tongue, and an OSA patient with the tongue adjusted to normal size. By performing the one step staggered fluid–structure interaction numerical simulations, we analyzed the collapse of the airway in these three models under the influence of cNEP. Restoring the tongue to normal size in the OSA patient significantly improves the airway critical closing pressure under cNEP treatment compared to the patient with enlarged tongue. The enlargement of the tongue in the OSA patient hindered the widening of the velopharyngeal airway under the action of cNEP. The numerical results reveal that cNEP treatment can effectively widen the laryngopharyngeal airway, thus providing a potential therapeutic option for OSA patients with laryngopharyngeal obstruction. Tongue enlargement in OSA patients is a critical factor influencing the efficacy of cNEP treatment. This study reveals the reasons for cNEP treatment failure in obese patients and the potential value of cNEP targeted therapy.</div></div>","PeriodicalId":55381,"journal":{"name":"Biocybernetics and Biomedical Engineering","volume":"45 3","pages":"Pages 590-599"},"PeriodicalIF":6.6,"publicationDate":"2025-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144911947","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-06-27DOI: 10.1016/j.bbe.2025.06.001
Marcin Sińczuk , Jacek Rogala , Piotr Bogorodzki
This study explores the importance of scanner-specific calibration measurements for Magnetic Resonance Spectroscopy Thermometry (MRST) in human brain temperature estimations. Data acquisition was conducted on a 3-T GE scanner. Calibration constants for the water-chemical shift were obtained using a temperature-controlled phantom containing an aqueous solution of N-acetyl aspartate (NAA), Creatine (Cr), and Choline (Cho), and data from three different research groups using the same metabolites. Temperatures were estimated utilizing correlation of water chemical shift with NAA, Cr and Cho. For data acquisition, commercially available single-voxel point-resolved spectroscopy (PRESS) sequences were used for calibrations and in vivo temperature estimations. Each sequence included spectras without (WU) and with (WS) water suppression. In vivo study consisted of two PRESS sequences, one before and one after extensive 30-minute fMRI task acquisition. Significant differences were found between absolute brain temperatures measured using scanner-specific calibrations and those from other researchers, varying from −0.68 °C to + 0.37 °C for NAA, −0.92 °C to 0.37 °C for Cr, and −0.78 °C to 0.7 °C for Cho. Each method reported a similar temperature decrease of −0.26 ∓ 0.03 °C between before and after fMRI measurements. These findings suggest that while absolute temperatures from non-scanner specific calibrations may be inaccurate, comparative estimates are valid.
{"title":"MRS thermometry – Importance of scanner-specific calibrations for accurate brain temperature estimations","authors":"Marcin Sińczuk , Jacek Rogala , Piotr Bogorodzki","doi":"10.1016/j.bbe.2025.06.001","DOIUrl":"10.1016/j.bbe.2025.06.001","url":null,"abstract":"<div><div>This study explores the importance of scanner-specific calibration measurements for Magnetic Resonance Spectroscopy Thermometry (MRST) in human brain temperature estimations. Data acquisition was conducted on a 3-T GE scanner. Calibration constants for the water-chemical shift were obtained using a temperature-controlled phantom containing an aqueous solution of N-acetyl aspartate (NAA), Creatine (Cr), and Choline (Cho), and data from three different research groups using the same metabolites. Temperatures were estimated utilizing correlation of water chemical shift with NAA, Cr and Cho. For data acquisition, commercially available single-voxel point-resolved spectroscopy (PRESS) sequences were used for calibrations and in vivo temperature estimations. Each sequence included spectras without (WU) and with (WS) water suppression. In vivo study consisted of two PRESS sequences, one before and one after extensive 30-minute fMRI task acquisition. Significant differences were found between absolute brain temperatures measured using scanner-specific calibrations and those from other researchers, varying from −0.68 °C to + 0.37 °C for NAA, −0.92 °C to 0.37 °C for Cr, and −0.78 °C to 0.7 °C for Cho. Each method reported a similar temperature decrease of −0.26 ∓ 0.03 °C between before and after fMRI measurements. These findings suggest that while absolute temperatures from non-scanner specific calibrations may be inaccurate, comparative estimates are valid.</div></div>","PeriodicalId":55381,"journal":{"name":"Biocybernetics and Biomedical Engineering","volume":"45 3","pages":"Pages 451-456"},"PeriodicalIF":5.3,"publicationDate":"2025-06-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144489482","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-06-27DOI: 10.1016/j.bbe.2025.06.002
Natalia Kowalczyk, Jacek Rumiński, Magdalena Mazur-Milecka
The COVID-19 pandemic has underscored the importance of wearing facial masks and monitoring respiratory health to prevent the spread of the virus. In this study, we developed a model for segmenting facial masks in thermal images. We applied the model to segment face masks in different conditions, including a person walking toward the observing camera. The segmented regions were further processed using different erosion masks to analyze the influence of the selected sources on the quality of the estimated respiratory signals. The Signal-to-Noise Ratio (SNR) was used as a quality measure. Additionally, the extracted respiratory signals were compared with two reference signals: binary signals generated by participants who signaled the inhalation phase and pressure signals measured with a respiratory belt. Our findings show a high level of concordance between the respiratory signals derived from the segmented mask region and those from the respiratory belt, validating the effectiveness of thermal imaging for capturing respiratory patterns. Notably, the signal-to-noise ratio (SNR) was higher for the segmented mask than the detection methods used in previous works. Specifically, for the mask segmentation task, the mean SNR improved by 4.3 compared to facial mask detection. The segmentation model achieved a mean Average Precision (mAP) of 0.992 for segmentation tasks and 0.857 mAP at the 50–95 % threshold using the Yolov8 “nano” architecture. This study underscores the potential of thermal imaging for non-invasive respiratory monitoring and highlights the explainability and accuracy of selecting the facial mask region for signal extraction.
{"title":"Improving the quality of respiratory signals extracted from the segmented mask area","authors":"Natalia Kowalczyk, Jacek Rumiński, Magdalena Mazur-Milecka","doi":"10.1016/j.bbe.2025.06.002","DOIUrl":"10.1016/j.bbe.2025.06.002","url":null,"abstract":"<div><div>The COVID-19 pandemic has underscored the importance of wearing facial masks and monitoring respiratory health to prevent the spread of the virus. In this study, we developed a model for segmenting facial masks in thermal images. We applied the model to segment face masks in different conditions, including a person walking toward the observing camera. The segmented regions were further processed using different erosion masks to analyze the influence of the selected sources on the quality of the estimated respiratory signals. The Signal-to-Noise Ratio (SNR) was used as a quality measure. Additionally, the extracted respiratory signals were compared with two reference signals: binary signals generated by participants who signaled the inhalation phase and pressure signals measured with a respiratory belt. Our findings show a high level of concordance between the respiratory signals derived from the segmented mask region and those from the respiratory belt, validating the effectiveness of thermal imaging for capturing respiratory patterns. Notably, the signal-to-noise ratio (SNR) was higher for the segmented mask than the detection methods used in previous works. Specifically, for the mask segmentation task, the mean SNR improved by 4.3 compared to facial mask detection. The segmentation model achieved a mean Average Precision (mAP) of 0.992 for segmentation tasks and 0.857 mAP at the 50–95 % threshold using the Yolov8 “nano” architecture. This study underscores the potential of thermal imaging for non-invasive respiratory monitoring and highlights the explainability and accuracy of selecting the facial mask region for signal extraction.</div></div>","PeriodicalId":55381,"journal":{"name":"Biocybernetics and Biomedical Engineering","volume":"45 3","pages":"Pages 457-468"},"PeriodicalIF":5.3,"publicationDate":"2025-06-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144489483","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-06-14DOI: 10.1016/j.bbe.2025.05.011
Víctor Gutiérrez-de Pablo , María Herrero-Tudela , Marina Sandonís-Fernández , Jesús Poza , Aarón Maturana-Candelas , Víctor Rodríguez-González , Miguel Ángel Tola-Arribas , Mónica Cano , Hideyuki Hoshi , Yoshihito Shigihara , Roberto Hornero , Carlos Gómez
Dementia and mild cognitive impairment (MCI) due to Alzheimer’s disease (AD) are neurological pathologies associated with disruptions in brain electromagnetic activity, typically studied using magnetoencephalography (MEG) and electroencephalography (EEG). To quantify diverse brain properties, different families of parameters can be computed from MEG and EEG (i.e., spectral, non-linear, morphological, functional connectivity, or network structure and organisation). However, studying these characteristics separately overlooks the complex nature of brain activity. Integrative frameworks can be useful to unveil the intricate neurophysiological fingerprint, as well as to characterise pathological conditions comprehensively. To that purpose, data fusion methodologies are crucial, despite their interpretational challenges. In this study, Machine Learning (ML) models were trained to discriminate between groups of severity, whereas the SHapley Additive eXplanations (SHAP) algorithm was afterwards utilised to assess the relevance of the input characteristics into the output classification. Three databases were analysed: MEG (55 healthy controls, HC, 42 MCI patients, and 86 AD patients), EEG1 (51 HC, 52 MCI, and 100 AD), and EEG2 (45 HC, 69 MCI, and 82 AD). The best results for the three-class classification problem were obtained by Gradient Boosting for the MEG database: 3-class Cohen’s kappa coefficient of 0.5452 and accuracy of 72.63 %. Afterwards, using SHAP on Gradient Boosting, it has been shown that spectral features were identified as highly relevant across all databases. Furthermore, morphology measures presented high relevance for the MEG database, whereas EEG1 and EEG2 databases showed functional connectivity and multiplex organisation measures, respectively, as relevant subgroups of parameters. Finally, commonly relevant features across databases were selected using SHAP to generate the neurophysiological fingerprints of AD and MCI. This study highlights the relevance of different MEG and EEG parameters in characterising neurological pathologies. The proposed framework, based on MEG and EEG, can be used to generate interpretable, robust, and accurate neurophysiological fingerprints of AD and MCI.
{"title":"Integrative and interpretable framework to unveil the neurophysiological fingerprint of Alzheimer’s disease and mild cognitive impairment: A machine learning-SHAP approach","authors":"Víctor Gutiérrez-de Pablo , María Herrero-Tudela , Marina Sandonís-Fernández , Jesús Poza , Aarón Maturana-Candelas , Víctor Rodríguez-González , Miguel Ángel Tola-Arribas , Mónica Cano , Hideyuki Hoshi , Yoshihito Shigihara , Roberto Hornero , Carlos Gómez","doi":"10.1016/j.bbe.2025.05.011","DOIUrl":"10.1016/j.bbe.2025.05.011","url":null,"abstract":"<div><div>Dementia and mild cognitive impairment (MCI) due to Alzheimer’s disease (AD) are neurological pathologies associated with disruptions in brain electromagnetic activity, typically studied using magnetoencephalography (MEG) and electroencephalography (EEG). To quantify diverse brain properties, different families of parameters can be computed from MEG and EEG (i.e., spectral, non-linear, morphological, functional connectivity, or network structure and organisation). However, studying these characteristics separately overlooks the complex nature of brain activity. Integrative frameworks can be useful to unveil the intricate neurophysiological fingerprint, as well as to characterise pathological conditions comprehensively. To that purpose, data fusion methodologies are crucial, despite their interpretational challenges. In this study, Machine Learning (ML) models were trained to discriminate between groups of severity, whereas the SHapley Additive eXplanations (SHAP) algorithm was afterwards utilised to assess the relevance of the input characteristics into the output classification. Three databases were analysed: MEG (55 healthy controls, HC, 42 MCI patients, and 86 AD patients), EEG1 (51 HC, 52 MCI, and 100 AD), and EEG2 (45 HC, 69 MCI, and 82 AD). The best results for the three-class classification problem were obtained by Gradient Boosting for the MEG database: 3-class Cohen’s kappa coefficient of 0.5452 and accuracy of 72.63 %. Afterwards, using SHAP on Gradient Boosting, it has been shown that spectral features were identified as highly relevant across all databases. Furthermore, morphology measures presented high relevance for the MEG database, whereas EEG1 and EEG2 databases showed functional connectivity and multiplex organisation measures, respectively, as relevant subgroups of parameters. Finally, commonly relevant features across databases were selected using SHAP to generate the neurophysiological fingerprints of AD and MCI. This study highlights the relevance of different MEG and EEG parameters in characterising neurological pathologies. The proposed framework, based on MEG and EEG, can be used to generate interpretable, robust, and accurate neurophysiological fingerprints of AD and MCI.</div></div>","PeriodicalId":55381,"journal":{"name":"Biocybernetics and Biomedical Engineering","volume":"45 3","pages":"Pages 438-450"},"PeriodicalIF":5.3,"publicationDate":"2025-06-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144279921","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Diffusion tensor imaging (DTI) was used to observe degeneration processes at the microstructural and functional levels in the brains of patients with Parkinson’s disease (PD). Two tensor-based unscented Kalman filter (UKF) was used for analyses of eight regions: the substantia nigra, putamen, caudate nucleus, globus pallidus, primary motor cortex, preprimary motor cortex, supplementary motor area (SMA), presupplementary motor area (pre-SMA), and whole brain of patients with PD (n = 14) and controls (n = 12). We analyzed eight DTI metrics in the entire brain and eight brain regions separately for each hemisphere using univariate and multivariate statistical analysis and their correlation with the clinical parameters. The most affected brain regions in patients with PD were the substantia nigra, pre-SMA, globus pallidus, and caudate nucleus. These results suggest that DTI is an adequate tool for evaluating structural and functional alterations, including inflammation, reduced fiber length, changes in neurite density, axonal growth, demyelination, and axonal damage or loss, in the studied brain regions of patients with PD. The results also revealed a generalized brain degeneration process. In conclusion, DTI can be applied for in vivo studies of the degenerative process and could be considered a complementary method in future studies to improve the accuracy of PD diagnosis.
{"title":"Diffusion tensor imaging technique for studying brain microstructural changes in Parkinson’s disease patients","authors":"Beata Toczylowska , Małgorzata Michałowska , Malgorzata Chalimoniuk , Piotr Ladyzynski , Leszek Krolicki , Urszula Fiszer","doi":"10.1016/j.bbe.2025.05.010","DOIUrl":"10.1016/j.bbe.2025.05.010","url":null,"abstract":"<div><div>Diffusion tensor imaging (DTI) was used to observe degeneration processes at the microstructural and functional levels in the brains of patients with Parkinson’s disease (PD). Two tensor-based unscented Kalman filter (UKF) was used for analyses of eight regions: the substantia nigra, putamen, caudate nucleus, globus pallidus, primary motor cortex, preprimary motor cortex, supplementary motor area (SMA), presupplementary motor area (pre-SMA), and whole brain of patients with PD (n = 14) and controls (n = 12). We analyzed eight DTI metrics in the entire brain and eight brain regions separately for each hemisphere using univariate and multivariate statistical analysis and their correlation with the clinical parameters. The most affected brain regions in patients with PD were the substantia nigra, pre-SMA, globus pallidus, and caudate nucleus. These results suggest that DTI is an adequate tool for evaluating structural and functional alterations, including inflammation, reduced fiber length, changes in neurite density, axonal growth, demyelination, and axonal damage or loss, in the studied brain regions of patients with PD. The results also revealed a generalized brain degeneration process. In conclusion, DTI can be applied for in vivo studies of the degenerative process and could be considered a complementary method in future studies to improve the accuracy of PD diagnosis.</div></div>","PeriodicalId":55381,"journal":{"name":"Biocybernetics and Biomedical Engineering","volume":"45 3","pages":"Pages 426-437"},"PeriodicalIF":5.3,"publicationDate":"2025-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144177851","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-05-28DOI: 10.1016/j.bbe.2025.05.005
Fatemeh Ghorbani, Andry Rakotonirainy, Mohammed Elhenawy
This study examined the temporal dynamics of face perception using event-related potentials (ERPs) to investigate how familiarity and repetition influence early and late stages of face processing. A generalised linear mixed-effects (GLME) model was employed to assess the amplitude and latency of the P100, N170, and N250 ERP components across three stimulus types (famous, non-famous, and scrambled faces), three repetition conditions (first presentation, immediate repeat, delayed repeat), and two brain hemispheres. The P100 component, associated with early visual processing, showed no significant modulation by stimulus familiarity or repetition, suggesting stable perceptual encoding across conditions. In contrast, N170 and N250 amplitudes were significantly affected by repetition, indicating enhanced neural responses during repeated exposure, particularly in the right hemisphere. Latency analyses revealed that N250 component was also sensitive to repetition timing, with delayed repetitions eliciting shorter response time, implying shifts in processing efficiency and memory engagement. Multivariate time-series decoding further demonstrated higher discriminability between scrambled and familiar faces compared to non-famous faces, particularly during first and delayed repeat conditions. Notably, decoding performance declined for immediate repeats, suggesting reduced neural differentiation during short-interval repetition. These findings provide new insights into how repetition and familiarity modulate the neural underpinnings of face perception, emphasizing the role of temporal dynamics and hemispheric specialization in face processing.
{"title":"The effect of familiarity and repetition on neural activity during visual face perception","authors":"Fatemeh Ghorbani, Andry Rakotonirainy, Mohammed Elhenawy","doi":"10.1016/j.bbe.2025.05.005","DOIUrl":"10.1016/j.bbe.2025.05.005","url":null,"abstract":"<div><div>This study examined the temporal dynamics of face perception using event-related potentials (ERPs) to investigate how familiarity and repetition influence early and late stages of face processing. A generalised linear mixed-effects (GLME) model was employed to assess the amplitude and latency of the P100, N170, and N250 ERP components across three stimulus types (famous, non-famous, and scrambled faces), three repetition conditions (first presentation, immediate repeat, delayed repeat), and two brain hemispheres. The P100 component, associated with early visual processing, showed no significant modulation by stimulus familiarity or repetition, suggesting stable perceptual encoding across conditions. In contrast, N170 and N250 amplitudes were significantly affected by repetition, indicating enhanced neural responses during repeated exposure, particularly in the right hemisphere. Latency analyses revealed that N250 component was also sensitive to repetition timing, with delayed repetitions eliciting shorter response time, implying shifts in processing efficiency and memory engagement. Multivariate time-series decoding further demonstrated higher discriminability between scrambled and familiar faces compared to non-famous faces, particularly during first and delayed repeat conditions. Notably, decoding performance declined for immediate repeats, suggesting reduced neural differentiation during short-interval repetition. These findings provide new insights into how repetition and familiarity modulate the neural underpinnings of face perception, emphasizing the role of temporal dynamics and hemispheric specialization in face processing.</div></div>","PeriodicalId":55381,"journal":{"name":"Biocybernetics and Biomedical Engineering","volume":"45 3","pages":"Pages 399-413"},"PeriodicalIF":5.3,"publicationDate":"2025-05-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144154937","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-05-28DOI: 10.1016/j.bbe.2025.04.004
Antonina Pater , Lukasz Roszkowiak , Krzysztof Siemion , Jakub Zak , Karol Deptuch , Anna Korzynska
Population screening in the form of cervical smears was introduced to reduce cervical cancer morbidity. However, the manual evaluation of cervical smears is a labour-intensive and meticulous task. This challenge has led to the development of various computer-aided cell identification methods aimed at improving the quality of evaluations and reducing the time required for slide analysis. These supportive tools for pathologists should adhere to the Bethesda classification system for cell types, facilitating integration with established clinical practices. Despite advances, the automatic classification of clustered cells in cervical smears remains a significant challenge for both automated and semiautomated methods that utilize classical image processing and deep learning techniques.
This study introduces a novel method for fragmenting clusters to improve the classification of clustered cells in digital images of Papanicolaou smears. The proposed method integrates explainable AI and marker-guided watershed segmentation ensuring both interpretability and versatility of an overall pipeline for diagnostician support in smear evaluation. Validation of the method was performed on a recently developed Papanicolaou cytology Bialystok dataset, as well as the publicly available CRIC dataset, demonstrating the method’s generalizability across different datasets.
The cell classification pipeline, enhanced by the proposed declustering method, improved the overall harmonic mean of recall and precision (F1 score) by 13.27 percentage points compared with the results obtained without this additional processing. The improvement in classifying the most critical cell type according to the Bethesda system (HSIL cell class) was even more significant, with an increase of 35.72 percentage points compared with classifying the entire cluster.
{"title":"Cytoplasm and nuclei as a basis for Bethesda cell cluster classification in cervical smears","authors":"Antonina Pater , Lukasz Roszkowiak , Krzysztof Siemion , Jakub Zak , Karol Deptuch , Anna Korzynska","doi":"10.1016/j.bbe.2025.04.004","DOIUrl":"10.1016/j.bbe.2025.04.004","url":null,"abstract":"<div><div>Population screening in the form of cervical smears was introduced to reduce cervical cancer morbidity. However, the manual evaluation of cervical smears is a labour-intensive and meticulous task. This challenge has led to the development of various computer-aided cell identification methods aimed at improving the quality of evaluations and reducing the time required for slide analysis. These supportive tools for pathologists should adhere to the Bethesda classification system for cell types, facilitating integration with established clinical practices. Despite advances, the automatic classification of clustered cells in cervical smears remains a significant challenge for both automated and semiautomated methods that utilize classical image processing and deep learning techniques.</div><div>This study introduces a novel method for fragmenting clusters to improve the classification of clustered cells in digital images of Papanicolaou smears. The proposed method integrates explainable AI and marker-guided watershed segmentation ensuring both interpretability and versatility of an overall pipeline for diagnostician support in smear evaluation. Validation of the method was performed on a recently developed Papanicolaou cytology Bialystok dataset, as well as the publicly available CRIC dataset, demonstrating the method’s generalizability across different datasets.</div><div>The cell classification pipeline, enhanced by the proposed declustering method, improved the overall harmonic mean of recall and precision (F1 score) by 13.27 percentage points compared with the results obtained without this additional processing. The improvement in classifying the most critical cell type according to the Bethesda system (HSIL cell class) was even more significant, with an increase of 35.72 percentage points compared with classifying the entire cluster.</div></div>","PeriodicalId":55381,"journal":{"name":"Biocybernetics and Biomedical Engineering","volume":"45 3","pages":"Pages 414-425"},"PeriodicalIF":5.3,"publicationDate":"2025-05-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144166979","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-05-24DOI: 10.1016/j.bbe.2025.05.006
Lei Sun , Xiang Geng , Yichong Wang , Daoyuan Wang , Xin Ma , Gangmin Ning
Diabetic foot is a kind of common complications in diabetic patients, manifesting as ischemic lesions in plantar tissue. The plantar tissue ischemia can be estimated by the perfusion of plantar blood flow. Currently, it lacks effective methods to monitor and assess plantar blood flow perfusion in practice and detailed information is not available. This study intends to address this problem from the view of mathematical modelling. It aims to provide quantitative approach to investigate the risk of ischemia in diabetic foot by simulating the hemodynamic status along the plantar vascular network. A simplified topological structure map of the main plantar arteries was established considering the anatomical structure and physical data of each vessel. Aided by bond graph technology, a Poiseuille’s law based hemodynamic model of the plantar vessels was constructed. The model simulated the blood flow distribution in the feet of healthy individuals and diabetic patients. Compared with healthy individuals, diabetic patients show increased plantar blood pressure and decreased blood supply. The reduction in blood flow on the lateral side of the foot is greater than on the medial side in diabetic patients. The study marks the first application of vascular modeling in the plantar region, providing a quantitative method for investigating foot hemodynamics in diabetic patients and offering new potential for assessing the progression of diabetic foot and the risk of tissue ulceration.
{"title":"Hemodynamic simulation for plantar arteries by anatomical structure based mathematical model","authors":"Lei Sun , Xiang Geng , Yichong Wang , Daoyuan Wang , Xin Ma , Gangmin Ning","doi":"10.1016/j.bbe.2025.05.006","DOIUrl":"10.1016/j.bbe.2025.05.006","url":null,"abstract":"<div><div>Diabetic foot is a kind of common complications in diabetic patients, manifesting as ischemic lesions in plantar tissue. The plantar tissue ischemia can be estimated by the perfusion of plantar blood flow. Currently, it lacks effective methods to monitor and assess plantar blood flow perfusion in practice and detailed information is not available. This study intends to address this problem from the view of mathematical modelling. It aims to provide quantitative approach to investigate the risk of ischemia in diabetic foot by simulating the hemodynamic status along the plantar vascular network. A simplified topological structure map of the main plantar arteries was established considering the anatomical structure and physical data of each vessel. Aided by bond graph technology, a Poiseuille’s law based hemodynamic model of the plantar vessels was constructed. The model simulated the blood flow distribution in the feet of healthy individuals and diabetic patients. Compared with healthy individuals, diabetic patients show increased plantar blood pressure and decreased blood supply. The reduction in blood flow on the lateral side of the foot is greater than on the medial side in diabetic patients. The study marks the first application of vascular modeling in the plantar region, providing a quantitative method for investigating foot hemodynamics in diabetic patients and offering new potential for assessing the progression of diabetic foot and the risk of tissue ulceration.</div></div>","PeriodicalId":55381,"journal":{"name":"Biocybernetics and Biomedical Engineering","volume":"45 3","pages":"Pages 390-398"},"PeriodicalIF":5.3,"publicationDate":"2025-05-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144131213","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-05-21DOI: 10.1016/j.bbe.2025.05.007
Judyta Sobiech, Teodor Buchner
The phenomenon of synchronization, which appears during the coupling of two or more oscillators, is often observed in nature and many fields of science. The case of several rhythms is particularly interesting, as the synchronizations may compete. In this work, we report the occurrence of the phase coupling between three biological rhythms: heart rate, respiration and locomotor rhythm. We analyzed its appearance during cycling on a stationary bike with stepwise increasing speed. Presence of synchronization was investigated using the synchrogram in order to detect phase-locking epochs.
Study group consisted of 21 young healthy volunteers (15 women and 6 men). Specific subgroups were also defined.
During exercise synchronization between all rhythms appeared in all subjects. Cardiolocomotor synchronization occurred significantly more often, which underlines the role of cardiac entrainment to phasic baroreflex input. In cardiolocomotor synchronization we observed a direct phase coupling (1:1) lasting at least 20 s. It occurred in 76 % of subjects, which manifests a dominant role of this type of synchronization. Additionally, we observed subjects, in whom epochs of simultaneous synchronization of all rhythms appeared.
Total duration of synchronization was significantly larger in men than in women and in less experienced cyclists. No significant difference in age or BMI was observed, however the group was relatively uniform. During the post-exercise epoch, cardiorespiratory synchronization appeared in nine subjects. Due to its frequent occurence, synchronization is a good candidate for an objective estimator of physical stress, applicable in diagnostic, therapeutical and lifestyle scenarios.
{"title":"Rhythmic physical exercise as a protocol to study synchronization phenomenon","authors":"Judyta Sobiech, Teodor Buchner","doi":"10.1016/j.bbe.2025.05.007","DOIUrl":"10.1016/j.bbe.2025.05.007","url":null,"abstract":"<div><div>The phenomenon of synchronization, which appears during the coupling of two or more oscillators, is often observed in nature and many fields of science. The case of several rhythms is particularly interesting, as the synchronizations may compete. In this work, we report the occurrence of the phase coupling between three biological rhythms: heart rate, respiration and locomotor rhythm. We analyzed its appearance during cycling on a stationary bike with stepwise increasing speed. Presence of synchronization was investigated using the synchrogram in order to detect phase-locking epochs.</div><div>Study group consisted of 21 young healthy volunteers (15 women and 6 men). Specific subgroups were also defined.</div><div>During exercise synchronization between all rhythms appeared in all subjects. Cardiolocomotor synchronization occurred significantly more often, which underlines the role of cardiac entrainment to phasic baroreflex input. In cardiolocomotor synchronization we observed a direct phase coupling (1:1) lasting at least 20 s. It occurred in 76 % of subjects, which manifests a dominant role of this type of synchronization. Additionally, we observed subjects, in whom epochs of simultaneous synchronization of all rhythms appeared.</div><div>Total duration of synchronization was significantly larger in men than in women and in less experienced cyclists. No significant difference in age or BMI was observed, however the group was relatively uniform. During the post-exercise epoch, cardiorespiratory synchronization appeared in nine subjects. Due to its frequent occurence, synchronization is a good candidate for an objective estimator of physical stress, applicable in diagnostic, therapeutical and lifestyle scenarios.</div></div>","PeriodicalId":55381,"journal":{"name":"Biocybernetics and Biomedical Engineering","volume":"45 3","pages":"Pages 380-389"},"PeriodicalIF":5.3,"publicationDate":"2025-05-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144107325","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Although the optic nerve head (ONH) has demonstrated a significant viscoelastic response to changes in intraocular pressure (IOP), existing computational models of ONH biomechanics have yet to fully account for these viscoelastic properties. In this study, we introduce and evaluate a mesh-free beam-in-solid coupling algorithm to model the complex viscoelastic behavior of anisotropic collagen fibers within a viscoelastic scleral matrix. We also incorporated viscoelastic formulations for the retina, lamina cribrosa, and optic nerve to provide a more comprehensive understanding of ONH tissue mechanics. We compared the biomechanics of the ONH resulting from hyperelastic and viscoelastic scleral material formulations using an eye-specific finite element model of the posterior human eye. This model integrates the detailed 3D microstructure of the load-bearing lamina cribrosa, including interspersed laminar neural tissues, as well as the heterogeneous, anisotropic behavior of the collagenous sclera and pia. The viscoelastic material properties were validated against published experimental tensile tests of human scleral and retinal tissue samples. Simulations of ONH biomechanical responses were conducted by applying changes in IOP and cerebrospinal fluid pressure (CSFP) typical of body position transitions, such as moving from sitting to supine, over a 250 ms period. In both simulations, the ONH tissues exhibited greater stresses and strains in the supine position compared to sitting, as anticipated. The laminar surface showed posterior deformation (+6 µm) during the transition from sitting to supine when using the hyperelastic material model, whereas it deformed anteriorly (−5.7 µm) with the viscoelastic model. Furthermore, the radial scleral canal expansion at the anterior laminar insertion was significantly smaller in the viscoelastic formulation (9 µm) compared to the hyperelastic formulation (19.8 µm). All results aligned with experimental observations. While the stresses, strains, and deformations remained within physiological ranges for both models, there were substantial differences between the two formulations, particularly in terms of deformation. Improving the accuracy of material formulations in ONH models is expected to enhance our understanding of ONH biomechanics. However, further experimental validation is needed to confirm these results and strengthen their applicability.
{"title":"Viscoelastic modeling of optic nerve head biomechanics: Effects of intraocular and cerebrospinal fluid pressure","authors":"Alireza Karimi , Reza Razaghi , Seyed Mohammadali Rahmati , J. Crawford Downs","doi":"10.1016/j.bbe.2025.05.008","DOIUrl":"10.1016/j.bbe.2025.05.008","url":null,"abstract":"<div><div>Although the optic nerve head (ONH) has demonstrated a significant viscoelastic response to changes in intraocular pressure (IOP), existing computational models of ONH biomechanics have yet to fully account for these viscoelastic properties. In this study, we introduce and evaluate a mesh-free beam-in-solid coupling algorithm to model the complex viscoelastic behavior of anisotropic collagen fibers within a viscoelastic scleral matrix. We also incorporated viscoelastic formulations for the retina, lamina cribrosa, and optic nerve to provide a more comprehensive understanding of ONH tissue mechanics. We compared the biomechanics of the ONH resulting from hyperelastic and viscoelastic scleral material formulations using an eye-specific finite element model of the posterior human eye. This model integrates the detailed 3D microstructure of the load-bearing lamina cribrosa, including interspersed laminar neural tissues, as well as the heterogeneous, anisotropic behavior of the collagenous sclera and pia. The viscoelastic material properties were validated against published experimental tensile tests of human scleral and retinal tissue samples. Simulations of ONH biomechanical responses were conducted by applying changes in IOP and cerebrospinal fluid pressure (CSFP) typical of body position transitions, such as moving from sitting to supine, over a 250 ms period. In both simulations, the ONH tissues exhibited greater stresses and strains in the supine position compared to sitting, as anticipated. The laminar surface showed posterior deformation (+6 µm) during the transition from sitting to supine when using the hyperelastic material model, whereas it deformed anteriorly (−5.7 µm) with the viscoelastic model. Furthermore, the radial scleral canal expansion at the anterior laminar insertion was significantly smaller in the viscoelastic formulation (9 µm) compared to the hyperelastic formulation (19.8 µm). All results aligned with experimental observations. While the stresses, strains, and deformations remained within physiological ranges for both models, there were substantial differences between the two formulations, particularly in terms of deformation. Improving the accuracy of material formulations in ONH models is expected to enhance our understanding of ONH biomechanics. However, further experimental validation is needed to confirm these results and strengthen their applicability.</div></div>","PeriodicalId":55381,"journal":{"name":"Biocybernetics and Biomedical Engineering","volume":"45 3","pages":"Pages 357-379"},"PeriodicalIF":5.3,"publicationDate":"2025-05-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144090192","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}