The present invention relates to the field of cardiac disease, in particular to the assessment of coronary vessels, in particular to determine the patterns of blockage or restriction to the blood flow through a coronary vessel. More particularly, the present invention relates to a computer-implemented method to quantify the extent of functional coronary artery disease. In addition, the invention provides a computer device for determining the functional pattern of coronary disease in a mammal. More particularly the invention provides methods to select a mammal suffering from coronary disease to benefit from a percutaneous coronary intervention.
Invasive functional assessment of coronary artery disease (CAD) has been regarded as the gatekeeper for revascularization in patients with chronic coronary syndromes. Guidelines advocate evaluating the reduction in coronary flow using pressure-derived indices to decide upon the need for revascularization. Intracoronary pressure measurements are typically performed in the distal segment of the coronary artery reflecting cumulative pressure losses along the epicardial vessel. Focal narrowing can be entirely responsible for the pressure drops; nonetheless, diffuse functional deterioration can be also observed outside angiographic stenotic regions contributing to the total decrease in coronary perfusion pressure. Coronary angiography remains to date the most utilized method to guide stent implantation. The length of the lesion can be quantified by quantitative coronary angiography (QCA) or alternatively, more precisely, using intravascular imaging. Both approaches aim to guide stent selection to cover the atherosclerotic plaque, restore epicardial conductance and improve myocardial perfusion. However, in almost a third of patients after an angiographically successful percutaneous coronary intervention (PCI), epicardial conductance remains suboptimal. In diffuse functional disease, PCI is of limited benefit in terms of coronary physiology whereas in focal CAD PCI usually restores epicardial conductance. Furthermore, patients with low fractional flow reserve (FFR) after percutaneous revascularization have been shown to be at an increased risk of adverse events compared to patients with high post-PCI FFR. Gain in epicardial conductance with PCI can be predicted by assessing the distribution of epicardial resistance. A pullback maneuver during intracoronary pressure measurements identifies the presence, location, magnitude and extent of pressure drops. Two factors, namely (i) the magnitude of FFR drops and (ii) extension of functional CAD are predictive of improvement in epicardial conductance after percutaneous revascularization. Thus, quantifying the extent of functional disease may have prognostic capability for post-PCI FFR.
In the present invention, we sought to quantify the mismatch in the extent of CAD between anatomical and functional evaluations and to assess the impact of functional-anatomical mismatch (FAM) on FFR after PCI.
According to a first aspect, there is provided a computer device for quantifying the extent of functional coronary artery disease (CAD) comprising a processor configured to:
This approach is that it is less vulnerable to artefacts in the pullback curves compared to prior art devices that do not make use of piece-wise linearization by means of an automated change point detection algorithm. It is clear that piece-wise linearization by means of an automated change point detection algorithm, by means of the parameters determined based on the linear segments is less sensitive to local artefacts, which are for example the result of temporary or local measurement errors, etc. Furthermore, when use is made of this improved quantification in the context of FAM, as explained in further detail below, which quantifies the mismatch between the anatomical and functional CAD length, and it also becomes possible to take into account the impact of residual pressure losses outside the treated region on post PCI physiology. Moreover, the FAM approach is based on the presence and length of disease rather than on the magnitude of pressure drops making this approach less influenced by the interaction in cases of serial lesions. Further this allows an improved assessment of the functional pattern of CAD which may improve patient selection for PCI by avoiding stenting lesions which don't result in sufficient post-PCI benefits, by reducing the risk of peri-procedural myocardial infarction and by increasing the chance of a net clinical benefit from revascularization. In this way patients with a negative FAM, i.e. having diffuse functional CAD, may be better treated with optimal medical therapy or coronary artery bypass grafting, and patients with a positive FAM may be better treated with PCI.
It is clear that said FFR data refers to said set of fractional flow reserve values.
According to an embodiment, there is provided a computer device, configured to:
According to an embodiment, there is provided a computer device, configured to:
According to an embodiment, there is provided a computer device, configured to:
According to an embodiment there is provided a computer device configured to quantify the extent of the functional coronary artery disease (CAD), or in other words the function lesion length, based on and/or such that it corresponds to:
According to an embodiment the FFR data comprises the set of FFR values obtained at different positions of a coronary vessel between the ostium and the most distal part of the coronary vessel.
According to an embodiment the set of FFR values corresponds to an FFR pullback curve.
According to an embodiment there is provided a computer device is configured to quantify the functional lesion length from analysis of the FFR pullback curve.
According to an embodiment the computer device is further configured to perform the quantification of the functional lesion length and/or the classification of the coronary vessel in segments after one or more of the following:
According to an embodiment, there is provided a computer device, wherein the computer device further comprises a display configured to display said healthy, focal and/or diffused diseased fragments, optionally on an image of the coronary artery, optionally wherein the displayed image of the coronary artery is a 2-dimensional image.
According to an embodiment, there is provided a computer device, wherein the automated change-point detection algorithm is configured to detect one or more change points in the set of FFR values, such that said change points each correspond to a position along the coronary vessel where an attribute of the set of FFR values changes, wherein:
According to an embodiment, there is provided a computer device, wherein:
According to an embodiment, there is provided a computer device, wherein the computer device is further configured to classify the coronary vessel such that:
According to an embodiment, the logistic regression model is configured to provide a binary separation.
According to an embodiment, the computer device is further configured to apply the logistic regression model in a two-steps approach, in which:
According to an embodiment, there is provided a computer device, wherein said automated change-points detection algorithm is configured to operate based on a penalized parametric global method.
According to an embodiment, there is provided a computer device, wherein the display is further configured to display the image of the coronary artery in a 2-dimensional image.
According to an embodiment, there is provided a computer device, further configured to obtain the set of FFR values from:
According to an embodiment, there is provided a computer device, wherein the computer device is further configured to predict the response to a percutaneous coronary intervention (PCI) by said quantifying of the extent of functional CAD, and/or wherein the computer device is further configured to quantify the extent of functional CAD as the sum of the lengths of the diseased segments.
According to an embodiment, there is provided a computer device, wherein the computer device is further configured to select a mammal suffering from coronary artery disease (CAD) to be eligible for a percutaneous coronary intervention (PCI) by said quantifying of the extent of functional CAD, and selecting a mammal when the extent of functional disease in the coronary artery is smaller than the extent of anatomical disease in the coronary artery; and/or wherein the computer device is further configured to calculate a Functional Anatomical Mismatch (FAM) as the difference between the extent of anatomical CAD and the extent of functional, thereby identifying two lesion endotypes: (1) functional CAD circumscribed within the anatomical CAD when FAM>0, and (2) functional CAD extending beyond the anatomical CAD when FAM<0.
According to an embodiment, there is provided a computer device, wherein the computer device (system) is configured to operate offline.
According to an embodiment, there is provided a computer device, wherein the computer device is configured to perform said automatic classification.
According to a second aspect, there is provided a computer-implemented method to quantify the extent of functional coronary artery disease (CAD) comprising the following steps: i) processing a set of fractional flow reserve (FFR) values obtained at different positions of a coronary vessel between the ostium and the most distal part of the coronary vessel, ii) classifying the coronary vessel in healthy segments, focal diseased segments and/or diffused diseased segments by carrying out a piece-wise linearization of said FFR data by applying an automated change-points detection algorithm, and optionally iii) displaying said healthy, focal and/or diffused diseased fragments on an image of the coronary artery, and optionally said automated change-points detection algorithm is based on a penalized parametric global method.
According to an embodiment, there is provided a computer-implemented method, wherein said method comprises the step of obtaining the set of FFR values from a pull-back curve, or 3-dimensional quantitative coronary angiography, or a CT scan, or intravascular imaging (e.g. optical coherence tomography (OCT) or intravascular ultrasound (IVUS), or the combination between coronary angiography and intravascular imaging or the combination of a CT scan and intravascular imaging.
According to a third aspect, there is provided a computer-implemented method for developing an automated classifier for use in the computer device according to the first aspect for performing the classification of the coronary vessel in healthy focal and/or diffused diseased segments and/or for use in the computer-implemented method according to the second aspect for performing the classification of the coronary vessel in healthy focal and/or diffused diseased segments, wherein:
The mismatch between anatomy and physiology regarding epicardial lesion severity has been widely recognized in the prior art. For example, in the FAME study, more than one-third of lesions with an angiographic 50% to 70% diameter stenosis demonstrated an FFR 50.80 whereas one-fifth of lesions with a 71% to 90% angiographic diameter stenosis demonstrated an FFR>0.80. Disconnection between anatomy and physiology goes beyond the assessment of lesion significance. The length of CAD also differs between anatomical and functional evaluations. In the present invention, we have determined the length of functional CAD lesion with the means of a specially developed automatic algorithm. Accordingly, our novel computer-implemented method shows that when the length of functional disease (in a coronary artery of a patient) is greater than its anatomical equivalent either derived from QCA or optical coherence tomography (OCT) then the FAM value is <0 (see
In one aspect the invention relates to a computer-implemented method to quantify the extent of functional coronary artery disease (CAD) comprising the following steps: i) processing a set of fractional flow reserve (FFR) values obtained at different positions of the coronary vessel between the ostium and the most distal part of the coronary vessel, ii) classifying the coronary vessel in healthy, focal and/or diffused diseased segments by carrying out a piece-wise linearization of said FFR data by applying an automated change-points detection algorithm and optionally iii) displaying said healthy, focal and/or diffused diseased fragments on a 2-dimensional image of the coronary artery.
In specific aspects in the computer-implemented method the FFR values are obtained from a pull-back curve, 3-dimensional quantitative coronary angiography, CT scan or OCT. It is thus clear that the set of FFR values, or in other words the FFR data or FFR pullback curve, according to an embodiment can be obtained as data that was measured, generated and/or recorded from measurements of suitable pressure sensors during an FFR pullback operation, or in other words FFR data obtained from pressure measurements in the coronary artery vessel, which is an invasive measurement. It is however clear that preferably the embodiment of the computer implemented method does not include the invasive step of making the pressure measurements in the coronary artery vessel, and preferably only processes data received as an input, resulting from such measurements. According to alternative embodiments, the set of FFR values, or in other words the FFR data or FFR pullback curve does not result from direct pressure measurements inside the coronary vessel but is calculated by means of computational fluid dynamics simulations applied to a 3D model of the coronary vessel as reconstructed from medical imaging, such as for example 3-dimensional quantitative coronary angiography, CT scan, OCT or IVUS. It is clear that according to such an embodiment, the FFR data can be obtained by means of non-invasive measurements, such as for example 3-dimensional quantitative coronary angiography, CT scan. When according to such embodiments, there is made use of invasive measurements, such as for example OCT or IVUS, it is clear that preferably the embodiment of the computer implemented method does not include the invasive step of making the measurements in the coronary artery vessel, and preferably only processes data received as an input, resulting from such measurements, and preferably the medical imaging data from these measurements, or a 3 dimensional model of the coronary vessel as reconstructed from such medical imaging data.
In another aspect an in vitro method is provided to predict the response to a percutaneous coronary intervention (PCI) by quantifying the extent of functional CAD.
In yet another aspect an in vitro method is provided to select a mammal suffering from coronary artery disease (CAD) to be eligible for a percutaneous coronary intervention (PCI) comprising the application of the computer-implemented method described herein and selecting a mammal when the extent of functional disease in the coronary artery is smaller than the extent of anatomical disease in the coronary artery.
In a specific aspect the method is an offline method.
In yet another aspect the method to quantify the extent of functional coronary artery disease is an automatic classification method.
In another aspect a computer device is provided for evaluating the functional pattern of coronary artery disease in a mammal, said computer device configured to process a set of FFR values obtained at different positions of the coronary vessel between the ostium and the most distal part of the coronary vessel and classifying the coronary vessel in focal and/or diffused diseased segments by carrying out a piece-wise linearization of said FFR data by applying an automated change-points detection algorithm.
According to a further aspect, there is provided a computer-implemented method to quantify the extent of functional coronary artery disease (CAD) comprising the following steps: i) processing a set of fractional flow reserve (FFR) values obtained at different positions of the coronary vessel between the ostium and the most distal part of the coronary vessel, ii) classifying the coronary vessel in healthy, focal and/or diffused diseased segments by carrying out a piece-wise linearization of said FFR data by applying an automated change-points detection algorithm and optionally iii) displaying said healthy, focal and/or diffused diseased fragments on an image of the coronary artery.
According to an embodiment, there is provided a computer-implemented method wherein said automated change-points detection algorithm is based on a penalized parametric global method.
According to an embodiment, there is provided a computer-implemented method wherein in step iii) the displayed image of the coronary artery is a 2-dimensional image.
According to an embodiment, there is provided a computer-implemented method wherein the FFR values are obtained from a pull-back curve or 3-dimensional quantitative coronary angiography or CT scan or intravascular imaging (e.g. optical coherence tomography (OCT) or intravascular ultrasound (IVUS) or the combination between coronary angiography and intravascular imaging or the combination of a CT scan and intravascular imaging.
According to a further aspect, there is provided a method to predict the response to a percutaneous coronary intervention (PCI) by quantifying the extent of functional CAD according to the previous aspect.
According to a further aspect, there is provided a method to select a mammal suffering from coronary artery disease (CAD) to be eligible for a percutaneous coronary intervention (PCI) comprising the application of the computer-implemented method according to a previous aspect and selecting a mammal when the extent of functional disease in the coronary artery is smaller than the extent of anatomical disease in the coronary artery.
According to an embodiment, there is provided a method according to a previous aspect wherein the method is an offline method.
According to an embodiment, there is provided a method according to a previous aspect wherein the method is an automatic classification method.
According to a further aspect, there is provided a computer device for evaluating the functional pattern of coronary artery disease in a mammal, said computer device configured to process a set of FFR values obtained at different positions of the coronary vessel between the ostium and the most distal part of the coronary vessel and classifying the coronary vessel in healthy, focal and/or diffused diseased segments by carrying out a piece-wise linearization of said FFR data by applying an automated change-points detection algorithm, optionally said automated change-points detection algorithm is based on a penalized parametric global method.
According to an embodiment, there is provided a computer device wherein the set of FFR values are obtained from a pull-back curve, or 3-dimensional quantitative coronary angiography, or a CT scan, or intravascular imaging (e.g. optical coherence tomography (OCT) or intravascular ultrasound (IVUS), or the combination between coronary angiography and intravascular imaging or the combination of a CT scan and intravascular imaging.
Exemplary embodiments will now be described, for example with reference to the following Figures. When in the context of these Figures, and/or the description, there is made reference to colors as present in the drawings as appended to the filing of this patent application. These colors have been translated into corresponding indications by means of different line styles in the following Figures.
The present invention will be described with respect to particular embodiments and with reference to certain drawings, but the invention is not limited thereto but only by the claims. Any reference signs in the claims shall not be construed as limiting the scope. The drawings described are only schematic and are non-limiting. In the drawings, the size of some of the elements may be exaggerated and not drawn on scale for illustrative purposes. Where the term “comprising” is used in the present description and claims, it does not exclude other elements or steps. Where an indefinite or definite article is used when referring to a singular noun e.g. “a” or “an”, “the”, this includes a plural of that noun unless something else is specifically stated. Furthermore, the terms first, second, third and the like in the description and in the claims, are used for distinguishing between similar elements and not necessarily for describing a sequential or chronological order. It is to be understood that the terms so used are interchangeable under appropriate circumstances and that the embodiments of the invention described herein are capable of operation in other sequences than described or illustrated herein. The following terms or definitions are provided solely to aid in the understanding of the invention. Unless specifically defined herein, all terms used herein have the same meaning as they would to one skilled in the art of the present invention. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art (e.g. in molecular biology, interventional cardiology fluid physics, biochemistry, and/or computational biology/biomechanics).
Functional complete revascularization of the coronary arteries has been associated with improved clinical outcomes after percutaneous coronary interventions (PCI). Nevertheless, in one third of patients, coronary perfusion pressure remains low even after a successful procedure. In the present invention we developed a computer-implemented method which was able to quantify the mismatch in the extent of CAD between anatomical and functional invasive evaluations. In the invention we have assessed the impact of this mismatch on post-PCI fractional flow reserve (FFR).
Accordingly, the invention provides in a first embodiment a computer-implemented method to quantify the extent of functional coronary artery disease (CAD) comprising the following steps: i) processing a set of fractional flow reserve (FFR) values obtained at different positions of the coronary vessel between the ostium and the most distal part of the coronary vessel, ii) classifying the coronary vessel in focal and/or diffused diseased segments by carrying out a piece-wise linearization of said FFR data by applying an automated change-points detection algorithm based on a penalized parametric global method and optionally, and iii) displaying said focal and/or diffused diseased fragments on a 2-dimensional image of the coronary artery.
In a particular embodiment the FFR values are obtained from a pull-back curve, 3-dimensional quantitative coronary angiography, CT scan or optical coherence tomography (OCT).
Automatic change-point detection algorithms are known in the art but have never been used in the context of cardiac disease. Several of such algorithms are known in the art. For example, the binary segmentation method proposed by Scott, A. J. and Knott, M. (1974) Biometrics, 30 (3):507-512 is a well-known changepoint search method. Yet another approach is based on the algorithm of Jackson B. et al (2005) IEEE, Signal Processing Letters, 12(2):105-108. In the present invention the method of Killick R et al (2012) J. Am. Stat. Assoc. 107, 1590-1598) has been used.
In yet another embodiment the invention provides an in vitro method for predicting the response to a percutaneous coronary intervention (PCI) by quantifying the extent of functional CAD.
In yet another embodiment the invention provides a method to select a mammal suffering from coronary artery disease (CAD) to be eligible for a percutaneous coronary intervention (PCI) comprising the application of the computer-implemented method according to the methods herein described and selecting a mammal when the extent of functional disease in the coronary artery is smaller than the extent of anatomical disease in the coronary artery.
In a particular embodiment the mammal is eligible for a PCI when the functional-anatomical mismatch (FAM) is higher than 0 (FAM=>0). This means that the length of the functional CAD is described within the anatomical defined CAD.
In another embodiment the mammal is not eligible for a PCI (or not eligible for having a successful PCI procedure) when the FAM is lower than 0 (FAM=<0). This means that the length of the functional CAD extends beyond the anatomical defined lesion. This means that a negative FAM value reflects to pressure losses outside the determined anatomical lesion.
In specific embodiments the methods are offline methods.
In other specific embodiments the methods are automatic classification methods.
In yet another embodiment the invention provides a computer device for evaluating the functional pattern of coronary artery disease in a mammal, said computer device configured to process a set of FFR values obtained at different positions of the coronary vessel between the ostium and the most distal part of the coronary vessel and classifying the coronary vessel in focal and/or diffused diseased segments by carrying out a piece-wise linearization of said FFR data by applying an automated change-points detection algorithm based on a penalized parametric global method.
In a specific embodiment the set of FFR values are obtained from a pull-back curve, 3-dimensional quantitative coronary angiography, CT scan or optical coherence tomography (OCT).
The methods of the invention can also be used in veterinary applications. Mammals comprise cats, dogs, horses, cows, goats, sheep and preferably human subjects (human patients).
In other specific embodiment the fractional flow reserve (FFR) data are obtained by a manual or motorized pullback device which device is attached to the pressure wire.
In yet another particular embodiment there is no need for a motorized pullback device but instead the FFR curve is obtained by a pressure wire comprising a multiple of built-in pressure sensors.
In another particular embodiment, the catheter is configured to obtain diagnostic information about the coronary vessel. In this respect, the catheter can include one or more sensors, transducers, and/or other monitoring elements configured to obtain the diagnostic information about the vessel. The diagnostic information includes one or more of pressure, flow (velocity), images (including images obtained using ultrasound (e.g., intravascular ultrasound—IVUS), optical coherence tomography (OCT), thermal, and/or other imaging techniques), temperature, and/or combinations thereof. These one or more sensors, transducers, and/or other monitoring elements are positioned less than 30 cm, less than 10 cm, less than 5 cm, less than 3 cm, less than 2 cm, and/or less than 1 cm from a distal tip of the catheter in some instances. In some instances, at least one of the one or more sensors, transducers, and/or other monitoring elements is positioned at the distal tip of the catheter. In another particular embodiment the catheter comprises at least one element configured to monitor pressure within the coronary vessel. The pressure monitoring element can take the form a piezo-resistive pressure sensor, a piezo-electric pressure sensor, a capacitive pressure sensor, an electromagnetic pressure sensor, an optical pressure sensor, and/or combinations thereof. In some instances, one or more features of the pressure monitoring element are implemented as a solid-state component manufactured using semiconductor and/or other suitable manufacturing techniques.
In yet another embodiment the catheter comprises a pressure wire (or a guide wire). Examples of commercially available guide wire products that include suitable pressure monitoring elements include, without limitation, the Prime Wire PRESTIGE® pressure guide wire, the Prime Wire® pressure guide wire, and the ComboWire® XT pressure and flow guide wire, each available from Volcano Corporation, as well as the Pressure Wire™ Certus guide wire and the Pressure Wire™ Aeris guide wire, each available from St. Jude Medical, Inc or COMET™ FFR pressure guidewire from Boston Scientific. The pressure wire is also configured to obtain diagnostic information about the coronary vessel. In some instances, the pressure wire is configured to obtain the same diagnostic information as the catheter. In other instances, the pressure wire is configured to obtain different diagnostic information than the catheter, which may include additional diagnostic information, less diagnostic information, and/or alternative diagnostic information. The diagnostic information obtained by the pressure wire includes one or more of pressure, flow (velocity), images (including images obtained using ultrasound (e.g. IVUS), OCT, thermal, and/or other imaging techniques), temperature, and/or combinations thereof.
Similar to the catheter the pressure wire also includes at least one element configured to monitor pressure within the vessel. The pressure monitoring element can take the form a piezo-resistive pressure sensor, a piezo-electric pressure sensor, a capacitive pressure sensor, an electromagnetic pressure sensor, an optical pressure sensor, and/or combinations thereof. In some instances, one or more features of the pressure monitoring element are implemented as a solid-state component manufactured using semiconductor and/or other suitable manufacturing techniques. In a particular embodiment the pressure wire can comprise multiple pressure sensors, e.g. at least 10, at least 20, at least 30, at least 40, at least 50, or more pressure sensors. It is clear that according to such embodiments of the pressure wire, the multiple pressure sensors are provided at different positions along the length of the pressure wire, and thus configured to, even when stationary, after being introduced into the coronary vessel up to the distal end of the coronary vessel, determine a plurality of pressure measurements at different positions along the length of the coronary vessel, or in other words at different positions between the ostium and the distal end of the coronary vessel.
In a particular embodiment the pressure wire is configured to monitor pressure within the vessel while being moved through the lumen of the vessel. In some instances, the pressure wire is configured to be moved through the lumen of the vessel and across the stenosis present in the vessel. In that regard, the pressure wire is positioned distal of the stenosis and moved proximally (i.e. pulled back) across the stenosis to a position proximal of the stenosis in some instances. Movement of the pressure wire can be controlled manually by medical personnel (e.g. hand of a surgeon) in some embodiments. In other preferred embodiments, movement of the pressure wire is controlled automatically by a movement control device (e.g. a pullback device, such as the Trak Back® II or Volcano R-100 Device available from Volcano Corporation). In that regard, the movement control device controls the movement of the pressure wire at a selectable and known speed (e.g. 5.0 mm/s, 2.0 mm/s, 1.0 mm/s, 0.5 mm/s, etc.) in some instances. Movement of the pressure wire through the vessel is continuous for each pullback, in some instances. In other instances, the pressure wire is moved step-wise through the vessel (i.e. repeatedly moved a fixed amount of distance and/or a fixed amount of time).
In yet another embodiment the invention provides a system for evaluating coronary artery disease in a patient under hyperaemic conditions, comprising i) a coronary catheter comprising a pressure sensor, said catheter further comprising a pressure wire comprising at least one pressure sensor, ii) a computing device in communication with the catheter and the pressure wire, the computing device configured to generate an FFR curve based on a multiple of FFR values (the latter are relative pressure measurements from pressures obtained over the total length of the coronary vessel relative to the pressure in the ostium), iii) said computer device comprising a computer algorithm which calculates the length (or the extent) of the a functional coronary disease based on the FFR curve, and the correlation with the anatomical coronary artery disease, the computer output displays an FAM value which informs an interventional cardiologist of a successful percutaneous coronary intervention based on the positive or negative value of the FAM value. In particular, when the FAM value is negative (i.e. <0) there is no likelihood of conducing a successful PCI on the patient.
In the present invention a “system” is equivalent to a “device” or an “apparatus”. It is clear that such a system, device and/or apparatus could comprise any suitable number of interconnected subsystems or components which could be housed in the same housing or in a plurality of different housings.
A computing device is generally representative of any device suitable for performing the processing and analysis techniques discussed within the present disclosure. In some embodiments, the computing device includes a processor, random access memory, and a storage medium. In that regard, in some particular instances the computing device is programmed to execute steps associated with the data acquisition and analysis described herein. Accordingly, it is understood that any step related to data acquisition, data processing, calculation of the FAM, instrument control, and/or other processing or control aspects of the present disclosure may be implemented by the computing device using corresponding instructions stored on or in a non-transitory computer readable medium accessible by the computing device. In some instances, the computing device is a console device. In some instances, the computing device is portable (e.g. handheld, on a rolling cart, etc.). Further, it is understood that in some instances the computing device comprises a plurality of computing devices. In that regard, it is particularly understood that the different processing and/or control aspects of the present disclosure may be implemented separately or within predefined groupings using a plurality of computing devices. Any divisions and/or combinations of the processing and/or control aspects described herein across multiple computing devices are within the scope of the present disclosure.
It is understood that any communication pathway between the catheter and the computing device may be utilized, including physical connections (including electrical, optical, and/or fluid connections), wireless connections, and/or combinations thereof. In that regard, it is understood that the connection is wireless in some instances. In some instances, the connection a communication link over a network (e.g. intranet, internet, telecommunications network, and/or other network). In that regard, it is understood that the computing device is positioned remote from an operating area where the catheter is being used in some instances. Options for the connection include a connection over a network which can facilitate communication between the catheter and the remote computing device regardless of whether the computing device is in an adjacent room, an adjacent building, or in a different state/country. Further, it is understood that the communication pathway between the catheter and the computing device is a secure connection in some instances. Further still, it is understood that, in some instances, the data communicated over one or more portions of the communication pathway between the catheter and the computing device is encrypted.
The present invention provides a computer device and a computer-implemented method for the quantification of the extension of functional coronary artery disease (CAD) in a mammal, such as a human patient. It is clear that the extension of the functional CAD corresponds to the functional length of the CAD. The method determines the mismatch in the extent of CAD between anatomical and physiological invasive evaluations based on angiography, intravascular imaging and intracoronary hyperemic pressure tracing pullbacks. In particular, the extent of functional disease derived from FFR data (such as FFR data derived from pullback curves) can be quantified using a specially developed algorithm provided herein. The clinical relevance of the methods provided is that the mismatch between the length of anatomical and functional CAD (i.e. abbreviated as FAM, either derived from QCA or OCT) predicts improvement in epicardial conductance after percutaneous revascularization. It is clear that the length or extent of the anatomical CAD is determined by means of a detection of a particular part of the vessel comprising a reduction of the diameter, or the lumen area of the coronary vessel, or any other suitable indicator of an anatomical diversion of the vessel which for example exceeds a predetermined threshold. It is clear that the length or extent of the anatomical CAD is correlated to the part of the vessel, which can be considered as anatomically diseased as its anatomy impacts the blood flow along the coronary vessel negatively. It is clear that the length or extent of the functional CAD is determined by means of a detection of particular linearized segments of the vessel correlating to a reduction in the FFR values, or any other suitable indicator of a pressure change, which exceeds a predetermined threshold, and thereby determines the extent or length of the vessel which can be considered as functionally diseased based on the fact that, irrespective of detectable anatomical indicators, the functionality of the coronary artery in these segments is negatively affected.
QCA is based on conventional angiography and identifies CAD length as the extent of the stenotic segment. It is clear that, as described above, this refers to the extent or length of the anatomical CAD, which according to this embodiment is determined by means of QCA. On the other hand, OCT, possessing higher spatial resolution, derives lesion length from the selection of proximal and distal reference cross-sections without atherosclerotic plaques. It is clear that, this refers to the extent or length of the anatomical CAD, which according to this embodiment is determined by means of OCT. Therefore, it is expected that embodiments with CAD anatomical length derived from OCT will be equal or longer than embodiments with the QCA-derived length. We observed that the anatomical length of CAD was shorter when derived from QCA compared to the one derived from OCT; still, FFR pullbacks derived CAD length was longer (
Pressure pullbacks can show two distinct functional CAD endotypes, namely predominant focal or diffuse. In the focal functional CAD, pressure drops are commonly restricted to anatomical stenosis. In this disease endotype, PCI restores epicardial conductance, results in higher post-PCI FFR, increases the likelihood of relieving patients from angina and is associated with improved clinical outcomes. In contrast, in patients with functional diffuse disease, PCI results in minor improvement in vessel physiology, low post-PCI FFR and higher likelihood of persistent angina. Several novel methods are available to assess the pattern of CAD aiming at predicting the results of PCI in terms of coronary physiology. The pullback pressure gradient (PPG) index (Coroventis Research, Uppsala, Sweden), instant wave-free ratio (iFR) co-registration system (Philips, Best, the Netherlands) and the FFRCT revascularisation planner (HeartFlow Inc, Redwood city, USA) are novel approaches that may further personalize clinical decision making and refine revascularization strategies in patients with chronic coronary syndromes. In the present invention, we developed a complementary approach to predict the response to PCI by quantifying the extent of functional CAD from FFR pullbacks. The larger the functional length of CAD, the lower the functional gain obtained with PCI and the lower the likelihood of functional revascularization. This approach is analogous to the PPG where millimeters with functional disease are calculated based on an FFR threshold. In other words, according to this embodiment, there is defined a threshold of for example an FFR drop ≥0.0015/mm for labeling the parts of the coronary vessel exhibiting FFR deterioration. It is clear that, according to this embodiment, similarly this threshold, defines the parts of the coronary vessel which do not exhibit FFR deterioration. In the context of the PPG index, as referred to above, the length or extent of the functional disease was derived from a pullback curve, by for example aggregating the length of all parts of the curve where the FFR drop, or in other words the FFR reduction was ≥0.0015/mm. According to the embodiment of the current approach, as described above, first the set of FFR values, for example representing an FFR pullback curve, are transformed by means of the piece-wise linearization by applying an automated change-points detection algorithm into a sequency of segments, such as for example healthy segments and diseased segments, for example comprising focal diseased segments, diffused diseased segments, or any other suitable diseased segments. With the current approach, the length of functional disease is computed based on an automated algorithm classifying the FFR curve segments as healthy or diseased. In other words, after processing the set of FFR values, representing for example an FFR pullback curve, the automated change-points detection algorithm, converts the set of FFR values, by means of piece-wise linearization, into a sequence of linear segments, which are classified as healthy segments or diseased segments. According to a preferred embodiment, such a segment is classified as a healthy segment, when the segment does not exhibit FFR deterioration, or in other words, when for example the FFR drop and segment length of the segment define a position in a coordinate system, in which the FFR drop is the y-axis and in which the segment length is the x-axis, which is above a predetermined first classification threshold function. It is clear that the FFR drop is defined as the difference between the FFR values at the distal and at the proximal point of the segment. According to the exemplary embodiment shown in
It is also evident that the information obtained regarding characteristics of the coronary artery disease, such as for example the measurement of the length of the anatomical and functional region, can be considered in addition to other representations of the lesion or stenosis and/or the vessel, such as e.g. IVUS, for example including virtual histology, OCT, ICE, Thermal, Infrared, flow, Doppler flow, and/or other vessel data-gathering modalities, to provide a more complete and/or accurate understanding of the vessel characteristics. For example, in some instances the information regarding characteristics of the lesion or stenosis and/or the vessel as obtained by the system of the invention are utilized to confirm information calculated or determined using one or more other vessel data-gathering modalities.
Finally, it is to be understood that although particular embodiments, specific configurations as well as materials and/or molecules, have been discussed herein for methods according to the present invention, various changes or modifications in form and detail may be made without departing from the scope of this invention. The following examples are provided to better illustrate particular embodiments, and they should not be considered limiting the application. The application is limited only by the claims.
Clinical characteristics of patients are shown in Table 1. Overall, 117 patients (131 vessels) were included: 71 patients (81 vessels) in the derivation cohort and 48 patients (50 vessels) in the validation cohorts. QCA and OCT anatomical lesion lengths were available for all (n=50) and part (n=36) of the validation cohort, respectively (Table 2). FFR motorized pullbacks pre and post PCI were available in all cases.
Anatomical, functional, and procedural characteristics of the derivation and validation cohort are presented in Table 2. From the FFR curves, 431 segments were extracted. In detail, 151 (observer 1) and 156 (observer 2) segments were visually assessed as healthy, 101 (observer 1) and 106 (observer 2) as focal disease, 146 (observer 1) and 147 (observer 2) as diffuse disease (
PCI was performed in 50 vessels included in the validation cohort. Pre-PCI FFR was 0.74 [0.67-0.77] and diameters stenosis was 53.0% [47.25-59.50]. Anatomical CAD length derived from QCA was 16.05 mm [11.40-22.05], anatomical CAD length derived from OCT was 28.0 mm [16.63-38.0] and functional CAD length was 67.12 mm [25.38-91.37] (p<0.001). No correlation emerged between the extent or length of anatomical CAD derived from QCA and the extent or length of functional CAD derived from FFR pullbacks (r=0.124, 95% CI 0.168 to 0.396, p=0.390,
Mean stent length was 27.45±11.52 mm. Mean post-PCI FFR was 0.86 [0.82-0.89]. An explanatory example visualizing vessels with positive and negative FAM that underwent PCI and post-PCI FFR measurement is presented in
Patients in whom functional disease was confined within the anatomical lesion (i.e. FAM≥0) had the strongest improvement in relative functional gain (FAMQCA≥0.701±0.235 vs. FAMQCA<0 0.441±0.225, p<0.001). FAM either derived from QCA or OCT predicted functional gain (FAMQCA AUC 0.84, 95% CI 0.71 to 0.93, p<0.001 and FAMOCT AUC 1.00, 95% CI 0.93 to 1.00, p<0.001). The best FAMQCA and FAMOCT cutoff values predicting 50% gain in epicardial conductance were −57.64 mm and −37.19 mm, respectively. Percent FFR drops within the anatomical lesions either derived from QCA or OCT were strongly correlated with functional gain (r=0.792, 95% CI 0.655 to 0.879, p<0.001 for QCA and r=0.789, 95% CI 0.615 to 0.890, p<0.001 for OCT,
Tables
This is a multicenter, prospective registry of patients undergoing clinically indicated coronary angiography in whom motorized FFR pullback evaluations were performed before PCI. Patients presenting with acute coronary syndromes, previous coronary artery bypass grafting, significant valvular disease, severe obstructive pulmonary disease or bronchial asthma, coronary ostial lesions, severe tortuosity, or severe calcification were excluded. Patients with adequate pressure tracings and pullback curves were included in this analysis. The study was approved by the Ethics Committee at each participating center. The study population is a combination of two prospective studies NCT03824600 and NCT03782688.
Angiographies were performed using a dedicated acquisition protocol. Two angiographic projections separated at least 30 degrees were obtained for each target lesion after the administration of intracoronary nitrates (
Examinations were performed using the OPTIS™ OCT systems (Abbott Vascular). OCT pullback at 36 mm/s were acquired before pre-dilation if feasible. OCT-derived anatomical lesion length was defined as the distance between the proximal and distal reference segments using the OCT automated lumen detection feature. Stent diameter selection was based on the distal reference mean external elastic lamina (EEL)-based diameters rounded down to the nearest available stent size (usually in 0.25 mm increments) to determine stent diameter. If the EEL could not be adequately visualized, the stent diameter is chosen using the mean lumen diameter at the distal reference rounded up to the next stent size. Optimization of the device for performed based on OCT at operator discretion.
Fractional flow reserve (FFR) measurements were performed with the Pressure Wire X (Abbott Vascular, Chicago, Il, USA) that was connected to a motorized pullback device at a speed of 1 mm/s (R 100, Philips Volcano, San Diego, Ca, USA). Pressure pullback measurements were acquired at a sampling frequency of 100 Hz. A continuous intravenous adenosine infusion was given at a dose of 140 mg/kg/min via a peripheral or central vein to obtain steady-state hyperemia for at least 2 min. The position of the pressure sensor was recorded with a contrast injection to identify the pullback initial position for co-registration purposes. In cases undergoing PCI, FFR measurements were repeated at the same anatomical location. FFR gain was defined as FFR post-minus FFR pre-PCI. If FFR drift (>0.03) was observed, the FFR pullback was repeated. For the FFR gain a FFR post and pre-PCI were determined as the ratio of the distal pressure or Pd at the most distal part of the coronary artery, with respect to the proximal or aortic pressure Pa at the ostium of the coronary artery.
The FFR curve along the vessel axis was reconstructed by applying a moving average filter with a window size of 10 s, followed by an infinite impulse response low pass elliptic filter (0.1 Hz cutoff frequency) for smoothing (
Such an embodiment is especially useful for quantifying the extent and/or patterns of coronary artery functional disease in a coronary vessel from a patient under hyperaemic conditions, wherein the patient is a mammal, for example a human. Under such conditions there can be generated pressure values that represent an FFR pullback curve during a pullback operation. Determining such an FFR pullback curve is done by determining FFR values from measurements of the movable pressure sensor, also referred to as distal pressure or Pd, with respect to a stationary pressure sensor, also referred to as proximal or aortic pressure Pa, during the pullback time period. It is clear that for example the stationary pressure sensor is positioned at the ostium of the vessel and that the movable pressure sensor, during the pullback time period is moved between a more distal part of the vessel, for example the most distal part of the vessel, or a part of the vessel distal of a suspected stenosis, stricture or lesion, and the ostium of the vessel. When such measurements are for example performed under hyperaemic conditions in a coronary artery, then these values that are determined based on these measurements during the pullback time period, during which the movable sensor is moved along the vessel, are referred to as FFR values of an FFR pullback curve, and are typically determined as the ratio of Pd/Pa, wherein Pd and Pa could for example be determined from the measured pressure values after any suitable form of pre-processing such as for example by means of a moving mean or average function which is configured to filter out the rhythmic and/or periodical component of the heartbeat cycle. It is thus clear that according to a preferred embodiment FFR could for example be defined as the ratio of mean or average distal coronary pressure, which is the pressure measured by the movable pressure sensor, and the mean or average aortic pressure, which for example is the pressure measured by the stationary pressure sensor, measured during, preferably maximal, hyperaemia that is preferably achieved through administration of a potent vasodilator such as for example adenosine, ATP or papaverine either by IV infusion or by intracoronary (IC) bolus injection. Under such conditions of hyperemia, by rendering myocardial microvascular resistance constant and minimal, the impact of disease in the epicardial conduit artery on myocardial blood flow is more advantageously separated out.
An automatic algorithm was developed for functional length quantification from FFR curves. The first step of the algorithm consisted in the piece-wise linearization of each FFR curve, see for example
Two cohorts were defined to develop and validate the part of the algorithm performing automatic FFR segments classification. The derivation cohort consisted of patients with CAD defined as distal FFR<0.90. For this cohort, only baseline (i.e. pre-PCI) FFR pullbacks were included. These were selected in a consecutive fashion from all patients included in the registry. The validation cohort included subsequent patients with CAD defined as a distal FFR 0.80 who underwent OCT-guided PCI and FFR measurement after stent implantation.
Two independent observers (observer 1: CaC; observer 2: SN) preliminarily adjudicated by visual inspection each one of the piece-wise linearized FFR curve segments belonging to the two cohorts as ‘healthy’, i.e. without FFR deterioration, or as ‘diseased’. Then, the two observers performed a further adjudication on ‘diseased’ segments, discriminating between ‘focal’ or ‘diffuse’, based on the presence of step-ups in the FFR pullback linearized curve.
The visual adjudication of the derivation cohort was used to develop the automatic classifier, based on a two-variables logistic regression. The two independent variables considered for the logistic regression were the length of the linearized segment and the associated FFR drop.
The detection of main changes in the distributional properties of FFR pullback curves was here addressed implementing a change points identification strategy. The implemented approach leads to a piece-wise linearization of FFR pullback curves based on a change points detection problem, where a change point is defined as a sample of the acquired FFR pullback curve at which an attribute of the curve suddenly changes. It is clear that a sample of the acquired FFR pullback curve, corresponds to a particular position along the part of the coronary vessel where the corresponding FFR pullback curve was generated. It is clear that a sudden change of the attribute of the FFR pullback curve, corresponds to an identifiable change, at that particular position in a relevant attribute of the FFR pullback curve, such as for example described in further detail below. As for example described below, the attribute of the FFR pullback curve, could for example be the average value and/or the slope along segments of the FFR pullback curve, or in other words subsets of the set of FFR values obtained at different positions of a coronary vessel between the ostium and the most distal part of the coronary vessel. Technically, a parametric global method detailed in Killick, R. et al (2012) J. Am. Stat. Assoc. 107, 1590-1598 and Lavielle M. (2005) Signal Processing 85, 1501-1510 were implemented here in MATLAB environment (MathWorks, Natick, MA, US) for FFR pullback change points identification. The steps of the implemented algorithm leading to a single change point detection are the following:
It is thus clear that the change point is configured to divide the FFR pullback curve in two segments, in which the change point defines the endpoint between these two segments.
The problem expressed by points 1-5 can be translated into an algorithm as explained in the followings. Given a generic FFR pullback curve, FFR=(FFR1, FFR2, . . . , FFRN), where FFRi is the FFR value at i-th sample of the curve and N the total number of samples, the problem consists in finding the k-th sample minimizing the cost function
J(k)=Σi=1k−1Δ(FFRi;χ([FFR1 . . . FFRk−1]))+Σi=kNΔ(FFRi;χ([FFRk . . . FFRN])) (1),
where χ is the empirical estimation of the statistical property of interest and Δ is the deviation measure. Since we are interested in highlighting changes in average value and slope along the FFR pullback curve, here a linear function was adopted as statistical property of interest. This is like to say that for a generic interval between points m and n along the FFR pullback curve (
A generic FFR pullback curve might have several change points, the number of change points being unknown a priori. Since adding change points decreases the residual error, the overfitting of the FFR curve is avoided by adding a penalty term which is a linear function of the number of change points to the cost function, which can be expressed as:
where kr and kC are the first and the last sample of the FFR pullback curve, respectively, C is the number of change points, and B is the fixed penalty term (set equal to 0.1 in this study). The minimization of the cost function was obtained implementing an algorithm based on dynamic programming with early abandonment Killick, R. et al (2012) J. Am. Stat. Assoc. 107, 1590-1598.
It is thus clear that the one or more change points are configured to divide the FFR pullback curve, in two or more segments, in which the change points define the endpoint between two segments, or in other words between two neighboring segments. It is further clear that said segments according to the embodiment shown correspond to a linear function between the bordering change points, or in other words correspond to linearized segments. It is clear the FFR pullback curve corresponds to a set of fractional flow reserve (FFR) values obtained at different positions of a coronary vessel between the ostium and the most distal part of the coronary vessel. It is thus clear that each segment extends between a proximal point of the segment and a distal point of the segment, and corresponds to a subset comprising the FFR values obtained at different positions between said proximal point of the segment and the distal point of the segment, of the set of FFR values obtained at different positions of a coronary vessel between the ostium and the most distal part of the coronary vessel. In other words, each linearized segment corresponds to a linearized subset of the FFR pullback curve, or in other words the corresponding set of FFR values, bordered by at least one change point. It is clear that the proximal point of the segment corresponds to a position closer to the ostium of the vessel of the set of FFR values. It is clear that the distal point of the segment corresponds to a position closer to the distal part of the vessel of the set of FFR values. It is further clear that the two or more segments, according to the embodiment shown, extend between the ostium and the most distal part of the coronary vessel in such a way that a first segment extends between the ostium and a first change point, and a last segment extends between the last change point and the most distal part of the coronary vessel. It is clear that in an embodiment in which there are two or more change points detected in the FFR pullback curve, there will be a corresponding sequence of one or more segments which extends between this first segment and the last segment.
Once the piece-wise linearization of FFR pullback curve has been carried out, each linearized segment of the curve was then characterized by two quantities (
On each segment a third quantity, the segment slope, was defined as the ratio between FFR drop and segment length.
Two independent observers (CaC, SN) adjudicated by visual inspection each segment of the piece-wise linearized of FFR pullback curve as belonging to one of five segment phenotypes: (1) healthy segment; (2) focal disease segment; (3) diffuse disease segment; (4) pressure recovery segment; (5) artifact.
The automatic adjudication of each segment of the piece-wise linearized of FFR pullback curves to one of the defined five segment phenotypes was based on the application of a logistic regression model on FFR pullback curves segments from the derivation cohort (81 Patients). Technically, the logistic regression model was built on segments that obtained the same adjudication from the two observers (216 segments), while segments manually classified as artefacts were not considered. Two independent variables were considered in the logistic regression model to discriminate the segments: FFR drop and segment length. It is clear that in this way suitable first and second classification threshold functions were derived for the logistic regression model to discriminate healthy segments from diseased segments, and focal diseased segments from diffuse diseased segments as mentioned above. It is clear that the first and second classification threshold functions are different. It is further clear that the first classification threshold function, according to the embodiment shown, operates on all the segments. It is further clear that the second classification threshold function, according to the embodiment shown, only operates on the segments that are classified as diseased by means of the first classification function. It is clear that any other suitable classification threshold functions might be derived, from any suitable data set, such as a suitable derivation cohort, comprising any suitable number of patients, from which any suitable number of observers, provides a suitable classification of the segments. According to an embodiment, the classification threshold functions, such as for example for classifying the segments by means of the logistic regression model, can be derived from any suitable supervised, data-driven approach. It is clear that according to alternative embodiments any suitable training dataset, comprising any suitable size could be used to determine such suitable classification threshold functions for a suitable model for automatically classifying the segments, for example based on FFR drop, segment length, slope or any other suitable parameter of the segments.
Since the logistic regression provides a binary separation, a two-steps approach was implemented for the piece-wise linearized FFR pullback curve segments automatic adjudication:
In the healthy segment phenotype, pressure recovery segments were identified as the ones meeting all the following criteria: they have a positive FFR drop, they are contiguous to a diseased segment, and they are shorter than 20 mm.
The performance of the automatic adjudication was evaluated by comparison with the manual adjudication by the two observers (CaC; SN) on piece-wise linearized FFR pullback curves belonging to the validation cohort (50 patients, 179 segments). The results of the automatic adjudication are reported in
The ability of the proposed automatic adjudication method in discriminating healthy, focal, and diffuse disease segments clearly emerges (
Sensitivity of FAMQCA and FAMOCT to serial lesions was tested with two approaches: (1) considering only the contiguous segments in the definition of the functional length and (2) excluding serial lesions from the analysis. When considering only the contiguous segments to define the functional length, both FAMQCA and FAMOCT were still correlated with FFR relative gain after PCI (r=0.606, 95% CI 0.387 to 0.760, p<0.001,
The automatic piece-wise linearization and classification of the FFR curve segments allowed to derive the extent of the functional disease, namely the functional length, which was expressed in millimeters. In detail, the functional length of disease for each coronary artery was obtained as the summation of the length of all linearized FFR curve segments classified as diseased by the algorithm. In the presence of serial or multiple lesions (i.e. functionally diseased segments separated by functionally healthy segments), the functional length was considered as the sum of all (i.e., contiguous and non-contiguous) diseased segments.
The difference between the anatomical and the functional length of CAD was defined as Functional Anatomical Mismatch (FAM) (
PCI was performed following standard of care guided by FFR and OCT, both executed before and after stent implantation. Intraprocedural PCI guidance or stent optimizations based on either physiology or imaging were left at operator's discretion. New generation DES were used in all cases. To quantify the impact of PCI, the relative functional gain was defined as post PCI FFR minus pre-PCI FFR divided by 1−pre-PCI FFR.
Continuous data are presented as mean (±SD) or median [25th-75th percentiles]. Categorical data are presented as counts and proportions (%). Differences were evaluated using the univariate Mann-Whitney non-parametric U test. Spearman's correlation coefficients were calculated to assess the relationship between FAM and post-PCI FFR. Agreement between observers was assessed by the intraclass correlation coefficient (ICC). The optimal cutoff values of FAM to predict relative functional gain were calculated using receiver operating characteristic (ROC) curves. The discriminant ability of FAM value to predict optimal post-PCI physiologic results was evaluated with area under curve (AUC). Optimal relative functional gain was defined as an increase in epicardial conductance greater than 50%.
Number | Date | Country | Kind |
---|---|---|---|
2020468.1 | Dec 2020 | GB | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2021/087477 | 12/23/2021 | WO |
Number | Date | Country | |
---|---|---|---|
20240130674 A1 | Apr 2024 | US |