Patent Description:
Cardiac fibrillation is one of the leading causes of morbidity and mortality in the developed world, where atrial fibrillation (AF) is the most common sustained arrhythmia (<NPL>), affecting approximately <NUM> million patients in the United States, and its prevalence increases with age (<NPL>). Although AF is not considered a life-threating arrhythmia, it affects the circulatory system of patients and their general health and quality of life; further, it poses a significant burden on the healthcare system (<NPL>): e.g., AF is a significant risk factor for stroke, with about <NUM>% of strokes occurring in people with AF (<NPL>). Unfortunately, currently available drug treatments for AF are less than satisfactory. An expert panel convened by the Heart, Lung and Blood Institute to resolve this issue determined that "AF detection by using emerging methods and technologies" constitutes an important future direction in preventing AF (<NPL>). To date, there is no effective clinical method for early prediction and precise identification of AF events.

Interestingly, while persistent AF is associated with an increased HRV on average (<NPL>), reduced HRV indices are observed just prior to arrhythmogenic events (<NPL>). Although the correlation between changes in heart beat complexity and prevalence of AF has been known for more than three decades (<NPL>; <NPL>; <NPL>), a method to detect HRV changes on line and in short period before AF event occurs does not exist.

<CIT> describes a method for medical premonitory event estimation includes one or more processors to perform operations. This method comprises: acquiring a first set of physiological information of a subject, and a second set of physiological information of the subject received during a second period of time; calculating a first and second risk scores associated with estimating a risk of a potential cardiac arrhythmia event for the subject based on applying the first and second sets of physiological information to one or more machine learning classifier models; providing at least the first and second risk scores associated with the potential cardiac arrhythmia event as a time changing series of risk scores; and classifying the first and second risk scores associated with estimating the risk of the potential cardiac arrhythmia event for the subject based on the one or more thresholds.

<CIT> describes a system for cardiac monitoring, that may be used to sense the presence of an intermittent arrhythmia in an individual. <CIT> discloses an apparatus and method for ambulatory, real-time detection of Atrial Fibrillation (AF) providing an overall accuracy that refers to detection of AF, irrespective of the duration of AF and beat-to-beat classification.

There is provided, in an embodiment, a method comprising receiving a plurality of data segments each representing human cardiac activity; calculating, from each of said data segments, a first value representing a distribution of a heart rate variability (HRV) parameter; acquiring a signal representing cardiac activity in a patient; deriving a second value representing distribution of an HRV parameter in said signal; and indicating an atrial fibrillation (AF) event in said patient, when a difference between said first and second values is above a predetermined threshold.

There is also provided in an embodient, a system comprising at least one hardware processor; and a non-transitory computer-readable storage medium having stored thereon program instructions, the program instructions executable by the at least one hardware processor to: receive a plurality of data segments each representing human cardiac activity, calculate, from each of said data segments, a first value representing a distribution of a heart rate variability (HRV) parameter, acquire a signal representing cardiac activity in a patient, derive a second value representing distribution of an HRV parameter in said signal, and indicate an atrial fibrillation (AF) event in said patient, when a difference between said first and second values is above a predetermined threshold.

There is further provided, in an embodiment, a computer program product comprising a non-transitory computer-readable storage medium having program instructions embodied therewith, the program instructions executable by at least one hardware processor to: receive a plurality of data segments each representing human cardiac activity; calculate, from each of said data segments, a first value representing a distribution of a heart rate variability (HRV) parameter; acquire a signal representing cardiac activity in a patient; derive a second value representing distribution of an HRV parameter in said signal; and indicate an atrial fibrillation (AF) event in said patient, when a difference between said first and second values is above a predetermined threshold.

In some embodiments, said plurality of data segments represent cardiac activity in a plurality of subjects.

In some embodiments, at least some of said data segments are acquired from said patient.

In some embodiments, said HRV parameter is an entropy parameter.

In some embodiments, said entropy parameter is modified Shannon entropy.

In some embodiments, said calculating comprises calculating a distribution of said HRV parameter in all said data segments.

In some embodiments, each of said data segments is acquired over a period of between <NUM> and <NUM> minutes.

In some embodiments, said signal is acquired over a period of between <NUM> and <NUM> minutes.

In some embodiments, the method further comprises training, and the program instructions are further executable to train, at a training stage, a machine learning classifier on a training set comprising: (i) said plurality of data segment, wherein at least some of said segments comprise cardiac activity reflecting an arrythmia event, and (ii) labels indicating the occurrence of an arrythmia event in said cardiac activity.

In some embodiments, said data segments are labelled with said labels.

In some embodiments, said indicating is performed by applying, at an inference stage, said machine learning classifier to said signal, to indicate an AF event in said patient.

In some embodiments, said indicating comprises at least one of detecting an ongoing AF event, and predicting an upcoming AF event.

In some embodiments, said upcoming AF event occurs between <NUM> and <NUM> minutes after said indicating.

In some embodiments, said detecting is based, at least in part, on calculating statistical estimators of said distribution.

In some embodiments, said statistical estimators comprise at least standard deviation and deviation from normal distribution.

In some embodiments, said deviation from normal distribution is based on the Kolmogorov-Smirnov statistic.

In some embodiments, the method further comprises applying, and the program instructions are further executable to apply, a clustering algorithm to said plurality of data segments, to cluster at least some of said data segments based, at least in part, on said HRV parameter. In some embodiments, said clustering algorithm is a spectral clustering algorithm.

A system and method are disclosed herein for detecting and/or predicting atrial fibrillation (AF) or other irregular regularity. In some embodiments, the present invention provides for real-time prediction or detection of AF or other irregular regularity by employing an algorithm based on pattern classification of the distribution of a modified Shannon entropy of beating intervals of heart beats intervals. The present method and system produce a real-time alarm when certain changes occur in the modified Shannon entropy of a patient.

The present invention is based on the finding that persistent AF is associated with an increased HRV on average (<NPL>). However, reduced HRV indices are observed just prior to events (<NPL>). Therefore, an atypical decrease in HRV might be predictive for an upcoming episode of AF, and thus a useful marker to be used in the detection, prevention, research and treatment of AF.

The present invention provides clinical and diagnostic tools for early prediction and detection of AF events by examining changes in HRV. Real-time analysis of HRV of a patient, based on short segments of heartbeat intervals, according to the present method may help to detect and prevent AF events, thereby reducing the probability of associated events, such as strokes. In some embodiments, the disclosed invention allows for prediction of AF events shortly (e.g., <NUM>. ) before they occur. Early prediction of AF events thus will allow to determine which hemodynamic, mechanical or electrical parameters change just before each episode, and help during ablation procedure.

The present invention can be implemented, for example, as an early detection and prediction system for atrial fibrillation based on continuous heart rate monitoring, as well as for early prediction of arrhythmia episodes in general. In some embodiments, heart beat segments can be acquired from a mobile device, a wearable device, a device that is used to monitor or diagnose medical conditions by medical experts, and/or based on a smart implantable cardiac defibrillator. Because the disclosed method provides for quantifying the HRV accurately based on relatively short segments of heart beats, it can be implemented based on, e.g., a remote ECG acquiring system. Thus, the present method can provide for efficient and cost-effective local and/or remote monitoring of multiple AF patients.

In some embodiments, the present disclosure can be implemented, for example, as a diagnostic system used for assessing the underlying etiology and/or the arrythmia profile of each patient's AF events, allowing personalized treatment.

In some embodiments, methods for predicting and/or detecting AF according to the present disclosure comprise a first step of acquiring and/or receiving one or more prerecorded baseline segments of cardiac activity comprising multiple heart beat intervals of a subject. The baseline segments comprise a plurality of heart beat intervals taken, e.g., at rest, when no arrythmia event was in progress. In some embodiments, the recording of the baseline signal is performed over a period of, e.g., <NUM> to <NUM> minutes. A modified Shannon entropy distribution and/or other HRV indices are may then be calculated for the baseline signal. In some embodiments, the present method further comprises calculating a difference in Shannon entropies between every two or more consecutive heart beat intervals of the baseline signal, wherein the calculation of the baseline distribution is based on the calculated differences in Shannon entropies.

In some embodiments, the modified Shannon Entropy is used to separate patients into groups by, e.g., Spectral Clustering.

In real time, a second, shorter, segment of the patient cardiac activity measurement is acquired (the "real-time segment" ), wherein a modified Shannon entropy distribution and/or other HRV indices are calculated for the real-time segment as well. The distributions of the baseline and real-time segments are compared to one another. In some embodiments, when a difference in distribution width between the baseline and second segments is above a predetermined threshold, an AF event is predicted and/or detected. In some embodiments, the recording of the real-time segment is performed over a period of <NUM> seconds to <NUM> minutes.

In some embodiments, the present disclosure may be implemented as a machine leanring classifier trained on a training set comprising a plurality of cardiac activity segments acquired from a population of subjects, wherein at an inference stage the classifier may be configured to predict and/or detect AF event in a real-time segment acquired from a monitored subject.

Accordingly, in some embodiments, a machine leanring classifier may be trained, e.g., on a training set comprising a plurality of recorded baseline segments measurements from, e.g., a group of subjects. In some embodiments, the baseline segments may comprise a plurality of heart beat intervals of each subject taken, e.g., at rest, when no arrythmia event was in progress. In some embodiments, each segment may comprise, e.g., <NUM> to <NUM> minutes of recorded cardiac activity. In some embodiments, an existing database of cardiac activity recordings form a plurality of subjects may comprise at least a portion of the training set. In some embodiments, a modified Shannon entropy distribution and/or other HRV indices are may be derived for each segment within the training set. In some embodiments, modified Shannon entropy may be derived, e.g., by calculating a difference in Shannon entropies between every two or more consecutive heart beat intervals of the baseline signal.

In some embodiments, the training set may be labelled with appropriate labels indicating the existence of an AF event within each segment. In some embodiments, the machine learning classifier may be trained on the training set to classify a segment of cardiac activity as indicative of a concurrent and/or imminent (e.g., within a time horizon of <NUM>-<NUM> mins. ) AF event. In some embodiments, a trained classifier of the present disclosure may be used to detect and or predict an AF event, based, at least in part, on a segment of cardiac activity measurement acquired in real time from a subject, wherein the segment may comprise a recording of cardiac activity over a period of <NUM> seconds to <NUM> minutes.

With reference to <FIG>, in some embodiments, an exemplary system <NUM> in accordance with the present disclosure comprises processing unit <NUM>, comprising at least one hardware processor. Processing unit <NUM> is communicatively connected to a non-transient computer readable storage device <NUM>. Processing unit <NUM> is configured to automatically control the operation of system <NUM>, based on one or more applications, sets of software instructions, and algorithms stored storage device <NUM>. Storage device <NUM> may be configured to have stored thereon digital representations of cardiac signals and/or measurements. In some embodiments, processing unit <NUM> is configured to receive and process one or more cardiac signals and/or measurements collected by one or more cardiac signal sensors <NUM>. System <NUM> may further comprise a communication module <NUM>, which may be configured to connect system <NUM> through a communication network, such as the Internet, a local area network, a wide area network, and/or a wireless network. For example, system <NUM> may be configured to receive cardiac signals from a remote acquisition unit over a network via communication module <NUM>. System <NUM> described herein is only an exemplary embodiment, and may have more or fewer components than shown, may combine two or more components, or a may have a different configuration or arrangement of the components. The various components of system <NUM> may be implemented in hardware, software, or a combination of both hardware and software. According to various other embodiments, processing unit <NUM> or processing tasks performed thereby may be implemented by a handheld or worn computing device such as, but not limited to, a smart phone, a tablet computer, a notepad computer, smart watch, and the like. In addition, aspects of the present system which can be implemented by computer program instructions, may be executed on a general-purpose computer, a special-purpose computer, or other programmable data processing apparatus.

<FIG> is a flowchart illustrating the functional steps of an exemplary prediction and/or detection machine leanring classifier <NUM>, according to some embodiments of the present disclosure. Classifier <NUM> is configured to calculate HRV over relatively short real-time segments of heart beat, based on pattern classification of distribution of the modified Shannon entropy of beating intervals. Classifier <NUM> works in two stages: a learning stage (steps <NUM>-<NUM>) and an inference stage (steps <NUM>-<NUM>).

At a learning stage, starting at <NUM>, a training set of baseline cardiac activity segments of a single or a plurality of subjects may be acquired or received. The baseline segments can be acquired using an ECG and/or similar acquiring methodology from the subjects at rest. In some embodiments, the baseline segments comprise a segment of ECG and/or other bio-signal measurement which can be used to calculate the heart beating intervals (e.g., saturated oxygen, pulse detector, blood pressure, etc.). In some embodiments, each segment is recorded over a time period that is sufficient to create a reliable distribution curve for statistics, for example, <NUM> minutes. In some embodiments, the recording of the baseline segments is performed over a period of <NUM> to <NUM> minutes where no arrhythmia events are present.

At <NUM>, modified Shannon entropy is calculated with respect to each segment in the training set. For example, the difference in Shannon entropy of a pair in each three consecutive beats (or more beats with N-<NUM> combination) is calculated (to reduce noise in the signal). The absolute value between each pair of beats is calculated, to provide a modified Shannon entropy value.

At <NUM>, the distribution of modified Shannon entropy values is calculated for the entire segment. In some embodiments, the distribution is normalized for the average modified Shannon entropy, and the shape of this probability density function is calculated. In some embodiments, the distribution of the modified Shannon entropy of the segment is represented by estimation of statistical parameters.

At <NUM>, a machine learning classifier is trained on the training set acquired and processed in steps <NUM>-<NUM>.

At an inference stage <NUM>, a real-time sample segment is taken over a relatively short period of time, e.g., <NUM>. In some embodiments, the real-time recording of is performed over a period of <NUM> seconds to <NUM> minutes.

At <NUM>, classifier <NUM> is applied to the real-time segment, to predict and/or detect an AF event, based, at least in part, on a distribution parameter of Shannon entropy and/or other HRV indices in the real-time segment.

At <NUM>, in case of significant differences from the baseline distribution, the real-time measurements and processing may be repeated, e.g., between <NUM> and <NUM> times.

In some embodiments, the present disclosure employs, at least in part, the following definitions and assumptions for predicting and/or detecting AF events. For example, in the case of a finite set of consecutive beat intervals from a certain subject (time intervals between consecutive heart beats): <MAT> The present modified Shannon entropy of order <NUM> may be defined as <MAT>.

And for a general order M where <MAT>, M > <NUM>: <MAT>.

It is further assumed that there exists a probability density function parametric family <MAT>.

It should be noted that θ might be a multi-dimensional parameter set, for which distribution there will be used the notion Y(θ). Satisfying that for every possible set D of consecutive beat intervals from the same subject exists a θ ∈ Θ, such that <MAT>.

For the probability density function fX(x, θ), there is defined a "distribution width" function: <MAT> where W(fX, (x, θ)) might be thought of as a generalized variance.

At the learning stage, the prediction algorithm works to estimate the characteristic θ for a certain subject in an arrhythmia-free time interval. For example, in the case of a learning set <MAT> satisfying <MAT> there is calculated an estimation for θ, for which the notion θ̂ will be used. For θ̂ there is calculated the estimated distribution width <MAT>.

In the real-time calculation stage, the prediction algorithm defines a "narrow" distribution as: <MAT> Where, T(w) is a function <MAT> used in threshold determination. One might think of T(w) as a function determining how low a value of W should be to satisfy P(w ≤ W) < α. The real-time calculation set is then: <MAT> satisfying <MAT>.

Accordingly, with every new beat, the prediction algorithm estimates the current Ŵ for n iterations. In some embodiments, n is a predetermined constant number, which can be chosen based on either heuristic methods, or through statistical analysis. If n or more consecutive real-time calculation sets are considered "narrow," it is considered as an "alarm. " The implementation of the present prediction algorithm requires determination of an estimator for the distribution width W and a threshold estimating function T.

In some embodiments, the present disclosure can be implemented, for example, as a diagnostic system used for assessing the underlying etiology and/or the arrythmia profile of each patient's AF events, allowing personalized treatment. In some embodiments, the modified Shannon Entropy is used to separate patients into groups by Spectral Clustering (see, e.g., <NPL>). In some embodiments, the spectral clustering is based on a Laplace Kernel distance metric: <MAT> where N is the length of the patient's entropy signal.

In some embodiments, the threshold estimating function, T, may be implemented using classic statistical estimators. In some embodiments, threshold estimating function, T, may be implemented using a machine learning algorithm.

In some embodiments, an exemplary detection algorithm according to some embodiments is based on the assumptions and definitions above, with the following additional definitions. For a set of N observations <MAT>, the empirical distribution function is defined as: <MAT>.

The estimators for the expectancy standard deviation are defined as the maximum likelihood estimators: <MAT> The estimated normal distribution is defined as: <MAT> where Φ is the cumulative distribution function of a normal standard random variable and the deviation from normality is defined as the Kolmogorov-Smirnov statistic: <MAT>.

Although under the present algorithm assumptions described above, the probability density function of the modified entropy exists, even if it is assumed that beat intervals are normally distributed, the use of a sum of absolute values in the modified entropy calculation might result in a different parametric family for the modified entropy. Moreover, the calculation of consecutive modified entropies used overlapping beat intervals, thus, one cannot assume Ŝi and Ŝi+<NUM> to be statistically independent. Thus, under different assumptions, different parametric families of distribution functions may be used.

This theoretical result was clearly observed in numerical attempts to fit a normal distribution density function to sets of modified entropy samples, even when calculated from signals taken from the healthy subjects. <FIG> illustrates a density function fitting calculated and plotted by a Matlab Distribution fitting toolbox, in which a discrepancy between the histogram and the best fitting is evident. Accordingly, instead of finding a parametric family and truly calculating the distribution width from the estimated probability density function, in some embodiments, the algorithm is configured to try to estimate the distribution width directly from the samples, some examples for such an estimator are described below. In some embodiments the algorithm directly assesses the need for an alarm from the beat intervals sequence. Nevertheless, in some embodiments, the full process of distribution parametric family determination, parameter estimation and threshold estimation is utilized.

For a set of N modified entropy samples <MAT>, first calculate the empirical cumulative distribution function: <MAT> where <NUM>A is the indicator function of event A. Then, define a tail <MAT>, and find the samples Ŝlow,Ŝhigh satisfying: <MAT> <MAT>.

The estimator to the distribution width is then: <MAT>.

The use of empirical distribution functions yielded an interesting observation. <FIG> shows multiple empirical distribution functions of <NUM>-minute sets of beat intervals from the same subject on the same figure. Difference can be noted between panel A, plotted using a recording from a subject with occasional AF episodes, and panel B, plotted using a recording from a healthy subject. In the healthy subject, it can qualitatively be seen that although some variability between the distribution of different sets exists, the variability is low and there are no distinct groups of curves. Conversely, in panel A, representing the AF patient's records, larger variability can be seen with <NUM> distinguishable groups of distribution width. Another, more subtle difference is the presence of some very narrow curves in the AF patient's results, which are absent in the healthy patient's results. These differences suggest that the empirical distribution width might hold the information needed to differentiate patients with and without AF.

A natural choice for a parameter to estimate distribution width is variance. The estimation of variance was made by the unbiased version of the maximum likelihood estimator: <MAT>.

To reduce the noise, a moving average filter may be used, which produce a smoother signal. <FIG> shows a time plot with the σ̂<NUM> estimator, while <FIG> shows the same time interval with the filtered estimator, wherein annotations of AF onset are marked by Xs. Higher values during arrhythmia episodes can be seen, with and the low values often preceding their onset.

An example of the threshold determination from a learning set <MAT>, comprised of beat intervals from a patient, may be:.

Using this threshold, the preliminary definition for an alarm is:<MAT>.

<FIG> shows a plot of the σ̂<NUM> estimator after a running average filter in a representative time interval from an AF subject. Annotations of AF onset are marked by Xs, and alarms are marked by Xs of lighter shade.

<FIG> shows the same interval from the previous figures, now with the alarms (using α = <MAT>) marked on the plot. <FIG> shows a similar plot of another time interval from the same subject.

As a point of comparison, <FIG> shows a similar plot from a healthy subject. Even though some fluctuating does occur in healthy subjects, extremely low values are rare, and therefore in most of healthy subjects, no alarms are produced.

Using the Laplace Kernel distance metric for the Spectral Clustering algorithm, we clustered a patients cohort of both AF and healthy patients into <NUM> groups. Than we calculated the mean Multi-Scale Entropy (see, <NPL>) measure for each group. In <FIG>, one can see the correlation between the <NUM> groups from a controlled rabbits' experiment (<FIG>, see, e.g., Rosenberg, A. (<NUM>), "Non-Invasive In-Vivo Analysis of Intrinsic Clock-Like Pacemaker Mechanisms: Decoupling Neural Input Using Heart Rate Variability Measurements" (master's thesis)) to that of the clustering (<FIG>). On further investigation, the clustering was found to group patients with similar pattern of arrhythmia.

In some embodiments, the clustering method herein is used to classify patients to groups to assist the personalization of medical treatment.

Atrial fibrillation is classically defined as ". a supraventricular arrhythmia characterized electrocardiographically by low-amplitude baseline oscillations (fibrillatory or f waves) and an irregularly irregular ventricular rhythm. " (see, <NPL>). Thus, in some embodiments, an algorithm identifying irregular irregularity of heart rate is used to detect atrial fibrillation episodes.

Heart rate irregularity show a strong correlation with the standard deviation of the modified Shannon entropy, but does not tell regular irregularity and irregular irregularity apart. The normality measure described above tell regular irregularity and irregular irregularity apart but is not specific for irregular irregularity and may give similar results in irregular irregularity and in normal sinus rhythm with complex dynamics. Thus, in some embodiments, both measures are used together to detect AF episodes. In <FIG> one may see that although both the standard deviation (denoted "var" in the legend) and the normality measure (denoted "ks" in the legend) can rise separately on non-AF times, they rise together only during the AF episodes.

Such a propagated signal may take any of a variety of forms, including, but not limited to, electromagnetic, optical, or any suitable combination thereof.

Computer program code for carrying out operations for aspects of the present invention may be written in any combination of one or more programming languages, including an object-oriented programming language such as Java, Smalltalk, C++ or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages.

Aspects of the present invention are described herein with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. These computer program instructions may be provided to a hardware processor of a general-purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.

The flowcharts and block diagrams in the Figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present invention.

Claim 1:
A system comprising: an electrocardiogram (ECG) acquiring system (<NUM>); at least one hardware processor (<NUM>); and a non-transitory computer-readable storage medium (<NUM>) having stored thereon program instructions, the program instructions executable by the at least one hardware processor to receive (<NUM>), from the ECG acquiring system, data segments representing human cardiac activity, wherein the program instructions are characterized by:
recording a baseline signal, comprising a plurality of data segments over a first period;
calculating (<NUM>) a difference in Shannon entropies between every two or more consecutive heartbeat intervals of the baseline signal, to obtain from each of said data segments a first value representing distribution of modified Shannon entropy ;
acquiring a real-time segment (<NUM>) representing cardiac activity in a patient over a second period, shorter than the first period;
calculating a second value representing distribution of modified Shannon entropy of the real-time segment;
calculating a difference between said first and second values; and
indicating an atrial fibrillation (AF) event in said patient when the calculated difference is above a predetermined threshold.