Hiroshi Nishiura

Comparative estimation of the reproduction number for pandemic influenza from daily case notification data

The reproduction number, , defined as the average number of secondary cases generated by a primary case, is a crucial quantity for identifying the intensity of interventions required to control an epidemic. Current estimates of the reproduction number for seasonal influenza show wide variation and, in particular, uncertainty bounds for for the pandemic strain from 1918 to 1919 have been obtained only in a few recent studies and are yet to be fully clarified. Here, we estimate using daily case notifications during the autumn wave of the influenza pandemic (Spanish flu) in the city of San Francisco, California, from 1918 to 1919. In order to elucidate the effects from adopting different estimation approaches, four different methods are used: estimation of using the early exponential-growth rate (Method 1), a simple susceptible–exposed–infectious–recovered (SEIR) model (Method 2), a more complex SEIR-type model that accounts for asymptomatic and hospitalized cases (Method 3), and a stochastic susceptible–infectious–removed (SIR) with Bayesian estimation (Method 4) that determines the effective reproduction number at a given time t. The first three methods fit the initial exponential-growth phase of the epidemic, which was explicitly determined by the goodness-of-fit test. Moreover, Method 3 was also fitted to the whole epidemic curve. Whereas the values of obtained using the first three methods based on the initial growth phase were estimated to be 2.98 (95% confidence interval (CI): 2.73, 3.25), 2.38 (2.16, 2.60) and 2.20 (1.55, 2.84), the third method with the entire epidemic curve yielded a value of 3.53 (3.45, 3.62). This larger value could be an overestimate since the goodness-of-fit to the initial exponential phase worsened when we fitted the model to the entire epidemic curve, and because the model is established as an autonomous system without time-varying assumptions. These estimates were shown to be robust to parameter uncertainties, but the theoretical exponential-growth approximation (Method 1) shows wide uncertainty. Method 4 provided a maximum-likelihood effective reproduction number 2.10 (1.21, 2.95) using the first 17 epidemic days, which is consistent with estimates obtained from the other methods and an estimate of 2.36 (2.07, 2.65) for the entire autumn wave. We conclude that the reproduction number for pandemic influenza (Spanish flu) at the city level can be robustly assessed to lie in the range of 2.0–3.0, in broad agreement with previous estimates using distinct data.

Transmission dynamics and control of Ebola virus disease (EVD): a review

The complex and unprecedented Ebola epidemic ongoing in West Africa has highlighted the need to review the epidemiological characteristics of Ebola Virus Disease (EVD) as well as our current understanding of the transmission dynamics and the effect of control interventions against Ebola transmission. Here we review key epidemiological data from past Ebola outbreaks and carry out a comparative review of mathematical models of the spread and control of Ebola in the context of past outbreaks and the ongoing epidemic in West Africa. We show that mathematical modeling offers useful insights into the risk of a major epidemic of EVD and the assessment of the impact of basic public health measures on disease spread. We also discuss the critical need to collect detailed epidemiological data in real-time during the course of an ongoing epidemic, carry out further studies to estimate the effectiveness of interventions during past outbreaks and the ongoing epidemic, and develop large-scale modeling studies to study the spread and control of viral hemorrhagic fevers in the context of the highly heterogeneous economic reality of African countries.

Pros and cons of estimating the reproduction number from early epidemic growth rate of influenza A (H1N1) 2009.

BackgroundIn many parts of the world, the exponential growth rate of infections during the initial epidemic phase has been used to make statistical inferences on the reproduction number, R, a summary measure of the transmission potential for the novel influenza A (H1N1) 2009. The growth rate at the initial stage of the epidemic in Japan led to estimates for R in the range 2.0 to 2.6, capturing the intensity of the initial outbreak among school-age children in May 2009.MethodsAn updated estimate of R that takes into account the epidemic data from 29 May to 14 July is provided. An age-structured renewal process is employed to capture the age-dependent transmission dynamics, jointly estimating the reproduction number, the age-dependent susceptibility and the relative contribution of imported cases to secondary transmission. Pitfalls in estimating epidemic growth rates are identified and used for scrutinizing and re-assessing the results of our earlier estimate of R.ResultsMaximum likelihood estimates of R using the data from 29 May to 14 July ranged from 1.21 to 1.35. The next-generation matrix, based on our age-structured model, predicts that only 17.5% of the population will experience infection by the end of the first pandemic wave. Our earlier estimate of R did not fully capture the population-wide epidemic in quantifying the next-generation matrix from the estimated growth rate during the initial stage of the pandemic in Japan.ConclusionsIn order to quantify R from the growth rate of cases, it is essential that the selected model captures the underlying transmission dynamics embedded in the data. Exploring additional epidemiological information will be useful for assessing the temporal dynamics. Although the simple concept of R is more easily grasped by the general public than that of the next-generation matrix, the matrix incorporating detailed information (e.g., age-specificity) is essential for reducing the levels of uncertainty in predictions and for assisting public health policymaking. Model-based prediction and policymaking are best described by sharing fundamental notions of heterogeneous risks of infection and death with non-experts to avoid potential confusion and/or possible misuse of modelling results.

Natural History of Dengue Virus (DENV)—1 and DENV—4 Infections: Reanalysis of Classic Studies

BACKGROUND The natural history of wild-type dengue virus (DENV) infections of humans, including incubation and infectious periods, requires further study. METHODS Two experimental studies in the Philippines of DENV-4 (1924-1925) and DENV-1 (1929-1930) were reexamined. The intrinsic incubation periods were fitted to log-normal distribution using the maximum likelihood method, and the infectious and extrinsic incubation periods were assessed by proportions of successful transmissions causing clinically apparent dengue. Correlations between the intrinsic incubation period and other variables and univariate associations between clinical severity and serotype were also examined. RESULTS Mean+/-SD incubation periods were 6.0+/-1.4 and 5.7+/-1.5 days for DENV-4 and DENV-1, respectively. Significant negative correlations were observed between the incubation period and duration of fever (r=-0.43 and -0.33). Even 1 and 2 days before the onset of fever, 80.0% (95% confidence interval [CI], 44.9%-100%) and 25.0% (CI, 0%-67.4%) of biting experiments caused clinically apparent dengue. DENV-1 infections resulted in a significantly longer duration of fever than DENV-4 infections (P<.01). CONCLUSIONS Incubation period was negatively correlated with disease severity, potentially reflecting a dose-response mechanism. The historical data provided useful details concerning serotype differences in the natural history of primary DENV infections.

Early transmission dynamics of Ebola virus disease (EVD), West Africa, March to August 2014

The effective reproduction number, Rt, of Ebola virus disease was estimated using country-specific data reported from Guinea, Liberia and Sierra Leone to the World Health Organization from March to August, 2014. Rt for the three countries lies consistently above 1.0 since June 2014. Country-specific Rt for Liberia and Sierra Leone have lied between 1.0 and 2.0. Rt<2 indicate that control could be attained by preventing over half of the secondary transmissions per primary case.

Surveillance of wild birds for avian influenza virus

TOC Summary: A targeted, hypothesis-based approach and local surveys over broad geographic areas are needed.

Excess Mortality Associated With Influenza A and B Virus in Hong Kong, 1998–2009

BACKGROUND Although deaths associated with laboratory-confirmed influenza virus infections are rare, the excess mortality burden of influenza estimated from statistical models may more reliably quantify the impact of influenza in a population. METHODS We applied age-specific multiple linear regression models to all-cause and cause-specific mortality rates in Hong Kong from 1998 through 2009. The differences between estimated mortality rates in the presence or absence of recorded influenza activity were used to estimate influenza-associated excess mortality. RESULTS The annual influenza-associated all-cause excess mortality rate was 11.1 (95% confidence interval [CI], 7.2-14.6) per 100,000 person-years. We estimated an average of 751 (95% CI, 488-990) excess deaths associated with influenza annually from 1998 through 2009, with 95% of the excess deaths occurring in persons aged ≥65 years. Most of the influenza-associated excess deaths were from respiratory (53%) and cardiovascular (18%) causes. Influenza A(H3N2) epidemics were associated with more excess deaths than influenza A(H1N1) or B during the study period. CONCLUSIONS Influenza was associated with a substantial number of excess deaths each year, mainly among the elderly, in Hong Kong in the past decade. The influenza-associated excess mortality rates were generally similar in Hong Kong and the United States.

Does Glycosylation as a modifier of Original Antigenic Sin explain the case age distribution and unusual toxicity in pandemic novel H1N1 influenza

BackgroundA pandemic novel H1N1 swine-origin influenza virus has emerged. Most recently the World Health Organization has announced that in a country-dependent fashion, up to 15% of cases may require hospitalization, often including respiratory support. It is now clear that healthy children and young adults are disproportionately affected, most unusually among those with severe respiratory disease without underlying conditions. One possible explanation for this case age distribution is the doctrine of Original Antigenic Sin, i.e., novel H1N1 may be antigenically similar to H1N1 viruses that circulated at an earlier time. Persons whose first exposure to influenza viruses was to such similar viruses would be relatively immune. However, this principle is not sufficient to explain the graded susceptibility between ages 20 and 60, the reduced susceptibility in children below age 10, and the unusual toxicity observed.MethodsWe collected case data from 11 countries, about 60% of all cases reported through mid-July 2009. We compared sequence data for the hemagglutinin of novel H1N1 with sequences of H1N1 viruses from 1918 to the present. We searched for sequence differences that imply loss of antigenicity either directly through amino acid substitution or by the appearance of sites for potential glycosylation proximal to sites known to be antigenic in humans. We also considered T-cell epitopes.ResultsIn our composite, over 75% of confirmed cases of novel H1N1 occurred in persons ≤ 30 years old, with peak incidence in the age range 10-19 years. Less than 3% of cases occurred in persons over 65, with a gradation in incidence between ages 20 and 60 years.The sequence data indicates that novel H1N1 is most similar to H1N1 viruses that circulated before 1943. Novel H1N1 lacks glycosylation sites on the globular head of hemagglutinin (HA1) near antigenic regions, a pattern shared with the 1918 pandemic strain and H1N1 viruses that circulated until the early 1940s. Later H1N1 viruses progressively added new glycosylation sites likely to shield antigenic epitopes, while T-cell epitopes were relatively unchanged.ConclusionsIn this evolutionary context, Original Antigenic Sin exposure should produce an immune response increasingly mismatched to novel H1N1 in progressively younger persons. We suggest that it is this mismatch that produces both the gradation in susceptibility and the unusual toxicity. Several murine studies suggest specific cell types as a likely basis of the unusual toxicity. These studies also point to widely available pharmaceutical agents as plausible candidates for mitigating the toxic effects. The principle of Original Antigenic Sin modified by glycosylation appears to explain both the case age distribution and the unusual toxicity pattern of the novel H1N1 pandemic. In addition, it suggests pharmaceutical agents for immediate investigation for mitigation potential, and provides strategic guidance for the distribution of pandemic mitigation resources of all types.

Aerosol transmission is an important mode of influenza A virus spread

Influenza A viruses are believed to spread between humans through contact, large respiratory droplets and small particle droplet nuclei (aerosols), but the relative importance of each of these modes of transmission is unclear. Volunteer studies suggest that infections via aerosol transmission may have a higher risk of febrile illness. Here we apply a mathematical model to data from randomized controlled trials of hand hygiene and surgical face masks in Hong Kong and Bangkok households. In these particular environments, inferences on the relative importance of modes of transmission are facilitated by information on the timing of secondary infections and apparent differences in clinical presentation of secondary infections resulting from aerosol transmission. We find that aerosol transmission accounts for approximately half of all transmission events. This implies that measures to reduce transmission by contact or large droplets may not be sufficient to control influenza A virus transmission in households.

Monocyte Chemotactic Factor in Rheumatoid Arthritis Synovial Tissue PROBABLY A CROSS-LINKED DERIVATIVE OF S19 RIBOSOMAL PROTEIN

The extracts of rheumatoid arthritis-synovial lesions from seven patients possessed a strong chemotactic activity for monocytes and a negligible one for polymorphonuclear leukocytes. These results are consistent with a prominent histological feature of the synovial lesion, the mononuclear cell predominant infiltration. The major monocyte chemotactic factor in the synovial tissue extracts was purified to a single protein peak in reverse phase high performance liquid chromatography with a C4 column. NH-terminal amino acid analysis of the initial 20 residues yielded a single sequence. Surprisingly, this sequence was completely identical to that of S19 ribosomal protein. The purified sample demonstrated two protein bands in SDS-polyacrylamide gel electrophoresis with apparent molecular masses of 34 and 68 kDa. These sizes were 2 and 4 times that of S19 ribosomal protein, suggesting that the chemotactic factor would be a dimer or tetramer of S19 ribosomal protein cross-linked by factor XIIIa. A recombinant human S19 ribosomal protein was prepared as a fusion protein with a maltose binding protein in Escherichia coli. After treatment with factor XIIIa, cross-linked recombinant S19 ribosomal protein exhibited the monocyte chemotactic activity, although the untreated recombinant protein did not.