Impact of regional heterogeneity on the severity of COVID-19

1. Introduction

Coronavirus disease 2019 (COVID-19), which is caused by the SARS-CoV-2 virus, has become a global health threat, imposing a substantial disease burden on our society.¹ Although it has fluctuated, the epidemic has continued.²

COVID-19 is more infectious than other respiratory tract infections,³ such as influenza, and is sometimes fatal to the elderly and those with underlying medical conditions. Considering these characteristics, the capacity of healthcare facilities for COVID-19 patients is an important factor in the management of newly infected COVID-19 patients.

When we consider countermeasures against COVID-19, the proportion of patients with severe disease is extremely important.^4–6 Mild and/or asymptomatic cases can be easily managed as they need no specific treatment. Conversely, severe cases often require intensive, multidisciplinary care. Furthermore, the duration of the disease in severe cases is generally longer than that of other viral pneumonias⁷ and therefore places a greater burden on healthcare capacity. It is important, therefore, to precisely recognize the factors associated with the severity of COVID-19. For example, new variants may contribute to the severity of COVID-19 cases.⁸ However, such variants are not the only determinants of severity; age, sex, past medical history, and many other factors have been attributed to the severity of COVID-19.^5,9,10 In addition, it is conceivable that there are other factors influencing the severity of COVID-19 aside from such biological aspects.

In Japan, the number and proportion of severe COVID-19 cases vary by prefecture. This phenomenon seems difficult to explain; nevertheless, there are differences in the management of COVID-19 cases among prefectures, and these differences may contribute to the capacity of healthcare facilities in each prefecture. Understanding the extent of such regional differences will help us tailor the countermeasures against COVID-19 more appropriately in each prefecture. The main objective of our study was to examine this regional heterogeneity quantitatively.

2. Methods

We included 27,865 cases registered between January 2020 and February 2021 in the COVID-19 Registry of Japan (COVIREGI-JP)¹¹. We used a Bayesian hierarchical model to construct a prefecture-specific random intercept to evaluate the effect of regional heterogeneity.

The capacity of medical facilities can be evaluated by the number of hospital beds; however, the differences in health policies against COVID-19 among local municipalities should also be taken into consideration as the planning of hospital bed distribution and the priority of hospitalization indication are influenced by the policy of each local municipality. Considering this, we constructed two models: one which included the number of beds secured for COVID-19 patients in each prefecture as a fixed effect (model B) and one which did not (model A).

The independent variable was the National Early Warning Score (NEWS)¹² of COVID-19 patients on the day of admission, and the explanatory variables for the fixed effect were sex, age of 65 years and older, and any risk factor of severe infection (cardiovascular disease, cerebrovascular disease, chronic respiratory disease, liver disease, diabetes, hypertension, obesity, dialysis or sever renal failure, and malignancy) for model A and the number of beds secured for COVID-19 patients in each prefecture for model B.

The independent variable, the NEWS score at the time of admission, was assumed to follow a Poisson distribution. We set four separate sampling chains, each consisting of 12,000 random samples, including 1,000 burnin. The sampling convergence was evaluated using the Gelman–Rubin statistics (R-hat below 1.1) and by visually inspecting a trace plot. All analyses were conducted using the R software version 4.0.5.¹³

3. Results

Table 1 presents the demographic characteristics of the hospitalized COVID-19 cases registered in COVIREGI-JP. Tables 2 and 3 present the posterior distribution of the model parameters.

All variables other than the number of beds included as fixed effects exhibited significant associations with the NEWS score at the time of admission in both models.

Figure 1 presents the choropleth graph of the median value of the posterior distribution of prefecture-specific random intercept in model A. Figure 2 presents the choropleth graph of the median value of the posterior distribution of prefecture-specific random intercept in model B.

• Download figure
• Open in new tab

Figure 1. Choropleth map of prefecture-specific random effect on the severity of COVID-19 cases at the time of admission in model A

Thick colors indicate higher risk of higher NEWS score on admission.

Black color (Tokushima prefecture) represents N/A as there are no registered patients in this prefecture.

• Download figure
• Open in new tab

Figure 2. Choropleth map of prefecture-specific random effect on the severity of COVID-19 cases at the time of admission in model B

Thick colors indicate higher risk of higher NEWS score on admission.

Black color (Tokushima prefecture) represents N/A as there are no registered patients in this prefecture.

Several prefectures exhibited a higher risk of higher NEWS scores upon admission, but these prefectures did not significantly change between model A and B. The Bayesian factor of model A over model B was 136988.360.

4. Discussion

The results revealed the existing regional heterogeneity in the severity of COVID-19 cases at the time of admission in Japan. Our results suggest that the health policy of each prefecture might have an impact on the severity of hospitalized COVID-19 cases on admission.

A wide variety was observed in the median values of the posterior distribution of random effects for the different prefectures. That of Kanagawa prefecture in model A was 0.340 [95% Credible Interval (CrI) 0.269, 0.410], and that of Iwate prefecture was -0.431 [95% CrI -0.593, -0.268], indicating a substantial discrepancy in the severity of hospitalized COVID-19 cases between these two prefectures.

The posterior distribution of random effects in model B exhibited a similar trend to that in model A. The coefficient of Kanagawa prefecture in model B was 0.306 [95% CrI 0.213, 0.402], and that of Iwate prefecture was -0.413 [95% CrI -0.579, -0.249]. Other prefectures demonstrated similar results in both models, suggesting that factors other than the number of beds secured for COVID-19 patients (e.g., prioritization policy for COVID-19 case admission and number of physicians) had an impact on the severity of COVID-19 on admission.

We should take various factors into consideration when we interpret our results. The potential factors associated with the variation of the severity on admission could be viewed from two directions, namely, the local epidemiology of COVID-19 and the COVID-19 healthcare policy in each prefecture. From the former viewpoint, a regional surge in the number of COVID-19 patients or a rapid increase in the severity of disease due to new SARS-CoV-2 variants may contribute to the change in the severity on admission.^8,14 However, our analysis did not include patients with new SARS-CoV-2 variants as the study period was limited to before March 2021. While it is possible that COVID-19 cases due to new variants carry a greater risk of serious illness, there was a substantial regional heterogeneity in the severity of COVID-19 cases at the time of admission, as previous studies have demonstrated.^15–17 Presumably, variants are contributing factors to the severity of COVID-19 cases, but it should be noted that other factors also have a non-negligible influence on the severity of disease.

The latter point is relevant to policymakers (and healthcare providers) in the municipalities, mainly related to the on-site indication for hospitalization that is operated in each prefecture. The policy implemented to manage the waves of infection could be affected by the total medical capacity in the area. The number of beds secured for COVID-19 patients is one of the major resources.^15,18 The number of beds per population in each prefecture was included as one of the fixed effects in model B, and the random effect of each prefecture did not drastically change. Also, the number of physicians eligible for COVID-19 management may have an impact as this factor could create a bottleneck for patient acceptance.¹⁹ However, it is difficult to precisely grasp the number of physicians who are engaging in COVID-19 patient management. This is because the specialty of the medical doctor is not a good proxy for the number of physicians engaging in COVID-19 management under such an emergency period.²⁰

The local government of each prefecture has established policy for the process of admission prioritization as a way of integrating the information from the regional medical situations, including the aforementioned factors, a strategy that may determine the variation in the severity on admission.^21,22 This means that there may be differences in whether a patient can be hospitalized and how long it takes for a patient to be hospitalized, even for clinically similar cases, depending on the prefecture they belong to.

Two types of strategy for symptomatic COVID-19 patients can be assumed. The following are simplified examples: 1) moderate to severe cases are admitted to hospital, whereas mild cases are managed at home or in a dedicated facility (severity on admission will increase), and 2) all patients, whether with mild or severe disease, are admitted to hospital (severity on admission will decrease). Both strategies have their own risks and benefits, and the local government of each prefecture has adopted each strategy depending on the height of each wave of the epidemic.

The strategy adopted should be evaluated by the outcomes of COVID-19 patients, damage to the medical system for non-COVID-19 diseases, and the impact on hospital and regional healthcare economy. A previous study in China demonstrated that the number of beds and medical staff per COVID-19 patient was negatively correlated with mortality.²³ Conversely, a recent study from Italy reported that greater access to hospitals did not affect the mortality of COVID-19 patients older than 80 years but affected that of patients younger than 80 years,²⁴ suggesting that the best approach may be to change the strategy according to age group or other properties of the patients. Evaluation of the regional policy could provide important insights to construct a better healthcare system against COVID-19 and other emerging infectious diseases.

Our analysis had several limitations. First, we examined data solely from Japan; thus, it is not appropriate to generalize our results globally. Nevertheless, it is easy to imagine that, even in other countries, there would be differences in health policies among local authorities and that they might have an impact on the management of COVID-19 cases. Second, we could not include other confounding factors, such as human resources and variance of SARS-CoV-2; these aspects should be addressed in future studies.

In conclusion, our analysis revealed a possible association between prefecture and an increased/decreased risk of severe COVID-19 infection at the time of admission. In addition, the number of beds secured for COVID-19 patients in each prefecture is not a single cause of such heterogeneity; therefore, other factors could be significant for the management of COVID-19 cases in Japan. Countermeasures against COVID-19 will be more appropriate if we take these insights into consideration.

Author contributions

Shinya Tsuzuki: Data Curation, Formal Analysis, Methodology, Visualization, Original Draft Preparation. Yusuke Asai: Data Curation, Methodology, Validation, Review and Editing. Nobuaki Matsunaga: Original Draft Preparation, Review and Editing. Haruhiko Ishioka: Original Draft Preparation, Review and Editing. Takayuki Akiyama: Data curation, Validation, Review and Editing. Norio Ohmagari: Conceptualization, Funding Acquisition, Project Administration, Review & Editing.

Data availability

The data that support the findings of this study are available upon request to the corresponding author. The data are not publicly available due to privacy or ethical restrictions.

Funding

This work was supported by a grant from the Health and Labor Sciences Research, “Research for risk assessment and implementation of crisis management functions for emerging and re-emerging infectious diseases.”

Conflict of interest

None declared.

Ethical approval

This study was approved by the National Center for Global Health and Medicine (NCGM) ethics review committee (NCGM-G-003494-0). Information regarding opting out of our study is available on the registry website.

Patient consent

Informed consent was obtained in the form of opt-out on the registry website.

Permission to reproduce material from other sources

Not applicable.