Next Article in Journal
Papillomavirus-like Particles in Equine Medicine
Next Article in Special Issue
Pending Reorganization of Hantaviridae to Include Only Completely Sequenced Viruses: A Call to Action
Previous Article in Journal
Optimal Conditions for In Vitro Assembly of Respiratory Syncytial Virus Nucleocapsid-like Particles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seroprevalence of Hantavirus in Forestry Workers, Northern France, 2019–2020

by
Catarina Krug
1,2,
Emma Rigaud
3,
Dieyenaba Siby-Diakite
4,
Laetitia Bénézet
1,
Pavlos Papadopoulos
4,
Henriette de Valk
1,
Gaëtan Deffontaines
3,
Alexandra Septfons
1,† and
Jean-Marc Reynes
4,*,†
1
Santé Publique France, 94410 Saint-Maurice, France
2
European Centre for Disease Prevention and Control (ECDC), 169 73 Solna, Sweden
3
Caisse Centrale de la Mutualité Sociale Agricole, 93000 Bobigny, France
4
Institut Pasteur, Université Paris Cité, 75015 Paris, France
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Viruses 2023, 15(2), 338; https://doi.org/10.3390/v15020338
Submission received: 16 December 2022 / Revised: 16 January 2023 / Accepted: 21 January 2023 / Published: 25 January 2023

Abstract

:
We aimed to estimate the seroprevalence of Puumala orthohantavirus (PUUV) among forestry workers in northern France, and to explore sociodemographic risk factors. We conducted a random cross-sectional seroprevalence survey among 1777 forestry workers in 2019–2020. The presence of immunoglobulin G against PUUV antigens in serum was assessed using enzyme-linked immunosorbent assay and confirmed using immunofluorescence assay. Poisson regression models were used to explore factors associated with seropositivity. Weighted seroprevalence was 5% (3–6) in northeastern France, 4% (2–6) in north central France, and 1% in two regions located in the center of the country (Auvergne and Limousin). There were no seropositive workers detected in northwestern France. Seropositivity was associated with age, sex, and cumulative seniority in the forestry sector. Seroprevalence was highest in known endemic areas of the northeast and lowest in the northwest. Nevertheless, we found serological evidence of PUUV infection in two regions located in the center of the country, suggesting circulation of the virus in these regions, previously thought to be non-endemic.

1. Introduction

Zoonotic hantaviruses are rodent-borne enveloped viruses with a tri-segmented, negative sense, single-stranded RNA genome, belonging to the Orthohantavirus genus (Subfamily Mammantavirinae, Family Hantaviridae, Order Bunyavirales). They are responsible in Eurasia for mild to severe hemorrhagic fever with renal syndrome (HFRS) [1].
Puumala orthohantavirus (PUUV) is the hantavirus species most commonly detected in patients in metropolitan France, whereas Seoul orthohantavirus (SEOV) and Tula orthohantavirus (TULV) are detected sporadically. As in some other European countries [1,2], the spatial distribution of HFRS cases is not homogenous in metropolitan France, with detection of human cases being restricted to the northeastern area, including the highly endemic regions of Nord-Pas-de-Calais, Picardie, Franche-Comté [3,4]. There were less than ten individuals with anti-hantavirus antibodies (probable SEOV cases) reported outside the northeastern area in the nineties [5]. The main risk factors for human PUUV infection include living near the forest, the presence of rodents in the dwellings, and leisure or professional activities in the forest [6].
In contrast to the spatial distribution of HFRS cases, the main natural hosts of these three hantavirus species are present in most regions of metropolitan France. Hantavirus hosts include the bank vole (Myodes glareolus) for PUUV, the brown rat (Rattus norvegicus) for SEOV, and the common vole (Microtus arvalis) for TULV. There are various hypotheses to explain this heterogeneous geographic distribution of PUUV, including environmental factors (influencing the rodent population dynamics and virus survival outside the host), rodent susceptibility (including virus replication and excretion), genetic characteristics of the viruses (favoring or not, its replication, its excretion, its survival outside the host), and human behaviors (contact with rodents or their excreta) [2,3,7]. PUUV prevalence studies in bank vole populations have only been conducted in the northeastern quarter of metropolitan France, in the endemic area, and its vicinity [8,9,10,11,12,13,14,15,16,17]. Bank voles have been found infected in these areas, with varying prevalence level of anti-PUUV antibodies according to time and space. There is no information outside this area.
Cases may be missed outside the northeastern quarter of metropolitan France because patients are not frequently tested for hantavirus [18]. The occurrence of the disease is not expected by clinicians and awareness of hantavirus diseases is likely to be low among clinicians in non-endemic areas [18].
There is little information on the presence of PUUV in French regions that are considered non-endemic. Given that forestry workers are at higher risk than the rest of the population due to potential exposure to rodents, particularly the bank vole, host of PUUV, and their excreta [19], conducting a study targeting forestry workers can help us document the distribution of hantavirus. Our primary aim was to estimate the overall and regional prevalence of antibodies against PUUV antigens among forestry workers in northern France. Our secondary aim was to compare seroprevalence between geographic areas, particularly northeast (considered to be endemic) and other northern areas (considered to be non-endemic); and to determine sociodemographic factors associated with hantavirus seropositivity.

2. Materials and Methods

2.1. Sample Size Calculations

The original study was designed to investigate the seroprevalence of Borrelia burgdorferi sl. The sample size for our study was therefore determined to answer this research objective. The number of required subjects was calculated for each of the following five geographic areas: northwestern France including Normandie (Haute-Normandie and Basse-Normandie), Bretagne and Pays-de-la-Loire regions; north central France including Hauts-de-France (Picardie, Nord-Pas-de-Calais), Île-de-France, and Centre regions; northeastern France including Alsace, Lorraine, Champagne-Ardenne, Bourgogne, and Franche-Comté regions; and Auvergne and Limousin regions, both located in the center of France. Assuming a B. burgdorferi sl seroprevalence of 5% in northwestern and north central France and of 15% in northeastern France [20,21,22], a sample of 2591 individuals would be required to achieve a minimum precision of 1.5% and a maximum precision of 4% in the estimation of B. burgdorferi sl seroprevalence. A total of 1778 forestry workers were included from May 2019 to March 2020. The number of participants was lower than required due to the COVID-19 pandemic as the lockdown led to the interruption of data collection. We conducted power calculations to estimate the precision of hantavirus seroprevalence estimation using the available sample (see Table S1).

2.2. Study Design, Participants, and Samples

We conducted a random cross-sectional seroprevalence survey in northern France from May 2019 to March 2020. The target population consisted of all forestry workers aged 18 or older, working at least once a week in a wooded environment or park covered by the French farmers’ insurance fund (Mutualité Sociale Agricole, MSA) of one of the following administrative regions: Normandie (Haute-Normandie and Basse-Normandie), Bretagne, Pays-de-la-Loire, Hauts-de-France (Picardie, Nord-Pas-de-Calais), Île-de-France, Centre, Alsace, Lorraine, Champagne-Ardenne, Bourgogne, Franche-Comté, Auvergne, and Limousin.
To recruit study participants, we used a stratified random sample. Strata were defined by region (as defined by affiliation to the local MSA insurance fund), occupational activity (silviculturist, forestry technician, woodcutter, etc.), and working status (self-employed, salaried, other). Within each stratum, simple random sampling was used to select individuals. The number of individuals sampled was proportional to the size of the strata.
Individuals were invited for a medical visit by their occupational health physician. After obtaining informed consent, participants were interviewed, and blood samples were collected. The questionnaire administered included sociodemographic information such as age and region of MSA insurance fund affiliation, as well as details on occupational activity including main profession, average weekly exposure to the forest, and cumulative seniority in forestry activities. The questionnaire also included questions on the main departments where forestry activities occurred during the last 12 months, as well as secondary departments or countries. All blood samples were anonymized at the time of collection by giving each a unique number. Blood samples were kept at room temperature (10–30 °C) and centrifuged (3200 rpm for 10 min, at 20–24 °C) within 48 h. Sera were frozen at −20 °C within 2 h after centrifugation. Serum was divided into five aliquots and stored at −30 °C until analyzed.

2.3. Serological Methods

Collected sera were screened for the presence of anti-PUUV immunoglobulin G (IgG) at 1:100 dilution using homemade enzyme-linked immunosorbent assay (ELISA) according to Rossi and colleagues’ method [23]. ELISA-positive sera were screened at 1:64 dilution using a homemade immunofluorescence assay (IFA) according to Niklasson and colleagues’ method [24], to confirm the presence of anti-PUUV IgG. Both assays were conducted under ISO 15,189 accreditation. The Hantavirus National Reference Center has estimated a sensitivity of 100% and a specificity of 94.8% for its ELISA and of 100% and 98.1% for its IFA for serum samples from patients taken 10 days after symptom onset (Jean-Marc Reynes, personal communication). A negative result indicated a lack of PUUV exposure while a positive result indicated hantavirus exposure but not necessarily specific to PUUV, due to serological cross-reaction with other hantaviruses [25].

2.4. Statistical Analysis

All statistical analyses were weighted. The weights used considered the sampling weight and were adjusted for non-response. Non-response was corrected by reweighting using the equal-quantile score method [26], based on socio-demographic and geographic data and the professional activity sector available in the sampling frame for all individuals. A calibration by raking ratio method was then applied, using the distributions by sex, age group, and professional activity sector in the target population, using the SAS macro Calmar [27].
We described population demographic characteristics with absolute and relative frequencies for categorical variables and with mean and standard deviation for continuous variables. We explored factors associated with hantavirus seropositivity by estimating seroprevalence ratios and 95% confidence intervals (CI) using weighted Poisson regression models with robust standard errors. The clustering of observations by geographic area was accounted for using a fixed area effect in all models. The outcome of interest was hantavirus seropositivity, and we used one model per factor of interest (the main exposure variable). Causal diagrams were made between hantavirus seropositivity and each factor of interest using DAGitty ([28,29]; diagrams not shown). Based on those causal diagrams, it was deemed reasonable to use unconditional models for the geographic area, age, and main profession (i.e., no important confounders were identified for those factors). Also, we identified age as a confounder of the association between seniority and hantavirus seropositivity; and main profession as a confounder of the association between average weekly exposure time in forest and hantavirus seropositivity. Nevertheless, the strong collinearity between age and seniority (Pearson correlation coefficient of 0.76), precluded a multivariable analysis testing both factors. Analyses were performed using Stata 14.

3. Results

Out of 4484 forestry workers randomly selected, 1778 participated in this study (response rate of 41% [30], and serological test results were available for 1777. Demographic characteristics of the study population are described in Table 1. The mean (range) age of participants was 44 (18–87) years, 94% (1714/1777) were men, and one-third of participants were affiliated to three French farmers’ insurance funds regions in northeastern France: Lorraine (13%, 202/1777), Champagne-Ardenne (12%, 148/1777) and Alsace (10%, 133/1777). As for the characteristics of the work performed, the majority of participants were working in one of four occupational activities: woodcutter (34%, 622/1768), forest technician or ranger (21%, 357/1768), forestry machine operator (18%, 319/1768), and silviculturist (17%, 289/1768). Seventy-eight percent (1399/1758) of participants were exposed to the forest, forest borders or edges, parks, and wooded gardens for longer than 20 h a week (Table 1).
Out of 1777 individuals, 50 (2.8%) had IgG against PUUV antigens in their serum. We estimated a weighted seroprevalence of 3.3% (95% CI: 2.5–4.3). Seroprevalence varied across age groups, being lowest among individuals aged 18–40 years (1%, 95% CI: 0.73–3) and highest among those older than 60 years (6%, 95% CI: 3–11). Specifically, seroprevalence was four times higher (95% CI: 1.5–10.2) among individuals older than 60 years than among those aged 18–40 years (Table 2). All workers with anti-PUUV IgG were male. In Figure 1, we observe a seroprevalence increasing from northwest to northeast France, with no seropositive workers detected in northwestern France, a low seroprevalence in Centre (1%, 95% CI: 0.25–3), Limousin (1%, 95% CI: 0.09–3) and Auvergne (1%, 95% CI: 0.33–4) regions, and highest seroprevalence in Hauts-de-France (7%, 95% CI: 4–12) and Franche-Comté regions (8%, 95% CI: 5–14; Table 2). Workers with anti-PUUV IgG in Auvergne and Centre reported having been only exposed in their own region, with the exception of one worker in Auvergne who also traveled to Limousin and Centre. The worker with anti-PUUV IgG in Limousin was not exposed in other French regions, but traveled to Eastern Europe/Balkans, hantavirus endemic areas. There was no association between traveling to neighboring countries, and having anti-PUUV IgG.
Seroprevalence appeared to be greater among forest technicians and forest rangers (6%, 95% CI: 4–9) than among woodcutters (3%, 95% CI: 2–5) or silviculturists (2%, 95% CI: 0.63–4), but differences were not statistically significant (Table 2). Seroprevalence increased with average weekly exposure time to the forest, forest borders or edges, parks, and wooded gardens, with a seroprevalence of 1% (95% CI: 0.32–4) among workers exposed for less than 10 hours a week and 4% (95% CI: 3–5) among workers exposed for longer than 20 h a week. Corresponding to a prevalence ratio of four (95% CI: 1.2–12.5) between the two categories after adjusting for the main profession. Seroprevalence was lower among forestry workers with one to ten years of work in the sector (1%, 95% CI: 0.50–2) than among those with over 20 years of work in the sector (5%, 95% CI: 3–8). Specifically, seroprevalence was five times higher (95% CI: 1.9–12.5) among individuals with 21–30 years of work than among those with 1–10 years of work in the sector (Table 2). There was no association between hunting, fishing, hiking, or harvesting, and having anti-PUUV IgG.

4. Discussion

Our primary aim was to estimate the overall and regional seroprevalence of the hantavirus among forestry workers in northern France. In our study, the overall seroprevalence for hantavirus was estimated at 3.3% (95% CI: 2.5–4.3), which aligns with a previous meta-analysis by Riccò and colleagues [19] that indicated a seroprevalence of 4.1% (95% CI: 2.7–6.1) among European forestry workers. Nevertheless, results may not be directly comparable to ours, because this metanalysis included data from various European countries (west, south, and east European countries) with heterogeneous hantavirus incidences, and a variety of laboratory methodologies.
Our secondary aim was to compare seroprevalence between geographic areas, particularly northeast and other northern areas. As anticipated, the seroprevalence among forestry workers was greatest in known endemic areas of northeastern and north central France, particularly, in the highly endemic regions of Franche-Comté and Hauts-de-France [18]. Champagne-Ardenne region is also a known endemic region. We observed lower seroprevalence in the Champagne-Ardenne region (4%) than a previous study conducted in the Ardennes (part of the Champagne-Ardenne region) by Penalba and colleagues [31], who described a 15% seroprevalence in 55 forestry workers [31]. No details on the laboratory methods used in their study were described, which makes it difficult to compare with our results. Île-de-France, an endemic area with a lower number of cases than other endemic areas, had a lower seroprevalence of 3%, in accordance with a previous study conducted in 1991–1992 among forestry workers [32].
We did not find evidence of the circulation of PUUV in northwestern France (Normandie, Bretagne, and Pays-de-la-Loire), which is considered non-endemic [18]. Specifically, given that there are around 1000 forestry workers in northwestern France, and that we found zero hantavirus workers with anti-PUUV IgG among the 396 tested, we estimate at a 95% confidence level that the real seroprevalence of hantavirus in that area is below 1.3% [33]. This is in agreement with the absence of detection of hantavirus cases by healthcare practitioners [18]. In contrast, we found evidence of hantavirus infection (three workers with anti-PUUV IgG) in the Auvergne and Limousin regions, located in the center of the country and previously considered to be non-endemic [18]. A landscape and regional environmental analysis of the spatial distribution of hantavirus human cases in Europe showed that HFRS cases should be expected in these areas [2]. Similarly, a previous study performed in the 80′s in the Auvergne region (located in the center of France) had 2.3% of more than 1000 farmers testing positive for IgG against Hantaan orthohantavirus or PUUV antigens by immunofluorescence (serum dilution starting at 1:16) [34]. One of the limitations of our study was that we were unable to estimate hantavirus seroprevalence across the whole of France. Now that our study found evidence of hantavirus infection in workers in an area previously considered to be non-endemic, we recommend closer monitoring of southern France regions, as the area of virus circulation could very well extend. Nevertheless, we cannot dismiss the hypothesis that the individuals with anti-PUUV IgG in Auvergne and Limousin did not recall being exposed to endemic regions or countries (recall bias).
Finally, we also aimed to determine sociodemographic factors associated with the prevalence of anti-PUUV IgG. Age was positively associated with the presence of anti-PUUV IgG similar to the results reported in other studies [35,36]. This is explained by the fact that exposure to hantavirus accumulates over the course of a lifetime, due to repeated exposure to rodents infected excreta, and by the long half-life of these antibodies [37]. Similar to the results reported in other studies [35,36], seroprevalence also increased with longer exposure time per week in the forest, forest borders or edges, parks, and wooded gardens and with cumulative seniority in the forestry sector. The proportion of women in our target population was low, which affected our ability to draw conclusions on the relationship between sex at birth and the presence of anti-PUUV IgG.
Forestry workers are at higher risk of interacting with rodent hosts than the general population. Thus, our results cannot be generalized to the overall French population. We can also hypothesize that forestry workers are more aware than the general French population of occupational-related infections and their preventive measures. Indeed, there are specific HFRS prevention campaigns targeting forestry workers, as it is part of the list of occupational diseases in France [38]. Moreover, a previous French study indicated that forestry workers had greater knowledge about tick-borne diseases and were more likely to protect themselves than the general French population [39], indicating that they do follow recommendations and have a better knowledge of occupational diseases than the general population. Since the original study was aimed at investigating the seroprevalence of B. burgdorferi sl, we did not collect information on the use of preventive measures for hantavirus infections (such as the use of masks while working), and we could not evaluate their relationship with hantavirus seroprevalence.
In our study, we limited selection bias by randomly selecting forestry workers and by adjusting for non-response according to differences in sociodemographic and professional characteristics. Nevertheless, we cannot exclude the possibility of residual selection bias due to differences in hantavirus exposure between participants and non-participants with the same socio-demographic and professional characteristics.

5. Conclusions

In conclusion, the regional seroprevalence of hantavirus antibodies among forestry workers was greatest in the known endemic areas of northeastern and north central France, and lowest in the non-endemic northwestern France. Nevertheless, we found indications of the circulation of hantavirus in the Auvergne and Limousin regions, located in the center of the country, even though they were previously thought to be non-endemic. These results may indicate some risk of local HFRS cases in these regions, and the need to raise awareness among general and hospital practitioners and among at-risk populations. Practitioners in non-endemic areas should still be aware of the possibility of HFRS cases, especially among patients returning from endemic areas. An improved understanding of hantavirus infection rates in reservoir host species and virus transmission to humans is still needed to fully understand the absence of infection in non-endemic areas. Rodents’ ecology is highly variable, not only at the geographical level, but also over time, due to complicated interactions with their environment [40].
Routine surveillance has the disadvantage that cases might be missed in newly endemic areas that are still thought to be non-endemic [18]. The information obtained from routine surveillance may thus be complemented by sporadic zero surveillance, a useful tool to assess the cumulative exposure to hantavirus and the geographic distribution of human exposure. A reinforcement of sensitization campaigns among occupational health workers on hantavirus infection could help early identification of cases in endemic and non-endemic areas.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15020338/s1, Table S1: Precision of hantavirus seroprevalence estimates, assuming a sample size of approximately 1781 individuals and a true seroprevalence lying between 0.5 and 10%.

Author Contributions

Conceptualization, C.K., H.d.V., A.S. and J.-M.R.; methodology, C.K., L.B. and A.S.; software, C.K., L.B. and A.S.; validation, H.d.V., A.S. and J.-M.R.; formal analysis, C.K.; investigation, C.K., E.R., D.S.-D., P.P., G.D., A.S. and J.-M.R.; resources, E.R., D.S.-D., L.B., P.P. and G.D.; data curation, C.K., L.B. and A.S.; writing—original draft preparation, C.K.; writing—review and editing, all authors; visualization, C.K.; supervision, A.S. and J.-M.R.; project administration, A.S.; funding acquisition, C.K., E.R., G.D., A.S. and J.-M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Caisse Centrale de la Mutualité Sociale Agricole, the French Agency for Food, Environmental and Occupational Health & Safety (ANSES), and Santé Publique France. The National Reference Center for hantavirus was financially supported by Santé Publique France.

Institutional Review Board Statement

The study was approved by the Committee for the Protection of Persons participating in medical research (CPP Ouest II, Angers; reference number: 18-09-13-58557 2018/66) and by the French national data protection authority (Commission Nationale de l’Informatique et des Libertés, DR-2020-011).

Informed Consent Statement

Informed consent was obtained from all study participants.

Data Availability Statement

Data is available upon request.

Acknowledgments

We thank the forestry workers who accepted to participate in this study, as well as occupational physicians and nurses.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Faber, M.; Krüger, D.H.; Auste, B.; Stark, K.; Hofmann, J.; Weiss, S. Molecular and epidemiological characteristics of human Puumala and Dobrava-Belgrade hantavirus infections, Germany, 2001 to 2017. Eurosurveillance 2019, 24, 180–675. [Google Scholar]
  2. Zeimes, C.B.; Quoilin, S.; Henttonen, H.; Lyytikäinen, O.; Vapalahti, O.; Reynes, J.-M.; Reusken, C.; Swart, A.N.; Vainio, K.; Hjertqvist, M.; et al. Landscape and regional environmental analysis of the spatial distribution of hantavirus human cases in Europe. Front. Public Health 2015, 3, 54. [Google Scholar] [PubMed] [Green Version]
  3. Reynes, J.-M.; Penalba, C.; Galempoix, J.-M. Hantavirus. EMC Mal. Infect. 2022, in press. [Google Scholar]
  4. Brun, A.; Greusard, M.; Reynes, J.-M.; Grenier, M.; Bamoulid, J.; Giraudoux, P.; Lepiller, Q.; Chirouze, C.; Bouiller, K.; Bailly, B. Description of an outbreak of hemorrhagic fever with renal syndrome in the southern Jura Mountains, France, in 2021. Infect. Dis. Now 2023. [Google Scholar] [CrossRef] [PubMed]
  5. Le Guenno, B. Les hantavirus. Med. Mal. Infect. 1997, 27, 703–710. [Google Scholar] [CrossRef]
  6. Crowcroft, N.; Infuso, A.; Ilef, D.; Le Guenno, B.; Desenclos, J.; Van Loock, F.; Clement, J. Risk factors for human hantavirus infection: Franco-Belgian collaborative case-control study during 1995-6 epidemic. BMJ 1999, 318, 1737–1738. [Google Scholar] [CrossRef] [Green Version]
  7. Madrières, S.; Tatard, C.; Murri, S.; Vulin, J.; Galan, M.; Piry, S.; Pulido, C.; Loiseau, A.; Artige, E.; Benoit, L.; et al. How Bank Vole-PUUV Interactions Influence the Eco-Evolutionary Processes Driving Nephropathia Epidemica Epidemiology—An Experimental and Genomic Approach. Pathogens 2020, 9, 789. [Google Scholar] [CrossRef]
  8. Augot, D.; Muller, D.; Demerson, J.; Boué, F.; Caillot, C.; Cliquet, F. Dynamics of Puumala virus infection in bank voles in Ardennes department (France). Pathol. Biol. 2006, 54, 572–577. [Google Scholar]
  9. Sauvage, F.; Penalba, C.; Vuillaume, P.; Boue, F.; Coudrier, D.; Pontier, D.; Artois, M. Puumala hantavirus infection in humans and in the reservoir host, Ardennes region, France. Emerg. Infect. Dis. 2002, 8, 1509. [Google Scholar] [CrossRef]
  10. Deter, J.; Bryja, J.; Chaval, Y.; Galan, M.; Henttonen, H.; Laakkonen, J.; Voutilainen, L.; Vapalahti, O.; Vaheri, A.; Salvador, A.R.; et al. Association between the DQA MHC class II gene and Puumala virus infection in Myodes glareolus, the bank vole. Infect. Genet. Evol. 2008, 8, 450–458. [Google Scholar] [CrossRef]
  11. Augot, D.; Sauvage, F.; Boue, F.; Bouloy, M.; Artois, M.; Demerson, J.; Combes, B.; Coudrier, D.; Zeller, H.; Cliquet, F.; et al. Spatial and temporal patterning of bank vole demography and the epidemiology of the Puumala hantavirus in northeastern France. Epidemiol. Infect. 2008, 136, 1638–1643. [Google Scholar] [CrossRef] [PubMed]
  12. Dupinay, T.; Pounder, K.C.; Ayral, F.; Laaberki, M.-H.; Marston, D.A.; Lacôte, S.; Rey, C.; Barbet, F.; Voller, K.; Nazaret, N.; et al. Detection and genetic characterization of Seoul virus from commensal brown rats in France. Virol. J. 2014, 11, 32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Castel, G.; Couteaudier, M.; Sauvage, F.; Pons, J.-B.; Murri, S.; Plyusnina, A.; Pontier, D.; Cosson, J.-F.; Plyusnin, A.; Marianneau, P.; et al. Complete genome and phylogeny of Puumala hantavirus isolates circulating in France. Viruses 2015, 7, 5476–5488. [Google Scholar] [CrossRef] [PubMed]
  14. Murri, S.; Madrières, S.; Tatard, C.; Piry, S.; Benoit, L.; Loiseau, A.; Pradel, J.; Artige, E.; Audiot, P.; Leménager, N.; et al. Detection and genetic characterization of Puumala Orthohantavirus S-segment in areas of France non-endemic for nephropathia epidemica. Pathogens 2020, 9, 721. [Google Scholar] [CrossRef]
  15. Monchatre-Leroy, E.; Murri, S.; Castel, G.; Calavas, D.; Boué, F.; Henaux, V.; Marianneau, P. First insights into Puumala orthohantavirus circulation in a rodent population in Alsace, France. Zoonoses Public Health 2018, 65, 540–551. [Google Scholar] [CrossRef]
  16. Dubois, A.; Castel, G.; Murri, S.; Pulido, C.; Pons, J.-B.; Benoit, L.; Loiseau, A.; Lakhdar, L.; Galan, M.; Marianneau, P.; et al. Bank vole immunoheterogeneity may limit Nephropatia Epidemica emergence in a French non-endemic region. Parasitology 2018, 145, 393–407. [Google Scholar] [CrossRef] [Green Version]
  17. Deter, J.; Chaval, Y.; Galan, M.; Gauffre, B.; Morand, S.; Henttonen, H.; Laakkonen, J.; Voutilainen, L.; Charbonnel, N.; Cosson, J.F. Kinship, dispersal and hantavirus transmission in bank and common voles. Arch. Virol. 2008, 153, 435–444. [Google Scholar] [CrossRef]
  18. Reynes, J.-M.; Carli, D.; Renaudin, B.; Fizet, A.; Bour, J.-B.; Brodard, V.; Cart-Tanneur, E.; Dewilde, A.; El Hamri, M.; Fleury, H.; et al. Surveillance des infections humaines par hantavirus en France métropolitaine, 2012–2016. Bull. Epidemiol. Hebd. 2017, 23, 492–499. [Google Scholar]
  19. Riccò, M.; Peruzzi, S.; Ranzieri, S.; Magnavita, N. Occupational hantavirus infections in agricultural and forestry workers: A systematic review and metanalysis. Viruses 2021, 13, 2150. [Google Scholar] [CrossRef]
  20. Christiann, F.; Rayet, P.; Ngueodjibaye, D.B.; Patey, O.; Godefroy, A.; Klein, J.; Lapègue, R.; Gargasson, J.-F.T.-L.; Godfroid, E.; Lafaix, C. Endemic level of Lyme borreliosis in a region of central France: A sero-epidemiologic examination involving blood donors. Eur. J. Epidemiol. 1997, 13, 361. [Google Scholar] [CrossRef]
  21. Rigaud, E.; Jaulhac, B.; Garcia-Bonnet, N.; Hunfeld, K.-P.; Femenia, F.; Huet, D.; Goulvestre, C.; Vaillant, V.; Deffontaines, G.; Abadia-Benoist, G. Seroprevalence of seven pathogens transmitted by the Ixodes ricinus tick in forestry workers in France. Clin. Microbiol. Infect. 2016, 22, 735. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Zhioua, E.; Rodhain, F.; Binet, P.; Perez-Eid, C. Prevalence of antibodies to Borrelia burgdorferi in forestry workers of Ile de France, France. Eur. J. Epidemiol. 1997, 13, 959–962. [Google Scholar] [CrossRef] [PubMed]
  23. Rossi, C.A.; Ksiazek, T.G. Enzyme linked Immunosorbent Assay. In Manual of Hemorrhagic Fever with Renal Syndrome and Hantavirus Pulmonary Syndrome; Asan Institute for Life Sciences: Seoul, Republic of Korea, 1999; pp. 87–91. [Google Scholar]
  24. Niklasson, B.; Lundkvist, A. Indirect immunofluorescent antibody test (IFT). In Manual of Hemorrhagic Fever with Renal Syndrome and Hantavirus Pulmonary Syndrome; Asan Institute for Life Sciences: Seoul, Republic of Korea, 1999; pp. 83–86. [Google Scholar]
  25. Chu, Y.K.; Rossi, C.; Leduc, J.W.; Lee, H.W.; Schmaljohn, C.S.; Dalrymple, J.M. Serological relationships among viruses in the Hantavirus genus, family Bunyaviridae. Virology 1994, 198, 196–204. [Google Scholar] [CrossRef] [PubMed]
  26. El Tinge, J.L.; Yansaneh, I.S. Diagnostics for formation of nonresponse adjustment cells, with an application to income nonresponse in the US Consumer Expenditure Survey. Surv. Methodol. 1997, 23, 33–40. [Google Scholar]
  27. Deville, J.-C.; Särndal, C.-E.; Sautory, O. Generalized raking procedures in survey sampling. J. Stat. Data Sci. Educ. 1993, 88, 1013–1020. [Google Scholar] [CrossRef]
  28. Dohoo, I.R.; Martin, W.; Stryhn, H.E. Veterinary Epidemiologic Research, 2nd ed.; VER Inc.: Charlottetown, PEI, Canada, 2009. [Google Scholar]
  29. Textor, J.; Van der Zander, B.; Gilthorpe, M.S.; Liśkiewicz, M.; Ellison, G.T. Robust causal inference using directed acyclic graphs: The R package ‘dagitty’. Int. J. Epidemiol. 2016, 45, 1887–1894. [Google Scholar] [CrossRef] [Green Version]
  30. The American Association for Public Opinion Research [AAPOR]. Standard Definitions: Final Dispositions of Case Codes and Outcome Rates for Surveys; AAPOR: Deerfield, IL, USA, 2011. [Google Scholar]
  31. Penalba, C.; Goupil, P.; Reska, A.; Reveil, J.; Le Guenno, B.; Camprasse, M. La fièvre hémorragique avec syndrome rénal dans les Ardennes. Étude épidémiologique du personnel de l’ONF. Bull. Epidemiol. Hebd. 1995, 5, 95. [Google Scholar]
  32. Kerouani, G.; Binet, P.; Coudrier, D.; Lallement, P.; Perez, C.; Rollin, P. Seroprevalence of hemorrhagic fever with renal syndrome in foresters of Ile-de-France. Presse Med. 1992, 21, 1636. [Google Scholar]
  33. Lwanga, S.K.; Lemeshow, S. Sample Size Determination in Health Studies: A Practical Manual; World Health Organization: Geneva, Switzerland, 1991. [Google Scholar]
  34. Rollin, P.; Laveran, H.; Gonzalez, J.P.; Coudrier, D.; Sureau, P.; Beytout, D. Fièvre hémorragique avec syndrome rénal dansle centre de la France: Données épidémiologiques. Presse Med. 1987, 16, 4. [Google Scholar]
  35. Bergstedt Oscarsson, K.; Brorstad, A.; Baudin, M.; Lindberg, A.; Forssén, A.; Evander, M.; Ahlm, C. Human Puumala hantavirus infection in northern Sweden; increased seroprevalence and association to risk and health factors. BMC Infect. Dis. 2016, 16, 566. [Google Scholar] [CrossRef] [Green Version]
  36. Penalba, C.; Galempoix, J.; Lanoux, P. Epidémiologie des infections à hantavirus en France. Med. Mal. Infect. 2001, 31, 272–284. [Google Scholar] [CrossRef]
  37. Vapalahti, O.; Lundkvist, Å.; Vaheri, A. Human immune response, host genetics, and severity of disease. In Hantaviruses. Current Topics in Microbiology and Immunology; Springer: Berlin/Heidelberg, Germany, 2001; pp. 153–169. [Google Scholar]
  38. Santé Publique France [SpF], Les Infections à Hantavirus dues au Virus Puumala en France. Available online: https://www.santepubliquefrance.fr/maladies-et-traumatismes/maladies-transmissibles-de-l-animal-a-l-homme/hantavirus/documents/depliant-flyer/les-infections-a-hantavirus-dues-au-virus-puumala-en-france (accessed on 26 October 2022).
  39. Septfons, A.; Figoni, J.; Gautier, A.; Soullier, N.; de Valk, H.; Desenclos, J.-C. Increased awareness and knowledge of Lyme Borreliosis and tick bite prevention among the general population in France: 2016 and 2019 health barometer survey. BMC Public Health 2021, 21, 1808. [Google Scholar] [CrossRef] [PubMed]
  40. Olsson, G.E.; Leirs, H.; Henttonen, H. Hantaviruses and their hosts in Europe: Reservoirs here and there, but not everywhere? Vector-Borne Zoonotic Dis. 2010, 10, 549–561. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Adjusted prevalence of immunoglobulin G against Puumala orthohantavirus antigens among forestry workers in northern France, 2019–2020.
Figure 1. Adjusted prevalence of immunoglobulin G against Puumala orthohantavirus antigens among forestry workers in northern France, 2019–2020.
Viruses 15 00338 g001
Table 1. Demographic features of forestry workers targeted in the Puumala orthohantavirus seroprevalence survey, northern France, 2019–2020.
Table 1. Demographic features of forestry workers targeted in the Puumala orthohantavirus seroprevalence survey, northern France, 2019–2020.
Demographic Features Weighted Distribution
% (n/N)
Age in years, mean (min–max) 44 (18–87)
Age group, years18–4041% (702/1777)
41–5023% (430/1777)
51–6028% (520/1777)
>608% (125/1777)
SexMale94% (1714/1777)
Female6% (63/1777)
Administrative region 1Alsace10% (133/1777)
Lorraine13% (202/1777)
Champagne-Ardenne12% (148/1777)
Bourgogne6% (44/1777)
Franche-Comté9% (145/1777)
Auvergne8% (164/1777)
Hauts-de-France (Picardie, Nord-Pas-de-Calais)9% (143/1777)
Île-de-France3% (33/1777)
Centre9% (197/1777)
Limousin5% (144/1777)
Normandie (Basse Normandie, Haute-Normandie)3% (81/1777)
Bretagne4% (135/1777)
Pays-de-la-Loire7% (208/1777)
Geographic area 2Northeast51% (672/1777)
North central24% (401/1777)
Northwest12% (396/1777)
Auvergne8% (164/1777)
Limousin5% (144/1777)
Occupational activitySilviculturist17% (289/1768)
Woodcutter34% (622/1768)
Forestry machine operator18% (319/1768)
Forest technician/forest ranger21% (357/1768)
Hunting technology/gamekeeper/fishery guardian8% (144/1768)
Gardener/landscaper1% (20/1768)
Other1% (17/1768)
Average weekly exposure time in the forest, forest borders or edges, parks, and wooded gardens exclusively, hours<1010% (160/1758)
10–2012% (199/1758)
>2078% (1399/1758)
Cumulative seniority in forestry activities, years1–1027% (483/1777)
11–2025% (448/1777)
21–3022% (406/1777)
31–4017% (317/1777)
>408% (123/1777)
1 Region of French farmers’ insurance fund affiliation. 2 Northeast includes Alsace, Lorraine, Champagne-Ardenne, Bourgogne, and Franche-Comté regions; North central includes Hauts-de-France, Île-de-France, and Centre regions; Northwest includes Normandie, Bretagne, and Pays-de-la-Loire regions. Note that Auvergne and Limousin regions are located in the center of France.
Table 2. Weighted prevalence and prevalence ratio of immunoglobulin G (IgG)against Puumala orthohantavirus (PUUV) antigens among forestry workers, northern France, 2019–2020.
Table 2. Weighted prevalence and prevalence ratio of immunoglobulin G (IgG)against Puumala orthohantavirus (PUUV) antigens among forestry workers, northern France, 2019–2020.
Demographic Features Number of Workers with Anti-PUUV IgG/Number of Workers TestedWeighted Seroprevalence, (95% CI) 3Weighted Prevalence Ratio 4p-Value for the Prevalence Ratio 5
Age group, years18–408/7021% (0.73–3)Reference-
41–5012/4304% (2–7)3.2 (1.4–7.6)0.008
51–6024/5205% (3–7)3.4 (1.6–7.1)0.001
>606/1256% (3–11)3.9 (1.5–10.2)0.005
SexMale50/17143% (3–5)NA
Female0/63-NA
Administrative region 1Alsace2/1332% (0.44–5)NA
Lorraine8/2024% (2–7)NA
Champagne-Ardenne7/1484% (2–9)NA
Bourgogne2/445% (1–17)NA
Franche-Comté13/1458% (5–14)NA
Auvergne2/1641% (0.33–4)NA
Hauts-de-France (Picardie, Nord-Pas-de-Calais)12/1437% (4–12)NA
Île-de-France1/333% (0.50–16)NA
Centre2/1971% (0.25–3)NA
Limousin1/1441% (0.09–3)NA
Normandie (Basse Normandie, Haute-Normandie)0/81-NA
Bretagne0/135-NA
Pays-de-la-Loire0/208-NA
Geographic area 2Northeast32/6725% (3–6)Reference-
North central15/4014% (2–6)0.77 (0.43–1.4)0.364
Northwest0/396---
Auvergne2/1641% (0.33–4)0.26 (0.07–0.96)0.044
Limousin1/1441% (0.09–3)0.11 (0.02–0.69)0.018
Occupational activitySilviculturist4/2892% (0.63–4)0.49 (0.17–1.4)0.190
Woodcutter15/6223% (2–5)Reference-
Forestry machine operator8/3193% (2–7)1.1 (0.50–2.6)0.763
Forest technician/forest ranger22/3576% (4–9)1.6 (0.87–3.1)0.126
Hunting technology/gamekeeper/fishery guardian0/144---
Gardener/landscaper0/20---
Other1/175% (0.74–23)1.5 (0.23–9.4)0.674
Average weekly exposure time in forest, forest borders or edges, parks, and wooded gardens exclusively, hours<102/1601% (0.32–4)Reference-
10–205/1993% (1–6)2.8 (0.71–10.8)0.141
>2042/13994% (3–5)3.8 (1.2–12.5)0.026
Cumulative seniority in forestry activities, years1–105/4831% (0.50–2)Reference-
11–206/4482% (0.80–4)1.6 (0.53–4.7)0.408
21–3016/4065% (3–8)4.9 (1.9–12.5)0.001
31–4014/3175% (3–8)4.2 (1.7–10.8)0.002
>409/1237% (4–13)5.4 (2.0–14.7)0.001
PUUV, Puumala orthohantavirus, IgG, immunoglobulin G. 1 Region of French farmers’ insurance fund affiliation. 2 Northeast includes Alsace, Lorraine, Champagne-Ardenne, Bourgogne, and Franche-Comté regions; North central includes Hauts-de-France, Île-de-France, and Centre regions; Northwest includes Normandie, Bretagne, and Pays-de-la-Loire regions. Note that Auvergne and Limousin regions are located in the center of France. 3 Weighted seroprevalence. Seroprevalence is not adjusted for the sensitivity and specificity of the diagnostic test. 4 Calculated using weighted univariable Poisson regression models with hantavirus seropositivity (yes/no) as the outcome and the variable on the left as the exposure. The association between hantavirus and average weekly exposure time in forests, forest borders or edges, parks, and wooded gardens was adjusted for occupational activity. 5 The overall p-values were 0.007 for age, 0.025 for the geographic area, 0.056 for occupational activity, 0.218 for average weekly exposure time in the forest, and 0.001 for cumulative seniority in the forestry sector.
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

Krug, C.; Rigaud, E.; Siby-Diakite, D.; Bénézet, L.; Papadopoulos, P.; de Valk, H.; Deffontaines, G.; Septfons, A.; Reynes, J.-M. Seroprevalence of Hantavirus in Forestry Workers, Northern France, 2019–2020. Viruses 2023, 15, 338. https://doi.org/10.3390/v15020338

AMA Style

Krug C, Rigaud E, Siby-Diakite D, Bénézet L, Papadopoulos P, de Valk H, Deffontaines G, Septfons A, Reynes J-M. Seroprevalence of Hantavirus in Forestry Workers, Northern France, 2019–2020. Viruses. 2023; 15(2):338. https://doi.org/10.3390/v15020338

Chicago/Turabian Style

Krug, Catarina, Emma Rigaud, Dieyenaba Siby-Diakite, Laetitia Bénézet, Pavlos Papadopoulos, Henriette de Valk, Gaëtan Deffontaines, Alexandra Septfons, and Jean-Marc Reynes. 2023. "Seroprevalence of Hantavirus in Forestry Workers, Northern France, 2019–2020" Viruses 15, no. 2: 338. https://doi.org/10.3390/v15020338

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