West Nile virus in Portugal

Literature reviews

We queried PUBMED for articles related to WNV and Portugal using the query “((west nile virus) AND (Portugal)) OR (vírus nilo AND Portugal)” in all fields; thus leaving the query non-specific enough to capture a wide range of studies, and searching for articles both in Portuguese and English. The query returned 42 hits as of 27th January 2020, from which 17 did not relate directly to the circulation of WNV or its vectors in the country. A second literature search was performed manually among Portuguese academic journals and PUBMED. References within the studies found in our two searches were also scanned for potential studies related to WNV circulation in Portugal. We also scanned the GenBank nucleotide database for existing WNV sequences from Portugal, searched with the query “((west nile virus) AND (Portugal)) OR ((WNV) AND (Portugal))” OR (vírus nilo AND Portugal). The query returned 70 hits as of 27th January 2021, of which 62 did not relate to WNV genomic regions.

We provide the results of our literature review in the form of tables in the main text. Studies found in the PUBMED and manual searches are grouped by content type in Table 1. All evidence related to the presence of WNV in Portugal across all studies found is summarized in Table 2 (also available in a permanent figshare repository [30]). Studies presenting time series of relevant mosquito species in Portugal are summarized in Table 3. All publicly available WNV genetic sequences obtained in Portugal are listed in Table 4.

Mosquito data

We used mosquito data collected by the REVIVE initiative in Portugal [27]. Data included mosquito counts (adult and immature stages) for Culex pipiens per county in Portugal during the period 2017-2019. The REVIVE data is in the public domain [27].

Transmission suitability

The basic reproduction number R₀ of WNV in the zoonotic reservoir can be summarised as the product of two factors: the number of adult female mosquitoes per host M, and the transmission potential of each adult female mosquito P (i.e. R₀=MP). Both M and P oscillate in time, driven by the influence of multiple eco-climatological factors. M, the ratio between the population size of mosquitoes V and hosts N (M =V / N), depends on the local effects of climate on mosquito survival and reproduction, affecting V [31], as well as on the reproduction and migratory behaviours of the avian species, affecting N [32]. For most epidemiological contexts and pathogens, high resolution data on either V or N is often lacking. The transmission potential P depends on viral, mosquito and host factors (e.g. incubation periods, lifespans, etc), some of which depend on local climate variations [31,33,34]. Since adequate geo-temporal data for V and N is generally difficult to obtain, many recent theoretical epidemiology studies resort to estimating the transmission potential of mosquito-borne viruses using ‘suitability indices’. These indices are functions of the parameters that compose either M or P or both, and generally depend on covariates such as climatic variables, vegetation indices, host demography, mosquito and disease incidence reports, etc [34–36]. While a suitability index is an incomplete measure of the real transmission potential of a mosquito-borne virus, much work has been done to demonstrate the potential of an index to capture the geo-temporal dynamics of the number of reported mosquito-borne infections in host species [34,37,38].

In this study, we assessed local WNV transmission potential by estimating one of such suitability indices - termed the index P (the transmission potential of each adult female mosquito from R₀=MP). The theory and practice of P has been previously described in full by Obolski et al., that also introduced the MVSE R-package capable of estimating P given some epidemiological priors and climatic data [34]. The index P has been used in several studies focusing on the population dynamics of the dengue, chikungunya, and Zika viruses (e.g. [34,37,39,40]). Recently, the index was also successfully used to estimate the transmission potential of WNV in Israel dependent on local relative humidity and temperature variables [41]. In that study, P was interpreted as a proxy for the risk of spillover to the human population given that it measured the WNV transmission potential by single adult female mosquitoes in the animal reservoir. In the current study, we use the same approach as applied to Israel, informed by the epidemiological priors applied in that study (see Table in main text of [41]), which relate to spp. Culex, WNV and a theoretical, average bird species.

To inform the estimation of the index P, climatic data was downloaded from the Copernicus platform at www.copernicus.eu. We used the available dataset “essential climate variables for assessment of climate variability” [42], which includes climatic variables at a time resolution of 1 month (1981-2019), and gridded spatial resolution of 0.25°x0.25°. In the current study we used MVSE R-package version 1.01r [43], which uses temperature, relative humidity and precipitation as input data. For each latitude-longitude Portuguese data point in the Copernicus dataset, we estimated the index P as described by Obolski et al. for Brazil, the Honduras and South America [34]. The resulting geotemporal estimations, as presented in the main text, are made available as a comma separated file in a permanent figshare repository [30].

Analysis of peak suitability timing and size

Peaks in the monthly index P time series 1981-2019 per county were found by the maximal value and month where it occurred, using a circular mean to calculate the average time point of the year for the north and south regions, independently. The distribution of values throughout the years were then compared between regions with a circular ANOVA yielding a p=0.005. Circular data was handled with the R-package ‘circular’ [44]. Mean peak values were calculated across counties for each year, using Welch’s t-test to compare north and south regions yielding a p=1 ×10⁻¹⁶.

Analysis of long-term climatic and suitability trends

We applied linear quantile mixed models to the estimated index P, temperature, precipitation and humidity between 1981 and 2019, for which the Copernicus climate dataset was available. Our data contained yearly median values of the mentioned variables for each Portuguese county. Since the data included highly irregular points, which might skew standard linear regression models, we chose to use the more robust quantile regression using the median (τ =0.5). Furthermore, as different regions have repeated measures, we employed a mixed-effects quantile regression (using the ‘lqmm’ R-package [45]). Our models contained a random intercept term for each region, a fixed effect dummy variable coding for north versus south, and a continuous variable representing the time unit (year). We tested the interaction between the year and north-south terms and found it statistically significant (p<0.01) in all models except when index P was the dependent variable (p>0.7) for which the term was dropped. All the other fixed effects were statistically significant for all models (p<0.001) and are summarized in Table 5.

Analysis of long-term suitability periodicity

Detrending and wavelet analysis followed the workflow of Damineli et al. [46]. Detrending with a Maximal Overlap Discrete Wavelet Transform was performed in a multiresolution analysis with the R-package ‘wavelets’ [47]. Significant periods of the continuous wavelets were assessed considering the spectrum of an autoregressive model of order 1, as a null model using the R-package ‘biwavelet’ [48].

Equine WNV surveillance data 2016-2020

Between 2016 and 2020 the virology laboratory of INIAV performed molecular and serology testing in 165 horses from several regions of the country, reacting to local suspected neuroinvasive disease. RNA was extracted by using the BioSprint 96 workstation with the MagAttract 96 cador pathogen kit, according to manufacturer’s instructions (Qiagen, Hilden, Germany). The samples were screened for WNV by real-time reverse transcriptase PCR (RT-qPCR) targeting the NS2A gene [49]. Equine sera were tested for WNV-IgM antibodies by capture ELISA (ID Screen West Nile IgM Capture ELISA, IDVET, Montepellier, France). Twenty three horses were diagnosed as WNV-positive on the basis of IgM positivity, with one horse being also positive by RT-qPCR (July 2020, Escalos de Cima-Castelo Branco, center Portugal). Four of eight horses that were followed for clinical progression died or were euthanized. These results are summarized in Table 2, including the time and region of sampling.

Evidence for West Nile virus circulation in Portugal

The earliest evidence of WNV circulation in Portugal came from serological surveys in the south of the country in 1966-1967, reporting high titres of hemagglutination-inhibiting antibodies and neutralizing capacity against WNV in animal serum [50,51]. Such findings led to several initiatives in 1969 to create institutional links between the national health and veterinary medicine services with the purpose of seeding future studies and raising awareness to the virus [52]. Some of the earlier initiatives included the capture and testing of animals (e.g. equines, rodents, birds, mosquitoes) and serological testing of human serum [51–54], eventually leading to the first viral isolation from the mosquito species Anopheles maculipennis [53].

After a period of three decades with seemingly no research output in the country, further evidence for the circulation of WNV in Portugal accrued from several studies. Between 1999 and 2002, Formosinho and colleagues reported serological data from several avian species sampled across natural reserves and from equines sampled in the district of Santarém [55]. A few studies in the period 2004-2010 also reported positive findings in several avian species [56], mosquitoes [57,58], and equines [56]. More recently, between 2015 and 2020, following reports of several clinically compatible equine outbreaks in the south of the country, a number of equines were found serologically positive to the virus in the Alentejo and Algarve [59,60]. Some of this evidence is reported for the first time in this study (see Methods for details, and Table 2 for summary).

Although such evidence supports a long history of local circulation in animals and mosquitoes, autochthonous WNV infections in humans in Portugal were only first reported in 2004. This was the case of two Irish bird-watchers found positive for both WNV IgG and IgM antibodies. They visited Ria Formosa, a marshland rich in endemic and migratory birds in the Algarve between the 26th June and 10th of July [61]. Two studies published thereafter were able to detect WNV by RT-PCR in pools of mosquitoes collected near the putative region of infection that year [57,58]. A third case was reported in July 2010 in the region of Setúbal, from a local rural resident with no history of international travel [62]. The resident tested positive for both WNV IgG and IgM antibodies in the weeks following the clinical diagnosis. Reactive mosquito collection performed around the patient’s residence between the 13-18th of July revealed species suited for WNV transmission (e.g. Culex pipiens), although all mosquito pools tested negative by RT-PCR [62]. Two months later, two equine cases were notified about 4 kilometers away from the reported human case. The most recent human infection was reported in a rural area of Almancil (Algarve) in 2015 in a patient testing positive for both WNV IgG and IgM antibodies and viral neutralization, with no record of international travel or history of vaccination against flaviviruses [63]. In August and September of the same year several outbreaks in equines were reported in Loulé (Algarve) with some animals testing positive by serology [60].

It is likely that the main reason why none of the human or equine reported infections in Portugal have so far been confirmed by RT-PCR is due to diagnosis and screening being performed well past the viraemic phase. This might also explain why obtaining sequencing data related to WNV Portuguese infections has not been successful. To date, a single study has been published on the local phylogenetics of WNV, with sequencing solely obtained from mosquito pools (Table 4) [58]. Such sequences, based on a short segment of the NS5 gene, demonstrated a close genetic relationship with past WNV lineages circulating in the Mediterranean basin (France, Morocco and Italy), but fell short in providing definite conclusions on viral diversity, phylodynamics and persistence in the country.

When summarized, the past evidence for the circulation of WNV in Portugal is mostly restricted to reports from the southern regions of the country, in particular the districts of Faro, Beja, Setúbal, Évora, Lisboa, Portalegre and Santarém (Figure 1, Table 2). Such reports include evidence from avian species, equines, mosquitoes and humans, achieved via different confirmation methods. It is unknown whether the observed WNV circulation in the southern regions, compared to the northern regions, stems from sampling bias, or from disparate ecological conditions dictating WNV epidemic activity. Sampling bias could be due to large differences in regional biodiversity and population sizes of sentinel species (mosquitoes, birds, equines) across the country. Indeed, there is evidence that the equine population is larger in the south [64], while the population sizes and diversity of birds in the north and south are largely unknown for many resident and migratory species. To the ecological point, climate is a major driver of the population dynamics of mosquito and avian species participating in the zoonotic cycle of WNV. According to the Köppen climate classification, Portugal presents a north-south divide in climate types, with the northern region classified as Temperate Mediterranean, and the southern region as Warm Mediterranean. A lack of past studies exploring the effects of climate on WNV epidemiology in Portugal has, however, not allowed for much speculation on its role in the seemingly restricted circulation of the virus to the south of the country. Considering the possible sampling bias and lack of the above mentioned ecological and genetic information, and with other factors still to be researched in the national context, there is currently no consensus on where and whether WNV circulates in Portugal sporadically or endemically.

• Download figure
• Open in new tab

Figure 1 Past evidence of WNV circulation in Portugal.

Summary of data sources (articles, national reports) on various types of evidence for the circulation of WNV in Portugal. Panel A presents the geographical region of each data source (latitude-longitude). Positioning is only approximate, given that some reports refer to regions of different areas. For data sources referring to large areas (e.g. districts), the capital city of the region was used. Not all sources are mapped, e.g. not reporting where positive samples were found. Mapped words represent the type of evidence reported, aggregated into two categories: molecular, for PCR, isolation; and serology, for HI, VNT, ELISA. Panel B shows a timeline of the number of data sources per year. The color scale (on the right) represents the different types of samples reported, and refers to both panels. Not all years between 1966 and 2020 had data sources. Detailed summary of data sources in Table 2

West Nile virus vectors in Portugal

To date, very few articles in the research literature have provided geo-temporal data on the population size and dynamics of mosquito species relevant for the transmission of WNV in Portugal (Table 3). We obtained REVIVE Culex pipiens data collected between 2017 and 2019 across Portuguese counties (Figure 2). The methods of collection (e.g. type of trap used), as well as resources and regional objectives, varied by year, as described in detail in the published reports [27]. Summarizing the data revealed that both sampling coverage and the total number of yearly collections per county varied in time, both for adult (Figures 2A-C) and immature (Figures 2E-G) forms of this mosquito species. Moreover, we noticed that even in the presence of many project-related factors affecting sampling, Culex pipiens was found to be present across the country (Figures 2DH). Since Culex pipiens is considered the main vector for the possible transmission of WNV in Portugal, these findings stand in contrast with the north-south divide of WNV circulation found in the data review from the literature and public databases (Figure 1, Table 2).

• Download figure
• Open in new tab

Figure 2 REVIVE Culex pipiens counts in Portugal (2017-2019).

Borders represent Portuguese counties, the level at which Culex pipiens data are reported. In panels A-C and E-G, counties are colored according to the yearly observations (transformed to lo g₁₀(counts+ 1)) of adult and immature mosquitoes (respectively). In panels D and H, counties are colored by the yearly mean observations (transformed to lo g₁₀(mean(counts)+1)) of adult and immature mosquitoes (respectively). Counties for which there was zero counts are coloured in white. Counties for which collections were not performed are in grey.

West Nile virus transmission potential in Portugal

In this study we use a climate-driven suitability index (termed index P) to estimate the geo-temporal transmission potential of WNV in Portugal (see Methods for details). This index is the theoretical transmission potential P of each adult female mosquito in the zoonotic reservoir (derived from the decomposition of R₀=MP). P can be interpreted as a proxy for the risk of spillover to the human and equine populations (as demonstrated recently for Israel [41]). Typically, the estimated seasonal variations in index P would be validated by quantifying geo-temporal correlations between time series of the index and suspected or confirmed infections in humans or equines. However, WNV infection time series do not exist for Portugal given the absence of reported chains of transmission (Figure 1, Table 2).

We thus first compared the only available time series data related to WNV transmission in Portugal, which included the newly estimated index P and regional Culex pipiens population sizes from the REVIVE initiative (Figure 3). Given the north-south disparity in past evidence for the circulation of WNV (Figure 1) and the widespread geographical distribution of Culex pipiens (Figure 2), we compared directly two of the most northern and southern districts of the country. Although the index P measures transmission per mosquito and not population size, we found the index to present peaks and troughs coinciding in time with those of the mosquito population sizes, both in the south (Figures 3AC) and north (Figures 3BD). Apart from the previously demonstrated potential for the index P to be correlated with WNV human infections [41], these results thus suggest that (at least in Portugal) the index can also be used to estimate the seasonal timing of variation in Culex pipiens population sizes.

• Download figure
• Open in new tab

Figure 3 Monthly REVIVE Culex pipiens counts in Portugal versus estimated transmission potential index P (2017-2019).

Panels A-D present the monthly time series for log10[1+mosquito population size] (grey bars) and mean (colored lines) and standard deviation (coloured areas) for four districts: two from the south (A for Beja in purple, C for Faro in purple) and two from the north (B for Bragança in red, D for Viana do Castelo in blue). The map panel presents all districts, coloured by their location. Panel E shows the mean index P for all districts 2009-2019 (North in grey and South in black). Panel F shows the seasonal index P per month from panel E independently of year. Panel G is the distribution of index P across the period 2009-2019. For the estimation of the index P, the monthly mean of the climate variables was used per district.

The patterns in Figures 3A-D further suggestde several characteristics of the typical season of WNV transmission in Portugal, as informed by both mosquito population size and estimated transmission potential. Seasons appeared to occur between the months of March and November of each year, often peaking in the summer and reaching minimum transmission potential values during the winter months. We calculated the peak timing for each region using the index P time series of all counties, and found that historically peaks occur in July, although on average ≈ 12 days later in the south compared to the north (t-test p<0.01, see Methods for details). When comparing the northern and southern districts for a longer time period (2009-2019, Figures 3EF), we found that seasonal variation in index P typically presented higher values during the summer months in the north compared to the south (≈ 35 % higher peaks, t-test p<0.001). However, compared to the northern districts, the southern districts presented wider seasonal waves, with transmission potential increasing earlier in the spring and remaining non-zero well into late autumn (Figure 3F). The higher values of P in the north, suggesting higher summer transmission potential compared to the south, contrasted with the existing evidence for WNV circulation (Figure 1). However, the long-term time distribution of index P seasonal waves in the north was skewed, with lower medians than in the south (Figure 3G), potentially supporting the north-south divide found in past evidence for WNV circulation in the country (Figure 1).

To explore in detail the regional differences in transmission potential across the country, we next estimated the index P per available latitude-longitude coordinate in our climate dataset for the time period between 2016 and 2019 (Figure 4). Yearly summary statistics per latitude-longitude coordinate revealed a north-south divide similar to the one presented by past evidence of WNV circulation in the country (Figure 1). In particular, we found the southern region of Portugal to universally present higher median yearly index P when compared to the north (Figures 4ACEG). However, the geographical ‘boundary’ between north and south, in terms of WNV transmission potential, changed over the years (compare e.g. Figures 4A to 4C). This output showed that, in accordance to previously described geo-temporal patterns for other mosquito-borne viruses [34] and WNV [41] in different countries, natural climate variations strongly influence both inter- and intra-yearly variation in transmission potential. When quantifying the number of months within a year for which the estimated transmission potential was not zero, we again found a clear division between the north and south (Figures 4BDFH). Regions in the south and south-west typically spent nine or more months of the year with non-zero transmission potential, contrasting to regions in the north and north-east which spent only half or less of the year with non-zero transmission potential. Both the higher yearly median and longer seasonal waves in the south compared to the north, could thus be seen to support the mentioned north-south divide in past evidence for WNV circulation in Portugal (Figure 1), hinting that the south may provide particularly favourable conditions for the virus.

• Download figure
• Open in new tab

Figure 4 Estimated WNV transmission potential across Portugal (2016-2019).

Maps present summary statistics of estimated index P (WNV transmission potential) at the spatial resolution available in the climate dataset. Panels A, C, E and G present the yearly median index P per coordinate according to the presented color scale on the right. Panels B, D, F and H present the number of months during each year for which the estimated index P per coordinate was above zero, according to the presented color scale on the right. Boundaries in dark grey represent districts.

Long-term WNV suitability and climatic trends in Portugal

Estimating transmission potential of a pathogen via the index P, while informed by climatic variables, provides means to reconstruct historical trends when long-term climatic data is available. We thus used the past 40 years of annual climate data per Portuguese county to explore long-time trends in climate variables and estimated WNV transmission potential, using linear quantile mixed models (see Methods for technical details).

We found that all four explored variables presented significant differences between the north and south regions of the country (Figure 5). Among these, precipitation had the smallest change over time, with a yearly decrease of ≈−1.2× 10⁻⁶ in the north, and increase of ≈ 4.8 × 10⁻⁶ in the south (Figure 5A). Following the Köppen climate classifications, temperature was historically warmer in the south compared to the north. There were similar yearly increases of ≈ 0.041 (north) and ≈ 0.039 (south) Celsius per year between 1981 and 2019, equating to cumulative changes of ≈ 1.55 and ≈ 1.44 degrees for the north and south, respectively (Figure 5B). As the country warmed, it also became dryer, with an accumulated decrease in relative humidity of ≈ 2.8 % in the north and ≈ 1.6 % in the south (Figure 5C). Such historical climate trends translated into an yearly increase in estimated transmission potential for both regions (Figure 5D). The index P had an yearly increase of ≈ 0.00097, accumulating a change of ≈ 0.038 between 1981 and 2019. During this period, the south had an absolute difference of ≈ 0.083 in estimated transmission potential compared to the north.

• Download figure
• Open in new tab

Figure 5 Long-term trends and periodicity of transmission potential of WNV and climate across different regions of Portugal.

Panels A, B, and C present the median of climatic variables among counties (points) within the north (blue) and south (orange) regions defined by the set of districts listed in the map subpanel (on the right). Panel D presents the same as A-C but for the index P (transmission potential). Panel E presents the cross-wavelet transform of the median yearly transmission potential between north and south time-series. Significant periodic components are circled in white (p−value <0.01), where the heat map denotes power with peaks indicated by black dots (wavelet ridges) and mean power throughout time shown in the right subplot.

As described in Figure 3, the northern and southern regions presented an evident yearly periodicity in climate-driven transmission potential, with cycles being stable throughout the years and in phase between the regions. To explore the possibility of periodicity beyond the 12 month cycles, we performed wavelet analysis on the north and south yearly median of transmission potential (see Methods for details). The analysis revealed longer periodic components in both regions, which were seen to change between 1981 and 2019 (Figure 5E). Before 2005, periodicity was mostly found in ≈ 8 year cycles, although a strong 2-3 year component was present around 1995. From 2005 onward, periodicity shifted to become strongly dominated not by one but two components at 2-3 and 4-5 years. The presence of such quasi-periods indicates a Moran effect [65], i.e. that climate synchronizes transmission potential across large geographical areas in Portugal beyond the classic intra-year, summer winter cycles. Finally, the significant shifts in the periodic components over the years do not support an historically stable influence of climate, but are rather suggestive that inter-year instability is an expected feature of WNV transmission potential in Portugal.