The claimed invention relates to Electrocardiogram-based analysis systems and methods. More particularly, the claimed invention relates to ECG-based differentiation of LQT1 and LQT2 mutation in patients.
The electrocardiogram (ECG) is based on the electrical activity of the heart muscle cells. In the resting stage, the inside of the cardiac cells has a negative charge compared to the outside of the cells. The resulting voltage difference between the internal and the external spaces of the cell membrane is called transmembrane potential. The discharging of this voltage is known as depolarization and is associated with the start of the contraction of the heart muscle cell fibers. After contraction of the ventricles, the heart muscle cells redevelop substantially the same voltage over the cell membrane. This recovery phase is called the repolarization process of the heart ventricles. An ECG measured from the skin surface measures a total electrical component created by the depolarization and repolarization of the heart's muscular cells.
The repolarization of the heart is made possible in part by ion channels within the myocardial cells of the heart which allow an ion current to redistribute charge. It is highly important that the regulation of the ion currents during the ventricular repolarization process occurs without interference, since a delay in this process or any other abnormalities can lead to a substantially increased risk for sudden cardiac death.
Medical professionals have used electrocardiograms (ECG's) to examine the ventricular repolarization period, also known as the QT interval, to check for elongation of the QT interval. In general, an elongated QT interval may be considered to be indicative of a delay in the ventricular repolarization process. While medications in some cases may be the cause of an elongated QT interval, for many people are predisposed to have an elongated QT interval due to one or more congenital mutations. Those patients having a congenital predisposition for an elongated QT interval are referred-to as having Long QT Syndrome, or LQTS. The clinical course and the precipitating risk factors in the congenital Long QT Syndrome (LQTS) are genotype specific. Among LQTS mutations, KvLQT1 (LQT1) and HERG (LQT2) mutations have the higher likelihood of recurrent cardiac events and their diagnosis is crucial to reduce lethal outcome.
The Long QT Syndrome (LQTS) is an inherited disease caused by genetically determined defects in trans-membrane ion channel subunit. LQTS patients are at high risk of sudden cardiac death due to the development of ventricular tachycardia degenerating in ventricular fibrillation and cardiac arrest. The prevalence of the syndrome may be expected to occur in 1 in 3000-5000 individuals per year in the United States. The number of cases of sudden cardiac death associated with the LQTS is unknown but among the 300,000 sudden cardiac deaths documented each year, one may expect around 2-5% having LQTS-related arrhythmic death. In the US, the syndrome remains an under diagnosed disorder because an estimated 10 to 15% of the LQTS gene carrier patients have a QT interval duration near normal values.
Seven mutant genes have been associated with LQTS: KvLQT1/minK, HERG, SCN5A, Ankyrin B, KCNE1, KCNE2 and KCNJ2. Most of these genes encode cardiac ion channels and their mutation leads to dysfunction of the ion current kinetics. Among the current 150 mutations identified in the seven LQTS genes, LQT1 and LQT2 represent the majority of cases (88%) whereas LQT3 account for (7%) and the others are very rare. Typically in these types of LQTS, the arrhythmias occur in conjunction with vigorous physical exercise and emotional stress. The mutation of the HERG (LQT2) gene decreases the rapidly activating delayed rectifier potassium (K+) current (IKr). KvLQT1 mutation (LQT1) is associated with blockade of the slowly activating delayed rectifier K+ currents (IKs). The inhibition of IKr/IKs ion currents is associated with a prolongation of the action potentials within the heart leading in general to a prolonged QT interval on the surface ECGs.
Once a patient is identified as having an elongated QT interval, there is still a clinical need for discriminating mutation-specific syndrome (LQT1 vs. LQT2 patients) because of the mutation-specific therapeutic strategies one may consider. In this case, the QT/QTc prolongation is not a useful marker, the QT interval is, in average, prolonged similarly within the two groups. Traditionally, genetic testing has been prescribed to determine which form of Long QT Syndrome is the cause. Unfortunately, genetic testing can cost thousands of dollars and is time consuming. The expense of the testing alone prevents many patients from receiving important information about the root cause of their elongated repolarization interval. Therefore, there exists a need for an economical, reliable, and quicker way to differentiate between patients with the LQT1 and LQT2 mutation.
A method differentiating LQT1 mutation from LQT2 mutation is disclosed. An ECG signal is obtained for a patient. At least a first ECG parameter and a second ECG parameter are determined from the ECG signal. A probability that the patient is an LQT1 carrier or an LQT2 carrier is determined based on a regression model which takes into account the first ECG parameter and the second ECG parameter.
A system for assessing repolarization abnormalities is also disclosed. The system has a processor configured to differentiate between LQT1 and LQT2 based on at least two ECG parameters from ECG data. The system also has a data input coupled to the processor and configured to provide the processor with the ECG data. The system further has a user interface coupled to either the processor or the data input. It is at least one goal of the claimed invention to provide an improved method and system for the identification and differentiation of LQT1 and LQT2 mutations in patients.
It will be appreciated that for purposes of clarity and where deemed appropriate, reference numerals have been repeated in the figures to indicate corresponding features, and that the various elements in the drawings have not necessarily been drawn to scale in order to better show the features.
A surface electrocardiogram (ECG) may be measured by an ECG capture device which can have one or more leads which are coupled to a person's body in various locations. The electrical activity occurring within individual cells throughout the heart produces a cardiac electrical vector which can be measured at the skin's surface by the ECG capture device leads. The signal registered at the skin's surface originates from many simultaneously propagating activation fronts at different locations, each of which affects the size of the total component. One type of ECG capture device is a twelve-lead signal device, although ECG capture devices of any number of leads may be used to gather a set of ECG signals for use in assessing repolarization abnormality.
Using various scalar and vector-based parameters determined from one or more ECG signals, it may be determined whether or not a patient with Long QT Syndrome has the LQT1 genetic mutation or the LQT2 genetic mutation as the cause of the Long QT Syndrome by using a regression model which takes into account at least two of the scalar or vector based ECG parameters.
A study has been done to show the viability of this approach as follows:
Experimental Results:
Population
The study population consists of 49 LQT1 and 25 LQT2 carrier genotyped patients from 26 LQT1 families and 19 LQT2 families. These patients are all enrolled in the International Registry for the LQTS in whom 12-lead Holter ECGs were recorded. The HERG and KvLQT1 mutations were identified in each subject using standard genetic tests.
ECG Recordings
Twenty-four hours 12-lead Holter ECGs were acquired using the H12 recorders from Mortara Instrument (Mortara Instrument, Milwaukee, Wis.). This equipment provides 24-hour Holter digital ECG signal at a sampling frequency of 180 Hz and with 16 bit amplitude resolution (2.5 μV).
Eight true leads are recorded and the remaining four leads (augmented limb leads aVR, aVL, aVF and lead III) were computed.
Measurement Technique
All measurements are based on median beats from 10 consecutive cardiac cycles throughout entire 24-hour Holter recording. Only beats with stable heart rate (HR) are taken into account. The HR-stability assessment was based on the computation of the average HR within the 10 beats. This set of beats was accepted if all beats met the following criterion:
We investigated electrocardiographic and vectorcardiographic parameters in order to better characterize and understand the abnormalities differentiating LQT1 from LQT2 patients. To exclude circadian rhythms influence, we focused our work on the diurnal period.
We measured classical repolarization ECG measurements as the QT interval, Q to Tpeak interval (QTpeak), Tpeak to Tend interval (TpTe) and the magnitude of the T-wave (Tmag) from lead II and lead V5, using our own developed software for Comprehensive Analysis of the repolarization Signal (COMPAS). The vectorcardiographic measurements were based on the Principal Components Analysis (PCA) of the repolarization segment defined between the J point and the point located 220 ms before the next R peak (the shorter one in the set of 10 continuous beats) in order to ensure that the analysis encompasses all components of the ventricular repolarization. Those skilled in the art are familiar with this type of PCA analysis. The PCA analysis was applied to the eight original leads.
Briefly, the method relies on the computation of the singular value decomposition (SVD) [B9], in which any matrix A (M×N) can be written as:
A=USVT
where U is a (M×N) column-orthogonal matrix, S is a diagonal matrix with positive or zero elements (the singular values), and V an (N×N) orthogonal matrix containing the right singular vectors.
Being A a matrix (M×8), having in columns the 8 leads and in row the M samples of the repolarization signal of one median cardiac beat, applying SVD, we get U, S and VT matrices; S will contain the 8 eigenvalues in descending order. Multiplying AV, provides the projection of our original data on the principal components.
The repolarization signal within the space defined by the three first components (ev1, ev2, ev3) will be called the T-loop signal, note that it was not constructed from the original orthogonal Frank leads. The evn signals are the eigenleads. The plane (ev1ev2) defines the preferential plane of the T-loop. Following the dipolar theory of electrocardiography, the T-loop represents the path followed by the cardiac vector (VECG) during the repolarization process of the heart ventricles.
PCA measurements were obtained from the COMPAS PCA analysis package, which offers standard PCA parameters: complexity of repolarization (λ2/λ1), T-loop planarity (λ3) and other T-wave and T-loop morphology parameters such as the right (αR) and left (αL) slopes of the T-wave, computed on ev1. It also includes more complex repolarization parameters, based on the morphology of the T-loop in the principal components, such as maximum vector (MV), the interval from beginning on QRS to the instant at MV (LQTa) and QT interval calculated on T-Loop (LQTe), where the end of the T-wave is identified as the point Te:
The end of the T-wave is identified from the T-loop in its preferential plane. The detection of the end of the T-loop is defined as the earliest point of the current beat to meet the two criteria: 1) to be inside the circle of radius equals 25% of MV and centre in the point identified as the projection of the beginning of QRS on the preferential plane and 2) to have a repolarization velocity inferior of 3 mV/s.
Statistical Analysis
The ECG and VCG parameters were analyzed using multivariate analysis involving binary logistic regression to design models selecting the most discriminating parameters in the subset of factors we have investigated. A best subset regression model was selected based on the Akaike Information Criterion. Averaged values of parameters were compared using non-parametric test (Kruskal-Wallis). A p-value<0.05 was considered statistically significant.
We use an analytical strategy in which a referential model based on classical ECG parameters is first implemented then, a second model including both classical and vectocardiographic parameters quantifying T-wave and T-loop morphology is investigated.
The comparison between models is based on the area under the curve (AUC) of the receiver operating characteristics (ROC). The ROC curves are obtained by varying the threshold on the probability of being LQT2 patient from 0 to 1, this probability for each patient is provided by the binary logistic model.
Results
The clinical characteristics of the study population and the values from the 12-lead ECG parameters for the two populations are summarized in
The role of heart rate (HR) on the repolarization measurements is fundamental regardless of the type of measurements one can consider. In the discrimination between LQT1 and LQT2 patients, the effect of HR cannot be neglected even if in average the HR is not significantly different between the two types of mutations (
Initially, we focused our investigation on the limited HR range in which the number of subjects was maximal for the two groups. We chose the range in which RR interval varies from 775 ms to 800 ms corresponding to an HR between 75 and 77.5 bpm. One LQT1 and two LQT2 subjects were excluded from the analysis because they did not have any cardiac beat in this HR range.
Based on median beat measurements from the Holter ECG recording,
As previously described the LQT2 patients present very low and flat T-waves. In our study, these results are confirmed by revealing different values between the two groups in terms of T-wave slopes and MV values. LQT1 have taller T-wave: 0.33±0.14 mV vs. 0.11±0.16 mV for LQT1 and LQT2 patients respectively (p<0.05) and steeper slopes. Right and left slopes of T-wave were almost three fold higher in LQT1 patients: 8.0±3.0 vs. 3.6±1.7 μV/ms in LQT1 and LQT2 for the left slope and 12.7±5.1 vs. 4.3±2.6 μV/ms for the right slope.
The right slope of T-wave is the best parameter for the discrimination between LQT1 and LQT2 patients, with AUC=0.95 and 87.5% specificity and 87% sensitivity. T-wave magnitude in lead V5 was the best classical ECG parameter, with AUC=0.85 and 73.9% specificity and 78.3% sensitivity.
Other vectorial parameters such as T-loop roundness and T-loop planarity were significantly higher in LQT2 than in LQT1, revealing a profound difference in the repolarization process affecting the overall orientation and electrical activity within the myocardium between the two mutations.
The morphology of the T-wave is different between males and females. We investigated the role of gender in the two LQTS populations. These results are summarized in
Statistical Model for the Discrimination Between LQT1 and LQT2 Patients: Multivariate Analysis
First, we designed a model based on clinical parameters (QT, and RR), then added all the scalar ECG parameters and finally introduced the vectorial parameters.
Clinical Model
This model has been implemented as a reference. It describes the level of discriminant power obtained when using our statistical strategy and when only QT and RR intervals are used, the currently used clinical electrocardiographic parameters. As expected, this model performed poorly (see.
Scalar Model
The parameters included in the model were: RR, QTpeak, TpTe, T-wave magnitude and QT from lead II and V5. Using the best subset model based on the AIC criteria, a model using 2 parameters (Tmag and TpTe from lead V5) was selected. The model provided a good discrimination of the two mutations: 95.8% specificity and 87.0% sensitivity.
Vectorial Model
The design of the computerized model was based on the following list of repolarization parameters: parameters from scalar model and λ2/λ1, λ3, right and left slope of T-wave from lead ev1, MV, Loop QT, Loop QTapex and 3 parameters from T-wave modeling with Hill equation. The best model with three parameters relied on the following ones: αR, QTpeak from Lead II and QT interval from the T-loop could discriminate the two groups with 95.8% specificity and 95.7% sensitivity, for LQT2 prediction.
Evaluation in a Wider HR Range
We evaluated the stability of our predictive models for larger HR intervals to partially assess the validity of the models. We fitted the current models with data from larger ranges of heart rate. We defined three HR ranges:
The ROC curves for the scalar model are presented in
The information contained in the morphology of the T-wave from both scalar ECG and vectorcardiographic signals are relevant in discriminating the patients according to their type of mutations (HERG vs. KvLQT1). In particular, we observed that HERG patients present lower values of MV (related to T-wave magnitude in all scalar leads), this can also be assessed using the left and right slopes of T-wave (flatter T-wave). T, and the TpTe interval is significantly longer than in LQT1 patients. These results are consistent with previous work.
T-wave magnitude and MV show very different heart rate dependency between mutations. In LQT2, there is almost no HR dependency, while in LQT1 T-wave magnitude has an inverted HR dependency, T-wave amplitude increases when HR decreases.
We believe that these differences between the two LQTS groups are linked to the abnormal kinetics of the delayed rectifier potassium currents of the myocardium cells. It is the efflux of potassium ions that produces the T-wave of the electrocardiogram. Agents that delay or prevent this movement of potassium will modify the appearance of the T-wave. More precisely, IKr is a current moving potassium ions out of the cardiac cells during the plateau phase of the action potential. One may expect to see changes on the surface ECG recordings in the late part of the T-wave.
Clinical Usefulness of Repolarization Morphology
Genotype identification by ECG is useful for stratifying molecular genetic studies. With 5 disease genes and 170 mutations already identified, it is very costly and time consuming to screen all known genes and mutational sites, limiting the application of genetic studies. With a typical ECG pattern, the suspected gene can be the initial target for testing, with a higher likelihood of rapid identification of the mutation. Such a strategy will significantly reduce time and costs, allowing more families to be genotyped and enhancing genotype-phenotype correlation studies. Furthermore, if therapeutic interventions based on specific genotype are shown to be effective, genotype identification by ECG could be helpful for therapeutic decision-making. These findings are applicable only to patients and families with an established clinical diagnosis of LQTS, from LQT1 and LQT2 mutations.
The multivariate analysis demonstrated that the use of very simple measurements such a the amplitude of the T-wave, the length of specific intervals such as TpTe provides excellent discrimination between LQT1 and LQT2 patients. Thus, the need for more sensitive quantifiers does not seem to bring a significant improvement (˜+8% sensitivity). It is noteworthy that in case of smaller inhibition of IKr current such as in drug induced LQTS such parameters may become more useful.
We evidenced very different behavior of T-wave magnitude between the two types of mutation across heart rate. Such observation may be very interesting and should be investigated further. One may imagine designing a simple test in which the T-wave magnitude and RR relationship may be modeled in order to identify a loss of relationship, according to our observation in LQT2 patients such loss seems to characterize patients with IKr-inhibition. Such test may be applicable to the congenital and drug-induced LQTS.
Identifying IKr or IKs-related abnormalities of the surface ECGs may find broader application than discriminating the type of congenital LQTS, in a more general population estimating the propensity of an individual to the drug-induced LQTS is a relevant question. The drug-induced QT prolongation and the development of lethal arrhythmias of an individual due to a compound affecting the kinetics of the ion currents of the myocardial cells, depends on numerous factors including the individual genetic make-up. Thus, an individual predisposition to certain compound could be assessed based on ECG factors or the development of electrocardiographic signs for a predisposition to cardiac events on certain drug could be evaluated. Although this test was not able to screen normal people from those with mutations, it is possible that similar models may be developed to screen normal patients from those with LQT1 and LQT2 in addition to being able to differentiate between LQT1 carriers and LQT2 carriers.
We demonstrated that LQT1 and LQT2 patients have specific phenotypic expressions of their mutation on the surface ECG and these expressions can be quantified using simple parameters to discriminate these two types of mutations. We evidenced different behaviors of T-wave magnitude between mutations across heart rate. Finally, we developed a binary logistic model computing the probability of a LQTS patient to carry an HERG mutation based on simple ECG parameters.
A data input 64 is coupled to the processor 62 and configured to provide the processor with ECG data. An ECG capture device 66 may optionally be coupled to the data input 64 to enable the live capture of ECG data. Examples of ECG capture devices include, but are not limited to, a twelve-lead ECG device, an eight-lead ECG device, a two lead ECG device, a Holter device, a bipolar ECG device, and a uni-polar ECG device. Similarly, a database 68 may optionally be coupled to the data input 64 to provide previously captured ECG signal data to the processor. Database 68 can be as simple as a memory device holding raw data or formatted files, or database 68 can be a complex relational database. Depending on the embodiment, none, one, or multiple databases 68 and/or ECG capture devices 66 may be coupled to the data input 64. The ECG capture device 66 may be coupled to the data input 64 by a wired connection, an optical connection, or by a wireless connection. Suitable examples of wireless connections may include, but are not limited to, RF connections using an 802.11x protocol or the Bluetooth® protocol. The ECG capture device 66 may be configured to transmit data to the data input 64 only during times which do not interfere with data measurement times of the ECG capture device 66. If interference between wireless transmission and the measurements being taken is not an issue, then transmission can occur at any desired time. Furthermore, in embodiments having a database 68, the processor 62 may be coupled to the database 68 for storing results or accessing data by bypassing the data input 64.
The system 60 also has a user interface 70 which may be coupled to either the processor 62 and/or the data input 64. The user interface 70 can be configured to display the ECG signal data, a statistical plot of LQT1 vs LQT2 determination based on the ECG signal data, and/or an overall determination of LQT1 vs LQT2 status with or without confidence indicators. The user interface 70 may also be configured to allow a user to select ECG signals from a database 68 coupled to the data input 64, or to start and stop collecting data from an ECG capture device 66 which is coupled to the data input 64.
LQT2 prediction: log(odds)=log(p2/(1/p2))=−7.19−15.70(Tmag)+0.10(TpTe) (Example 1)
LQT2 prediction: log(odds)=log(p2/(1−p2))=11.7−1.35(RS T)−0.32QTp+0.25(T−loop QT) (Example 2)
LQT2 prediction: log(odds)=11.73−1352(αR on ev1)+0.2480(Loop Offset QT)−0.3163(QTpeak from Lead II) (Example 3)
It should be understood that the LQT2 prediction examples above could be converted to LQT1 probabilities instead.
Other embodiments may include a filtering 50 step which may include statistical combinations of multiple beats from the ECG signals. As a non-limiting example, a median beat may be created from a number of consecutive beats from each lead. In some embodiments, one or more leading beats may be discarded. In other embodiments, one or more trailing beats may be discarded. In further embodiments, only beats with a stable heart rate may be taken into account. An example of a suitable definition of beats with a stable heart rate is when the heart rate for a given beat varies less than ten percent in beats of the previous two minutes. In other embodiments other percentages, time-frames, and definitions of a stable heart rate may be used without deviating from the scope of the claimed invention.
The advantages of a method and system to differentiate between LQT1 and LQT2 carriers have been discussed herein. Embodiments discussed have been described by way of example in this specification. It will be apparent to those skilled in the art that the forgoing detailed disclosure is intended to be presented by way of example only, and is not limiting. Various alterations, improvements, and modifications will occur and are intended to those skilled in the art, though not expressly stated herein. These alterations, improvements, and modifications are intended to be suggested hereby, and are within the spirit and the scope of the claimed invention. Additionally, the recited order of processing elements or sequences, or the use of numbers, letters, or other designations therefore, is not intended to limit the claims to any order, except as may be specified in the claims. Accordingly, the invention is limited only by the following claims and equivalents thereto.
This application claims priority to U.S. provisional patent application 60/780,854, entitled “Quantitative Repolarization Patterns Identifying KvLQT1 and HERG Mutation in Patients with Long QT Syndrome”, and filed Mar. 10, 2006. The entire 60/780,854 patent application is hereby officially incorporated by reference.
This invention was made with government support under grant HL068226 and HL033843 awarded by NIH. The government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
20080033313 A1 | Feb 2008 | US |
Number | Date | Country | |
---|---|---|---|
60780854 | Mar 2006 | US |