Journal of Infectious Diseases and Treatment Open Access

  • ISSN: 2472-1093
  • Journal h-index: 7
  • Journal CiteScore: 1.06
  • Journal Impact Factor: 0.77
  • Average acceptance to publication time (5-7 days)
  • Average article processing time (30-45 days) Less than 5 volumes 30 days
    8 - 9 volumes 40 days
    10 and more volumes 45 days

Original Research Article - (2022) Volume 8, Issue 8

Meteorological Factors Facilitated Increased Legionnaires' Disease Notifications and Mortality Pre-COVID-19 and COVID-19 in Hong Kong
Yin Zhang1, Frank Y. Chen1*, Geoffrey K.F. Tso1 and Lin Yang2
 
1Department of Management Sciences, City University of Hong Kong, Hong Kong
2Department of Health and Social Sciences, The Hong Kong Polytechnic University, Hong Kong
 
*Correspondence: Frank Y. Chen, Department of Management Sciences, City University of Hong Kong, Hong Kong, Email:

Received: 27-Sep-2022, Manuscript No. IPJIDT-22-14441; Editor assigned: 29-Sep-2022, Pre QC No. IPJIDT-22-14441 (PQ); Reviewed: 13-Oct-2022, QC No. IPJIDT-22-14441; Revised: 18-Oct-2022, Manuscript No. IPJIDT-22-14441 (R); Published: 25-Oct-2022, DOI: 10.36648/2394-9988-8.8.40

Abstract

Annual notification number of Legionnaires’ disease (LD) increased from 66 cases in 2015 to 104 cases in 2020 in Hong Kong, with a case fatality rate (CFR) of 17.3% in 2020. This study investigated possible meteorological factors associated with confirmed LD cases and deaths. The negative binomial regression model was used for monthly LD confirmed cases and the Poisson regression model was used for monthly LD mortality cases. A decision tree model was adopted to identify the thresholds of meteorological factors if any. The confirmed LD cases were highest in July while the LD mortality was highest in August. A negative binomial regression model confirmed that maximum air temperature (P=0.081, RR=1.3, 95% CI=0.97-1.75) and rainfall (P=0.009, RR=1.02, 95% CI=1.0004-1.003) were positively correlated with the increase in confirmed LD cases. A Poisson regression model confirmed that maximum air temperature (P=0.034, RR=1.71, 95% CI=1.04-2.80) was positively correlated with the surge in LD mortality. When rainfall exceeds 78.8 mm, the public should be alerted of the increased risk of contracting LD. Temperature and rainfall were positively associated with the confirmed LD cases in Hong Kong.

Keywords

Legionnaires’ disease; Meteorological causes; Negative binomial regression; Poisson regression; Decision tree

INTRODUCTION

Legionnaires’ disease (LD) ranked fourth (17.3%) in terms of Case Fatality Rate (CFR) among the 13 notifiable diseases with all other causes (AOC) in Hong Kong in 2020 (Supplementary Table S1). The associations of environmental factors with confirmed Legionnaires’ disease (LD) cases have been reported in the literature; Tanimoto, et al. found apparent evidence of Legionella spp. bacteria being transmitted via water mist in hot spring spas among the ageing populations in Japan [1]. Lam, et al. found that Legionella spp. bacteria widely colonized in cooling towers in Singapore, which has abundant high rise buildings [2]. With the COVID-19 lockdown and subsequent reopening, Rhoads, et al. expressed concern that extended water stagnation in buildings would increase the number of confirmed LD cases, and Palazzolo, et al. reported that an Italian male aged 40 years who worked as a restaurant dishwasher contracted LD 3 days after the COVID-19 lockdown restriction was eased in Rome [3,4]. O’Connor, et al. developed a PCR test kit to detect bacteria in the Legionella genus in water with a single test. Hozalski, et al. suggested that flushing of stagnant water after the COVID-19 shutdown could reduce the Legionella infection risk, and Lancaster, et al. encouraged building management personnel to monitor water stagnancy before reopening to reduce morbidity and mortality [5-7].

The literature reports clinical epidemiological evidence of LD and COVID-19 pneumonia coinfection, with Legionella spp. being the main species isolated in Europe, the United States (US), and Canada [8-12]. Dey, et al. reported LD and COVID-19 coinfection in patients, similar to that of “1918 Spanish flu” bacterial coinfection. Lai, et al. warned that up to 50% of patients with LD and COVID-19 coinfection do not survive [13,14].

Other works in the literature identified climatic variables and seasons associated with the number of confirmed LD cases. Ricketts, et al. explored quarterly and annual data and found some evidence suggesting that the effect of the interaction of relative humidity and temperature on the increased number of confirmed LD cases in England was highly significant. Precipitation increased the confirmed LD cases in Taiwan. Brandsema, et al. found an association of warm and wet weather with an increased number of confirmed LD cases in the summer of 2010 in the Netherlands [15-17]. Oda, et al. reported a male aged 62 years contracting LD after the Heavy Rain Event of July 2018 in Japan [18]. Temperature, rainfall, and atmospheric pressure have been associated with the numbers of confirmed LD cases in Denmark, Germany, Italy and the Netherlands [19]. The average temperature was significantly associated with the number of confirmed LD cases in South Korea [20]. Simmering, et al. analyzed LD cases in the middle Atlantic region and confirmed that LD cases were more common in summer in the US [21]. Falconi, et al. observed a seasonal pattern from mid-September to mid-August in the number of confirmed LD cases after 2003 in the US [22]. Ozeki, et al. found two seasonal peaks of confirmed LD cases between June and November 2005-2009 in both Tokyo and Saitama, Japan [23].

Previous literature focused on the association of the number of confirmed LD cases with environmental factors, epidemiological factors, or meteorological factors in Asia, Europe, or North America. Accordingly, our research aims to fill the knowledge gap and analyze the number of confirmed LD cases and deaths in association with potential meteorological factors in HK.

Methods

Data

We obtained

• Daily death micro-data sets from the Census and Statistical Department of HK from 2015-2020,

• Monthly meteorological variables from the HK Observatory Department from 2015-2021 and

• Monthly number of confirmed LD cases from the Department of Health from January 2015 to November 2021 [24-25]. We retrieved LD deaths by using the International Classification of Disease, 10th revision (ICD-10) A48.1 in the death registry data from the Census and Statistics Department of Hong Kong. All methods were performed in accordance with the relevant guidelines and regulations.

Statistical Analysis

We employed negative binomial regression to model our raw monthly confirmed LD case as the dependent variable considering raw meteorological factors as independent variables since we found that the raw monthly confirmed LD cases are over dispersed counts variables where the mean of the confirmed LD cases does not equal the variance confirmed LD cases. While, we used Poisson regression to model raw monthly LD deaths with raw meteorological variables, because the mean and variance of the LD deaths are almost the same. Finally, we constructed a decision tree to identify possible weather thresholds for public policy decision making regarding LD.

Results

Confirmed Cases Mortality

LD ranked the fourth (17.3%) in terms of CFR among the 13 notifiable diseases with all other cause (AOC) mortality data in HK in 2020. The CFR of LD more than doubled (101.9%) in HK from 2019-2020 (Figure 1), when the COVID-19 outbreak began.

infectious-fatality

Figure 1: Annual change in the case fatality rate of the 6 leading causes of death from 2018-2020 in Hong Kong.

A total of 591 confirmed LD cases were reported in HK from January 2015 to November 2021. We found that the confirmed LD cases increased in HK during this time (Figure 2).

infectious-ld

Figure 2: The bar chart shows the accumulated monthly number of confirmed LD cases. We can see that July had the largest accumulated number of confirmed LD cases in HK, 2015-2020.

June to August had the largest number of confirmed LD cases (40%), and the summer month of July accounted for the largest number of accumulated confirmed LD cases (15%).

The annual rate of confirmed LD cases increased from 0.88 to 1.4 per 100,000 standard populations from 2015 to 2018 and remains 1.4 PTSP in 2018, 2019 and 2020, respectively. An abnormal surge in the number of monthly confirmed LD cases (21 cases) in February 2020 is shown in Supplementary Figure S1. This was the month immediately after the first case of COVID-19 case was reported in HK in January 2020.

53 LD deaths (case fatality rate 10.1%) were identified in HK from 2015 to 2020. We found that the annual increase in LD CFR doubled from 2019 to 2020 in HK after the COVID-19 outbreak. Leung, et al. found a more than 4-fold increase in the number of confirmed LD cases from 2005 to 2015 [26]. August had the largest number of accumulated LD deaths in our study period (Figure 3).

infectious-death

Figure 3: The bar chart shows the accumulated monthly number of LD-related death cases. We can see that August had the largest number of accumulated LD-related deaths in HK, 2015-2020.

Descriptive Analysis of Confirmed LD Cases and Mortality with Monthly Weather Variables

We now examined the associations of confirmed LD cases and mortality in HK with the city’s meteorological variables, described below. Supplementary Table S2 shows the monthly number of confirmed LD cases with meteorological variables covering 83 months in HK from 2015-2021. Supplementary Table S3 shows monthly LD mortality with meteorological variables covering 72 months in HK from 2015-2020. Since the LD mortality data we retrieved were daily data, we stratified the LD-related mortality into months per year. We aimed to determine whether the number of confirmed LD confirmed cases and related deaths in HK were associated with monthly meteorological variables, including maximum air temperature (in °C), average air temperature (in °C), relative humidity and rainfall (in mm)

Negative Binomial Regression for LD Confirmed Cases

Previous studies used Poisson regression to model confirmed LD cases. On the basis of their approach, we conducted the goodness of fit and Wald’s statistics tests, which rejected the null hypothesis (Supplementary Appendix D). Then, we constructed a negative binomial regression model, and the results were confirmed by the t test, log-likelihood ratio test, and χ2 tests (Supplementary Appendix E).

Using the number of confirmed LD cases in the 83 month dataset from January 2015 to November 2021, we constructed a negative binomial regression model to study the marginal effects of air temperature, relative humidity, and rainfall on confirmed numbers. Table 1 shows that the maximum air temperature and rainfall amount were statistically associated, but the relative humidity and average air temperature were not [22,27,28].

Variable Coefficient S.E. Pr(>z) RR 95% CI
Intercept 0.478 1.625 0.769 1.61 0.07-39
Maximum air temperature 0.264 0.151 0.081* 1.3 0.97-1.75
Average air temperature -0.231 0.154 0.135 0.79 0.60-1.075
Relative humidity 0.493 1.871 0.792 0.61 0.02-23
Rainfall 0.002 0.001 0.009*** 1.02 1.0004-1.003

Table 1: Negative binomial model results for the number of confirmed LD cases associated with monthly weather variables. CI, confidence interval; RR, relative risk; S.E., Standard Error; 95% CI, 95% Confidence Interval of relative risk. *p<0.1, **p<0.05, ***p<0.01. N=591 number of confirmed LD cases.

The relative risk associated with the maximum air was 1.3 while the relative risk associated with the rainfall amount was 1.02. This confirmed that monthly maximum air temperature was the most important factor positively correlated to the number of monthly confirmed LD cases, followed by rainfall. To examine the effect of the month of July, we defined the month of July of every year as a categorical variable (true/false, true if the confirmed case occurred in July, false if otherwise). The July effect was statistically significant (P=0.025), as shown in Supplementary This quantitatively confirmed our visual finding that July had the largest number of accumulated monthly confirmed LD cases, as shown in Figure 3. Furthermore, to examine the February 2020 spike in confirmed LD cases, we defined February 2020 as a categorical variable (true/false, true if the confirmed case occurred in February 2020, false if otherwise). Similarly, the February 2020 spike was statistically significant (P=0.016) (Table S4-S7).

The relative risk associated with the February 2020 spike was 3.27 (95% CI 1.24-8.58). This further confirmed the abnormally large number of confirmed LD cases in February 2020 amid the COVID-19 pandemic, as shown in Supplementary Figure S1.

Poisson Regression for LD Mortality

The 72 month monthly LD-related mortality data had a Poisson distribution because the mean (0.74) almost equaled the variance (0.79). We also confirmed the hypothesis by the goodness of fit test (Supplementary Appendix F). Poisson regression (Supplementary Appendix G) was applied to study the marginal effects of air temperatures, relative humidity, and rainfall on mortality numbers. Table 2 shows that the factors of maximum air temperature, average air temperature, and rainfall amount were statistically significant, but relative humidity was not. The relative risk associated with the maximum air temperature equaled 1.71. This confirmed that monthly maximum air temperature was the most important factor positively related to the number of monthly LD deaths (Table S8 -S9).

Variable Coefficient S.E. Pr(>z) RR 95% CI
Intercept -6.585 2.622 0.012** 0.001 8×10-6-0.24
Maximum air temperature 0.536 0.252 0.034* 1.71 1.04-2.80
Average air temperature -0.424 0.253 0.093* 0.65 0.40-1.07
Relative humidity 3.14 2.397 0.190 23.11 0.21-2537
Rainfall -0.002 0.001 0.035** 0.997 0.9957-0.9998

Table 2: Poisson model results for the number of LD-related deaths associated with monthly weather variables. CI, confidence interval; RR, relative risk; S.E.,Standard Error; 95% CI, 95% Confidence Interval of relative risk. *p<0.1, **p<0.05, ***p<0.01. N=53 number of LD-related deaths.

Decision Management for LD

Considering that maximum monthly air temperature and rainfall amount were important factors that were positively associated with the number of monthly confirmed LD cases, we built a decision tree model [29] (Supplementary Appendix G). When rainfall reaches the “danger” zone, i.e., when the amount of rainfall is above 78.8 (in mm) and the Gini impurity of rainfall equals 3.4%, public messaging should alert the general public and healthcare providers to the risk of new-onset LD (Supplementary Figure S2). We found that a large number of confirmed LD cases occurred when rainfall was higher than or equal to 334.8 (in mm) on the right side of the tree in Supplementary Figure S3.

Discussion

Our retrospective study analyzed a total of 591 confirmed LD cases in HK from 2015-2021 and 53 LD-related deaths in HK from 2015- 2020. The LD case fatality rate remained high, at over 10%. An abnormal increase in the number of confirmed LD cases occurred in February 2020, which was the month immediately following that of the first reported COVID-19 case in HK (January 2020). An increase in the LD notification number was observed from 2015 to 2020 (157%), and the annual increase of CFR of LD doubled (101.9%), from 2019 to 2020, after the COVID-19 outbreak. Overall, June to August had the largest number of confirmed LD cases (40%), which is in accordance with that of Ricketts, et al., 15 and the summer month of July accounted for the largest number of accumulated confirmed LD cases (15%) from 1 January 2015 to 30 November 2021. In addition to HK, North America, Europe, and Australia also experienced increases in the number of confirmed LD cases amid the COVID-19 pandemic [30]. We also found that the number of confirmed LD cases and related deaths had a lag of approximately one month, from July to August, which is in accordance with Riddell, et al. who stated that the incubation period of Legionella spp. is 2-10 days. We also found that the LD-related death characteristics by gender and age are highly consistent with that of the confirmed LD cases in Leung, et al. [26].

As reported in the previous study [17], warm weather is an important factor in LD development [22]. Regarding the number of monthly confirmed LD cases in HK from January 2015-November 2021, our study data quantitatively supplements that of Leung, et al. in terms of the number of confirmed LD cases in HK (Supplementary Table S4). Furthermore, we are the first to report that the February 2020 spike was positively correlated with an abnormally large number of confirmed LD cases in HK in 2020; February 2020 was the month immediately following that of the first reported case of COVID-19 in HK (January 2020). In recent years, there has been increasing interest in studying the meteorological factors potentially associated with the increase in the number of confirmed LD cases. Unexpectedly, our study found that in HK, relative humidity (P>0.05) was not significantly associated with the number of confirmed LD cases or related deaths; these results are in contrast with those of Ricketts, et al. who found a significant association with humidity in England and Wales, but in accordance with those of Chen, et al. in Taiwan [15,16]. Temperature and amount of rainfall, on the other hand, are related to the confirmed LD cases in HK, which is in accordance with studies in Scotland, Netherlands, Taiwan, Denmark, and Seoul Korea [17,19,20,26,27].

Strong seasonal patterns were observed in the monthly data for maximum air temperature, average air temperature, relative humidity, and rainfall, though the two air temperature variables were highly correlated visually, which is in accordance with Ozeki, et al. Our initial simple analysis using 83 month average (83 MA) and relative risk found that the probability of contracting LD was inconclusive; for example, the relative risks associated with the monthly maximum air temperature above the 83 MA (26.9°C, RR=0.98, 95% CI=0.88-1.07), the monthly average air temperature above the 83 MA (24.2°C, RR=1.03, 95% CI=0.93-1.13), relative humidity and rainfall were inconclusive. Likewise, we found that the probabilities of dying from LD were increased by 72% and 62% when exposed to a monthly maximum air temperature above the 72 MA (26.7°C, RR=1.72, 95% CI=1.06-2.79) and a monthly average air temperature above the 72 MA (24.05°C, RR=1.62, 95% CI=1.014-2.595); however, the results of relative humidity and rainfall were inconclusive (RR=0.65, 95% CI=0.416- 1.018; RR=1.03, 95% CI=0.658-1.620), respectively (Supplementary Appendix H). Therefore, we further used more appropriate statistical models to analyze the raw number of confirmed LD cases (instead of transformed variable) as the response variable and raw weather variables (instead of transformed variables) as the explanatory variables to achieve conclusive results. Our study is the first to use negative binomial regression to model the number of confirmed LD cases and confirm that maximum air temperature is a robust risk factor positively related to the number of monthly confirmed LD cases (P=0.081, RR=1.3, 95% CI=0.71- 1.75). Rainfall was also positively related to the number of monthly confirmed LD cases with a relative risk marginally higher than one (P=0.009, RR=1.02, 95% CI=1.0004-1.003). Furthermore, our Poisson regression model confirmed that maximum air temperature was positively related to the number of LD-related deaths (P=0.034, RR=1.71, 95% CI=1.04-2.80) [31,32].

From a public healthcare management policy perspective, our retrospective study suggests that the local government should inform the general public more about LD risk in HK. Frequent disinfection of water pipes in buildings during COVID-19 lockdowns should be implemented before reopening to prevent LD. Testing for LD and COVID-19 coinfection are also suggested to avoid missing a potential LD case in a COVID-19 patient. These are ongoing public health concerns. Rainfall exceeding 78.8 (in mm) is likely to increase the likelihood of contracting LD in HK, which is in accordance with Oda, et al. [18]. We hypothesized that during the hot summer months (June-August) and rainy days in Hong Kong, more people tend to stay indoors longer with the air conditioning turned on in the office or at home in a typical summer month [29]. This could potentially result in people inhaling contaminated water aerosols released from the air conditioning [33]. As such, this could potentially increase the risk of contracting and dying from LD during the hot summer months. However, this hypothesis should be tested by future study. Weather related risk alerts should be disseminated from the HK Observatory Department to the general public and healthcare professionals when the weather risk threshold identified in our study is reached to prevent LD in HK; government should also educate people on regular air conditioning maintenance during hot summers to minimize confirmed LD cases and deaths in HK.

Conclusion

The annual case fatality rate of LD doubled from 2019 to 2020 in HK after the COVID-19 outbreak. The risk of contracting LD in July (P=0.025, RR=1.81) was 81% higher than the risk in any other month. Negative binomial regression confirmed that maximum air temperature (P=0.081, RR=1.3) and rainfall (P=0.009, RR=1.02) were positively related to the number of confirmed LD cases. Poisson regression found that maximum air temperature (P=0.034, RR=1.71) was positively correlated to the number of LD-related deaths. From a public policy perspective, the general public and healthcare professionals should be alerted when the amount of rainfall increases above 78.8 (in mm), which will increase the likelihood of contracting LD in HK because people would be highly likely or likely to stay indoors longer during rainy days and hot summers with air conditioning on for potential inhalation of contaminated air aerosols.

Author Contributions Statement

Y.Z. conceived and designed the analysis; contributed analysis tools; conducted the analysis; drafted the paper and revised the paper. F.C. conceived and designed the analysis; reviewed and revised the paper; and managed the project. G.T. conceived and designed the analysis; collected data; reviewed and revised the paper. L.Y. reviewed and revised the paper; checked compliance. All authors reviewed the manuscript.

Declaration of Competing Interest

The authors declare that there are no competing interests.

Data Availability Declaration

Some examples of data availability statements are given below: The data that support the findings of this study are available:

• From [the Census and Statistical Department, Hong Kong, SAR, China]: Daily mortality per person micro data sets from the Census and Statistical department in HK 2015-2020,. Restrictions apply to the availability of these data, which were used under license for this study. Such data can be obtained via a direct request to the Census Department (with a certain fee).

• From the Department of Health: Monthly legionaries’ disease numbers from January 2015 to November 2021 are available. The number of notifiable infectious diseases by month is available from the Department of Health, The Government of the Hong Kong Special Administrative Region. January 2015 to December 2021.

Declaration of Funding Sources

This research did not receive any specific grant from funding agencies in the public, commercial, or not for profit sectors.

Acknowledgements

We thank the Census and Statistical Department, Government of Hong Kong Special Administrative Region, for providing the micro data sets used in our research. We thank American Journal Experts (AJE) for English language editing on this paper.

REFERENCES

Citation: Zhang Y, Chen FY, Tso GKF, Yang L (2022) Meteorological Factors Facilitated Increased Legionnaires’ Disease Notificationsand Mortality Pre-COVID-19 and COVID-19 in Hong Kong. J Infect Dis treat. 8:40.

Copyright: © 2022 Zhang Y, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited