Abstract

BackgroundThe first pandemic season of A/H3N2 influenza virus (1968/1969) resulted in significant mortality in the United States, but it was the second pandemic season of A/H3N2 influenza virus (1969/1970) that caused the majority of deaths in England. We further explored the global pattern of mortality caused by the pandemic during this period

MethodsWe estimated the influenza-related excess mortality in 6 countries (United States, Canada, England and Wales, France, Japan, and Australia) using national vital statistics by age for 1967–1978. Geographical and temporal pandemic patterns in mortality were compared with the genetic drift of the influenza viruses by analyzing hemagglutinin and neuraminidase sequences from GenBank

ResultsIn North America, the majority of influenza-related deaths in 1968/1969 and 1969/1970 occurred during the first pandemic season (United States, 70%; Canada, 54%). Conversely, in Europe and Asia, the pattern was reversed: 70% of deaths occurred during the second pandemic season. The second pandemic season coincided with a drift in the neuraminidase antigen

ConclusionWe found a consistent pattern of mortality being delayed until the second pandemic season of A/H3N2 circulation in Europe and Asia. We hypothesize that this phenomenon may be explained by higher preexisting neuraminidase immunity (from the A/H2N2 era) in Europe and Asia than in North America, combined with a subsequent drift in the neuraminidase antigen during 1969/1970

Annual influenza epidemics are sustained in the human population through gradual mutations in hemagglutinin and neuraminidase, the surface antigens of the virus. The genetic makeup of the influenza virus allows frequent minor drifts every 2–5 years in response to selection pressure to evade human immunity [1]. Rarely, reassortment between human and nonhuman viruses results in larger shifts, in which a new virus subtype emerges and replaces the previously circulating one [2]. A new pandemic virus rapidly invades the human population with partial or no immunity and may cause severe illness worldwide [3, 4]. Although the impact of influenza is not always higher during pandemics than during interpandemic periods, a shift in the age distribution of mortality toward younger age groups distinguishes pandemic from epidemic impact [5, 6]

The influenza virus responsible for the last pandemic, A/Hong Kong/68 (A/H3N2), was first isolated in Hong Kong in July 1968 [7]. The new A/H3N2 virus exhibited a shift in hemagglutinin but not in neuraminidase; it replaced A/H2N2 viruses that had been circulating in all countries since 1957. Despite rapid and extensive spread by international air travel [7, 8], the impact of the new virus was not the same in all geographical regions. A marked increase in mortality occurred in the United States during the first pandemic season (1968/1969), especially in persons <65 years old, but was not seen elsewhere [7]. Conversely, in England, the second pandemic season (1969/1970) of A/H3N2 virus proved to be more severe than the first [9, 10]

The reasons for the delayed severe impact in England are still not understood [4, 11]. Such a delay is counterintuitive, since a novel virus introduced in a susceptible population should demonstrate decreasing impact over time as immunity increases [9, 12]. Here, we analyze monthly mortality data on 6 countries (on 4 continents) and review published morbidity and virological studies to extend the current understanding of the Hong Kong A/H3N2 pandemic. To explain the epidemiological patterns, we investigate the genetic sequences of influenza surface antigens. Finally, we discuss the implications of these patterns for pandemic preparedness

Methods

A detailed description of the data sources and analytic approach is given in the Appendix

Data Sources

Mortality and population dataMonthly age-specific data on pneumonia and influenza (P&I) and all-cause mortality were compiled from vital statistics on the United States, Canada, England and Wales (referred to as “England” for simplicity), France, Japan, and Australia for 1967–1979, which includes the 2 pandemic seasons, 1968/1969 and 1969/1970. No details on age were available for Japan. Before analysis, we calculated the monthly incidence per 100,000 persons for both mortality outcomes

Review of the literature for virological surveillance and morbidity dataWe compiled literature reports, mostly from World Health Organization sources, to compare the date of the first isolation of the pandemic strain and the timing of the 2 pandemic seasons in each country [13–19 ] (figure 1). To compare influenza morbidity patterns, we analyzed published weekly time series of influenza virological isolates and summary estimates of morbidity. We restricted the analysis to follow-up studies that covered both pandemic seasons [9, 17, 20–23 ]. We also searched the available published literature for general-population influenza seroepidemiological surveys that reported the prevalence of new influenza infections [9, 20–24 ]

Figure 1

Timing of the first 2 pandemic seasons of A/H3N2 influenza virus in 6 countries. Shown are isolation, sporadic, and epidemic activities, as defined by virus surveillance reports [13–19 ]

Figure 1

Timing of the first 2 pandemic seasons of A/H3N2 influenza virus in 6 countries. Shown are isolation, sporadic, and epidemic activities, as defined by virus surveillance reports [13–19 ]

Influenza genetic sequencesTemporal and/or geographical differences in the antigens of circulating viruses may explain intercountry differences in patterns of mortality. For viruses collected during 1966–1975, we analyzed all hemagglutinin and neuraminidase genetic sequences published in GenBank (http://www.ncbi.nlm.nih.gov/GenBank/) and in the Influenza Sequence Database (http://www.flu.lanl.gov) [25]. We chose this period in order to include strains that circulated just before the pandemic until the first major postpandemic antigenic change (A/Victoria/75) [11]. We analyzed 38 sequences for the hemagglutinin H3 gene and 53 sequences for the neuraminidase N2 gene

Analytic Approach

To evaluate the impact of the 2 pandemic seasons in the 6 countries, we calculated P&I and all-cause seasonal excess mortality and adjusted for age and nondemographic factors

Estimating P&I and all-cause seasonal excess mortality P&I and all-cause seasonal excess mortality were computed for 2 age categories (all ages and <65 years old) as the increase in mortality above a seasonal baseline, during epidemic months, for both the 1968/1969 and the 1969/1970 pandemic seasons. We used a Serfling-type regression model [26] and identified influenza epidemic months from the time series of deaths specifically attributed to influenza

We adjusted excess mortality using the mid-1969 population of England as the reference population, to account for differences in age structure between countries. To compensate for nondemographic differences (e.g., access to health care and coding for underlying cause of death [27, 28]), we calculated the percentage increase in mortality as the excess mortality divided by the baseline mortality during winter (expected mortality) for both P&I and all-cause seasonal excess mortality. This measure was introduced by Serfling, for comparison of influenza mortality in different age groups [29]

Exposure to influenza during the last A/H2N2 epidemic season (1967/1968) may have reduced the impact of the A/H3N2 pandemic by granting protection against neuraminidase [3, 11]. We also analyzed the age-standardized excess mortality and corresponding percentage increase for this epidemic season

Proportion of deaths, illnesses, and infections during each pandemic season (1968/1969 and 1969/1970)To estimate the relative mortality during each pandemic season, we summed the excess mortality for 1968/1969 and 1969/1970 and computed the proportion occurring during each pandemic season separately for each country. We compared the resulting estimates with those of previous studies of mortality [5, 9, 24, 30–35 ]. Similarly, we computed the proportion of clinical illnesses and infections identified by serologic tests during each pandemic season from the available published literature [9, 17, 20–24 ]

Results

Timing of the First 2 A/H3N2 Pandemic Seasons

In most countries, the pandemic strain was isolated shortly after its appearance in Hong Kong in July 1968 (figure 1), with several months of sporadic activity before the epidemic took off. There was no documented cocirculation of A/H3N2 and A/H2N2 viruses in the same geographical area; the last reported isolation of A/H2N2 was in August 1968 in Australia [14, 36]

Mortality

Epidemic pattern and impact of the first 2 A/H3N2 pandemic seasons (all ages)The epidemic patterns of the 2 pandemic seasons were different in the 6 countries studied (figure 2). The P&I mortality time series reveals a large epidemic in the United States in 1968/1969, followed by a milder one in 1969/1970, late in the winter season. In Canada, the 2 epidemic patterns were similar in amplitude and timing. In the other 4 countries, entirely different patterns emerged. The first epidemic was mild, followed by a much more intense epidemic the next season. Similar patterns were observed from all-cause mortality time series (data not shown)

Figure 2

Influenza epidemic mortality curves for the first 2 pandemic seasons of A/H3N2 virus in 6 countries. Pneumonia and influenza (P&I) monthly mortality time series/100,000 persons are represented for the United States, England and Wales, Canada, France, Japan, and Australia. Note the geographical differences in the relative mortality of each pandemic season. First pandemic season, 1968/1969; second pandemic season, 1969/1970. In Australia, the first pandemic season was in 1969, and the second pandemic season was in 1970 (figure 1)

Figure 2

Influenza epidemic mortality curves for the first 2 pandemic seasons of A/H3N2 virus in 6 countries. Pneumonia and influenza (P&I) monthly mortality time series/100,000 persons are represented for the United States, England and Wales, Canada, France, Japan, and Australia. Note the geographical differences in the relative mortality of each pandemic season. First pandemic season, 1968/1969; second pandemic season, 1969/1970. In Australia, the first pandemic season was in 1969, and the second pandemic season was in 1970 (figure 1)

Of the sum of P&I excess deaths that occurred during 1968/1969 and 1969/1970, in the United States, 70% occurred during the first pandemic season; in Canada, 54%; in Japan, 32%; in England, 23%; in Australia, 22%; and, in France, 15% (table 1). The same trend was obtained for all-cause excess mortality, which supports a 1-year delay in the major impact in England, France, Japan, and Australia

Table 1

Influenza mortality impact of the last A/H2N2 epidemic season and the first 2 A/H3N2 pandemic seasons

Table 1

Influenza mortality impact of the last A/H2N2 epidemic season and the first 2 A/H3N2 pandemic seasons

Age standardization and computation of the percentage increase reduced but did not eliminate the intercountry differences in the sum of excess deaths that occurred during 1968/1969 and 1969/1970, for both P&I and all causes (table 1). For instance, P&I excess mortality was 2–3-fold higher in England and France than elsewhere. All-cause excess mortality was 2–5-fold lower in North America than elsewhere, and the corresponding percentage increase remained lower, suggesting that the total pandemic mortality was substantially lower in North America than in other countries

Last epidemic season of the A/H2N2 era (1967/1968) Similar intercountry differences were found for the last A/H2N2 epidemic season (table 1). P&I excess mortality was 2–3-fold higher in England and France than elsewhere. All-cause excess mortality was lower in North America than in other countries. The percentage increase in P&I and all-cause mortality was consistently lower in North America than elsewhere, supporting a lower impact for the last A/H2N2 epidemic season in this region. Interestingly, in all countries, all-cause excess mortality during the last A/H2N2 epidemic season was greater than or equal to that during the major A/H3N2 pandemic season

Pandemic “age shift” and impact in persons <65 years old (1968/1969 and 1969/1970)For all countries studied, the proportion of excess P&I mortality in persons <65 years old increased substantially (by 2.2–4.6-fold) during the first A/H3N2 pandemic season, compared with that during the last A/H2N2 epidemic season (table 1). This signature age shift of a novel virus subtype, taken together with virus surveillance data that indicated exclusive A/H3N2 circulation, clearly supports pandemic virus activity in all countries throughout the first A/H3N2 season

Intercountry differences in the relative impact of each pandemic season reported for all ages also existed for persons <65 years old (table 2). In this age group, the majority of P&I and all-cause mortality occurred during the first pandemic season in North America but during the second pandemic season in England, France, and Australia. Regarding the sum of deaths during the 2 pandemic seasons, all of our mortality measures showed a lower impact in North America in persons <65 years old than in Europe and Asia

Table 2

Influenza mortality of the first 2 A/H3N2 pandemic seasons in persons <65 years old

Table 2

Influenza mortality of the first 2 A/H3N2 pandemic seasons in persons <65 years old

Geographically distinct pandemic patternsFor the Hong Kong A/H3N2 pandemic, 2 geographically distinct mortality patterns emerged. The North American pattern (United States and Canada) was characterized by a first pandemic season that was more severe than the second, although, the impact of each pandemic season was more balanced in Canada than in the United States. There were minimal geographical variations in local pandemic patterns across the United States (Appendix). By contrast, local pandemic patterns were heterogeneous in Canada, which explains the more-balanced impact in the national analysis

The “smoldering” pattern in Europe and Asia (England, France, Japan, and Australia) was characterized by a second pandemic season 2–5 times more severe than the first. The last A/H2N2 epidemic season was less severe in North America than elsewhere

Review of the Literature for Proportion of Deaths, Illnesses, and Infections during Each Pandemic Season (1968/1969 and 1969/1970)

On average, across all mortality studies, including the present one, in the United States, an estimated 70% of influenza-related deaths occurred during the first pandemic season; in Canada, 57%; in England, 30%; and, in Australia, 31% (table 3). These studies gave fairly consistent estimates, although they did not use the same approach for estimating excess mortality. No study other than ours included France and Japan

Table 3

Relative impact of the first 2 influenza A/H3N2 pandemic seasons on mortality, morbidity, and serological infections in the United States, Canada, England, and Australia

Table 3

Relative impact of the first 2 influenza A/H3N2 pandemic seasons on mortality, morbidity, and serological infections in the United States, Canada, England, and Australia

Morbidity time-series analysis revealed a pattern similar to the mortality time series in the United States, England, and Australia, in terms of amplitude and timing (figure 3A–3C). In these 3 countries, the relative morbidity and mortality of each pandemic season were very close (table 3), ruling out potential differences in case fatalities. In all countries, high morbidity was associated with high mortality

Figure 3

Influenza epidemic morbidity curves for the first 2 pandemic seasons of A/H3N2 virus. A Weekly no. of influenza A isolates in England [9]. B Weekly no. of influenza-like illnesses (ILIs) in Michigan [20–23 ]. C Weekly no. of influenza A isolates in New South Wales, Australia [17]

Figure 3

Influenza epidemic morbidity curves for the first 2 pandemic seasons of A/H3N2 virus. A Weekly no. of influenza A isolates in England [9]. B Weekly no. of influenza-like illnesses (ILIs) in Michigan [20–23 ]. C Weekly no. of influenza A isolates in New South Wales, Australia [17]

General-population seroepidemiological surveys reported that the majority of influenza infections occurred during the first pandemic season in the United States and Canada, closely matching the pattern of deaths and clinical illnesses (table 3). Conversely, in England, patterns of serological infections were discordant with mortality and morbidity patterns. Infections during the first pandemic season were nearly as frequent as were those during the second pandemic season, whereas the large majority of deaths and clinical illnesses occurred during the second pandemic season. These discordances imply that asymptomatic (subclinical) infections were frequent during the first pandemic season in England. Similarly, in the other countries in Europe and Asia that were studied, high rates of asymptomatic infections during the first pandemic season were reported [15, 17, 37]. From the combined analysis of mortality, morbidity, and serological studies, we concluded that asymptomatic infections were frequent during the first pandemic season in Europe and Asia and were less frequent in North America

Phylogenetic Studies

We searched for temporal and/or geographical differences in the genes of influenza surface antigens. The hemagglutinin H3 gene sequences available from the 2 pandemic seasons clustered in a single group in the hemagglutinin tree (n=13 for 1968/1969–1969/1970) (figure 4A). A drift in this gene appeared in 1971/1972, which was later than our period of interest. Taken together with the single antigenic cluster defined by Smith et al. for the hemagglutinin between July 1968 and late 1970 [1], the single genetic cluster found here suggests that there was no significant evolution in this antigen during the 2 pandemic seasons

Figure 4

Phylogenetic relationships of the hemagglutinin H3 gene (1968–75) (A) and the neuraminidase N2 gene (1967–72) (B) of influenza virus. For the hemagglutinin gene, black dots indicate the antigenic clusters (A) as defined elsewhere [1]. For the neuraminidase gene, black dots denote the genetic clusters corresponding to the 1968/1969 and 1969/1970 sequences (B); clusters were supported by bootstrap values >90%. Braces indicate earlier or later clusters. Italic font indicates neuraminidase genes from A/H2N2 viruses; boldface font indicates neuraminidase genes from A/H3N2 viruses. Hemagglutinin and neuraminidase sequences (nt 49–960 and 1–1310, respectively) were aligned by use of the Clustal W program of MacVector (version 7.2.2). Phylogenetic analysis was performed by use of the neighbor-joining method, with Kimura’s 2-parameter distance (1000 bootstrap replicates). A/Hong-Kong/1/68 was used as the root for the hemagglutinin tree, and A/Singapore/1/57 was used as the root for the neuraminidase tree. The first A/H3N2 virus isolated was A/Hong Kong/1/68

Figure 4

Phylogenetic relationships of the hemagglutinin H3 gene (1968–75) (A) and the neuraminidase N2 gene (1967–72) (B) of influenza virus. For the hemagglutinin gene, black dots indicate the antigenic clusters (A) as defined elsewhere [1]. For the neuraminidase gene, black dots denote the genetic clusters corresponding to the 1968/1969 and 1969/1970 sequences (B); clusters were supported by bootstrap values >90%. Braces indicate earlier or later clusters. Italic font indicates neuraminidase genes from A/H2N2 viruses; boldface font indicates neuraminidase genes from A/H3N2 viruses. Hemagglutinin and neuraminidase sequences (nt 49–960 and 1–1310, respectively) were aligned by use of the Clustal W program of MacVector (version 7.2.2). Phylogenetic analysis was performed by use of the neighbor-joining method, with Kimura’s 2-parameter distance (1000 bootstrap replicates). A/Hong-Kong/1/68 was used as the root for the hemagglutinin tree, and A/Singapore/1/57 was used as the root for the neuraminidase tree. The first A/H3N2 virus isolated was A/Hong Kong/1/68

Phylogenetic analysis of the neuraminidase N2 gene highlighted similarities in this gene between late A/H2N2 and early A/H3N2 viruses (figure 4B). The neuraminidase gene of A/H2N2 viruses from 1967/1968 formed 2 distinct clusters (clusters I and II in figure 4B). The neuraminidase gene of A/H3N2 viruses from 1968/1969 formed a single cluster, genetically close to the A/H2N2 cluster I. Although few sequences were available for A/H3N2 viruses from 1969/1970, they all formed a distinct cluster comprising strains from England and Canada, genetically close to the A/H2N2 cluster II. For A/H3N2 viruses, the 1968/1969 genetic cluster differed from the 1969/1970 genetic cluster by 29 nt. The clusters were supported by bootstrap values >90%. The neuraminidase genes of more-recent A/H3N2 viruses were derived from the 1969/1970 genetic cluster. In addition, we predicted and compared the protein sequence of neuraminidase. Similar to the nucleotide analysis, the protein analysis revealed 2 distinct clusters for 1968/1969 and 1969/1970. The clusters differed by 11 amino acid changes, 3 of them located in characterized antigenic sites [38]. These findings are in agreement with those of a recent study showing that neuraminidase had an increased rate of change immediately after the emergence of A/H3N2 viruses, compared with that during the more-recent interpandemic period [38]. These results highlight temporal differences in the neuraminidases of influenza viruses circulating during the 2 pandemic seasons, whereas the hemagglutinins remained unchanged

Discussion

The present study is based on data from 6 countries and demonstrates distinct differences in mortality patterns for the Hong Kong A/H3N2 pandemic. The United States and Canada displayed the expected pattern—high mortality during the first pandemic season (1968/1969), when the emerging virus first circulated, followed by a relatively mild second pandemic season (1969/1970). This pattern was observed only in North America. For the 4 other countries in Europe and Asia that were studied, we identified an opposite mortality pattern: more than two-thirds of all influenza-related deaths that occurred during 1968/1969–1969/1970 occurred during the second pandemic season. Although this smoldering pattern had previously been reported for England [9], the present study is the first multinational study to carefully compare and quantify the mortality and temporal pattern of the Hong Kong A/H3N2 pandemic. Published mortality time series for at least 8 additional countries in Europe and elsewhere support the idea that the 1-year delay in mortality might be the most common experience in continents other than North America [7, 31, 39–41 ]

We next entertained a possible hypothesis to explain the unexpected smoldering pattern in Europe and Asia (figure 5) that combines (1) the effect of geographical differences in preexisting immunity to neuraminidase at the time of emergence of A/H3N2 (remaining from the A/H2N2 era) and (2) the effect of genetic drift in the neuraminidase antigen during 1969/1970. Preexisting immunity to the neuraminidase N2 gene may have contributed to the differential impact of the first pandemic season in the 2 regions. It is believed that exposure to A/H2N2 viruses during 1967/1968 attenuated the impact of the 1968 pandemic [3, 11]. A study showed that individuals infected by A/H2N2 during 1967/1968 were protected against influenza infection during the first A/H3N2 pandemic season, and infections in persons with prior antineuraminidase antibodies were more frequently asymptomatic, compared with those in persons without prior antineuraminidase antibodies [20]

Figure 5

Hypotheses for explaining the differences in influenza pandemic patterns between Europe and Asia and North America. The genetic evolution of hemagglutinin and neuraminidase are indicated from the first pandemic season of A/H2N2 viruses (1957/1958) to the second pandemic season of A/H3N2 viruses (1969/1970). * and †, successive drifts; •, virus circulation and impact in each region (•, low; ••, high)

Figure 5

Hypotheses for explaining the differences in influenza pandemic patterns between Europe and Asia and North America. The genetic evolution of hemagglutinin and neuraminidase are indicated from the first pandemic season of A/H2N2 viruses (1957/1958) to the second pandemic season of A/H3N2 viruses (1969/1970). * and †, successive drifts; •, virus circulation and impact in each region (•, low; ••, high)

The smoldering pattern in Europe and Asia is consistent with high preexisting immunity to neuraminidase. The greater mortality of the last A/H2N2 epidemic season in this region suggests greater exposure to late A/H2N2 viruses (figure 5). Hence, immunity to neuraminidase may have been greater in Europe and Asia during the fall of 1968, when the A/H3N2 virus started to circulate. Furthermore, the neuraminidase of late A/H2N2 viruses circulating in Europe and Asia may have been closer to that of early A/H3N2 viruses [14, 16, 42]

A drift in the neuraminidase antigen found in our phylogenetic analysis coincided with the severe second A/H3N2 pandemic season in Europe and Asia (figure 5). This analysis is in agreement with antigenic studies conducted by the World Influenza Center (London, UK) in Europe and Japan. In 1968/1969, influenza viruses were antigenically similar to the first A/H3N2 virus identified in Hong Kong, in both hemagglutinin and neuraminidase [42, 43]. In 1969/1970, the neuraminidase antigen had changed, but the hemagglutinin antigen had not [41–43 ]. Accordingly, Lindstrom et al. recently found that 2 distinct neuraminidase lineages circulated during 1968/1969 and 1969/1970 and hypothesized that such genetic diversity could be due to multiple reassortments with A/H2N2 viruses [38]

We were unfortunately unable to study geographical heterogeneity in the neuraminidase gene, since very few North American viral sequences from 1968/1969 to 1969/1970 were available in the public domain. Rapid fluxes in international populations [8] and the available sequences would suggest that the same viruses circulated in North America and in Europe and Asia

Assuming that the same viruses circulated in all continents studied, differences in prior immunity followed by a drift in the neuraminidase antigen could explain the intercountry differences in mortality patterns. What sets North America apart epidemiologically is the severity of the first pandemic season, which included frequent influenza illness and death (figure 5). Hence, at the end of the first pandemic season, a substantial proportion of the North American population had antibodies to the novel hemagglutinin H3 gene. Since the hemagglutinin did not change between the 2 pandemic seasons, these antibodies would have protected this population during the second pandemic season

One of the other factors we considered was immunity to the hemagglutinin H3 gene in the very elderly, which remained from childhood exposure to H3-like viruses before 1892 [10, 32, 44–46 ]. However, such immunity would not account for the intercountry differences in mortality found in persons <65 years old. Furthermore, influenza vaccination coverage was too limited during this time to explain these differences [9]. Finally, environmental differences and weather could have played a role—but this explanation is unlikely, given the geographical heterogeneity of the countries studied [47]. The A/H3N2 pandemic strain appeared early during the winter season in all countries. Certainly, weather did not hamper the full potential for causing widespread epidemics and mortality in every country during the 1968/1969 pandemic season

Despite the geographical differences in the timing of the major mortality, the A/H3N2 pandemic was relatively mild in all countries, compared with surrounding severe epidemics [48], including the last A/H2N2 epidemic season. The mildness of the 1968 pandemic is perhaps not entirely unexpected, considering preexisting immunity to the neuraminidase antigen in all age groups [3, 11] and to the hemagglutinin antigen in the elderly [10, 32, 44–46 ]

It has been reported that, for the 1950s and 1960s, both P&I excess mortality and all-cause crude excess mortality (unadjusted) were systematically 2–3-fold higher in Europe than in North America [27, 28, 31]. Adjusting for differences in age structure and mortality outside influenza periods (data not shown) reduced but did not eliminate these differences. Therefore, we studied the percentage increase in mortality above the winter baseline, which further reduced intercountry differences [29]. A more thorough investigation of the intercountry differences in absolute mortality was beyond the scope of the present study; factors such as population density, demographics, climate, country size, and health care systems could affect influenza transmission and disease outcome and, in turn, mortality [27, 28, 31, 49]. In the present study, the smoldering pattern relies on the unbiased comparison of the relative impact of each pandemic season

The present study was based on analysis of both P&I (as an underlying cause of death) and all-cause mortality. Analysis of P&I data could potentially have been biased as a result of intercountry differences in coding of the underlying cause of death and the transition to a new system of classification of deaths around 1968 (Appendix). However, because our analysis of all-cause mortality data was not affected by such possible bias and because it produced similar results in terms of temporal and geographical differences in patterns, nosological bias was not an issue in the present study

Another limitation of the present study was the poor temporal and geographical resolution of influenza genetic sequences from this pandemic, both in terms of the overall number of sequences and the lack of systematic sampling (outliers were preferentially sequenced [1]). This bias tends to preclude the use of phylogenetic studies to derive the epidemiological importance of each strain. However, in this work, analyses at the amino acid level and antigenic studies support the significance of the drift in the neuraminidase antigen [38, 41–43 ]

Why should we care about studying the geographical and temporal impact of pandemic influenza? The pattern in Europe and Asia during the 1968/1969–1969/1970 period suggests a possible favorable opportunity for pandemic response should a future pandemic be like the one in 1968. The smoldering first pandemic season observed in 4 of the 6 countries suggests that a pandemic vaccine available 1 year after the emergence of the new subtype could have prevented the majority of deaths and illnesses associated with the emerging pandemic strain. The exact role that neuraminidase played in driving the mortality of the 1968 pandemic remains to be confirmed, but it seems that the 1968 pattern, unique in its conservation of the neuraminidase antigen, allowed ample time to produce and distribute a pandemic vaccine. In the meantime, vaccination targeted against neuraminidase could be complementary to treatment and prophylaxis by antiviral agents, which would probably be available in short supply in most countries during a pandemic situation [50]. Unfortunately, given our limited experience with pandemic influenza—only 3 pandemics occurred during the 20th century—our ability to predict the likelihood of a 1968-like pandemic in the future is quite limited

Multinational Influenza Seasonal Mortality Study (MISMS) Group

MISMS group members (other than the authors already listed) are as follows: J. Dushoff (Princeton University, Princeton), D. M. Fleming (Birmingham Research Unit of the Royal College of General Practitioners, Birmingham, UK), A. W. Hampson (World Health Organization Collaborating Centre for Reference and Research on Influenza, Melbourne, Australia), N. Sugaya (Keiyu Hospital, Yokohama, Japan), T. W. Tam (Health Canada, Ottawa, Canada), A. J. Valleron (Institut National de la Santé et de la Recherche Médicale, Paris, France), and J. Watkins (University of Wales College of Medicine, Cardiff, UK)

Acknowledgments

We thank O. T. Bjørnstad, B. T. Grenfell, D. L. Smith, F. E. McKenzie, and J. Dushoff for fruitful discussions on this study

References

1.
Smith
DJ
Lapedes
AS
de Jong
JC
, et al. 
Mapping the antigenic and genetic evolution of influenza virus
Science
2004
, vol. 
305
 (pg. 
371
-
6
)
2.
Nicholson
K
Hay
A
Textbook of influenza
1998
Oxford
Blackwell
3.
Kilbourne
ED
Perspectives on pandemics: a research agenda
J Infect Dis
1997
, vol. 
176
 
Suppl 1
(pg. 
S29
-
31
)
4.
Nguyen‐Van‐Tam
JS
Hampson
AW
The epidemiology and clinical impact of pandemic influenza
Vaccine
2003
, vol. 
21
 (pg. 
1762
-
8
)
5.
Noble
G
Bear
A
Epidemiological and clinical aspects of influenza
Basic and applied influenza research
1982
Boca Raton
CRC Press
6.
Simonsen
L
Clarke
MJ
Schonberger
LB
Arden
NH
Cox
NJ
Fukuda
K
Pandemic versus epidemic influenza mortality: a pattern of changing age distribution
J Infect Dis
1998
, vol. 
178
 (pg. 
53
-
60
)
7.
Cockburn
WC
Delon
PJ
Ferreira
W
Origin and progress of the 1968–69 Hong Kong influenza epidemic
Bull World Health Organ
1969
, vol. 
41
 (pg. 
345
-
8
)
8.
Longini
IM
Jr
Fine
PE
Thacker
SB
Predicting the global spread of new infectious agents
Am J Epidemiol
1986
, vol. 
123
 (pg. 
383
-
91
)
9.
Miller
DL
Pereira
MS
Clarke
M
Epidemiology of the Hong Kong‐68 variant of influenza A2 in Britain
Br Med J
1971
, vol. 
1
 (pg. 
475
-
9
)
10.
Stuart‐Harris
CH
Pandemic influenza: an unresolved problem in prevention
J Infect Dis
1970
, vol. 
122
 (pg. 
108
-
15
)
11.
Stuart‐Harris
C
Epidemiology of influenza in man
Br Med Bull
1979
, vol. 
35
 (pg. 
3
-
8
)
12.
Cox
NJ
Subbarao
K
Global epidemiology of influenza: past and present
Annu Rev Med
2000
, vol. 
51
 (pg. 
407
-
21
)
13.
Sharrar
RG
National influenza experience in the USA, 1968–69
Bull World Health Organ
1969
, vol. 
41
 (pg. 
361
-
6
)
14.
Cameron
AS
Roder
DM
Esterman
AJ
Moore
BW
Mortality from influenza and allied infections in South Australia during 1968–1981
Med J Aust
1985
, vol. 
142
 (pg. 
14
-
7
)
15.
Fukumi
H
Summary report on Hong Kong influenza in Japan
Bull World Health Organ
1969
, vol. 
41
 (pg. 
353
-
9
)
16.
Roden
AT
National experience with Hong Kong influenza in the United Kingdom, 1968–69
Bull World Health Organ
1969
, vol. 
41
 (pg. 
375
-
80
)
17.
Ferris
AA
Influenza in Melbourne, Australia, 1969
Bull World Health Organ
1969
, vol. 
41
 pg. 
399
 
18.
Influenza A2/Hong‐Kong/68 in the world
Wkly Epidemiol Rec
1969
, vol. 
44
 (pg. 
187
-
8
)
19.
Influenza—worldwide 1969–70
MMWR Morb Mortal Wkly Rep
1970
, vol. 
19
 (pg. 
105
-
6
)
20.
Monto
AS
Kendal
AP
Effect of neuraminidase antibody on Hong Kong influenza
Lancet
1973
, vol. 
1
 (pg. 
623
-
5
)
21.
Monto
AS
Davenport
FM
Napier
JA
Francis
T
Jr
Modification of an outbreak of influenza in Tecumseh, Michigan by vaccination of schoolchildren
J Infect Dis
1970
, vol. 
122
 (pg. 
16
-
25
)
22.
Monto
AS
Kioumehr
F
The Tecumseh study of respiratory illness. IX. Occurrence of influenza in the community, 1966–1971
Am J Epidemiol
1975
, vol. 
102
 (pg. 
553
-
63
)
23.
Monto
AS
Cavallaro
JJ
The Tecumseh study of respiratory illness. II. Patterns of occurrence of infection with respiratory pathogens, 1965–1969
Am J Epidemiol
1971
, vol. 
94
 (pg. 
280
-
9
)
24.
Davignon
L
Williams
MC
Surveillance of influenza in Quebec 1968–71
Can J Public Health
1974
, vol. 
65
 (pg. 
284
-
91
)
25.
Macken
C
Lu
H
Goodman
J
Boykin
L
Osterhaus
AD
Cox
N
Hampson
AW
The value of a database in surveillance and vaccine selection
Options for the control of influenza IV
2001
Amsterdam
Elsevier
(pg. 
103
-
6
)
26.
Serfling
R
Methods for current statistical analysis of excess pneumonia‐influenza deaths
Public Health Rep
1963
, vol. 
78
 (pg. 
494
-
506
)
27.
Langmuir
AD
Housworth
J
A critical evaluation of influenza surveillance
Bull World Health Organ
1969
, vol. 
41
 (pg. 
393
-
8
)
28.
World Health Organization
Mortality at ages 65 and older—United States, Canada, and western Europe
Stat Bull
1971
, vol. 
52
 (pg. 
8
-
11
)
29.
Serfling
RE
Sherman
IL
Houseworth
WJ
Excess pneumonia‐influenza mortality by age and sex in three major influenza A2 epidemics, United States, 1957–58, 1960 and 1963
Am J Epidemiol
1967
, vol. 
86
 (pg. 
433
-
41
)
30.
Stroup
DF
Thacker
SB
Herndon
JL
Application of multiple time series analysis to the estimation of pneumonia and influenza mortality by age 1962–1983
Stat Med
1988
, vol. 
7
 (pg. 
1045
-
59
)
31.
Assaad
F
Cockburn
WC
Sundaresan
TK
Use of excess mortality from respiratory diseases in the study of influenza
Bull World Health Organ
1973
, vol. 
49
 (pg. 
219
-
33
)
32.
Simonsen
L
Reichert
TA
Miller
MA
Kawaoka
Y
The virtues of antigenic sin: consequences of pandemic recycling on influenza‐associated mortality
Options for the control of influenza V
2004
Amsterdam
Elsevier
(pg. 
791
-
4
)
33.
Housworth
WJ
Spoon
MM
The age distribution of excess mortality during A2 Hong Kong influenza epidemics compared with earlier A2 outbreaks
Am J Epidemiol
1971
, vol. 
94
 (pg. 
348
-
50
)
34.
Alling
DW
Blackwelder
WC
Stuart‐Harris
CH
A study of excess mortality during influenza epidemics in the United States, 1968–1976
Am J Epidemiol
1981
, vol. 
113
 (pg. 
30
-
43
)
35.
Clifford
RE
Smith
JW
Tillett
HE
Wherry
PJ
Excess mortality associated with influenza in England and Wales
Int J Epidemiol
1977
, vol. 
6
 (pg. 
115
-
28
)
36.
Influenza in the world
Wkly Epidemiol Rec
1969
, vol. 
44
 (pg. 
12
-
27
)
37.
Sohier
R
Henry
M
Epidemiological data on Hong Kong influenza in France
Bull World Health Organ
1969
, vol. 
41
 (pg. 
402
-
4
)
38.
Lindstrom
SE
Cox
NJ
Klimov
A
Genetic analysis of human H2N2 and early H3N2 influenza viruses, 1957–1972: evidence for genetic divergence and multiple reassortment events
Virology
2004
, vol. 
328
 (pg. 
101
-
19
)
39.
Rocchi
G
Ragona
G
De Felici
A
Muzzi
A
Epidemiological evaluation of influenza in Italy
Bull World Health Organ
1974
, vol. 
50
 (pg. 
401
-
6
)
40.
Sprenger
MJ
Mulder
PG
Beyer
WE
Van Strik
R
Masurel
N
Impact of influenza on mortality in relation to age and underlying disease, 1967–1989
Int J Epidemiol
1993
, vol. 
22
 (pg. 
334
-
40
)
41.
Werner
J
Schudrowitz
C
Kohler
H
Antigenic variation of neuraminidase of human type A influenza (H3N2) viruses isolated in Berlin (West)
Zentralbl Bakteriol [Orig A]
1975
, vol. 
233
 (pg. 
440
-
6
)
42.
Schild
GC
Henry‐Aymard
M
Pereira
MS
, et al. 
Antigenic variation in current human type A influenza viruses: antigenic characteristics of the variants and their geographic distribution
Bull World Health Organ
1973
, vol. 
48
 (pg. 
269
-
78
)
43.
Aymard‐Henry
M
Schild
GC
Newman
RW
Regamey
RH
Perkins
FT
Antigenic variation in the neuraminidases of human influenza A strains: evidence of antigenic shift in neuraminidase predating changes in hemagglutinin
Influenza vaccine for men and horses
1973
Basel
Karger
(pg. 
28
-
30
)
44.
Davenport
FM
Minuse
E
Hennessy
AV
Francis
T
Jr
Interpretations of influenza antibody patterns of man
Bull World Health Organ
1969
, vol. 
41
 (pg. 
453
-
60
)
45.
Masurel
N
Marine
WM
Recycling of Asian and Hong Kong influenza A virus hemagglutinins in man
Am J Epidemiol
1973
, vol. 
97
 (pg. 
44
-
9
)
46.
Fukumi
H
Interpretation of influenza antibody patterns in man: existence and significance of Hong Kong antibody in old people prior to the Hong Kong influenza epidemic
Bull World Health Organ
1969
, vol. 
41
 (pg. 
469
-
73
)
47.
IRI/LDEO Climate Data Library
 
Available at: http://ingrid.ldeo.columbia.edu. Accessed 5 August 2004
48.
Simonsen
L
Olson
D
Viboud
C
Miller
M
Knobler
S
Oberholtzer
K
Pandemic influenza and mortality: past evidence and projections for the future
Forum on microbial threats. Pandemic influenza: assessing capabilities for prevention and response
2004
Washington, DC
Institute of Medicine, The National Academy of Sciences
49.
Reichert
TA
Sharma
A
Prado
S
Osterhaus
AD
Cox
N
Hampson
AW
A global pattern for influenza activity
Options for the control of influenza IV
2001
Amsterdam
Elsevier
(pg. 
87
-
94
)
50.
Kaiser
J
Influenza: girding for disaster. Facing down pandemic flu, the world’s defenses are weak
Science
2004
, vol. 
306
 (pg. 
394
-
7
)
51.
Gill
PW
Babbage
NF
Gunton
PE
Flower
W
Garrett
DA
Did the Asian virus protect us from Hong Kong influenza?
Med J Aust
1971
, vol. 
2
 (pg. 
53
-
4
)
52.
Simonsen
L
Reichert
TA
Viboud
C
Blackwelder
WC
Miller
MA
Impact of influenza vaccination on seasonal mortality in the US elderly
Arch Intern Med
2005
, vol. 
165
 (pg. 
265
-
72
)
53.
Foy
HM
Cooney
MK
Allan
I
Kenny
GE
Rates of pneumonia during influenza epidemics in Seattle, 1964 to 1975
JAMA
1979
, vol. 
241
 (pg. 
253
-
8
)

Figures and Tables

Table A1

Average intercountry differences in age-adjusted excess mortality related to influenza

Table A1

Average intercountry differences in age-adjusted excess mortality related to influenza

Table A2

Change in summer pneumonia and influenza (P&I) mortality at the transition of International Classification of Diseases (ICD) and in later years, by country (1967–1978)

Table A2

Change in summer pneumonia and influenza (P&I) mortality at the transition of International Classification of Diseases (ICD) and in later years, by country (1967–1978)

Financial support: National Institutes of Health

R.F.G. is under contract to Fogarty International Center

Study group members are listed after the text