Method and algorithm for defining the pathologic state from a plurality of intrinsically and extrinsically derived signals
The present invention describes methods and algorithms for processing a plurality of clinically relevant signals/data/genotypes intrinsic to a given patient and/or signals/data derived from external diagnostic equipment. The intrinsic signals can be acquired from device-based sensors and analogous extrinsic data can be obtained from imaging equipment. These signals/data are input into software algorithms that use digital signal processing in order to output informational data sets of clinical and technical relevance after comparisons are made to patients with access to this technology whose outcome under varying treatments is known. This data is used to define prognosis, guide clinical decision-making, make treatment suggestions, and direct programming of cardiac devices. Evaluation of an implanted device using such technology will confirm appropriately programmed settings within such a device. Application of these methodologies and algorithms to analyze signals obtained by an implanted device allows for closed loop programming of such a device based on multiple parameters using digital signal processing. These technologies will also allow for open connectivity between external imaging equipment and implanted devices as to derive a plethora of diagnostic and prognostic information. Comparisons of large data sets of intrinsically and extrinsically derived parameters are used to develop a translation function that relates device based indices to analogous, conventional indices that are currently state of the art. Multiple indices reflecting a plurality of diagnostic information are input into such software algorithms (intrinsically and/or extrinsically derived) in order to generate prognostic information and present therapeutic options based on evidence based medicine.
This patent application claims priority to provisional patent applications No. 60/647,102 filed Jan. 26, 2005; and 60/660,101 filed Mar. 9, 2005, incorporated herein by reference.
BACKGROUND OF THE INVENTION1. Field of the Invention
This application pertains to a system and apparatus for determining the status of a patient and more particularly to a system and apparatus wherein measurements from the patient are combined with external information such as statistical information collected from other patients and results of diagnostic testing to obtain a diagnosis, treatment options and prognostic information for the patient quickly and accurately using evidence based medicine.
2. Description of the Prior Art
Major advances are occurring in the development of imaging modalities and implantable technologies capable of diagnosing and treating a variety of pathophysiology. Thus, the clinician has numerous diagnostic tests at his or her disposal and an option of therapeutic regimens. Coupled with the wealth information from these tests, there is a need to expeditiously evaluate the results of a plethora of diagnostic tests, and compare the data to historical data sets. Treatment algorithms designed to assess such informational data sets provide better direct therapy (evidence-based medicine).
To a large extent, practice decisions based on anecdotal data and individual studies have dictated medical practice to date. However, large scale population studies and the development of registries along with digitization of acquired diagnostic data obtained from these studies provide more powerful statistical analyses of patient outcome by comparing differing therapeutic modalities. The availability of such data will depend on a means for device-device communication and the evolution of wide range digitization of medical records and of data derived from different diagnostic equipment. Though the majority of the algorithms and examples described herein are in reference to advances in cardiovascular medicine and genetics, the inventions described have broad range application to any medical field.
The advent of digital signal processing (DSP) in evolving technologies (e.g. implantable pacemaker/defibrillators) will allow for informational data sets to be available in discrete numeric format, processed by high-speed microprocessors, and incorporated into device-based software algorithms. Such processing will allow a user to relate, in a complementary fashion, clinically relevant data obtained from both an imaging apparatus and implanted device as to derive a composite of information for diagnostic purposes and for optimizing patient management. The application of DSP will enable calculation algorithms to perform several processing operations simultaneously.
By way of example, software algorithms for readily evaluating a number of variables descriptive of cardiac performance and electromechanical dysynchrony (physiological properties) are assimilated within an implanted device or downloaded between an implanted device and a separate apparatus/extrinsic diagnostic equipment (e.g. echocardiography machine, cardiac MRI). A composite of physiological properties obtained from the implanted device is digitized (if necessary) and compared to normal and pathologic values as to generate physiological descriptors that have a numerical score. Models constructed to predict probability of outcome from the data obtained are implemented for such comparisons (Selker et al. Patient specific predictions of outcomes in myocardial infarction for real-time emergency use: a thrombolytic predictive instrument. Ann Intern med 1997; 127: 538-56). These physiologic descriptors are then input into software algorithms as to produce an informational data set of clinical and technical relevance. This informational data set (IDS) is output in form of an easily interpretable set of recommendations or prognostic data for the clinician. In one embodiment, such IDS can be downloaded into removable digital storage media or other media compatible with implanted device software and incorporated into an electronic medical record (EMR). In a preferred embodiment, the IDS is available to closed loop control systems that direct device based therapies (e.g. Cardiac Resynchronization Therapy, CRT).
The following references provide background information for the present application and illustrate the state of the art. All these references are incorporated by reference.
U.S. Pat. Nos. 6,804,559, 6,795,732, 6,792,308, 6,816,301, 6,572,560, 6,070,100, 6,725,091, 6,628,988, 6,740,033, 5,971,931, 5,833,623, 6,826,509, 6,805,667, 6,574,511, 6,418,346, 5,549,650
Published US patent applications: 20040176810, 20030083702, 20020026103, 20040111127, 20020072784, 20030216620, 20040167587, 20050182447, 20050043895
REFERENCES IN PEER-REVIEWED JOURNALS
- Meluzin Jaroslav, Novak Miroslav, Mullerova Jolana, Krejci Jan, Hude Petr, Eisenberger Martin, Dusek Ladislav, Dvorak, Ivo, Spinarova Lenka, A fast and simple echocardiographic method of determination of the optimal atrioventricular delay in patients after biventricular stimulation. PACE, 2004. 27: p. 58-64.
- Perego Giovanni B, Chianca Roberto, Facchini Mario, Frattola Alessandra, Balla Eva, Zucchi Stefania, Cavaglia Sergio, Vicini Ilaria, Negretto Marco, Osculati Giuseppe, Simultaneous vs. sequential biventricular pacing in dilated cardiomyopathy: an acute hemodynamic study. The European Journal of Heart Failure, 2003. 5: p. 305-313.
- Ritter P, Padeletti L, Gillio-Meina L, Gaggini G., Determination of the optimal atrioventricular delay in DDD pacing. Europace, 1999. 1: p. 126-130.
- Van Gelder Berry M, Bracke Frank A, Meijer Albert, Lakerveld Lex J M, Pijis Nico H J, Effect of optimizing the VV interval on left ventricular contractility in cardiac resynchronization therapy. Am J Cardiol, 2004. 93: p. 1500-1503.
- Yu C M, Lin H, Zhang Q., High prevalence of left ventricular systolic and diastolic asynchrony in patients with congestive heart failure and normal QRS duration. Heart, 2003. 89: p. 54-60.
- Pappone C, Augello G, Rosanio S, et al. First Human Chronic Experience with Cardiac Contractility Modulation by Nonexcitatory Electrical Currents for Treating Systolic Heart Failure: Mid-Term Safety and Efficacy Results from a Multicenter Study. J Cardiovasc Electrophysiol, 2004; 15, 418-422.
- Dipla K, Mattiello J A, Margulies K B et al. The sarcoplasmic reticulum and the sodium/calcium exchanger both contribute to the calcium transient of failing human ventricular myocytes. Circ Res 1999; 84: 435-444.
- Burkoff D, Shemer I, Feizen B, et al. Electric currents applied during the refractory period can modulate cardiac contractility in vitro and in vivo. Heart Fail Rev 2001; 6: 27-34.
- Pappone C, Rosanio S, Burkoff D, et al. Cardiac Contractility Modulation by Electric Currents Applied During the Refractory Period in Patients with Heart Failure Secondary to Ischemic or Idiopathic Dilated Cardiomyopathy. Am J Cardiol 2002; 90: 1307-1313.
- Willems R, Sipido K R. Nonexcitatory Stimulation as a novel treatment for heart failure: cause for excitement? European Heart Journal 2004. 25: 626-628.
- Padeletti L, Barold S S. Digital Technology for Cardiac Pacing. Am J Cardiolo 2005; 95: 479-482.
- Thomas J D, Greenberg N L, Garcia M J. Digital echocardiography 2002: now is the time. J Am Soc Echocardiograpy 2002; 15: 831-8.
- Feignebaum H. Digital echocardiography [review]. Am J Cardiology 2000; 86: 2G-3G.
The present invention relates to acquiring a plurality of diagnostic information based on intrinsic properties of a given patient and extrinsic information derived from diagnostic testing performed on the patient, and combining the available data to derive prognostic information about a specific pathologic state using evidence based medicine. Recommendations for therapy delivered intrinsically via an implanted device or extrinsically via various therapeutic modalities are made with such analysis algorithms. Some of the examples used herein are focused mainly on implantable cardiac rhythm management (CRM) devices and other similar devices for diagnosing or treating heart failure. These types of devices generate a multitude of measurements diagnostic of a specific pathologic state, based on software algorithms incorporated into device-based platforms to obtain a composite of relevant physiological data from implanted sensors/transducers. The algorithms generate informational data sets for diagnostic/monitoring purposes with the object of guiding physician management and/or programming of an implanted device via a closed loop control system. Such informational data sets may be used, for example, to guide titration of pharmaceutical therapies, diagnose myocardial ischemia or determine a patient's candidacy for coronary revascularization or valve replacement surgery. This technology is also capable of incorporating downloadable indices/data from extrinsic diagnostic equipment into implanted devices/programmers and vice versa at periodic intervals via removable digital media (e.g. removable hard drive, magnetic-optical disc) or wireless telemetry (e.g. Bluetooth). Bi-directional communication of this data is used to confirm adequate functioning of an implanted device and verify diagnoses made with extrinsic equipment (cross-verification). This data can be examined to confirm response to specific therapeutic modalities such as cardiac resynchronization therapy (CRT), cardiac contractility modulation (CCM) or alternate Non-Pharmacologic Inotropic Therapy (NPIT).
In an alternate application, these algorithms are used in the field of genetic medicine through a similar means, combining data reflective of properties intrinsic to an individual in conjunction with environmentally based or extrinsic characteristics of that individual in context of known predictors of a given disease state derived from large population studies.
Diagnostic imaging modalities are moving toward digitization with standardization of storage format adherent to standard models (Digital Imaging and Communication in Medicine). Thus, there is a need to formulate composite indices which can be digitally processed and input into fast software algorithms for processing. Extrinsically derived composite indices can be evaluated from time to time and compared to similar indices generated from an implanted device (intrinsic), assimilated into an implanted device's existing diagnostic data. Translation of mathematical indices derived from signals acquired by implanted device sensors into conventional indices commonly used with external imaging modalities will facilitate device-device communication.
Through the utilization of this invention, a means of translation between analogous intrinsic/extrinsic indices is developed from the acquired pooled data after significant numbers of patients gain access to these technologies. The concepts underlying this invention are extended to imaging technologies including, but not limited to, echocardiography, magnetic resonance imaging, PET scans, nuclear imaging or computed tomography, and other diagnostic/laboratory tests. The extrinsically derived composite data can be downloaded into an implanted device for storage (e.g. medical record keeping) and available to the clinician in combination with similarly derived device based indices. The combined informational data set (intrinsic and extrinsic) can then more accurately provide prognostic information and generate recommendations for various therapies (e.g. using neural networks). Comparisons between extrinsic and intrinsic diagnostic data can be used to confirm diagnoses (e.g. presence of myocardial ischemia). A method and means for correlating intrinsically and extrinsically derived data with a translation function is developed once large numbers of patients have access to these technologies and patient outcome under varying clinical circumstances is determined. Digitization and standardization of storage formats will facilitate the application of such a translation function for derivation of analogous intrinsic and extrinsic indices using a universal mathematical language that will provide the clinician with prognostic information and treatment suggestions.
By way of example, this technology is generally described in relation to its application to implanted devices in patients with congestive heart failure as this particular field is in a state of rapid technologic evolution. Such implanted devices can be expected to have lead based and non-lead based sensors and transducers capable of acquiring data that is analogous to the diagnostic information obtained from extrinsic imaging modalities. A brief review and examples of such technologies are included herein along with related references. This review serves to familiarize the reader with the mechanical and physiological principles helpful for understanding this invention by correlating device based sensor data to analogous data acquired with echocardiography.
Patients referred for consideration of CRT implantation undergo echocardiographic or other imaging analysis. Current ultrasound technologies allow for a means of evaluating changes in regional volume, differential myocardial motion and contractility. Temporal frame rates for evaluating such subtle differences in timing are currently under 10 milliseconds. Available equipment manufactured by a number of companies such as General Electric, Philips, TomTec and Siemens are capable of measuring changes in both global and regional volumes within the cardiac chambers during the cardiac cycle. Analysis of tissue Doppler data allows for calculation of tissue velocities, myocardial strain and strain rate derived from the spatial gradient of tissue velocity (strain rate equation). This can also be performed by defining relative locations of ultrasound reflectors in two or three dimensional space over time using a technique referred to as speckle tracking (GE). Such data is currently processed at fast enough rates as to generate real time parametric imaging. As improvements in processing and microprocessor robustness occur (Moore's law), one can expect a number of analyses (unique and redundant) will be performed simultaneously in digital format and be able to be compiled into informational data sets that can be used for diagnostic purposes. This also holds true for implantable CRM and heart failure devices.
Recent advances in cardiac MRI have allowed temporal frame rates that approach the level of ultrasound imaging. These frame rates are adequate enough to allow for regional strain analyses that define degree and location of electromechanical dysynchrony during two cardiac cycles. Magnetic resonance tagging, phase contrast MRI, and stimulated MRI (DENSE techniques) are MR methodologies capable of characterizing myocardial motion and strain. Newer techniques such as use of harmonic phase analysis (HARP) vastly improve upon imaging and processing time and allow for more rapid strain analysis. In the Fourier domain, spatial modulation of magnetization (SPAMM) tagged MR data will have multiple spectral peaks. Any off-origin spectral peak will relate to tissue motion. Spatial derivatives of harmonic phase can be computed for any pixels as to derive strain measurements. Such techniques allow for acquisition of a number of cardiac performance parameters and measurements of electromechanical timing (dysynchrony). The details of such innovations are described in the cardiac MRI literature.
An index of dysynchrony or dysynchrony index, DI, and cardiac performance parameters, CP, may be derived by any and all of these imaging modalities as well as within an implanted device. The DIs and CPs allow for a means of quantifying and localizing electromechanical dysynchrony and quantifying systolic and diastolic cardiac performance. A number of DIs and CPs have been defined in the literature and these, as well as, several novel DIs and CPs developed by the inventor (see below and cited references) may be used for such analysis. The user can implement fast software algorithms using DSP to compile informational data sets (IDS) descriptive of cardiac performance (CP) and electromechanical dysynchrony.
The present invention provides for a rapid interpretation of physiological properties (DI, CP) as to generate relevant information, IDS, to the physician who is considering delivering a specific therapy to the patient (e.g. programming a CRT device, altering medical therapies, ordering diagnostic tests). Any available index or index to be defined may be incorporated into algorithms for performing such tasks and those mentioned in this application are by way of example. Demographic information and genotype can also be input into these algorithms.
Non-Pharmacologic Inotropic Therapy—Cardiac Contractility Modulation
CRT devices are well known in the field and to those experienced in the art. Non-excitatory stimulation, NES, is an evolving technology that holds promise to improve heart failure symptoms and may be less well known to the reader. Thus, a description of this technology and references are incorporated herein for the sake of completeness.
Devices capable of delivering non-excitatory stimulation to myocardium during absolute refractory periods (Impulse Dynamics) have been demonstrated to improve cardiac contractility by a mechanism termed cardiac contractility modulation (CCM). Initial experience with CCM devices suggests patients with congestive heart failure have improvements in symptoms and may be expected to have an improved outcome but such technology is expensive and the battery longevity of the device is often less than 6 months. These devices function by augmenting intracellular calcium concentration in regions of stimulated myocardium. This results in improvements in global cardiac performance for reasons that are currently not well understood. References that relate to this technology are cited below and can be reviewed if necessary. The inventor describes novel indices related to global tethering effects illustrative of how CCM affects cardiac performance and explains theories and potential mechanisms by which CCM improves global cardiac function with only local current delivery. Indices related to global tethering can be derived using this technology and evaluated by a physician to determine which patients may benefit from such expensive technology. The same indices can be monitored as part of a closed loop system within an implanted device as to confirm benefit and direct how often such therapy should be delivered as to reduce battery expenditure.
Application to Digital Devices—Interfacing/Telemetry
Currently implanted pacemakers/defibrillators which rely on analog technology will be evolving into devices that universally transform all analog signals into digital format. This will allow data compression and storage of extensive amounts of information in digital libraries within the device. Transfer of relevant clinical data (informational data set) between imaging equipment and implanted devices through digital technology can then be easily accomplished. This transfer may be by direct connection or wireless telemetry sent to implanted devices and device programmers. Security measures to ensure compliance with HIPAA regulations may require the use of fixed network addresses that place restrictions on access to such data. The clinical utility of the methods and algorithms described herein will be best realized when digital processing operations occurring within implanted devices and extrinsic diagnostic equipment can occur expeditiously and in a standardized format.
BRIEF DESCRIPTION OF THE DRAWINGS
By way of example, the algorithms described herein serve to guide programming of a CRT device and program the frequency of NPIT as part of a closed loop system, as well as, identify the presence of myocardial ischemia. The dysynchrony indices will be generally referred to as DI and any cardiac performance parameter as CP. These physiological properties are well known in the art as established in the above-referenced published U.S. patent applications. The examples described are included as to illustrate how the inventive algorithms function, though the application of such algorithms are by no means limited to these clinical scenarios. The DI's and CP's can be either intrinsically acquired or extrinsically derived (Dix, CPx) as denoted in
Referring to
Analysis of physiological properties obtained from large-scale clinical trials with a variety of imaging technologies (e.g. ultrasound, MRI) defines specific characteristics. For purposes of explanation, most of the explanations in this application will focus on properties that distinguish heart failure patients with varying forms of electromechanical dysynchrony and differences in cardiac performance and genetic markers for disease. A number of ongoing studies are examining these characteristics and how they relate to patient outcome and response to varying therapies.
In one application of the invention, these characteristics can be defined using monitoring technologies within implanted cardiac devices. These characteristics are compiled from numerous patients and evaluated as pooled data to correlate clinical outcome/prognosis based on various CPs, DIs, genotypes, relevant demographics and patient diagnostic information.
By way of example, these characteristics can determine how to program interval timing in a CRT. These characteristics can also direct the duration and frequency of activation of novel pacing modalities such as CCM or alternate NPIT. Characteristics derived from external monitoring equipment, (e.g. CPx and DIx), can be incorporated into reference/template data in digital format for comparison to real time acquired DIx, CPx, using predictor algorithms in the Analysis Apparatus, 200 along with intrinsically derived indices obtained from the implanted device itself (e.g. CP and DI). These physiological properties are assigned a numerical value along a graduated scale and output from 200 as physiological descriptors, 210.
In a preferred embodiment, the range of the numerical scale is based on similarly acquired data gathered from population studies, preferably using monitoring techniques capable of transmitting such information via wireless telemetry to a central processing center/data bank and back to the specific monitoring equipment (e.g. implanted device, echocardiography machine). The resulting numerically scored data, 210, is input to the Characterization Apparatus, 300, in
At step 110 the data is compiled, and if in analog form, is digitized in 100 and then input to 200 where the physiological descriptors are scored as to characterize the detected pathophysiological features of the patient using multivariate statistical techniques such as Discriminant Analysis. Other statistical techniques may be used to perform an analysis of these physiological descriptors. By way of example, Discriminant Analysis will compare measured DI. DIx and CP, CPx values related to regional and global systolic/diastolic function and attributes or assigns a binary numerical value to a range of characteristics based on comparisons to normal template data acquired from population studies. As an example, an arbitrary characteristic A may be defined as a function of three-impedance based measurements; variation in the integral under regionally derived dynamic impedance waveforms, Z (t) dt, during the systolic ejection phase, (Discriminant characteristic X), the standard deviation in peak regional impedance from multiple electrode combinations (Discriminant characteristic Y) and the standard deviation of peak positive dZ/dt during systole in multiple vectors (Discriminant characteristic Z). These characteristics are automatically obtained by device based sensors, accrued and digitized within the Detection Apparatus, 100. In one embodiment, an operator can have the option to review any automatically obtained characteristic (e.g. impedance waveform) before entry into 200, to confirm that an adequate signal is present for data processing. Analysis apparatus 200 assigns to characteristic A a specific value as described below. The values of A reflect a composite of findings that are weighted based on level of relevance to characteristic A and results of pooled data (i.e., known results of extrinsic diagnostic tests and intrinsic sensor based data) obtained from patients implanted with devices with this capability (this data will have been telemetered to the central processing center and then telemetered back to the implanted device). Thus, characteristic A may equal a1 when Discriminant characteristics X, Y and Z fall into a specific range of values. The relative importance of Discriminant characteristics X, Y and Z for Characteristic A is represented by values x, y and z. In this example x may be weighted at 50 percent (x=0.50), y and z at 25 percent, and
A=a1 when Discriminant value L<=xX+yY+zZ
By way of example, A equal to a2 or a3 is found when the morphology of regionally derived impedance waveforms (X) septal and lateral peak Z (Y) and regional comparisons of dZ/dt in more than one vector (Z) are consistent with dilated, non-coronary CM, and values in excess of a5 are found when values are consistent with multiple infarctions (more heterogeneous properties). Values of a4 and a5 are seen under circumstances consistent with a single infarct zone. The higher the number following a, the more likely the cardiomyopathy (CM) is a result of multiple regional infarctions rather than non-coronary CM (a1). Changes in any or all of these characteristics over time or at increased heart rates can indicate a new myocardial infarction or myocardial ischemia, respectively, and prompt a clinician to perform angiography.
Referring to
In this case example, C is equal to MRa (step 360 ensues) when Discriminant Analysis of echocardiographic data demonstrate that the MR flow exceeds a certain value. The direction of MR can be detected based on extrinisic echocardiographic imaging and noted to be eccentric rather than central. Such a finding is consistent with electromechanical dysynchrony and when taken in context of other variables at step 370, suggests that such a patient will respond to CRT therapy and LV pacing can be activated if available in the implanted device. C=MRx indicates that mitral valve surgery should be considered. Other complementary data (intrinsic or extrinsic) can be evaluated before such a recommendation is made. Some examples include; larger effective regurgitant orifice areas determined by ultrasonic color Doppler, and increases in LV end-diastolic dimensions or assessment of MR from implanted device based sensors.
Again the value of C is assigned in Analysis Apparatus 200. In 200 the degree of MR is given a numerical value based on comparisons of MR severity derived from normal and abnormal patients along a graduated scale. Thus, in this particular patient, echocardiographic assessment of MR determines that it is severe (e.g. EROA value near maximal) and the value of MR is graded numerically as 9 out of 10, where 10 is most severe and 0 being normal with no MR present). Peak endocardial acceleration, Acceleration index, at query box DV, has characteristics that are consistent with multiple myocardial segments reaching peak myocardial acceleration at significantly different times during the cardiac cycle (septal and lateral regions). In fact, the delays in peak myocardial acceleration are more pronounced and values of peak myocardial acceleration are lower at increased heart rates, suggesting that these territories have viable, dysynchronous myocardium. Such a patient would benefit from revascularization. Step 370 then ensues, and at query box E, the patient is labeled as someone with multi-segmental abnormalities in anterior, lateral and inferior walls based on a composite of echocardiographic and acceleration indices. As conditions DV=ALI and E=Vali are met, the Characterization Apparatus will define the patient as having significant dysynchrony in anterior, lateral and inferior territories with impaired wall motion but viable myocardium. The composite of this data is communicated to the clinician in any format. By way of example, when the values of A, B, C, DV, E fall within a specific range, a message is communicated in display box 500: “this patient has been identified as one who will have a positive outcome from revascularization and mitral valve replacement, based on patients with similar clinical findings this patients' prognosis will improve with revascularization and mitral valve surgery if other factors do not increase operative risk. Individual outcomes may vary.”
The ability of the system to generate such data is dependent on existing data sets/statistical probabilities of outcome (prognostic information) with different permutations of parameters (e.g. A-E) that are compared in 400 to findings in other patients who's outcome has been determined and stored in the CPC. Neural networks or other techniques are applied in 400 to update and ‘teach’ the system at periodic intervals from data made available from the CPC.
After revascularization (e.g. coronary artery bypass surgery), changes in the parameters may or may not still indicate that the patient would benefit from CRT. Post-operative assessment of intrinsic device based indices in conjunction with extrinsic imaging techniques such as tissue Doppler echocardiography can be used to determine further management. In a preferred embodiment, the extrinsic and intrinsic data acquired can be analyzed in a complementary fashion, compared to previous case examples with similar findings in the central processing center and the system can provide the clinician with statistical data detailing expected response rates and prognosis with and without CRT.
Though future generation pacing systems may include LV pacing capabilities, activation of resynchronization therapy may only appropriate at specific times (e.g. exercise induced conduction abnormalities, changes in electromechanical properties with remodeling/reverse remodeling) and thus, CRT activation can be implemented when needed. This will save costs to the system and battery longevity. When activated, the system can make real time recommendations about programming of interval timing or perform these automatically as part of a closed loop system. Similarly, activation of regionally delivered NPIT will occur at the appropriate times and to the appropriate leads when directed by such a closed loop system as described in more detail below. Activation of NPIT will occur when any cardiac performance index or indices indicate(s) impairments in cardiac performance. By way of example, if measurements of inotropy using implanted lead based accelerometer data showed a reduction in systolic cardiac performance the system will activate NPIT. In a preferred embodiment, activation of NPIT is based on assessment of a plurality of parameters as denoted in
The apparatus predictor algorithms need not rely on a single index but will incorporate redundancy techniques as to improve specificity much in the same way an operator may define hypokinetic myocardium if such findings are found in more than one imaging plane on a conventional 2 dimensional echocardiogram. This is also important as inadequate signal to noise ratios or the lack of specific sensors may compromise the systems ability to define certain characteristics. Thus, the system can rely on a plurality of signals, intrinsic and/or extrinsic in its analysis. This case is by way of an example and illustrates how any number of redundant or repetitive steps occurs with fast software algorithms.
In the example described herein this patient will be defined as someone likely to respond to revascularization surgery and possibly mitral valve replacement and CRT. This result can be output in any format (digital or analog), numeric or as written language. Neural networks or comparable means can analyze the available data, draw comparisons between the individual case and a library/data bank of other cases stored in the central processing center with outcome information prior to reporting any observations and before making recommendations in 400. The apparatus/algorithms may specifically report such characteristics via a graphical display and/or printed out by monitor/printer 500 as described above. This will allow an operator to review and manually input data based on subjective assessment of the acquired data (step 510) before final processing in 300 (
As the current invention allows for acquisition of data acquired from large population studies with entry of the data into a central data bank, the algorithms employed and recommendations generated (IDS) apply evidence-based medicine. Initially, such technology will be used for monitoring purposes. As more patients have access to equipment with such technologic capability, the more predictive the IDS will be for determining prognosis. Once validated with large scale clinical trials, formal recommendations can be relayed to the practitioner based on the predictor algorithms described herein. When possible, the device can perform automatic programming with closed loop control systems. Though the device based algorithms have the capacity to perform these functions, such sequential, advancements in the application of such technology can be achieved by software downloads (device upgrades) at the appropriate times (comprehensive central data bank, FDA approval).
It is readily apparent that such a technologic advancement holds promise to reduce health care costs and improve patient outcome as patients most likely to benefit from specific therapies will have access to such treatment and those patients unlikely to benefit from specific therapies will avoid the potential risks/health care costs of such treatment.
Outcome analysis performed within the device and/or a central processing center gathers and processes such data via wireless telemetry from extrinsic diagnostic equipment and/or implanted devices with such capability. Thus, at pre-programmed intervals such intrinsically and extrinsically derived data sends acquired data in numeric format (e.g. characteristics A-E) to the central processing center (
Any and all data received by the central processing center (CPC) is stored in a data bank, and tracked for all patients who have such data entered. The data can be from implanted devices (intrinsic) or external diagnostic equipment (extrinsic). Patient specific therapies and outcome are input to the CPC. The outcome data can be in the form of any index of cardiac performance or hard endpoints such as morbidity or mortality data. The data is processed via neural networks or other technique in 400 with the output data being any recommendations for patient treatment communicated to the user in 500 (e.g. recommendation for valve replacement) or activation of CRT or NPIT in an implanted device.
In this fashion, predictor algorithms in the CPC or within the implanted device/extrinsic diagnostic equipment can generate prognostic data for individual patients based on population studies using evidence based medicine. The pooled data is used to help predict outcome and guide treatment. The centrally processed pooled data can be made available to expert panels elected by the American Heart Association, various regulatory committees such as the Food and Drug Administration, and Centers for Medicare and Medicaid Services (CMS).
Non-Excitatory Stimulation and Global Tethering/Left Ventricular Torsion
The normal heart contracts and relaxes as a helix. Embryologic development of the heart is in a helical fashion. Viewing the heart from apex to base, this rotation is counterclockwise during systole and clockwise during diastole. Such helical rotational properties/torsion significantly contribute to cardiac performance, though this is currently not recognized in both academic and clinical circles. Such properties are lost in patients with dysynchronous contractile patterns and cardiomyopathy. Diastolic uncoiling and systolic coiling of the heart results in a suction effect for diastolic inflow and a twisting, wringing effect for systolic output. Pathologic myocardium and/or conduction abnormalities alter this pattern and lead to diastolic and systolic dysfunction.
The relative contribution of LV torsion to cardiac performance is significantly underestimated. Patients with severe impairments in contractility but absent heart failure symptoms may have retained cardiac performance as LV torsion is preserved. Likewise, patients who appear to have preserved contractility and advanced heart failure symptoms might be expected to have impairments in LV torsion. Electromechanical activation patterns are a third factor contributing to overall cardiac performance. Thus, the inventor recognized that three major mechanical factors, myocardial thickening, electromechanical timing and left ventricular torsion contribute to cardiac performance, patient symptomatology, and clinical outcome (
As mentioned above, initial studies in patients with systolic heart failure who have implanted devices capable of delivering non-excitatory electrical current have demonstrated improved systolic performance and reductions in heart failure symptomatology.
The inventor believes that CCM or other NPIT function by increasing regional contractility, and thus, increasing local systolic strain and peak strain rate but not necessarily in remote myocardial segments. NPIT increases systolic longitudinal displacement of regional (e.g. septal) segments but properties of distant territories (e.g. lateral) are effected as well. This occurs as a result of global tethering; tethering effects that affect far-field tissue velocity/motion despite regional current delivery (e.g. only in the interventricular septum). By augmenting strain and motion of the interventricular septum, NPIT will in turn increase the systolic rotation of the left ventricle and augment diastolic recoil/counter-rotation. These effects will increase LV torsion and improve both systolic and diastolic cardiac performance (
Development of multi-site pacing systems capable of NPIT promises to further enhance the benefits of such technology and avoid any potential for increases in dysynchrony resulting from regional non-excitatory stimulation. In patients with electromechanical dysynchrony, CRT-NPIT pacing systems can improve cardiac performance and synchronize myocardial activation. Because of the anisotropic properties of the myopathic heart, an individuals needs for resynchronization therapy and non-pharmacologic inotropic therapy may vary from myocardial segment to segment, patient to patient and time to time. In order to maximize a given patient's response to such device based therapies, specific algorithms and control systems will be necessary.
Multiple novel rotational/torsion parameters can be derived that describe the helical properties of the heart during the cardiac cycle and are ideal for analysis of the effects of NPIT. For illustrative purposes,
Patients with organized atrial tachyarrhythmias (e.g. atrial flutter) lack the normal A″ (
In the extrinsic apparatus designed to evaluate a plurality of signals, algorithms for identifying patients likely to respond to NPIT can utilize input physiological properties including septal contractility, septal and distal segmental longitudinal displacement, regional and global strain analyses and LV torsion into Analysis Apparatus 200 for numerical scoring (step 110 in
For explanatory purposes, we will describe the method and means by which such an algorithm functions using non-invasive imaging echocardiographic techniques and then illustrate how these methodologies are applied to be intrinsic to device based software. By way of example, referring to
In the current invention, the algorithms described above, are employed within a CRT and/or NPIT device in an analogous manner. Regional trans-myocardial changes in impedance are surrogate to strain measurements reflective of contractility and lead-based accelerometers/piezoelectric transducers are surrogate to measurements of rotation/torsion, longitudinal displacement. Referring to
The algorithms detailed herein are described using physiological indices derived from implanted cardiac devices and extrinsic imaging modalities as this field is in a state of rapid evolution. The principles underlying this invention can be applied to other fields that are similarly evolving at a rapid pace. Identification of genetic markers for predicting disease/outcome is one such field.
In an alternate embodiment, the concepts herein are extended to algorithms that incorporate data intrinsic to the body such as chromosomal properties related to genetic markers of disease (e.g. coronary artery disease, cardiomyopathy, long QT syndrome, oncogenes), circulating lipoproteins, serum markers; and extrinsic to the body such as results of various diagnostic tests (demographic data, phenotype, results of imaging modalities, pathologic findings from biopsies). This data can be input to 100 along with CP and DI and is relevant as an increasing number of genetic mutations have been identified that are predictive of disease states (e.g. progressive cardiomyopathy) and sudden cardiac death (Villard E, Duboscq-Bidot L, Charron P, et al. Eur Heart J 2005; 26:794-803). Treatment options including device based therapies, pharmacologic treatment and delivery of genetic material into pathologic tissue may be more appropriately chosen when the invention herein is widely available.
By way of example, analysis algorithms are being developed for detection of genomic properties in human cancer cell lines. These methods are currently employed using intrinsic or extrinsic markers for disease. Greater predictive power can be expected by combining both intrinsic and extrinsic data/signals.
Daruwala et al. has developed a versatile statistical algorithm to detect genome copy number variation (Daruwala R S, Rudra A, Ostrer H, et al. A versatile statistical analysis algorithm to detect genome copy number variation: Proceedings of the National Academy of Sciences of the United States of America, 2004. November 16; 101 (46): 16292-7). This algorithm analyzes genomic data obtained from a variety of array technologies allowing for detection of regions of the genome with altered copy number. The algorithm combines data from multiple ‘intrinsic’ sources and thus facilitates the discovery of genes and markers important for predicting risk of developing cancer. It implements a priorless maximum a posteriori estimator and dynamic programming implementation to increase computational robustness. This specific technique and others like it are being utilized to predict or detect cancer.
Algorithms that implement genetic markers for identifying patients prone to developing coronary artery disease are also being developed (Jones, Karen Ann et al., Atherosclerosis-associated Genes; patent application 20020015950). This algorithm is ultimately employed to treat or prevent “a disease associated with the altered expression of a gene that is coexpressed with one or more known atherosclerosis-associated genes in a subject in need, the method comprising the step of administering to the subject in need the pharmaceutical composition . . . in an amount effective for treating or preventing the disease”.
The health care costs associated with widespread genetic testing and ethical issues associated with using genetic predictors of disease for clinical decision making (e.g. pharmaceutical/genetic therapy, prophylactic surgery) are tremendous. Thus, a more accurate means for predicting outcome will lead to development of guidelines which determine who in the general population should be tested and which treatment options are preferred. This will reduce health care costs, reduce the frequency of and risks of unnecessary therapies/surgeries and prevent the emotional trauma from positive results in patients with an otherwise low risk for developing a specific pathologic state.
Risk assessment tools exist for predicting probability of developing cancer. One example is the Breast Cancer Risk Assessment Tool (http://brca.nci.nih.gow/brc/). This is a risk stratification tool that relies on extrinsic data and not genetic markers. Certain models/calculation methods have been developed such as the Gail Model and NSABP model (http://www.halls.md/breast/riskcom.htm). The construct of these models employed evidence based medicine. Genetic data will have more predictive power and would be weighed more heavily in the algorithm described herein. Nonetheless, combining extrinsic data such as that used in the Gail Model and intrinsic data that relates to genetic markers for disease will result in greater predictive power than either methodology alone.
Referring to
As described above, several factors can be used to derive a numeric score for properties A-F. A can relate to genomic copy number variation and presence of other genomic traits (e.g. EVI1, MYC genes). B can relate to presence of and type of disease-associated germlne mutations (e.g. BRCA1, BRCA2). C can be a genotypic descriptor related to likelihood of genetic expression.
Extrinsic data can include demographic data (history of tobacco or alcohol intake), results of imaging studies (e.g. breast sonograms, mammograms, MRI), pathologic findings from biopsies and family history of cancer. Extrinsic physiological descriptors can represent, D, family history, E, results of breast MRI, and F, biopsy findings. By way of example, D can be scored based on heritage (e.g. Ashknenazi descent), number of immediate (first-degree) and second-degree relatives that have been diagnosed with specific cancer types, the life-expectancy of such relatives, response of familial cancer to therapy in a weighted fashion (e.g. Discriminant Analysis). E can relate to MRI properties such as the anisotropy of breast tissue density, presence of specific findings associated with increased risk of breast cancer. F can relate to pre-malignant findings (e.g. dysplasia, papillomatous tissue, ductal carcinoma in situ) on breast biopsies performed on a given individual or even an individual's family members.
Specific physiological indices are entered into the algorithm using a discrete numeric format along a graduated scale as described above. Physiological descriptors such as presence of BRCA mutations (descriptor B) and biopsy pathology data (descriptor F) will be weighed more heavily than other less predictive descriptors (e.g. using Discriminant analysis). The output informational data set provided to the physician will guide decision-making by providing an estimation of probable outcome (e.g. percent risk of developing invasive breast cancer, percent risk of dying from breast cancer) and the effect of various therapeutic regimens (e.g. increased surveillance, treatment with Tamoxifen, future genetic therapies, prophylactic mastectomy/oophorectomy) for modifying such an outcome. The central processing center, CPC, and pooled data bank would ideally be located where chromosomal analysis is performed (i.e. Myriad laboratories) as few centers are currently available for such genetic testing. It would be mandatory for any patients having such testing to have all demographic data, family history and other descriptors sent to the CPC. This data in turn would be communicated into 400 for processing. In this example, apparatus 300 is an separate apparatus or contained within extrinsic diagnostic equipment such as a breast MRI machine or mammography unit.
Thus, the spirit and scope of this invention relates to a variety of pathologic states and not just those related to the field of cardiology.
Translation Function
When this invention is applied to implanted cardiac device technology, mathematical methods for correlating the various indices can be used as to allow translation of device based data to imaging based data and vice versa. Novel mathematical units that are assigned a value along a graduated scale can apply to both intrinsically and extrinsically derived data. These units can be derived as to develop a universal mathematical language that is applicable to both device-based and non-device based sensors/transducers. By way of example, a number of patients who have had an impedance based dysynchrony parameter, CPPz, acquired/derived from an implanted device, and a multi-dimensional echocardiographic measurements of dysynchrony, SMDI, derived can be compared (though not necessary for understanding this concept, SMDI represents the Systolic Myocardial Dysynchrony Integral, a triple integration in polar coordinates of the Curved Anatomic M-mode myocardial dimension, coursing around a 360 arc about the heart, between base to apex, endocardial to epicardial surfaces, where the limits of integration are between the time of aortic valve opening and time of end diastolic motion of all myocardial segments; the derivation of SMDI can be found in the published patent application Ser. No. 10/779,162). Such a comparison can include patient demographic data, clinical variables and patient outcome. Regression or other statistical analysis of extrinsic and intrinsic parameters is used to correlate intrinsic and extrinsic parameters. The relationship need not be linear but may be exponential, curvilinear, determined by a derived differential equation or otherwise correlated in some other fashion. By way of example, SMDI and CPPz can have a linear relationship. Thus, SMDI and CPPz are directly proportional along a graduated scale between normal and pathologic values (
A translation function can, in a more simplified approach, be used to translate mathematical units derived from one type of implanted sensor into conventionally recognized units. For example, displacement data derived from lead based accelerometers to indices used with tissue Doppler or speckle tracking techniques (e.g. longitudinal mm distance, radians rotated). This can be defined for an individual patient based on simultaneously acquired data with, for example, an ultrasonic interface connected ('162) or by using a look up table developed from pooled data (registry) derived from numerous patients who have had both intrinsic and extrinsic data derived. Regression analysis or other technique can be used to correlate any intrinsically and extrinsically derived index (
Referring to
(CPPZ)=m*(SMDI)
A linear relationship between intrinsic impedance parameter CPPz and extrinsic ultrasound parameter SMDI is present. Slope m relates how CPPz and SMDI changes as a function of one another.
This relationship can be derived as a function of time during the cardiac cycle and in this example SMDI becomes an instantaneous measurement of cardiac dysynchrony, MDI, and CPPz an instantaneous impedance based measurement of dysynchrony which may be directly proportionate to each other. (Cardiac motion/deformation can be symmetric and synchronous at one point during the cardiac cycle (e.g. end diastole) and not at another time frame (e.g. end-systole)).
integral Z(t) dt α MDI(t) dt
To further elaborate on the concepts described, the reader is referred to
In
In these examples, the peak displacement seen at end-systole is noted to vary under different circumstances. If the same patient has NPIT delivered to lateral leads, the resulting curves may become more disparate. Simultaneous delivery of NPIT to septal and lateral electrodes in this patient may also lead to less congruence in displacement curves (
Circumferentially located lead based accelerometers placed about a CS lead will enable an implanted system to derive rotational indices (e.g. peak systolic radial displacement, time of peak rotational velocity) and thus the system will be capable of evaluating all three elements of cardiac performance, The Three T's; thickening, timing and torsion (
Referring to
In the latter algorithms described, a plurality of signals is derived from one transducer/sensor type without need for multiple sensor technology. The accrued data represents more specific information and translation between intrinsic and extrinsic indices will be relatively simple. Translation can be performed by comparing external monitoring data and device based data acquired from a single individual (e.g. manually downloaded, via an interface or wireless telemetry). Translation can also be performed for large samples of patients (from a centralized database/registry) in a more generalized fashion. When multiple sensor technology is implemented, as detailed in the algorithms described previously, the resulting informational data set has broader application and prognostic information and clinical recommendations can be made using neural network techniques, for example. Thus, mathematical indices derived using multiple sensor technology can be incorporated into algorithms that generate prognostic information. The resulting informational data sets will have more widespread ramifications than the specific information provided by an individual index (e.g. displacement) which is best suited for optimizing device function through closed loop control systems. Device software algorithms can implement more complex, device-based, closed loop control systems using multiple sensor technology if the costs to the system and processing times are not excessive.
Closed Loop System Function Confirmation
In an alternate embodiment, this technology can confirm the appropriate functioning of closed loop systems within implanted devices from time to time. Similar physiological descriptors may be acquired by implanted devices and external diagnostic equipment and constitute multiple complementary data sets. As an example, myocardial acceleration can be derived using tissue Doppler techniques or implanted lead-based accelerometers (
Data acquired by any imaging modality can be made available for entry into an apparatus that performs such complex algorithms in an efficacious manner. In this fashion, a physician can render diagnoses and make management decisions based on a number of indices simultaneously rather than relying on independent physiological properties. Digital signal processing allows for imaging data derived from alternate imaging modalities to be compiled in a central processing unit and analyzed in context of similarly derived parameters acquired from implanted devices. Wireless data transfer can be via a number of means such as the Internet, removable digital media or telecommunication. The gathering of voluminous data banks/libraries/registries of data acquired from external monitoring equipment, implanted devices/sensors coupled with patient demographic data and clinical outcomes will lead us to practice evidence based medicine. The derived data will be used to develop treatment guidelines and clinical recommendations appropriate for each patient as an individual. In turn, this technology will reduce health care costs and benefit society as a whole.
The spirit and scope of the concepts described herein are applicable to a variety of diagnostic modalities, monitoring systems, disciplines, implanted devices and a variety of indices derived by external and internal means. This information can be communicated to a physician any number of ways and sent via wireless telemetry to an electronic medical record. Derivation of universal mathematical indices/units of expression will facilitate translation and bi-directional communication of diagnostic data between an implanted device and external diagnostic equipment and confirm functionality of closed loop control systems and extrinsic diagnostic imaging modalities (cross-verification).
Data Communication
Parametric imaging techniques can be used to illustrate intrinsically derived data such as time of peak impedance/acceleration/displacement or device based indices of cardiac performance. In a preferred embodiment, intrinsically derived indices that are analogous to extrinsically derived data are displayed using color encoded pixels superimposed on images of the heart generated, for example, by echocardiography. By way of example, regional differences in time of peak impedance or time of peak displacement derived from implanted sensors can be displayed in similar fashion to tissue Doppler derived indices portrayed on state of the art echocardiography machines.
Many other modifications may be made to the invention without departing from its scope as defined in the appended claims.
Claims
1. A method of generating cardiac information comprising:
- comparing pooled diagnostic data stored in a centralized data bank with data collected from a plurality of similarly acquired data from individual patients including sensor-derived data from within the patients' implanted cardiac rhythm management device; and
- deriving information from said comparison, said information being representative of the pathologic state of the individual patient.
2. The method of claim 1 further comprising the providing of one or more of a prognosis, recommended treatment options and programming of said cardiac rhythm management device.
3. The method of claim 1 wherein said comparing and deriving is performed in stand-alone equipment.
4. The method of claim 1 wherein said comparing and deriving is performed in a patient's cardiac management system.
5. The method of claim 1 further comprising using said information to confirm appropriate functioning of cardiac rhythm management devices at periodic intervals.
6. A method of generating programming parameters for implantable devices comprising:
- combining a plurality of one or more of patient demographic information, genetic characteristics and results of diagnostic testing from any given patient to generate patient specific prognostic information based on statistical comparisons to similar data acquired from patients whose outcome under varying treatments is stored in a central processing center.
7. A method of generating programming parameters for cardiac rhythm management devices comprising:
- collecting analogous data from cardiac rhythm management device-based sensors and external diagnostic equipment; and
- using said data to obtain said parameters, said parameters being selected to control the frequency and duration of pacing modalities selected to increase cardiac contractility.
Type: Application
Filed: Jan 19, 2006
Publication Date: Jul 27, 2006
Inventor: Stuart Schecter (Great Neck, NY)
Application Number: 11/334,935
International Classification: A61N 1/08 (20060101); A61N 1/372 (20060101);
