ITree: a user-driven tool for interactive decision-making with classification trees.

MOTIVATION: ITree is an intuitive web tool for the manual, semi-automatic, and automatic induction of decision trees. It enables interactive modifications of tree structures and incorporates Relative Expression Analysis for detecting complex patterns in high-throughput molecular data. This makes ITree a versatile tool for both research and education in biomedical data analysis. RESULTS: The tool allows users to instantly see the effects of modifications on decision trees, with updates to predictions and statistics displayed in real time, facilitating a deeper understanding of data classification processes. AVAILABILITY AND IMPLEMENTATION: Available online at https://itree.wi.pb.edu.pl. Source code and documentation are hosted on GitHub at https://github.com/hsokolowski/iTree and in supplement.

GPU-accelerated lung CT segmentation based on level sets and texture analysis.

This paper presents a novel semi-automatic method for lung segmentation in thoracic CT datasets. The fully three-dimensional algorithm is based on a level set representation of an active surface and integrates texture features to improve its robustness. The method's performance is enhanced by the graphics processing unit (GPU) acceleration. The segmentation process starts with a manual initialisation of 2D contours on a few representative slices of the analysed volume. Next, the starting regions for the active surface are generated according to the probability maps of texture features. The active surface is then evolved to give the final segmentation result. The recent implementation employs features based on grey-level co-occurrence matrices and Gabor filters. The algorithm was evaluated on real medical imaging data from the LCTCS 2017 challenge. The results were also compared with the outcomes of other segmentation methods. The proposed approach provided high segmentation accuracy while offering very competitive performance.

A comparison of different machine-learning techniques for the selection of a panel of metabolites allowing early detection of brain tumors.

Metabolomics combined with machine learning methods (MLMs), is a powerful tool for searching novel diagnostic panels. This study was intended to use targeted plasma metabolomics and advanced MLMs to develop strategies for diagnosing brain tumors. Measurement of 188 metabolites was performed on plasma samples collected from 95 patients with gliomas (grade I-IV), 70 with meningioma, and 71 healthy individuals as a control group. Four predictive models to diagnose glioma were prepared using 10 MLMs and a conventional approach. Based on the cross-validation results of the created models, the F1-scores were calculated, then obtained values were compared. Subsequently, the best algorithm was applied to perform five comparisons involving gliomas, meningiomas, and controls. The best results were obtained using the newly developed hybrid evolutionary heterogeneous decision tree (EvoHDTree) algorithm, which was validated using Leave-One-Out Cross-Validation, resulting in an F1-score for all comparisons in the range of 0.476-0.948 and the area under the ROC curves ranging from 0.660 to 0.873. Brain tumor diagnostic panels were constructed with unique metabolites, which reduces the likelihood of misdiagnosis. This study proposes a novel interdisciplinary method for brain tumor diagnosis based on metabolomics and EvoHDTree, exhibiting significant predictive coefficients.

Exploring microRNAs as predictive biomarkers for type 2 diabetes mellitus remission after sleeve gastrectomy: A pilot study.

OBJECTIVE: This study aimed to evaluate microRNAs (miRNAs) as predictive biomarkers for type 2 diabetes (T2D) remission 12 months after sleeve gastrectomy (SG). METHODS: A total of 179 serum miRNAs were profiled, and 26 clinical variables were collected from 46 patients. Two patients were later excluded because of hemolysis, and six patients with unclear remission status were set aside to evaluate the prediction models. The remaining 38 patients were included for model building. Variable selection was done using different approaches, including Least Absolute Shrinkage and Selection Operator (LASSO). Prediction models were then developed using LASSO and assessed in the validation set. RESULTS: A total of 26 out of 38 patients achieved T2D remission 12 months after SG. The prediction model with only clinical variables misclassified two patients, which were correctly classified using miRNAs. Two miRNA-only models achieved an accuracy of one but performed poorly for the validation set. The best miRNA model was a mixed model (accuracy: 0.974) containing four miRNAs (hsa-miR-32-5p, hsa-miR-382-5p, hsa-miR-1-3p, and hsa-miR-21-5p) and four clinical variables (T2D medication, sex, age, and fasting blood glucose). These miRNAs are involved in pathways related to obesity and insulin resistance. CONCLUSIONS: This study suggests that four serum miRNAs might be predictive biomarkers for T2D remission 12 months after SG, but further validation studies are needed.

Towards a patient-specific hepatic arterial modeling for microspheres distribution optimization in SIRT protocol.

Selective internal radiation therapy (SIRT) using Yttrium-90 loaded glass microspheres injected in the hepatic artery is an emerging, minimally invasive therapy of liver cancer. A personalized intervention can lead to high concentration dose in the tumor, while sparing the surrounding parenchyma. We propose a computational model for patient-specific simulation of entire hepatic arterial tree, based on liver, tumors, and arteries segmentation on patient's tomography. Segmentation of hepatic arteries down to a diameter of 0.5 mm is semi-automatically performed on 3D cone-beam CT angiography. The liver and tumors are extracted from CT-scan at portal phase by an active surface method. Once the images are registered through an automatic multimodal registration, extracted data are used to initialize a numerical model simulating liver vascular network. The model creates successive bifurcations from given principal vessels, observing Poiseuille's and matter conservation laws. Simulations provide a coherent reconstruction of global hepatic arterial tree until vessel diameter of 0.05 mm. Microspheres distribution under simple hypotheses is also quantified, depending on injection site. The patient-specific character of this model may allow a personalized numerical approximation of microspheres final distribution, opening the way to clinical optimization of catheter placement for tumor targeting.

In silico modeling of magnetic resonance flow imaging in complex vascular networks.

The paper presents a computational model of magnetic resonance (MR) flow imaging. The model consists of three components. The first component is used to generate complex vascular structures, while the second one provides blood flow characteristics in the generated vascular structures by the lattice Boltzmann method. The third component makes use of the generated vascular structures and flow characteristics to simulate MR flow imaging. To meet computational demands, parallel algorithms are applied in all the components. The proposed approach is verified in three stages. In the first stage, experimental validation is performed by an in vitro phantom. Then, the simulation possibilities of the model are shown. Flow and MR flow imaging in complex vascular structures are presented and evaluated. Finally, the computational performance is tested. Results show that the model is able to reproduce flow behavior in large vascular networks in a relatively short time. Moreover, simulated MR flow images are in accordance with the theoretical considerations and experimental images. The proposed approach is the first such an integrative solution in literature. Moreover, compared to previous works on flow and MR flow imaging, this approach distinguishes itself by its computational efficiency. Such a connection of anatomy, physiology and image formation in a single computer tool could provide an in silico solution to improving our understanding of the processes involved, either considered together or separately.

Evolutionary approach for relative gene expression algorithms.

A Relative Expression Analysis (RXA) uses ordering relationships in a small collection of genes and is successfully applied to classiffication using microarray data. As checking all possible subsets of genes is computationally infeasible, the RXA algorithms require feature selection and multiple restrictive assumptions. Our main contribution is a specialized evolutionary algorithm (EA) for top-scoring pairs called EvoTSP which allows finding more advanced gene relations. We managed to unify the major variants of relative expression algorithms through EA and introduce weights to the top-scoring pairs. Experimental validation of EvoTSP on public available microarray datasets showed that the proposed solution significantly outperforms in terms of accuracy other relative expression algorithms and allows exploring much larger solution space.

Multi-test decision tree and its application to microarray data classification.

OBJECTIVE: The desirable property of tools used to investigate biological data is easy to understand models and predictive decisions. Decision trees are particularly promising in this regard due to their comprehensible nature that resembles the hierarchical process of human decision making. However, existing algorithms for learning decision trees have tendency to underfit gene expression data. The main aim of this work is to improve the performance and stability of decision trees with only a small increase in their complexity. METHODS: We propose a multi-test decision tree (MTDT); our main contribution is the application of several univariate tests in each non-terminal node of the decision tree. We also search for alternative, lower-ranked features in order to obtain more stable and reliable predictions. RESULTS: Experimental validation was performed on several real-life gene expression datasets. Comparison results with eight classifiers show that MTDT has a statistically significantly higher accuracy than popular decision tree classifiers, and it was highly competitive with ensemble learning algorithms. The proposed solution managed to outperform its baseline algorithm on 14 datasets by an average 6%. A study performed on one of the datasets showed that the discovered genes used in the MTDT classification model are supported by biological evidence in the literature. CONCLUSION: This paper introduces a new type of decision tree which is more suitable for solving biological problems. MTDTs are relatively easy to analyze and much more powerful in modeling high dimensional microarray data than their popular counterparts.

Computational modeling of MR flow imaging by the lattice Boltzmann method and Bloch equation.

In this work, a computational model of magnetic resonance (MR) flow imaging is proposed. The first model component provides fluid dynamics maps by applying the lattice Boltzmann method. The second one uses the flow maps and couples MR imaging (MRI) modeling with a new magnetization transport algorithm based on the Eulerian coordinate approach. MRI modeling is based on the discrete time solution of the Bloch equation by analytical local magnetization transformations (exponential scaling and rotations). Model is validated by comparison of experimental and simulated MR images in two three-dimensional geometries (straight and U-bend tubes) with steady flow under comparable conditions. Two-dimensional geometries, presented in literature, were also tested. In both cases, a good agreement is observed. Quantitative analysis shows in detail the model accuracy. Computational time is noticeably lower to prior works. These results demonstrate that the discrete time solution of Bloch equation coupled with the new magnetization transport algorithm naturally incorporates flow influence in MRI modeling. As a result, in the proposed model, no additional mechanism (unlike in prior works) is needed to consider flow artifacts, which implies its easy extensibility. In combination with its low computational complexity and efficient implementation, the model could have a potential application in study of flow disturbances (in MRI) in various conditions and in different geometries.

Vascular system modeling in parallel environment - distributed and shared memory approaches.

This paper presents two approaches in parallel modeling of vascular system development in internal organs. In the first approach, new parts of tissue are distributed among processors and each processor is responsible for perfusing its assigned parts of tissue to all vascular trees. Communication between processors is accomplished by passing messages, and therefore, this algorithm is perfectly suited for distributed memory architectures. The second approach is designed for shared memory machines. It parallelizes the perfusion process during which individual processing units perform calculations concerning different vascular trees. The experimental results, performed on a computing cluster and multicore machines, show that both algorithms provide a significant speedup.

Multiscale model of liver DCE-MRI towards a better understanding of tumor complexity.

The use of quantitative imaging for the characterization of hepatic tumors in magnetic resonance imaging (MRI) can improve the diagnosis and therefore the treatment of these life-threatening tumors. However, image parameters remain difficult to interpret because they result from a mixture of complex processes related to pathophysiology and to acquisition. These processes occur at variable spatial and temporal scales. We propose a multiscale model of liver dynamic contrast-enhanced (DCE) MRI in order to better understand the tumor complexity in images. Our design couples a model of the organ (tissue and vasculature) with a model of the image acquisition. At the macroscopic scale, vascular trees take a prominent place. Regarding the formation of MRI images, we propose a distributed model of parenchymal biodistribution of extracellular contrast agents. Model parameters can be adapted to simulate the tumor development. The sensitivity of the multiscale model of liver DCE-MRI was studied through observations of the influence of two physiological parameters involved in carcinogenesis (arterial flow and capillary permeability) on its outputs (MRI images at arterial and portal phases). Finally, images were simulated for a set of parameters corresponding to the five stages of hepatocarcinogenesis (from regenerative nodules to poorly differentiated HepatoCellular Carcinoma).

Physiological modeling of tumor-affected renal circulation.

One way of gaining insight into what can be observed in medical images is through physiological modeling. For instance, anatomical and functional modifications occur in the organ during the appearance and the growth of a tumor. Some of these changes concern the vascularization. We propose a computational model of tumor-affected renal circulation that represents the local heterogeneity of different parts of the kidney (cortex, medulla). We present a simulation of vascular modifications related to vessel structure, geometry, density, and blood flow in case of renal cell carcinoma. We also use our model to simulate computed tomography scans of a kidney affected by the renal cell carcinoma, at two acquisition times after injection of a contrast product. This framework, based on a physiological model of the organ and physical model of medical image acquisition, offers an opportunity to help radiologists in their diagnostic tasks. This includes the possibility of linking image descriptors with physiological perturbations and markers of pathological processes.

Coupling texture analysis and physiological modeling for liver dynamic MRI interpretation.

We coupled our physiological model of the liver, to a MRI simulator (SIMRI) in order to find image markers of the tumor growth. Some pathological modifications related to the development of Hepatocellular carcinoma are simulated (flows, permeability, vascular density). Corresponding images simulated at typical acquisition phases (arterial, portal) are compared to real images. The evolution of some textural features with arterial flow is also presented.

Simulation of biphasic CT findings in hepatic cellular carcinoma by a two-level physiological model.

In this paper, we present a two-level physiological model that is able to reflect morphology and function of vascular networks, in clinical images. Our approach results from the combination of a macroscopic model, providing simulation of the growth and pathological modifications of vascular network, and a microvascular model, based on compartmental approach, which simulates blood and contrast medium transfer through capillary walls. The two-level model is applied to generate biphasic computed tomography of hepatocellular carcinoma. A contrast-enhanced sequence of simulated images is acquired, and enhancement curves extracted from normal and tumoral regions are compared to curves obtained from in vivo images. The model offers the potential of finding early indicators of disease in clinical vascular images.

[Evaluation of adiponectin and TNFalpha genes expression in women with gestational diabetes. Preliminary results]. / Ocena ekspresji genów adiponektyny i TNFalpha w tkance tluszczowej i lozysku kobiet z cukrzyca ciazowa, badania wstepne.

OBJECTIVES: Recent studies suggest the essential role of different cytokines realised from adipose tissue in pathogenesis of gestational diabetes. The aim of the study was evaluation of adiponectin (diabetes development protective factor) and TNFalpha (one of the most important insulin resistance mediator) genes expression in maternal visceral and subcutaneous adipose tissue as well as placental tissue. MATERIAL AND METHODS: The study group consists of patients with gestational diabetes, healthy pregnant glucose tolerant women represented the control group. Tissue samples--placental tissue, visceral and subcutaneous adipose tissue--were obtained from patients who were undergoing ceasarean section. RT-PCR technique was performed for evaluation of adoponectin and TNF alpha genes expression. RESULTS: Our study showed decreased adiponectin gene expression and increased TNFalpha gene expression in visceral adipose tissue in pregnant women with gestational diabetes. An expression of adiponectin gene in diabetic placental tissue was not observed. We also noticed slight increase of TNF alpha gene expressionin placenta of diabetic cases. CONCLUSIONS: Decreased adiponectin and increased TNFalpha genes expression in adipose tissue of pregnant women with gestational diabetes seems to play a significant role in insulin resistance appearance and can lead to development of diabetes in pregnancy.

Physiologically based modeling of 3-D vascular networks and CT scan angiography.

In this paper, a model-based approach to medical image analysis is presented. It is aimed at understanding the influence of the physiological (related to tissue) and physical (related to image modality) processes underlying the image content. This methodology is exemplified by modeling first, the liver and its vascular network, and second, the standard computed tomography (CT) scan acquisition. After a brief survey on vascular modeling literature, a new method, aimed at the generation of growing three-dimensional vascular structures perfusing the tissue, is described. A solution is proposed in order to avoid intersections among vessels belonging to arterial and/or venous trees, which are physiologically connected. Then it is shown how the propagation of contrast material leads to simulate time-dependent sequences of enhanced liver CT slices.

Fast algorithm for 3-D vascular tree modeling.

In this short paper, accelerated three-dimensional computer simulations of vascular trees development, preserving physiological and haemodynamic features, are reported. The new computation schemes deal: (i). with the geometrical optimization of each newly created bifurcation; and (ii). with the recalculation of blood pressures and radii of vessels in the whole tree. A significant decrease of the computation time is obtained by replacing the global optimization by the fast updating algorithm allowing more complex structure to be simulated. A comparison between the new algorithms and the previous one is illustrated through the hepatic arterial tree.

Hepatic tumor enhancement in computed tomography: combined models of liver perfusion and dynamic imaging.

The objective of this study is to show how computational modeling can be used to increase our understanding of liver enhancement in dynamic computer tomography. It relies on two models: (1). a vascular model, based on physiological rules, is used to generate the 3D hepatic vascular network; (2). the physical process of CT acquisition allows to synthesize timed-stamped series of images, aimed at tracking the propagation of a contrast material through the vessel network and the parenchyma. The coupled models are used to simulate the enhancement of a hyper-vascular tumor at different acquisition times, showing a maximum conspicuity during the arterial phase.