MTBI Mathematical and Theoretical Biology Institute Modeling

MTBI Mathematical and Theoretical Biology Institute

Modeling Influenza Epidemics and Related Issues Carlos Castillo-Chavez Department of Biological Statistics and Computational Biology Department of Theoretical and Applied Mechanics Current Collaborators (Direct or Indirect) Miriam Nuno, Cornell University Zhilan Feng, Purdue University Maia Martcheva, Polytechnic University Laura Jones, Cornell University

Influenza virus Picture

Historical Overview History suggests that influenza pandemics have probably happened during at least the last four centuries. During the 20 th century, three pandemics and several "pandemic scares" occurred. 1918: Spanish Flu The Spanish Influenza pandemic is the catastrophe against which all modern pandemics are measured. It is estimated that approximately 20 to 40 percent of the worldwide population became ill and that over 20 million people died. Between September 1918 and April 1919, approximately 500, 000 deaths from the flu occurred in the U. S. alone. Many people died from this very quickly. Some people who felt well in the morning became sick by noon, and were dead by nightfall. Those who did not succumb to the disease within the first few days often died of complications from the flu (such as pneumonia) caused by bacteria. One of the most unusual aspects of the Spanish flu was its ability to kill young adults. The reasons for this remain uncertain. With the Spanish flu, mortality rates were high among healthy adults as well as the usual high-risk groups. The attack rate and mortality was highest among adults 20 to 50 years old. The severity of that virus has not been seen again.

Historical Overview (Taken from Web Site) History suggests that influenza pandemics have probably happened during at least the last four centuries. 20 th century pandemics and "pandemic scares. " 1918: Spanish Flu It is estimated that approximately 20 to 40 percent of the worldwide population became ill and that over 20 million people died. Between September 1918 and April 1919, approximately 500, 000 deaths from the flu occurred in the U. S. alone. Some people who felt well in the morning became sick by noon, and were dead by nightfall. Those who did not succumb to the disease within the first few days often died of complications from the flu (such as pneumonia) caused by bacteria. With the Spanish flu, mortality rates were high among healthy adults as well as the usual high-risk groups. The attack rate and mortality was highest among adults 20 to 50 years old. The severity of that virus has not been seen again.

1957: Asian Flu In February 1957, the Asian influenza pandemic was first identified in the Far East. Immunity to this strain was rare in people less than 65 years of age, and a pandemic was predicted. The 1957 pandemic virus was quickly identified, due to advances in scientific technology. Vaccine was available in limited supply by August 1957. The virus came to the U. S. quietly, with a series of small outbreaks over the summer of 1957. When U. S. children went back to school in the fall, they spread the disease in classrooms and brought it home to their families. Infection rates were highest among school children, young adults, and pregnant women in October 1957. Most influenza-and pneumonia-related deaths occurred between September 1957 and March 1958. The elderly had the highest rates of death. During January and February 1958, there was another wave of illness among the elderly. This is an example of the potential "second wave" of infections that can develop during a pandemic. The disease infects one group of people first, infections appear to decrease and then infections increase in a different part of the population. Although the Asian flu pandemic was not as devastating as the Spanish flu, about 69, 800 people in the U. S. died.

1968: Hong Kong Flu In early 1968, the Hong Kong influenza pandemic was first detected in Hong Kong. The first cases in the U. S. were detected as early as September of that year, but illness did not become widespread in the U. S. until December. Deaths from this virus peaked in December 1968 and January 1969. Those over the age of 65 were most likely to die. The same virus returned in 1970 and 1972. The number of deaths between September 1968 and March 1969 for this pandemic was 33, 800, making it the mildest pandemic in the 20 th century.

1976: Swine Flu Scare When a novel virus was first identified at Fort Dix, it was labeled the "killer flu. " Experts were extremely concerned because the virus was thought to be related to the Spanish flu virus of 1918. The concern that a major pandemic could sweep across the world led to a mass vaccination campaign in the United States. In fact, the virus--later named "swine flu"--never moved outside the Fort Dix area. Research on the virus later showed that if it had spread, it would probably have been much less deadly than the Spanish flu.

1977: Russian Flu Scare In May 1977, influenza A/H 1 N 1 viruses isolated in northern China, spread rapidly, and caused epidemic disease in children and young adults (< 23 years) worldwide. The 1977 virus was similar to other A/H 1 N 1 viruses that had circulated prior to 1957. (In 1957, the A/H 1 N 1 virus was replaced by the new A/H 2 N 2 viruses). Because of the timing of the appearance of these viruses, persons born before 1957 were likely to have been exposed to A/H 1 N 1 viruses and to have developed immunity against A/H 1 N 1 viruses. Therefore, when the A/H 1 N 1 reappeared in 1977, many people over the age of 23 had some protection against the virus and it was primarily younger people who became ill from A/H 1 N 1 infections. By January 1978, the virus had spread around the world, including the United States. Because illness occurred primarily in children, this event was not considered a true pandemic. Vaccine containing this virus was not produced in time for the 1977 -78 season, but the virus was included in the 1978 -79 vaccine.

1997: Avian Flu Scare The most recent pandemic "scares" occurred in 1997 and 1999. In 1997, at least a few hundred people became infected with the avian A/H 5 N 1 flu virus in Hong Kong and 18 people were hospitalized. Six of the hospitalized persons died. This virus was different because it moved directly from chickens to people, rather than having been altered by infecting pigs as an intermediate host. In addition, many of the most severe illnesses occurred in young adults similar to illnesses caused by the 1918 Spanish flu virus. To prevent the spread of this virus, all chickens (approximately 1. 5 million) in Hong Kong were slaughtered. The avian flu did not easily spread from one person to another, and after the poultry slaughter, no new human infections were found. In 1999, another novel avian flu virus – A/H 9 N 2 – was found that caused illnesses in two children in Hong Kong. Although both of these viruses have not gone on to start pandemics, their continued presence in birds, their ability to infect humans, and the ability of influenza viruses to change and become more transmissible among people is an ongoing concern.

Influenza and Cross-Immunity OUTLINE Ø Epidemiology of Influenza Ø Ø Ø Ø Cross-Immunity Two strain models (CHALL) Age structure Results Limitations Example Conclusions Future Work

Motivation ØResearchers have explored the possible mechanism(s) underlying the recurrence of epidemics and persistence of co-circulating virus strains of influenza types between pandemics. ØIt is believed that interactions between human populations with host (animal) populations such as swine, duck, equine may be responsible for the observed epidemic and pandemic outbreaks.

Evidence of Cross-Immunity ØStudy in 1979 shows that only 3% of individuals previously exposed to strain A/Hong Kong/68 or A/England/72 were infected with a similar subtype H 3 N 2, while 23% were infected without a previous exposure. ØFrequency of detection of antibody-positive sera between 1977 and 1978 changed from 0% to 38% for young people. (reappearance of H 1 N 1 subtype). ØHouston study indicates that no cross-reactive immunity exist between subtypes H 1 N 1 and H 3 N 2.

What is Cross-Immunity? ØThere are three subtypes of influenza A: H 1 N 1, H 2 N 2, and H 3 N 2 ØInfection with an influenza subtype A strain may provide cross protection against other antigenically similar circulating strains.

Cross Immunity n n n -coefficient of cross-immunity Relative reduction on susceptibility due to prior exposure to a related strain. =0, represents total cross-immunity =1, represents no cross-immunity 0< <1, represents partial cross-immunity >1, represents immune deficiency

EXPERIMENTAL EVIDENCE (1) 1974 Study < 3% With Prior (2) A/HONG KONG/68 (H 3 N 2) OR (3) A PRIOR A/ENGLAND/72 (H 3 N 2) (4) GOT (5) A/PORT CHALMERS/73 (6) Vs. 23% With NO Prior Experience. (7) (2) 1976 Appearance of (8) A/VICTORIA/75 (H 3 N 2) (9) Relative Frequency of First Infected/Previously Infected (By Another Strain of H 3 N 2( Was Approx. 41%).

EXPERIMENTAL EVIDENCE (3) 1977 Co-Circulating H 1 N 2 Strains (4) Individuals bon before 1952 “GOT” a strain of H 1 N 1 (5) DETECTION OF ANTIBODY-POSITIVE SERA (6) YOUNG: Changed From 0% to 9%. (7) OLDER: Did not change (remained at 9%). (8) (4) 1982 (Glezen) No Cross-Immunity Between Subtypes (H 1 N 1 & H 3 N 2)

Facts Experimental results indicate that Cross-immunity shares the following features (Couch and Kasel, 1983): ØExhibits subtype specificity. ØExhibits cross-reactivity to variants within a subtype, but with reduced cross-reactivity for variants that are antigenically distant from the initial variant. ØExhibits a duration of at least five to eight years. ØBe able to account for the observation that resistance to reinfection with H 1 N 1 may last 20 years

First Approaches ØIn 1975 epidemiological interference of virus populations was introduced [Dietz]. ØIn 1989 age-structure, proportionate mixing and crossimmunity are studied [Castillo-Chavez, et. al]. ØIn 1989 interactions between human and animal host populations are studied as a source of recombinants in strains and cross-immunity.

Models with age structure

SIR Model with Age Structure q s(t, a) : Density of susceptible individuals with age a at time t. q i(t, a) : Density of infectious individuals with age a at time t. q r(t, a) : Density of recovered individuals with age a at time t. # of susceptible individuals with ages in (a 1 , a 2) at time t # of infectious individuals with ages in (a 1 , a 2) at time t # of recovered individuals with ages in (a 1 , a 2) at time t

Parameters q q q : recruitment/birth rate. (a): age-specific probability of becoming infected. c(a): age-specific per-capita contact rate. (a): age-specific per-capita mortality rate. (a): age-specific per-capita recovery rate.

Equations

Initial and Boundary Conditions

Demographic Steady State n(t, a): density of individual with age a at time t n(t, a) satisfies the Mackendrick Equation We assume that the total population density has reached this demographic steady state.

Mixing Rules 1. p(t, a, a`) 0 2. 3. 4. 5. Proportionate mixing:

Stability of Disease-free Steady State

Characteristic Equation The characteristic equation has a unique real solution which is dominant, that is, the real part of rest solutions is less than This dominant solution.

R 0<1, Disease-free State Is Stable q The characteristic equation has a unique dominant real solution. That is, the real part of all other solutions is less than this dominant solution; q The dominant solution is negative iff R 0<1; q Whenever R 0<1, the disease-free steady state is locally asymptotically stable.

Endemic Steady States

Endemic Steady States Moving forward, one can formally solve for the steady states. The existence of endemic steady states is determined by the roots of the following equation: f(B*) is a decreasing function of B* with f( )=0. R 0>1, there exists a unique endemic (e. g. non trivial) steady states; R 0<1, an endemic steady state does not exist.

Cross Immunity n n n -coefficient of cross-immunity Relative reduction on susceptibility due to prior exposure to a related strain. =0, represents total cross-immunity =1, represents no cross-immunity 0< <1, represents partial cross-immunity >1, represents immune deficiency

Transfer Diagram: X Y 1 Z 1 Y 2 V 2 Z 2 V 1 W NOTATION: X Susceptible Individuals Yi Individuals Infected Strain i. Zi Individuals Recovered Strain i. But Susceptible to Other Strain. Vi Individuals Infected by Strain I but Recovered From the Other Strain. W Recovered Individuals.

Hopf-Bifurcation point ì > 0 a <a. C 1 ï (a 1 ) = í < 0 a 1 > a 1 C w 2 ï = Hopf - Bifurcation î a 1 I 1

Hopf -Bifurcation 2 Strain 2 E Strain 1 2 Strain 2 is present E 0 E 1 1 Strain 1 is present *We have oscillations for a range of the quarantine period

Results ØAge-structure is sufficient to drive sustained oscillations in a multiple strain model[Castillo-Chavez, Hethcote, Andreasen, Liu and Levin]. ØFor a heterogeneous population and age-dependent mortality, cross-immunity provides an explanation to the observed recurrence of strains [Castillo-Chavez, et. al]. ØUnder the assumption of cross-immunity and multiple strains, quarantine periods support sustained oscillations [Castillo-Chavez, et. al].

A model on the role of cross-immunity and quarantine ØIn the following model we explore crossimmunity and quarantine periods as the mechanisms driving the observed epidemic outbreaks. ØMultiple outbreaks appear and grow in number as we increase the levels of cross protection for the circulating type A influenza strains.

m. I m. S S b S(I + V 1 1 1 ) A 1 2 d. I I m. Q 1 g. I 1 1 Q d. I 2 2 2 m. Q a Q 2 2 2 V 2 g. I R 1 2 2 1 1 V ) m. V 1 1 g. V 1 1 W 2 2 A m. V 2 g V 2 b s. R (I + V A m. R 2 1 1 2 I 2 m. R b s R (I + V ) 2 2 1 R 1 A m. I 1 a. Q 1 b S(I + V ) 2 Q 2 m. W 2

The Model b S(I + V ) b S(I + V d. S = dt A A ( d. I b S I + V ) (m = - + g + d )Ii dt A i i i j j j )- m i i i d. Qi = d i Ii - ( m + a i ) Q i dt (I + V d. R = g I + a Q - m. R - b s. R dt A (I + V ) ( d. V = b s. R - m +g )V dt A i j i i i i j i d. W = g V + g V - m. W dt i i i j j i j ) S

Hopf-Bifurcation point ì > 0 a <a. C 1 ï (a 1 ) = í < 0 a 1 > a 1 C w 2 ï = Hopf - Bifurcation î a 1 I 1

Hopf -Bifurcation 2 Strain 2 E Strain 1 2 Strain 2 is present E 0 E 1 1 Strain 1 is present *We have oscillations for a range of the quarantine period

Conclusions n n Cross-immunity allows for the incorporation of immune system level phenomena into population models. Population structure is critical in the study of communicable diseases. The impact of reservoir of infection needs further study. Identification of mechanisms capable of generating sustained oscillations is important.

Future Work n n The study of influenza and other communicable diseases at global scales is critical (several Russian researchers and Longini have done some important past work). The impact of continued population growth and environmental changes also must be considered in the study of disease spread.