Advertisement
Congenital: Fontan| Volume 157, ISSUE 3, P1146-1155, March 2019

The first cohort of prospective Fontan surgical planning patients with follow-up data: How accurate is surgical planning?

Open ArchivePublished:December 11, 2018DOI:https://doi.org/10.1016/j.jtcvs.2018.11.102

      Abstract

      Objective

      Fontan surgical planning is an image-based, collaborative effort, which is hypothesized to result in improved patient outcomes. A common motivation for Fontan surgical planning is the progression (or concern for progression) of pulmonary arteriovenous malformations. The purpose of this study was to evaluate the accuracy of surgical planning predictions, specifically hepatic flow distribution (HFD), a known factor in pulmonary arteriovenous malformation progression, and identify methodological improvements needed to increase prediction accuracy.

      Methods

      Twelve single-ventricle patients who were enrolled in a surgical planning protocol for Fontan surgery with pre- and postoperative cardiac imaging were included in this study. Computational fluid dynamics were used to compare HFD in the surgical planning prediction and actual postoperative conditions.

      Results

      Overall, HFD prediction error was 17 ± 13%. This error was similar between surgery types (15 ± 18% and 18 ± 10% for revisions vs Fontan completions respectively; P = .73), but was significantly lower (6 ± 7%; P = .05) for hepatic to azygous shunts. Y-grafts and extracardiac conduits showed a strong correlation between prediction error and discrepancies in graft insertion points (r = 0.99; P < .001). Improving postoperative anatomy prediction significantly reduced overall HFD prediction error to 9 ± 6% (P = .03).

      Conclusions

      Although Fontan surgical planning can offer accurate HFD predictions for specific graft types, methodological improvements are needed to increase overall accuracy. Specifically, improving postoperative anatomy prediction was shown to be an important target for future work. Future efforts and refinements to the surgical planning process will benefit from an improved understanding of the current state and will rely heavily on increased follow-up data.

      Key Words

      Abbreviations and Acronyms:

      CMR (cardiac magnetic resonance), HFD (hepatic flow distribution), IVC (inferior vena cava), PAVM (pulmonary arteriovenous malformation), TCPC (total cavopulmonary connection)
      Figure thumbnail fx1
      Comparison of predicted and postoperative hepatic flow distribution (indicated by percentage) for a representative patient. The overlay shows predicted (yellow) and postoperative (purple) anatomies.
      Fontan surgical planning accurately predicts hepatic flow distribution for specific graft types and will benefit from improved anatomical predictions.
      Fontan surgical planning can provide accurate predictions of hepatic flow distribution for specific graft types. Anatomical prediction was shown to be a key methodological shortcoming in the surgical planning process. With continued improvements, surgical planning might be a useful tool to avoid pulmonary arteriovenous malformations in Fontan patients.
      See Commentary on page 1156
      Fontan surgical planning is an image-based, collaborative effort, which is hypothesized to result in improved patient outcomes. The general process has been detailed in multiple publications and has been implemented for select cases over the past decade.
      • Trusty P.M.
      • Slesnick T.C.
      • Wei Z.A.
      • Rossignac J.
      • Kanter K.R.
      • Fogel M.A.
      • et al.
      Fontan surgical planning: previous accomplishments, current challenges, and future directions.
      • Fogel M.A.
      • Khiabani R.H.
      • Yoganathan A.
      Imaging for preintervention planning pre- and post-fontan procedures.
      • Slesnick T.C.
      • Yoganathan A.P.
      Computational modeling of Fontan physiology: at the crossroads of pediatric cardiology and biomedical engineering.
      • de Zélicourt D.A.
      • Marsden A.
      • Fogel M.A.
      • Yoganathan A.P.
      Imaging and patient-specific simulations for the Fontan surgery: current methodologies and clinical applications.
      Despite its use in dozens of cases, few follow-up data have been available to evaluate the accuracy of surgical planning predictions.
      A common motivation for Fontan surgical planning is the progression (or concern for progression) of pulmonary arteriovenous malformations (PAVMs). PAVMs are extremely rare in the average population (2/100,000), but are much more common in Fontan patients in whom a poorly designed total cavopulmonary connection (TCPC) might lead to unbalanced hepatic flow distribution (HFD), a known factor in PAVM formation/progression.
      • Duncan B.W.
      • Desai S.
      Pulmonary arteriovenous malformations after cavopulmonary anastomosis.
      • Khurshid I.
      • Downie G.H.
      Pulmonary arteriovenous malformation.
      • Srivastava D.
      • Preminger T.
      • Lock J.E.
      • Mandell V.
      • Keane J.F.
      • Mayer J.E.
      • et al.
      Hepatic venous blood and the development of pulmonary arteriovenous malformations in congenital heart disease.
      • Bernstein H.S.
      • Brook M.M.
      • Silverman N.H.
      • Bristow J.
      Development of pulmonary arteriovenous fistulae in children after cavopulmonary shunt.
      • Shinohara T.
      • Yokoyama T.
      Pulmonary arteriovenous malformation in patients with total cavopulmonary shunt: what role does lack of hepatic venous blood flow to the lungs play?.
      • McElhinney D.B.
      • Marshall A.C.
      • Lang P.
      • Lock J.E.
      • Mayer J.E.
      Creation of a brachial arteriovenous fistula for treatment of pulmonary arteriovenous malformations after cavopulmonary anastomosis.
      • Vettukattil J.J.
      Pathogenesis of pulmonary arteriovenous malformations: role of hepatopulmonary interactions.
      This type of surgical planning can involve Fontan revision cases (a previous Fontan surgery resulted in PAVMs and therefore must be revised) as well as Fontan completion cases (stage 2-3, often performed in patients with more complex anatomies and hence the need for added insight to determine the best surgical strategy). In either case, the purpose of Fontan surgical planning is to determine which TCPC design will result in balanced HFD. The goal to optimize HFD has led to a variety of Fontan connections.
      Although limited to case studies and relatively short-term follow-up data, several previous studies have provided a preliminary understanding of surgical planning prediction accuracy.
      • Fogel M.A.
      • Khiabani R.H.
      • Yoganathan A.
      Imaging for preintervention planning pre- and post-fontan procedures.
      Sundareswaran and colleagues published the first use of surgical planning to correct PAVMs (determined by an increase in oxygen saturation) in 2009.
      • Sundareswaran K.S.
      • de Zélicourt D.
      • Sharma S.
      • Spray T.L.
      • Rossignac J.
      • Sotiropoulos F.
      • et al.
      Correction of pulmonary arteriovenous malformation using image-based surgical planning.
      Unfortunately, this single-patient case report included no postoperative imaging data to compare the predicted and postoperative HFD. Haggerty and colleagues provided the most thorough study to date, which included only 4 patients with short follow-up times (2 patients less than 1 month).
      • Haggerty C.M.
      • de Zélicourt D.A.
      • Restrepo M.
      • Rossignac J.
      • Spray T.L.
      • Kanter K.R.
      • et al.
      Comparing pre- and post-operative fontan hemodynamic simulations: implications for the reliability of surgical planning.
      In their study HFD and graft resistance were compared and showed “sufficient agreement” between the predicted and postoperative states. However, the use of steady flow conditions and fixed outlet flow splits, as well as the small sample size, limit the generalizability of these findings.
      The purpose of this study was to evaluate the accuracy of Fontan surgical planning predictions and identify which methodological improvements are needed to improve prediction accuracy. Specifically, in this study we focused on the prediction of HFD. Pre- and postoperative imaging data from patients enrolled in the surgical planning process were used to address these questions. This study offers, to our knowledge, the first assessment of prediction accuracy and methodological shortcomings for Fontan revisions and Fontan completion cases using longer-term follow-up data not previously available.

      Methods

      Patient Selection

      A total of 12 single-ventricle patients were included in this study. All patient data were received from the Children's Hospital of Philadelphia and Children's Healthcare of Atlanta under institutional review board approval (H17434 and H09279, respectively) with a waiver of consent. Inclusion criteria were that the patient: (1) was enrolled in the surgical planning process before surgery; (2) had pre- and postoperative cardiac magnetic resonance (CMR) and phase contrast CMR imaging; and (3) imaging quality was sufficient for accurate anatomical segmentation as well as flow segmentation at every TCPC inlet/outlet. This resulted in Fontan revision cases (n = 5) and Fontan completion cases (n = 7). All patients in this study were enrolled in surgical planning because of PAVM development or the concern for PAVM development due to atypical vasculature or clinical history. Clinical data including age, sex, body surface area, imaging and surgery dates, and diagnosis were obtained for each patient.

      Magnetic Resonance Imaging Acquisition

      All CMR imaging was performed with a Siemens 1.5 T magnetic resonance imaging system (Siemens Medical Solutions, Malvern, Pa). Patients were scanned supine, head first in the scanner with electrocardiogram leads placed. After localizers were obtained, a stack of contiguous, static, diastolic steady state free precession images were obtained from the diaphragm to thoracic inlet to assess anatomy and provide inputs for computational fluid dynamics modeling. Slice thickness was generally 3-4 mm and in plane resolution was 1 × 1 mm.
      Through plane, retrospectively gated, phase contrast magnetic resonance imaging was used to assess flow in the cavae, branch pulmonary arteries, and across the aortic valve. Inferior vena cava (IVC) flow was measured supra-hepatically. Velocity encoding was generally 150 cm/s for the aorta and 60 cm/s for the other vessels (superior vena cava, Fontan baffle, right pulmonary artery [RPA], and left pulmonary artery [LPA]). Slice thickness was generally 3 mm with in-plane resolution of 1 × 1 mm. The number of phases was a function of the heart rate and ranged from 20 to 30. Postoperative CMR images were acquired at the follow-up times specified in Table 1.
      Table 1Demographic and surgical data
      PatientSexDiagnosisSurgery typePAVMs
      Data indicate the lung that contained the malformations.
      Surgical optionAge at surgery, yearsAge at follow-up, yearsFollow-up, monthsHFD prediction, % LPAPostoperative HFD, % LPAHFD error
      1FH, PARevisionLeftHep to AZ4.712.897.860564
      2FHRevisionLeftY-graft19.019.00.338457
      3MH, DRevisionRightHep to AZ11.615.040.41001000
      4FURevisionRightY-graft12.712.81.6533221
      5MH, PARevisionLeftECC17.518.39.0267145
      6MHLHSFontan completionECC1.34.335.9175134
      7MH, HLHSFontan completionY-graft2.62.60.2274821
      8FPA, THFontan completionY-graft3.08.363.7607717
      9MPA, DILVFontan completionECC2.22.20.373807
      10FHFontan completionHep to AZ1.11.10.351565
      11FH, DORVFontan completionHep to AZ3.24.211.7718716
      12FH, HLHSFontan completionHep to Inn1.41.52.1482523
      PAVM, Pulmonary arteriovenous malformation; HFD, hepatic flow distribution; LPA, left pulmonary artery; F, female; H, heterotaxy; PA, pulmonary atresia; Hep to AZ, hepatic to azygous shunt; M, male; D, dextrocardia; U, unbalanced canal; ECC, extracardiac conduit; HLHS, hypoplastic left heart syndrome; TH, tricuspid hypoplasia; DILV, double-inlet left ventricle; DORV, double-outlet right ventricle; Hep to Inn, hepatic to innominate vein.
      Data indicate the lung that contained the malformations.

      Anatomic Reconstruction and Blood Flow Segmentation

      Patient-specific anatomies were reconstructed from axial CMR images using methods previously developed.
      • Martins T.
      • Barbara A.
      • Silva G.
      • Faria T.
      • Dalava B.
      • Silva J.
      InVesalius: three-dimensional medical reconstruction software.
      • de Moraes
      • Thiago F.
      • Amorim P.H.
      InVesalius-an open-source imaging application.
      Geomagic Studio (Geomagic Inc, Research Triangle Park, NC) was used to fit a surface around the reconstructed point-cloud and export the surface for mesh generation. Patient-specific blood flow waveforms were segmented from phase contrast magnetic resonance images for all vessels of interest using previously validated methods.
      • Heiberg E.
      • Sjögren J.
      • Ugander M.
      • Carlsson M.
      • Engblom H.
      • Arheden H.
      Design and validation of Segment - freely available software for cardiovascular image analysis.
      • Bidhult S.L.
      • Carlsson M.
      • Steding-Ehrenborg K.
      • Arheden H.
      • Heiberg E.
      A new method for vessel segmentation based on a priori input from medical expertise in cine phase-contrast magnetic resonance imaging.
      Pulmonary flow distribution was calculated as:
      PFD=QLPAQLPA+QRPA


      where Q is flow rate. Changes in flows from pre- to postoperative were calculated as:
      FlowChange=|PostoperativeQPreoperativeQ| to show absolute magnitude changes in flow rates.

      Surgical Planning Prediction

      For complex Fontan cases, surgical planning has been used to compare fundamentally different surgical options (extracardiac conduit vs Y-graft vs hepatic to azygous, etc). In addition, this process is used to determine which location on the pulmonary arteries the graft should be anastomosed. SURGEM III, a solid modeling software designed specifically for Fontan surgical planning was used to generate the surgical planning anatomy prediction.
      • Luffel M.
      • Sati M.
      • Rossignac J.
      • Yoganathan A.P.
      • Haggerty C.M.
      • Restrepo M.
      • et al.
      SURGEM: a solid modeling tool for planning and optimizing pediatric heart surgeries.
      First, the preoperative TCPC (for Fontan revisions) or bidirectional Glenn (for Fontan completion) and hepatics were imported. With input from the respective clinician, the desired surgical option was then created and exported as a surface mesh. Preoperative flow waveforms reconstructed from phase contrast–magnetic resonance images were directly used as the “predicted” flow waveforms. This technique makes the simplifying assumption that postoperative flows will be identical to preoperative flows. It is important to mention that during the surgical planning process, multiple surgical options are created. However, for the sake of comparing surgical planning predictions with postoperative data, only the prediction of the actual implemented surgical option is considered.

      Computational Fluid Dynamics

      The 3-D anatomy (from either the preoperative scan, surgical planning prediction, or postoperative scan) was imported into ANSYS Workbench (ANSYS Inc, Canonsburg, Pa), where vessel extensions of length 10 × (vessel diameter) were added to overcome entrance effects and establish an appropriate velocity profile. A polyhedral mesh of approximately DIVC/20-mm elements was used to achieve mesh-independent results, where DIVC is the diameter of the IVC.
      • Wei Z.A.
      • Tree M.
      • Trusty P.M.
      • Wu W.
      • Singh-Gryzbon S.
      • Yoganathan A.
      The advantages of viscous dissipation rate over simplified power loss as a fontan hemodynamic metric.
      All simulations were performed using ANSYS Fluent (Release 17.1; ANSYS Inc), which is a finite volume pressure-based Navier–Stokes solver. Blood was modeled as a single-phase Newtonian fluid (μ = 0.04 g/cm/s; ρ = 1.06 g/cm3). The appropriate patient-specific blood flow waveforms extracted from phase contrast--magnetic resonance imaging were used as boundary conditions for each TCPC inlet and outlet.
      • Wei Z.A.
      • Tree M.
      • Trusty P.M.
      • Wu W.
      • Singh-Gryzbon S.
      • Yoganathan A.
      The advantages of viscous dissipation rate over simplified power loss as a fontan hemodynamic metric.
      • Wei Z.A.
      • Trusty P.M.
      • Tree M.
      • Haggerty C.M.
      • Tang E.
      • Fogel M.
      • Yoganathan A.P.
      Can time-averaged flow boundary conditions be used to meet the clinical timeline for Fontan surgical planning?.
      Twenty cardiac cycles were simulated to overcome transition effects and achieve period stability, using the final cycle for data analysis.
      To investigate the potential accuracy of surgical planning if methodological improvements allowed for more accurate anatomy and flow predictions, 2 additional simulations were run representing either “improved” anatomy or flow predictions. This resulted in 4 simulations for each patient:
      • 1.
        Predicted (simulation uses predicted anatomy and predicted flows);
      • 2.
        Actual postoperative (simulation uses postoperative anatomy and postoperative flows);
      • 3.
        “Improved” anatomy prediction (simulation uses actual postoperative anatomy and predicted flows). The postoperative anatomy represents a perfect anatomical prediction from the surgical planning process; and
      • 4.
        “Improved” flow prediction (simulation uses predicted anatomy and actual postoperative flows). The postoperative flows represent a perfect flow prediction from the surgical planning process.
      In addition, preoperative simulations (preoperative anatomy and preoperative flows) were also run for Fontan revision cases. This allowed an investigation of the relationship between preoperative HFD and PAVM progression.

      HFD

      HFD was quantified by seeding massless particles at the IVC and calculating the total flux of particles leaving the left and right pulmonary arteries. HFD was defined as:
      HFD=θLPAθLPA+θRPA


      where θ is the total flux of particles throughout a cardiac cycle. The error in HFD prediction is defined as:
      ErrorinHFDprediction=|HFDpost-opHFDpredicted|.


      Anatomy Comparison

      To compare the predicted and postoperative anatomies, the TCPCs were first registered to account for differences in imaging coordinate systems. A mesh comparison software (CloudCompare, version 2.10) was then used to quantify average and maximum deviations between the surfaces of the 2 TCPCs. This was done for the full TCPC and the graft alone. Graft insertion offset was calculated by measuring the distance between the predicted and postoperative anastomosis locations (distance between center points of each anastomosis). For Y-grafts, the largest insertion offset of the 2 branches was used.

      Statistical Analysis

      SPSS version 25 (IBM Corp, Armonk, NY) was used for statistical analyses. The Shapiro–Wilk test was used to determine normality for each parameter. Pearson and Spearman correlations were used to investigate bivariate correlations between HFD prediction error and clinical, hemodynamic, and anatomic parameters for linear and monotonic relationships, respectively. Depending on normality, either a Wilcoxon rank sum test or a 2-sample t test was used to test for equal medians between the revision and Fontan completion groups, as well as between various surgical connection types. A paired sample t test was used to test for differences between surgical planning methodologies (current vs improved anatomy/flow). Statistical significance was determined using P < .05. Values are shown as average ± standard deviation (median, interquartile range).

      Results

      Clinical Data

      Clinical and surgical data are given in Table 1. The cohort consisted of 5 Fontan revisions and 7 Fontan completion cases. Implemented surgical options included 4 hepatic to azygous shunts, 4 Y-grafts, 3 traditional extracardiac conduits, and 1 hepatic to innominate vein connection. Average follow-up time was 22 ± 32 [6, 39] months.

      Revisions Versus Fontan Completions

      Age at surgery and age at follow-up were significantly different between the revision (13.1 ± 5.7 [12.7, 10.1] and 15.6 ± 2.9 [15.0, 5.8] years) and Fontan completion (2.1 ± 0.9 [2.2, 1.7] and 3.5 ± 2.5 [2.6, 2.8] years) cases respectively (P < .001 for both; Table 2). Follow-up time was not significantly different between the revision and Fontan completion cases (30 ± 41 [9, 68] and 16 ± 25 [2, 64] months respectively; P = .49, Table 2). No significant differences in pre- to postoperative changes in flow rates were seen between the revision and Fontan completion cases (Table 2). Additionally, no flow rates (grouped according to surgery type or vessel) showed consistent directionality in flow rate changes. Significant differences in geometric variations (between the predicted and actual postoperative anatomy) were seen between the revision and Fontan completion cases in terms of the TCPC as a whole and the graft alone (Table 2).
      Table 2Fontan revision and Fontan completion comparison
      VariableRevisionFontan completionP value
      Age at surgery, y13.1 ± 5.7 (12.7-10.1)2.1 ± 0.9 (2.2-1.7)<.001
      Age at follow-up, y15.6 ± 2.9 (15.0-5.8)3.5 ± 2.5 (2.6-2.8)<.001
      Follow-up time, mo30 ± 41 (9-68)16 ± 25 (2-64).492
      HFD prediction error15 ± 18 (7-31)18 ± 10 (17-16).795
      IVC flow change, L/min0.22 ± 0.16 (0.15-0.32)0.25 ± 0.19 (0.17-0.33).782
      SVC flow change, L/min0.35 ± 0.13 (0.4-0.24)0.44 ± 0.32 (0.34-0.56).573
      AZ flow change, L/min0.61 ± 0.45 (0.63-0.82)0.20 ± 0.17 (0.2-0.32).104
      LPA flow change, L/min0.57 ± 0.19 (0.57-0.32)0.48 ± 0.37 (0.35-0.67).638
      RPA flow change, L/min0.80 ± 0.60 (0.91-1.13)0.25 ± 0.22 (0.14-0.39).052
      IVC flow change, %42 ± 42 (35-62)41 ± 26 (32-35).955
      SVC flow change, %36 ± 21 (35-40)48 ± 29 (52-44).462
      AZ flow change, %49 ± 30 (51-58)40 ± 28 (48-52).656
      LPA flow change, %38 ± 21 (30-41)59 ± 35 (64-38).274
      RPA flow change, %62 ± 40 (55-58)41 ± 53 (19-32).471
      PFD change, %9 ± 9 (4-14)14 ± 6 (13-12).261
      TCPC deviation, mm3.3 ± 0.8 (3-1.4)1.6 ± 0.5 (1.5-0.9).001
      TCPC max deviation, mm17.6 ± 2.0 (19-3.5)9.2 ± 3.7 (9.4-4.1).004
      Graft deviation, mm5.6 ± 2.8 (6.3-5.1)2.8 ± 1.9 (2.7-3.0).062
      Graft max deviation, mm15.3 ± 4.4 (15.5-8.3)9.3 ± 4.5 (9.4-6.2).046
      Graft insertion offset, mm13.5 ± 8.3 (15.0-16.1)5.8 ± 4.5 (6.8-7.5).063
      Data are presented as mean ± standard deviation (interquartile range). Change represents absolute difference between preoperative and postoperative flows. HFD, Hepatic flow distribution; IVC, inferior vena cava; SVC, superior vena cava; AZ, azygous vein; LPA, left pulmonary artery; RPA, right pulmonary artery; PFD, pulmonary flow distribution; TCPC, total cavopulmonary connection.

      HFD Prediction Error

      The predicted and postoperative HFD can be found in Table 1. Overall, the HFD prediction error was 17 ± 13 [17, 17]%, and was not significantly different between revisions (15 ± 18 [7, 31]%) and Fontan completion (18 ± 10 [17, 16]%) cases (P = .73; Table 2). Fontan completion predictions underestimated HFD in 6 of 7 cases, whereas revisions were evenly split between overestimations and underestimations. Computational fluid dynamics results comparing the predicted and postoperative streamlines for all Fontan revisions and Fontan completion cases are shown in Figures 1 and 2, respectively. In addition, Figure 1 shows the preoperative HFD for all Fontan revisions, confirming a lack of hepatic flow to the lung with PAVMs. Overall, no significant correlations were found between HFD prediction error and age at surgery, age at follow-up, or follow-up time. Moderate correlations were seen between the percent change in IVC flow rate (r = 0.60; P = .04) and the change in pulmonary flow distribution (r = 0.60; P = .04) from pre- to postoperative with HFD prediction error.
      Figure thumbnail gr1
      Figure 1Streamline comparison for the preoperative (Pre-Op), predicted, and postoperative (Post-Op) states for all Fontan revision cases. Patients 1, 2, and 5 were diagnosed with left-sided pulmonary arteriovenous malformations (PAVMs), and patients 3 and 4 with right-sided PAVMs. PAVMs regressed in each case when the revision resulted in increased hepatic flow to the affected lung. Hepatic flow distribution is noted as the percent of hepatic flow distribution to the left pulmonary artery (LPA). SVC, Superior vena cava; IN, innominate vein; RPA, right pulmonary artery; AZ, azygous vein; FC, Fontan conduit; LSVC, left superior vena cava.
      Figure thumbnail gr2
      Figure 2Streamline comparison between the predicted and postoperative (Post-Op) states for all Fontan completion patients. Hepatic flow distribution is noted as the percent of hepatic flow distribution to the left pulmonary artery (LPA). The overlay column shows a comparison of the predicted (yellow) and Post-Op (purple) anatomies. SVC, Superior vena cava; RPA, right pulmonary artery; FC, Fontan conduit; AZ, azygous vein; LSVC, left superior vena cava; IN, innominate vein.

      Connection Types

      HFD prediction error was associated with surgical connection type. A comparison of HFD prediction error between graft types is shown in Figure 3. Hepatic to azygous shunts had significantly lower prediction errors than other connection types (6 ± 7 [5, 12]% vs 22 ± 13 [21, 22]% respectively; P = .05; Figure 3, B). In addition, a strong, positive correlation was seen between HFD prediction error and graft insertion offset within the Y-graft and extracardiac conduit groups (r = 0.99; P < .001; Figure 4, D). Example cases of low, moderate, and high graft insertion offsets are shown in Figure 4, A-C.
      Figure thumbnail gr3
      Figure 3Effect of connection type on hepatic flow distribution (HFD) prediction error. A, High variability in prediction error was observed within and across graft types. B, Hepatic to azygous (Hep to AZ) shunts showed significantly lower HFD prediction errors than other connection types. ECC, Extracardiac conduit.
      Figure thumbnail gr4
      Figure 4Relationship between hepatic flow distribution (HFD) prediction error and graft insertion offset for extracardiac conduit and Y-graft Fontan connections. Graft insertion offset describes the distance between the predicted and implemented graft insertion locations. Representative cases are shown for (A) low, (B) moderate, and (C) high graft insertion offsets. A strong correlation (D) was seen between prediction error and offset for these connection types. The overlay figures (panels A-C) show a comparison of the predicted (blue) and postoperative (red) total cavopulmonary connections on the left, and show a color map of the offset between the predicted and postoperative grafts on the right for each representative case. All cases use the same color scale.

      Potential Methodological Improvements

      The HFD prediction error associated with the “improved” postoperative anatomy and flow predictions and the current methodology is shown in Table 3. Although a reduced or identical prediction error was seen for most patients using either methodological improvement (8 of 12 and 7 of 12 for the improved anatomy and improved flow scenarios, respectively), a more substantial reduction in error was seen when the improved anatomy was used. The current HFD prediction error (17 ± 13 [17, 17]%) was significantly reduced by improving anatomy prediction (9 ± 6 [9, 12]%; P = .03, paired sample t test), but remained nearly the same when only an improved flow prediction was used (18 ± 17 [14, 24]%; P = .73, paired sample t test). When the 2 potential methodological improvements were compared, improved postoperative anatomy prediction resulted in equivalent or more accurate HFD predictions for 9 of 12 patients compared with improved flow prediction.
      Table 3Comparison of HFD prediction errors between the current surgical planning process and potential methodological improvements
      PatientSurgery typeVessels presentConnection typeHFD prediction error
      CurrentImproved anatomy predictionImproved flow prediction
      1RevisionAZ, INNHep to AZ4611
      2RevisionAZY-graft7131
      3RevisionAZ, LSVCHep to AZ000
      4RevisionY-graft211029
      5RevisionAZECC451658
      6Fontan completionECC341818
      7Fontan completionAZY-graft21738
      8Fontan completionAZY-graft17614
      9Fontan completionLSVCECC7135
      10Fontan completionAZHep to AZ523
      11Fontan completionAZ, LSVCHep to AZ16323
      12Fontan completionAZ, INNHep to Inn231813
      Average17 ± 13 (17-17)9 ± 6 (9-12)
      Statistically significant difference (paired sample t test) from the current method.
      18 ± 17 (14-24)
      The presence of any nonstandard vessels (in addition to the inferior/superior vena cava and left/right pulmonary artery) is indicated in the vessels present column. HFD, Hepatic flow distribution; AZ, azygous shunt; INN, innominate vein; Hep to AZ, hepatic to azygous shunt; LSVC, left superior vena cava; ECC, extracardiac conduit; Hep to Inn, hepatic to innominate connection.
      Statistically significant difference (paired sample t test) from the current method.

      Discussion

      Previously limited by a lack of postoperative data, the current study offers, to our knowledge, the first assessment of prospective Fontan surgical planning accuracy for Fontan revisions and Fontan completions using medium-term postoperative data. This study incorporates a unique data set resulting from more than a decade of surgical planning experience and provides a methodological assessment necessary for the improvement of surgical planning accuracy.
      Although an exact “cutoff” for HFD to prevent PAVMs is currently unknown, and might vary between patients, the Fontan revision results (Table 1, Figure 1) confirm a lack of hepatic flow to the lung with PAVMs in each case. Furthermore, PAVMs regressed in each case when the revision resulted in increased hepatic flow to the affected lung. In combination with previous studies, these results emphasize the importance of achieving a balanced HFD through appropriate TCPC design.
      In this study, HFD prediction error averaged 17 ± 13 [17, 17]% across all patients. Prediction error was similar between revision and Fontan completion cases, but differed across connection types. Intuitively, hepatic to azygous connections are more robust to variations in surgical implementation because all hepatic flow will join the azygous flow and then travel through the entire azygous vein before interacting with other flows regardless of exact placement of the shunt. No colliding flows from multiple vessels are present locally in hepatic to azygous connections, in contrast with Y-graft and extracardiac conduit connections where slight offsets might substantially change the interactions between various inlets and therefore stray from predicted results (Figure 4, D).
      • Restrepo M.
      • Luffel M.
      • Sebring J.
      • Kanter K.
      • del Nido P.
      • Veneziani A.
      • et al.
      Surgical planning of the total cavopulmonary connection: robustness analysis.
      Therefore, hepatic to azygous predictions were quite accurate, whereas predictions for other connection types were more varied.

      Methodologic Improvements

      Capitalizing on the available postoperative data, various simulations were run using the postoperative anatomy or flows as a surrogate for an “improved” prediction. Although it is unlikely that anatomy or flow prediction techniques will ever produce exact matches to postoperative outcomes, this analysis is instructive by showing the full potential of surgical planning accuracy if methodological improvements in either of these areas offered extremely accurate predictions. To reiterate, improvements in anatomy prediction led to a significant (P = .03) reduction in HFD prediction error. Interestingly, improvements in flow prediction did not result in similar error reduction (Table 3). These findings stress 2 critical points: (1) postoperative anatomy prediction is a primary factor in HFD prediction; and (2) anatomy prediction methods must improve for Fontan surgical planning to provide more accurate HFD predictions.
      Again, improved flow prediction alone did not result in more accurate HFD predictions (P = .73; Table 3). In general, HFD is primarily driven by graft placement.
      • Tang E.
      • Restrepo M.
      • Haggerty C.M.
      • Mirabella L.
      • Bethel J.
      • Whitehead K.K.
      • et al.
      Geometric characterization of patient-specific total cavopulmonary connections and its relationship to hemodynamics.
      In a complex connection such as the TCPC, relatively small offsets in graft placement and angulation can largely alter the collisions and interactions among flows from various vessels.
      • Tang E.
      • Restrepo M.
      • Haggerty C.M.
      • Mirabella L.
      • Bethel J.
      • Whitehead K.K.
      • et al.
      Geometric characterization of patient-specific total cavopulmonary connections and its relationship to hemodynamics.
      These variations can affect the preferential streaming of inlet flows including hepatic flow, which in turn will determine HFD. Naturally, severe changes in individual flow rates can affect HFD prediction; however, this was not observed in this cohort.
      Though it is common to use indexed flow rates (normalized for body surface area) in pediatric studies involving changes over time, raw flow rates are shown in this study (Table 2) to emphasize the changes in actual input to the surgical planning process. Although an indexed flow rate might remain constant over a several-year follow-up, the actual flow rate (and therefore the flow rate that needs to be predicted) does not. These raw data better represent how boundary conditions for the surgical planning process change over time.
      Multiple methods exist to predict postoperative Fontan anatomies and flows, ranging from simple (basic computer-aided design software and using preoperative flows as the “predicted” postoperative flows) to more sophisticated (designated surgical planning software and lumped parameter modeling) methods.
      • Luffel M.
      • Sati M.
      • Rossignac J.
      • Yoganathan A.P.
      • Haggerty C.M.
      • Restrepo M.
      • et al.
      SURGEM: a solid modeling tool for planning and optimizing pediatric heart surgeries.
      • Marsden A.
      • Esmaily-Moghadam M.
      Multiscale modeling of cardiovascular flows for clinical decision support.
      As the methods have progressed in complexity, anatomy prediction methods have become faster (software is designed specifically for Fontan surgical planning) and flow prediction methods have become slower (more complex calculations and “full body” modeling). Meeting the clinical timeline for most surgical planning cases requires accelerated analysis.
      • Trusty P.M.
      • Slesnick T.C.
      • Wei Z.A.
      • Rossignac J.
      • Kanter K.R.
      • Fogel M.A.
      • et al.
      Fontan surgical planning: previous accomplishments, current challenges, and future directions.
      Conveniently, the present results indicate that anatomy prediction is a primary shortcoming, which can hopefully be improved without lengthening the surgical planning process.
      Accurate anatomy prediction involves predicting a viable surgical option and accurately implementing that option. Modeling a viable surgical option is heavily dependent on high-quality imaging data, clinician involvement, and inclusion of relevant anatomical landmarks. Current methods include the heart, aorta, and pulmonary circulation, but future efforts could potentially use additional organs and the process of chest closure. When a surgeon has selected the surgical option to implement, closely replicating that option in vivo might be challenging. Little intraoperative guidance is currently offered as part of the surgical planning process. Some efforts have explored 3-D printing and augmented/virtual reality as planning/guidance tools, but further refinements are needed.
      • Sodian R.
      • Weber S.
      • Markert M.
      • Loeff M.
      • Lueth T.
      • Weis F.C.
      • et al.
      Pediatric cardiac transplantation: three-dimensional printing of anatomic models for surgical planning of heart transplantation in patients with univentricular heart.
      • Ventola C.L.
      Medical applications for 3D printing: current and projected uses.
      • McCloy R.
      • Stone R.
      Science, medicine, and the future. Virtual reality in surgery.
      • Seymour N.E.
      • Gallagher A.G.
      • Roman S.A.
      • O’Brien M.K.
      • Bansal V.K.
      • Andersen D.K.
      • et al.
      Virtual reality training improves operating room performance: results of a randomized, double-blinded study.
      In addition, growth is another difficult factor to model that might be necessary to improve surgical planning results.
      • Restrepo M.
      • Mirabella L.
      • Tang E.
      • Haggerty C.M.
      • Khiabani R.H.
      • Fynn-Thompson F.
      • et al.
      Fontan pathway growth: a quantitative evaluation of lateral tunnel and extracardiac cavopulmonary connections using serial cardiac magnetic resonance.
      • Restrepo M.
      • Tang E.
      • Haggerty C.M.
      • Khiabani R.H.
      • Mirabella L.
      • Bethel J.
      • et al.
      Energetic implications of vessel growth and flow changes over time in Fontan patients.
      It is possible that “variations” between a predicted and postoperative anatomy are because of growth rather than “imperfect” surgical implementation.
      Finally, patient 3 is a unique case that deserves attention. As shown in Figure 1, this patient had 100% HFD to the left lung in the preoperative state. This led to right-sided PAVMs and the need for a Fontan revision. Also shown in Figure 1, the surgical option implemented during the Fontan revision similarly resulted in 100% HFD to the left lung. This was predicted during the surgical planning process but implemented nonetheless because of other concerns. Importantly, not all surgical options for patient 3 resulted in 100% HFD to the left lung. This case illustrates the importance of being able to predict which options will perform poorly in addition to which ones perform well. This case also shows that surgical planning predictions are not always the main determinant for decision-making in current clinical practice. Assessing how surgical planning affects clinical decision-making and patient outcomes is an important next step when surgical planning predictions are known to offer accurate results.
      If Fontan surgical planning is to be used in clinical practice, the importance and necessity of follow-up data and validation studies such as this cannot be overstated. Understanding the current accuracy and methodological shortcomings is imperative to correctly use the results and progress the surgical planning process. Future efforts and refinements to the surgical planning process will greatly benefit from an improved understanding of the current state of the art (Video 1).
      Figure thumbnail fx2
      Video 1Discussion of surgical planning predictions, accuracy, and implications by Dr Timothy Slesnick. Video available at: https://www.jtcvs.org/article/S0022-5223(18)33252-5/fulltext.

      Limitations

      Although important conclusions can be drawn from this study, a more complete understanding of surgical planning accuracy and methodological needs with increased statistical power will require a substantial amount of data most likely through a multicenter study. Additionally, the current study includes a broad range of follow-up times and offers only a “snapshot” of the patients’ postoperative hemodynamics. If available, the inclusion of serial data in a similar study could offer a better understanding of hemodynamic changes over time. The predicted results in this study are representative of the specific surgical planning process used. Results might vary with other prediction techniques. However, this study used one of the most advanced anatomy prediction methods and we conclude that postoperative anatomy prediction is a limiting factor, which is unlikely to change on the basis of prediction technique. Finally, 2 experienced surgeons were involved in these surgical planning cases. Although we saw no differences in prediction accuracy between the 2 surgeons, it is possible that results might vary on the basis of the surgeon involved. Finally, the results from this study do not indicate how the surgical strategy changed because of the surgical planning process. A surgical strategy was not developed before the surgical planning process for each case to determine how surgical planning influenced the original plan. Therefore, these results report prediction accuracy and do not quantify the effect on clinical decision-making.

      Conclusions

      Overall, HFD prediction error was 17 ± 13%. This error was similar between Fontan revisions and Fontan completions, but varied across surgical connection types. Although Fontan surgical planning can offer accurate HFD predictions for specific graft types, methodological improvements are needed to increase overall accuracy. Specifically, improving postoperative anatomy prediction was shown to be an important target for future efforts that would substantially improve flow field modeling, and therefore increase HFD prediction accuracy. Future efforts and refinements to the surgical planning process will greatly benefit from an improved understanding of the current state of the art and will rely heavily on increased follow-up data.

      Conflict of Interest Statement

      Mark A. Fogel reports a modest research grant (Siemens). All other authors have nothing to disclose with regard to commercial support.
      The authors acknowledge the use of ANSYS software, which was provided through an Academic Partnership between ANSYS, Inc, and the Cardiovascular Fluid Mechanics Lab at the Georgia Institute of Technology. We also acknowledge Dr Jarek Rossignac from the College of Computing at Georgia Tech for his contributions to the development of the SURGEM III platform.

      Supplementary Data

      References

        • Trusty P.M.
        • Slesnick T.C.
        • Wei Z.A.
        • Rossignac J.
        • Kanter K.R.
        • Fogel M.A.
        • et al.
        Fontan surgical planning: previous accomplishments, current challenges, and future directions.
        J Cardiovasc Transl Res. 2018; 11: 133-144
        • Fogel M.A.
        • Khiabani R.H.
        • Yoganathan A.
        Imaging for preintervention planning pre- and post-fontan procedures.
        Circ Cardiovasc Imaging. 2013; 6: 1092-1101
        • Slesnick T.C.
        • Yoganathan A.P.
        Computational modeling of Fontan physiology: at the crossroads of pediatric cardiology and biomedical engineering.
        Int J Cardiovasc Imaging. 2014; 30: 1073-1084
        • de Zélicourt D.A.
        • Marsden A.
        • Fogel M.A.
        • Yoganathan A.P.
        Imaging and patient-specific simulations for the Fontan surgery: current methodologies and clinical applications.
        Prog Pediatr Cardiol. 2010; 30: 31-44
        • Duncan B.W.
        • Desai S.
        Pulmonary arteriovenous malformations after cavopulmonary anastomosis.
        Ann Thorac Surg. 2003; 76: 1759-1766
        • Khurshid I.
        • Downie G.H.
        Pulmonary arteriovenous malformation.
        Postgrad Med J. 2002; 78: 191-197
        • Srivastava D.
        • Preminger T.
        • Lock J.E.
        • Mandell V.
        • Keane J.F.
        • Mayer J.E.
        • et al.
        Hepatic venous blood and the development of pulmonary arteriovenous malformations in congenital heart disease.
        Circulation. 1995; 92: 1217-1222
        • Bernstein H.S.
        • Brook M.M.
        • Silverman N.H.
        • Bristow J.
        Development of pulmonary arteriovenous fistulae in children after cavopulmonary shunt.
        Circulation. 1995; 92: II309-II314
        • Shinohara T.
        • Yokoyama T.
        Pulmonary arteriovenous malformation in patients with total cavopulmonary shunt: what role does lack of hepatic venous blood flow to the lungs play?.
        Pediatr Cardiol. 2001; 22: 343-346
        • McElhinney D.B.
        • Marshall A.C.
        • Lang P.
        • Lock J.E.
        • Mayer J.E.
        Creation of a brachial arteriovenous fistula for treatment of pulmonary arteriovenous malformations after cavopulmonary anastomosis.
        Ann Thorac Surg. 2005; 80: 1604-1609
        • Vettukattil J.J.
        Pathogenesis of pulmonary arteriovenous malformations: role of hepatopulmonary interactions.
        Heart. 2002; 88: 561-563
        • Sundareswaran K.S.
        • de Zélicourt D.
        • Sharma S.
        • Spray T.L.
        • Rossignac J.
        • Sotiropoulos F.
        • et al.
        Correction of pulmonary arteriovenous malformation using image-based surgical planning.
        JACC Cardiovasc Imaging. 2009; 2: 1024-1030
        • Haggerty C.M.
        • de Zélicourt D.A.
        • Restrepo M.
        • Rossignac J.
        • Spray T.L.
        • Kanter K.R.
        • et al.
        Comparing pre- and post-operative fontan hemodynamic simulations: implications for the reliability of surgical planning.
        Ann Biomed Eng. 2012; 40: 1-13
        • Martins T.
        • Barbara A.
        • Silva G.
        • Faria T.
        • Dalava B.
        • Silva J.
        InVesalius: three-dimensional medical reconstruction software.
        Virtual Rapid Manuf. 2007; 3: 135-141
        • de Moraes
        • Thiago F.
        • Amorim P.H.
        InVesalius-an open-source imaging application.
        Comput Vis Med Image Process. 2011; 3: 405-407
        • Heiberg E.
        • Sjögren J.
        • Ugander M.
        • Carlsson M.
        • Engblom H.
        • Arheden H.
        Design and validation of Segment - freely available software for cardiovascular image analysis.
        BMC Med Imaging. 2010; 10: 1
        • Bidhult S.L.
        • Carlsson M.
        • Steding-Ehrenborg K.
        • Arheden H.
        • Heiberg E.
        A new method for vessel segmentation based on a priori input from medical expertise in cine phase-contrast magnetic resonance imaging.
        J Cardiovasc Magn Reson. 2014; 16: P355
        • Luffel M.
        • Sati M.
        • Rossignac J.
        • Yoganathan A.P.
        • Haggerty C.M.
        • Restrepo M.
        • et al.
        SURGEM: a solid modeling tool for planning and optimizing pediatric heart surgeries.
        Comput Des. 2015; 70: 3-12
        • Wei Z.A.
        • Tree M.
        • Trusty P.M.
        • Wu W.
        • Singh-Gryzbon S.
        • Yoganathan A.
        The advantages of viscous dissipation rate over simplified power loss as a fontan hemodynamic metric.
        Ann Biomed Eng. 2018; 46: 404-416
        • Wei Z.A.
        • Trusty P.M.
        • Tree M.
        • Haggerty C.M.
        • Tang E.
        • Fogel M.
        • Yoganathan A.P.
        Can time-averaged flow boundary conditions be used to meet the clinical timeline for Fontan surgical planning?.
        J Biomech. 2017; 50: 172-179
        • Restrepo M.
        • Luffel M.
        • Sebring J.
        • Kanter K.
        • del Nido P.
        • Veneziani A.
        • et al.
        Surgical planning of the total cavopulmonary connection: robustness analysis.
        Ann Biomed Eng. 2015; 43: 1321-1334
        • Tang E.
        • Restrepo M.
        • Haggerty C.M.
        • Mirabella L.
        • Bethel J.
        • Whitehead K.K.
        • et al.
        Geometric characterization of patient-specific total cavopulmonary connections and its relationship to hemodynamics.
        JACC Cardiovasc Imaging. 2014; 7: 215-224
        • Marsden A.
        • Esmaily-Moghadam M.
        Multiscale modeling of cardiovascular flows for clinical decision support.
        Appl Mech Rev. 2015; 67: 1-11
        • Sodian R.
        • Weber S.
        • Markert M.
        • Loeff M.
        • Lueth T.
        • Weis F.C.
        • et al.
        Pediatric cardiac transplantation: three-dimensional printing of anatomic models for surgical planning of heart transplantation in patients with univentricular heart.
        J Thorac Cardiovasc Surg. 2008; 136: 1098-1099
        • Ventola C.L.
        Medical applications for 3D printing: current and projected uses.
        P T. 2014; 39: 704-711
        • McCloy R.
        • Stone R.
        Science, medicine, and the future. Virtual reality in surgery.
        BMJ. 2001; 323: 912-915
        • Seymour N.E.
        • Gallagher A.G.
        • Roman S.A.
        • O’Brien M.K.
        • Bansal V.K.
        • Andersen D.K.
        • et al.
        Virtual reality training improves operating room performance: results of a randomized, double-blinded study.
        Ann Surg. 2002; 236: 454-458
        • Restrepo M.
        • Mirabella L.
        • Tang E.
        • Haggerty C.M.
        • Khiabani R.H.
        • Fynn-Thompson F.
        • et al.
        Fontan pathway growth: a quantitative evaluation of lateral tunnel and extracardiac cavopulmonary connections using serial cardiac magnetic resonance.
        Ann Thorac Surg. 2014; 97: 916-922
        • Restrepo M.
        • Tang E.
        • Haggerty C.M.
        • Khiabani R.H.
        • Mirabella L.
        • Bethel J.
        • et al.
        Energetic implications of vessel growth and flow changes over time in Fontan patients.
        Ann Thorac Surg. 2015; 99: 163-170

      Linked Article

      • Commentary: Individualized surgery for patients with single-ventricle physiology
        The Journal of Thoracic and Cardiovascular SurgeryVol. 158Issue 5
        • Preview
          Six decades of treating patients with single-ventricle physiology has taught us that the complexity and heterogeneity of the Fontan circulation does not lend itself to a “one size fits-all” strategy. Furthermore, we have learned that the technical details associated with the Fontan-Kreutzer (Fontan) operation can have a lasting impact on a physiologic single ventricle. The development of pulmonary arteriovenous malformations (PAVMs) is an important example of how the technical aspects of a Fontan operation—or similar total caval pulmonary connections—can influence long-term outcomes.
        • Full-Text
        • PDF
        Open Archive
      • Another step toward intelligent surgical planning
        The Journal of Thoracic and Cardiovascular SurgeryVol. 157Issue 3
        • Preview
          The applications of computational fluid dynamics (CFD) to cardiovascular medicine and surgery just keep getting better. Trusty and colleagues1 examine a way that CFD could be used to plan “complex” Fontan operations. These patients are those whose anatomy does not “guarantee” a relatively even distribution of hepatic blood flow (HBF) to both lungs with a conventional extracardiac conduit or lateral tunnel. The typical example would be an interrupted inferior vena cava with azygous continuation to a left superior vena cava and a right-sided hepatic veins-to-right pulmonary artery conduit.
        • Full-Text
        • PDF
        Open Archive