Long-term risk prediction is a priority for the prevention of atrial fibrillation (AF). P wave indices are electrocardiographic measurements describing atrial conduction. The role of P wave indices in the prospective determination of AF and mortality risk has had limited assessment. We quantified by digital caliper the P wave indices of maximum duration and dispersion in 1,550 Framingham Heart Study participants ≥60 years old (58% women) from single-channel electrocardiograms recorded from 1968 through 1971. We examined the association of selected P wave indices and long-term outcomes using Cox proportional hazards regression incorporating age, gender, body mass index, systolic blood pressure, treatment for hypertension, significant murmur, heart failure, and PR interval. Over a median follow-up of 15.8 years (range 0 to 38.7), 359 participants developed AF and 1,525 died. Multivariable-adjusted hazard ratios (HRs) per SD increase in maximum P wave duration were 1.15 (95% confidence interval [CI] 0.90 to 1.47, p = 0.27) for AF and 1.02 (95% CI 0.96 to 1.08, p = 0.18) for mortality. The upper 5% of P wave maximum duration had a multivariable-adjusted HR of 2.51 (95% CI 1.13 to 5.57, p = 0.024) for AF and an HR of 1.11 (95% CI 0.87 to 1.40, p = 0.20) for mortality. We found no significant associations between P wave dispersion with incidence of AF or mortality. In conclusion, maximum P wave duration at the upper fifth percentile was associated with long-term AF risk in an elderly community-based cohort. P wave duration is an electrocardiographic endophenotype for AF.
We sought to assess prospectively the association between P wave indices and clinical outcomes in the Framingham Heart Study (FHS). Specifically, we hypothesized that prolonged P wave duration and dispersion were associated with atrial fibrillation (AF) risk. We used a digital caliper to determine P wave duration and dispersion from single-channel electrocardiograms (ECGs) and then examined their relation to the outcomes of AF, heart failure (HF), and mortality in a prospective analysis.
Methods
The FHS is a longitudinal community-based study designed to investigate risk factors for cardiovascular disease. The study was initiated in 1948 with the enrollment of an original cohort. The present study used baseline data from the original cohort’s 11th biennial examination, which occurred from 1968 through 1971 and included 2,955 participants. Participants at this examination underwent comprehensive phenotypic assessment that included physical examination, clinical history, laboratory assessment including cholesterol levels, and single-channel ECG. Subjects <60 years of age (n = 1,323), with a history of AF (n = 29), or with illegible ECGs for measurement acquisition as described later (n = 53) were excluded from the present analysis. No subjects had pacemaker implantation or missing covariate data. After exclusions, 1,550 participants (905 women) were eligible for inclusion. All participants provided written informed consent and the study was approved by Boston University Medical Center’s institutional review board.
Participants at index examinations underwent comprehensive clinical assessments including physical examination, laboratory assessment, and review of medications, clinical history, and interim clinical events. Participants had a single-channel ECG recorded sequentially on lined paper at 25 mm/s and 0.1 mV/mm (Hewlett-Packard, Palo Alto, California) sequentially for each of the 12 electrocardiographic leads as described previously. ECGs were converted to tagged image format using commercially available software for electrocardiographic measurement modified for the present study (Rigel 1.7.4, AMPS, New York, New York; available at: http://www.amps-llc.com ). The software uses a computer mouse for digital caliper manipulation, fourfold magnification, and automated measurement recording to avoid transcription error. P wave indices such as area and terminal force were not determined in the present study owing to the design of the software.
Measurements consisted of identification of the initiation of P wave onset and its termination at P wave offset. P wave onset was determined as a positive deflection after the T wave preceding the QRS waveform and deviating from the isoelectric line comprised of the TP interval. P wave offset was determined as the return to isoelectric baseline immediately before the QRS waveform. An analogous protocol was employed for negative deflections consistent with P waves in lead V 1 . Up to 3 consecutive P waves were measured in each of the 12 leads to account for physiologic variability. Leads with ambiguous P wave onset or offset were not included for measurement. All measurements were determined by a single investigator (J.W.M.). Previous assessment of quality control had demonstrated an intrarater intraclass correlation coefficients of r = 0.80 for maximum P wave duration and r = 0.82 for P wave dispersion (defined later). Of note, intraclass correlation coefficients of 0.50, 0.75, and 1.00 suggest moderate, excellent, and perfect correlations, respectively. The PR interval was quantified independent of the present investigation as described previously.
P wave indices employed included maximum P wave duration, minimum P wave duration, and P wave dispersion. Maximum P wave duration was defined as the longest and minimum as the shortest duration P wave from all leads. P wave dispersion was defined as the difference between the maximum and minimum P wave durations. P wave indices were derived from mean measurements and measurements were recorded in milliseconds.
FHS participants receive ongoing surveillance for cardiovascular events and mortality; records are retrieved from outside medical facilities and hospitalizations. Adjudication for cardiovascular events, including HF, is performed by 3 FHS investigators who review pertinent data from all available hospitalizations, ambulatory medical records, and FHS examinations. ECGs are reviewed by a study cardiologist and classified as AF or atrial flutter, as previously documented. As in previous analyses from the FHS, participants with tracings showing atrial flutter were combined with those showing AF.
Descriptive statistics were generated for the 1,550 subjects included in the analysis. Continuous variables were described using means ± SDs and categorical variables using percentages. Mean P wave duration was computed for each lead by summing the P wave durations for each lead and dividing by the number of leads measured, followed by determination of maximum and minimum P waves defined by magnitude of P wave duration across the 12 leads. Correlations between the P wave indices were calculated.
We constructed multivariable Cox proportional hazards regression models to examine the association between P wave indices (PR interval, maximum P wave duration, and P wave dispersion) and the separate outcomes of AF, HF, and all-cause mortality. We determined the hazard ratio (HR) and 95% confidence interval (CI) relative to a 1-SD increase in each of the P wave indices. In secondary analyses we stratified P wave indices by >95th percentile (i.e., 5% longest measurements) compared to those in ≤95th percentile. All analyses were initially adjusted for age and gender. Multivariable adjustment was performed using age, gender, body mass index, systolic blood pressure, treatment for hypertension, significant murmur (≥3/6 systolic murmur or presence of any diastolic murmur), and HF. These variables were derived from a previous Framingham AF risk score analysis. HF was diagnosed by established validated criteria. Cases of prevalent HF were excluded from analyses examining the outcome of HF. For maximum P wave duration, multivariable adjustment also included PR interval. We confirmed the proportionality of hazards assumption separately for each of the 3 outcomes.
All analyses were performed using SAS 9.1.3 (SAS Institute, Cary, North Carolina). Statistical significance was determined by a 2-sided p value <0.05.
Results
Of 1,550 subjects included in the present investigation, 58% were women and mean age was 68.2 ± 5.5 years. Median follow-up was 15.8 years (range 0.0 to 38.1). Six percent had a previous myocardial infarction and 2% prevalent HF. Further baseline characteristics of participants and distribution of the P wave indices are presented in Table 1 .
Entire Cohort | Men | Women | |
---|---|---|---|
(n = 1,550) | (n = 645) | (n = 905) | |
Clinical characteristics | |||
Age (years) | 68.2 ± 5.5 | 68.0 ± 5.4 | 68.3 ± 5.6 |
Women | 58% | — | — |
Body mass index (kg/m 2 ) | 25.6 ± 3.9 | 25.7 ± 3.2 | 25.5 ± 4.3 |
Systolic blood pressure (mm Hg) | 146 ± 22 | 143 ± 22 | 148 ± 23 |
Treatment for hypertension | 20% | 14% | 24% |
Significant precordial murmur ⁎ | 5% | 4% | 5% |
Prevalent heart failure | 2% | 2% | 2% |
Myocardial infarction | 6% | 10% | 3% |
Heart rate (beats/min) | 77 ± 13 | 74 ± 13 | 79 ± 13 |
Electrocardiographic left ventricular hypertrophy | 4% | 4% | 3% |
QRS interval (ms) | 77 ± 15 | 80 ± 16 | 75 ± 14 |
P wave indices | |||
P wave duration (ms) | 86 ± 9 | 87 ± 9 | 86 ± 9 |
Maximum P wave duration (ms) | 105 ± 12 | 108 ± 12 | 104 ± 12 |
P wave dispersion (ms) | 41 ± 11 | 43 ± 12 | 39 ± 11 |
PR interval (ms) | 163 ± 23 | 166 ± 25 | 160 ± 22 |
⁎ Defined as systolic murmur grade >3 of 6 or any diastolic murmur auscultated by a physician.
Correlations between P wave indices assessed are presented in Table 2 . Maximum P wave duration and P wave dispersion were strongly correlated (r = 0.74). PR interval was moderately correlated with maximum P wave duration (r = 0.34) and less so with P wave dispersion (r = 0.17).
Maximum P Wave Duration | P Wave Dispersion | PR Interval | |
---|---|---|---|
Maximum P wave duration | 1.00 | ||
P wave dispersion | 0.74 | 1.00 | |
PR interval | 0.34 | 0.17 | 1.00 |
By the conclusion of follow-up, 359 participants developed AF, 337 developed HF, and 1,525 died; these categories were nonexclusive. Age- and gender-adjusted incidence rates per 10,000 person-years were stratified at the 95th percentile of measurements and are presented in Table 3 . The AF incidences per 10,000 person-years were 143 (95% CI 117 to 166) at ≤95th percentile and 228 (95% CI 126 to 297) at >95th percentile for maximum P wave duration. Incidence rates per 10,000 years for mortality were similar across the strata of maximum P wave duration (324 for ≤95th percentile, 95% CI 307 to 338, and 351 for >95th percentile, 95% CI 300 to 387). P wave dispersion for ≤95th percentile had a higher incidence of AF, HF, and mortality compared to those values in >95th percentile.
P Wave Indices | Atrial Fibrillation | Heart Failure | All-Cause Mortality |
---|---|---|---|
Maximum P wave duration | |||
≤95th percentile | 143 (117–166) | 119 (95–141) | 324 (307–338) |
>95th percentile | 228 (126–297) | 108 (41–163) | 351 (300–387) |
P wave dispersion | |||
≤95th percentile | 149 (122–173) | 122 (97–144) | 325 (309–340) |
>95th percentile | 105 (37–161) | 69 (19–114) | 319 (265–358) |
PR interval | |||
≤95th percentile | 139 (113–162) | 118 (94–140) | 323 (306–337) |
>95th percentile | 224 (147–281) | 127 (64–178) | 346 (308–375) |
Hazards for associations between P wave indices and outcomes are presented in Table 4 . Gender- and age-adjusted HR per SD increase in maximum P wave duration was 1.12 (95% CI 1.00 to 1.25, p = 0.047) for AF and multivariable-adjusted HR was 1.15 (95% CI 0.90 to 1.47, p = 0.08). When stratified at the 95th percentile in a secondary analysis, the upper 5% of maximum P wave duration had an HR of 2.19 (95% CI 1.46 to 3.30, p <0.001) for AF in gender- and age-adjusted analysis and an HR of 2.51 (95% CI 1.13 to 5.57, p = 0.024) in multivariable-adjusted analysis. In comparison, PR interval stratified at the 95th percentile conveyed an HR of 1.65 (95% CI 1.20 to 2.29, p = 0.002) for AF in age- and gender-adjusted analysis and an HR of 1.59 (95% CI 0.83 to 3.03, p = 0.16) in multivariable-adjusted analysis. P wave maximum duration was not significantly associated with mortality when examined by SD increase. After stratification at the 95th percentile, the upper 5% of maximum P wave duration had an HR of 1.28 (95% CI 1.01 to 1.61, p = 0.04) for all-cause mortality in age- and gender-adjusted analysis and an HR of 1.11 (95% CI 0.87 to 1.40, p = 0.20) in multivariable-adjusted analysis. In contrast, PR interval remained significantly associated with mortality in all analyses.
P wave Indices | Atrial Fibrillation | All-Cause Mortality | ||
---|---|---|---|---|
HR (95% CI) | p Value | HR (95% CI) | p Value | |
Maximum P wave duration | ||||
Age and gender adjusted, 1-SD increase ⁎ | 1.12 (1.00–1.25) | 0.047 | 1.04 (0.99–1.10) | 0.13 |
Multivariable adjusted, 1-SD increase † | 1.15 (0.90–1.47) | 0.27 | 1.02 (0.96–1.08) | 0.18 |
Age and gender adjusted, >95th percentile ⁎ | 2.19 (1.46–3.30) | <0.001 | 1.28 (1.01–1.61) | 0.04 |
Multivariable-adjusted, >95th percentile † | 2.51 (1.13–5.57) | 0.024 | 1.11 (0.87–1.40) | 0.20 |
P wave dispersion | ||||
Age and gender adjusted, 1-SD increase ⁎ | 1.00 (0.90–1.11) | 1.00 | 1.03 (0.98–1.08) | 0.32 |
Multivariable adjusted, >95th percentile † | 0.84 (0.20–3.54) | 0.81 | 0.93 (0.74–1.17) | 0.55 |
PR interval | ||||
Age and gender adjusted ⁎ | 1.21 (1.08–1.34) | <0.001 | 1.07 (1.02–1.13) | 0.001 |
Multivariable adjusted † | 1.15 (0.90–1.47) | 0.27 | 1.07 (1.01–1.12) | 0.012 |
Age and gender adjusted, >95th percentile ⁎ | 1.65 (1.20–2.29) | 0.002 | 1.27 (1.07–1.50) | 0.006 |
Multivariable adjusted, >95th percentile † | 1.59 (0.83–3.03) | 0.16 | 1.24 (1.05–1.47) | 0.011 |