U.S. Patent Documents
U.S. Pat. No. 4,417,306 November, 1983 Citron et al.
U.S. Pat. No. 5,419,338 May, 1995 Sarma et al.
U.S. Pat. No. 5,437,285 August, 1995 Verrier et al.
U.S. Pat. No. 5,560,368 October, 1996 Berger.
U.S. Pat. No. 5,560,370 October, 1996 Verrier et al.
U.S. Pat. No. 6,132,381 October, 2000 Forbes et al.
U.S. Pat. No. 6,324,423 November, 2001 Callahan and Shell
U.S. Pat. No. 6,577,894 June, 2003 Callahan and Shell
Other References
The present invention relates to measuring cardiac function intervals.
It is known that alteration of the QT, QTc or RR interval on the electrocardiogram may be a marker for sudden death(1-15). Measurements of the QT interval are generally taken from a 12-lead electrocardiogram where one to three heart beats are analyzed either individually or averaged(16;17). The 12-lead electrocardiogram provides only point-in-time data reflecting approximately 17 seconds of time that is required to inscribe a 12-lead ECG. The QT interval duration is dynamic, however, and can vary by upwards of 100 msec in a twenty four hour period(18-26). Thus the measurement of either a single or a few 12-lead ECGs sampled during 24 hours will miss the beat-to-beat dynamicity data that is inherent in the changes that occur. The dynamic data reflecting changes in ECG intervals is captured by longer recordings, generally 24 hours of continuous ECG data, referred to as either 24 hour ambulatory ECG (AECG) or Holter Monitoring(18;27-36). Heretofore, beat-to-beat ECG data, both short and long-term recordings, has been averaged due primarily to constraints in computing power. Unfortunately, averaging minimizes the understanding of the beat-to-beat variability inherent in QT interval data. Moreover, methods to analyze large data sets of cardiac intervals have been incomplete. For example the methods for beat-to-beat binning of QT and QTc intervals described by Callahan and Shell where limited to analysis of only outliers(37), calculating the % of beats that exceed a certain threshold. The disclosures by Shell and Callahan do not teach a method to analyze central tendency, variance, kurtosis or other statistical properties of the histogram as appropriate for Gausian or non-Gaussian distributions.
Increases in the QT and QTc interval measurements on a 12-lead Electrocardiogram (ECG) are associated with an increased risk of cardiac dysrhythmias and sudden cardiac death. See, for example, Algra(38), Schwartz(39-41) and Sawicki(42;43). The increased QTc interval length is associated with an increased risk of sudden death from all causes. The prolongation of the QTc interval induced by pharmaceuticals has been associated with Torsade de Pointes and sudden death; the pharmaceutical induction of prolonged QTc intervals has formed the basis for removal of pharmaceuticals from the market. There is, however, no readily agreed upon method to measure the dynamic changes in the QTc interval, particularly for long term recordings of the ECG.
While the resting 12-lead electrocardiogram may provide important spatial information regarding the status of ventricular repolarization, the use of a single 12-lead ECG measured randomly in time may disregard potentially important prognostic data regarding the dynamicity, temporal relationships, and circadian rhythms of the QT interval.
It is known that the QT interval may undergo significant changes over both the short and long term due to circadian rhythms. See, for example, Yi, et al(44) who teach the association between circadian rhythm and sudden death associated with acute myocardial infarction. See also, for example, Callahan and Shell who describe a method to assess circadian changes in the QT interval.
It is known that the QTc interval may undergo significant changes over both the shorter and longer term due to autonomic control. See, for example, Cappatto et al, Browne et al(45), and Kautzner, et al(46;47), demonstrated the relationship between sympathetic and vagal tone on the QT and QTc interval.
Thus, a single 12-lead ECG taken at a given point in time may provide misleading and inaccurate cardiac risk data. Therefore, analysis of the QT interval for an entire 24-hour period, reflecting circadian and autonomic changes, may provide additional information regarding the risk of sudden death not available on the single, random 12-lead ECG.
It is now possible to measure the QT interval on 24-hour Holter (AECG) recordings(18;29;31;48-64) These measurements have generally been reported as averages over short time periods, typically between about 15 seconds and about five minutes, for example Molnar et al(65-67) or Yanaga, et al(68). The use of averaged QT measurements may obscure significant short-term variations in the QT intervals. Conversely, beat-to-beat measurements retain the natural variability data that may be important for calculating a patient's risk of dysrhythmia and sudden death.
More recently beat-to-beat QT interval measurements have been used but methods to analyze the beat-to-beat changes have been incomplete.
Although beat-to-beat variability of the QT interval has been described by Berger and others (69), little is known regarding normal ranges in variability and measures of the QT interval over a 24-hour period using beat-to-beat measurements.
Molnar and colleagues published a study that gives some indication of the dynamic range of the QT intervals using five minute averages and not beat-to-beat measurements(70). They reported a mean maximum QTc interval of 495 ms for normal subjects using 24-hour ambulatory monitoring. They also showed a mean intra-subject change of 95 ms. Molnar further reported six normal female subjects as having a maximum mean QTc interval measurement of more than 500 ms. These mean maximum measurements were taken over a five-minute period.
The use of average QTc measurements obscures the dynamicity of individual beats. Measurements of central tendency, skewness and shape of histograms have not been used extensively to describe the relationship of QT and QTc measurements in histograms representing beat-to-beat QT, QTc or RR intervals. These measurements may be important to give an overall picture of the status of the subject.
It is an objective of the present invention, in a preferred embodiment, to enable the assessment of the QT and QTc intervals and other cardiac function intervals on a beat-to-beat basis, providing a composite histogram of the individual beats with QT and QTc intervals.
It is another objective of the present invention, in a preferred embodiment, to enable the measurement and assessment of the QT and QTc intervals and other cardiac function intervals over an extended period of time, including not only periods of time greater than about one minute but also periods of time lasting at least 24 hours and even longer, in some cases.
In accordance with the present invention, in a preferred embodiment, this and other objectives are achieved by providing a method for analyzing beat-to-beat QT intervals from high-resolution Ambulatory Electrocardiographic monitoring (AECG) to detect the frequency distribution in a continuous AECG recording. Beat-to-beat QT and RR intervals may be measured to calculate beat-to-beat QTc. In a preferred embodiment, a composite of the entire frequency distribution of QT and QTc intervals taken from a set of observations with a common characteristic may be examined. Moreover, a composite of one characteristic may be statistically compared to a composite with other characteristics, including statistical methods that do not assume a normal distribution of the histogram.
The present invention, in a preferred embodiment, provides a method to analyze beat-to-beat QT data, stratify the data according to a time-series bin-array, and create a composite of multiple histograms. This method and apparatus may be applicable to a wide variety of different subjects including, for example, normal subjects, subjects with the Inherited Long QT syndrome (ILQTS), and subjects exposed to drug titration. The statistical characteristics of a normal subject group can be compared to either a second group, or to individuals who have taken a drug, have potential congenital heart disease, have been exposed to an environmental toxin, or have a disease which could cause prolongation of the QTc interval such as diabetes mellitus.
Further objects, advantages and other features of the present invention will be apparent to those skilled in the art upon reading the disclosure set forth herein.
The following detailed description is of the best presently contemplated mode of carrying out the invention. This description is not to be taken in a limiting sense, but is made merely for the purpose of illustrating the general principles of the invention. The scope of the invention is best defined by the appended claims.
In a preferred embodiment, standard 24-hour AECG recordings may be obtained using any commercially available Holter cassette tape recording device. An example of this type of device is a Reynolds Medical Tracker II (Reynolds Medical, Hertford UK) recorder. Also in the preferred embodiment, standard 24-hour AECG recordings can be obtained using commercially available digital Holter recorders that have a sufficient sample rate to allow detection and measurement of the cardiac intervals. For QT interval analysis, the sample rate can be between 128 and 2000 samples per seconds. The preferred embodiment is the use of a sample rate of at least of at least 1000 samples per second. These digital recorders must also be compatible with a Holter playback system that can produce beat-to beat interval measurements. An example of this type of device is the Reynolds Medical LifeCard CF recorder (Reynolds Medical, Hertford, UK). An example of the compatible Holter playback system is the Reynolds Medical Pathfinder 700 series. (Reynolds Medical, Hertford, UK). These recorders and playback systems are commercially available and need not be modified.
Analog signals from the Holter cassette recordings may be digitized at 12-bit or higher resolution using a Holter playback system that has the ability to perform interval measurements. The Reynolds Medical Pathfinder 700 series Holter analyzer is an example of this type of equipment (Reynolds Medical, Hertford, UK).
Using the digitized file of the electrocardiogram, a QT interval analysis may be accomplished in the following manner: The onset of a Q-wave (Qb) may be defined and a cursor may be placed at this point. The end of a T-wave (Te) may be defined and a second cursor may be placed at this point. The data from the digital file may then be replayed at 60-times normal time, while the cursors on the Qb and Te points may be monitored for stability. If either cursor wavers from the Qb or Te points, the cursors may be replaced and the affected portion of the data may be reanalyzed. The QT interval may be defined as the time difference between the time points at Qb and Te. The QT intervals may be measured for the entire AECG recording on a beat-to-beat basis. Other analysis systems that display digital data may be used.
The peak of an R-wave may be detected and a third cursor may be placed (Rp). Accordingly, each QT interval may be matched with the preceding R--R interval. For a 24-hour recording, this may result in approximately 100,000 beats for which a QT interval and an R--R interval may be defined. The data may then be output to a high-speed computer for post-analysis processing.
In the examples described herein use was made of AECG recordings from normal volunteers, subjects treated with placebo and subjects on-treatment in a drug treatment study, and recordings from subjects with inherited Long QT Syndrome (ILQT). These recordings help to demonstrate the potential effectiveness creating composite curves in accordance with the present invention.
In the examples described herein QTc was calculated by removing a time-series of the QT and preceding R--R intervals to a high-speed computer with both a fast processor and adequate disk storage space. For each QT interval, a QTc may be calculated using a variety of correction factors for the QT interval including Bazett's correction formula, Fridericia correction formula and linear correction formula.
The QT and QTc intervals may be individually placed in the bins according to their measurement as described in Shell and Callahan. In a preferred embodiment the composite curves are constructed by software programs that generate a time series of approximately 100,000 data points long of RR/QT/QTc triplets for each patient. Then the QTc data for each patient is binned in a histogram for that patient, finally, software is used to merge many patients'data into a composite data set (a“population”) and to take means and standard deviations of this population (assuming normalcy of the data). Finally, more the data thus aggregated into two or more populations can then be compared, again using a combination of software and procedures as described in Press et al(71), against each other to check for statistical difference between these two or more populations.
These aggregated population curves can then be used as a template for comparison against an single patient's binned histogram to determine what population (e.g. normal, inherited disorder, or drug-induced damaged) this particular patient belongs. The current embodiment assumes normal distributions, but this is not intrinsic to the method and more sophisticated distribution-distinguishing numerical analysis and statistics is declared here as well. The numeric procedures used are commonly described in Press, et al.
In this example, the 24 hour ambulatory ECG from a normal subject was analyzed. The 24-hour ECG was digitized. The QT and RR interval was determined for each beat using the Reynolds's analyzer. All extra beats were eliminated. All beats with prolonged QT intervals were inspected and artifact was eliminated. The QT and RR files were then used to construct a Histogram of QT and QTc intervals. The histogram was constructed with 10 msec intervals. The histogram of QTc intervals is depicted in
PROLONGED QTc Intervals in Inherited Long QT Syndrome.
The Inherited Long QT Syndrome is a genetic defect of the heart's ion channels. The patients with Inherited Long QT Syndrome are known to have intermittent prolonged QTc intervals. Often, however, many of the heart beats of patients with inherited long QT syndrome are within the normal range and the identification of these patients cannot be made from a single conventional 12 lead ECG. Since these patients, often children, die suddenly, failure to detect the presence of the abnormal gene can lead to sudden death of the infant, child or young adult, an unnecessary death since treatment is available to prevent such sudden death. In this example., a child with a known gene defect underwent 24 hour ambulatory monitoring. The ECG was digitized and the QT and RR intervals defined. The QT and QTc histograms are depicted in
Patients with Drug Induced Long QT Interval.
Many drugs can prolong the QTc interval and the drug induced prolonged QTc interval is associated with an increased incidence of sudden death. Many drugs have been removed from the market because they prolong the QTc interval. Cisapride is a drug that can prolong the QTc interval. In
Comparison of the QTc Interval Using Composite Curves in Normal Subjects to those with Inherited Long QT Syndrome and Drug Therapy
In
Comparison of the QTc Interval Using Composite Curves in a Subject Before and After Drug Intervention
Since one of the important uses of this methodology is to compare the QTc interval before and after the use of a pharmaceutical that could prolong the QTc interval, we compared a group of patients before and after the administration of a pharmaceutical. The 19 patients had 24 hour ambulatory ECG monitoring before and after the administration of drug. The composite curves before and after the administration of drug are depicted in
We then compared the means in these two groups using a standard T-test and analysis of variance. The p-value for the difference was less than 0.000008. The creation of the composite curves allows definitive differentiation of the three groups.
Comparison of a Single Individual to a Composite Set of Data
Frequently, one is confronted with the problem of defining if a set of QTc data is derived from a set of normal data. In this example, a single individual with Inherited Long QT Syndrome was compared to a set of normal subjects (
In a preferred embodiment, the present invention represents a new method for quantifying the QT and QTc interval measurements over a period of time. The invention allows a quantitative comparison of two or more sets of QT or QTc intervals. For example, the invention allows comparison of a group of patients before and after a drug. The method described allows application of a variety of statistical methods to define whether two or more sets of intervals are different from one another.
In a preferred embodiment, the method and apparatus may make use of high-speed computer processors, and large capacity data-storage media. In a preferred embodiment a 1 GHz Pentium IV processor with an 80-gigabyte hard drive may be used to analyze and store the large data files. Several custom-built software programs are used to generate a time series of approximately 100,000 data points of RR/QT/QTc measurements for each patient. Then the QTc data for each patient is binned in a histogram for that particular patient, finally, a combination of software and procedures are used to merge many patients'data into a composite data set (a“population”) and to take means and standard deviations of this population (assuming normalcy of the data). Finally, the data thus aggregated into two or more populations can then be compared, again using a combination of custom software and procedures, against each other to check for statistical difference between these two or more populations.
These aggregated population curves can then be used as a template for comparison against an single patient's binned histogram to determine what population (e.g. normal, inherited disorder, or drug-induced damaged) this particular patient belongs. The current embodiment assumes normal distributions, but this is not intrinsic to the method and more sophisticated distribution distinguishing numerical analysis and statistics is declared here as well.
Composite QTc histogram measurements in accordance with the present invention allows for a quantitative assessment of the number of specified intervals, such as QT and QTc in a 24-hour AECG recording.
The present invention, in a preferred embodiment, is directed to a method for the quantification of beat-to-beat QT and QTc interval measurements from ambulatory electrocardiographic recordings.
A QT binning technique in accordance with the present invention may be used to provide information about the effects of a pharmaceutical. For instance, in the example illustrated in
Although the preferred embodiment of the present invention has been described herein with respect to measurement and analysis of the QT interval, it will be recognized that a method in accordance with the present invention may also be useful in the measurement and analysis of a wide variety of other ECG and related biologically significant intervals
In a preferred embodiment, the method takes discreet measurements and discreet intervals and places them into a time series bin or an amplitude series bin. For example, all of the RR intervals in a sample could be selected and coded according to their length and then placed into bins. Each bin could be characterized by a frequency. The same analysis could be performed using any interval on the electrocardiogram.
The presently disclosed embodiments are to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims, rather than the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein.