Contrasted Epidemiological Patterns of West Nile Virus Lineages 1 and 2 Infections in France from 2015 to 2019

Since 2015, annual West Nile virus (WNV) outbreaks of varying intensities have been reported in France. Recent intensification of enzootic WNV circulation was observed in the South of France with most horse cases detected in 2015 (n = 49), 2018 (n = 13), and 2019 (n = 13). A WNV lineage 1 strain was isolated from a horse suffering from West Nile neuro-invasive disease (WNND) during the 2015 episode in the Camargue area. A breaking point in WNV epidemiology was achieved in 2018, when WNV lineage 2 emerged in Southeastern areas. This virus most probably originated from WNV spread from Northern Italy and caused WNND in humans and the death of diurnal raptors. WNV lineage 2 emergence was associated with the most important human WNV epidemics identified so far in France (n = 26, including seven WNND cases and two infections in blood and organ donors). Two other major findings were the detection of WNV in areas with no or limited history of WNV circulation (Alpes-Maritimes in 2018, Corsica in 2018–2019, and Var in 2019) and distinct spatial distribution of human and horse WNV cases. These new data reinforce the necessity to enhance French WNV surveillance to better anticipate future WNV epidemics and epizootics and to improve the safety of blood and organ donations.


Comparison of WNV Seasonal Patterns in
In France during the last five years, the most important equine WNV outbreaks were reported in 2015. Of the 49 cases reported in horses in 2015, 41 presented neuroinvasive forms, three febrile forms, and five asymptomatics were identified thanks to serosurveys implemented in identified WNV transmission foci [11]. There was no equine case noticed in 2016, while one asymptomatic horse was detected in 2017 in the vicinity of WNV human cases. A total of 13 and nine neuroinvasive cases were reported in 2018 and 2019 respectively, while four additional febrile forms were reported in 2019 ( Figure 1). Natural death or euthanasia occurred in 14.6% (6/41) and 15.4% (2/13) of the horses with was detected in 2017 in the vicinity of WNV human cases. A total of 13 and nine neuroinvasive cases were reported in 2018 and 2019 respectively, while four additional febrile forms were reported in 2019 ( Figure 1). Natural death or euthanasia occurred in 14.6% (6/41) and 15.4% (2/13) of the horses with West Nile neuroinvasive disease (WNND) in 2015 and 2018 respectively whereas no equine death was reported in 2019. These percentages are lower than usually described in the literature [30]. For the same period, for humans, the highest number of autochthonous cases was reported in 2018 with 26 laboratory-diagnosed human cases, including seven WNND, 18 febrile and one asymptomatic form. Interestingly, one blood donor, symptomatic a few days after the donation and one organ asymptomatic donor were tested positive for WNV the same year [31]. The total number of cases in 2018 represented a 27-fold increase compared with the 2015-2017 transmission seasons during which one, none, and two febrile cases were reported in 2015, 2016, and 2017 respectively. In 2019, one febrile and one WNND were reported ( Figure 1). WNND have generally occurred more frequently in horses than in humans in France so far.
For the first time in France since the implementation of an integrated WNV surveillance system, WNV infections were also reported through WNV surveillance in the avifauna in 2018. In total, four raptors (two northern goshawks (Accipiter gentilis), one common buzzard (Buteo buteo), and one long-eared owl (Asio otus) were diagnosed WNV positive in September and October 2018 in Corsica (owl) and Alpes-Maritimes (diurnal raptors) by the SAGIR network (a French network dedicated to wildlife disease surveillance). All these wild birds were found alive and suffered from serious nervous disorders.

Shifts in Temporal and Spatial Distribution of WNV Cases
In 2015, the onset of the first equine case was reported on week 33 (starting the 11 August) with a peak on week 38 (14 September to 20 September) and the last case was notified on week 44 (Figure 2a). For the same period, for humans, the highest number of autochthonous cases was reported in 2018 with 26 laboratory-diagnosed human cases, including seven WNND, 18 febrile and one asymptomatic form. Interestingly, one blood donor, symptomatic a few days after the donation and one organ asymptomatic donor were tested positive for WNV the same year [31]. The total number of cases in 2018 represented a 27-fold increase compared with the 2015-2017 transmission seasons during which one, none, and two febrile cases were reported in 2015, 2016, and 2017 respectively. In 2019, one febrile and one WNND were reported ( Figure 1). WNND have generally occurred more frequently in horses than in humans in France so far.
For the first time in France since the implementation of an integrated WNV surveillance system, WNV infections were also reported through WNV surveillance in the avifauna in 2018. In total, four raptors (two northern goshawks (Accipiter gentilis), one common buzzard (Buteo buteo), and one long-eared owl (Asio otus) were diagnosed WNV positive in September and October 2018 in Corsica (owl) and Alpes-Maritimes (diurnal raptors) by the SAGIR network (a French network dedicated to wildlife disease surveillance). All these wild birds were found alive and suffered from serious nervous disorders.

Shifts in Temporal and Spatial Distribution of WNV Cases
In 2015, the onset of the first equine case was reported on week 33 (starting the 11 August) with a peak on week 38 (14 September to 20 September) and the last case was notified on week 44 (Figure 2a). Each of the 49 equine cases were identified in three departments surrounding the Camargue area, a large region delineated by the Rhone delta and characterized by high biological and environmental diversity in South-Eastern France. A total of 33 confirmed cases corresponding to 26 distinct outbreaks were located in Bouches-du-Rhône department, 15 confirmed cases (12 outbreaks) in Gard department, and one in Hérault department. Only one WNF human case was confirmed later in the season (2 October 2015, onset of symptoms 27 September) in Gard department and one mosquito pool corresponding to Culex pipiens mosquitoes was found positive in the same area on 11 September 2015 (Figure 3c,d). Each of the 49 equine cases were identified in three departments surrounding the Camargue area, a large region delineated by the Rhone delta and characterized by high biological and environmental diversity in South-Eastern France. A total of 33 confirmed cases corresponding to 26 distinct outbreaks were located in Bouches-du-Rhône department, 15 confirmed cases (12 outbreaks) in Gard department, and one in Hérault department. Only one WNF human case was confirmed later in the season (2 October 2015, onset of symptoms 27 September) in Gard department and one mosquito pool corresponding to Culex pipiens mosquitoes was found positive in the same area on 11 September 2015 (Figure 3c,d).
In 2017, two human cases were diagnosed on the 21 August and 4 September in Alpes-Maritimes department. It was the first time in France that WNV was reported in this area. Following these cases, a serosurvey was carried out on 151 equids from a horse center located in the vicinity of the second WNV human case. Only one asymptomatic horse (1/151; CI 0%−2.9%) was found WNV-IgM positive (Figure 3e,f).
The WNV transmission season started much earlier and finished later in 2018 compared to the 2015-2017 period. Even if horses are particularly sensitive to the infection and can be used as indicators of virus circulation [32], autochthonous human cases were diagnosed before the occurrence of horse cases from 19 July 2018, while symptoms onset dated back to early July (week 27, 2-8 July). The number of human cases peaked on week 33 and the last case was notified week 46 (12)(13)(14)(15)(16)(17)(18) (Figure 2b). Four French departments around the Mediterranean Sea reported WNV infections (i.e., Alpes-Maritimes, Bouches-du-Rhône, Vaucluse, and Eastern Pyrenees). Alpes-Maritimes department was the first area reporting WNV cases and a cluster of 21 cases were located exclusively in this area (Figure 3h). Three raptors were found WNV positive later on in the season in the same area (in Nice and Antibes, Figure 3g). For the first time in 2018, an outbreak occurred in the south of Corsica with the first human case diagnosed at week 32 (6-12 August) and a second one at week 39 (24)(25)(26)(27)(28)(29)(30). Concomitantly, the onset of WNV equine outbreaks was reported on week 35 (starting 27 August) in North Corsica (Bastia) and four horses were found positive in 2018 in Corsica ( Figure 3g). Finally, one long-eared owl (Asio Otus) was also diagnosed positive in South Corsica. Interestingly, the other 2018 equine cases were located in two departments in the Camargue area, namely Bouches-du-Rhône and Hérault, already affected by WNV equine outbreaks in 2015 ( Figure 3g). The notification of WNV cases in equids started later than in humans and at a comparable period than in 2015 (week 33, 11 August 2015 versus week 35, 29 August 2018).
WNV activity in 2019 was lower than in 2018. It was characterized by a circulation of the virus in the Camargue region with WNV infected horses (n = 11) and by the resurgence of WNV disease in Corsica (n = 2 horses). Moreover, one WNF and one WNND human cases were detected in Var, a department with sporadic WNV transmission to humans and horses identified since 2003 [33]. Pathogens 2020, 9, x FOR PEER REVIEW 5 of 16 In 2017, two human cases were diagnosed on the 21 August and 4 September in Alpes-Maritimes department. It was the first time in France that WNV was reported in this area. Following these cases, Recent changes in the temporal and spatial distribution of French WNV cases can be highlighted from 2015-2019 data analysis. Specifically, a multiplication of circulation foci have been reported Pathogens 2020, 9, 908 6 of 16 during the last three years, with the emergence of WNV and recent description of clinical cases in the departments of Alpes-Maritimes, Var, and French Corsica island alongside the usual enzootic WNV circulation in the Camargue area during most of the period (2015, 2018-2019). WNV emergence in Alpes-Maritimes was associated with an increase of reported WNV cases in humans and birds but not in equids, while the distinct spatial distribution of human (mostly in Alpes-Maritimes) and horse (mostly in Camargue) WNV cases have been observed in Southern France these last years.
The intensity of WNV circulation and transmission is shaped in part by mosquito vector abundance and is influenced by biotic and abiotic factors favorable for mosquito proliferation [34,35]. Mosquito abundance was found to be significantly lower in Hérault than in the other two departments for the period 2015 to 2019 but the number of traps was six for this department compared to eight for the Bouches-du-Rhône and Gard (see Section 4.2 material and methods). The analysis of Culex pipiens mosquito abundance in the Camargue area (Bouches-Du-Rhône, Gard, and Hérault) during the vector season (June to October) indicated that the total number of trapped mosquitoes was significantly higher in 2018 than in 2016, 2017, and 2019 (p values ≤ 0.002). The difference observed between 2015 and 2018 was not found significant ( Figure 4b and Table 1). Moreover, since the importance of WNV outbreaks in Europe was found to be strongly correlated with the length of the mosquito proliferation season, early abundance of mosquitoes in June was compared in 2015-2019. The vector season started earlier in 2018 than in other years, as the number of mosquitoes trapped in June was significantly higher in 2018 than in 2015, 2016, 2017, and 2019 (p values < 0.0007) ( Table 1). The abundance ratio for June 2015 was 0.20, thus corresponding to an abundance five times higher in June 2018 than in June 2015, after controlling for the effect of the department (Figure 4a, Table 1).
Pathogens 2020, 9, x FOR PEER REVIEW 7 of 16      (Table S1). 9-14, and eight amino acid substitutions affecting different viral genes were observed between WNV-Akela/France/2015 and older French lineage 1 isolates and between lineage 2 WNV-7025/France/2018, and the closely genetically related WNV-Cremona4/Italy/2014 respectively (Tables S2 and S3). No amino acid substitutions correspond to established WNV molecular virulence determinants or to positively selected codons [37,38].   [36]. The optimal tree generated using the Neighbor-Joining method with the sum of branch length = 1,28228091 is shown. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) are shown next to the branches. The tree is drawn to scale, with branch lengths expressed in the same units as for the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Jukes-Cantor method and are in the units of the number of base substitutions per site. The analysis involved 50 nucleotide sequences, including Japanese Encephalitis Virus as an outgroup. All positions containing gaps and missing data were eliminated. There were a total of 10,475 positions in the final dataset.

Discussion
In France, WNV caused outbreaks involving several human and horse cases in the beginning of the 1960s before it disappeared for 35 years [7]. During the 2000-2008 period, four episodes of WNV transmission were reported in France. Only WNV lineage 1 was reported in France during this period, with closely genetically related WNV isolates belonging to the Western Mediterranean clade identified in the Camargue area in 2000 (in horses) [39] and in 2004 (in birds) ( Figure 5) [40]. It re-emerged following a cyclical and hardly predictable pattern and was mostly limited to Camargue, a high-risk area for WNV circulation due to high concentration of wetlands, mosquitoes, wild birds, and horses [41][42][43]. The same spatial Camargue location of equine outbreaks was pointed out in 2015, 2018, and 2019. The enhanced abundance of WNV competent mosquitoes in 2015 and 2018 and an earlier vector season in the Camargue area in 2018 have been identified as potential risk factors for higher WNV transmission. Moreover, the phylogenetic analysis of one WNV strain identified recently in Camargue supports WNV enzootic transmission in this region as it revealed that the sequence of WNV isolated from a confirmed equine case in 2015 is close to the lineage 1 strain that circulated in France in 2004 [44,45]. The percentage of nucleotide identity between French WNV lineage 1 isolates (>98.3%, Table S1) is coherent with the mean evolutionary rate of the European WNV strains (3.7 × 10 −4 substitutions/site/year) [46]. Most of the mutations distinguishing the viral isolates were synonymous and homogenously distributed along the viral genome, which suggests that the genetic evolution of French WNV strains arose through a strong and local diversifying selection.
In 2018, WNV lineage 2 belonging to the Central and Eastern European clade (CEC) as defined by Ziegler et al. [47] was isolated for the first time from wild raptors, which have been shown to be particularly susceptible to WNV neuro-invasive infections, in Alpes-Maritimes in France [18,48,49]. WNV lineage 2 emergence in France was associated with exceptional WNV activity and lineage 2 spread in Western and Northern-most territories (Germany) in Europe this same year [25]. During this year, WNV infections in Europe increased dramatically compared to previous transmission seasons. From June to November 2018, a large part of Europe faced a period of unusually hot weather that led to record-breaking temperatures [50]. Like in France, European WNV infections started earlier in 2018 than in previous years. Indeed, the first WNF cases were reported on 31 May (week 22) in Greece which is the earliest disease onset compared with previous years [51]. At the end of 2018, a total of 1503 human infected confirmed cases were reported in 11 countries of the European Union (with almost 92% of cases coming from Italy, Greece, Romania, and Hungary) [52]. This number exceeded the cumulative number of WNV reported infections of the seven previous years [53]. The highest increase compared to previous transmission season was observed in Bulgaria (15 fold) followed by France (13.5 fold), and Italy (10.9 fold) [53]. During the 2018 transmission season, reports from the ECDC also indicated a high transmission among horses with 285 outbreaks reported by European member states as follows: 149 in Italy, 91 in Hungary, 15 in Greece, and 13 in France representing an increase of 30% in comparison with the number of outbreaks in 2017 [53].
In particular, in 2018, there was a large WNV lineage 2 outbreak in Northern Italy, including the Piemonte regions. WNV lineage 2 circulation was first documented in Italy in 2011 and, since then, has settled in Northern Italy at least since 2013 [14,54]. Recent phylogenetic analysis [55] revealed that two Italian lineage 2 strains, namely clade A and clade B diverged between 2010 and 2012 from a central region of the Po Valley. Clade A spread towards Northeastern Italy and apparently became extinct in 2013-2014, whereas clade B spread north-west reaching the most Western regions of Italy. Such a WNV short distance introduction via infected birds coming from a neighboring country has been usually hypothesized [18,47]. According to the high percentage of nucleotide homology (99.76%) between the Italian 2014 clade B and French 2018 WNV strains (Table S1), we hypothesize that WNV lineage 2 gradually spread from Northern Italy to South-Eastern France in 2018 or 2017, considering that WNV lineage 2 was isolated from wild birds in Alpes-Maritimes in 2018 and that human cases were already reported in the same area in 2017. Moreover, Histidine instead of Proline residues are found at position 249 in the helicase part of NS3 in recent Italian and French lineage 2 strains. The role of this genetic modification in the modulation of WNV pathogenicity in mammalian and bird hosts has been regularly debated [56,57]. While a specific non-synonymous substitution, Lys2114Arg, was identified in several WNV lineage 2 isolates obtained in Germany in 2018 and could be associated with an increased WNV fitness as such a mutation was not evidenced in the French 2018 lineage 2 strains.
In 2018, most human cases in mainland France were reported in areas with WNV-related bird mortality, which is consistent with the positive relationship between WND human cases and seroprevalence level in passerine birds, previously demonstrated at the European level [58]. The absence of horse cases before the onset of human cases could be explained by the very low density of horses in the Alpes-Maritimes area and by the location of 13 out of 26 human cases in the urban city of Nice [59]. The increase observed in 2018 in WNV human cases in France could result from differing virulence or transmission properties for humans and for horses of WNV lineages 1 and 2 and from varying animal and human densities in areas reporting WNV infections in 2018 (Alpes-Maritimes with densely populated urbanized areas, while most horse cases were reported in a natural wetland, the Camargue area); the first hypotheses (with more transmission to and more cases in humans associated with WNV lineage 2 infections) would deserve more attention but are currently not supported by the literature.
We also document in 2018 the detection of WNV for the first time in French Corsica Island. This finding is not surprising as a serosurvey carried out in Corsica in 2014 highlighted that 9.4% of horses presented WNV antibodies. Among these positive horses, 66.6% were native from the island, indicative of WNV local circulation [60]. The identification of WNV clinical cases in humans, horses, and birds in Corsica further documents recent and active circulation of the virus. Nevertheless, the identification of the causative lineage could not be achieved as a low viral load was evidenced in clinical specimens collected on the island (one raptor, a long eared-owl) and as both lineages 1 and 2 were described recently in Sardinia and Italy [61].
Fewer outbreaks were reported in 2019 than in 2018 in France and Europe, but a changing epidemiological pattern of WNV circulation can be anticipated in France in the coming years. Indeed, the introduction of WNV lineage 2 in Hungary in 2004 was followed by strain adaptation and limited activity in 2005-2007 while extensive spread of the virus was reported from 2008 [18]. Moreover, a remarkable extension of the distribution area of WNV lineage 2 has been evidenced in 2018. Specifically equine cases and mortality on resident wild and captive birds were detected for the first time in Eastern and Southeastern Germany [25]. This introduction was followed in Germany by an increase of equine WNV outbreaks in 2019 and the reporting of the first five confirmed mosquito-borne autochthonous human cases [47]. Another important finding during the 2018 transmission season relates to WNV genome detection in one blood donor for the first time in France by the French blood establishment [31]. Interestingly this donor had spent time in Alpes-Maritimes before the occurrence of the first animal or human WNV cases. These new data on WNV spread in South France and on positive WNV screenings in blood products and organ transplants highlight the necessity to strengthen WNV integrated surveillance in France in order to primarily secure human blood, cells, or organ products. WNV surveillance has been mainly supported by clinical surveillance programs focusing on the analysis of moribund or dead birds and of horses and human patients with neuroinvasive signs, which may lack sensitivity and fail to detect low-level circulation. A combination of clinical event-based surveillance activities and active monitoring of WNV enzootic transmission through the regular monitoring of seroconversions in sentinel and/or resident birds and horses or through mosquito trapping and WNV screening would enhance the chance to early detect WNV transmission [62].
Finally, an enhanced transmission of WNV in Southeastern France in 2018 paralleled an unusually high number of outbreaks of another Culex-borne flavivirus, Usutu, and in most French metropolitan territories [63]. Such findings emphasize the need of unraveling the virological, ecological, and climatic factors responsible for Culex-borne flavivirus emergence in France and Europe [64,65].

Samples
Main organizational aspects of the French West Nile virus surveillance system in animals, humans, and vectors have been described previously in the article of Bahuon et al. [11]. Briefly, the surveillance is based on clinical case definition (human and equine) or criteria for dead or sick birds' reports, collection, and testing. No routine indicator-based surveillance is implemented on the animal population. Diagnostic specimens are from 1/suspect human cases, as well as Culex mosquito populations sampled in affected areas once the viral circulation has been confirmed. They are analyzed by the National Reference center (NRC) for arboviruses (IRBA-Armed Forces Biomedical Research Institute) 2/each horse and avian suspect cases confirmed by the National Reference laboratory (NRL) on West Nile virus (Anses, Animal Health Laboratory, Maisons-Alfort) [12]. Suspect West Nile cases correspond to human patients over 15 years old and equids presenting with fever (≥38.5 • C) and symptoms of viral meningitis or encephalitis; wild or captive birds (raptors, corvids, and turdids more specifically) found dead, and individuals displaying neurological symptoms during the surveillance period (1 June to end of November) in the at-risk area (i.e., counties in the Mediterranean area). Moreover wild bird surveillance has been extended to departments considered, according to a statistical model, at an increased risk of WNV transmission, and located along the Mediterranean Sea and the Rhone River in South Eastern France, as well as in Bas-Rhin in North Eastern France since 2019 [58].

Mosquito Collection
Mosquitoes were collected weekly from mid-May to late October, corresponding to the mosquito season in the Rhône Delta, Camargue. CDC-like traps (John W. Hock Company, Gainesville, FL, USA) were used without light and were baited with carbon-dioxide dry ice (−80 • C). The trapping network was composed of 8 traps in the departments of Gard, 8 traps in Bouches-du-Rhône, and 6 traps in Hérault departments. Mosquitoes were stored in the fridge, killed, and identified with identification morphological keys.

Serology
Blood samples were collected in dry tubes, allowed to clot, and centrifuged at 1500 rpm for 10 min and stored at +4 • C during 1 month at most or at −20 • C for long term archiving.
For equine suspected cases reported to the French NRL, sera were first screened for anti-WNV antibodies by competition ELISA (ID Screen West Nile competition kit, IDVet Company, Montpellier, France) in local veterinary laboratories. Then IgG positive sera were further analyzed by M-antibody capture ELISA for IgM detection (ID screen West Nile IgM capture, IDVet company, Montpellier, France) in local veterinary laboratories and confirmed at the NRL. Analysis and interpretation of ELISAs were performed according to the manufacturer's instructions. In the event of IgM positive screening, the first samples collected during WNV outbreaks were confirmed by microneutralization test (MNT) as described in Beck et al. [66]. A confirmed case was therefore defined as a clinical suspected horse with at least a positive IgM ELISA test.
For human WNV diagnosis, sera and cerebrospinal fluid (CSF) were tested by in-house ELISAs (indirect IgG and MAC-ELISAs) using precipitated and inactivated virus. A case of WNV infection is confirmed with the presence of IgM in CSF and/or IgM and IgG in sera and anti-WNV neutralizing antibodies [67].

Real-Time RT-PCR
Brain of horses and birds, EDTA blood, and CSF suspected to be infected with WNV were stored at −80 • C until analysis. Brains were grinded in Dulbecco modified Eagle's minimal essential medium (DMEM) with ceramic beads (MP Biomedicals, Illkirch, France) and FastPrep ribolyzer in BSL3 facilities. A total of 560 µL of Lysis buffer from the QIAamp Viral RNA kit (Qiagen, Hilden, Germany) were added to 140 µL of grinded material before RNA extraction with the automate QIAcube. Human samples (EDTA blood and CSF) were processed the same way. Every RNA extracts were subjected to real time (rt) RT-PCR following the protocol described earlier by Linke et al. [68].

Virus Isolation
One milliliter of brain homogenates of WNV rtRT-PCR positive wild birds and horses was prepared in DMEM culture medium and inoculated on T25 flask that had been seeded with Vero NK cells (ATCC: CCL81™), 24 h earlier and washed with DMEM before inoculation. After 1h 30 of incubation at 37 • C with 5% CO 2 , cells were washed twice with phosphate buffered saline (PBS), and complete medium (DMEM+ 1% penicillin-streptomycin+ 1% sodium pyruvate + 5% fetal calf serum) was added. The cells were observed each day from 3 days to 7 days post infection (pi). As soon as cytopathic effects (CPE) were detected, the supernatant was collected, stored at −80 • C, and RNA extracts subjected to rtRT-PCR to confirm WNV detection. Primary isolation was followed by a passage on Aedes albopictus (C6/36) (ATCC ® CRL1660™) cell line. A total of 200µL−1 mL of Vero cell supernatants was added to T25 flask that had been seeded with C6/36 24 h earlier and washed with Leibowitz L15 medium before inoculation. After 1h 30 of incubation at 28 • C without CO 2 , 6 mL of Leibowitz L15 media + 1% penicillin-streptomycin+ 1% sodium pyruvate + 1% L Glutamin+ 10% fetal calf serum were added. CPEs were not systematically observed in C6/36 cells and supernatants were collected on day 7 post-infection at the latest and tested as described above. This protocol is adapted from the OIE Manual of Diagnostic Tests and vaccines for Terrestrial Animals [69].

Nucleotide Sequencing and Sequence Analysis
Sequencing libraries were prepared from genomic RNAs extracted from virus isolate (2015) or from organ homogenates (2018) and whole-genome sequencing data were obtained as previously described (Ion Torrent sequencing and assembly with CLC Genomics Workbench for Genbank accession number MT863559, 2015 [70] or with bwa for Genbank accession numbers MT863560-1, 2018 [71]). Multiple alignment of the nucleotide sequences was performed using the ClustalW algorithm and phylogenetic analysis was performed using the Neighbor-Joining and Maximum Likelihood methods in MEGA7 [36].

Statistical Analysis
We used negative binomial generalized linear models to analyze the mosquito trapping data. The dependent variable was the number of trapped Cx. pipiens, and the independent variables were the department (Bouches-du-Rhône-reference class, Gard, or Hérault) and the year (2015-2019, the reference class being 2018). Two models were separately fitted: One for the yearly total number of trapped Cx. pipiens, and the other for the number of Cx. pipiens trapped in June. Statistical analyses were performed using R 3.6.1 [72].