Evaluation of models of sequestration flow in coronary arteries-Physiology versus anatomy?

BACKGROUND

Myocardial ischaemia results from insufficient coronary blood flow. Computed virtual fractional flow reserve (vFFR) allows quantification of proportional flow loss without the need for invasive pressure-wire testing. In the current study, we describe a novel, conductivity model of side branch flow, referred to as 'leak'. This leak model is a function of taper and local pressure, the latter of which may change radically when focal disease is present. This builds upon previous techniques, which either ignore side branch flow, or rely purely on anatomical factors. This study aimed to describe a new, conductivity model of side branch flow and compare this with established anatomical models.

METHODS AND RESULTS

The novel technique was used to quantify vFFR, distal absolute flow (Qd) and microvascular resistance (CMVR) in 325 idealised 1D models of coronary arteries, modelled from invasive clinical data. Outputs were compared to an established anatomical model of flow. The conductivity model correlated and agreed with the reference model for vFFR (r = 0.895, p < 0.0001; ＋0.02, 95% CI 0.00 to ＋ 0.22), Qd (r = 0.959, p < 0.0001; -5.2 mL/min, 95% CI -52.2 to ＋13.0) and CMVR (r = 0.624, p < 0.0001; ＋50 Woods Units, 95% CI -325 to ＋2549).

CONCLUSION

Agreement between the two techniques was closest for vFFR, with greater proportional differences seen for Qd and CMVR. The conductivity function assumes vessel taper was optimised for the healthy state and that CMVR was not affected by local disease. The latter may be addressed with further refinement of the technique or inferred from complementary image data. The conductivity technique may represent a refinement of current techniques for modelling coronary side-branch flow. Further work is needed to validate the technique against invasive clinical data.

Effect of turbulence and viscosity models on wall shear stress derived biomarkers for aorta simulations

Ascending aorta simulations provide insight into patient-specific hemodynamic conditions. Numerous studies have assessed fluid biomarkers which show a potential to aid clinicians in the diagnosis process. Unfortunately, there exists a large disparity in the computational methodology used to model turbulence and viscosity. Recognizing this disparity, some authors focused on analysing the influence of either the turbulence or viscosity models on the biomarkers in order to quantify the importance of these model choices. However, no analysis has yet been done on their combined effect. In order to fully understand and quantify the effect of the computational methodology, an assessment of the combined effect of turbulence and viscosity model choice was performed. Our results show that (1) non-Newtonian viscosity has greater impact (2.9-5.0%) on wall shear stress than Large Eddy Simulation turbulence modelling (0.1-1.4%), (2) the contribution of non-Newtonian viscosity is amplified when combined with a subgrid-scale turbulence model, (3) wall shear stress is underestimated when considering Newtonian viscosity by 2.9-5.0% and (4) cycle-to-cycle variability can impact the results as much as the numerical model if insufficient cycles are performed. These results demonstrate that, when assessing the effect of computational methodologies, the resultant combined effect of the different modelling assumptions differs from the aggregated effect of the isolated modifications. Accurate aortic flow modelling requires non-Newtonian viscosity and Large Eddy Simulation turbulence modelling.

A reduced complexity ECG imaging model for regularized inversion optimization

The resolution of the inverse problem of electrocardiography represents a major interest in the diagnosis and catheter-based therapy of cardiac arrhythmia. In this context, the ability to simulate several cardiac electrical behaviors was crucial for evaluating and comparing the performance of inversion methods. For this application, existing models are either too complex or do not produce realistic cardiac patterns. In this work, a low-resolution heart-torso model generating realistic whole heart cardiac mappings and electrocardiograms in healthy and pathological cases is designed. This model was built upon a simplified heart-torso geometry and implements the monodomain formalism by using the finite element method. In addition, a model reduction step through a sensitivity analysis was proposed where parameters were identified using an evolutionary optimization approach. Finally, the study illustrates the usefulness of the proposed model by comparing the performance of different variants of Tikhonov-based inversion methods for the determination of the regularization parameter in healthy, ischemic and ventricular tachycardia scenarios. First, results of the sensitivity analysis show that among 58 parameters only 25 are influent. Note also that the level of influence of the parameters depends on the heart region. Besides, the synthesized electrocardiograms globally present the same characteristic shape compared to the reference once with a correlation value that reaches 88%. Regarding inverse problem, results highlight that only Robust Generalized Cross Validation and Discrepancy Principle provide best performance, with a quasi-perfect success rate for both, and a respective relative error, between the generated electrocardiograms to the reference one, of 0.75 and 0.62.

Calibration of the Mechanical Boundary Conditions for a Patient-Specific Thoracic Aorta Model Including the Heart Motion Effect

OBJECTIVE

We propose a procedure for calibrating 4 parameters governing the mechanical boundary conditions (BCs) of a thoracic aorta (TA) model derived from one patient with ascending aortic aneurysm. The BCs reproduce the visco-elastic structural support provided by the soft tissue and the spine and allow for the inclusion of the heart motion effect.

METHODS

We first segment the TA from magnetic resonance imaging (MRI) angiography and derive the heart motion by tracking the aortic annulus from cine-MRI. A rigid-wall fluid-dynamic simulation is performed to derive the time-varying wall pressure field. We build the finite element model considering patient-specific material properties and imposing the derived pressure field and the motion at the annulus boundary. The calibration, which involves the zero-pressure state computation, is based on purely structural simulations. After obtaining the vessel boundaries from the cine-MRI sequences, an iterative procedure is performed to minimize the distance between them and the corresponding boundaries derived from the deformed structural model. A strongly-coupled fluid-structure interaction (FSI) analysis is finally performed with the tuned parameters and compared to the purely structural simulation.

RESULTS AND CONCLUSION

The calibration with structural simulations allows to reduce maximum and mean distances between image-derived and simulation-derived boundaries from 8.64 mm to 6.37 mm and from 2.24 mm to 1.83 mm, respectively. The maximum root mean square error between the deformed structural and FSI surface meshes is 0.19 mm. This procedure may prove crucial for increasing the model fidelity in replicating the real aortic root kinematics.

Segmentation of the aorta in systolic phase from 4D flow MRI: multi-atlas vs. deep learning

OBJECTIVE

In the management of the aortic aneurysm, 4D flow magnetic resonance Imaging provides valuable information for the computation of new biomarkers using computational fluid dynamics (CFD). However, accurate segmentation of the aorta is required. Thus, our objective is to evaluate the performance of two automatic segmentation methods on the calculation of aortic wall pressure.

METHODS

Automatic segmentation of the aorta was performed with methods based on deep learning and multi-atlas using the systolic phase in the 4D flow MRI magnitude image of 36 patients. Using mesh morphing, isotopological meshes were generated, and CFD was performed to calculate the aortic wall pressure. Node-to-node comparisons of the pressure results were made to identify the most robust automatic method respect to the pressures obtained with a manually segmented model.

RESULTS

Deep learning approach presented the best segmentation performance with a mean Dice similarity coefficient and a mean Hausdorff distance (HD) equal to 0.92+/- 0.02 and 21.02+/- 24.20 mm, respectively. At the global level HD is affected by the performance in the abdominal aorta. Locally, this distance decreases to 9.41+/- 3.45 and 5.82+/- 6.23 for the ascending and descending thoracic aorta, respectively. Moreover, with respect to the pressures from the manual segmentations, the differences in the pressures computed from deep learning were lower than those computed from multi-atlas method.

CONCLUSION

To reduce biases in the calculation of aortic wall pressure, accurate segmentation is needed, particularly in regions with high blood flow velocities. Thus, the deep learning segmen-tation method should be preferred.

Computer-aided shape features extraction and regression models for predicting the ascending aortic aneurysm growth rate

OBJECTIVE

ascending aortic aneurysm growth prediction is still challenging in clinics. In this study, we evaluate and compare the ability of local and global shape features to predict the ascending aortic aneurysm growth.

MATERIAL AND METHODS

70 patients with aneurysm, for which two 3D acquisitions were available, are included. Following segmentation, three local shape features are computed: (1) the ratio between maximum diameter and length of the ascending aorta centerline, (2) the ratio between the length of external and internal lines on the ascending aorta and (3) the tortuosity of the ascending tract. By exploiting longitudinal data, the aneurysm growth rate is derived. Using radial basis function mesh morphing, iso-topological surface meshes are created. Statistical shape analysis is performed through unsupervised principal component analysis (PCA) and supervised partial least squares (PLS). Two types of global shape features are identified: three PCA-derived and three PLS-based shape modes. Three regression models are set for growth prediction: two based on gaussian support vector machine using local and PCA-derived global shape features; the third is a PLS linear regression model based on the related global shape features. The prediction results are assessed and the aortic shapes most prone to growth are identified.

RESULTS

the prediction root mean square error from leave-one-out cross-validation is: 0.112 mm/month, 0.083 mm/month and 0.066 mm/month for local, PCA-based and PLS-derived shape features, respectively. Aneurysms close to the root with a large initial diameter report faster growth.

CONCLUSION

global shape features might provide an important contribution for predicting the aneurysm growth.

Prediction of the 3D shape of the L1 vertebral body from adjacent vertebrae

The aim of treatments of vertebral fractures is the anatomical reduction to restore the physiological biomechanics of the spine and the stabilization of the fracture to allow bone healing. However, the three-dimensional shape of the fractured vertebral body before the fracture is unknown in the clinical setting. Information about the pre-fracture vertebral body shape could help surgeons to select the optimal treatment. The goal of this study was to develop and validate a method based on Singular Value Decomposition (SVD) to predict the shape of the vertebral body of L1 from the shapes of T12 and L2. The geometry of the vertebral bodies of T12, L1 and L2 vertebrae of 40 patients were extracted from CT scans available from the VerSe2020 open-access dataset. Surface triangular meshes of each vertebra were morphed onto a template mesh. The set of vectors with the node coordinates of the morphed T12, L1 and L2 were compressed with SVD and used to build a system of linear equations. This system was used to solve a minimization problem and to reconstruct the shape of L1. A leave-one-out cross-validation was performed. Moreover, the approach was tested against an independent dataset with large osteophytes. The results of the study show a good prediction of the shape of the vertebral body of L1 from the shapes of the two adjacent vertebrae (mean error equal to 0.51 ± 0.11 mm on average, Hausdorff distance equal to 2.11 ± 0.56 mm on average), compared to current CT resolution typically used in the operating room. The error was slightly higher for patients presenting large osteophytes or severe bone degeneration (mean error equal to 0.65 ± 0.10 mm, Hausdorff distance equal to 3.54 ± 1.03 mm). The accuracy of the prediction was significantly better than approximating the shape of the vertebral body of L1 by the shape of T12 or L2. This approach could be used in the future to improve the pre-planning of spine surgeries to treat vertebral fractures.

Patient-specific finite element simulation of peripheral artery percutaneous transluminal angioplasty to evaluate the procedure outcome without stent implantation

The purpose of this work is to present a patient-specific (PS) modeling approach for simulating percutaneous transluminal angioplasty (PTA) endovascular treatment and assessing the balloon sizing influence on short-term outcomes in peripheral arteries, i.e. without stent implantation. Two 3D PS stenosed femoral artery models, one with a dominant calcified atherosclerosis while the other with a lipidic plaque, were generated from pre-operative computed tomography angiography images. Elastoplastic constitutive laws were implemented within the plaque and artery models. Implicit finite element method (FEM) was used to simulate the balloon inflation and deflation for different sizings. Besides vessel strains, results were mainly evaluated in terms of the elastic recoil ratio (ERR) and lumen gain ratio (LGR) attained immediately after PTA. Higher LGR values were shown within the stenosed region of the lipidic patient. Simulated results also showed a direct and quantified correlation between balloon sizing and LGR and ERR for both patients after PTA, with a more significant influence on the lumen gain. The max principal strain values in the outer arterial wall increased at higher balloon sizes during inflation as well, with higher rates of increase when the plaque was calcified. Results show that our model could serve in finding a compromise for each stenosis type: maximizing the achieved lumen gain after PTA, but at the same time without damaging the arterial tissue. The proposed methodology can serve as a step toward a clinical decision support system to improve angioplasty balloon sizing selection prior to the surgery.

"Neuro-Fiber Mapping": An Original Concept of Spinal Cord Neural Network Spatial Targeting Using Live Electrostimulation Mapping to (Re-)Explore the Conus Medullaris Anatomy

Spinal cord (SC) anatomy is often assimilated to a morphologically encapsulated neural entity, but its functional anatomy remains only partially understood. We hypothesized that it could be possible to re-explore SC neural networks by performing live electrostimulation mapping, based on "super-selective" spinal cord stimulation (SCS), originally designed as a therapeutical tool to address chronic refractory pain. As a starting point, we initiated a systematic SCS lead programming approach using live electrostimulation mapping on a chronic refractory perineal pain patient, previously implanted with multicolumn SCS at the level of the conus medullaris (T12-L1). It appeared possible to (re-)explore the classical anatomy of the conus medullaris using statistical correlations of paresthesia coverage mappings, resulting from 165 different electrical configurations tested. We highlighted that sacral dermatomes were not only located more medially but also deeper than lumbar dermatomes at the level of the conus medullaris, in contrast with classical anatomical descriptions of SC somatotopical organization. As we were finally able to find a morphofunctional description of "Philippe-Gombault's triangle" in 19th-century historical textbooks of neuroanatomy, remarkably matching these conclusions, the concept of "neuro-fiber mapping" was introduced.

Validation of a novel numerical model to predict regionalised blood flow in the coronary arteries

Background

Ischaemic heart disease results from insufficient coronary blood flow. Direct measurement of absolute flow (mL/min) is feasible, but has not entered routine clinical practice in most catheterisation laboratories. Interventional cardiologists therefore rely on surrogate markers of flow. Recently, we described a computational fluid dynamics (CFD) method for predicting flow that differentiates inlet, side branch and outlet flows during angiography. In the current study, we evaluate a new method that regionalises flow along the length of the artery.

Methods and Results

Three-dimensional coronary anatomy was reconstructed from angiograms from 20 patients with chronic coronary syndrome. All flows were computed using CFD by applying the pressure gradient to the reconstructed geometry. Side branch flow was modelled as a porous wall boundary. Side branch flow magnitude was based on morphometric scaling laws with two models: a homogenous model with flow loss along the entire arterial length; and a regionalised model with flow proportional to local taper. Flow results were validated against invasive measurements of flow by continuous infusion thermodilution (Coroventis™, Abbott). Both methods quantified flow relative to the invasive measures: homogenous (r 0.47, P 0.006; zero bias; 95% CI -168 to +168 mL/min); regionalised method (r 0.43, P 0.013; zero bias; 95% CI -175 to +175 mL/min).

Conclusions

During angiography and pressure-wire assessment, coronary flow can now be regionalised and differentiated at the inlet, outlet and side branches. The effect of epicardial disease on agreement suggests the model may be best targeted at cases with a stenosis close to side branches.

Assessment of shape-based features ability to predict the ascending aortic aneurysm growth

The current guidelines for the ascending aortic aneurysm (AsAA) treatment recommend surgery mainly according to the maximum diameter assessment. This criterion has already proven to be often inefficient in identifying patients at high risk of aneurysm growth and rupture. In this study, we propose a method to compute a set of local shape features that, in addition to the maximum diameter D, are intended to improve the classification performances for the ascending aortic aneurysm growth risk assessment. Apart from D, these are the ratio DCR between D and the length of the ascending aorta centerline, the ratio EILR between the length of the external and the internal lines and the tortuosity T. 50 patients with two 3D acquisitions at least 6 months apart were segmented and the growth rate (GR) with the shape features related to the first exam computed. The correlation between them has been investigated. After, the dataset was divided into two classes according to the growth rate value. We used six different classifiers with input data exclusively from the first exam to predict the class to which each patient belonged. A first classification was performed using only D and a second with all the shape features together. The performances have been evaluated by computing accuracy, sensitivity, specificity, area under the receiver operating characteristic curve (AUROC) and positive (negative) likelihood ratio LHR+ (LHR-). A positive correlation was observed between growth rate and DCR (r = 0.511, p = 1.3e-4) and between GR and EILR (r = 0.472, p = 2.7e-4). Overall, the classifiers based on the four metrics outperformed the same ones based only on D. Among the diameter-based classifiers, k-nearest neighbours (KNN) reported the best accuracy (86%), sensitivity (55.6%), AUROC (0.74), LHR+ (7.62) and LHR- (0.48). Concerning the classifiers based on the four shape features, we obtained the best accuracy (94%), sensitivity (66.7%), specificity (100%), AUROC (0.94), LHR+ (+∞) and LHR- (0.33) with support vector machine (SVM). This demonstrates how automatic shape features detection combined with risk classification criteria could be crucial in planning the follow-up of patients with ascending aortic aneurysm and in predicting the possible dangerous progression of the disease.

Refining Our Understanding of the Flow Through Coronary Artery Branches; Revisiting Murray's Law in Human Epicardial Coronary Arteries

Background: Quantification of coronary blood flow is used to evaluate coronary artery disease, but our understanding of flow through branched systems is poor. Murray’s law defines coronary morphometric scaling, the relationship between flow (Q) and vessel diameter (D) and is the basis for minimum lumen area targets when intervening on bifurcation lesions. Murray’s original law (Q α D^P) dictates that the exponent (P) is 3.0, whilst constant blood velocity throughout the system would suggest an exponent of 2.0. In human coronary arteries, the value of Murray’s exponent remains unknown.

Aim: To establish the exponent in Murray’s power law relationship that best reproduces coronary blood flows (Q) and microvascular resistances (Rmicro) in a bifurcating coronary tree.

Methods and Results: We screened 48 cases, and were able to evaluate inlet Q and Rmicro in 27 branched coronary arteries, taken from 20 patients, using a novel computational fluid dynamics (CFD) model which reconstructs 3D coronary anatomy from angiography and uses pressure-wire measurements to compute Q and Rmicro distribution in the main- and side-branches. Outputs were validated against invasive measurements using a Rayflow™ catheter. A Murray’s power law exponent of 2.15 produced the strongest correlation and closest agreement with inlet Q (zero bias, r = 0.47, p = 0.006) and an exponent of 2.38 produced the strongest correlation and closest agreement with Rmicro (zero bias, r = 0.66, p = 0.0001).

Conclusions: The optimal power law exponents for Q and Rmicro were not 3.0, as dictated by Murray’s Law, but 2.15 and 2.38 respectively. These data will be useful in assessing patient-specific coronary physiology and tailoring revascularisation decisions.

Development of Digital Twins to Optimize Trauma Surgery and Postoperative Management. A Case Study Focusing on Tibial Plateau Fracture

Background and context: Surgical procedures are evolving toward less invasive and more tailored approaches to consider the specific pathology, morphology, and life habits of a patient. However, these new surgical methods require thorough preoperative planning and an advanced understanding of biomechanical behaviors. In this sense, patient-specific modeling is developing in the form of digital twins to help personalized clinical decision-making.

Purpose: This study presents a patient-specific finite element model approach, focusing on tibial plateau fractures, to enhance biomechanical knowledge to optimize surgical trauma procedures and improve decision-making in postoperative management.

Study design: This is a level 5 study.

Methods: We used a postoperative 3D X-ray image of a patient who suffered from depression and separation of the lateral tibial plateau. The surgeon stabilized the fracture with polymethyl methacrylate cement injection and bi-cortical screw osteosynthesis. A digital twin of the patient’s fracture was created by segmentation. From the digital twin, four stabilization methods were modeled including two screw lengths, whether or not, to inject PMMA cement. The four stabilization methods were associated with three bone healing conditions resulting in twelve scenarios. Mechanical strength, stress distribution, interfragmentary strains, and fragment kinematics were assessed by applying the maximum load during gait. Repeated fracture risks were evaluated regarding to the volume of bone with stress above the local yield strength and regarding to the interfragmentary strains.

Results: Stress distribution analysis highlighted the mechanical contribution of cement injection and the favorable mechanical response of uni-cortical screw compared to bi-cortical screw. Evaluation of repeated fracture risks for this clinical case showed fracture instability for two of the twelve simulated scenarios.

Conclusion: This study presents a patient-specific finite element modeling workflow to assess the biomechanical behaviors associated with different stabilization methods of tibial plateau fractures. Strength and interfragmentary strains were evaluated to quantify the mechanical effects of surgical procedures. We evaluate repeated fracture risks and provide data for postoperative management.

Influence of balloon design, plaque material composition, and balloon sizing on acute post angioplasty outcomes: An implicit finite element analysis

In this work we propose a generic modeling approach for simulating percutaneous transluminal angioplasty (PTA) endovascular treatment, and evaluating the influence of balloon design, plaque composition, and balloon sizing on acute post-procedural outcomes right after PTA, without stent implantation. Clinically-used PTA balloons were classified into two categories according to their compliance characteristics, and were modeled correspondingly. Self-defined elastoplastic constitutive laws were implemented within the plaque and artery models, after calibration based on experimental and clinical data. Finite element method (FEM) implicit solver was used to simulate balloon inflation and deflation. Besides balloon profile at max inflation, results are mainly assessed in terms of the elastic recoil ratio (ERR) and lumen gain ratio (LGR) obtained immediately after PTA. No variations in ERR nor LGR values were detected when the balloon design changed, despite the differences observed in their profile at max inflation. Moreover, LGR and ERR inversely varied with the augmentation of calcification level within the plaque (-11% vs. +4% respectively, from fully lipidic to fully calcified plaque). Furthermore, results showed a direct correlation between balloon sizing and LGR and ERR, with noticeably higher rates of change for LGR (+18% and +2% for LGR and ERR respectively for a calcified plaque and a balloon pressure increasing from 10 to 14 atm). However a larger LGR comes with a higher risk of arterial rupture. This proposed methodology opens the way for evaluation of angioplasty balloon selections towards clinical procedure optimization.

Patient-Specific Finite Element Models of Posterior Pedicle Screw Fixation: Effect of Screw's Size and Geometry

Pedicle screw fixation is extensively performed to treat spine injuries or diseases and it is common for thoracolumbar fractures. Post-operative complications may arise from this surgery leading to back pain or revisions. Finite element (FE) models could be used to predict the outcomes of surgeries but should be verified when both simplified and realistic designs of screws are used. The aim of this study was to generate patient-specific Computed Tomography (CT)-based FE models of human vertebrae with two pedicle screws, verify the models, and use them to evaluate the effect of the screws' size and geometry on the mechanical properties of the screws-vertebra structure. FE models of the lumbar vertebra implanted with two pedicle screws were created from anonymized CT-scans of three patients. Compressive loads were applied to the head of the screws. The mesh size was optimized for realistic and simplified geometry of the screws with a mesh refinement study. Finally, the optimal mesh size was used to evaluate the sensitivity of the model to changes in screw's size (diameter and length) and geometry (realistic or simplified). For both simplified and realistic models, element sizes of 0.6 mm in the screw and 1.0 mm in the bone allowed to obtain relative differences of approximately 5% or lower. Changes in screw's length resulted in 4-10% differences in maximum deflection, 1-6% differences in peak stress in the screws, 10-22% differences in mean strain in the bone around the screw; changes in screw's diameter resulted in 28-36% differences in maximum deflection, 6-27% differences in peak stress in the screws, and 30-47% differences in mean strain in the bone around the screw. The maximum deflection predicted with realistic or simplified screws correlated very well (R ² = 0.99). The peak stress in screws with realistic or simplified design correlated well (R ² = 0.82) but simplified models underestimated the peak stress. In conclusion, the results showed that the diameter of the screw has a major role on the mechanics of the screw-vertebral structure for each patient. Simplified screws can be used to estimate the mechanical properties of the implanted vertebrae, but the systematic underestimation of the peak stress should be considered when interpreting the results from the FE analyses.

Numerical simulation of fenestrated graft deployment: Anticipation of stent graft and vascular structure adequacy

Fenestrated endovascular aneurism repair (FEVAR) is a minimally invasive technique, and its success depends on the adequacy of the correspondence between the visceral arteries ostia and position of the fenestrations of the stent graft (SG) during its deployment in juxtarenal aneurisms. However, the fenestration position is generally determined from a preoperative computerised tomography (CT) scan, without considering the vascular deformation induced by the insertion of the endovascular tools. Catheterisation difficulties may occur during clinical procedures. Accordingly, the objective of this work is to present an initial proof of concept aimed at anticipating and optimising the position of the fenestrations, while considering the vascular deformation induced by the insertion of the endovascular tools. The proposed method relies on the finite element method to simulate the SG deployment in a vascular structure (VS), and considers the vascular deformation induced by the tools. After determining the optimal simulation parameters for a patient-specific case, the robustness of the method is demonstrated on six other representative anatomies. The simulated SG is also compared with post-deployment CT observations, and demonstrates good adequacy. The results show that the numerically corrected fenestration positions, as determined from the simulated results following the insertion of the endovascular tools, deviate from those of the standard plan (as determined from the preoperative CT scan). This indicates that the SG-VS adequacy could be improved via simulation-based planning, to anticipate potential catheterisation difficulties.

Machine-Learning based model order reduction of a biomechanical model of the human tongue

BACKGROUND AND OBJECTIVES

This paper presents the results of a Machine-Learning based Model Order Reduction (MOR) method applied to a complex 3D Finite Element (FE) biomechanical model of the human tongue, in order to create a Digital Twin Model (DTM) that enables real-time simulations. The DTM is designed for future inclusion in a computer assisted protocol for tongue surgery planning.

METHODS

The proposed method uses an "a posteriori" MOR that allows, from a limited number of simulations with the FE model, to predict in real time mechanical responses of the human tongue to muscle activations.

RESULTS

The MOR method is evaluated for simulations associated with separate single tongue muscle activations. It is shown to be able to account with a sub-millimetric spatial accuracy for the non-linear dynamical behavior of the tongue model observed in these simulations.

CONCLUSION

Further evaluations of the MOR method will include tongue movements induced by multiple muscle activations. At this stage our MOR method offers promising perspectives for the use of the tongue model in a clinical context to predict the impact of tongue surgery on tongue mobility. As a long term application, this DTM of the tongue could be used to predict the functional consequences of the surgery in terms of speech production and swallowing.

The influence of angioplasty balloon sizing on acute post-procedural outcomes: a Finite Element Analysis

Atherosclerosis is one of the most common vascular pathologies in the world. Among the most commonly performed endovascular treatments, percutaneous transluminal angioplasty (PTA) has been showing significantly positive clinical outcomes. Due to the complex geometries, material properties and interactions that characterize PTA procedures, finite element analyses of acute angioplasty balloon deployment are limited. In this work, finite element method (FEM) was used to simulate the inflation and deflation of a semi-compliant balloon within the 3D model of a stenosed artery with two different plaque types (lipid and calcified). Self-defined constitutive models for the balloon and the plaque were developed based on experimental and literature data respectively. Balloon deployment was simulated at three different inflation pressures (10, 12 and 14 atm) within the two plaque types. Balloon sizing influence on the arterial elastic recoil obtained immediately after PTA was then investigated. The simulated results show that calcified plaques may lead to higher elastic recoil ratios compared to lipid stenosis, when the same balloon inflation pressures are applied. Also, elastic recoil increases for higher balloon inflation pressure independent of the plaque type. These findings open the way for a data-driven assessment of angioplasty balloon sizing selection and clinical procedures optimization.Clinical Relevance- The FE model developed in this work aims at providing quantitative evaluation of recoil after balloon angioplasty. It may be useful for both manufacturers and clinicians to improve efficiency of angioplasty balloon device design and sizing selection with respect to plaque geometry and constitution, consequently enhancing clinical outcomes.

PRIMAGE project: predictive in silico multiscale analytics to support childhood cancer personalised evaluation empowered by imaging biomarkers

PRIMAGE is one of the largest and more ambitious research projects dealing with medical imaging, artificial intelligence and cancer treatment in children. It is a 4-year European Commission-financed project that has 16 European partners in the consortium, including the European Society for Paediatric Oncology, two imaging biobanks, and three prominent European paediatric oncology units. The project is constructed as an observational in silico study involving high-quality anonymised datasets (imaging, clinical, molecular, and genetics) for the training and validation of machine learning and multiscale algorithms. The open cloud-based platform will offer precise clinical assistance for phenotyping (diagnosis), treatment allocation (prediction), and patient endpoints (prognosis), based on the use of imaging biomarkers, tumour growth simulation, advanced visualisation of confidence scores, and machine-learning approaches. The decision support prototype will be constructed and validated on two paediatric cancers: neuroblastoma and diffuse intrinsic pontine glioma. External validation will be performed on data recruited from independent collaborative centres. Final results will be available for the scientific community at the end of the project, and ready for translation to other malignant solid tumours.

Patient-specific simulation of guidewire deformation during transcatheter aortic valve implantation

Transcatheter aortic valve implantation is a recent mini-invasive procedure to implant an aortic valve prosthesis. Prosthesis positioning in transcatheter aortic valve implantation appears as an important aspect for the success of the intervention. Accordingly, we developed a patient-specific finite element framework to predict the insertion of the stiff guidewire, used to position the aortic valve. We simulated the guidewire insertion for 2 patients based on their pre-operative CT scans. The model was designed to primarily predict the position and the angle of the guidewires in the aortic valve, and the results were successfully compared with intraoperative images. The present paper describes extensively the numerical model, which was solved by using the ANSYS software with an implicit resolution scheme, as well as the stabilization techniques which were used to overcome numerical instabilities. We performed sensitivity analysis on the properties of the guidewire (curvature angle, curvature radius, and stiffness) and the conditions of insertion (insertion force and orientation). We also explored the influence of the model parameters. The accuracy of the model was quantitatively evaluated as the distance and the angle difference between the simulated guidewires and the intraoperative ones. A good agreement was obtained between the model predictions and intraoperative views available for 2 patient cases. In conclusion, we showed that the shape of the guidewire in the aortic valve was mainly determined by the geometry of the patient's aorta and by the conditions of insertion (insertion force and orientation).

Simulation of nonlinear transient elastography: finite element model for the propagation of shear waves in homogeneous soft tissues

In this study, visco-hyperelastic Landau's model, which is widely used in acoustical physic field, is introduced into a finite element formulation. It is designed to model the nonlinear behaviour of finite amplitude shear waves in soft solids, typically, in biological tissues. This law is used in finite element models based on elastography, experiments reported in Jacob et al, the simulations results show a good agreement with the experimental study: It is observed in both that a plane shear wave generates only odd harmonics and a nonplane wave generates both odd and even harmonics in the spectral domain. In the second part, a parametric study is performed to analyse the influence of different factors on the generation of odd harmonics of plane wave. A quantitative relation is fitted between the odd harmonic amplitudes and the non-linear elastic parameter of Landau's model, which provides a practical guideline to identify the non-linearity of homogeneous tissues using elastography experiment.

Personalized modeling for real-time pressure ulcer prevention in sitting posture

Ischial pressure ulcer is an important risk for every paraplegic person and a major public health issue. Pressure ulcers appear following excessive compression of buttock's soft tissues by bony structures, and particularly in ischial and sacral bones. Current prevention techniques are mainly based on daily skin inspection to spot red patches or injuries. Nevertheless, most pressure ulcers occur internally and are difficult to detect early. Estimating internal strains within soft tissues could help to evaluate the risk of pressure ulcer. A subject-specific biomechanical model could be used to assess internal strains from measured skin surface pressures. However, a realistic 3D non-linear Finite Element buttock model, with different layers of tissue materials for skin, fat and muscles, requires somewhere between minutes and hours to compute, therefore forbidding its use in a real-time daily prevention context. In this article, we propose to optimize these computations by using a reduced order modeling technique (ROM) based on proper orthogonal decompositions of the pressure and strain fields coupled with a machine learning method. ROM allows strains to be evaluated inside the model interactively (i.e. in less than a second) for any pressure field measured below the buttocks. In our case, with only 19 modes of variation of pressure patterns, an error divergence of one percent is observed compared to the full scale simulation for evaluating the strain field. This reduced model could therefore be the first step towards interactive pressure ulcer prevention in a daily set-up.

Finite element simulation of the insertion of guidewires during an EVAR procedure: example of a complex patient case, a first step toward patient-specific parameterized models

Deformations of the vascular structure due to the insertion of tools during endovascular treatment of aneurysms of the abdominal aorta, unless properly anticipated during the preoperative planning phase, may be the source of intraoperative or postoperative complications. We propose here an explicit finite element simulation method which enables one to predict such deformations. This method is based on a mechanical model of the vascular structure which takes into account the nonlinear behavior of the arterial wall, the prestressing effect induced by the blood pressure and the mechanical support of the surrounding organs and structures. An analysis of the model sensitivity to the parameters used to represent this environment is done. This allows determining the parameters that have the largest influence on the quality of the prediction and also provides realistic values for each of them as no experimental data are available in the literature. Moreover, for the first time, the results are compared with 3D intraoperative data. This is done for a patient-specific case with a complex anatomy in order to assess the feasibility of the method. Finally, the predictive capability of the simulation is evaluated on a group of nine patients. The error between the final simulated and intraoperatively measured tool positions was 2.1 mm after the calibration phase on one patient. It results in a 4.6 ± 2.5 mm in average error for the blind evaluation on nine patients.

Enhancement of Peripheral Stents Reliability: Developing Interactive Procedure Planning by Means of Numerical Simulations and Clinical Software Development

Around 20% of the population over 60-years-old have peripheral arterial disease. Symptoms can cause major lifestyle limitation through pain on walking (intermittent claudication). Progressive disease can also lead to critical limb ischemia, which is the major cause of amputation for adults. The implant of a stent is the endovascular therapy of election for the treatment of peripheral arterial disease. This, however, can increase the risk of fatigue fracture of the implanted device, under in vivo loading conditions. RT3S (real time simulations for safer vascular stenting) is a scientific project involving both universities and private companies, for a funded budget of more than $4.5 million (www.rt3s.eu). The project aims at providing physicians with interactive clinical software for an enhanced pre-operative planning of femoral artery stenting including a prediction of in vivo stent fracture risk, computed by means of numerical simulations. Software-developing companies developed an interactive application allowing physicians to upload the patient's clinical images and plan the stenting procedure, estimating the device fatigue fracture risk, for a specific case of stent and anatomical geometry/boundary conditions. Thanks to the application of innovative response surface techniques, results obtained with the numerical model (requiring 8-day running time) can be exploited quasi real time by the software user. Outcomes of the project will be also exploited for designing activities of next-generation stent devices.

Finite-element-based matching of pre- and intraoperative data for image-guided endovascular aneurysm repair

Endovascular repair of abdominal aortic aneurysms is a well-established technique throughout the medical and surgical communities. Although increasingly indicated, this technique does have some limitations. Because intervention is commonly performed under fluoroscopic control, 2-D visualization of the aneurysm requires the injection of a contrast agent. The projective nature of this imaging modality inevitably leads to topographic errors, and does not give information on arterial wall quality at the time of deployment. A specially adapted intraoperative navigation interface could increase deployment accuracy and reveal such information, which preoperative 3-D imaging might otherwise provide. One difficulty is the precise matching of preoperative data (images and models) and intraoperative observations affected by anatomical deformations due to tool-tissue interactions. Our proposed solution involves a finite-element-based preoperative simulation of tool-tissue interactions, its adaptive tuning regarding patient specific data, and the matching with intraoperative data. The biomechanical model was first tuned on a group of ten patients and assessed on a second group of eight patients.

Evaluation of the generality and accuracy of a new mesh morphing procedure for the human femur

Various papers described mesh morphing techniques for computational biomechanics, but none of them provided a quantitative assessment of generality, robustness, automation, and accuracy in predicting strains. This study aims to quantitatively evaluate the performance of a novel mesh-morphing algorithm. A mesh-morphing algorithm based on radial-basis functions and on manual selection of corresponding landmarks on template and target was developed. The periosteal geometries of 100 femurs were derived from a computed tomography scan database and used to test the algorithm generality in producing finite element (FE) morphed meshes. A published benchmark, consisting of eight femurs for which in vitro strain measurements and standard FE model strain prediction accuracy were available, was used to assess the accuracy of morphed FE models in predicting strains. Relevant parameters were identified to test the algorithm robustness to operative conditions. Time and effort needed were evaluated to define the algorithm degree of automation. Morphing was successful for 95% of the specimens, with mesh quality indicators comparable to those of standard FE meshes. Accuracy of the morphed meshes in predicting strains was good (R(2)>0.9, RMSE%<10%) and not statistically different from the standard meshes (p-value=0.1083). The algorithm was robust to inter- and intra-operator variability, target geometry refinement (p-value>0.05) and partially to the number of landmark used. Producing a morphed mesh starting from the triangularized geometry of the specimen requires on average 10 min. The proposed method is general, robust, automated, and accurate enough to be used in bone FE modelling from diagnostic data, and prospectively in applications such as statistical shape modelling.

Physiological changes in maternal cortisol do not alter expression of growth-related genes in the ovine placenta

OBJECTIVES

The aim of this study was to investigate the effect of cortisol on growth-related genes in the ovine placenta.

STUDY DESIGN

Ewes carrying singleton pregnancies were operated on between 112 and 116 days of gestation (115 ± 0.4, term = 147 days) and randomly assigned into three groups: six control animals, five ewes that were administered cortisol by continuous intravenous infusion (1 mg/kg/day, high cortisol), and five ewes that were adrenalectomized and replaced with 0.5-0.6 mg cortisol/kg/day and 3 μg aldosterone/kg/day to produce cortisol concentrations equivalent to pre-pregnancy values (low cortisol). At necropsy (130 ± 0.2 days of gestation), placental tissue was frozen and stored at -80 °C for mRNA analysis.

MAIN OUTCOME MEASURES

To assess potential molecular mechanisms by which cortisol alters placental structure and function and fetal growth.

RESULTS

Cortisol levels did not significantly affect 11β-hydroxysteroid dehydrogenase 1 and 2 enzymes, glucocorticoid receptor, mineralocorticoid receptor and angiotensin II receptor, type 1 (AT1R) expression levels. Gene expression levels of AT2R were increased in the high cortisol group for type B placentomes. There was little effect of cortisol on the insulin-like growth factor (IGF) axis. There was significantly more IGF-I mRNA in B versus A type and more IGFBP-2 mRNA in B and C type versus A type placentomes regardless of treatment (p < 0.05).

CONCLUSIONS

These data suggest that cortisol increases placental AT2R expression at high concentrations whereas it has little effect on the placental IGF axis.

Statistical Shape Modeling of Femurs Using Morphing and Principal Component Analysis

In this paper, we describe a method for automatically building a statistical shape model by applying a morphing method and a principal component analysis (PCA) to a large database of femurs. One of the major challenges in building a shape model from a training data set of 3D objects is the determination of the correspondence between different shapes. In our work, we solve this problem by using a morphing method. The morphing method consists of deforming the same template mesh over a large database of femur geometries, which results in isotopological meshes and one to one correspondences; i.e., the resulting meshes have the same number of nodes, the same number of elements, and the same connectivity in all morphed meshes. By applying the morphing-based registration followed by PCA to a large database of femurs, we demonstrate that the method can be used to derive a low dimensional representation of the main variabilities of the femur geometry.

Output characterization of a self-pulsating and aperiodic optical fiber source based on cascaded regeneration

We investigate the properties of a self-pulsating fiber cavity based on cascaded regeneration. The mechanisms that govern the number of oscillating pulses in the cavity, the pulse peak power, the pulse width, the wavelength tunability as well as the generation of sub-picosecond pulses are identified, analyzed and quantified. We find that the described self-pulsating cavity enables the oscillation of quasi transform-limited pulses with a pulsewidth of 4.8 ps at 1540.0 nm when using 0.4 nm non-Gaussian bandpass filters. Sub-picosecond pulses with an autocorrelation width of 471 fs are generated from the same self-pulsating source with modified bandpass filters and the addition of a chromatic dispersion compensator. The number of eigenpulses that oscillate simultaneously in the cavity can be adjusted from 0 up to 29,500 with proper cavity adjustment. This source has dual-wavelength output and can be tuned throughout the gain band of the amplifiers.

Fluid/structure interaction applied to the simulation of Abdominal Aortic Aneurysms

Aneurysms are a local dilatation of a vessel wall of at least twice the normal diameter (commonly accepted definition). They are asymptomatic and rupture is often lethal. Thus, prediction of rupture is an important stake. Aiming at a diagnosis tool relying on patient specific data and general physiological values, we created a virtual aneurysm model based on real imaging data. Fluid/structure interaction (FSI) simulations were made to compute the displacement and stress for the wall. For the fluid, the only in vivo measures used were for the inlet velocity. The mandatory output boundary condition has been implemented with the first order Windkessel model equations. Structure has been much more complicated to handle because of the association of a realistic geometry (no symmetry) and a full fluid/structure interaction approach. We used surface elements to stabilize the structure and to model surrounding organs. Validation parameters are the displacement, the Von Mises stress and the pressure profile at the outlet. The main difference with other studies relies on the association of all these elements in order to prepare industrial applications as the main goal of this study was to build an automated tool easy to use by people who are not experts in numerical simulation.

Investigation of self-phase modulation based optical regeneration in single mode As2Se3 chalcogenide glass fiber

We investigate the feasibility of all-optical regeneration based on self-phase modulation in single mode As2Se3 chalcogenide fiber. By combining the chalcogenide fiber with a bandpass filter, we achieve a near step-like power transfer function with no pulse distortion. The device is shown to operate with 5.8 ps duration pulses, thus demonstrating the feasibility of this device operating with high bit-rate data signals. These results are achieved with pulse peak powers <10 W in a fully passive device, including only 2.8 m of chalcogenide fiber. We obtain an excellent agreement between theory and experiment and show that both the high nonlinearity of the chalcogenide glass along with its high normal dispersion near 1550 nm enables a significant device length reduction in comparison with silica-based devices, without compromise on the performance. We find that even for only a few meters of fiber, the large normal dispersion of the chalcogenide glass inhibits spectral oscillations that would appear with self-phase modulation alone. We measure the two photon absorption attenuation coefficient and find that it advantageously affects the device transfer function.

In-band OSNR and chromatic dispersion monitoring using a fibre optical parametric amplifier

This paper presents an all-optical, in-band optical signal-to-noise ratio (OSNR) and chromatic dispersion monitor. We demonstrate monitoring over the 1 nm bandwidth of our signal, which is a 10 GHz pulse train of 8.8 ps pulses. The monitor output power, as measured on a slow detector, has a 1.9 dB dynamic range when the signal OSNR is varied by 20 dB, and a 1.6 dB dynamic range when +/- 150 ps/nm of chromatic dispersion is applied. Cascaded four-wave mixing occurring in the optical parametric amplifier provides the nonlinear power transfer function responsible for the monitoring. An analysis using the signals' probability density functions show that the nonlinear power transfer function provides preferential gain to clean undispersed pulses when compared to noisy and/or dispersed pulses. Our analysis includes a consideration of the applicability of the device to high duty cycle systems, and simulations on monitoring of a 40 Gb/s pulse train with a 50% duty cycle.

Osmoregulation of vasopressin secretion via activation of neurohypophysial nerve terminals glycine receptors by glial taurine

Osmotic regulation of supraoptic nucleus (SON) neuron activity depends in part on activation of neuronal glycine receptors (GlyRs), most probably by taurine released from adjacent astrocytes. In the neurohypophysis in which the axons of SON neurons terminate, taurine is also concentrated in and osmo-dependently released by pituicytes, the specialized glial cells ensheathing nerve terminals. We now show that taurine release from isolated neurohypophyses is enhanced by hypo-osmotic and decreased by hyper-osmotic stimulation. The high osmosensitivity is shown by the significant increase on only 3.3% reduction in osmolarity. Inhibition of taurine release by 5-nitro-2-(3-phenylpropylamino)benzoic acid, niflumic acid, and 4,4'-diisothiocyanatostilbene-2,2'-disulfonic acid suggests the involvement of volume-sensitive anion channels. On purified neurohypophysial nerve endings, activation of strychnine-sensitive GlyRs by taurine or glycine primarily inhibits the high K(+)-induced rise in [Ca(2+)](i) and subsequent release of vasopressin. Expression of GlyRs in vasopressin and oxytocin terminals is confirmed by immunohistochemistry. Their implication in the osmoregulation of neurohormone secretion was assessed on isolated whole neurohypophyses. A 6.6% hypo-osmotic stimulus reduces by half the depolarization-evoked vasopressin secretion, an inhibition totally prevented by strychnine. Most importantly, depletion of taurine by a taurine transport inhibitor also abolishes the osmo-dependent inhibition of vasopressin release. Therefore, in the neurohypophysis, an osmoregulatory system involving pituicytes, taurine, and GlyRs is operating to control Ca(2+) influx in and neurohormone release from nerve terminals. This elucidates the functional role of glial taurine in the neurohypophysis, reveals the expression of GlyRs on axon terminals, and further defines the role of glial cells in the regulation of neuroendocrine function.

Importance of attributable risk in monitoring adverse events after immunization: hepatitis B vaccination in children

OBJECTIVES

Most vaccine safety data present only the postvaccination incidence of all adverse events rather than an estimate of attributable risk. This study sought to illustrate the difference between the 2 estimates with data from a hepatitis B immunization program.

METHODS

The incidence of health problems occurring before and after each dose of hepatitis B vaccine in a cohort of 1130 children were compared.

RESULTS

Although 47.5% of all children reported an adverse event during the 4 weeks following each of the 3 doses, adverse events attributable to immunization occurred in only 10.6% of children.

CONCLUSIONS

Postimmunization incidence systematically overestimates the risk of adverse events. Estimating actual attributable risk is necessary to avoid false beliefs regarding immunization.

A single dose of an effective whole cell pertussis vaccine does not significantly increase protection in children primed with a less effective vaccine

We evaluated if a single dose of a protective whole cell pertussis vaccine given before school entry to children primed with a less effective vaccine would increase their protection. A school cohort including 3876 students and a family cohort including 162 children were assessed. Although there was a trend toward increased protection. the better vaccine did not provide a significant improvement. These results suggest that a single dose of an effective vaccine given to children primed with a less effective one does not raise the protection to at level similar to that provided by three doses of the better vaccine.