Next Article in Journal
EBV-Associated Hub Genes as Potential Biomarkers for Predicting the Prognosis of Nasopharyngeal Carcinoma
Next Article in Special Issue
Influenza A and B Viruses in Fine Aerosols of Exhaled Breath Samples from Patients in Tropical Singapore
Previous Article in Journal
Hepatitis B Virus Epsilon (ε) RNA Element: Dynamic Regulator of Viral Replication and Attractive Therapeutic Target
Previous Article in Special Issue
The Experience of Testing for Coronavirus Disease (COVID-19) at a Single Diagnostic Center in Paraguay before the Introduction of Vaccination
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear and Multidelayed Effects of Meteorological Drivers on Human Respiratory Syncytial Virus Infection in Japan

1
Division of International Health (Public Health), Graduate School of Medical and Dental Sciences, Niigata University, Niigata 951-8510, Japan
2
Japan Society for the Promotion of Science, Tokyo 102-0083, Japan
3
College of Health and Medicine, School of Medicine, University of Tasmania, Hobart 7000, Australia
*
Author to whom correspondence should be addressed.
Viruses 2023, 15(9), 1914; https://doi.org/10.3390/v15091914
Submission received: 9 August 2023 / Revised: 7 September 2023 / Accepted: 11 September 2023 / Published: 12 September 2023
(This article belongs to the Special Issue Animal and Human Respiratory Viruses—Causes of the Next Pandemic?)

Abstract

:
In this study, we aimed to characterize the nonlinear and multidelayed effects of multiple meteorological drivers on human respiratory syncytial virus (HRSV) infection epidemics in Japan. The prefecture-specific weekly time-series of the number of newly confirmed HRSV infection cases and multiple meteorological variables were collected for 47 Japanese prefectures from 1 January 2014 to 31 December 2019. We combined standard time-series generalized linear models with distributed lag nonlinear models to determine the exposure–lag–response association between the incidence relative risks (IRRs) of HRSV infection and its meteorological drivers. Pooling the 2-week cumulative estimates showed that overall high ambient temperatures (22.7 °C at the 75th percentile compared to 16.3 °C) and high relative humidity (76.4% at the 75th percentile compared to 70.4%) were associated with higher HRSV infection incidence (IRR for ambient temperature 1.068, 95% confidence interval [CI], 1.056–1.079; IRR for relative humidity 1.045, 95% CI, 1.032–1.059). Precipitation revealed a positive association trend, and for wind speed, clear evidence of a negative association was found. Our findings provide a basic picture of the seasonality of HRSV transmission and its nonlinear association with multiple meteorological drivers in the pre-HRSV-vaccination and pre-coronavirus disease 2019 (COVID-19) era in Japan.

1. Introduction

The human respiratory syncytial virus (HRSV) infects the human respiratory tract, causing clinically severe pneumonia in young children and bronchitis in infants [1]. Globally, HRSV-induced acute lower respiratory tract infections (ALRIs) are responsible for approximately 70,000 deaths in children under the age of 5 years and hospitalization of approximately 3.4 million people [2,3]. The magnitude of the global disease burden caused by HRSV has recently been recognized to affect not only infants and young children but also the elderly (≥65 years), and it remains a significant public health concern [4,5]. The transmission dynamics of HRSV infection epidemics are seasonally driven in temperate regions, with most infections occurring during annual autumn–winter (September–February) seasons; however, a seasonal shift in the HRSV infection epidemics to the summer–autumn (June–November) seasons between 2017 and 2019 has been reported in Japan [6,7,8,9]. For instance, epidemics of HRSV infection in the United States and 15 European countries, where the climatic environment is like that of Japan, are also concentrated in autumn–winter seasons [6,10]. It should be noted that off-season HRSV infection epidemics were observed during the coronavirus disease 2019 (COVID-19) pandemic in Japan (i.e., no major epidemic in 2020, and an unusually high number of cases reported in 2021) [11,12].
The transmission dynamics of HRSV are shaped by a web of biological and nonbiological factors, and it has been suggested that in addition to intrinsic transmission factors such as innate immunity, viral survival, and antigenic variation, extrinsic seasonal climate-driven factors influence the transmission dynamics [13,14]. Although previous studies have reported partially conflicting results, several phenomenological trends in the relative role of multiple meteorological conditions (e.g., mean ambient temperature, relative humidity, precipitation, and wind speed) in HRSV transmission have emerged. For instance, two systematic reviews concluded that high ambient temperature and relative humidity were associated with HRSV transmission in tropical and subtropical regions, while in temperate regions, low ambient temperature and high relative humidity were associated with the occurrence of several HRSV-caused epidemics [15,16]. In other epidemiological studies, ambient temperature has been found to have a negative association with HRSV infection incidence in European countries [17,18], the United States [19,20], and Mainland China [21,22,23,24]. Meanwhile, relative humidity has been reported to have a positive association with HRSV infection incidence in European countries and Hong Kong [17,25,26]. Overall, few epidemiological studies have described the relationship between precipitation and wind speed and HRSV infection incidence; however, an association between peak precipitation and HRSV infection incidence has been observed in several tropical regions [27,28,29,30,31], and a moderate negative correlation between wind speed and its incidence has been observed in Kenya [32]. Indeed, the differences in these reports could partially be explained by the variation in the spatial scales of analysis, application of different statistical methods, and consideration of various potential confounders. Epidemiologic studies assessing the relationship between climate variability and transmission could have methodologic limitations that may introduce bias and limit causal inferences [33,34,35].
Although previous studies have gradually revealed the associations between the climate variability and the transmission dynamics of HRSV, the effect of meteorological drivers on the timing and intensity of its epidemics is poorly understood. Generally, the effects of climate variability and disease transmission are complex, nonlinear, and often shape multiple, delayed events, limiting statistical inferences based on traditional deterministic linear modelling and single-delayed effects frameworks [36]. While these traditional methods aid in the forecasting of HRSV infection epidemics, current modelling frameworks may oversimplify the nature of relationships between meteorological drivers and disease dynamics. The time-varying distributed lag nonlinear models (DLNMs) developed in the last decade have been relatively underused, despite being well suited for time-series statistical modelling of infectious disease dynamics [37,38]. In recent years, these models have begun to be adapted to climate-sensitive vector-borne diseases (e.g., dengue fever, malaria, and West Nile virus) [39,40,41,42,43,44], but their application for respiratory virus infections, including HRSV infections, has been quite limited. In the present study, we aimed to explore the nonlinear and multidelayed effects of multiple meteorological drivers (i.e., mean ambient temperature, relative humidity, precipitation, and wind speed) on the seasonality of HRSV infections using weekly surveillance data across 47 prefectures in Japan over the 2014–2019 period by designing a cutting-edge time-series statistical model. These exposure–lag–response associations can reveal how meteorological drivers affect HRSV infection risk in the weeks leading up to an outbreak. Furthermore, if the explanatory powers of these meteorological drivers can be determined, the seasonality and transmission of these viruses in different environments can be better understood, which would contribute to the development of early warning systems to respond in a timely manner to annual HRSV epidemic surges and public health strategies and interventions. By characterizing the relative role of these meteorological drivers on HRSV infection dynamics, we attempted to gain an insight into the seasonal mechanisms of the disease, thereby rendering basic information to disentangle its extrinsic environmental drivers.

2. Materials and Methods

2.1. Design Setting

This time-series statistical modeling analysis investigated the potential nonlinear and multilagged relationships of HRSV transmission with multiple meteorological drivers (i.e., mean ambient temperature, relative humidity, rainfall, and wind speed) in Japan, using an ecological study design. In this study, we systematically collected and analyzed epidemiological data for the number of weekly newly confirmed HRSV infection cases and meteorological drivers (i.e., mean ambient temperature, relative humidity, rainfall, and wind speed) from 1 January 2014 to 31 December 2019 (between the 1st week of 2014 and 52nd week of 2019) across 47 prefectures in Japan, covering the entire country. We selected these locations based on the geographical diversity and availability of data across the study period. It should be noted that the transmission dynamics of HRSV have drastically changed since 2020 owing to the COVID-19 pandemic; hence, this period was not included in the analysis [11,12]. The post-COVID-19 pandemic period will be evaluated in a future project. Japan is located at latitudes of approximately 26–43° N and longitudes of approximately 127–141° E in the Western Pacific Region and comprises a total of 47 Japanese prefectures (Figure S1). Given their geoclimatic characteristics, 46 prefectures, excluding Okinawa, are classified as temperate regions, and Okinawa as a subtropical region with four distinct seasons in Japan; summers are hot and humid (June–August), while winters are cold and dry (December–February).

2.2. Empirical Datasets

2.2.1. National HRSV Surveillance Data

The HRSV infection epidemiological data used in this study were obtained from the Infectious Disease Weekly Report (IDWR), sourced from the National Epidemiological Surveillance of Infectious Diseases (NESID) data published by the National Institute of Infectious Diseases, Japan (NIID), under the authority of the Ministry of Health, Labor and Welfare, Japan (MHLW) [45]. The MHLW manages approximately 3000 pediatric sentinel sites (i.e., hospitals and clinics) in Japan and reports the number of patients diagnosed with HRSV infection on a weekly basis to the prefecture or municipal public health sectors in Japan [46,47]. A confirmed case of HRSV infection is defined as a positive result in a rapid diagnostic test (RDT) using a test kit licensed in Japan or a laboratory confirmation such as virus isolation or antibody titer increase in paired sera according to the MHLW guidelines [48]. The number of sentinels assigned to each public health service area is based on population size: a public health center with <30,000 individuals is assigned one sentinel, a center with 30,000–75,000 individuals, two sentinels, and that with >75,000 individuals, three or more sentinels, as determined by the following formula: 3  +  (population − 75,000)/50,000 [9]. These sentinel sites forward clinical data to approximately 60 prefectural or municipal public health sectors, and the data are electronically reported to the NIID; the number of HRSV infection cases is released weekly on the NIID website. In the present study, we extracted the number of weekly newly confirmed HRSV infection cases across all 47 prefectures in Japan reported in weeks 1–52 from 2014 to 2019, from the NESID database.

2.2.2. Meteorological Data

The Automated Meteorological Data Acquisition System (AMeDAS) developed by the Japan Meteorological Agency is a high-resolution surface observation network for investigating meteorological conditions in Japan. In the present study, we retrieved prefecture-specific daily time-series of meteorological data across 47 prefectures in Japan collected at the AMeDAS over the study period [49]. The weekly meteorological data, including mean ambient temperature in degrees Celsius (°C), relative humidity in percent (%), precipitation in millimeters (mm), and wind speed in meters per second (m/s), published by the website were calculated from the daily records and utilized as independent variables in the time-series statistical models presented here. These meteorological variables were chosen according to their availability and use in previous HRSV modelling studies [8,11,17,19,50,51,52]. Meteorological data collected from meteorological observatories (meteorological stations) situated in the prefectural capitals were utilized for each prefecture.

2.3. Statistical Analysis

2.3.1. Descriptive Statistics

To determine the key characteristics of the multitudes of time-series datasets used in this study, we described the probability distributions of the number of weekly newly confirmed HRSV infection cases and meteorological drivers (i.e., mean ambient temperature, relative humidity, precipitation, and wind speed) across all 47 prefectures in Japan during the study period by utilizing the following descriptive statistics: mean, standard deviation (SD), minimum (Min), 25th percentile (P25), 50th percentile (P50), 75th percentile (P75), and maximum (Max).

2.3.2. Construction of the Time-Series Statistical Model

To establish a robust and reliable time-series statistical model, multiple stages were incorporated into it [11,53,54]. Prior to constructing the model, we checked the probability distribution of the dependent variables and number of weekly newly confirmed HRSV infection cases (the normality of probability distribution was assessed by the Shapiro–Wilk test) (Figure S2), followed by an assessment of the relationships (e.g., linearity) between number of weekly newly confirmed HRSV infection cases and each independent variable. Generally, by assessing the transformation of dependent variables in time-series-driven systems, these variables can be used as stationary effects, often improving forecasting accuracy. All dependent variables and independent variables included in the statistical models were assessed for multicollinearity using pairwise Spearman’s rank-order cross-correlation coefficient (ρ). If the variables were found to be strongly linearly correlated (cut-off of |ρ| > 0.8), the variable with the largest mean absolute statistical correlation with the other independent variables was eliminated [55]. In the preliminary analysis, no independent variables showing strong statistical linear correlations were observed (Table S1).
In the present study, we fitted standard time-series multivariate generalized linear models (GLMs) with a quasi-Poisson distribution family and logarithmic-link function, allowing for overdispersion in the observational data, by adding the time-varying DLNMs to simultaneously describe nonlinear and multidelayed dependencies between the incidence relative risks (IRRs) of HRSV infection and meteorological drivers. Briefly, a DLNM was designed for each targeted meteorological driver and added to the regression equation of GLM. This method was used as a primary model to disentangle the underlying complex association between the number of weekly newly confirmed HRSV infection cases and different multiple meteorological drivers (i.e., mean ambient temperature, relative humidity, precipitation, and wind speed) as main exposures. In a DLNM, a cross-basis function that combines two tensor product of basis-functions is introduced: one represents the dependent variable’s probability distribution in the independent variable and lagging dimensions to simultaneously assess the lag effect and the other represents the nonlinear effects of the exposure driver [37,38,56].
The cross-basis term for time-varying DLNMs acts as basis function predictor in two dimensions: exposure and lag spaces. Specifically, we modeled exposure–response associations using a natural cubic spline with three degrees of freedom (df) and modeled the lag–response association using a natural cubic spline with three df. To adjust for the possible harvesting and misalignment from the delayed heterogeneity of meteorological drivers of HRSV infection dynamics, based on multiple accumulated studies investigating the incubation period (i.e., approximately 2–8 days) and reporting delays of HRSV [8,17,19,50,51,52,57,58], we considered temporal lags (i.e., delays in potential effect) of up to 2 weeks as the default lag structure for the cross-basis function of each independent variable related to the dependent variable. These temporal time lags are useful for describing the biologically and physically plausible time lags in the population dynamics, natural history of the HRSV host reservoir, and subsequent incubation period before disease notification.
The general algebraic definition of the time-series statistical models is formulated as follows:
log [ ( y i , t | r i , t ) ] : = α 0 + s β s f g x i , t , s ; d f + u h g z i , t , u ; θ + e t + p t + s t + log r i , t 1 + log O [ t ] + ε [ i , t ]
where y [ i , t ] is the outcome time-series; r [ i , t ] is the expected time-series of the number of weekly newly confirmed HRSV infection cases in prefecture i on week t; the term α 0 corresponds to the overall intercept; s β s f g x i , t , s ; d f denotes the cross-basis function of DLNMs with exposure and multilagged effects modelled by a natural cubic spline function and a linear function of multiple meteorological drivers (i.e., mean ambient temperature, relative humidity, precipitation, and wind speed) in prefecture i in week t, respectively. We also modeled baseline risk along with shared long-term seasonal variations and cycles and short-term trends by incorporating natural cubic splines of time (7 degrees of freedom (df) per year) as term u h g z i , t , u ; θ , year as term e [ t ] , and number of public holidays per week as term p [ t ] , with the fixed effects variables as possible confounders [11,53,54,59,60,61]. s [ t ] denotes prefectural characteristics or regional variable indicators in prefecture i. Furthermore, the autocorrelation of residuals in the case of infectious disease was pathogen-specific and needed to be accounted for; therefore, autoregressive terms at order one (as term r [ i , t 1 ] ) was incorporated into the statistical models, accounting for potential serial correlation [53,54,59,60,61]. The teem log O [ t ] denotes the logarithm of the yearly population (per 100,000) by prefecture as the offset term [62]. The use of these population estimates for our denominator in the models allowed for the exponentiation of each coefficient to be expressed as an IRR, which is an intuitive representation of the association of variables with the increased risk of HRSV infection incidence. The term ε [ i , t ] indicate errors. To quantify the total contribution, independent effects, and relative importance of multiple meteorological drivers, we included all variables in the same model [53,63]. By including all variables of interest in the same regression equation, we were able to strengthen the interpretation of the effects as independent or additive, based on accumulated empirical knowledge.
To test the sensitivity of the results to the modeling choices described above, we repeated the analysis by varying the df of the natural cubic spline of time from 7 df per year to 3 df or 11 df per year. We also performed a sensitivity analysis of the observed effect on the weeks of lags accounted for in the statistical model by modifying the length of the lag period from 2–3 and 4 weeks. We quantified the estimates as cumulative-week (i.e., 0–2 weeks) and single-week (i.e., 0, 1, and 2 weeks) IRRs together with the 95% confidence intervals (CIs) at the 25th and 75th percentiles of all meteorological drivers to determine the strengths of the associations. The reference value for ambient temperature and relative humidity were set as the median (50th percentile), while the precipitation and wind speed were set as 0.0 mm and 0.0 m/s, respectively. Additionally, we plotted the exposure–response curves to describe the overall cumulative associations. Statistical significance was considered at a p-value of <0.05 (i.e., type I error), on a two-tailed test. All analyses were performed using STATA version 15.1 statistical software (Stata Corp, College Station, TX, USA) and R statistical programming software version 4.1.0 (R Foundation for Statistical Computing, Vienna, Austria) using the “dlnm” [64].

2.4. Ethical Considerations

The present ecological modeling study analyzed publicly available data in Japan. As such, the epidemiological datasets utilized in this study were de-identified and fully anonymized in advance, and the analysis of publicly available data with no identifying information did not require ethical approval. The present study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).

3. Results

3.1. Descriptive Analysis

Overall, 721,709 newly confirmed HRSV infection cases were observed during the 313-week study period across all 47 prefectures in Japan (Figure S1). The mean number of the weekly newly confirmed cases, i.e., the average for all prefectures and weeks during the study period, was 49 (range, 0–892) (Table 1). The weekly mean ambient temperature, relative humidity, precipitation, and wind speed also had wide ranges (−6.7–32.2 °C, 30.8–97.2%, 0.0–94.1 mm, and 0.9–11.3 m/s, respectively).

3.2. Assessing the Nonlinear Effects of Meteorological Drivers on HRSV Infection Incidence

The pooled overall cumulative relationships between multiple meteorological drivers (i.e., mean ambient temperature, relative humidity, precipitation, and wind speed) and HRSV infection incidence are illustrated in Figure 1. Notably, we observed a nonlinear relationship between mean ambient temperature and HRSV infection incidence (Figure 1A). In particular, for the range of approximately 10.0–15.0 °C, we found a linear inverse mean ambient temperature–incidence relationship, with low mean ambient temperatures associated with increased HSRV infection incidence. In contrast, the cumulative IRR increased when the mean ambient temperature was dramatically higher than approximately 20.0 °C. Specifically, the corresponding cumulative IRRs were 1.023 (95% CI, 1.008–1.039) at 8.1 °C (25th percentile), and 1.068 (95% CI, 1.056–1.079) at 22.7 °C (75th percentile), respectively, with reference to the IRR at 16.3 °C (Table 1). A J-shaped association was observed between relative humidity and HRSV infection incidence (Figure 1B). Relative humidity exhibited a positive and almost linear relationship with HRSV infection incidence above approximately 70.0%. The 2-week cumulative IRRs were 0.997 (95% CI, 0.971–1.003) at a relative humidity of 63.1% (25th percentile) and 1.045 (95% CI, 1.032–1.059) at a relative humidity of 76.4% (75th percentile) with reference to the IRR at 70.4% (Table 1). The overall relationship between precipitation and HRSV infection incidence was weak (Figure 1C). An almost nonlinear positive association trend with incidence was observed across the entire precipitation range. For instance, the 2-week cumulative IRRs were 1.027 (95% CI, 1.015–1.038) at a wind speed of 0.7 mm (25th percentile) and 1.064 (95% CI, 1.032–1.097) at 6.2 mm (75th percentile) with reference to that at 0.0 mm (Table 1). A clear inverse association was observed between wind speed and HRSV infection incidence (Figure 1D). Specifically, the cumulative IRR decreased rapidly with increasing wind speed across the entire range. The 2-week cumulative IRRs were 0.793 (95% CI, 0.697–0.903) at a wind speed of 2.2 m/s (25th percentile) and 0.775 (95% CI, 0.678–0.885) at 3.4 m/s (75th percentile), with reference to that at 0.0 m/s (Table 1).

3.3. Sensitivity Analysis

The main findings described above were confirmed by repeating the series of sensitivity analyses utilizing alternative DLNM specifications. For all meteorological drivers (i.e., mean ambient temperature, relative humidity, precipitation, and wind speed), changing the incorporation of the df of the natural cubic spline of time from 7 df per year to 3 df (Figure S3) or 11 df per year (Figure S4) revealed that the observed risk effect shape was substantially robust over the different parameterizations. Generally, although we observed wider 95% CIs after slightly increasing the df for each of the meteorological drivers, the estimated shapes of the exposure-response functions remained largely similar. Further, for all meteorological drivers, on changing the incorporation from 0–2 lag weeks to 0–3 lag weeks (Figure S5) and 0–4 lag weeks (Figure S6), the observed shapes of the incidence curves were found to be consistent. Generally, model uncertainties tend to increase when more weeks are included. Taken together, the series of sensitivity analyses confirmed the robustness of the findings of the main analysis.

4. Discussion

In the present study, we retrospectively modelled the nonlinear and multidelayed temporal effects of multiple meteorological drivers (i.e., mean ambient temperature, relative humidity, precipitation, and wind speed) on HRSV seasonality using multiseason data across all 47 prefectures in Japan. To our knowledge, this study presents one of the most comprehensive assessments of the effects of meteorological drivers on HRSV infection seasonality across a large gradient of climate conditions related to latitudes and longitudes. Overall, we found evidence of a modest nonmonotonic association between these meteorological drivers and HRSV infection incidence over lags of 0–2 weeks, contributing to approximately 14.8% of the variation in IRR over the study period. More specifically, an increase in mean ambient temperature from 8.1 °C (25th percentile) to 22.7 °C (75th percentile) was associated with an increase of approximately 4.5% in the cumulative risk of HRSV infection incidence over a 2-week period. The other meteorological drivers (relative humidity, precipitation, and wind speed) explained approximately 4.8%, 3.7%, and 1.8% of the variation, respectively. Despite these simplified assumptions, our findings objectively reveal that the temporal driving patterns of HRSV infection could be partially explained by these meteorological drivers. In addition to highlighting the relative role played by climate variability, our findings indicate that there may be an additional underlying mechanism involved in shaping the time-dependent transmission heterogenicity of HRSV that remains to be identified. By identifying causal relationships between meteorological drivers and HRSV infection epidemic patterns, we can build predictive models of the epidemic that would help inform guidelines for timing of monoprophylaxis, such as palivizumab, and HRSV vaccination when it becomes available [65,66].
Our results are fairly consistent with the results of previous studies that investigated the broad mechanistic principles underlying the association between meteorological factors and HRSV infection seasonality. For instance, the observation that moderately low ambient temperatures lead to high viral transmission rates has been made in many previous studies. In particular, the negative association between mean ambient temperature and HRSV infection incidence, which was observed in this present study, is biologically valid and consistent with the results of previous studies conducted in mainland China and Brazil [23,67]. Moreover, laboratory studies and animal experiments have described that low ambient temperature plays a potential role in modulating the viability and stability of respiratory viruses by affecting the properties of viral surface proteins and lipid membrane and proportion of droplet nuclei [68,69,70,71]. Lower ambient temperatures may also enhance virus susceptibility by triggering changes in human physiology. In contrast, a linear increase in the cumulative IRR as the mean ambient temperature increased to above approximately 20.0 °C was observed in this study, similar to the form of nonlinear association identified in two studies that were recently conducted in Singapore [50,51]. For instance, HRSV infection epidemics in tropical and subtropical regions (e.g., Hong Kong, Singapore, Malaysia, and Colombia) have been associated with high ambient temperatures [72]. Recently, it has also been observed that the HRSV infection seasonality has changed from autumn–winter (September–February) to summer–autumn (June–November) in Japan after 2017, which is generally a temperate region [6,7,8,9]. A notable study using national Japanese surveillance and meteorological data examined the causes of summer–autumn epidemic conditions and found that the interaction between ambient temperature (≥28.2 °C) and relative humidity (≥79.0%) had a 1.92-fold (95% CI, 1.60–2.23) marginal effect on HRSV infection incidence [9]. Another study by Yusuf et al. that examined HRSV infection incidence throughout a year in nine cities that had different climates and geographical locations suggested that the incidence was related to ambient temperature in a bimodal fashion, with greater disease incidences at ambient temperatures of 24.0–30.0 °C and 2.0–6.0 °C [73]. These accumulated epidemiological findings suggest that HRSVs are partially active at relatively high and low ambient temperatures. Although the underlying cause for this complex form of association is unclear, our study, which captured data over multiple seasons across climate zones in Japan, found that ambient temperature had a U-shaped association with incidence. Our findings provide an important perspective on how the variation in ambient temperature between regions may affect HRSV infection epidemics in Japan.
There is ample evidence of a positive association between humidity (including relative and/or absolute) and HRSV infection epidemics, generally consistent with our results. Studies conducted in several regions, including Mexico, Spain, and Italy, have shown an association between high HRSV infection incidence with high humidity [18,25,73]. Indeed, our results are consistent with the modeling results of a study conducted in Singapore that found an inverse J-shaped association between humidity and HRSV infection epidemics [51]. One potential mechanism by which high humidity may facilitate HRSV transmission is through increased viral survival; as HRSV is transmitted via droplets and respiratory secretions, it tends to be rapidly inactivated in small aerosols at low humidity levels and higher HRSV stability in large-particle aerosols at higher humidity levels [18,51,73,74].
To our knowledge, little is known about the effects of precipitation and wind speed on HRSV dynamics; however, our findings show that high levels of precipitation and low levels of wind speeds are weakly associated with increased HRSV infection risk in Japan. Previously, several studies in Italy, Singapore, Philippine, Thailand, and Brazil have already reported that high precipitation is associated with increased HRSV infection incidence, suggesting that the epidemic seasonality may show a similar pattern to precipitation [51,52,67,75,76]. Indeed, the peak in HRSV infection epidemics coincides with the rainy season in several tropical locations [27,28,29,30,31]. In contrast, limited studies in Singapore and Kenya have found wind speed to be negatively associated with the HRSV infection incidence [32,51]. Generally, some of the previous literature suggests that low wind speeds can prolong the presence of infected aerosol particles in the air, thereby facilitating the transmission of respiratory viruses [77]. Also, higher wind speeds can disperse droplets and decrease the concentration of infected aerosol particles [77,78]. Overall, these epidemiological findings were partially consistent with the trends in our results; however, further studies are needed to assess the relative role of these drivers in the HRSV transmission.
Overall, a possible another interpretation for our observations is that ambient meteorological drivers are associated with the amount of time spent indoors and outdoors, thereby facilitating transmission. Time spent indoors in closed or confined environments could theoretically facilitate the transmission of respiratory viral infections such as HRSV. Indoor congestion has often been hypothesized as a primary contributor to the surge in HRSV infection cases during colder seasons in temperate regions [79,80]. Conversely, in subtropical regions, the prevalence of hot and humid climatic conditions during summer drives individuals to seek shelter in densely populated, air-conditioned settings, thereby creating a favorable environment for viral transmission [81]. Notably, most meteorological conditions vary between indoor and outdoor environments and the form of associations also differs (e.g., influenza and COVID-19) [33,82,83]. For instance, these meteorological drivers may be a proxy for another variable that affects transmission (e.g., human behavior and living environment), or there may be a complex relationship between climate variability and disease transmission (e.g., nonlinear or threshold effects). Indeed, these ambiguities may render complex nonlinear associations more plausible. Despite the regional proximity and similar climates, the epidemiological patterns of HRSV infection in the Western Pacific and South-East Asia regions, including Japan, appear to exhibit subtle differences, suggesting that climatic factors alone do not entirely dictate epidemic seasonality. Indeed, one of our significant concerns is the lack of studies on the association between ambient meteorological drivers, human behavior, and HRSV infection epidemics. Thus, detailed, multicountry, multicity studies with different epidemiological contexts and long-term studies are therefore needed to further disentangle complex interactions and present robust evidence.
There are several caveats when interpreting the results of this study. First, our study was based on the analysis of secondary data and is regarded as a type of ecological study in causal inference, which is prone to unobserved confounding [8,11,53,54]. In particular, we did not have access to data on the different genetic strains of HRSV that emerged over time in the different prefectures; therefore, we did not explicitly model the effects of HRSV subtypes (A and B), which could have helped to explain the differences in the HRSV transmission seasonality in Japan [84]. Furthermore, increased contact structures (i.e., contact patterns and rates) among older students during school terms have been also hypothesized to play an important role in HRSV seasonality [85,86]. These factors influence the accurate modeling of the true extent of disease infections and transmission within populations; however, disease surveillance data allow researchers to determine temporal trends in disease dynamics that can be further modeled and used in public health decision making. In particular, the influence of individual exposure on infection risk and causality at the individual level cannot be inferred over a large geographical area of Japan. Second, our time-series statistical model only accounted for meteorological drivers; other important drivers such as human behavioral changes could also drive viral transmission, although some of these were associated with the selected meteorological drivers (e.g., indoor crowding in hot/cold weather) [24,79,80,81]. Since the beginning of the COVID-19 pandemic, a near cessation of HRSV transmission during periods of movement restrictions and physical distancing has been observed. However, on/off-season HRSV infection epidemics have been reported in multiple countries, which may be attributed to the relaxation of nonpharmaceutical interventions and an increase in population susceptibility [87,88]. As all our data were collected prior to the COVID-19 pandemic, our current modeling framework was not devised to incorporate these factors. Fourth, our study may be limited in its generalizability due to the collection of data from a specific geographical area (47 Japanese prefectures) and a limited period (2014–2019), thus precluding the representation of different epidemic periods across diverse regions. For similar future studies, it may be beneficial to consider increasing the number of geographical regions, expanding the set of indicators, and including additional data, such as individual-level data. Fifth, although we only utilized meteorological drivers such as mean ambient temperature, relative humidity, precipitation, and wind speed as explanatory variables in our time-series statistical models, other environmental confounders such as diurnal temperature ranges, sunshine hours, and ultraviolet levels may also contribute to HRSV transmission [50,51,89,90,91]. In particular, some air pollutants (e.g., particulate matter of aerodynamic diameter of less than 2.5 µm, nitrogen dioxide, sulfur dioxide, and carbon monoxide) have been widely considered in epidemiological studies of HRSV and further modelling studies are needed to quantify the contribution of these potential drivers to the transmission dynamics of HRSV across Japan and to understand the interrelationships [50]. Seventh, because our study focused on pooled associations between meteorological drivers and seasonal variations in incidence, we did not fully characterize the association between the heterogeneous transmission of HRSV in specific Japanese prefectures, which was beyond the scope of this study. For instance, the occurrence in local HRSV infection epidemics in remote prefectures of Japan should be considered, as the epidemiological effect on HRSV transmission may differ among regions [11,46]. To identify the actual causes of the heterogeneity among the prefectures, we need to validate these causes by formulating sophisticated spatiotemporal hierarchical modeling frameworks (e.g., approximate Bayesian inference using integrated nested Laplace approximation) and DLNMs [92] while simultaneously accounting for spatial heterogeneity and autocorrelation or a two-stage time-series design involving a downscaling procedure [37,56]. Moreover, the application of techniques such as analysis of partial differential equation systems and numerical simulation may also allow for more flexible modelling of exposure–lag–response relationships [93,94,95,96]. Indeed, this our method does not provide site-specific information on meteorological relationships with HRSV infection incidence. Therefore, our present results reflect an overall average trend during the study period in Japan that should be interpreted with caution. Finally, although this study successfully described the space-varying, nonlinear, and multidelayed associations between the meteorological drivers and HRSV infection seasonality across Japan, there remain some future studies that could expand our current knowledge. Particularly, the extent to which the future long-term risk of HRSV infection may change under different current and future climate change scenarios has not been predicted and is not sufficiently clear. For instance, harmonized input data and common simulation protocols from the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP) can provide useful insights into the possible consequences of several key climate change options open to us (e.g., dengue, malaria, and vibriosis) [97].

5. Conclusions

Disease transmission is a dynamic interplay involving a network of contributing factors. Here, using empirical epidemiological data, we characterized the nonlinear and multidelayed temporal associations between the meteorological drivers and HRSV infection seasonality in the pre-HRSV-vaccination and pre-COVID-19 era across Japan. As a result, our findings indicate that hot and humid weather, high levels of precipitation, and low levels of wind speeds are associated with increased HRSV infection risk in Japan. Note that most HRSV and other seasonal respiratory viruses actually transmit indoors, rather than outdoors; therefore, the relevance of these outdoor climate correlations maybe most useful as a surrogate for indoor environments at the time of such HRSV incidence and transmission events. Our findings align, to an extent, with available mechanistic explanations and reports from the previous literature from other regions of the world. The meteorological drivers investigated in our method may help researchers to partially explain differences in the strength of seasonality across the different regions of Japan, in contrast with previous methods that assume a linear relationship of seasonality with disease risk and changing exposures. More broadly, this study may suggest to policymakers the need for public health strategies and interventions that are flexibly adapted to the climatic conditions in different regions of Japan to mitigate HRSV transmission (e.g., increase the intensity and frequency of interventions in periods of hot and humid weather, high levels of precipitation, and low levels of wind speeds). More importantly, substantively, with several prophylactic vaccine candidates (e.g., prefusion-stabilized F protein vaccine) on the horizon, our present findings on the HRSV infection seasonality in Japan could help to predict the onset of the HRSV season and may inform the timing of immunization as well as seasonal health service planning by pediatric clinical and laboratory diagnostic teams to prevent severe outcomes in high-risk groups [65,66]. Future analyses should consider how the associations addressed in this study can inform the forecasting of the spatiotemporal distribution and seasonality of HRSV under anticipated climate regimes, immunization, and COVID-19-related disruptions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15091914/s1, Figure S1: The geographic distribution of the 47 Japanese prefectures and their locations; Figure S2: Probability distribution of weekly newly confirmed HRSV cases across all the included prefectures and days; Figure S3: Assessing the pooled nonlinear association of the IRRs of HRSV incidence with meteorological variables; Figure S4: Assessing the pooled nonlinear association of the IRRs of HRSV incidence with meteorological variables; Figure S5: Assessing the pooled nonlinear association of the IRRs of HRSV incidence with meteorological variables; Figure S6: Assessing the pooled nonlinear association of the IRRs of HRSV incidence with meteorological variables; Table S1: Spearman’s rank-order linear correlation matrix between weekly newly confirmed HRSV cases and meteorological variables.

Author Contributions

Conceptualization, K.W.; methodology, K.W.; software, K.W.; validation, K.W.; formal analysis, K.W.; investigation, K.W.; resources, K.W.; data curation, K.W.; writing—original draft preparation, K.W.; writing—review and editing, I.S.K. and R.S.; visualization, K.W.; supervision, I.S.K. and R.S.; project administration, K.W.; funding acquisition, K.W. and R.S. All authors have read and agreed to the published version of the manuscript.

Funding

K.W. received funding from the Grants-in-Aid for Scientific Research (KAKENHI) of the Japan Society for the Promotion of Science (JSPS) (22J23183), the Community Medical Research Grant of the Niigata City Medical Association (GC03220213), and the Tsukada Medical Research Grant (grant number not available). R.S. received funding from the Japan Initiative for Global Research Network on Infectious Diseases (J-GRID) by the Japan Agency for Medical Research and Development (AMED) (15fm0108009h0001–21wm0125005h0002), KAKENHI by the JSPS (21K10414), the Health and Labor Sciences Research Grants, Ministry of Health, Labor and Welfare, Japan (H30-Shinkougyousei-Shitei-002 and H30-Shinkougyousei-Shitei-004), and the Niigata Prefecture Coronavirus Infectious Disease Control Research and Human Resources Development Support Fund (grant number not available). The funder of the study had no role in the study design, data collection, data analysis, data interpretation, or writing of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

An anonymized dataset that enables replication of the analysis is publicly available and is available from the corresponding author upon request.

Acknowledgments

We also thank all members of the local governments, public health centers and institutes, and the National Institute of Infectious Diseases, Japan, for their significant efforts in surveillance, laboratory testing, epidemiological investigations, and data collection.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Piedimonte, G.; Perez, M.K. Respiratory syncytial virus infection and bronchiolitis. Pediatr. Rev. 2014, 35, 519–530. [Google Scholar] [CrossRef]
  2. Nair, H.; Nokes, D.J.; Gessner, B.D.; Dherani, M.; Madhi, S.A.; Singleton, R.J.; O’Brien, K.L.; Roca, A.; Wright, P.F.; Bruce, N.; 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, 1545–1555. [Google Scholar] [CrossRef] [PubMed]
  3. Hirve, S.; Crawford, N.; Palekar, R.; Zhang, W. Clinical characteristics, predictors, and performance of case definition-Interim results from the WHO global respiratory syncytial virus surveillance pilot. Influenza Other Respir. Viruses 2020, 14, 647–657. [Google Scholar] [CrossRef]
  4. Wang, X.; Li, Y.; Deloria-Knoll, M.; Madhi, S.A.; Cohen, C.; Ali, A.; Basnet, S.; Bassat, Q.; Brooks, W.A.; Chittaganpitch, M.; et al. Global burden of acute lower respiratory infection associated with human metapneumovirus in children under 5 years in 2018: A systematic review and modelling study. Lancet Glob. Health 2021, 9, e33–e43. [Google Scholar] [CrossRef]
  5. Mazur, N.I.; Higgins, D.; Nunes, M.C.; Melero, J.A.; Langedijk, A.C.; Horsley, N.; Buchholz, U.J.; Openshaw, P.J.; McLellan, J.S.; Englund, J.A.; et al. The respiratory syncytial virus vaccine landscape: Lessons from the graveyard and promising candidates. Lancet Infect. Dis. 2018, 18, e295–e311. [Google Scholar] [CrossRef]
  6. Miyama, T.; Iritani, N.; Nishio, T.; Ukai, T.; Satsuki, Y.; Miyata, H.; Shintani, A.; Hiroi, S.; Motomura, K.; Kobayashi, K. Seasonal shift in epidemics of respiratory syncytial virus infection in Japan. Epidemiol. Infect. 2021, 149, e55. [Google Scholar] [CrossRef]
  7. 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. [Google Scholar] [CrossRef]
  8. Wagatsuma, K.; Koolhof, I.S.; Shobugawa, Y.; Saito, R. 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. PLoS ONE 2021, 16, e0248932. [Google Scholar] [CrossRef]
  9. Shobugawa, Y.; Takeuchi, T.; Hibino, A.; Hassan, M.R.; Yagami, R.; Kondo, H.; Odagiri, T.; Saito, R. Occurrence of human respiratory syncytial virus in summer in Japan. Epidemiol. Infect. 2017, 145, 272–284. [Google Scholar] [CrossRef]
  10. Broberg, E.K.; Waris, M.; Johansen, K.; Snacken, R.; Penttinen, P. Seasonality and geographical spread of respiratory syncytial virus epidemics in 15 European countries, 2010 to 2016. Eurosurveillance 2018, 23, 17-00284. [Google Scholar] [CrossRef]
  11. Wagatsuma, K.; Koolhof, I.S.; Shobugawa, Y.; Saito, R. Decreased human respiratory syncytial virus activity during the COVID-19 pandemic in Japan: An ecological time-series analysis. BMC Infect. Dis. 2021, 21, 734. [Google Scholar] [CrossRef] [PubMed]
  12. Ujiie, M.; Tsuzuki, S.; Nakamoto, T.; Iwamoto, N. Resurgence of Respiratory Syncytial Virus Infections during COVID-19 Pandemic, Tokyo, Japan. Emerg. Infect. Dis. 2021, 27, 2969–2970. [Google Scholar] [CrossRef] [PubMed]
  13. Florin, T.A.; Plint, A.C.; Zorc, J.J. Viral bronchiolitis. Lancet 2017, 389, 211–224. [Google Scholar] [CrossRef]
  14. Sloan, C.; Moore, M.L.; Hartert, T. Impact of pollution, climate, and sociodemographic factors on spatiotemporal dynamics of seasonal respiratory viruses. Clin. Transl. Sci. 2011, 4, 48–54. [Google Scholar] [CrossRef]
  15. Tang, J.W.; Loh, T.P. Correlations between climate factors and incidence--a contributor to RSV seasonality. Rev. Med. Virol. 2014, 24, 15–34. [Google Scholar] [CrossRef]
  16. Obando-Pacheco, P.; Justicia-Grande, A.J.; Rivero-Calle, I.; Rodríguez-Tenreiro, C.; Sly, P.; Ramilo, O.; Mejías, A.; Baraldi, E.; Papadopoulos, N.G.; Nair, H.; et al. Respiratory Syncytial Virus Seasonality: A Global Overview. J. Infect. Dis. 2018, 217, 1356–1364. [Google Scholar] [CrossRef]
  17. Li, Y.; Wang, X.; Broberg, E.K.; Campbell, H.; Nair, H. Seasonality of respiratory syncytial virus and its association with meteorological factors in 13 European countries, week 40 2010 to week 39 2019. Eurosurveillance 2022, 27, 2100619. [Google Scholar] [CrossRef]
  18. Nenna, R.; Evangelisti, M.; Frassanito, A.; Scagnolari, C.; Pierangeli, A.; Antonelli, G.; Nicolai, A.; Arima, S.; Moretti, C.; Papoff, P.; et al. Respiratory syncytial virus bronchiolitis, weather conditions and air pollution in an Italian urban area: An observational study. Environ. Res. 2017, 158, 188–193. [Google Scholar] [CrossRef]
  19. Baker, R.E.; Mahmud, A.S.; Wagner, C.E.; Yang, W.; Pitzer, V.E.; Viboud, C.; Vecchi, G.A.; Metcalf, C.J.E.; Grenfell, B.T. Epidemic dynamics of respiratory syncytial virus in current and future climates. Nat. Commun. 2019, 10, 5512. [Google Scholar] [CrossRef]
  20. Pitzer, V.E.; Viboud, C.; Alonso, W.J.; Wilcox, T.; Metcalf, C.J.; Steiner, C.A.; Haynes, A.K.; Grenfell, B.T. Environmental drivers of the spatiotemporal dynamics of respiratory syncytial virus in the United States. PLoS Pathog. 2015, 11, e1004591. [Google Scholar] [CrossRef] [PubMed]
  21. Yu, J.; Liu, C.; Xiao, Y.; Xiang, Z.; Zhou, H.; Chen, L.; Shen, K.; Xie, Z.; Ren, L.; Wang, J. Respiratory Syncytial Virus Seasonality, Beijing, China, 2007–2015. Emerg. Infect. Dis. 2019, 25, 1127–1135. [Google Scholar] [CrossRef] [PubMed]
  22. Zhang, H.; Wen, S.; Zheng, J.; Chen, X.; Lv, F.; Liu, L. Meteorological factors affecting respiratory syncytial virus infection: A time-series analysis. Pediatr. Pulmonol. 2020, 55, 713–718. [Google Scholar] [CrossRef] [PubMed]
  23. Ye, Q.; Fu, J.F.; Mao, J.H.; Shang, S.Q. Haze is a risk factor contributing to the rapid spread of respiratory syncytial virus in children. Environ. Sci. Pollut. Res. Int. 2016, 23, 20178–20185. [Google Scholar] [CrossRef] [PubMed]
  24. Xu, B.; Wang, J.; Li, Z.; Xu, C.; Liao, Y.; Hu, M.; Yang, J.; Lai, S.; Wang, L.; Yang, W. Seasonal association between viral causes of hospitalised acute lower respiratory infections and meteorological factors in China: A retrospective study. Lancet Planet. Health 2021, 5, e154–e163. [Google Scholar] [CrossRef]
  25. Lapeña, S.; Robles, M.B.; Castañón, L.; Martínez, J.P.; Reguero, S.; Alonso, M.P.; Fernández, I. Climatic factors and lower respiratory tract infection due to respiratory syncytial virus in hospitalised infants in northern Spain. Eur. J. Epidemiol. 2005, 20, 271–276. [Google Scholar] [CrossRef]
  26. Tang, J.W.; Lai, F.Y.; Wong, F.; Hon, K.L. Incidence of common respiratory viral infections related to climate factors in hospitalized children in Hong Kong. Epidemiol. Infect. 2010, 138, 226–235. [Google Scholar] [CrossRef]
  27. Fall, A.; Dia, N.; Cisse, E.H.A.; Kiori, D.E.; Sarr, F.D.; Sy, S.; Goudiaby, D.; Richard, V.; Niang, M.N. Epidemiology and Molecular Characterization of Human Respiratory Syncytial Virus in Senegal after Four Consecutive Years of Surveillance, 2012–2015. PLoS ONE 2016, 11, e0157163. [Google Scholar] [CrossRef]
  28. Rodriguez-Martinez, C.E.; Sossa-Briceño, M.P.; Acuña-Cordero, R. Relationship between meteorological conditions and respiratory syncytial virus in a tropical country. Epidemiol. Infect. 2015, 143, 2679–2686. [Google Scholar] [CrossRef]
  29. Djelantik, I.G.; Gessner, B.D.; Soewignjo, S.; Steinhoff, M.; Sutanto, A.; Widjaya, A.; Linehan, M.; Moniaga, V. Incidence and clinical features of hospitalization because of respiratory syncytial virus lower respiratory illness among children less than two years of age in a rural Asian setting. Pediatr. Infect. Dis. J. 2003, 22, 150–157. [Google Scholar] [CrossRef]
  30. Chan, P.K.; Sung, R.Y.; Fung, K.S.; Hui, M.; Chik, K.W.; Adeyemi-Doro, F.A.; Cheng, A.F. Epidemiology of respiratory syncytial virus infection among paediatric patients in Hong Kong: Seasonality and disease impact. Epidemiol. Infect. 1999, 123, 257–262. [Google Scholar] [CrossRef] [PubMed]
  31. Omer, S.B.; Sutanto, A.; Sarwo, H.; Linehan, M.; Djelantik, I.G.; Mercer, D.; Moniaga, V.; Moulton, L.H.; Widjaya, A.; Muljati, P.; et al. Climatic, temporal, and geographic characteristics of respiratory syncytial virus disease in a tropical island population. Epidemiol. Infect. 2008, 136, 1319–1327. [Google Scholar] [CrossRef]
  32. Nyoka, R.; Omony, J.; Mwalili, S.M.; Achia, T.N.O.; 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, e0178323. [Google Scholar] [CrossRef]
  33. Nottmeyer, L.; Armstrong, B.; Lowe, R.; Abbott, S.; Meakin, S.; O’Reilly, K.; von Borries, R.; Schneider, R.; Royé, D.; Hashizume, M.; et al. The association of COVID-19 incidence with temperature, humidity, and UV radiation—A global multi-city analysis. Sci. Total. Environ. 2022, 854, 158636. [Google Scholar] [CrossRef]
  34. Weaver, A.K.; Head, J.R.; Gould, C.F.; Carlton, E.J.; Remais, J.V. Environmental Factors Influencing COVID-19 Incidence and Severity. Annu. Rev. Public Health 2022, 43, 271–291. [Google Scholar] [CrossRef]
  35. Villeneuve, P.J.; Goldberg, M.S. Methodological Considerations for Epidemiological Studies of Air Pollution and the SARS and COVID-19 Coronavirus Outbreaks. Environ. Health Perspect. 2020, 128, 95001. [Google Scholar] [CrossRef] [PubMed]
  36. Rocklöv, J.; Dubrow, R. Climate change: An enduring challenge for vector-borne disease prevention and control. Nat. Immunol. 2020, 21, 479–483. [Google Scholar] [CrossRef] [PubMed]
  37. Gasparrini, A.; Guo, Y.; Hashizume, M.; Lavigne, E.; Zanobetti, A.; Schwartz, J.; Tobias, A.; Tong, S.; Rocklöv, J.; Forsberg, B.; et al. Mortality risk attributable to high and low ambient temperature: A multicountry observational study. Lancet 2015, 386, 369–375. [Google Scholar] [CrossRef]
  38. Gasparrini, A.; Armstrong, B.; Kenward, M.G. Distributed lag non-linear models. Stat. Med. 2010, 29, 2224–2234. [Google Scholar] [CrossRef] [PubMed]
  39. Liyanage, P.; Tozan, Y.; Overgaard, H.J.; Aravinda Tissera, H.; Rocklöv, J. Effect of El Niño-Southern Oscillation and local weather on Aedes dvector activity from 2010 to 2018 in Kalutara district, Sri Lanka: A two-stage hierarchical analysis. Lancet Planet. Health 2022, 6, e577–e585. [Google Scholar] [CrossRef]
  40. Chen, Y.; Li, N.; Lourenço, J.; Wang, L.; Cazelles, B.; Dong, L.; Li, B.; Liu, Y.; Jit, M.; Bosse, N.I.; et al. Measuring the effects of COVID-19-related disruption on dengue transmission in southeast Asia and Latin America: A statistical modelling study. Lancet Infect. Dis. 2022, 22, 657–667. [Google Scholar] [CrossRef]
  41. Kim, Y.; Ratnam, J.V.; Doi, T.; Morioka, Y.; Behera, S.; Tsuzuki, A.; Minakawa, N.; Sweijd, N.; Kruger, P.; Maharaj, R.; et al. Malaria predictions based on seasonal climate forecasts in South Africa: A time series distributed lag nonlinear model. Sci. Rep. 2019, 9, 17882. [Google Scholar] [CrossRef] [PubMed]
  42. Kim, Y.M.; Park, J.W.; Cheong, H.K. Estimated effect of climatic variables on the transmission of Plasmodium vivax malaria in the Republic of Korea. Environ. Health Perspect. 2012, 120, 1314–1319. [Google Scholar] [CrossRef]
  43. Moirano, G.; Gasparrini, A.; Acquaotta, F.; Fratianni, S.; Merletti, F.; Maule, M.; Richiardi, L. West Nile Virus infection in Northern Italy: Case-crossover study on the short-term effect of climatic parameters. Environ. Res. 2018, 167, 544–549. [Google Scholar] [CrossRef] [PubMed]
  44. Soh, S.; Aik, J. The abundance of Culex mosquito vectors for West Nile Virus and other flaviviruses: A time-series analysis of rainfall and temperature dependence in Singapore. Sci. Total Environ. 2021, 754, 142420. [Google Scholar] [CrossRef] [PubMed]
  45. National Institute of Infectious Diseases, Japan (NIID). National Epidemiological Surveillance of Infectious Diseases (NESID) Infectious Diseases Weekly Report (IDWR). Available online: https://www.niid.go.jp/niid/ja/idwr.html (accessed on 1 February 2023).
  46. Hibino, A.; Saito, R.; Taniguchi, K.; Zaraket, H.; Shobugawa, Y.; Matsui, T.; Suzuki, H. Molecular epidemiology of human respiratory syncytial virus among children in Japan during three seasons and hospitalization risk of genotype ON1. PLoS ONE 2018, 13, e0192085. [Google Scholar] [CrossRef]
  47. Zaraket, H.; Saito, R. Japanese Surveillance Systems and Treatment for Influenza. Curr. Treat. Options Infect. Dis. 2016, 8, 311–328. [Google Scholar] [CrossRef]
  48. Ministry of Health, Labour and Welfare, Japan (MHLW). Notification of Physicians and Veterinarians Based on the Infectious Diseases Law (Human Respiratory Syncytial Virus). Available online: http://www.mhlw.go.jp/bunya/kenkou/kekkaku-kansenshou11/01-05-15.html (accessed on 1 February 2023).
  49. Japan Meteorological Agency. Meteorological Data Search. Available online: https://www.data.jma.go.jp/obd/stats/etrn/index.php (accessed on 19 July 2023).
  50. Lee, M.H.; Mailepessov, D.; Yahya, K.; Loo, L.H.; Maiwald, M.; Aik, J. Air quality, meteorological variability and pediatric respiratory syncytial virus infections in Singapore. Sci. Rep. 2023, 13, 1001. [Google Scholar] [CrossRef]
  51. Ali, S.T.; Tam, C.C.; Cowling, B.J.; Yeo, K.T.; Yung, C.F. Meteorological drivers of respiratory syncytial virus infections in Singapore. Sci. Rep. 2020, 10, 20469. [Google Scholar] [CrossRef]
  52. Thongpan, I.; Vongpunsawad, S.; Poovorawan, Y. Respiratory syncytial virus infection trend is associated with meteorological factors. Sci. Rep. 2020, 10, 10931. [Google Scholar] [CrossRef]
  53. Wagatsuma, K.; Koolhof, I.S.; Saito, R. The Relative Roles of Ambient Temperature and Mobility Patterns in Shaping the Transmission Heterogeneity of SARS-CoV-2 in Japan. Viruses 2022, 14, 2232. [Google Scholar] [CrossRef]
  54. Wagatsuma, K.; Koolhof, I.S.; Saito, R. Was the Reduction in Seasonal Influenza Transmission during 2020 Attributable to Non-Pharmaceutical Interventions to Contain Coronavirus Disease 2019 (COVID-19) in Japan? Viruses 2022, 14, 1417. [Google Scholar] [CrossRef] [PubMed]
  55. Chan, Y.H. 2003 Biostatistics 104: Correlational analysis. Singap. Med. J. 2003, 44, 614–619. [Google Scholar]
  56. Gasparrini, A.; Masselot, P.; Scortichini, M.; Schneider, R.; Mistry, M.N.; Sera, F.; Macintyre, H.L.; Phalkey, R.; Vicedo-Cabrera, A.M. Small-area assessment of temperature-related mortality risks in England and Wales: A case time series analysis. Lancet Planet. Health 2022, 6, e557–e564. [Google Scholar] [CrossRef]
  57. Lessler, J.; Reich, N.G.; Brookmeyer, R.; Perl, T.M.; Nelson, K.E.; Cummings, D.A. Incubation periods of acute respiratory viral infections: A systematic review. Lancet Infect. Dis. 2009, 9, 291–300. [Google Scholar] [CrossRef] [PubMed]
  58. Taniguchi, K.; Hashimoto, S.; Kawado, M.; Murakami, Y.; Izumida, M.; Ohta, A.; Tada, Y.; Shigematsu, M.; Yasui, Y.; Nagai, M. Overview of infectious disease surveillance system in Japan, 1999–2005. J. Epidemiol. 2007, 17, S3–S13. [Google Scholar] [CrossRef]
  59. Bhaskaran, K.; Gasparrini, A.; Hajat, S.; Smeeth, L.; Armstrong, B. Time series regression studies in environmental epidemiology. Int. J. Epidemiol. 2013, 42, 1187–1195. [Google Scholar] [CrossRef]
  60. Imai, C.; Hashizume, M. A systematic review of methodology: Time series regression analysis for environmental factors and infectious diseases. Trop. Med. Health 2015, 43, 1–9. [Google Scholar] [CrossRef]
  61. Imai, C.; Armstrong, B.; Chalabi, Z.; Mangtani, P.; Hashizume, M. Time series regression model for infectious disease and weather. Environ. Res. 2015, 142, 319–327. [Google Scholar] [CrossRef]
  62. Prefecture/City/Town/Village. Population Data. Available online: https://uub.jp/rnk/rnk.cgi?T=p&S=j&B=20201001 (accessed on 1 January 2023).
  63. Ma, Y.; Pei, S.; Shaman, J.; Dubrow, R.; Chen, K. Role of meteorological factors in the transmission of SARS-CoV-2 in the United States. Nat. Commun. 2021, 12, 3602. [Google Scholar] [CrossRef]
  64. Gasparrini, A. Distributed Lag Linear and Non-Linear Models in R: The Package dlnm. J. Stat. Softw. 2011, 43, 1–20. [Google Scholar]
  65. Papi, A.; Ison, M.G.; Langley, J.M.; Lee, D.-G.; Leroux-Roels, I.; Martinon-Torres, F.; Schwarz, T.F.; van Zyl-Smit, R.N.; Campora, L.; Dezutter, N.; et al. Respiratory Syncytial Virus Prefusion F Protein Vaccine in Older Adults. N. Engl. J. Med. 2023, 388, 595–608. [Google Scholar] [CrossRef]
  66. Mazur, N.I.; Terstappen, J.; Baral, R.; Bardají, A.; Beutels, P.; Buchholz, U.J.; Cohen, C.; Crowe, J.E., Jr.; Cutland, C.L.; Eckert, L.; et al. Respiratory syncytial virus prevention within reach: The vaccine and monoclonal antibody landscape. Lancet Infect. Dis. 2023, 23, e2–e21. [Google Scholar] [CrossRef] [PubMed]
  67. Lopes, G.P.; Amorim, Í.P.S.; Melo, B.O.; Maramaldo, C.E.C.; Bomfim, M.R.Q.; Lima Neto, L.G.; Alves, M.S.; Silva, F.B.; Soeiro-Pereira, P.V.; Falcai, A. Identification and seasonality of rhinovirus and respiratory syncytial virus in asthmatic children in tropical climate. Biosci. Rep. 2020, 40, BSR20200634. [Google Scholar] [CrossRef]
  68. Polozov, I.V.; Bezrukov, L.; Gawrisch, K.; Zimmerberg, J. Progressive ordering with decreasing temperature of the phospholipids of influenza virus. Nat. Chem. Biol. 2008, 4, 248–255. [Google Scholar] [CrossRef] [PubMed]
  69. Chan, K.H.; Peiris, J.S.; Lam, S.Y.; Poon, L.L.; Yuen, K.Y.; Seto, W.H. The Effects of Temperature and Relative Humidity on the Viability of the SARS Coronavirus. Adv. Virol. 2011, 2011, 734690. [Google Scholar] [CrossRef] [PubMed]
  70. Shaman, J.; Kohn, M. Absolute humidity modulates influenza survival, transmission, and seasonality. Proc. Natl. Acad. Sci. USA 2009, 106, 3243–3248. [Google Scholar] [CrossRef]
  71. Lowen, A.C.; Mubareka, S.; Steel, J.; Palese, P. Influenza virus transmission is dependent on relative humidity and temperature. PLoS Pathog. 2007, 3, 1470–1476. [Google Scholar] [CrossRef]
  72. Bloom-Feshbach, K.; Alonso, W.J.; Charu, V.; Tamerius, J.; Simonsen, L.; Miller, M.A.; Viboud, C. Latitudinal variations in seasonal activity of influenza and respiratory syncytial virus (RSV): A global comparative review. PLoS ONE 2013, 8, e54445. [Google Scholar] [CrossRef]
  73. Yusuf, S.; Piedimonte, G.; Auais, A.; Demmler, G.; Krishnan, S.; Van Caeseele, P.; Singleton, R.; Broor, S.; Parveen, S.; Avendano, L.; et al. The relationship of meteorological conditions to the epidemic activity of respiratory syncytial virus. Epidemiol. Infect. 2007, 135, 1077–1090. [Google Scholar] [CrossRef]
  74. Hall, C.B.; Douglas, R.G., Jr.; Geiman, J.M. Possible transmission by fomites of respiratory syncytial virus. J. Infect. Dis. 1980, 141, 98–102. [Google Scholar] [CrossRef]
  75. Vandini, S.; Corvaglia, L.; Alessandroni, R.; Aquilano, G.; Marsico, C.; Spinelli, M.; Lanari, M.; Faldella, G. Respiratory syncytial virus infection in infants and correlation with meteorological factors and air pollutants. Ital. J. Pediatr. 2013, 39, 1. [Google Scholar] [CrossRef]
  76. Paynter, S.; Yakob, L.; Simões, E.A.; Lucero, M.G.; Tallo, V.; Nohynek, H.; Ware, R.S.; Weinstein, P.; Williams, G.; Sly, P.D. Using mathematical transmission modelling to investigate drivers of respiratory syncytial virus seasonality in children in the Philippines. PLoS ONE 2014, 9, e90094. [Google Scholar] [CrossRef] [PubMed]
  77. Coccia, M. How do low wind speeds and high levels of air pollution support the spread of COVID-19? Atmos. Pollut. Res. 2021, 12, 437–445. [Google Scholar] [CrossRef] [PubMed]
  78. Yuan, J.; Wu, Y.; Jing, W.; Liu, J.; Du, M.; Wang, Y.; Liu, M. Association between meteorological factors and daily new cases of COVID-19 in 188 countries: A time series analysis. Sci. Total Environ. 2021, 780, 146538. [Google Scholar] [CrossRef] [PubMed]
  79. Murray, E.L.; Klein, M.; Brondi, L.; McGowan, J.E., Jr.; van Mels, C.; Brooks, W.A.; Kleinbaum, D.; Goswami, D.; Ryan, P.B.; Bridges, C.B. Rainfall, household crowding, and acute respiratory infections in the tropics. Epidemiol. Infect. 2012, 140, 78–86. [Google Scholar] [CrossRef] [PubMed]
  80. Dowell, S.F. Seasonality—Still confusing. Epidemiol. Infect. 2012, 140, 87–90. [Google Scholar] [CrossRef]
  81. Tang, J.W. The effect of environmental parameters on the survival of airborne infectious agents. J. R. Soc. Interface 2009, 6 (Suppl. S6), S737–S746. [Google Scholar] [CrossRef]
  82. Lei, H.; Yang, M.; Dong, Z.; Hu, K.; Chen, T.; Yang, L.; Zhang, N.; Duan, X.; Yang, S.; Wang, D.; et al. Indoor relative humidity shapes influenza seasonality in temperate and subtropical climates in China. Int. J. Infect. Dis. 2023, 126, 54–63. [Google Scholar] [CrossRef]
  83. Verheyen, C.A.; Bourouiba, L. Associations between indoor relative humidity and global COVID-19 outcomes. J. R. Soc. Interface 2022, 19, 20210865. [Google Scholar] [CrossRef]
  84. White, L.J.; Waris, M.; Cane, P.A.; Nokes, D.J.; Medley, G.F. The transmission dynamics of groups A and B human respiratory syncytial virus (hRSV) in England & Wales and Finland: Seasonality and cross-protection. Epidemiol. Infect. 2005, 133, 279–289. [Google Scholar] [CrossRef]
  85. Waris, M.; White, L.J. Seasonality of respiratory syncytial virus infection. Clin. Infect. Dis. 2006, 43, 541. [Google Scholar] [CrossRef] [PubMed]
  86. Paiva, T.M.; Ishida, M.A.; Benega, M.A.; Constantino, C.R.; Silva, D.B.; Santos, K.C.; Oliveira, M.I.; Barbosa, H.A.; Carvalhanas, T.R.; Schuck-Paim, C.; et al. Shift in the timing of respiratory syncytial virus circulation in a subtropical megalopolis: Implications for immunoprophylaxis. J. Med. Virol. 2012, 84, 1825–1830. [Google Scholar] [CrossRef] [PubMed]
  87. Foley, D.A.; Yeoh, D.K.; Minney-Smith, C.A.; Martin, A.C.; Mace, A.O.; Sikazwe, C.T.; Le, H.; Levy, A.; Moore, H.C.; Blyth, C.C. The Interseasonal Resurgence of Respiratory Syncytial Virus in Australian Children Following the Reduction of Coronavirus Disease 2019-Related Public Health Measures. Clin. Infect. Dis. 2021, 73, e2829–e2830. [Google Scholar] [CrossRef]
  88. Casalegno, J.S.; Ploin, D.; Cantais, A.; Masson, E.; Bard, E.; Valette, M.; Fanget, R.; Targe, S.C.; Myar-Dury, A.F.; Doret-Dion, M.; et al. Characteristics of the delayed respiratory syncytial virus epidemic, 2020/2021, Rhône Loire, France. Eurosurveillance 2021, 26, 2100630. [Google Scholar] [CrossRef] [PubMed]
  89. Onozuka, D. The influence of diurnal temperature range on the incidence of respiratory syncytial virus in Japan. Epidemiol. Infect. 2015, 143, 813–820. [Google Scholar] [CrossRef] [PubMed]
  90. Boogaard, I.; van Oosten, M.; van Rijt, L.S.; Muskens, F.; Kimman, T.G.; Lambrecht, B.N.; Buisman, A.M. Respiratory syncytial virus differentially activates murine myeloid and plasmacytoid dendritic cells. Immunology 2007, 122, 65–72. [Google Scholar] [CrossRef] [PubMed]
  91. Paynter, S.; Weinstein, P.; Ware, R.S.; Lucero, M.G.; Tallo, V.; Nohynek, H.; Barnett, A.G.; Skelly, C.; Simões, E.A.; Sly, P.D.; et al. Sunshine, rainfall, humidity and child pneumonia in the tropics: Time-series analyses. Epidemiol. Infect. 2013, 141, 1328–1336. [Google Scholar] [CrossRef]
  92. Lowe, R.; Lee, S.A.; O’Reilly, K.M.; Brady, O.J.; Bastos, L.; Carrasco-Escobar, G.; de Castro Catão, R.; Colón-González, F.J.; Barcellos, C.; Carvalho, M.S.; et al. Combined effects of hydrometeorological hazards and urbanisation on dengue risk in Brazil: A spatiotemporal modelling study. Lancet Planet. Health 2021, 5, e209–e219. [Google Scholar] [CrossRef]
  93. Alesemi, M.; Shahrani, J.S.A.; Iqbal, N.; Shah, R.; Nonlaopon, K. Analysis and Numerical Simulation of System of Fractional Partial Differential Equations with Non-Singular Kernel Operators. Symmetry 2023, 15, 233. [Google Scholar] [CrossRef]
  94. Botmart, T.; Naeem, M.; Shah, R.; Iqbal, N. Fractional View Analysis of Emden-Fowler Equations with the Help of Analytical Method. Symmetry 2022, 14, 2168. [Google Scholar] [CrossRef]
  95. Yasmin, H.; Aljahdaly, N.H.; Saeed, A.M.; Shah, R. Investigating Families of Soliton Solutions for the Complex Structured Coupled Fractional Biswas–Arshed Model in Birefringent Fibers Using a Novel Analytical Technique. Fractal Fract. 2023, 7, 491. [Google Scholar] [CrossRef]
  96. Yasmin, H.; Aljahdaly, N.H.; Saeed, A.M.; Shah, R. Probing Families of Optical Soliton Solutions in Fractional Perturbed Radhakrishnan–Kundu–Lakshmanan Model with Improved Versions of Extended Direct Algebraic Method. Fractal Fract. 2023, 7, 512. [Google Scholar] [CrossRef]
  97. The Lancet Planetary Health. Modelling health futures. Lancet Planet. Health 2021, 5, e395. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Assessing the pooled nonlinear association of IRRs of HRSV incidence with meteorological variables. (A) Overall association of the 2-week cumulative risk of the percent change in the estimated human respiratory syncytial virus (HRSV) infection incidence with weekly mean ambient temperature (unit: °C). (B) Overall association of the 2-week cumulative risk of the percent change in the estimated HRSV infection incidence with weekly relative humidity (unit: %). (C) Overall association of the 2-week cumulative risk of the percent change in the estimated HRSV infection incidence with weekly precipitation (unit: mm). (D) Overall association of the 2-week cumulative risk of the percent change in the estimated HRSV infection incidence with daily weekly wind speed (unit: m/s). The present study covers the period between 1 January 2014 and 29 November 2019 (between the 1st week of 2014 and 52nd week of 2019) across all 47 prefectures in Japan. The red, yellow, blue, and green lines represent the estimated incidence relative risks (IRRs) of HRSV infection, with the shaded bands representing the 95% confidence intervals (CIs). The corresponding reference values are 16.3 °C (A), 70.2% (B), 0.0 mm (C), and 0.0 m/s (D). The relevant selected estimates (i.e., IRR and 95% CIs) for this figure are shown in Table 2.
Figure 1. Assessing the pooled nonlinear association of IRRs of HRSV incidence with meteorological variables. (A) Overall association of the 2-week cumulative risk of the percent change in the estimated human respiratory syncytial virus (HRSV) infection incidence with weekly mean ambient temperature (unit: °C). (B) Overall association of the 2-week cumulative risk of the percent change in the estimated HRSV infection incidence with weekly relative humidity (unit: %). (C) Overall association of the 2-week cumulative risk of the percent change in the estimated HRSV infection incidence with weekly precipitation (unit: mm). (D) Overall association of the 2-week cumulative risk of the percent change in the estimated HRSV infection incidence with daily weekly wind speed (unit: m/s). The present study covers the period between 1 January 2014 and 29 November 2019 (between the 1st week of 2014 and 52nd week of 2019) across all 47 prefectures in Japan. The red, yellow, blue, and green lines represent the estimated incidence relative risks (IRRs) of HRSV infection, with the shaded bands representing the 95% confidence intervals (CIs). The corresponding reference values are 16.3 °C (A), 70.2% (B), 0.0 mm (C), and 0.0 m/s (D). The relevant selected estimates (i.e., IRR and 95% CIs) for this figure are shown in Table 2.
Viruses 15 01914 g001
Table 1. Descriptive statistics for the number of weekly newly confirmed HRSV cases and meteorological variables.
Table 1. Descriptive statistics for the number of weekly newly confirmed HRSV cases and meteorological variables.
Potential DriversMeanSDMinP25P50P75Max
Weekly newly confirmed cases4976072258892
Mean ambient temperature (°C)15.78.3–6.78.116.322.732.2
Relative humidity (%)69.59.530.863.170.276.497.2
Precipitation (mm)4.76.40.00.72.96.294.1
Wind speed (m/s)2.90.90.92.22.83.411.3
Abbreviations: SD, standard deviation; Min: minimum; P25, 25th percentile; P25, 50th percentile; P75, 75th percentile; Max: maximum. Notes: the present study covers the period between 1 January 2014 and 29 November 2019 (between the 1st week of 2014 and 52nd week of 2019) across all 47 prefectures in Japan.
Table 2. Forecasting the specific IRRs of nonlinear associations between the HRSV incidence and meteorological variables.
Table 2. Forecasting the specific IRRs of nonlinear associations between the HRSV incidence and meteorological variables.
Potential DriversLag (Weeks)
0120–2
IRR (95% CI)IRR (95% CI)IRR (95% CI)IRR (95% CI)
Mean ambient temperature (°C)
8.1 °C1.053
(1.015, 1.093)
1.117
(1.070, 1.166)
0.869
(0.837, 0.903)
1.023
(1.008, 1.039)
22.7 °C0.936
(0.910, 0.964)
1.037
(1.001, 1.073)
1.099
(1.067, 1.131)
1.068
(1.056, 1.079)
Relative humidity (%)
63.1%0.998
(0.986, 1.011)
0.979
(0.966, 0.991)
1.009
(0.997, 1.022)
0.997
(0.971, 1.003)
76.4%1.007
(0.997, 1.017)
1.039
(1.029, 1.050)
0.998
(0.988, 1.008)
1.045
(1.032, 1.059)
Precipitation (mm)
0.7 mm1.002
(0.995, 1.008)
1.012
(1.005, 1.019)
1.015
(1.005, 1.019)
1.027
(1.015, 1.038)
6.2 mm0.983
(0.964, 1.002)
1.033
(1.014, 1.053)
1.047
(1.027, 1.067)
1.064
(1.032, 1.097)
Wind speed (m/s)
2.2 m/s0.972
(0.877, 1.078)
0.887
(0.799, 0.985)
0.919
(0.829, 1.019)
0.793
(0.697, 0.903)
3.4 m/s0.976
(0.875, 1.088)
0.864
(0.774, 0.965)
0.918
(0.824, 1.022)
0.775
(0.678, 0.885)
Abbreviations: IRR, incidence relative risk; CI, confidence interval; HRSV, human respiratory virus. Notes: the present study covers the period between 1 January 2014 and 29 November 2019 (between the 1st week of 2014 and 52nd week of 2019) across all 47 prefectures in Japan. Associations between the HRSV infection incidence and mean ambient temperature (unit: °C), relative humidity (unit: %), precipitation (unit: mm), and wind speed (unit: m/s) are described as IRRs with 95% confidence intervals (CIs) with reference to the IRRs at 16.3 °C, 70.2%, 0.0 mm, and 0.0 m/s, respectively; 8.1 °C and 22.7 °C correspond to the 25th and 75th percentiles of mean ambient temperature, respectively; 63.1% and 76.4% correspond to the 25th and 75th percentiles of relative humidity, respectively; 0.7 mm and 6.2 mm correspond to the 25th and 75th percentiles of precipitation, respectively; 2.2 m/s and 3.4 m/s correspond to the 25th and 75th percentiles of wind speed, respectively.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wagatsuma, K.; Koolhof, I.S.; Saito, R. Nonlinear and Multidelayed Effects of Meteorological Drivers on Human Respiratory Syncytial Virus Infection in Japan. Viruses 2023, 15, 1914. https://doi.org/10.3390/v15091914

AMA Style

Wagatsuma K, Koolhof IS, Saito R. Nonlinear and Multidelayed Effects of Meteorological Drivers on Human Respiratory Syncytial Virus Infection in Japan. Viruses. 2023; 15(9):1914. https://doi.org/10.3390/v15091914

Chicago/Turabian Style

Wagatsuma, Keita, Iain S. Koolhof, and Reiko Saito. 2023. "Nonlinear and Multidelayed Effects of Meteorological Drivers on Human Respiratory Syncytial Virus Infection in Japan" Viruses 15, no. 9: 1914. https://doi.org/10.3390/v15091914

APA Style

Wagatsuma, K., Koolhof, I. S., & Saito, R. (2023). Nonlinear and Multidelayed Effects of Meteorological Drivers on Human Respiratory Syncytial Virus Infection in Japan. Viruses, 15(9), 1914. https://doi.org/10.3390/v15091914

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop