The influence of systemic-to-pulmonary arterial shunts and peripheral vasculatures in univentricular circulations: Focus on coronary perfusion and aortic arch hemodynamics through computational multi-domain modeling

Initial palliation for univentricular hearts can be achieved via a systemic-to-pulmonary shunt (SPS). SPS configurations differ depending on the proximal anastomosis location, which might lead to dissimilar coronary and upper body perfusions. Mathematical modeling can be used to explore the local and global hemodynamic effects of the SPSs. In literature there are few patient-specific models of SPS that specifically address the influence of both the local and peripheral vasculature. In this study, multi-domain models of univentricular circulations were developed to investigate local hemodynamics and flow distribution in the presence of two shunt configurations. We also analyzed the relative impact of local and peripheral vascular resistances on coronary perfusion and flows through the upper aortic branches.A two-step approach was followed. First, two patient-specific models were based on clinical data collected from univentricular patients having different shunts and peripheral vasculatures. Each model coupled a three-dimensional representation of SPS, aortic arch (AA) and pulmonary arteries, with a lumped parameter model (LPM) of peripheral vasculature closing the circulatory loop. Then, two additional models of hypothetical subjects were created by coupling each customized LPM with the other patient’s three-dimensional anatomy.Flow rates and pressures predicted by the patient-specific models revealed overall agreement with clinical data. Differences in the local hemodynamics were seen during diastole between the two models. Varying the three-dimensional models, while keeping an identical LPM, led to comparable flow distribution through the AA, suggesting that peripheral vasculatures have a dominant effect on local hemodynamics with respect to the shunt configuration.     

3D Rapid Prototyping for Otolaryngology—Head and Neck Surgery: Applications in Image-Guidance,Surgical Simulation and Patient-Specific Modeling

The aim of this study was to demonstrate the role of advanced fabrication technology across a broad spectrum of head and neck surgical procedures, including applications in endoscopic sinus surgery, skull base surgery, and maxillofacial reconstruction. The initial case studies demonstrated three applications of rapid prototyping technology are in head and neck surgery: i) a mono-material paranasal sinus phantom for endoscopy training ii) a multi-material skull base simulator and iii) 3D patient-specific mandible templates. Digital processing of these phantoms is based on real patient or cadaveric 3D images such as CT or MRI data. Three endoscopic sinus surgeons examined the realism of the endoscopist training phantom. One experienced endoscopic skull base surgeon conducted advanced sinus procedures on the high-fidelity multi-material skull base simulator. Ten patients participated in a prospective clinical study examining patient-specific modeling for mandibular reconstructive surgery. Qualitative feedback to assess the realism of the endoscopy training phantom and high-fidelity multi-material phantom was acquired. Conformance comparisons using assessments from the blinded reconstructive surgeons measured the geometric performance between intra-operative and pre-operative reconstruction mandible plates. Both the endoscopy training phantom and the high-fidelity multi-material phantom received positive feedback on the realistic structure of the phantom models. Results suggested further improvement on the soft tissue structure of the phantom models is necessary. In the patient-specific mandible template study, the pre-operative plates were judged by two blinded surgeons as providing optimal conformance in 7 out of 10 cases. No statistical differences were found in plate fabrication time and conformance, with pre-operative plating providing the advantage of reducing time spent in the operation room. The applicability of common model design and fabrication techniques across a variety of otolaryngological sub-specialties suggests an emerging role for rapid prototyping technology in surgical education, procedure simulation, and clinical practice.     

Numerical modeling of hemodynamics scenarios of patient-specific coronary artery bypass grafts

A fast computational framework is devised to the study of several configurations of patient-specific coronary artery bypass grafts. This is especially useful to perform a sensitivity analysis of the hemodynamics for different flow conditions occurring in native coronary arteries and bypass grafts, the investigation of the progression of the coronary artery disease and the choice of the most appropriate surgical procedure. A complete pipeline, from the acquisition of patient-specific medical images to fast parameterized computational simulations, is proposed. Complex surgical configurations employed in the clinical practice, such as Y-grafts and sequential grafts, are studied. A virtual surgery platform based on model reduction of unsteady Navier–Stokes equations for blood dynamics is proposed to carry out sensitivity analyses in a very rapid and reliable way. A specialized geometrical parameterization is employed to compare the effect of stenosis and anastomosis variation on the outcome of the surgery in several relevant cases.     

Personalized Computer Simulation of Diastolic Function in Heart Failure

The search for a parameter representing left ventricular relaxation from non-invasive and invasive diagnostic tools has been extensive, since heart failure (HF) with preserved ejection fraction (HF-pEF) is a global health problem. We explore here the feasibility using patient-specific cardiac computer modeling to capture diastolic parameters in patients suffering from different degrees of systolic HF. Fifty eight patients with idiopathic dilated cardiomyopathy have undergone thorough clinical evaluation, including cardiac magnetic resonance imaging (MRI), heart catheterization, echocardiography, and cardiac biomarker assessment. A previously-introduced framework for creating multi-scale patient-specific cardiac models has been applied on all these patients. Novel parameters, such as global stiffness factor and maximum left ventricular active stress, representing cardiac active and passive tissue properties have been computed for all patients. Invasive pressure measurements from heart catheterization were then used to evaluate ventricular relaxation using the time constant of isovolumic relaxation Tau (s). Parameters from heart catheterization and the multi-scale model have been evaluated and compared to patient clinical presentation. The model parameter global stiffness factor, representing diastolic passive tissue properties, is correlated signif-icantly across the patient population with s. This study shows that multi-modal cardiac models can successfully capture diastolic (dys) function, a prerequisite for future clinical trials on HF-pEF.     

Patient-specific modeling of dyssynchronous heart failure: A case study

The development and clinical use of patient-specific models of the heart is now a feasible goal. Models have the potential to aid in diagnosis and support decision-making in clinical cardiology. Several groups are now working on developing multi-scale models of the heart for understanding therapeutic mechanisms and better predicting clinical outcomes of interventions such as cardiac resynchronization therapy. Here we describe the methodology for generating a patient-specific model of the failing heart with a myocardial infarct and left ventricular bundle branch block. We discuss some of the remaining challenges in developing reliable patient-specific models of cardiac electromechanical activity, and identify some of the main areas for focusing future research efforts. Key challenges include: efficiently generating accurate patient-specific geometric meshes and mapping regional myofiber architecture to them; modeling electrical activation patterns based on cellular alterations in human heart failure, and estimating regional tissue conductivities based on clinically available electrocardiographic recordings; estimating unloaded ventricular reference geometry and material properties for biomechanical simulations; and parameterizing systemic models of circulatory dynamics from available hemodynamic measurements.     

Orbital and maxillofacial computer aided surgery: patient-specific finite element models to predict surgical outcomes

This paper addresses an important issue raised for the clinical relevance of Computer-Assisted Surgical applications, namely the methodology used to automatically build patient-specific finite element (FE) models of anatomical structures. From this perspective, a method is proposed, based on a technique called the mesh-matching method, followed by a process that corrects mesh irregularities. The mesh-matching algorithm generates patient-specific volume meshes from an existing generic model. The mesh regularization process is based on the Jacobian matrix transform related to the FE reference element and the current element.This method for generating patient-specific FE models is first applied to computer-assisted maxillofacial surgery, and more precisely, to the FE elastic modelling of patient facial soft tissues. For each patient, the planned bone osteotomies (mandible, maxilla, chin) are used as boundary conditions to deform the FE face model, in order to predict the aesthetic outcome of the surgery. Seven FE patient-specific models were successfully generated by our method. For one patient, the prediction of the FE model is qualitatively compared with the patient's post-operative appearance, measured from a computer tomography scan. Then, our methodology is applied to computer-assisted orbital surgery. It is, therefore, evaluated for the generation of 11 patient-specific FE poroelastic models of the orbital soft tissues. These models are used to predict the consequences of the surgical decompression of the orbit. More precisely, an average law is extrapolated from the simulations carried out for each patient model. This law links the size of the osteotomy (i.e. the surgical gesture) and the backward displacement of the eyeball (the consequence of the surgical gesture).     

Towards a Complete In Silico Assessment of the Outcome of Cochlear Implantation Surgery

Cochlear implantation (CI) surgery is a very successful technique, performed on more than 300,000 people worldwide. However, since the challenge resides in obtaining an accurate surgical planning, computational models are considered to provide such accurate tools. They allow us to plan and simulate beforehand surgical procedures in order to maximally optimize surgery outcomes, and consequently provide valuable information to guide pre-operative decisions. The aim of this work is to develop and validate computational tools to completely assess the patient-specific functional outcome of the CI surgery. A complete automatic framework was developed to create and assess computationally CI models, focusing on the neural response of the auditory nerve fibers (ANF) induced by the electrical stimulation of the implant. The framework was applied to evaluate the effects of ANF degeneration and electrode intra-cochlear position on nerve activation. Results indicate that the intra-cochlear positioning of the electrode has a strong effect on the global performance of the CI. Lateral insertion provides better neural responses in case of peripheral process degeneration, and it is recommended, together with optimized intensity levels, in order to preserve the internal structures. Overall, the developed automatic framework provides an insight into the global performance of the implant in a patient-specific way. This enables to further optimize the functional performance and helps to select the best CI configuration and treatment strategy for a given patient.     

Finite element simulation of spinal deformities correction by in situ contouring technique

Biomechanical models have been proposed in order to simulate the surgical correction of spinal deformities. With these models, different surgical correction techniques have been examined: distraction and rod rotation. The purpose of this study was to simulate another surgical correction technique: the in situ contouring technique. In this way, a comprehensive three-dimensional Finite Element (FE) model with patient-specific geometry and patient-specific mechanical properties was used. The simulation of the surgery took into account elasto-plastic behavior of the rod and multiple moments loading and unloading representing the surgical maneuvers. The simulations of two clinical cases of hyperkyphosis and scoliosis were coherent with the surgeon's experience. Moreover, the results of simulation were compared to post-operative 3D measurements. The mean differences were under 5 degrees for vertebral rotations and 5 mm for spinal lines. These simulations open the way for future predictive tools for surgical planning.     

Improving surgical outcome following the Norwood procedure

Background

The Norwood procedure consists of three palliative operations, performed in neonates with hypoplastic left heart syndrome. Especially the first stage (Norwood I) is associated with the highest mortality rates in paediatric cardiac surgery (up to 25%). During surgery, the aorta is reconstructed and a systemic-to-pulmonary shunt is applied. Originally the modified Blalock-Taussig shunt was used, but recently the right-ventricle-to-pulmonary-artery shunt is increasingly being employed. We reviewed the results of our operative strategy, where an individualised choice of shunt is made. Furthermore, attempts to reduce interstage mortality (between Norwood I and II) were assessed.

Methods

All neonates who underwent Norwood stage I palliation from August 2004 until November 2010 were included in this retrospective analysis. Mortality rates and management strategies were compared.

Results

Thirty-six patients were available for analysis. Overall 30-day mortality was 5.6% (2 patients) and interstage mortality after discharge was 14% (5 patients). In 2006, a novel clinical protocol was introduced, aimed at reduction of mortality during the interstage period. This resulted in reduction of interstage mortality from 23% to 9% (3 of 13 infants, versus 2 of 23), with a cumulative survival of 82% (maximum follow-up 4 years).

Conclusion

Early surgical results following the Norwood procedure using an individualised shunt choice are favourable.     

Predicting the subject-specific primary stability of cementless implants during pre-operative planning: preliminary validation of subject-specific finite-element models

Pre-operative planning help the surgeon in taking the proper clinical decision. The ultimate goal of this work is to develop numerical models that allow the surgeon to estimate the primary stability during the pre-operative planning session. The present study was aimed to validate finite-element (FE) models accounting for patient and prosthetic size and position as planned by the surgeon. For this purpose, the FE model of a cadaveric femur was generated starting from the CT scan and the anatomical position of a cementless stem derived by a skilled surgeon using a pre-operative CT-based planning simulation software. In-vitro experimental measurements were used as benchmark problem to validate the bone-implant relative micromotions predicted by the patient-specific FE model. A maximum torque in internal rotation of 11.4 Nm was applied to the proximal part of the hip stem. The error on the maximum predicted micromotion was 12% of the peak micromotion measured experimentally. The average error over the entire range of applied torques was only 7% of peak measurement. Hence, the present study confirms that it is possible to accurately predict the level of primary stability achieved for cementless stems using numerical models that account for patient specificity and surgical variability.     

Finite element simulation of spinal deformities correction by in situ contouring technique

Biomechanical models have been proposed in order to simulate the surgical correction of spinal deformities. With these models, different surgical correction techniques have been examined: distraction and rod rotation. The purpose of this study was to simulate another surgical correction technique: the in situ contouring technique. In this way, a comprehensive three-dimensional Finite Element (FE) model with patient-specific geometry and patient-specific mechanical properties was used. The simulation of the surgery took into account elasto–plastic behavior of the rod and multiple moments loading and unloading representing the surgical maneuvers. The simulations of two clinical cases of hyperkyphosis and scoliosis were coherent with the surgeon's experience. Moreover, the results of simulation were compared to post-operative 3D measurements. The mean differences were under 5° for vertebral rotations and 5 mm for spinal lines. These simulations open the way for future predictive tools for surgical planning.     

Can time-averaged flow boundary conditions be used to meet the clinical timeline for Fontan surgical planning?

Cardiovascular simulations have great potential as a clinical tool for planning and evaluating patient-specific treatment strategies for those suffering from congenital heart diseases, specifically Fontan patients. However, several bottlenecks have delayed wider deployment of the simulations for clinical use; the main obstacle is simulation cost. Currently, time-averaged clinical flow measurements are utilized as numerical boundary conditions (BCs) in order to reduce the computational power and time needed to offer surgical planning within a clinical time frame. Nevertheless, pulsatile blood flow is observed in vivo, and its significant impact on numerical simulations has been demonstrated. Therefore, it is imperative to carry out a comprehensive study analyzing the sensitivity of using time-averaged BCs. In this study, sensitivity is evaluated based on the discrepancies between hemodynamic metrics calculated using time-averaged and pulsatile BCs; smaller discrepancies indicate less sensitivity.The current study incorporates a comparison between 3D patient-specific CFD simulations using both the time-averaged and pulsatile BCs for 101 Fontan patients. The sensitivity analysis involves two clinically important hemodynamic metrics: hepatic flow distribution (HFD) and indexed power loss (iPL). Paired demographic group comparisons revealed that HFD sensitivity is significantly different between single and bilateral superior vena cava cohorts but no other demographic discrepancies were observed for HFD or iPL. Multivariate regression analyses show that the best predictors for sensitivity involve flow pulsatilities, time-averaged flow rates, and geometric characteristics of the Fontan connection. These predictors provide patient-specific guidelines to determine the effectiveness of analyzing patient-specific surgical options with time-averaged BCs within a clinical time frame.     

Uncertainty quantification in cerebral circulation simulations focusing on the collateral flow: Surrogate model approach with machine learning

Collateral circulation in the circle of Willis (CoW), closely associated with disease mechanisms and treatment outcomes, can be effectively investigated using one-dimensional–zero-dimensional hemodynamic simulations. As the entire cardiovascular system is considered in the simulation, it captures the systemic effects of local arterial changes, thus reproducing collateral circulation that reflects biological phenomena. The simulation facilitates rapid assessment of clinically relevant hemodynamic quantities under patient-specific conditions by incorporating clinical data. During patient-specific simulations, the impact of clinical data uncertainty on the simulated quantities should be quantified to obtain reliable results. However, as uncertainty quantification (UQ) is time-consuming and computationally expensive, its implementation in time-sensitive clinical applications is considered impractical. Therefore, we constructed a surrogate model based on machine learning using simulation data. The model accurately predicts the flow rate and pressure in the CoW in a few milliseconds. This reduced computation time enables the UQ execution with 100 000 predictions in a few minutes on a single CPU core and in less than a minute on a GPU. We performed UQ to predict the risk of cerebral hyperperfusion (CH), a life-threatening condition that can occur after carotid artery stenosis surgery if collateral circulation fails to function appropriately. We predicted the statistics of the postoperative flow rate increase in the CoW, which is a measure of CH, considering the uncertainties of arterial diameters, stenosis parameters, and flow rates measured using the patients’ clinical data. A sensitivity analysis was performed to clarify the impact of each uncertain parameter on the flow rate increase. Results indicated that CH occurred when two conditions were satisfied simultaneously: severe stenosis and when arteries of small diameter serve as the collateral pathway to the cerebral artery on the stenosis side. These findings elucidate the biological aspects of cerebral circulation in terms of the relationship between collateral flow and CH.     

Virtual stenting workflow with vessel-specific initialization and adaptive expansion for neurovascular stents and flow diverters

Endovascular intervention using traditional neurovascular stents and densely braided flow diverters (FDs) have become the preferred treatment strategies for traditionally challenging intracranial aneurysms. Modeling stent and FD deployment in patient-specific aneurysms and its flow modification results prior to the actual intervention can potentially predict the patient outcome and treatment optimization. We present a clinically focused, streamlined virtual stenting workflow that efficiently simulates stent and FD treatment in patient-specific aneurysms based on expanding a simplex mesh structure. The simplex mesh is generated using an innovative vessel-specific initialization technique, which uses the patient’s parent artery diameter to identify the initial position of the simplex mesh inside the artery. A novel adaptive expansion algorithm enables the acceleration of deployment process by adjusting the expansion forces based on the distance of the simplex mesh from the parent vessel. The virtual stenting workflow was tested by modeling the treatment of two patient-specific aneurysms using the Enterprise stent and the Pipeline Embolization Device (commercial FD). Both devices were deployed in the aneurysm models in a few seconds. Computational fluid dynamics analyses of pre- and post-treatment aneurysmal hemodynamics show flow reduction in the aneurysmal sac in treated aneurysms, with the FD diverting more flow than the Enterprise stent. The test results show that this workflow can rapidly simulate clinical deployment of stents and FDs, hence paving the way for its future clinical implementation.     

An eFTD-VP framework for efficiently generating patient-specific anatomically detailed facial soft tissue FE mesh for craniomaxillofacial surgery simulation

Accurate surgical planning and prediction of craniomaxillofacial surgery outcome requires simulation of soft tissue changes following osteotomy. This can only be achieved by using an anatomically detailed facial soft tissue model. The current state-of-the-art of model generation is not appropriate to clinical applications due to the time-intensive nature of manual segmentation and volumetric mesh generation. The conventional patient-specific finite element (FE) mesh generation methods are to deform a template FE mesh to match the shape of a patient based on registration. However, these methods commonly produce element distortion. Additionally, the mesh density for patients depends on that of the template model. It could not be adjusted to conduct mesh density sensitivity analysis. In this study, we propose a new framework of patient-specific facial soft tissue FE mesh generation. The goal of the developed method is to efficiently generate a high-quality patient-specific hexahedral FE mesh with adjustable mesh density while preserving the accuracy in anatomical structure correspondence. Our FE mesh is generated by eFace template deformation followed by volumetric parametrization. First, the patient-specific anatomically detailed facial soft tissue model (including skin, mucosa, and muscles) is generated by deforming an eFace template model. The adaptation of the eFace template model is achieved by using a hybrid landmark-based morphing and dense surface fitting approach followed by a thin-plate spline interpolation. Then, high-quality hexahedral mesh is constructed by using volumetric parameterization. The user can control the resolution of hexahedron mesh to best reflect clinicians’ need. Our approach was validated using 30 patient models and 4 visible human datasets. The generated patient-specific FE mesh showed high surface matching accuracy, element quality, and internal structure matching accuracy. They can be directly and effectively used for clinical simulation of facial soft tissue change.     

Laporovaginal surgery in cervical cancer: a Croatian experience

With correct staging a large number of patients with cervical cancer FIGO stages IA2 and IB can be spared of unnecessary radiation therapy by laparoscopic assisted vaginal radical hysterectomy (LAVRH) as an option of radical surgical treatment in such patients. The development of laparovaginal surgery, indication and contraindication were presented. Also, the surgical technique was described in detail. Fifty-two patients were followed up in 2003 after LAVRH or open surgery, performed in our single center. Only 5 (14%) patients died from cervical cancer within 3 years following the treatment. They were all clinical stage IB treated with open surgery. There were 4 (11%) complications following treatment and they were all in patients with clinical stage IB, also treated with open surgery. There was no complication in LAVRH treated patients. The results and complications of the sole Croatian center performing LAVRH or open surgery in patients with cervical cancer FIGO stages IA and IB were similar to those in centers across the world.     

Personalized Computational Modeling of Mitral Valve Prolapse: Virtual Leaflet Resection

Posterior leaflet prolapse following chordal elongation or rupture is one of the primary valvular diseases in patients with degenerative mitral valves (MVs). Quadrangular resection followed by ring annuloplasty is a reliable and reproducible surgical repair technique for treatment of posterior leaflet prolapse. Virtual MV repair simulation of leaflet resection in association with patient-specific 3D echocardiographic data can provide quantitative biomechanical and physiologic characteristics of pre- and post-resection MV function. We have developed a solid personalized computational simulation protocol to perform virtual MV repair using standard clinical guidelines of posterior leaflet resection with annuloplasty ring implantation. A virtual MV model was created using 3D echocardiographic data of a patient with posterior chordal rupture and severe mitral regurgitation. A quadrangle-shaped leaflet portion in the prolapsed posterior leaflet was removed, and virtual plication and suturing were performed. An annuloplasty ring of proper size was reconstructed and virtual ring annuloplasty was performed by superimposing the ring and the mitral annulus. Following the quadrangular resection and ring annuloplasty simulations, patient-specific annular motion and physiologic transvalvular pressure gradient were implemented and dynamic finite element simulation of MV function was performed. The pre-resection MV demonstrated a substantial lack of leaflet coaptation which directly correlated with the severe mitral regurgitation. Excessive stress concentration was found along the free marginal edge of the posterior leaflet involving the chordal rupture. Following the virtual resection and ring annuloplasty, the severity of the posterior leaflet prolapse markedly decreased. Excessive stress concentration disappeared over both anterior and posterior leaflets, and complete leaflet coaptation was effectively restored. This novel personalized virtual MV repair strategy has great potential to help with preoperative selection of the patient-specific optimal MV repair techniques, allow innovative surgical planning to expect improved efficacy of MV repair with more predictable outcomes, and ultimately provide more effective medical care for the patient.     

Schwannomas of the peroneal nerves: Clinical and functional results of surgical treatment

Objectives:Peroneal nerves Schwannomas are rare benign tumors. Literature is still poor of studies about clinical and functional outcomes after surgical treatment. We evaluated the pre-operative presentation of the disease and assessed clinical and functional outcomes after surgery.Methods:We collected all the cases of peroneal nerves’ neurinoma treated surgically between June 2016 and June 2020. We analyzed each patients’ personal data and carried out accurate clinical examinations before and after surgery. MRI was performed both pre-operatively and post-operatively.Results:We reported 9 cases of peroneal nerves schwannomas: five arising from the common peroneal nerve and four arising from the deep or superficial branches alone. Their mean size was 22.6 mm. Each patient showed sensation deficits before surgery; pre-operative MRC score was 4.2. Pre-Operative MSTS and LEFS scores were 23.6 and 64.4. Surgery was successful in each case. No local recurrence nor major complication occurred. Tumor size was significantly associated with both diagnostic delay and development of pre-operative deficits. Surgery was proven to be globally successful: post-operative evaluations highlighted a marked reduction of neurological signs and overall functional limitations.Conclusions:Surgical treatment at early stages of the disease represents a reliable and relatively safe therapeutic option.