It has been 50 years since Francis Fontan pioneered the operation that today bears his name. Initially designed for patients with tricuspid atresia, this procedure is now offered for a vast array of congenital cardiac lesions when a circulation with 2 ventricles cannot be achieved. As a result of technical advances and improvements in patient selection and perioperative management, survival has steadily increased, and it is estimated that patients operated on today may hope for a 30-year survival of >80%. Up to 70 000 patients may be alive worldwide today with Fontan circulation, and this population is expected to double in the next 20 years. In the absence of a subpulmonary ventricle, Fontan circulation is characterized by chronically elevated systemic venous pressures and decreased cardiac output. The addition of this acquired abnormal circulation to innate abnormalities associated with single-ventricle congenital heart disease exposes these patients to a variety of complications. Circulatory failure, ventricular dysfunction, atrioventricular valve regurgitation, arrhythmia, protein-losing enteropathy, and plastic bronchitis are potential complications of the Fontan circulation. Abnormalities in body composition, bone structure, and growth have been detected. Liver fibrosis and renal dysfunction are common and may progress over time. Cognitive, neuropsychological, and behavioral deficits are highly prevalent. As a testimony to the success of the current strategy of care, the proportion of adults with Fontan circulation is increasing. Healthcare providers are ill-prepared to tackle these challenges, as well as specific needs such as contraception and pregnancy in female patients. The role of therapies such as cardiovascular drugs to prevent and treat complications, heart transplantation, and mechanical circulatory support remains undetermined. There is a clear need for consensus on how best to follow up patients with Fontan circulation and to treat their complications. This American Heart Association statement summarizes the current state of knowledge on the Fontan circulation and its consequences. A proposed surveillance testing toolkit provides recommendations for a range of acceptable approaches to follow-up care for the patient with Fontan circulation. Gaps in knowledge and areas for future focus of investigation are highlighted, with the objective of laying the groundwork for creating a normal quality and duration of life for these unique individuals.

Background-In the Fontan circulation, pulmonary and systemic vascular resistances are in series. The implications of this unique arrangement on infradiaphragmatic venous physiology are poorly understood. Methods and Results-We studied the effects of respiration and gravity on infradiaphragmatic venous flows in 20 normal healthy volunteers (control) and 48 Fontan patients (atriopulmonary connection [APC] nϭ15, total cavopulmonary connection [TCPC] nϭ30). Hepatic venous (HV), subhepatic inferior vena caval (IVC), and portal venous (PV) flow rates were measured with Doppler ultrasonography during inspiration and expiration in both the supine and upright positions. The inspiratory-to-expiratory flow rate ratio was calculated to reflect the effect of respiration, and the supine-to-upright flow rate ratio was calculated to assess the effect of gravity. HV flow depended heavily on inspiration in TCPC compared with both control and APC subjects (inspiratory-to-expiratory flow rate ratio 3.4, 1.7, and 1.6, respectively; PϽ0.0001). Normal PV flow was higher in expiration, but this effect was lost in TCPC and APC patients (inspiratory-to-expiratory flow rate ratio 0.8, 1.0, and 1.1, respectively; Pϭ0.01). The respiratory influence on IVC flow was the same in all groups. Gravity decreased HV flow more in APC than in TCPC patients (supine-to-upright flow rate ratio 3.2 versus 2.1, respectively; PϽ0.04) but reduced PV flow equally in all groups. Conclusions-Gravity and respiration have important influences on infradiaphragmatic venous return in Fontan patients.Although gravity exerts a significant detrimental effect on lower body venous return, which is more marked in APC than in TCPC patients, the beneficial effects of respiration in TCPC patients are mediated primarily by an increase in HV flow. These effects may have important short-and long-term implications for the hemodynamics of the Fontan circulation.

In patients who are in functionally poorer condition after the Fontan operation, portal venous flow loses normal expiratory augmentation and adverse gravity influence is enhanced. These suboptimal flow dynamics, coupled with higher splanchnic venous pressures and lower transhepatic venous pressure gradients, suggest that hepatic sinusoids are congested, acting as "open tubes." Transhepatic gradient loss is incrementally worse with higher caval pressures. These observations may be responsible for late gastrointestinal problems in patients who have had the Fontan operation.

International audienceSimulation divergence due to backflow is a common, but not fully addressed, problem in three-dimensional simulations of blood flow in the large vessels. Because backflow is a naturally occurring physiologic phenomenon, careful treatment is necessary to realistically model backflow without artificially altering the local flow dynamics. In this study, we quantitatively compare three available methods for treatment of outlets to prevent backflow divergence in finite element Navier-Stokes solvers. The methods examined are (1) adding a stabilization term to the boundary nodes formulation, (2) constraining the velocity to be normal to the outlet, and (3) using Lagrange multipliers to constrain the velocity profile at all or some of the outlets. A modification to the stabilization method is also discussed. Three model problems, a short and long cylinder with an expansion, a right-angle bend, and a patient-specific aorta model, are used to evaluate and quantitatively compare these methods. Detailed comparisons are made to evaluate robustness, stability characteristics, impact on local and global flow physics, computational cost, implementation effort, and ease-of-use. The results show that the stabilization method offers a promising alternative to previous methods, with reduced effect on both local and global hemodynamics, improved stability, little-to-no increase in computational cost, and elimination of the need for tunable parameters

BackgroundMedical image analysis in clinical practice is commonly carried out on 2D image data, without fully exploiting the detailed 3D anatomical information that is provided by modern non-invasive medical imaging techniques. In this paper, a statistical shape analysis method is presented, which enables the extraction of 3D anatomical shape features from cardiovascular magnetic resonance (CMR) image data, with no need for manual landmarking. The method was applied to repaired aortic coarctation arches that present complex shapes, with the aim of capturing shape features as biomarkers of potential functional relevance. The method is presented from the user-perspective and is evaluated by comparing results with traditional morphometric measurements.MethodsSteps required to set up the statistical shape modelling analyses, from pre-processing of the CMR images to parameter setting and strategies to account for size differences and outliers, are described in detail. The anatomical mean shape of 20 aortic arches post-aortic coarctation repair (CoA) was computed based on surface models reconstructed from CMR data. By analysing transformations that deform the mean shape towards each of the individual patient’s anatomy, shape patterns related to differences in body surface area (BSA) and ejection fraction (EF) were extracted. The resulting shape vectors, describing shape features in 3D, were compared with traditionally measured 2D and 3D morphometric parameters.ResultsThe computed 3D mean shape was close to population mean values of geometric shape descriptors and visually integrated characteristic shape features associated with our population of CoA shapes. After removing size effects due to differences in body surface area (BSA) between patients, distinct 3D shape features of the aortic arch correlated significantly with EF (r = 0.521, p = .022) and were well in agreement with trends as shown by traditional shape descriptors.ConclusionsThe suggested method has the potential to discover previously unknown 3D shape biomarkers from medical imaging data. Thus, it could contribute to improving diagnosis and risk stratification in complex cardiac disease.Electronic supplementary materialThe online version of this article (doi:10.1186/s12880-016-0142-z) contains supplementary material, which is available to authorized users.

The objective of this work is to perform a virtual planning of surgical repairs in patients with congenital heart diseases-to test the predictive capability of a closed-loop multi-scale model. As a first step, we reproduced the pre-operative state of a specific patient with a univentricular circulation and a bidirectional cavopulmonary anastomosis (BCPA), starting from the patient's clinical data. Namely, by adopting a closed-loop multi-scale approach, the boundary conditions at the inlet and outlet sections of the three-dimensional model were automatically calculated by a lumped parameter network. Successively, we simulated three alternative surgical designs of the total cavopulmonary connection (TCPC). In particular, a T-junction of the venae cavae to the pulmonary arteries (T-TCPC), a design with an offset between the venae cavae (O-TCPC) and a Y-graft design (Y-TCPC) were compared. A multi-scale closed-loop model consisting of a lumped parameter network representing the whole circulation and a patient-specific three-dimensional finite volume model of the BCPA with detailed pulmonary anatomy was built. The three TCPC alternatives were investigated in terms of energetics and haemodynamics. Effects of exercise were also investigated. Results showed that the pre-operative caval flows should not be used as boundary conditions in post-operative simulations owing to changes in the flow waveforms post-operatively. The multi-scale approach is a possible solution to overcome this incongruence. Power losses of the Y-TCPC were lower than all other TCPC models both at rest and under exercise conditions and it distributed the inferior vena cava flow evenly to both lungs. Further work is needed to correlate results from these simulations with clinical outcomes.*Author for correspondence (giancarlo.pennati@polimi.it).One contribution of 11 to a Theme Issue 'Towards the virtual physiological human: mathematical and computational case studies'.This journal is

A closed-loop lumped parameter model of blood circulation is considered for single-ventricle shunt physiology. Its parameters are estimated by an inverse problem based on patient-specific haemodynamics measurements. As opposed to a black-box approach, maximizing the number of parameters that are related to physically measurable quantities motivates the present model. Heart chambers are described by a single-fibre mechanics model, and valve function is modelled with smooth opening and closure. A model for valve prolapse leading to valve regurgitation is proposed. The method of data assimilation, in particular the unscented Kalman filter, is used to estimate the model parameters from time-varying clinical measurements. This method takes into account both the uncertainty in prior knowledge related to the parameters and the uncertainty associated with the clinical measurements. Two patient-specific cases - one without regurgitation and one with atrioventricular valve regurgitation - are presented. Pulmonary and systemic circulation parameters are successfully estimated, without assumptions on their relationships. Parameters governing the behaviour of heart chambers and valves are either fixed based on biomechanics, or estimated. Results of the inverse problem are validated qualitatively through clinical measurements or clinical estimates that were not included in the parameter estimation procedure. The model and the estimation method are shown to successfully capture patient-specific clinical observations, even with regurgitation, such as the double peaked nature of valvular flows and anomalies in electrocardiogram readings. Lastly, biomechanical implications of the results are discussed.

The adoption of simulation tools to predict surgical outcomes is increasingly leading to questions about the variability of these predictions in the presence of uncertainty associated with the input clinical data. In the present study, we propose a methodology for full propagation of uncertainty from clinical data to model results that, unlike deterministic simulation, enables estimation of the confidence associated with model predictions. We illustrate this problem in a virtual stage II single ventricle palliation surgery example. First, probability density functions (PDFs) of right pulmonary artery (PA) flow split ratio and average pulmonary pressure are determined from clinical measurements, complemented by literature data. Starting from a zero-dimensional semi-empirical approximation, Bayesian parameter estimation is used to find the distributions of boundary conditions that produce the expected PA flow split and average pressure PDFs as pre-operative model results. To reduce computational cost, this inverse problem is solved using a Kriging approximant. Second, uncertainties in the boundary conditions are propagated to simulation predictions. Sparse grid stochastic collocation is employed to statistically characterize model predictions of post-operative hemodynamics in models with and without PA stenosis. The results quantify the statistical variability in virtual surgery predictions, allowing for placement of confidence intervals on simulation outputs.