Hospitalizations for varicella in children and adolescents in a referral hospital in Hong Kong, 2004 to 2008: A time series study
© Chan et al; licensee BioMed Central Ltd. 2011
Received: 6 August 2010
Accepted: 23 May 2011
Published: 23 May 2011
Varicella accounts for significant morbidities and remains a public health issue worldwide. Climatic factors have been shown to associate with the incidence and transmission of various infectious diseases. We describe the epidemiology of varicella in paediatric patients hospitalized at a tertiary referral hospital in Hong Kong from 2004 to 2008, and to explore the possible association between the occurrence of varicella infection and various climatic factors.
The hospital discharge database of Princess Margaret Hospital was retrospectively analyzed for admissions associated with varicella from 2004 to 2008. Meteorological data were obtained from the monthly meteorological reports from the Hong Kong Observatory website. Time series analysis was performed with Poisson regression using a Generalized Estimating Equation (GEE) approach.
During the study period, 598 children were hospitalized for varicella. The mean age on admission was 57.6 months, and the mean duration of hospitalization was 3.7 days. The overall complication rate was 47%. The mean monthly relative humidity, especially in cool seasons, was inversely correlated with the monthly varicella cases of the same month.
Varicella can lead to serious complications and prolonged hospitalization, even in previously healthy children. Lower relative humidity in cool seasons is associated with higher number of paediatric varicella hospital admissions. These findings are useful for a better understanding of the pattern of paediatric varicella hospitalization in Hong Kong.
KeywordsVaricella Chickenpox Time series study Child humidity cool season
Varicella, or chickenpox, is an infectious disease caused by the ubiquitous varicella-zoster virus (VZV). It used to be considered a relatively benign communicable disease of childhood . However, serious complications can occur, including secondary bacterial skin and soft tissue infections, cerebellitis, encephalitis, pneumonia and coagulopathy [1, 2]. Hospitalization rates due to chickenpox are considerably high in developed countries, especially among children. The reported chickenpox complication rates range from 40.7% to 83.3% of children hospitalized with the condition in various studies [3–5]. Furthermore, a mortality rate of 2-3 per 100,000 affected persons has been reported . Therefore, varicella remains an important public health issue worldwide. Despite the availability of a safe and effective varicella vaccine since 1986, there is no recommendation for universal vaccination against varicella in Hong Kong [6–8].
It has been reported that the epidemiology of varicella appears to vary among different geographic regions, climatic belts, population densities, and degrees of socioeconomic development [9–11]. A number of studies have been conducted to investigate the possible association between the occurrence of varicella infection and various environmental factors, such as temperature, rainfall and humidity. Wu et al has shown that season and temperature are significantly related to varicella incidence in Taiwan . Another study performed by Kokaze et al demonstrated that the annual temperature variation affected the seasonal variations of varicella incidence in Japan . In Hong Kong, the rates and characteristics of chickenpox in children have not been clearly described.
The objective of this study was to investigate the local epidemiology and seasonal trend of varicella in hospitalized children in Hong Kong. Furthermore, the possible underlying association between varicella incidence and various climatic factors was also explored using time series methodology.
A retrospective study was conducted to analyse the clinical information of all children and adolescents aged below 18 years who were hospitalized for varicella at Princess Margaret Hospital, a tertiary referral hospital for paediatric infectious diseases, between 1 January 2004 and 31 December 2008. Cases were identified by reviewing hospital discharge records retrieved using the relevant International Classification of Diseases 9th Edition (ICD-9) principal or secondary discharge diagnostic codes for varicella or chickenpox, as in other epidemiological studies [2, 3]. The inclusion criteria were: children aged below eighteen years; principal or secondary discharge diagnostic codes being either varicella or chickenpox; and patients being admitted for varicella to Princess Margaret Hospital during the period from 1 January 2004 to 31 December 2008. All hospital admissions, not limited to emergency admissions, were included for time series analysis. Cases of herpes zoster were excluded.
Hospital discharge records for all the cases fulfilling the inclusion criteria were retrieved. Relevant demographic, clinical and admission data were extracted from the hospital discharge records. Relevant demographic data included age and gender of the patients. Relevant clinical data included the presence or absence of associated complications and the type of complications of the patients. Important admission data included the month and year of hospital admission and the length of hospitalization. The diagnosis of encephalitis was based on clinical features (decreased consciousness and seizure or focal neurologic deficit), cerebrospinal fluid pleocytosis, electroencephalographic changes and neuroimaging findings . The diagnosis of cerebellitis was based on clinical findings of ataxia and other cerebellar signs. The diagnosis of pneumonia is based on clinical, laboratory and radiological features . The diagnosis of skin and soft tissue infections was based on clinical signs (cellulitis, skin abscess, myositis or necrotizing fasciitis), and surgical scarlet fever was based on clinical features (sandpaper-like erythema, strawberry tongue and Pastia's lines).
Meteorological data including monthly mean air temperature, monthly mean relative humidity and monthly total rainfall were obtained from the overall monthly reports of meteorological observations during the study period from the website of Hong Kong Observatory .
Age of the patients (months) and length of hospitalization (days) were expressed as means, medians and ranges. Comparison among groups was performed using the chi-square or Fisher's exact test, where appropriate. Mann-Whitney test was used for non-normally distributed data.
Spearman rank correlation tests were performed to examine the relationship between monthly varicella incidence and climatic variables, including monthly mean temperature, monthly mean relative humidity and monthly rainfall, with a lag of zero to two months. The weather variables found to be correlated with varicella incidence were included in univariate and multivariate regression models. The monthly varicella incidence was modeled using a generalized estimating equations (GEE) approach, with a Poisson distribution. This model enables both specification of an over-dispersion term and a first-order autoregressive structure that accounts for the autocorrelation of monthly numbers of varicella cases. Besides the meteorological factors, the autoregression term ylag (number of varicella cases at one month lag) was also included in the model. As GEE are not a full likelihood-modeling method, the quasi-likelihood based information criterion (QICu) developed by Pan was computed to select the most parsimonious model . When using QICu to compare two models, the model with the smaller statistic was preferred. Model with fewest parameters was preferred.
A p-value < 0.05 will be considered statistically significant. Data were analysed using SPSS software (version 16.0; SPSS Inc, Chicago, IL, USA) and SAS software (version 9 for Windows; SAS Institute, Inc., Cary, North Carolina).
Hospitalization figures for current study
There were a total of 328 male (54.8%) and 270 female (45.2%) patients hospitalized for varicella, with a slight male predominance at a ratio (1.2 : 1). The mean age was 57.6 months (range 1-204 months). 64 patients (10.7%) were aged 6 months or under, 335 (56.0%) were aged 7 months to 5 years, 126 (21.1%) were aged above 5 years to 10 years, and 73 (12.2%) were aged above 10 years.
Length of hospitalizatio
The overall mean length of hospitalization was 3.7 days (range 1-27 days). It was noted that the length of hospitalization was significantly longer in patients with complicated varicella compared to those without complications, with a mean of 4.3 days (range 1-24 days) vs 3.1 days (range 1-27 days), respectively (p < 0.001).
Complications of varicella
Among all cases, 281 out of 598 (47%) patients developed complications of varicella. It was found that skin and soft tissue infections were the most commonly occurred complication, accounting 121 out of 281 patients (43.1%) who developed complicated varicella. Ninety-nine patients (35.2%) had surgical scarlet fever. Neurological complications were observed in 51 patients (18.1%), including one case of encephalitis, one case of aseptic meningitis, four cases of cerebellitis, five cases of afebrile seizure and 40 cases of febrile seizure. Pneumonia was identified in 23 patients (8.2%). Nine patients (3.2%) developed other complications, including clinical sepsis, immune thrombocytopenic purpura and Henoch-Schonlein purpura in 4, 3 and 2 patients, respectively. There were five patients with severe complicated varicella requiring intensive care but there was no fatal case. The mean age of patients with complicated varicella was 51.6 months (range 1-204 months), which was not significantly different from those without complications (mean age of 62.9 months, range 1-204 months, p = 0.183).
Periodicity and climatic factors variation
Coefficients of Spearman correlation analysis between varicella cases (2004-2008) at Princess Margaret Hospital and climatic variables (* P < 0.05)
Result for best fitting model with GEE analysis
Where Y, ylag and Humid0 stand for number of varicella cases in that month, number of varicella cases in the previous one month (at lag of one) and monthly mean relative humidity in the same month. The monthly mean relative humidity of the same month was found to be a predictor for the monthly varicella cases.
We also conducted a stratified analysis by two seasons (October-March defined as cool seasons and the other months hot season) to see how the effect of humidity varies. Humidity remained a significant predictor of varicella in the cool season model (Beta = -3.02, 95% confidence interval -5.47, -0.58, p = 0.0152), while it was not significant in the hot season model.
This is the first time series study exploring the possible assoication between paediatric varicella hospitalization and weather factors in Hong Kong. Our dataset is relevant because the involved hospital, Princess Margaret Hospital, is a tertiary referral hospital for paediatric infectious diseases. The involved hospital received not only patients in its locality, but also patients transferred from other public hospitals. Age-specific hospitalization rates vary between studies and generally peak in children aged 1-4 years as they did in the present study [3, 4, 18]. The complication rates reported in literature among children hospitalized for varicella infection range from 40.7% to 83.3% [3–5]. In the present study, the complication rate was 47%. The distribution of different complications in this study was similar to that reported in other paediatric series, with skin and neurological complications occurring most commonly [4, 5, 14, 18].
Time series analysis has been used extensively to study the effect of weather factors on various infectious diseases [19–23]. In the current study, Poisson regression with generalized estimating equations (GEE) was used to examine the potential effect of various weather factors on the incidence of chickenpox in paediatric patients hospitalized in a tertiary hospital in Hong Kong during the period of 2004 to 2008. Monthly mean relative humidity, especially in cool seasons, was found to be inversely correlated with the number of varicella cases in the same month.
It has been suggested that the transmission potential of the VZV virus might be adversely affected by a combination of high ambient temperatures and humidity in tropical regions . A review of data from Southeast Asia showed that the peak incidence of varicella infection occurred during cooler months and in cooler, more temperate regions . Outbreaks of varicella appear to be more common in the cooler months of the year in India and in Thailand; incidence rates are highest in the temperate northern region of the country [26–28]. Data from studies outside Southeast Asia also support the notion of reduced VZV transmission in hot, humid climates . For example, Maretic and Cooray found that the number of chickenpox admissions to a regional hospital in the cool, dry season more than doubled compared to the hot monsoon period in their study of a Sri Lankan population . Lolekha et al reported that the age-adjusted seroprevalence was significantly higher in the cooler than in the warmer regions in Thailand . Wu et al has shown that season and temperature are significantly related to varicella incidence in Taiwan . Another study performed by Kokaze et al demonstrated that the annual temperature variation affected the seasonal variations of varicella incidence in Japan .
Varicella patients excreted the virus from their respiratory tract or vesicles and disseminated the virus to the environment via an aerosol route . VZV can be transmitted via airborne spread . However, the exact mechanism for the potential association between relative humidity and varicella incidence and transmission is not known. One could hypothesize that with a lower relative humidity, the density of air particles would be lower and hence the VZV could have stayed longer in the air during its airborne transmission to a long distance. Another possibility could be that in a dry environment with more dry skin, the patients with chickenpox could suffer from excessive itchiness where more skin scratching would facilitate more viral spread to the environment. In addition, the more common respiratory symptoms such as cough and sputum production in the dry season could enhance the aerosol spread of the VZV. However, all these hypotheses warrant further studies to explore the mechanism underlying the association between relative humidity and varicella infection.
Our study may have some limitations. First, it is based on hospitalization discharge records, and therefore children with mild disease and complications may have been excluded due to the absence of hospitalization. The potential climatic association with varicella might have been dampened by this selective sample. Second, using discharge diagnoses codes for data collection might result in errors due to the possibility of miscoding. But studies investigating the accuracy of ICD-9 codes for various diseases in different countries suggested that this is a useful tool for epidemiological studies [4, 34, 35]. In addition, chickenpox is considered to be a easily recognisable disease [2, 36]. Third, as the weather fluctuation in Hong Kong is relatively small, the potential effect of weather factors might have been better captured by decreasing the data interval, for instance, a weekly or daily basis, or increasing the study period. It should be acknowledged that varicella transmission is multifactorial [37–39]. Apart from potential weather attributes, other environmental and host factors may also affect the incidence and transmission of the disease. Further studies to incorporate more variables such as age, geographic location, and frequency of contact expressed by population density and school contact patterns are warranted in the future.
In conclusion, varicella is a disease that can provoke serious complications and long hospital stays, even in healthy children. Lower relative humidity, especially in the cool seasons, is associated with higher number of paediatric varicella hospital admissions. These findings are useful for a better understanding of the pattern of paediatric varicella hospitalization in Hong Kong.
Conflict of interest statement
The authors declare that they have no competing interests.
We would like to thank Dr Christine Wong, Dr SK Chuang and Dr Thomas Tsang of the Centre for Health Protection, Department of Heath for their kind support and assistance with retrieval of notification figures for chickenpox. We would also like to thank Prof SS Lee of the School of Public Health and Primary Care, the Chinese University of Hong Kong for helpful discussions.
- Boelle PY, Hanslik T: Varicella in non-immune persons: incidence, hospitalization and mortality rates. Epidemiol Infect. 2002, 129: 599-606. 10.1017/S0950268802007720.PubMed CentralView ArticlePubMed
- Marchetto S, de Benedictis FM, de Martino M, Versace A, Chiappini E, Bertaine C, Osimani P, Cordiali R, Gabiano C, Galli L: Epidemiology of hospital admissions for chickenpox in children: an Italian multicentre study in the pre-vaccine era. Acta Paediatrica. 2007, 96: 1490-3. 10.1111/j.1651-2227.2007.00465.x.View ArticlePubMed
- Galil K, Brown C, Lin F, Seward J: Hospitalizations for varicella in the United States, 1988 to 1999. Pediatr Infect Dis J. 2002, 21: 931-5. 10.1097/00006454-200210000-00009.View ArticlePubMed
- Lin F, Hadler JL: Epidemiology of primary varicella and herpes zoster hospitalizations: the pre-varicella vaccine era. J Infect Dis. 2000, 181: 1897-905. 10.1086/315492.View ArticlePubMed
- Carapetis JR, Russell DM, Curtis N: The burden and cost of hospitalised varicella and zoster in Australian children. Vaccine. 2004, 23: 755-61. 10.1016/j.vaccine.2004.07.025.View ArticlePubMed
- Tseng HF, Tan HF, Chang CK, Wang LY, Yang SE, Liau MY, Pan MJ, Chen CY, Huang SP: A seroepidemiology study of varicella among children aged 0-12 years in Taiwan. Southeast Asian Journal of Tropical Medicine and Public Health. 2005, 36: 1201-7.PubMed
- Heininger U, Seward JF: Varicella. Lancet. 2006, 368: 1365-76. 10.1016/S0140-6736(06)69561-5.View ArticlePubMed
- Nguyen HQ, Jumaan AO, Seward JF: Decline in mortality due to varicella after implementation of varicella vaccination in the United States. N Engl J Med. 2005, 352: 450-8. 10.1056/NEJMoa042271.View ArticlePubMed
- Fairley CK, Miller E: Varicella-zoster virus epidemiology: a changing scene?. J Infect Dis. 1996, 174 (Suppl 3): S314-9.View ArticlePubMed
- Yawn B, Yawn R, Lydick E: Community impact of childhood varicella infection. J Pediatr. 1997, 130: 759-65. 10.1016/S0022-3476(97)80019-4.View ArticlePubMed
- Gershon AA, Takahashi M, White CJ: Varicella vaccine. Vaccines. Edited by: Plotkin SA, Orenstein WA. 1999, Philadephia: Saunders, 475-507.
- Wu PY, Li YC, Wu HD: Risk factors for chickenpox incidence in Taiwan from a large-scale computerized database. International J of Dermat. 2007, 46: 362-366. 10.1111/j.1365-4632.2006.03050.x.View Article
- Kokaze A, Yoshida M, Sekine Y, Ishikawa M, Kurokochi T, Uchida Y, Matsunaga N, Takashima Y: The magnitude of variation in temperature within a year has an effect on the seasonal variations of chickenpox incidence in Japan. Epidemiol Infect. 2001, 126: 269-277.PubMed CentralPubMed
- Gnann JW: Varicella-zoster virus: atypical presentations and unusual complications. J Infect Dis. 2002, 186 (Suppl 1): S91-8.View ArticlePubMed
- Kong Hong Observatory. [http://www.hko.gov.hk]
- Pan W: Akaike's information criterion in generalized estimating equations. Biometrics. 2001, 57: 120-125. 10.1111/j.0006-341X.2001.00120.x.View ArticlePubMed
- Centre for Health Protection, Department of health, Hong Kong. [http://www.chp.gov.hk]
- Ziebold C, von Kries R, Lang R, Weigl J, Schmitt HJ: Severe complications of varicella in previously healthy children in Germany: a 1-year survey. Pediatrics. 2001, 108: e79-10.1542/peds.108.1.79.View ArticlePubMed
- Lu L, Lin H, Tian L, Yang W, Sun J, Liu Q: Time series analysis of dengue fever and weather in Guangzhou, China. BMC Public Health. 2009, 9: 395-10.1186/1471-2458-9-395.PubMed CentralView ArticlePubMed
- Luz PM, Mendes BVM, Codeco CT, Struchiner CJ, Galvani AP: Time Series Analysis of Dengue Incidence in Rio de Janeiro, Brazil. Am J Trop Med Hyg. 2008, 79: 933-939.PubMed
- Tian LW, Bi Y, Ho SC, Liu WJ, Liang S, Goggins WB, Chan EYY, Zhou S, Sung JJY: One-year delayed effect of fog on malaria transmission: a time-series analysis in the rain forest area of Mengla County, south-west China. Malar J. 2008, 7: 110-10.1186/1475-2875-7-110.PubMed CentralView ArticlePubMed
- Tong S, Hu W: Climate variation and incidence of Ross river virus in Cairns, Australia: a time-series analysis. Environ Health Perspect. 2001, 109 (12): 1271-1273. 10.1289/ehp.011091271.PubMed CentralView ArticlePubMed
- Hu W, Tong S, Mengersen K, Connell D: Weather variability and the incidence of cryptosporidiosis: comparison of time series poisson regression and SARIMA models. Ann Epidemiol. 2007, 17: 679-688. 10.1016/j.annepidem.2007.03.020.View ArticlePubMed
- Garnett GP, Cox MJ, Bundy DAP, Didier JM, Catherine St J: The age of infection with varicella-zoster virus in St Lucia, West Indies. Epidemiol Infect. 1993, 110: 361-372. 10.1017/S0950268800068308.PubMed CentralView ArticlePubMed
- Lee BW: Review of varicella zoster seroepidemiology in India and Southeast Asia. Trop Med Int Health. 1998, 3: 886-890. 10.1046/j.1365-3156.1998.00316.x.View ArticlePubMed
- White E: Chickenpox in Kerala. Indian Journal of Public Health. 1978, 22: 141-151.PubMed
- Basu RN, Jezek Z, Ward NA: The eradication of smallpox from India. 1979, World Health Organization, New Delhi, 255-270.
- Venkitaraman AR, John TJ: The epidemiology of varicella in staff and students of a hospital in the tropics. International Journal of Epidemiology. 1984, 13: 502-505. 10.1093/ije/13.4.502.View ArticlePubMed
- Lee BW: Review of varicella zoster seroepidemiology in India and South-east Asia. Tropical Medicine and Internation Health. 1998, 3 (11): 886-890. 10.1046/j.1365-3156.1998.00316.x.View Article
- Maretic Z, Cooray MPM: Comparisons between chickenpox in a tropical and a European country. Journal of Tropical Medicine and Hygiene. 1963, 66: 311-315.PubMed
- Lolekha S, Tanthiphabha W, Sornchai P, Kosuwan P, Sutra S, Warachit B, Chupupprakarn S, Hutagalung Y, Weil J, Bock HL: Effect of climatic factors and population density on varicella. Am J Trop Med Hyg. 2001, 64 (3, 4): 131-136.PubMed
- Suzuki K, Yoshikawa T, Ihira M, Ohashi M, Suga S, Asano Y: Spread of varicellazoster virus DNA to the environment from varicella patients who were treated with oral acyclovir. Pediatr Int. 2003, 45 (4): 458-60. 10.1046/j.1442-200X.2003.01746.x.View ArticlePubMed
- Asano Y, Yoshikawa T, Ihira M, Furukawa H, Suzuki K, Suga S: Spread of varicellazoster virus DNA to family members and environments from siblings with varicella in a household. Pediatrics. 1999, 103 (5): e61-10.1542/peds.103.5.e61.View ArticlePubMed
- Shinefield H, Black S, Williams WR, Marchant C, Reisinger K, Stewart T, Meissner HC, Guerrero J, Klopfer SO, Xu J, Schodel F, Kuter BJ: Dose Selection Study Group for Proquad. Dose-response study of a quadrivalent measles, mumps, rubella and varicella vaccine in healthy children. Pediatr Infect Dis J. 2005, 24: 670-5. 10.1097/01.inf.0000172901.29621.e9.View ArticlePubMed
- Heininger U, Desgrandchamps D, Schaad UB: Seroprevalence of varicella-zoster virus IgG antibodies in Swiss children during the first 16 months of age. Vaccine. 2006, 24: 3258-60. 10.1016/j.vaccine.2006.01.026.View ArticlePubMed
- Wise RP, Salive ME, Braun MM, Mootrey GT, Seward JF, Rider LG, Krause PR: Postlicensure safety surveillance for varicella vaccine. JAMA. 2000, 284: 1271-9. 10.1001/jama.284.10.1271.View ArticlePubMed
- Nelson WE, Behrman RE, Kliegman RM: Streptococcal infections. Textbook of Pediatrics. Edited by: James Todd. 1996, Philadelphia, PA: W.B. Saunders, 750-751. 15
- Brisson M, Edmunds WJ, Law B, Gay NJ, Walld R, Brownell M, Roos L, De Serres G: Epidemiology of varicella zoster virus infection in Canada and the United Kingdom. Epidemiol Infect. 2001, 127: 305-314.PubMed CentralPubMed
- O'Grady KA, Merianos A, Patel M, Gilbert L: High seroprevalence of antibodies to varicella zoster virus in adult women in a tropical climate. Tropical Medicine and International health. 2000, 5 (10): 732-736. 10.1046/j.1365-3156.2000.00622.x.View ArticlePubMed
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/11/366/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.