Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Med., 10 October 2025

Sec. Family Medicine and Primary Care

Volume 12 - 2025 | https://doi.org/10.3389/fmed.2025.1643828

This article is part of the Research TopicEmerging Global Population Health Risks: From Epidemiological Perspectives Volume IIView all 6 articles

Influence of meteorological conditions on herpes zoster occurrence: a retrospective cohort study

Micha&#x; Ochal
&#x;Michał Ochal1*Katarzyna Gli&#x;ska Lewczuk&#x;Katarzyna Glińska Lewczuk2Ewa Draga&#x;ska&#x;Ewa Dragańska2Iwona Cymes&#x;Iwona Cymes2Jerzy Romaszko&#x;Jerzy Romaszko1
  • 1Department of Family Medicine and Infectious Diseases, Collegium Medicum, School of Medicine, University of Warmia and Mazury in Olsztyn, Olsztyn, Poland
  • 2Department of Water Management and Climatology, Faculty of Agriculture and Forestry, University of Warmia and Mazury in Olsztyn, Olsztyn, Poland

Introduction: Herpes zoster is caused by the reactivation of the varicella zoster virus, typically affecting older adults and immunocompromised individuals. Although various studies have examined potential triggers, the influence of meteorological conditions on herpes zoster reporting remains unclear. We aimed to investigate whether meteorological and biometeorological factors, including the Universal Thermal Climate Index (UTCI) and the H index (air cooling power), have an effect on herpes zoster occurrence.

Methodology: We performed a retrospective analysis of medical records for herpes zoster cases (ICD-10 code B02) from 2009 to 2023 at two medical facilities in northeastern Poland. Local meteorological data were obtained from regional stations. Statistical analyses accounted for patient age, seasonality, and thermal conditions (as measured by UTCI).

Results and discussion: No seasonal pattern in herpes zoster reporting was observed, and there were no statistically significant correlations between basic meteorological variables (ambient temperature, relative humidity, wind speed, etc.) and the number of cases. However, results showed that under strong heat stress conditions (UTCI class 8), the risk of herpes zoster was approximately 20–25% higher compared to thermally neutral days.

Conclusion: Herpes zoster occurrence does not exhibit a seasonal pattern. Composite indices like the UTCI and the H index can be useful tools for predicting the risk of herpes zoster. The relative risk of herpes zoster is lower during cold stress conditions and higher during heat stress conditions.

Introduction

Herpes zoster (shingles) is caused by reactivation of the varicella zoster virus (VZV), which remains latent in sensory ganglia after a primary varicella (chickenpox) infection (1). Primary VZV infection can be especially severe in adults and immunocompromised individuals (2), but after recovery the virus adopts a dormant state instead of being eliminated from the body (3). Virus reactivation into herpes zoster is often precipitated by a weakened immune system due to factors such as aging, stress, chronic disease, or immunosuppressive therapy (46). In particular, cell-mediated immunity wanes with age, reducing the body’s ability to contain the virus and thereby increasing the risk of reactivation. A crucial role in maintaining VZV latency is played by CD8 + T lymphocytes. These immune cells provide surveillance by secreting cytokines such as interferon-gamma (IFN-γ), which enhance neuronal resistance to viral replication (7). Additionally, pro-inflammatory cytokines (interleukin-1, interleukin-6, and tumor necrosis factor-α) are implicated in VZV reactivation: elevated levels of these mediators in response to stress or infection can disrupt the neuronal environment, promoting viral replication (8). Clinically, herpes zoster manifests as a painful vesicular rash with fluid-filled blisters that eventually rupture, crust, and heal. The rash is usually limited to one side of the body, reflecting the reactivation of the virus in a specific sensory ganglion (dermatomal distribution) (9). The virus in the lesions is transmissible; contact with the vesicles can cause varicella in susceptible individuals. A common complication of herpes zoster is postherpetic neuralgia (PHN), a chronic neuropathic pain that can persist for months or even years after the rash resolves (10). The risk of PHN is particularly high in older patients (10).

Considerably less is known about the influence of environmental factors especially weather conditions—on the reactivation of VZV (11). For example, one study found that vitamin D supplementation in hemodialysis patients was associated with a lower occurrence of herpes zoster, suggesting that vitamin D deficiency (potentially linked to reduced sun exposure in winter) might increase reactivation risk (12). In contrast, an analysis by Berlinberg et al. (USA) noted a peak in herpes zoster cases during the summer months (13), implying a possible seasonal effect. Other studies have explored specific meteorological factors: Choi et al. (South Korea) reported that higher ambient temperatures were associated with increased zoster risk, consistent with a heat-stress effect on immunity (14), whereas Yang et al. (China, Shanghai) observed more frequent herpes zoster cases in winter, when colder temperatures and reduced sunlight (hence lower vitamin D levels) might weaken immune defenses (15). Notably, Kawai et al. (USA) found that higher ambient ultraviolet radiation (UVR) exposure was associated with an increased risk of herpes zoster in men, suggesting that UVR-induced immunosuppression may play a role in VZV reactivation (16). It should be noted that these studies quantified thermal stress using single meteorological parameters, such as mean air temperature or UV radiation, rather than composite indices. Such methodological differences may contribute to the heterogeneity of reported findings and highlight the advantage of using the UTCI in our analysis. These inconsistencies across studies highlight that the impact of meteorological factors on herpes zoster remains unclear and potentially context-dependent. Differences in study populations, geographic regions, climate types, and methodologies may explain why some studies found seasonal or thermal effects while others did not. This unresolved scientific picture underscores the need for further research, particularly in regions with moderate continental climates, where such analyses remain limited. Understanding whether and how weather conditions affect herpes zoster reporting is important for public health, as this knowledge could inform new prevention strategies for high-risk groups (17). Moreover, with global climate change leading to more frequent extreme weather events, clarifying the impact of meteorological factors on health is becoming increasingly crucial (18). In Poland, recent decades have witnessed a growing occurrence of extreme weather events, reflecting broader climate variability in Central Europe (19). Although the pan-European heatwave of 2003 had relatively limited direct mortality impact in Poland, the record warm summers of 2010 and several other heatwaves (1992, 1994, 2006, 2014, 2015) did significantly affect public health (19). Poland has also experienced devastating flood events, including the “Millennium Flood” of 1997 and the widespread inundations in 2010, both associated with intense and prolonged precipitation (20). Additionally, meteorological droughts have been documented in Polish river catchments between 1981 and 2010, impacting water resources and agricultural productivity (21). These events illustrate how, despite a temperate continental climate, Poland already contends with climate extremes that are relevant to interpreting meteorological influences on health outcomes.

In the context of herpes zoster, shifting climate patterns may alter any seasonal occurrence trends and could even exacerbate disease occurrence under specific conditions (22, 23). In light of these considerations, the aim of this study was to analyze the effect of meteorological and biometeorological conditions on the occurrence of herpes zoster. The concept of ‘biometeorological’ conditions refers to the interdisciplinary field of biometeorology, which investigates how weather and climate interact with human physiology and health, and provides the framework for applying indices such as the UTCI (24) or air cooling power (H) (25). This work specifically evaluates whether composite weather indices (such as the UTCI) are associated with changes in herpes zoster occurrence, thereby providing a foundation for future research toward more effective preventive strategies.

Materials and methods

Study population

In this retrospective study, data were collected from two medical facilities in northeastern Poland: the Non-Public Health Care Centre “Pantamed” in Olsztyn and the Władysław Biegański County Hospital in Iława (including its Night and Holiday Emergency Health Care Unit, NiŚOZ). In Poland, Night and Holiday Emergency Health Care Units (NiŚOZ) provide urgent outpatient care during nights, weekends, and public holidays, functioning similarly to urgent care centers in other healthcare systems (26). The two sites are about 70 km apart. Pantamed serves approximately 13,000 registered patients (providing ~60,000 consultations per year), while the Iława hospital’s NiŚOZ unit serves a county population of about 90,000, with 15,595 patient admissions to NiŚOZ recorded during the study period (27). From these facilities’ records, we extracted all diagnosed cases of herpes zoster (ICD-10 code B02) in adults for the years 2009–2023 (28). Information on comorbidities and other diagnoses was available in medical records but was not systematically included as covariates in the present analysis. For each case, the date of diagnosis, patient age, sex, and relevant medical history (e.g., prior visits and diagnoses) were recorded. We excluded repeat entries for the same patient’s zoster episode within a 30-day period to avoid double-counting. In total, 1,102 unique adult cases of herpes zoster were analyzed.

Climate characteristics

Meteorological data for the period 2009–2023 were obtained from the Institute of Meteorology and Water Management station in Olsztyn (Poland). According to the standards of the World Meteorological Organization (WMO, 2017) and the Polish Institute of Meteorology and Water Management (IMGW-PIB, 2015), synoptic stations provide measurements representative of the broader geographical area; therefore, the Olsztyn station data were considered reliable for both Olsztyn and Iława (29, 30). According to the updated Köppen–Geiger climate classification, this region has a cold continental climate without a dry season and with a warm summer (Dfb) (31). The daily meteorological parameters collected included air temperature (daily mean, maximum, and minimum), wind speed, relative humidity, and atmospheric pressure. We used these data (specifically the measurements at 12:00 UTC each day for temperature, humidity, vapor pressure, cloud cover, and wind) to calculate the UTCI for each day. It is expressed in °C and represents a composite indicator derived from a human heat balance model that integrates meteorological variables to quantify cold and heat stress (32). UTCI has previously been linked to adverse health effects, including increased mortality, cardiovascular hospitalizations, and, in emerging evidence, infectious diseases (33). In particular, Di Napoli et al. (Europe-wide study) demonstrated across Europe that UTCI values above moderate to strong stress levels (approximately 26–32 °C) were associated with elevated death counts, underscoring its relevance in capturing heat-related health risks (33).

The UTCI is expressed in °C and reflects the integrated effect of weather conditions on the human body’s thermal balance (34). Based on UTCI values, thermal stress levels are categorized into classes. Cold stress conditions occur at low temperatures combined with wind (UTCI classes 1–5, corresponding to ≤9.0 °C), while heat stress occurs at high temperatures with intense solar radiation (heat stress classes 7–10, >26.0 °C). Conditions in the intermediate range (UTCI 9.1–26.0 °C, class 6) indicate thermal comfort (no significant heat or cold stress) (24, 35). In our dataset, the extreme UTCI classes 1, 9, and 10 did not occur during the study period and were thus excluded from analysis. In addition to UTCI, we calculated the H index (air cooling power, expressed in W/m2) for each day at 12:00. This index quantifies the rate of heat loss from a unit surface of the human body due to convection and conduction under the influence of air temperature and wind speed (31). We applied two classes of the cooling power index: cold (H > 1,260 W/m2) and very hot (H < 210 W/m2). Both UTCI and H index values were computed using the BioKlima software (36).

Statistical analysis

Analyses of distribution and group differences were performed using Statistica v13.0 (Dell Inc., Tulsa, USA). Data distribution within groups was first assessed with the Shapiro–Wilk test for normality (p < 0.05 indicating significance). Homogeneity of variances was evaluated by Levene’s test (p < 0.05). Since most variables did not meet the assumption of normal distribution and homoscedasticity, we employed non-parametric methods. To compare differences among groups (e.g., monthly case numbers, age groups, or classes based on UTCI and H) we applied the Kruskal–Wallis one-way analysis of variance. When significant overall differences were found, Dunn’s post hoc test with Bonferroni correction was applied for multiple pairwise comparisons. A two-tailed p-value <0.05 was considered statistically significant for all tests. In addition, Spearman’s rank correlation coefficients were computed to assess the relationships between the number of herpes zoster cases and individual meteorological variables (e.g., temperature, humidity, precipitation, wind) as well as the UTCI and H. This correlation analysis was used to gage the strength and direction of any linear associations between weather factors and zoster occurrence.

To examine how the occurrence of herpes zoster relates to multiple meteorological variables simultaneously, we performed a multivariate analysis using Principal Coordinate Analysis (PCoA) in XLSTAT by Lumivero (version 2023.2.0) (37). PCoA was appliedto visualize similarities or dissimilarities in the dataset. Prior to the analysis, variables measured on different scales were standardized and square-root transformed to reduce skewness. We constructed a Euclidean distance matrix, and PCoA was then carried out to identify principal axes of variation. PCoA provided the percentage of variance explained by each axis, allowing us to evaluate how well the reduced-dimensional representation captured the structure of the original multivariate data.

Furthermore, we calculated the relative risk (RR) of herpes zoster under different thermal conditions defined by UTCI. Specifically, we compared the risk of herpes zoster during periods of cold stress and heat stress with the risk under thermoneutral conditions [UTCI class 6 (9–26 °C)] which served as the reference category for the total cohort and separately for patients aged >65 years. In addition to RR values with corresponding 95% CIs, we also report Z-scores and p-values obtained from pairwise comparisons. We considered an RR statistically significant if the 95% confidence interval (CI) did not include 1. RR calculations were performed using PAST software (version 4.03). We chose to use RR (a risk ratio) in this context, rather than an odds ratio, because our data encompassed the entire population at risk in the study period rather than a case–control sample. Using RR provides a more direct measure of the effect size of thermal stress on disease occurrence in this retrospective cohort analysis.

Results

During the 2009–2023 period, herpes zoster predominantly affected older adults. The occurrence of herpes zoster increased significantly with age, with the highest rates observed in the oldest age group (>65 years) (Figure 1).

Figure 1
Bar graph showing the mean daily number of HZ cases by age class. The age groups are 18-35 years, 36-45 years, 46-55 years, 56-65 years, and over 65 years. The over 65 years group has the highest mean, followed by 56-65 years, 18-35 years, 46-55 years, and 36-45 years. Error bars represent variability, with letters indicating statistical significance between groups.

Figure 1. Average number of admissions for herpes zoster by age class. Values with different superscripts significantly differ according to the Kruskal–Wallis test, followed by Dunn’s post hoc test (p < 0.05).

This group was significantly more affected than all younger age groups (Dunn’s post hoc test, p < 0.0001). Patients aged 56–65 years also showed a higher reporting compared to younger adults, though not as pronounced as the >65 group.

No significant seasonal pattern in herpes zoster occurrence was detected. While the average number of cases appeared slightly higher in the summer months, the variation was not statistically significant (Kruskal–Wallis test, p > 0.05; Figure 2).

Figure 2
Strip plot showing Universal Thermal Climate Index (UTCI) and mean daily herpes zoster cases across months. Green lines represent UTCI and red crosses indicate mean temperature. Black dots with error bars display mean and standard error of herpes zoster cases, with values ranging from 0.15 to 0.30 cases per day. Temperature on left axis spans negative thirty degrees Celsius to forty degrees Celsius.

Figure 2. Monthly distribution of the average number of herpes zoster admissions per day (black dots) and whiskers (±SE) on the background of daily UTCI strip plots (green).

For example, differences in occurrence between February and April or between August and November were not significant, even after adjusting for the different number of days in each month (data not shown). Consistent with the lack of clear seasonality, correlation analysis showed no meaningful association between the number of herpes zoster cases and the meteorological variables examined. None of the Spearman’s correlation coefficients between case counts and factors such as temperature, humidity, precipitation, or wind speed reached statistical significance. Even the correlation with the composite thermal index (UTCI) was essentially zero (r = 0.023, p > 0.05), indicating a negligible effect of thermal comfort on herpes zoster reporting.

The multivariate PCoA provided further insight (Figure 3). In the ordination, the occurrence of herpes zoster (for the total cohort) aligned more closely with certain variables such as the number of days with heat discomfort (Hhot) and relative humidity (RH), suggesting a positive association with these factors along the first principal coordinate (F1). In contrast, temperature-related variables (mean, maximum, and minimum air temperature, as well as the UTCI) plotted in the opposite direction from herpes zoster occurance on the F1 axis, indicating little to no direct association between simple temperature measures and the occurrence of the disease.

Figure 3
Bar chart illustrating the mean daily number of herpes zoster cases per day across different Universal Thermal Climate Index (UTCI) classes. Two bars per class represent the total cohort and those older than sixty-five years. Higher cases are observed in higher UTCI classes, particularly UTCI eight, with a noticeable increase in those over sixty-five.

Figure 3. Principal Coordinate Analysis (PCoA) plot showing the relationship between weather variables, the Universal Thermal Climate Index (UTCI), H index, and the occurrence of herpes zoster. The axes (F1 and F2) represent the primary dimensions of variability, with F1 capturing 43.74% and F2 capturing 16.27% of the variance. The scree plot indicates the eigenvalues and cumulative variability explained by each axis, with F1 and F2 together accounting for 60.01% of the variability in the dataset.

We next examined herpes zoster reporting across categories of thermal conditions as defined by UTCI, comparing the entire study population to the subset of older patients (>65 years) (Figure 4). We observed that extreme heat stress conditions were associated with higher numbers of cases. In both the overall population and in seniors, the average daily number of admissions was lowest under moderate thermal conditions (UTCI class 6) and increased markedly under heat stress conditions (peaking at UTCI classes 7 and 8). The >65 years group had fewer cases per day than the total population at each UTCI category, but showed a similar upward trend across UTCI classes. This indicates that very hot weather conditions correspond to an increase in herpes zoster cases in all adults, although the absolute occurance remains lower in the elderly subgroup.

Figure 4
Bar and line graph depicting Relative Risk (RR) and the number of herpes zoster (HZ) cases per Universal Thermal Climate Index (UTCI) class. The bar graph shows RR for the total cohort and those over 65. Line graphs indicate HZ cases for all patients and those over 65. RR peaks at UTCI 6, while HZ cases for all patients are highest at UTCI 6, and for patients over 65 at UTCI 5.

Figure 4. Average daily number of herpes zoster admissions in relation to UTCI class in the general population and in the group of patients aged >65. The lighter bars correspond to the total cohort, while the darker bars represent patients aged >65.

In addition to case counts, we calculated the relative risk of herpes zoster for different UTCI categories (Figure 5; Supplementary Table S2). In the general population, the relative risk of zoster increased steadily with warmer UTCI classes, reaching the highest values under the most extreme heat stress (UTCI class 8). For the >65 age group, a similar trend of rising RR with higher UTCI was observed, though the RR values for seniors were consistently lower than those for the total population at corresponding classes. In both the overall and >65 groups, the relative risk was markedly elevated in heat stress conditions compared to neutral conditions. Concordantly, the number of herpes zoster admissions per UTCI class was highest at the extreme heat category (UTCI class 8) in both populations. Thus, extreme thermal conditions were associated with an increased risk and frequency of herpes zoster, with the effect somewhat attenuated among older individuals. For example, the risk of herpes zoster was approximately 20–25% higher on days classified as strong heat stress (UTCI class 8) compared to thermally neutral days, whereas it was slightly lower under milder cold stress conditions (UTCI class 3), relative to neutral. Notably, these meteorological associations emerged despite the lack of any overt seasonal pattern in reporting, indicating that herpes zoster occurrence is not governed by season per se but can be modulated by specific thermal conditions.

Figure 5
Principal Coordinate Analysis plot with axes F1 and F2. Data points include air pressure, Herpes zoster, temperature metrics, wind speed, and relative humidity. Inset scree plot shows eigenvalues and cumulative variability across axes F1 to F10, with a red line indicating growth in cumulative variability.

Figure 5. Relationship between herpes zoster admissions and Universal Thermal Climate Index (UTCI) classes. Bars represent the relative risk (RR) of herpes zoster admissions for the total cohort (grey bars) and patients aged >65 years (red dashed line). The dashed green line represents the number of admissions per UTCI class for the total cohort, while the dotted red line represents the number of admissions for seniors aged >65 years. The graph did not show UTCI 2 due to a lack of patients >65 years.

Discussion

This study examined the relationship between weather conditions and herpes zoster occurrence, with a particular focus on using the composite index UTCI to capture thermal stress. In the overall analysis using straightforward meteorological variables (temperature, humidity, etc.), we did not find any significant association with herpes zoster case numbers, nor any clear seasonal trend. However, when we categorized days by UTCI-defined thermal stress levels, a distinct pattern emerged: periods of heat stress corresponded to higher rates of herpes zoster reactivation, whereas periods of cold stress were associated with lower rates. This contrast suggests that while individual weather parameters alone may not show an effect, the combined effect captured by UTCI can reveal an influence of extreme thermal conditions on herpes zoster occurrence. The multivariate PCoA supports this interpretation, as it showed herpes zoster cases clustering with indices of heat stress in the multidimensional analysis.

Unlike using temperature alone, the UTCI incorporates multiple meteorological factors (temperature, humidity, wind, and solar radiation) to assess human thermal stress, making it a more precise indicator of weather conditions relevant to health outcomes (38). Our findings underscore the value of such a comprehensive index: even when no single weather variable showed a strong effect, the integrated thermal load (as represented by UTCI) was able to distinguish conditions that predispose to viral reactivation. In contrast, the H index did not provide additional insights in our dataset, likely because UTCI integrates multiple meteorological variables and therefore offers a more comprehensive and physiologically relevant measure of thermal stress in this context. The analysis of H index (Supplementary Table S3) nevertheless indicated a significantly increased risk of herpes zoster under very hot conditions, although overall its explanatory value was lower than that of the UTCI. Notably, days meeting the “very hot” H-class criterion were relatively infrequent in our series, which limits the precision of the risk estimate. This finding should therefore be interpreted with caution and corroborated in larger datasets and settings with more frequent extreme heat exposure.

The UTCI has been previously linked to infectious diseases. In Poland, lower UTCI values, indicating stronger cold stress, correlated with higher consultation rates for respiratory infections (32). This supports the hypothesis that thermal stress may influence immune function and susceptibility to viral reactivation, including varicella-zoster virus.

Extreme heat is a well-recognized hazard for public health, and its importance is growing in the context of global climate change (33). For instance, the severe European heatwave of 2003 was associated with around 70,000 excess deaths, predominantly among the elderly (39). Heatwaves not only cause direct heat-related illness and mortality, but they can also aggravate chronic conditions such as cardiovascular and respiratory diseases (40). With climate change, the frequency and intensity of heat extremes are rising; the proportion of heat-related deaths in people over 65 has increased significantly (by an estimated 85% between 2000 and 2021) (41).

Physiological stress from extreme temperatures can influence the immune system and has been implicated in the reactivation of other latent viruses (42). For example, heat stress has been observed to precipitate reactivation of herpes simplex virus (HSV-1) in susceptible individuals (43). A comprehensive review by Kennedy et al. (UK) compared HSV-1 and VZV, demonstrating that both viruses share key mechanisms of latency and reactivation, such as persistence in sensory ganglia and reliance on host immune control. This parallel supports the biological plausibility that similar triggers, such as thermal stress, may also contribute to VZV reactivation. In the case of VZV, our results suggest that elevated ambient temperatures might act as a trigger that disturbs the balance between the virus and the host’s nervous system, allowing the virus to emerge from latency (44). Figures 4, 5 similarly demonstrate that HZ occurrence increases with higher UTCI classes, underscoring a significant effect of thermal stress on VZV reactivation. Marcus et al. (USA) reported that in vitro, lowering the incubation temperature from 37 °C to 34 °C markedly increased the reactivation of VZV [36]. Notably, a skin surface temperature around 34 °C can occur under relatively mild heat stress (for example, on a warm day when a person is lightly clothed, corresponding to UTCI classes 7–8, 26–38 °C). Although more extreme heat stress conditions (beyond class 8) were not encountered in our region during the study period, the occurrence of high UTCI values – and the associated physiological heat strain – was clearly linked to increased herpes zoster cases in our analysis (see Figure 3).

Our findings in Figure 2 and Supplementary Table S1 indicate no meaningful relationship between ambient temperature alone and HZ cases, implying temperature by itself is not sufficient to trigger reactivation. The slightly lower number of cases observed in February (relative to some longer months like August or November) is likely explained by February’s shorter duration, and this difference was only on the borderline of statistical significance; thus, it should be interpreted with caution. Overall, an analysis based purely on calendar months or single weather variables could miss important nuances. In contrast, using the UTCI, which combines multiple meteorological factors, provided a clearer picture of risk: it seems that it is the composite effect of temperature together with humidity, wind, and solar radiation (and even behavioral factors like clothing insulation) that contributes to changes in HZ occurrence. UTCI thus appears to be a more effective predictor of HZ reactivation risk than temperature alone.

Our results are in line with several studies that found no strong seasonality in herpes zoster occurrence. Investigations in Australia and South Korea, for instance, reported no significant seasonal pattern in herpes zoster cases, whether examining hospitalizations or emergency department visits (45, 46). On the other hand, Wu et al. (Taiwan, China) noted that there was generally no pronounced seasonality, the number of cases did peak in the hottest month (August), (47). Notably, the climate in the heavily populated region of Taipei (subtropical, Köppen Cfa) features summer temperatures around 28 °C, which corresponds to a moderate heat stress level (approximately UTCI class 7) (31). This observation from Taiwan, China aligns with our finding that conditions of heat stress can increase HZ risk.

The lack of a traditional seasonal pattern in our data supports the notion that it is not the season or the calendar month per se, but rather the actual thermal stress from temperature extremes that drives the risk of disease. Focusing solely on calendar-defined seasons can be misleading, especially as climate change leads to atypical weather (warmer winters, unseasonal heatwaves) that disrupt usual seasonal norms (48). Our findings support the notion that HZ occurrence is not strictly seasonal but is influenced by meteorological conditions. It should be noted, however, that host factors—particularly age—remain fundamental determinants of risk (49). In our study, older adults had higher baseline occurrence of HZ (Figure 1). Even so, when considering our age-specific results alongside the UTCI analysis (Figures 4, 5), we see that thermal stress influences both younger and older adults, highlighting UTCI as a useful predictor of risk across age groups.

These results suggest that it may be worth incorporating complex climatic indices, like the UTCI, into public health planning for HZ prevention. If certain weather conditions indeed heighten the risk of zoster, health authorities could consider timing preventive measures (such as vaccination campaigns) to precede periods of high thermal stress. Vaccination is an effective strategy to reduce the risk of HZ reactivation, and aligning vaccine administration with forecasts of elevated UTCI values might enhance its protective impact (50). To date, vaccination programs are generally organized around fixed schedules or known seasonal patterns for other diseases (for example, annual influenza campaigns) (51, 52), and HZ immunization is primarily recommended based on age or immune status rather than climate considerations (53). This paradigm was evident during the initial rollout of COVID-19 vaccines, where prioritization was based on vulnerability and risk rather than environmental factors (54). Currently, weather conditions are rarely taken into account for vaccinations; HZ vaccines are administered irrespective of temperature or season, focusing mainly on the patient’s age (55). Our findings highlight a potential benefit in considering meteorological risk factors in HZ prevention strategies, especially in the face of changing climate patterns.

In this 2009–2023 study, we observed that days characterized by higher thermal stress (UTCI) coincided with increased herpes zoster risk, indicating that biometeorological conditions may influence VZV reactivation. These results align with evidence from Taiwan and Hefei, China, as well as South Korea, where warmer temperatures were linked to higher herpes zoster reporting or health-care utilization (11, 14, 17). Conversely, data from Shanghai reported winter predominance, highlighting that the direction and magnitude of temperature-related effects may vary by climatic zone, population characteristics, and exposure metrics (15). Notably, our use of the integrative UTCI adds physiological relevance compared with single-parameter indices. Taken together, the literature and our findings support the biological plausibility that heat-related stress can destabilize host–virus equilibrium and facilitate reactivation in susceptible individuals. Given the absence of patient-level confounders (e.g., vaccination, immunosuppression) and the use of regional exposure, these results should be considered hypothesis-generating and call for larger, multicenter studies with finer exposure assessment and richer clinical data.

Limitations

This study has several limitations. First, because the data were drawn from only two medical facilities, the findings might not be fully generalizable to the entire region or other populations. These two facilities were chosen because they provided complete and standardized electronic medical records for the entire study period, but this restriction, while ensuring high data quality, inevitably limits the generalizability of the findings to broader populations. Second, the meteorological conditions were assigned based on regional weather station data, which may not capture each individual’s exact exposure (microclimates or personal environmental conditions could differ). Third, as a retrospective analysis of medical records, the study relies on the accuracy and consistency of those records; any differences in diagnostic or reporting practices between the two facilities could affect the results. Additionally, we did not have data on certain personal factors such as patients’ immune status, co-morbid conditions, or stress levels, which could influence the likelihood of VZV reactivation. Moreover, vaccination status against varicella or herpes zoster was not available in our dataset and therefore could not be controlled for as a potential confounder. The shingles vaccine became available in Poland in March 2023, but its impact can be disregarded in our retrospective analysis since nationwide access only began on April 1, 2025, when it was added to the list of reimbursed drugs for individuals over 65 years by the National Health Fund (NFZ) (56). Another limitation is the use of regional meteorological data rather than individual-level exposure measures. Such an ecological approach may introduce bias, as personal behaviors (e.g., time spent outdoors, air conditioning use) and microclimatic variations could differ substantially across patients, leading to potential exposure misclassification. Another potential source of classification bias arises from the use of ICD-10 codes, which do not allow full distinction between recurrent or prolonged cases. To reduce this risk, repeated consultations occurring within 30 days were excluded from the analysis, but some residual misclassification cannot be entirely ruled out. UCTI may overestimate cold stress in populations that adapt through protective clothing or heating, while at the same time underestimating true exposure variability when based on single daily measurements. Future research should consider incorporating these clinical variables and possibly a broader range of locations to provide a more comprehensive understanding of how weather and other factors interact to affect herpes zoster risk.

Conclusion

Our findings indicate that thermal stress, particularly heat stress captured by the UTCI, may play a role in shaping herpes zoster occurrence, even in the absence of clear seasonality. This work adds to the emerging evidence linking climatic conditions with infectious diseases, offering a new perspective on herpes zoster epidemiology. While the study is exploratory and constrained by its retrospective design, limited geographic scope, and absence of certain individual-level variables, it highlights important avenues for future research. Expanding analyses to larger, multicenter cohorts and integrating data on vaccination, immune function, and behavioral exposures will be essential to confirm and refine these associations. Ultimately, a deeper understanding of the interplay between climate variability and herpes zoster risk may inform both preventive strategies and public health planning in the era of global climate change.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by Ethics Committee of the University of Warmia and Mazury in Olsztyn (approval no. WMIL-KB/90/2024). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

MO: Writing – review & editing, Formal analysis, Writing – original draft, Data curation, Methodology, Investigation, Conceptualization. KG: Validation, Writing – review & editing. ED: Data curation, Writing – review & editing. IC: Writing – review & editing, Data curation. JR: Writing – review & editing, Project administration, Formal analysis.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2025.1643828/full#supplementary-material

References

1. Gershon, AA, and Gershon, MD. Pathogenesis and current approaches to control of varicella-zoster virus infections. Clin Microbiol Rev. (2013) 26:728–43. doi: 10.1128/CMR.00052-13

PubMed Abstract | Crossref Full Text | Google Scholar

2. Heininger, U, and Seward, JF. Varicella. Lancet. (2006) 368:1365–76. doi: 10.1016/S0140-6736(06)69561-5

PubMed Abstract | Crossref Full Text | Google Scholar

3. Cohen, JI. Clinical practice: Herpes zoster. N Engl J Med. (2013) 369:255–63. doi: 10.1056/NEJMcp1302674

PubMed Abstract | Crossref Full Text | Google Scholar

4. Oxman, MN, and Levin, MJShingles Prevention Study G. Vaccination against herpes zoster and Postherpetic neuralgia. J Infect Dis. (2008) 197:S228–36. doi: 10.1086/522159

PubMed Abstract | Crossref Full Text | Google Scholar

5. Weinberg, A, and Levin, MJ. VZV T cell-mediated immunity. Curr Top Microbiol Immunol. (2010) 342:341–57. doi: 10.1007/82_2010_31

PubMed Abstract | Crossref Full Text | Google Scholar

6. Levin, MJ. Immune senescence and vaccines to prevent herpes zoster in older persons. Curr Opin Immunol. (2012) 24:494–500. doi: 10.1016/j.coi.2012.06.002

PubMed Abstract | Crossref Full Text | Google Scholar

7. Laing, KJ, Ouwendijk, WJD, Koelle, DM, and Verjans, G. Immunobiology of varicella-zoster virus infection. J Infect Dis. (2018) 218:S68–74. doi: 10.1093/infdis/jiy403

PubMed Abstract | Crossref Full Text | Google Scholar

8. Peng, Q, Guo, X, Luo, Y, Wang, G, Zhong, L, Zhu, J, et al. Dynamic immune landscape and VZV-specific T cell responses in patients with herpes zoster and Postherpetic neuralgia. Front Immunol. (2022) 13:887892. doi: 10.3389/fimmu.2022.887892

PubMed Abstract | Crossref Full Text | Google Scholar

9. Johnson, RW, and Rice, AS. Clinical practice. Postherpetic neuralgia. N Engl J Med. (2014) 371:1526–33. doi: 10.1056/NEJMcp1403062

PubMed Abstract | Crossref Full Text | Google Scholar

10. Baron, R, Binder, A, and Wasner, G. Neuropathic pain: diagnosis, pathophysiological mechanisms, and treatment. Lancet Neurol. (2010) 9:807–19. doi: 10.1016/S1474-4422(10)70143-5

PubMed Abstract | Crossref Full Text | Google Scholar

11. Lv, X, Fang, X, Qian, T, Cai, Y, Gao, P, Chen, H, et al. Association between meteorological factors and outpatient visits for herpes zoster in Hefei, China: a time-series analysis. Int J Environ Res Public Health. (2023) 20:2097. doi: 10.3390/ijerph20032097

PubMed Abstract | Crossref Full Text | Google Scholar

12. Lin, LY, Bhate, K, Forbes, H, Smeeth, L, Warren-Gash, C, and Langan, S. Vitamin D deficiency or supplementation and the risk of human herpesvirus infections or reactivation: a systematic review protocol. BMJ Open. (2019) 9:e031867. doi: 10.1136/bmjopen-2019-031867

PubMed Abstract | Crossref Full Text | Google Scholar

13. Berlinberg, EJ, Kim, E, Deiner, MS, Patterson, C, Porco, TC, and Acharya, NR. Seasonality of herpes zoster and herpes zoster ophthalmicus. J Clin Virol. (2020) 126:104306. doi: 10.1016/j.jcv.2020.104306

PubMed Abstract | Crossref Full Text | Google Scholar

14. Choi, YJ, Lim, YH, Lee, KS, and Hong, YC. Elevation of ambient temperature is associated with an increased risk of herpes zoster: a time-series analysis. Sci Rep. (2019) 9:12254. doi: 10.1038/s41598-019-48673-5

PubMed Abstract | Crossref Full Text | Google Scholar

15. Yang, Y, Chen, R, Xu, J, Li, Q, Xu, X, Ha, S, et al. The effects of ambient temperature on outpatient visits for varicella and herpes zoster in Shanghai, China: a time-series study. J Am Acad Dermatol. (2015) 73:660–5. doi: 10.1016/j.jaad.2015.07.015

PubMed Abstract | Crossref Full Text | Google Scholar

16. Kawai, K, VoPham, T, Drucker, A, Curhan, SG, and Curhan, GC. Ultraviolet radiation exposure and the risk of herpes zoster in three prospective cohort studies. Mayo Clin Proc. (2020) 95:283–92. doi: 10.1016/j.mayocp.2019.08.022

PubMed Abstract | Crossref Full Text | Google Scholar

17. Lai, SW, Liao, KF, Kuo, YH, Lin, CL, Liu, CS, Hwang, BF, et al. The impacts of ambient temperature and ultraviolet radiation on the incidence of herpes zoster: an ecological study in Taiwan. Int J Clin Pract. (2021) 75:e13854. doi: 10.1111/ijcp.13854

PubMed Abstract | Crossref Full Text | Google Scholar

18. Kawai, K, Yawn, BP, Wollan, P, and Harpaz, R. Increasing incidence of herpes zoster over a 60-year period from a population-based study. Clin Infect Dis. (2016) 63:221–6. doi: 10.1093/cid/ciw296

PubMed Abstract | Crossref Full Text | Google Scholar

19. Graczyk, D, Kundzewicz, ZW, Choryński, A, Førland, EJ, Pińskwar, I, and Szwed, M. Heat-related mortality during hot summers in polish cities. Theor Appl Climatol. (2019) 136:1259–73. doi: 10.1007/s00704-018-2554-x

Crossref Full Text | Google Scholar

20. Pińskwar, I, Choryński, A, Graczyk, D, and Kundzewicz, ZW. Observed changes in extreme precipitation in Poland: 1991–2015 versus 1961–1990. Theor Appl Climatol. (2019) 135:773–87. doi: 10.1007/s00704-018-2372-1

Crossref Full Text | Google Scholar

21. Kubiak-Wójcicka, K, and Bąk, B. Monitoring of meteorological and hydrological droughts in the Vistula basin (Poland). Environ Monit Assess. (2018) 190:691. doi: 10.1007/s10661-018-7058-8

PubMed Abstract | Crossref Full Text | Google Scholar

22. Kawai, K, and Yawn, BP. Risk factors for herpes zoster: a systematic review and Meta-analysis. Mayo Clin Proc. (2017) 92:1806–21. doi: 10.1016/j.mayocp.2017.10.009

PubMed Abstract | Crossref Full Text | Google Scholar

23. Weitzman, D, Shavit, O, Stein, M, Cohen, R, Chodick, G, and Shalev, V. A population based study of the epidemiology of herpes zoster and its complications. J Infect. (2013) 67:463–9. doi: 10.1016/j.jinf.2013.06.016

PubMed Abstract | Crossref Full Text | Google Scholar

24. Blazejczyk, K, Epstein, Y, Jendritzky, G, Staiger, H, and Tinz, B. Comparison of UTCI to selected thermal indices. Int J Biometeorol. (2012) 56:515–35. doi: 10.1007/s00484-011-0453-2

PubMed Abstract | Crossref Full Text | Google Scholar

25. Malcheva, K, and Gocheva, A. Thermal comfort indices for the cold half-year in Sofia. Bulg J Meteorol Hydrol. (2014) 19:16–25. doi: 10.13140/RG.2.1.1769.4325

Crossref Full Text | Google Scholar

26. National Fund Health. Medical assistance at night, weekends and holidays. (2023). Available online at: https://ippez.pl/en/medical-assistance-at-night-weekends-and-holidays/. (Accessed September 01, 2025)

Google Scholar

27. Poland, S. Population. Size and structure and vital statistics in Poland by territorial division in 2023 (2025). Available online at: https://stat.gov.pl/en/topics/population/population/population-size-and-structure-and-vital-statistics-in-poland-by-territorial-division-in-2023-as-of-31-december,3,35.html. (Accessed September 01, 2025).

Google Scholar

28. Sundstrom, K, Weibull, CE, Soderberg-Lofdal, K, Bergstrom, T, Sparen, P, and Arnheim-Dahlstrom, L. Incidence of herpes zoster and associated events including stroke--a population-based cohort study. BMC Infect Dis. (2015) 15:488. doi: 10.1186/s12879-015-1170-y

PubMed Abstract | Crossref Full Text | Google Scholar

29. Kieloch, B, Kamińska, M, Lech, H, Górka, A, Grzelak-Agaciak, E, Ragin, G, et al. Instrukcja dla stacji meteorologicznych. Warszawa: Instytut Meteorologii i Gospodarki Wodnej – Państwowy Instytut Badawczy (IMGW-PIB) (2015). 442 p.

Google Scholar

30. World Meteorological Organization. Guide to meteorological instruments and methods of observation. WMO-No. 8th ed. Geneva: World Meteorological Organization (2017).

Google Scholar

31. Peel, MC, Finlayson, BL, and McMahon, TA. Updated world map of the Köppen-Geiger climate classification. Hydrol Earth Syst Sci. (2007) 11:1633–44. doi: 10.5194/hess-11-1633-2007

Crossref Full Text | Google Scholar

32. Romaszko, J, Skutecki, R, Bocheński, M, Cymes, I, Dragańska, E, Jastrzębski, P, et al. Applicability of the Universal Thermal Climate Index for predicting the outbreaks of respiratory tract infections: a mathematical modeling approach. Int J Biometeorol. (2019) 63:1231–41. doi: 10.1007/s00484-019-01740-y

PubMed Abstract | Crossref Full Text | Google Scholar

33. Di Napoli, C, Pappenberger, F, and Cloke, HL. Assessing heat-related health risk in Europe via the Universal Thermal Climate Index (UTCI). Int J Biometeorol. (2018) 62:1155–65. doi: 10.1007/s00484-018-1518-2

PubMed Abstract | Crossref Full Text | Google Scholar

34. Romaszko, J, Draganska, E, Jalali, R, Cymes, I, and Glinska-Lewczuk, K. Universal climate thermal index as a prognostic tool in medical science in the context of climate change: a systematic review. Sci Total Environ. (2022) 828:154492. doi: 10.1016/j.scitotenv.2022.154492

PubMed Abstract | Crossref Full Text | Google Scholar

35. Bröde, P, Fiala, D, Blazejczyk, K, Holmér, I, Jendritzky, G, Kampmann, B, et al. Deriving the operational procedure for the Universal Thermal Climate Index (UTCI). Int J Biometeorol. (2011) 56:481–94. doi: 10.1007/s00484-011-0454-1

PubMed Abstract | Crossref Full Text | Google Scholar

36. Błażejczyk, K, and Bröde, P. Principles of the new Universal Thermal Climate Index (UTCI) and its application to bioclimatic research in European scale. Miscellanea Geogr. (2010) 14:91–102. doi: 10.2478/mgrsd-2010-0009

Crossref Full Text | Google Scholar

37. Lumivero ed. XLSTAT. 1st ed Lumivero. Denver: Lumivero, (2025).

Google Scholar

38. Ghada, W, Estrella, N, Ankerst, DP, and Menzel, A. Universal Thermal Climate Index associations with mortality, hospital admissions, and road accidents in Bavaria. PLoS One. (2021) 16:e0259086. doi: 10.1371/journal.pone.0259086

PubMed Abstract | Crossref Full Text | Google Scholar

39. Bouchama, A. The 2003 European heat wave. Intensive Care Med. (2004) 30:1–3. doi: 10.1007/s00134-003-2062-y

PubMed Abstract | Crossref Full Text | Google Scholar

40. Ponjoan, A, Blanch, J, Alves-Cabratosa, L, Marti-Lluch, R, Comas-Cufi, M, Parramon, D, et al. Effects of extreme temperatures on cardiovascular emergency hospitalizations in a Mediterranean region: a self-controlled case series study. Environ Health. (2017) 16:32. doi: 10.1186/s12940-017-0238-0

PubMed Abstract | Crossref Full Text | Google Scholar

41. Climate Change, Heat and Health: World Health Organization (2024). Available online at: https://www.who.int/news-room/fact-sheets/detail/climate-change-heat-and-health. (Accessed September 01, 2025).

Google Scholar

42. Huang, W, Xie, P, Xu, M, Li, P, and Zao, G. The influence of stress factors on the reactivation of latent herpes simplex virus type 1 in infected mice. Cell Biochem Biophys. (2011) 61:115–22. doi: 10.1007/s12013-011-9167-7

PubMed Abstract | Crossref Full Text | Google Scholar

43. Jones, C. Intimate relationship between stress and human alpha-herpes virus 1 (HSV-1) reactivation from latency. Curr Clin Microbiol Rep. (2023) 10:236–45. doi: 10.1007/s40588-023-00202-9

PubMed Abstract | Crossref Full Text | Google Scholar

44. Thellman, NM, and Triezenberg, SJ. Herpes simplex virus establishment, maintenance, and reactivation: in vitro modeling of latency. Pathogens. (2017) 6:28. doi: 10.3390/pathogens6030028

PubMed Abstract | Crossref Full Text | Google Scholar

45. Jardine, A, Conaty, SJ, and Vally, H. Herpes zoster in Australia: evidence of increase in incidence in adults attributable to varicella immunization? Epidemiol Infect. (2011) 139:658–65. doi: 10.1017/S0950268810001949

PubMed Abstract | Crossref Full Text | Google Scholar

46. Kim, YJ, Lee, CN, Lim, CY, Jeon, WS, and Park, YM. Population-based study of the epidemiology of herpes zoster in Korea. J Korean Med Sci. (2014) 29:1706–10. doi: 10.3346/jkms.2014.29.12.1706

PubMed Abstract | Crossref Full Text | Google Scholar

47. Wu, PY, Wu, HD, Chou, TC, and Sung, FC. Varicella vaccination alters the chronological trends of herpes zoster and varicella. PLoS One. (2013) 8:e77709. doi: 10.1371/journal.pone.0077709

PubMed Abstract | Crossref Full Text | Google Scholar

48. Agency USEP. Climate Change Indicators: Seasonal Temperature. (2023). Available online at: https://www.epa.gov/climate-indicators/climate-change-indicators-seasonal-temperature. (Accessed September 01, 2025).

Google Scholar

49. Steinmann, M, Lampe, D, Grosser, J, Schmidt, J, Hohoff, ML, Fischer, A, et al. Risk factors for herpes zoster infections: a systematic review and meta-analysis unveiling common trends and heterogeneity patterns. Infection. (2024) 52:1009–26. doi: 10.1007/s15010-023-02156-y

PubMed Abstract | Crossref Full Text | Google Scholar

50. Vogelsang, EM, and Polonijo, AN. Social determinants of shingles vaccination in the United States. J Gerontol B Psychol Sci Soc Sci. (2022) 77:407–12. doi: 10.1093/geronb/gbab074

PubMed Abstract | Crossref Full Text | Google Scholar

51. Weinberger, B. Vaccination of older adults: influenza, pneumococcal disease, herpes zoster, COVID-19 and beyond. Immun Ageing. (2021) 18:38. doi: 10.1186/s12979-021-00249-6

PubMed Abstract | Crossref Full Text | Google Scholar

52. Rogers, CJ, Bahr, KO, and Benjamin, SM. Attitudes and barriers associated with seasonal influenza vaccination uptake among public health students; a cross-sectional study. BMC Public Health. (2018) 18:1131. doi: 10.1186/s12889-018-6041-1

PubMed Abstract | Crossref Full Text | Google Scholar

53. Bradley, S. Vaccinations for older adults: an update. Curr Geriatr Rep. (2014) 3:66–72. doi: 10.1007/s13670-014-0076-7

Crossref Full Text | Google Scholar

54. Grana, C, Ghosn, L, Evrenoglou, T, Jarde, A, Minozzi, S, Bergman, H, et al. Efficacy and safety of COVID-19 vaccines. Cochrane Database Syst Rev. (2022) 12:CD015477. doi: 10.1002/14651858.CD015477

PubMed Abstract | Crossref Full Text | Google Scholar

55. Cunningham, AL, Lal, H, Kovac, M, Chlibek, R, Hwang, SJ, Diez-Domingo, J, et al. Efficacy of the herpes zoster subunit vaccine in adults 70 years of age or older. N Engl J Med. (2016) 375:1019–32. doi: 10.1056/NEJMoa1603800

PubMed Abstract | Crossref Full Text | Google Scholar

56. PZH–PIB (National Institute of Public Health) NIoPHN (National Institute of Hygiene). Ministerstwo Zdrowia obejmuje refundacją kolejne szczepionki (The Ministry of Health extends reimbursement to additional vaccines). (2025). Available online at: https://szczepienia.pzh.gov.pl/ministerstwo-zdrowia-obejmuje-refundacja-kolejne-szczepionki/. (Accessed September 01, 2025).

Google Scholar

Keywords: herpes zoster, climate, Universal Thermal Climate Index (UTCI), temperature, weather, epidemiology

Citation: Ochal M, Glińska Lewczuk K, Dragańska E, Cymes I and Romaszko J (2025) Influence of meteorological conditions on herpes zoster occurrence: a retrospective cohort study. Front. Med. 12:1643828. doi: 10.3389/fmed.2025.1643828

Received: 09 June 2025; Accepted: 25 September 2025;
Published: 10 October 2025.

Edited by:

Phillip R. Kramer, Texas A&M University, United States

Reviewed by:

Mate Magyar, Semmelweis University, Hungary
Carlos Cifuentes-González, Tan Tock Seng Hospital, Singapore
Pilar Cárdenas-Soriano, Hospital Universitario Fundación Alcorcón, Spain

Copyright © 2025 Ochal, Glińska Lewczuk, Dragańska, Cymes and Romaszko. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Michał Ochal, b2NoYWxtaWNoYWw3QGdtYWlsLmNvbQ==

ORCID: Michał Ochal, orcid.org/0000-0003-0484-6778
Katarzyna Glińska-Lewczuk, orcid.org/0000-0002-7845-6353
Ewa Dragańska, orcid.org/0000-0002-8099-6359
Iwona Cymes, orcid.org/0000-0001-6055-8266
Jerzy Romaszko, orcid.org/0000-0002-7929-6822

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.