A method for detecting myocardial abnormality by using a current-ratio map calculated from an exercise-induced magnetocardiogram

The magnetocardiogram

The magnetic field produced by the heart's electrical activity is called the magnetocardiogram (MCG). The first 20 years of MCG research established most of the concepts, instrumentation, and computational algorithms in the field. Additional insights into fundamental mechanisms of biomagnetism were gained by studying isolated hearts or even isolated pieces of cardiac tissue. Much effort has gone into calculating the MCG using computer models, including solving the inverse problem of deducing the bioelectric sources from biomagnetic measurements. Recently, most magnetocardiographic research has focused on clinical applications, driven in part by new technologies to measure weak biomagnetic fields.

A movable unshielded magnetocardiography system

Magnetocardiography (MCG), which uses high-sensitivity magnetometers to record magnetic field signals generated by electrical activity in the heart, is a noninvasive method for evaluating heart diseases such as arrhythmia and ischemia. The MCG measurements usually require the participant keeping still in a magnetically shielded room due to the immovable sensor and noisy external environments. These requirements limit MCG applications, such as exercise MCG tests and long-term MCG observations, which are useful for early detections of heart diseases. Here, we introduce a movable MCG system that can clearly record MCG signals of freely behaving participants in an unshielded environment. On the basis of optically pumped magnetometers with a sensitivity of 140 fT/Hz1/2, we successfully demonstrated the resting MCG and the exercise MCG tests. Our method is promising to realize a practical movable multichannel unshielded MCG system that nearly sets no limits to participants and brings another kind of insight into the medical diagnosis of heart disease.

Assessing heart disease using a novel magnetocardiography device

The aim of this paper is to present the use of a portable, unshielded magnetocardiograph (MCG) and identify key characteristics of MCG scans that could be used in future studies to identify parameters that are sensitive to cardiac pathology. We recruited 50 patients with confirmed myocardial infarction (MI) within the past 12 weeks and 46 volunteers with no history of cardiac disease. A set of 38 parameters were extracted from MCG features including both signals from the sensor array and from magnetic images obtained from the device and principal component analysis was used to concentrate the information contained in these parameters into uncorrelated predictors. Linear fits of these parameters were then used to examine the ability of MCG to distinguish between sub-groups of patients. In the fist instance, the primary aim of this study was to ensure that MCG has a basic ability to separate a highly polarised patient group (young controls from post infarction patients) and to identify parameters that could be used in future studies to build a formal diagnostic tool kit. Parameters that parameterised left ventricular ejection fraction (LVEF) were identified and an example is presented to show differential low and high ejection fractions.

Clinical application of magnetocardiography

Magnetocardiography is a noninvasive contactless method to measure the magnetic field generated by the same ionic currents that create the electrocardiogram. The time course of magnetocardiographic and electrocardiographic signals are similar. However, compared with surface potential recordings, multichannel magnetocardiographic mapping (MMCG) is a faster and contactless method for 3D imaging and localization of cardiac electrophysiologic phenomena with higher spatial and temporal resolution. For more than a decade, MMCG has been mostly confined to magnetically shielded rooms and considered to be at most an interesting matter for research activity. Nevertheless, an increasing number of papers have documented that magnetocardiography can also be useful to improve diagnostic accuracy. Most recently, the development of standardized instrumentations for unshielded MMCG, and its ease of use and reliability even in emergency rooms has triggered a new interest from clinicians for magnetocardiography, leading to several new installations of unshielded systems worldwide. In this review, clinical applications of magnetocardiography are summarized, focusing on major milestones, recent results of multicenter clinical trials and indicators of future developments.

Subtraction magnetocardiogram for detecting coronary heart disease

BACKGROUND

A large-scale magnetocardiogram (MCG) database was produced, and standard MCG waveforms of healthy patients were calculated by using this database. It was clarified that the standard MCG waveforms are formed with the same shape and current distribution in healthy patients. A new subtraction method for detecting abnormal ST-T waveforms in coronary heart disease (CHD) patients by using the standard MCG waveform was developed.

METHODS

We used MCGs of 56 CHD patients (63 ± 3 years old) and 101 age-matched normal control patients (65 ± 5 years old). To construct a subtracted ST-T waveform, we used standard MCG waveforms produced from 464 normal MCGs (male: 268, female: 196). The standard MCG waveforms were subtracted from each subject's measured MCGs, which were shortened or lengthened and normalized to adjust to the data length and magnitude of the standard waveform. We evaluated the maximum amplitude and maximum current-arrow magnitude of the subtracted ST-T waveform.

RESULTS

The maximum magnetic field, maximum magnitude of current arrows, and maximum magnitude of total current vector increased according to the number of coronary artery lesions. The sensitivity and specificity of detecting CHD and normal control patients were 74.6% and 84.1%, respectively.

CONCLUSIONS

The subtraction MCG method can be used to detect CHD with high accuracy, namely, sensitivity of 74.6% and specificity of 84.1% (in the case of maximum amplitude of total current vector). Furthermore, the subtraction MCG magnitude and its current distribution can reflect the expanse of the ischemic lesion area and the progress from ischemia to myocardial infarction.

Repolarization spatial-time current abnormalities in patients with coronary heart disease

BACKGROUND

Magnetocardiography (MCG) is a new technique for visualizing a current distribution in the myocardium. In recent years, current distribution parameters (CDPs) have been developed based on the distribution. The CDPs reflect spatial-time current abnormalities in patients with coronary heart disease (CHD). However, the criteria and scoring method of the abnormalities using CDPs are still controversial.

METHOD

We measured MCG signals for 101 normal controls and 56 CHD patients (single-, double-, and triple-vessel diseases) using a MCG system. The CDPs (maximum current vector [MCV], total current vector [TCV], current integral map, and current rotation) during ventricular repolarization were analyzed. To evaluate the CDPs that are effective in distinguishing between normal controls and CHD patients, the areas under the receiver operating characteristic curve (A(z)) are calculated. Furthermore, the total scores ("0" to "4") of four CDPs with high A(z) values are also calculated.

RESULTS

MCV and TCV angles at the T-wave peak had the highest A(z) value. Furthermore, TCV angular differences between the ST-T segment also had high A(z) values. Using the four CDPs, the averaged total score for patients with triple-vessel disease was the highest ("2.67") compared to the other groups (normal controls: 0.53). Furthermore, based on the assumption that subjects with a total score over "1" were suspected of having CHD, sensitivity and specificity were 85.7% and 74.3%, respectively.

CONCLUSION

We concluded that the score and criteria using MCV and TCV during repolarization in CHD patients can reflect lesion areas and time changes of electrical activation dispersion due to ischemia.

Space-time database for standardization of adult magnetocardiogram-making standard MCG parameters

BACKGROUND

The magnetocardiogram (MCG) is a promising medical tool for detecting and visualizing abnormal cardiac electrical activation in heart-disease patients. However, there is no large-scale MCG database of healthy subjects, and there is little knowledge of gender- and age-related influences on MCG data.

METHODS AND RESULTS

We obtained MCG data from 869 subjects (554 men, 315 women) using a conventional 64-channel MCG system, which covers the whole heart. Electrocardiogram (ECG) data were also obtained; 464 people (268 men, 196 women) were identified as a normal group using ECG data. Time intervals (PQ, QRS, QT, and QTc), current distributions (maximum current vector (MCV), and the total current vector (TCV)) of MCG data of the 464 normal subjects were analyzed to obtain basic MCG parameters. Although mean values of PQ and QRS intervals of the male subjects were slightly longer than those of the female subjects, no intervals were correlated with gender or age. The correlation between PQ intervals of ECG and those of MCG was better than the correlation between QRS and QT intervals of ECG and those of MCG. Both MCV and TCV angles were much smaller than the electrical-axis angle in ECG. Although TCVs of the QRS and T waves were stable, the women's mean T-wave-TCV angles significantly increased with age. The maximum amplitude of the P wave was about 1.7 pT, and the value of the QRS complex was about 20-25 pT. Moreover, the T-wave amplitude decreases with age.

CONCLUSION

The MCG standard space-time parameters determined here provide a normal range for MCG parameters.

Detection of patients with coronary artery disease using cardiac magnetic field mapping at rest

We studied the use of cardiac magnetic field mapping to detect patients with CAD without subjecting them to stress. Fifty-nine healthy control subjects and 101 patients with CAD without previous MI were included. The optimal positions for detecting CAD were located in the left superior parasternal and in the inferior midsternal area. Values for ST slope, ST shift, T peak amplitude, ST-T integral, and magnetic field map orientation differed significantly between the 2 groups. Three parameters together in a multivariate analysis yielded a sensitivity of 84% and a specificity of 83% in distinguishing patients with CAD from control subjects. We suggest that cardiac magnetic field mapping is a promising technique to identify patients with CAD.

Visualization of cardiac dipole using a current density map: detection of cardiac current undetectable by electrocardiography using magnetocardiography

A close relationship exists between electric current and the magnetic field. However, electricity and magnetism have different physical characteristics, and magnetocardiography (MCG) may provide information on cardiac current that is difficult to obtain by electrocardiography (ECG). In the present study, we investigated the issue of whether the current density map method, in which cardiac current is estimated from the magnetic gradient, facilitates the visualization of cardiac current undetectable by ECG. The subjects were 50 healthy adults (N group), 40 patients with left ventricular overloading (LVO group), 15 patients with right ventricular overloading (RVO group), 10 patients with an old inferior myocardial infarction (OMI group), and 30 patients with diabetes mellitus (DM group). MCGs were recorded with a second derivative superconducting quantum interference device (SQUID) gradiometer using liquid helium. Isopotential maps and current density maps from unipolar precordial ECG leads and MCGs, respectively, were prepared, and the cardiac electric current was examined. The current density map at the ventricular depolarization phase showed one peak of current density in the N group. However, in the OMI group, the current density map showed multiple peaks of current density areas. In the RVO group, two peaks of current densities were detected at the right superior region and left thoracic region and these two diploles appeared to be from the right and left ventricular derived cardiac currents, respectively. Moreover, there was a significant correlation between the magnitude of the current density from the right ventricle and the systolic pulmonary arterial pressure. The current density map at the ventricular repolarization phase in the N group showed only a single current source. However, abnormal current sources in the current density maps were frequently detected even in patients showing no abnormalities on isopotential maps in the LVO, DM, and OMI groups. The findings herein suggest that opposing dipoles of the ventricular depolarization and repolarization vectors were summed and evaluated as a single dipole in the electrocardiogram. However, MCG facilitated the detection of multiple dipoles because of its superior spatial resolution as well as difference in physical properties between magnetic and electrical fields. Our results suggest that MCG with a current density map is useful for detecting cardiac current undetectable by ECG in an early stage.

Pseudo current density maps of electrophysiological heart, nerve or brain function and their physical basis

Background

In recent years the visualization of biomagnetic measurement data by so-called pseudo current density maps or Hosaka-Cohen (HC) transformations became popular.

Methods

The physical basis of these intuitive maps is clarified by means of analytically solvable problems.

Results

Examples in magnetocardiography, magnetoencephalography and magnetoneurography demonstrate the usefulness of this method.

Conclusion

Hardware realizations of the HC-transformation and some similar transformations are discussed which could advantageously support cross-platform comparability of biomagnetic measurements.

Reconstruction of action potential of repolarization in patients with congenital long-QT syndrome

A method for reconstructing an action potential during the repolarization period was developed. This method uses a current distribution-plotted as a current-arrow map (CAM)--calculated using magnetocardiogram (MCG) signals. The current arrows are summarized during the QRS complex period and subtracted during the ST-T wave period in order to reconstruct the action-potential waveform. To ensure the similarity between a real action potential and the reconstructed action potential using CAM, a monophasic action potential (MAP) and an MCG of the same patient with type-I long-QT syndrome were measured. Although the MAP had one notch that was associated with early afterdepolarization (EAD), the reconstructed action potential had two large and small notches. The small notch timing agreed with the occurrence of the EAD in the MAP. On the other hand, the initiation time of an abnormal current distribution coincides with the appearance timing of the first large notch, and its end time coincides with that of the second small notch. These results suggest that a simple reconstruction method using a CAM based on MCG data can provide a similar action-potential waveform to a MAP waveform without having to introduce a catheter.

Identifying patterns of spatial current dispersion that characterise and separate the Brugada syndrome and complete right-bundle branch block

The aim of the study was to detect patterns of spatial-current distribution in the late QRS and early ST-segments that distinguish Brugada-syndrome cases from complete right-bundle branch block (CRBBB). Magnetocardiograms (MCGs) were recorded from Brugada-syndrome patients (n = 6), CRBBB patients (n = 4) and the members of a control group (n = 33). The current distributions at six time points from Q-onset were estimated by producing current-arrow maps (CAMs). The angle of the current arrow of maximum amplitude at each time point was calculated. In the Brugada cases, the characteristic ST elevation was seen above the upper right chest, and abnormal currents appeared to be present in the right-ventricular outflow tract (RVOT). The angles of the abnormal arrows were -78 degrees +/- 51 degrees at 100 ms and -50 degrees +/- 61 degrees at 110 ms. In the cases of CRBBB, wide S- and R-waves were recorded above the upper right and lower right chest, respectively. The angles of the abnormal arrows for CRBBB were 152 degrees +/- 19 degrees at 100 ms, 159 degrees +/- 20 degrees at 110 ms, and 157 degrees +/- 19 degrees at 120 ms. The findings suggest that an abnormal current from the RVOT to the upper left chest may be a feature of the Brugada syndrome, and that the direction of this current is completely different from that seen in CRBBB.

Comparison of Automatic Repolarization Measurement Techniques in the Normal Magnetocardiogram

Multichannel MCG noninvasively measures cardiac magnetic field strength from many sites at the body surface, potentially providing useful regional information about ventricular repolarization. Previous work on ECGs has shown that automatic techniques for repolarization measurement are better than manual measurement at discriminating patients with cardiac conditions from normal subjects. Although automatic repolarization measurement techniques have been quantified for ECGs, no comparative data exists for the MCG. In this study four different automatic repolarization (QT) interval techniques for detecting T wave end in the MCG were compared. The influence of MCG filtering on the automatic algorithms was also quantified. MCGs were obtained at 49 sites over the heart from 23 normal subjects. Automatic measurements of the repolarization (QT) interval were made following the addition of different high pass (0.25, 0.5, 1 Hz) and low pass (100, 60, 40, 30 Hz) filters. There were consistent differences between automatic techniques in the unfiltered data amounting to greatest mean difference of 52.3 ms. Low pass filtering significantly increased the automatic repolarization (QT) interval relative to unfiltered measurement by 6.5 (3.2) ms (mean SD) for 100 Hz, 6.0 (3.0) ms for 60 Hz, 8.1 (3.2) ms for 40 Hz, and 8.8 (3.1) ms for 30 Hz across all techniques. High pass filtering significantly decreased the value by -2.6 (6.0) ms for 0.25 Hz, -5.5 (5.3) ms for 0.5 Hz, and -17.1 (7.8) ms for 1 Hz. Automatic measurements of repolarization (QT) in the MCG differ between techniques and are influenced by filtering. These effects should be considered when comparing results.

Classifying cases of fetal Wolff-Parkinson-White syndrome by estimating the accessory pathway from fetal magnetocardiograms

The paper presents an evaluation of the possibility of using fetal magnetocardiogram (FMCG) signals to estimate and classify the accessory pathway in fetal Wolff-Parkinson-White (WPW) syndrome. The FMCG signals of two fetuses with WPW syndrome (type A) were detected using a 64-channel superconducting quantum-interference device system. An average across the cycles of these signals was taken to obtain clear WPW signals. To determine the direction and position of the accessory pathway in a fetal heart accurately, the accessory pathway and activated pathway at the peak of the QRS complex thus obtained were estimated for each fetus, using a single-dipole model. The phase angle (about 90 degrees) between the equivalent current dipoles (ECDs) was the same for both fetuses. This angle suggested that the accessory pathway is in the left side of the heart, i.e. that the pathway exists in the left ventricle, which indicates type A WPW syndrome. Identification of the position of the accessory pathway in a fetus with WPW syndrome from the angle between the ECD of the accessory pathway and the ECD of the peak in the QRS complex was thus demonstrated.

Magnetocardiographic indices of left ventricular hypertrophy

OBJECTIVE

We tested the hypothesis that multichannel magnetocardiographic (MCG) mapping can detect and quantify the degree of left ventricular hypertrophy (LVH).

DESIGN

A cross-sectional study.

SETTING

Helsinki University Central Hospital, a tertiary referral center.

PARTICIPANTS

Forty-two patients with pressure overload induced LVH by gender-specific echocardiographic criteria (LVH group), and 12 healthy middle-aged controls.

MAIN OUTCOME MEASURES

MCG QRS-T area integrals and QRS-T angle in magnetic field maps in relation to echocardiographic LVH as well as left ventricular (LV) mass and structure. Conventional 12-lead electrocardiographic (ECG) LVH indices (Sokolow-Lyon voltage, Cornell voltage, Cornell voltage duration product) were assessed for comparison.

RESULTS

MCG QRS- and T-wave integrals provided complementary information of echocardiographic LV mass. Their combination, the QRS-T integral, and the QRS-T angle were increased in patients with LVH and, in those patients, correlated significantly with LV mass indexed to body surface area (r = 0.455;P = 0.002 and r= 0.379; P= 0.013, respectively). A QRS-T integral 16000 fT.s had identical sensitivity of 62% at 92% specificity as the gender-adjusted Cornell voltage duration product of 240 micro V.s for the detection of LVH.

CONCLUSIONS

The MCG method can detect patients with LVH and also quantify the degree of LVH in patients with increased LV mass.

Detection of spatial repolarization abnormalities in patients with LQT1 and LQT2 forms of congenital long-QT syndrome

The aim of this study is to detect the spatial current dispersion that appears in the T-wave of patients with congenital long-QT syndrome (LQTS). To observe this dispersion, magnetocardiograms (MCGs)--which have a high spatial resolution--of LQT1 patients (n = 7), LQT2 patients (n = 9) and a control group (n = 33) were recorded. The dispersion was evaluated by plotting current-arrow maps (CAMs) calculated from the MCG signals. In the case of LQT1, abnormal current arrows in the CAMs appeared above the inferior part of the heart in two LQT1 patients with a long corrected QT interval (QTc) (>0.6), and the current direction was from the left (origin side) to the right ventricular muscle (110 degrees). In six out of nine LQT2 patients, abnormal current arrows with angles below 20 degrees were observed above the right inferior part or lower septum; the current direction was from the right (origin side) to the left ventricular muscle. However, in the case of the LQT2 patients, the QTc values did not correlate with the abnormal current. These findings suggest that the origin of abnormal repolarization in LQT1 is the left ventricular muscle and the origin of that in LQT2 is the right ventricular muscle or lower septum. The estimation of the origin in LQTS patients can provide important information such as the risk factor of sudden death.

Abnormal auditory neural networks in patients with right hemispheric infarction, chronic dizziness, and moyamoya disease: a magnetoencephalogram study

The purpose of this study was to determine whether the auditory cortex is sensitive to cortical insults and to determine the specificity of the insults in three clinical situations with different cortical involvement. Auditory-evoked magnetic fields of ten normal subjects, 8 patients with right hemispheric infarction, 11 with chronic dizziness, and 2 with moyamoya disease were measured. To analyze the abnormality of auditory neural networks, the magnitude ratio and the angle difference (Deltatheta) between response vectors, which were determined from maximum current arrows corresponding to the N100m peak for contralateral and ipsilateral stimuli were used. A normal range of the parameters was defined so that abnormal values could be determined. Of the three parameters, Deltatheta was the most sensitive: 4 patients with right hemispheric infarction, 4 with chronic dizziness, and 1 with moyamoya disease had abnormal Deltatheta. The electrical activity in the patients with such abnormal Deltathetas had a circular current pattern. These findings suggest that right infarction lesions sometime affect the left auditory neural network, dizziness is caused by abnormal neural networks between the vestibular cortical area and the auditory cortex or by an imbalance between left and right auditory-cortex activities, and moyamoya-disease patients have almost normal auditory-electrical activity.

Detection of atrial-flutter and atrial-fibrillation waveforms by fetal magnetocardiogram

Two cases of fetal tachycardia are reported: atrial flutter and fibrillation. The waveforms from each case were detected by fetal magnetocardiograms (FMCGs) using a 64-channel superconducting quantum interference device (SQUID) system. Because the magnitude of supraventricular arrhythmia signals is very weak, two subtraction methods were used to detect the fetal MCG waveforms: subtraction of the maternal MCG signal, and subtraction of the fetal ORS complex signal. It was found that atrial-flutter waveforms showed a cyclic pattern and that atrial-fibrillation waveforms showed f-waves with a random atrial rhythm. Fast Fourier transform analysis determined the main frequency of the atrial flutter to be about 7Hz, and the frequency distribution of atrial fibrillation consisted of small, broad peaks. To visualise the current pattern, current-arrow maps, which simplify the observation of pseudo-current patterns in fetal hearts, of the averaged atrial flutter and fibrillation waveforms were produced. The map of the atrial flutter had a circular pattern, indicating a re-entry circuit, and the map of the atrial fibrillation indicated one wavelet, which was produced by a micro-re-entry circuit. It is thus concluded that an FMCG can detect supraventricular arrhythmia, which can be characterised by re-entry circuits, in fetuses.

Prenatal diagnosis of QT prolongation by fetal magnetocardiogram--use of QRS and T-wave current-arrow maps

To determine the T wave of a fetal magnetocardiogram (FMCG), we have evaluated the T/QRS ratio and obtained current-arrow maps that indicate weak currents. We measured FMCG signals for 52 normal fetuses and two abnormal fetuses with prolonged QT waves by using three superconducting quantum interference device (SQUID) systems: a nine-channel system, a 12-channel vector system and a 64-channel system. The T/QRS ratio was calculated for all the normal fetuses from the maximum magnitudes of the QRS complex and the T wave. Current-arrow maps of the QRS complex (R wave) and T wave were obtained by using the 64-channel system, and the phase differences of the total-current vectors were calculated by using the current-arrow maps. The results showed that the T/QRS ratio had a wide variability of 0.35 for the normal fetuses. However, the magnitude of the prolonged T wave was as weak as the detection limit of the SQUID magnetometer. Although the T/QRS ratios for the fetuses with QT prolongation were within the normal range (< 0.35), the weak magnitude of the prolonged T wave could be evaluated. On the other hand, by comparing the current-arrow maps of the R and T waves for the normal fetuses, we found that the maximum-current arrows were indicated as either in the same direction or in opposite directions. These patterns could be identified clearly by the phase differences. Very weak prolonged T waves for the two abnormal fetuses could be determined by using these current-arrow maps and phase differences. Consequently, although the T/QRS ratios of FMCG signals have a wide distribution, we have concluded that the current-arrow map and phase difference can be used to determine the T wave of an FMCG signal.

A method for detecting myocardial abnormality by using a total current-vector calculated from ST-segment deviation of a magnetocardiogram signal

A simple method to determine the state of ischaemia or fibrosis of myocardial cells has been developed. This method uses the ST wave of 64-channel magnetocardiogram (MCG) signals to calculate three parameters from the current-arrow map of the normal component signal of the MCG. One parameter is a total current vector that is obtained through summation of all current arrows. Another is a variance current vector calculated from the differential vector of two total current vectors at different times. The third is a flatness factor between the magnitude of the total current vector and the variance current vector. The three parameters are independent of the distance between the heart and the gradiometers. We measured the MCG signals of 29 healthy subjects, twenty patients with coronary artery disease (ten with previous myocardial infarction (MI) and ten with angina pectoris (AP)), and eight patients with cardiomyopathy (four with hypertrophic cardiomyopathy (HCM), three with dilated cardiomyopathy (DCM), and one with restrictive cardiomyopathy (RCM)). With our method, none of the healthy subjects tested positive for myocardial abnormalities, while 80% of the MI patients, 50% of the AP patients, and 100% of the cardiomyopathy patients tested positive. Although further testing is needed, we feel this simple technique enables easy diagnosis of myocardial damage.