Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Shifts in the epidemic season of human respiratory syncytial virus associated with inbound overseas travelers and meteorological conditions in Japan, 2014–2017: An ecological study

  • Keita Wagatsuma ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Writing – original draft

    waga@med.niigata-u.ac.jp

    Affiliation Division of International Health (Public Health), Graduate School of Medical and Dental Sciences, Niigata University, Niigata, Japan

  • Iain S. Koolhof,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation College of Health and Medicine, School of Medicine, University of Tasmania, Hobart, Tasmania, Australia

  • Yugo Shobugawa,

    Roles Conceptualization, Investigation, Writing – review & editing

    Affiliation Department of Active Aging, Graduate School of Medical and Dental Sciences, Niigata University, Niigata, Japan

  • Reiko Saito

    Roles Conceptualization, Investigation, Project administration, Writing – review & editing

    Affiliations Division of International Health (Public Health), Graduate School of Medical and Dental Sciences, Niigata University, Niigata, Japan, Infectious Diseases Research Center of Niigata University in Myanmar (IDRC), Graduate School of Medical and Dental Sciences, Niigata University, Niigata, Japan

Abstract

Few studies have examined the effects of inbound overseas travelers and meteorological conditions on the shift in human respiratory syncytial virus (HRSV) season in Japan. This study aims to test whether the number of inbound overseas travelers and meteorological conditions are associated with the onset week of HRSV epidemic season. The estimation of onset week for 46 prefectures (except for Okinawa prefecture) in Japan for 4-year period (2014–2017) was obtained from previous papers based on the national surveillance data. We obtained data on the yearly number of inbound overseas travelers and meteorological (yearly mean temperature and relative humidity) conditions from Japan National Tourism Organization (JNTO) and Japan Meteorological Agency (JMA), respectively. Multi-level mixed-effects linear regression analysis showed that every 1 person (per 100,000 population) increase in number of overall inbound overseas travelers led to an earlier onset week of HRSV epidemic season in the year by 0.02 week (coefficient –0.02; P<0.01). Higher mean temperature and higher relative humidity were also found to contribute to an earlier onset week by 0.30 week (coefficient –0.30; P<0.05) and 0.18 week (coefficient –0.18; P<0.01), respectively. Additionally, models that included the number of travelers from individual countries (Taiwan, South Korea, and China) except Australia showed that both the number of travelers from each country and meteorological conditions contributed to an earlier onset week. Our analysis showed the earlier onset week of HRSV epidemic season in Japan is associated with increased number of inbound overseas travelers, higher mean temperature, and relative humidity. The impact of international travelers on seasonality of HRSV can be further extended to investigations on the changes of various respiratory infectious diseases especially after the coronavirus disease 2019 (COVID-19) pandemic.

Introduction

Human respiratory syncytial virus (HRSV) is an infection of the respiratory tract which causes clinically severe pneumonia in children and bronchiolitis in infants [1]. Globally, acute lower respiratory tract infections caused by HRSV lead to approximately 70,000 annual deaths in children under the age of 5 years and approximately 3.4 million people require hospitalization [2]. Epidemics of HRSV are seasonal, typically having an epidemic during winter in temperate zones as seen in Japan [3]. The HRSV epidemic season in Japan normally begins around October, peaks in December-January, and ends in March-April. However, a shift in the HRSV epidemic season to summer periods (from June-August) has been reported in Japan more recently [4]. Although many previous studies looked at the association between HRSV cases and meteorological conditions [3,511], there may well be other factors that determine the onset week of HRSV epidemic season.

In recent years, with a large number of travelers coming into Japan, the likelihood for contact the local population have with foreign travelers has increased [12]. Until the start of coronavirus disease 2019 (COVID-19) pandemic in 2020, the number of inbound overseas travelers increased in the last several decades and was expected to continue to increase in the future. Several studies have suggested the possibility of increased transmission of infectious diseases owing to greater international travel and arrivals [1316], suggesting that the inflow of infected travelers may be associated with the epidemics. Indeed, we recently reported that early cases of COVID-19 in January-February 2020 in Japan were concentrated in prefectures where the numbers of foreign travelers were high [17].

This study aims to test whether the number of inbound overseas travelers and meteorological conditions are associated with the onset week of HRSV epidemic season in Japan. We included the number of inbound overseas travelers as a parameter as well as meteorological conditions, to examine the associations with the onset week of HRSV epidemic season during 2014–2017.

Materials and methods

Study design

This study was designed to evaluate the ecological association between meteorological (mean temperature and relative humidity) conditions and the number of inbound travelers on the onset week of HRSV epidemic season in Japan over a 4-year period (2014–2017).

Study area

The study area is Japan, which consists of 47 prefectures. Of the 47 prefectures, 46 prefectures, excluding Okinawa prefecture, belong to the temperate zones [3]. Okinawa belongs to the subtropical region where the seasonality of HRSV is different [18]. Therefore, data for this study, only included the 46 prefectures which have temperate climates and excluded the subtropical Okinawa prefecture.

HRSV national surveillance data in Japan

HRSV epidemiological data, available from the Infectious Disease Weekly Report (IDWR), sourced from the National Epidemiological Surveillance of Infectious Diseases (NESID), published by National Institute of Infectious Diseases (NIID) under the Ministry of Health, Welfare and Labour in Japan (MHLW) was used for the analysis in this paper [19]. The MHLW designates approximately 3,000 paediatric sentinel sites (i.e., hospitals and clinics) in Japan, which report numbers of patients diagnosed as HRSV infection at weekly basis to the prefecture or municipal public health sectors in Japan [18]. An HRSV case is defined by a positive rapid diagnostic test (RDT) kits licenced in Japan, or a laboratory confirmation such as virus isolation or antibody rise in paired sera according to the MHLW guidelines [20]. The prefectural data is reported to the NIID, and the collective numbers of HRSV cases per prefecture are released weekly through its website [19]. However, the information at the individual paediatric sentinel sites where each HRSV case count is reported is anonymised and not made public. Therefore, we could not specify whether these NIID data came from the same sentinel sites during the study period. Still we assume these changes of sentinel sites were minor and did not affect the quality of data. A sentinel site is located roughly over 30,000 populations according to the Infectious Disease Control Law in Japan and the prefecture government are responsible for maintaining the sites so that it does not affect the quality of data [21].

In this study, we extracted the total number of HRSV cases in each prefecture (excluding Okinawa prefecture) reported during the 1st week of 2014 and the 52nd week of 2018 from the NESID data by weekly basis. Since the number of sentinel sites differed by prefecture largely due to population size in the area, then the numbers of HRSV cases per sentinel site (numbers of patient visits due to HRSV cases per hospital/clinic) was calculated from the total number of cases divided by the number of sentinel sites. Using the number of weekly HRSV cases per sentinel site, we created the epidemic curve of HRSV in Japan from 2014 to 2017.

Definition of onset week of HRSV epidemic season

We used the onset week of HRSV epidemic season in 46 prefectures for 4-year period estimated by Yamagami et al. in the other paper [4]. They developed a new algorithm to estimate the start of the HRSV epidemic season based on the IDWR data, since NIID does not set the definition on the start of HRSV epidemic season in Japan. In their paper, in determining the epidemic starting point of the HRSV epidemic season in each prefecture, they make estimates that take into account both the number of HRSV reports included within the epidemic period relative to the total number of HRSV reports (capture rate) and the number of HRSV reports per week within the epidemic period (HRSV-reports/w). They calculate both these capture rates and the HRSV-reports/w search index to determine the starting point of the epidemic cycle (see below Eqs (13)).

In the following, Yamagami et al. [4] define the parameters for A(x) denotes HRSV-reports/w (the number of HRSV reported per week); B(x) denotes capture rate (the number of HRSV reports included within onset-trough periods relative to the total number of HRSV reports in the dataset); Index(x) denotes the search index. The capture rate A(x) and HRSV-reports/w (B(x)) were converted into the same scale and summed; x is a percentile value ranging from 0 to 100 in integer units; s is the number of epidemic cycles in the dataset; Onset(x) is the ordinal number of each onset week within the dataset from 1 to x in integer units; Trough is the ordinal number of each trough-week within the dataset from 0 to x in integer units; k is the number of HRSV reports in a given week.

(1)(2)(3)

The procedure for setting up the onset week of HRSV epidemic season by Yamagami et al. [4] was defined as follows. First, they transposed the number of HRSV reports in each prefecture to the percentile rank and selected the lowest percentile rank that was greater than the percentile value. Secondly, they selected the 101st percentile rank in each prefecture for the 101st percentile value. Third, they defined the reporting week corresponding to the intersection of these percentile ranks and the increase in the slope of the sine curve as the onset week. Finally, they applied this generation procedure to all 47 prefectures in Japan to estimate the onset week at prefectural level. The validity of the determined onset was confirmed by calculating the capture rate. The validity of the model was subsequently tested using data from 47 prefectures in the previous study [4].

In this study, the onset weeks of HRSV epidemic season by prefecture over the 4-year period from 2014 to 2017 were extracted from their paper [4] and used in the analysis.

Dependent variable

The dependent variable for the models is the onset week of HRSV epidemic season, ranged from weeks 1–52, by each prefecture [4].

Explanatory variables

Meteorological data.

The meteorological data published by the Japan Meteorological Agency (JMA) [22] was used as explanatory variables for the models in this paper. Monthly and yearly meteorological data observed at a meteorological observatory situated at the prefectural capital city was used for each respective prefecture. We extracted monthly and yearly mean temperature and relative humidity from 2014 to 2017. There were two meteorological observatories with missing data on relative humidity. Therefore, we selected the observatories sites with the nearest distance from the prefectural capital by substitution. Using these monthly meteorological data in each prefecture, mean temperature and relative humidity over the whole of Japan have been calculated.

Inbound overseas travelers.

As a surrogate of the monthly and yearly number of overall inbound overseas travelers by prefecture, the daily number of foreign visitors accommodated (e.g., hotels) in each prefecture from 2014 to 2017 were sourced from the travel statistics survey of the Japan National Tourism Organization (JNTO) [23]. In addition, we used yearly number of travelers by country from the top three countries (China, South Korea, and Taiwan) with the most travelers coming into Japan. In addition, Australia was used as the country with the largest number of travelers from the southern hemisphere from 2014 to 2017.

Ethical considerations

This ecological study meets the ethical and regulatory guidelines, including adherence to the legal requirements of the study country. The collection of HRSV data under the NESID is legally defined and implemented under the Infectious Disease Control Law of the Japanese government. Through this scheme, the requirement for informed consents from the patients were legally waived. The MHLW publishes these data as a web-based open source data [19]. Similarly, the meteorological data and the inbound overseas traveler’s data are also open to the public from the respective governmental website [22,23].

Statistical analysis

The statistical analyses were conducted in two steps. Firstly, before the main analysis, the association between each variable included in this study, onset week of HRSV epidemic season, mean temperature, relative humidity, number of overall inbound overseas travelers per 100,000 population by each prefecture in Japan, number of inbound overseas travelers from Taiwan, South Korea, China, and Australia per 100,000 populations by each prefecture in Japan, was calculated using Spearman’s rank test to evaluate the correlation between variables (N = 4×46 = 184). Subsequently, variables that showed strong multicollinearity were excluded in the regression model.

Secondly, multi-level mixed-effects linear regression model was performed to adjust for the effects of confounding variables (N = 4×46 = 184). We conceptualised our analysis with a multi-level structure consisting of factors nested within prefectures (prefectural-level) to account for the impact of year repeated measures. We fitted the data using a multi-level linear regression procedure with a mixed effect model, adjusting for both individual and prefectural levels as fixed effects. The parameters of the model were chosen to be random effects estimators using the maximum likelihood. The onset week of HRSV epidemic season by each year (continuous) was the dependent variable, and yearly number of inbound travelers (continuous) and meteorological (yearly mean temperature and relative humidity) conditions (continuous) were included as explanatory variables. Furthermore, the model was adjusted by inputting year variables (years 2014, 2015, 2016, and 2017) (dummy variables) as a covariate to control seasonal variation. Our first analysis involved the estimation of a null model (Model 1). The null model allows us to decompose the variance in the onset week of HRSV epidemic season to determine whether it is due to variation at the prefectural or individual level. Model 2 represents the number of overall inbound travelers per 100,000 population. Models 3, 4, 5, and 6 consider the numbers of inbound travelers per 100,000 population from Taiwan, South Korea, China, and Australia, respectively. The regression coefficients and P-values were calculated in each multi-level analysis. Goodness-of-fit statistics of each model was evaluated by calculating the Akaike’s Information Criterion values (AIC). The intra-class correlation (ICC) was calculated to evaluate the similarity within group.

All tests were two-sided and P-values less than 0.05 were considered statistically significant. All data analyses were performed using STATA software, Ver. 15.0 (Stata Corp LP College Station, TX, USA).

Results

Seasonality trends of HRSV cases, mean temperature, relative humidity, and inbound overseas travelers throughout Japan, 2014–2017

HRSV activity in Japan usually peaks in the winter season (from December-January). The peak epidemic of HRSV was observed in December in 2014 and 2015, but was changed to October in 2016, and to September in 2017 (Fig 1A, S1 Dataset). Correspondingly, the annual onset week of the HRSV epidemic season occurred earlier in the later year. The median onset weeks of the HRSV epidemic season were around 36.5 weeks (September) in 2014, 36.0 weeks (September) in 2015, 35 weeks (late August) in 2016, and 30 weeks (beginning of July) in 2017, respectively. The median onset week of HRSV epidemic season was approximately 6.0 weeks (1.5 months) earlier in 2017 than in 2014. On the other hand, the number of inbound overseas travelers increased year on year, doubled from 13,413,467 in 2014 to 28,691,073 in 2017 (S2 Dataset). However, the monthly number of inbound overseas travelers from all countries during 2014 and 2017 did not show clear seasonal trend (Fig 1B and S2 Dataset). The number of inbound overseas travelers from Taiwan, China and South Korea remained almost steady through the year, while Australia is observed to have a higher number of inbound travelers in the winter months of January and December (Fig 1C–1F, and S3S6 Datasets). This pattern did not differ during the four years. The mean temperature peaks at approximately 27°C in August of each year and the relative humidity is around 60–70% throughout the year in Japan (Fig 1G and 1H, S7 and S8 Datasets). Mean temperature and relative humidity did not vary significantly throughout the 4-year period.

thumbnail
Fig 1. Trends of HRSV cases, mean temperature, relative humidity, and inbound overseas travelers throughout Japan, 2014–2017.

(A) Epidemic curve of HRSV in Japan from 2014 to 2017. (B) Monthly number of overall inbound overseas travelers from all countries during 2014 and 2017, (C) Taiwan, (D) South Korea, (E) China, (F) Australia, (G) Monthly mean temperature (°C) throughout Japan from 2014 to 2017, (H) Monthly average relative humidity (%) from 2014 to 2017. The orange, yellow, green and brown lines represent 2014, 2015, 2016 and 2017, respectively.

https://doi.org/10.1371/journal.pone.0248932.g001

Spearman’s correlation analysis of variables included in the study

Table 1 shows the matrix of Spearman’s correlation coefficients [ρs] for variables included in the study for 4-year period (S9 Dataset). Onset week of HRSV epidemic season and mean temperature showed significant inverse correlation (Spearman’s ρ = –0.17). Onset week of HRSV epidemic season was also significant having an inverse correlation with the number of overall inbound travelers per 100,000 population and number of inbound travelers from individual country per 100,000 population (Spearman’s ρ = –0.20, –0.16, –0.32, –0.16 and –0.15, respectively). These results suggested the increased temperature and the higher number of overseas travelers are associated to earlier onset of HRSV season, respectively.

thumbnail
Table 1. Spearman’s correlation of variables included in the study (N = 184).

https://doi.org/10.1371/journal.pone.0248932.t001

Association between meteorological conditions and number of inbound travelers on onset week of HRSV epidemic season in Japan by multi-level analysis

The multi-level analysis showed that mean temperature, relative humidity, and the number of inbound travelers were negatively associated with the onset week of the HRSV epidemic season (Table 2, S9 Dataset). The null model with no predictors (Model 1) revealed no significant evidence of inter-prefectural variation (prefectural-level variance = 0.61, ICC = 0.01). In Model 2, which included the number of overall inbound travelers, we found that every 1 degree increases in mean temperature led to an earlier onset week within the year by 0.30 week (Model 2: coefficient –0.30; P<0.05). In addition, every 1.0% increase in relative humidity and every 1 person (per 100,000 populations) increase in the number of overall inbound travelers were associated with an earlier onset week of the HRSV epidemic season by 0.18 week and 0.02 week respectively (Model 2: coefficient –0.18; P<0.01 and coefficient –0.02; P<0.01). Even after adjusting for meteorological variables and year variables, travelers’ visits were independently associated with an earlier onset week of HRSV epidemic season. An increase in the number of inbound travelers from individual countries was associated with an earlier onset week of HRSV epidemic season such as Taiwan (Model 3: coefficient –0.07; P<0.01), South Korea (Model 4: coefficient –0.25; P<0.01), and China (Model 5: coefficient –0.09; P<0.05). On the other hand, an increase in the number of inbound travelers from Australia (Model 6: coefficient –0.60, P = 0.05) was marginally associated with an earlier onset week of the HRSV epidemic season. The onset week of HRSV epidemic was significantly earlier in 2017 and 2016 than in 2014 in all models. To this end, onset week of HRSV epidemic season was earlier in all models from 2014 to 2017.

thumbnail
Table 2. Regression coefficients with P-values for onset week of HRSV epidemic season: The results of multi-level mixed-effects linear regression analysis (N = 184).

https://doi.org/10.1371/journal.pone.0248932.t002

Discussion

In this study, we examined the associations between meteorological (mean temperature and relative humidity) conditions as well as the number of inbound overseas travelers, and the onset week of the HRSV epidemic season in Japan. Our findings suggest that these variables have a significant association with the onset week of HRSV epidemic season in Japan, and an increase of one unit in temperature, relative humidity and inbound travelers, are associated with 0.3, 0.2, 0.02 week earlier onset week in the year, respectively by multi-level analysis. To our knowledge, this is the first study to evaluate the effects of inbound overseas travelers and meteorological conditions on the shift in HRSV epidemic season in Japan.

The HRSV epidemic period historically has been between September and December in Japan [18], however, recently the onset week of HRSV epidemic season has begun to shift, becoming earlier in the year since 2016 [4]. Studies in the past have focused on the epidemic peak and the size (e.g. number of HRSV cases) [2426]. Although the number of HRSV cases is an important indicator, it is more useful to predict the onset week before the number of infections to increase. In a clinical setting, it is crucial to identify the start of an epidemic to allow paediatricians to assess whether or when to initiate palivizumab (anti-HRSV antibody) administration to prevent possible severe cases especially in high-risk children [27]. Our analysis suggests that not only meteorological conditions but also inbound overseas travelers may be associated with changing in the timing of HRSV epidemics, highlighting the need to consider these impacts in the timing of palivizumab administration. In April 2019, the Japan Pediatric Society published the consensus guidelines on the use of Palivizumab in Japan [28]. This guideline states that due to the shift of seasonality of HRSV in recent years in Japan, in order to increase the efficacy of palivizumab, it is necessary to raise the serum antibody titer to the level required for prophylaxis by the start of the onset week of HRSV epidemic season, and that pediatricians and others should lead a study in each prefecture to consider the month to start palivizumab and the duration and frequency of administration. In particular, they reported the shift of seasonality of HRSV in recent years in Japan, but they mentioned possible association with meteorological conditions as contributing factors. Therefore, our findings may bring a new information to the guideline that HRSV season may come earlier where they have more overseas travelers in the local area.

First, we examined the association between onset week of HRSV epidemic season and meteorological conditions. In previous studies, many researchers have reported on the association between meteorological conditions and HRSV transmission [3,511,13,24,29,30]. Most of the studies showed that higher temperature and humidity were associated with greater HRSV transmission in tropical and subtropical regions [3,69,11]. In contrast, in temperate zones, lower temperatures and higher humidity were associated with greater HRSV transmission [5]. In our results, we showed that higher temperature and/or higher humidity were associated with an earlier onset week of HRSV epidemic season. We previously reported that the number of HRSV cases in summer (from June-August) in Japan increased with higher temperature and higher relative humidity during 2007–2014 [26]. Our findings here also indicate that meteorological conditions have an association with HRSV transmission, corresponding with an early onset in the epidemic season. We observed a shift in the timing of the onset week of HRSV epidemic season, whereby, the epidemic shifted to occurring earlier in the year during summer. Global warming and consequences from climate change may affect the onset timing of HRSV epidemics. Meteorological condition in countries like Japan (i.e., temperate zones), which are near to the subtropical zone, might be affected by such subtle changes in climate.

Secondly, we considered the association between onset week of HRSV epidemic season and inbound overseas travelers to Japan. As the number of inbound overseas travelers has increased in recent years, there is a possibility that person-to-person close contact has increased, especially in densely populated cities such as urban areas (e.g., Tokyo, Osaka, Aichi, Fukuoka, and Hokkaido). An inflow of overseas travelers may affect the HRSV epidemic season owing to the infection route of HRSV having human-to-human transmission via droplets and contact with increased travelers resulting in greater contact rates [31]. In influenza, the travel has been reported to have an impact on the spread of infectious diseases [15,16], similarly, measles and novel coronavirus disease 2019 (COVID-19) have also been reported to be increased in transmission resulting from a greater traffic of people [13,14,17]. It is possible that the inflow of travelers infected with HRSV may be associated with the HRSV epidemic season in Japan based on what was observed with the other infectious diseases. The results of this study suggest that the inflow of overseas travelers may have important impact on the outbreak of HRSV in Japan. However, there are few studies looking at HRSV transmission and the association with population behavioural and movement dynamics in Japan. In our study, after controlling for mean temperature, relative humidity, and year, we showed that the onset week of HRSV epidemic season was earlier with a greater number of inbound overseas travelers into Japan. Therefore, the onset week of HRSV epidemic season may have shifted due to the inflow of people from areas with different epidemic periods compared to that in Japan.

In this analysis, we expected that the impact on the onset week of HRSV epidemic season might be affected because of seasonality differences of HRSV in Taiwan and Australia compared to Japan. Epidemics of the HRSV often occur from June-September in Taiwan [29,32,33] and from June-September in temperate climate areas such as Sydney and Melbourne in Australia [8,11,34,35]. In contrast, HRSV epidemic tends occur from November-January in Korea and China as also typically seen before 2016 in Japan [30,3639]. Thus, we expected that Taiwan and Australia would have a greater effect on causing an earlier onset week of HRSV epidemic season in Japan than Korea and China. However, all countries except Australia showed association with earlier onset of HRSV epidemic season in Japan. The marginal effect of Australia was presumably from the differences between the peak seasons for travelers to Japan (January and December) and that for HRSV in Australia (June-September). Further studies are needed to investigate the effect seasonality by the original location of where the travelers’ have departed from and the timing of the HRSV epidemic season.

Several limitations of this study need to be considered. First, since this is an ecological study for limited years, the longitudinal monitoring and evaluation may be leading to a different result. Secondly, the travelers’ data in this study does not take into account which local area of the country travelers have come from. The country areas are very large in China and Australia and have multiple climate zones. Therefore, the HRSV epidemic periods in each of these countries differ from area to area domestically, and the impact on the onset week of HRSV epidemic season may also differ. Thirdly, this study focused on onset week and did not analyse the peak week and size. Fourthly, this study does not consider the lag-time of inbound overseas travelers and meteorological conditions. The reason for this is that in the present analysis, we used 46 prefectures multiplied by 4 seasons (184 units), and used an annual bases for inbound overseas travelers and meteorological conditions instead of weekly bases. Therefore, a more detailed time-series level of analysis needs to be considered in our future studies. Finally, this study does not take into account on the numbers of outgoing Japanese travelers or changes in the domestic human mobility (e.g., inflow-outflow population between prefectures) over time. This may be a contributing factor in the changes of HRSV seasonality.

In conclusion, the earlier onset week of HRSV epidemic season in Japan is ecologically associated with an increased number of inbound overseas travelers, mean temperature, and relative humidity. We highlight the hypothesis that inbound overseas travelers may contribute to the shift in the HRSV epidemic season in addition to changes in meteorological conditions. The current findings provide useful information for public health decision making in prevention and the preparation for HRSV outbreaks in Japan. After the emergence of COVID-19, respiratory viral infections underwent significant changes. The number of HRSV cases has dramatically dropped and virtually no epidemics were observed except for small outbreaks in southernmost prefectures in Japan, such as Okinawa in November and Kagoshima in October in 2020 [19]. The dramatic decline of overseas travelers, strengthened personal hygiene behaviours and social infection control measurements may have contributed to such a decline of the HRSV cases in 2020 [4045]. Thus it is important that future studies include post- 2020 data to assess how the COVID-19 emergence changes the onset week of HRSV epidemic season. In particular, taking into account the impact of travel restrictions on the number of inbound overseas travelers, it may be possible to identify changes in the indirect impact on the onset week of HRSV epidemic season. Further evaluation and monitoring toward the change of respiratory infectious diseases after the COVID-19 pandemic should be jointly investigated at a global level by international coordinated schemes.

Supporting information

S1 Dataset. The weekly number of HRSV cases per sentinel site in Japan.

https://doi.org/10.1371/journal.pone.0248932.s001

(XLSX)

S2 Dataset. The monthly number of overall inbound overseas travelers from all countries to Japan.

https://doi.org/10.1371/journal.pone.0248932.s002

(XLSX)

S3 Dataset. The monthly number of inbound overseas travelers from Taiwan to Japan.

https://doi.org/10.1371/journal.pone.0248932.s003

(XLSX)

S4 Dataset. The monthly number of inbound overseas travelers from South Korea to Japan.

https://doi.org/10.1371/journal.pone.0248932.s004

(XLSX)

S5 Dataset. The monthly number of inbound overseas travelers from China to Japan.

https://doi.org/10.1371/journal.pone.0248932.s005

(XLSX)

S6 Dataset. The monthly number of inbound overseas travelers from Australia to Japan.

https://doi.org/10.1371/journal.pone.0248932.s006

(XLSX)

S7 Dataset. Dataset for the monthly mean temperature throughout Japan.

https://doi.org/10.1371/journal.pone.0248932.s007

(XLSX)

S8 Dataset. Dataset for the monthly average relative humidity throughout Japan.

https://doi.org/10.1371/journal.pone.0248932.s008

(XLSX)

S9 Dataset. Dataset of variables for Spearman’s correlation analysis and multi-level analysis.

https://doi.org/10.1371/journal.pone.0248932.s009

(XLSX)

Acknowledgments

We thank local governments, public health centers and institutes, and the National Institute of Infectious Diseases, Japan, for surveillance, laboratory testing, epidemiological investigations, and data collection. We thank Dr. Naohito Tanabe from the University of Niigata Prefecture in Japan for his useful comments. We would also like to thank Editage (www.editage.com) for English language editing.

References

  1. 1. Piedimonte G, Perez MK. Respiratory syncytial virus infection and bronchiolitis. Pediatr Rev. 2014;35(12):519–30. pmid:25452661.
  2. 2. Nair H, Nokes DJ, Gessner BD, Dherani M, Madhi SA, Singleton RJ, et al. Global burden of acute lower respiratory infections due to respiratory syncytial virus in young children: a systematic review and meta-analysis. Lancet. 2010;375(9725):1545–55. pmid:20399493.
  3. 3. Tang JW, Loh TP. Correlations between climate factors and incidence—a contributor to RSV seasonality. Rev Med Virol. 2014;24(1):15–34. pmid:24421259.
  4. 4. Yamagami H, Kimura H, Hashimoto T, Kusakawa I, Kusuda S. Detection of the Onset of the Epidemic Period of Respiratory Syncytial Virus Infection in Japan. Front Public Health. 2019;7:39. pmid:30931290.
  5. 5. Price RHM, Graham C, Ramalingam S. Association between viral seasonality and meteorological factors. Sci Rep. 2019;9(1):929. pmid:30700747.
  6. 6. Chen Z, Zhu Y, Wang Y, Zhou W, Yan Y, Zhu C, et al. Association of meteorological factors with childhood viral acute respiratory infections in subtropical china: An analysis over 11 years. Archives of Virology. 2014;159(4):631–9. pmid:24114148.
  7. 7. Tang JW, Lai FY, Wong F, Hon KL. Incidence of common respiratory viral infections related to climate factors in hospitalized children in Hong Kong. Epidemiol Infect. 2010;138(2):226–35. pmid:19631018.
  8. 8. Morley C, Grimwood K, Maloney S, Ware RS. Meteorological factors and respiratory syncytial virus seasonality in subtropical Australia. Epidemiol Infect. 2018;146(6):757–62. pmid:29560833.
  9. 9. Tian DD, Jiang R, Chen XJ, Ye Q. Meteorological factors on the incidence of MP and RSV pneumonia in children. PLoS One. 2017;12(3):e0173409. pmid:28282391.
  10. 10. Haynes AK, Manangan AP, Iwane MK, Sturm-Ramirez K, Homaira N, Brooks WA, et al. Respiratory syncytial virus circulation in seven countries with Global Disease Detection Regional Centers. J Infect Dis. 2013;208:246–54. pmid:24265484.
  11. 11. Paynter S, Ware RS, Sly PD, Weinstein P, Williams G. Respiratory syncytial virus seasonality in tropical Australia. Aust N Z J Public Health. 2015;39(1):8–10. pmid:25648729.
  12. 12. Henderson JC. Destination Development: Trends in Japan’s Inbound Tourism. Int J Tourism Research. 2017;19(1):89–98.
  13. 13. Cliff A, Haggett P. Time, travel and infection. Br Med Bull. 2004;69:89–99. pmid:15226199.
  14. 14. Sasaki A, Suzuki H, Sakai T, Sato M, Shobugawa Y, Saito R. Measles outbreaks in high schools closely associated with sporting events in Niigata, Japan. J Infect. 2007;55(2):179–83. pmid:17560656.
  15. 15. Belderok SM, Rimmelzwaan GF, van den Hoek A, Sonder GJ. Effect of travel on influenza epidemiology. Emerg Infect Dis. 2013;19(6):925–31. pmid:23735636.
  16. 16. Htwe KTZ, Dapat C, Shobugawa Y, Odagiri T, Hibino A, Kondo H, et al. Phylogeographic analysis of human influenza A and B viruses in Myanmar, 2010–2015. PLoS One. 2019;14(1):e0210550. pmid:30629691.
  17. 17. Wagatsuma K, Phyu WW, Osada H, W. Tang J, Saito R. Geographic correlation between the number of COVID-19 cases and the number of overseas travelers in Japan, Jan-Feb, 2020. Jpn J Infect Dis. 2020 (in press). pmid:32863355.
  18. 18. Hibino A, Saito R, Taniguchi K, Zaraket H, Shobugawa Y, Matsui T, et al. Molecular epidemiology of human respiratory syncytial virus among children in Japan during three seasons and hospitalization risk of genotype ON1. PLoS One. 2018;13(1):e0192085. pmid:29377949.
  19. 19. National Institute of Infectious Disease (NIID). Infectious Diseases Weekly Report (IDWR). Available from: https://www.niid.go.jp/niid/ja/idwr.html. (accessed Mar 1, 2021).
  20. 20. The Ministry of Health, Labor and Welfare, Japan (MHLW). Notification of physicians and veterinarians under the Infectious Diseases Act (RS virus). Available from: http://www.mhlw.go.jp/bunya/kenkou/kekkaku-kansenshou11/01-05-15.html. (accessed Aug 1, 2020).
  21. 21. National Institute of Infectious Diseases (NIID). National Epidemiological Surveillance for Infectious Diseases (NESID) in Japan. Available from: https://www.niid.go.jp/niid/images/epi/nesid/nesid_ja.pdf. (accessed Mar 1, 2021).
  22. 22. Japan Meteorological Agency (JMA). Meteorological statistics. Available from: http://www.data.jma.go.jp/obd/stats/etrn/index.php. (accessed Aug 2, 2020).
  23. 23. Japan National Tourism Organization (JNTO). Japan tourism statistics. Available from: https://statistics.jnto.go.jp/en/. (accessed Aug 3, 2020).
  24. 24. Nyoka R, Omony J, Mwalili SM, Achia TNO, Gichangi A, Mwambi H. Effect of climate on incidence of respiratory syncytial virus infections in a refugee camp in Kenya: A non-Gaussian time-series analysis. PLoS One. 2017;12(6):e0178323. pmid:28570627.
  25. 25. Bloom-Feshbach K, Alonso WJ, Charu V, Tamerius J, Simonsen L, Miller MA, et al. Latitudinal variations in seasonal activity of influenza and respiratory syncytial virus (RSV): a global comparative review. PLoS One. 2013;8(2):e54445. pmid:23457451.
  26. 26. Shobugawa Y, Takeuchi T, Hibino A, Hassan MR, Yagami R, Kondo H, et al. Occurrence of human respiratory syncytial virus in summer in Japan. Epidemiol Infect. 2017;145(2):272–84. pmid:27682641.
  27. 27. The IMpact-RSV Study Group. Palivizumab, a humanized respiratory syncytial virus monoclonal antibody, reduces hospitalization from respiratory syncytial virus infection in high-risk infants. Pediatrics. 1998;102:531–7. pmid:9738173.
  28. 28. Japan Pediatric Society. Consensus guidelines for the use of palivizumab in Japan. Available from: http://www.jpeds.or.jp/uploads/files/20190402palivizumabGL.pdf. (accessed Jan 8, 2021).
  29. 29. Lee JT, Chang LY, Wang LC, Kao CL, Shao PL, Lu CY, et al. Epidemiology of respiratory syncytial virus infection innorthern Taiwan, 2001–2005—seasonality, clinicalcharacteristics, and disease burden. J Microbiol Immunol Infect. 2007;40(4):293–301. pmid:17712463.
  30. 30. Yu J, Liu C, Xiao Y, Xiang Z, Zhou H, Chen L, et al. Respiratory Syncytial Virus Seasonality, Beijing, China, 2007–2015. Emerg Infect Dis. 2019;25(6):1127–35. pmid:31107230.
  31. 31. Bont L. Nosocomial RSV infection control and outbreak management. Paediatr Respir Rev. 2009;10:16–17. pmid:19651394.
  32. 32. Shek LPC, Lee BW. Epidemiology and seasonality of respiratory tract virus infections in the tropics. Paediatr Respir Rev. 2003;4(2):105–11. pmid:12758047.
  33. 33. Hsu CH, Lin CY, Chi H, Chang JH, Hung HY, Kao HA, et al. Prolonged seasonality of respiratory syncytial virus infection among preterm infants in a subtropical climate. PLoS One. 2014;9(10):e110166. pmid:25333346.
  34. 34. Di Giallonardo F, Kok J, Fernandez M, Carter I, Geoghegan JL, Dwyer DE, et al. Evolution of Human Respiratory Syncytial Virus (RSV) over Multiple Seasons in New South Wales, Australia. Viruses. 2018;10(9):476. pmid:30200580.
  35. 35. Hogan AB, Anderssen RS, Davis S, Moore HC, Lim FJ, Fathima P, et al. Time series analysis of RSV and bronchiolitis seasonality in temperate and tropical Western Australia. Epidemics. 2016;16:49–55. pmid:27294794.
  36. 36. Kim CK, Choi J, Callaway Z, Kim HB, Chung JY, Koh YY, et al. Clinical and epidemiological comparison of human metapneumovirus and respiratory syncytial virus in seoul, Korea, 2003–2008. J Korean Med Sci. 2010;25(3):342–7. pmid:20191030.
  37. 37. Ge X, Guo Y, Chen J, Hu R, Feng X. Epidemiology and Seasonality of Respiratory Viruses Detected from Children with Respiratory Tract Infections in Wuxi, East China. Med Sci Monit. 2018;24:1856–62. pmid:29599424.
  38. 38. Liu P, Xu M, He L, Su L, Wang A, Fu P, et al. Epidemiology of Respiratory Pathogens in Children with Lower Respiratory Tract Infections in Shanghai, China, from 2013 to 2015. Jpn J Infect Dis. 2018;71(1):39–44. pmid:29279451.
  39. 39. Kwon YS, Park SH, Kim MA, Kim HJ, Park JS, Lee MY, et al. Risk of mortality associated with respiratory syncytial virus and influenza infection in adults. BMC Infect Dis. 2017;17(1):785. pmid:29262784.
  40. 40. Britton PN, Hu N, Saravanos G, Shrapnel J, Davis J, Snelling T, et al. COVID-19 public health measures and respiratory syncytial virus. The Lancet Child & Adolescent Health. 2020;4(11):e42–e3. pmid:32956616.
  41. 41. Baker RE, Park SW, Yang W, Vecchi GA, Metcalf CJE, Grenfell BT. The impact of COVID-19 nonpharmaceutical interventions on the future dynamics of endemic infections. Proc Natl Acad Sci U S A. 2020;117(48):30547–53. pmid:33168723.
  42. 42. Nolen LD, Seeman S, Bruden D, Klejka J, Desnoyers C, Tiesinga J, et al. Impact of Social Distancing and Travel Restrictions on non-COVID-19 Respiratory Hospital Admissions in Young Children in Rural Alaska. Clin Infect Dis. 2020 Sep 5;ciaa1328. pmid:32888007.
  43. 43. Yeoh DK, Foley DA, Minney-Smith CA, Martin AC, Mace AO, Sikazwe CT, et al. The impact of COVID-19 public health measures on detections of influenza and respiratory syncytial virus in children during the 2020 Australian winter. Clin Infect Dis. 2020 Sep 28;ciaa1475. pmid:32986804.
  44. 44. Trenholme A, Webb R, Lawrence S, Arrol S, Taylor S, Ameratunga S, et al. COVID-19 and Infant Hospitalizations for Seasonal Respiratory Virus Infections, New Zealand, 2020. Emerg Infect Dis. 2021;27(2). pmid:33263515.
  45. 45. Kuitunen I, Artama M, Makela L, Backman K, Heiskanen-Kosma T, Renko M. Effect of Social Distancing Due to the COVID-19 Pandemic on the Incidence of Viral Respiratory Tract Infections in Children in Finland During Early 2020. Pediatr Infect Dis J. 2020;39(12):e423–e7. pmid:32773660.