Next Article in Journal
California Serogroup Viruses in a Changing Canadian Arctic: A Review
Next Article in Special Issue
The Evolution of Highly Pathogenic Avian Influenza A (H5) in Poultry in Nigeria, 2021–2022
Previous Article in Journal
Genetic Evolution and Variation of Human Adenovirus Serotype 31 Epidemic Strains in Beijing, China, during 2010–2022
Previous Article in Special Issue
Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses, 1988–2021
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 15
Issue 6
10.3390/v15061241
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Remote Sensing and Ecological Variables Related to Influenza A Prevalence and Subtype Diversity in Wild Birds in the Lluta Wetland of Northern Chile

by
Soledad Ruiz
1,2,
Pablo Galdames
2,
Cecilia Baumberger
2,
Maria Antonieta Gonzalez
2,
Camila Rojas
2,
Cristobal Oyarzun
2,
Katherinne Orozco
2,
Cristian Mattar
3,
Pamela Freiden
4,
Bridgette Sharp
4,
Stacey Schultz-Cherry
4,
Christopher Hamilton-West
2,* and
Pedro Jimenez-Bluhm
1,*
1
Escuela de Medicina Veterinaria, Facultad de Ciencias Biológicas, Facultad de Medicina, Facultad de Agronomía e Ingeniería Forestal, Pontificia Universidad Católica de Chile, Santiago 7820436, Chile
2
Departamento de Medicina Preventiva Animal, Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Santiago 8330111, Chile
3
Laboratory for Analysis of the Biosphere (LAB), Universidad de Chile, Santiago 8330111, Chile
4
Department of Infectious Diseases, St. Jude Children’s Research Hospital, Memphis, TN 38105, USA
*
Authors to whom correspondence should be addressed.
Viruses 2023, 15(6), 1241; https://doi.org/10.3390/v15061241
Submission received: 3 May 2023 / Revised: 19 May 2023 / Accepted: 23 May 2023 / Published: 25 May 2023
(This article belongs to the Special Issue Drivers of Evolution of Animal RNA Viruses, Volume II)
Browse Figures
Versions Notes

Abstract

:
The Lluta River is the northernmost coastal wetland in Chile, representing a unique ecosystem and an important source of water in the extremely arid Atacama Desert. During peak season, the wetland is home to more than 150 species of wild birds and is the first stopover point for many migratory species that arrive in the country along the Pacific migratory route, thereby representing a priority site for avian influenza virus (AIV) surveillance in Chile. The aim of this study was to determine the prevalence of influenza A virus (IAV) in the Lluta River wetland, identify subtype diversity, and evaluate ecological and environmental factors that drive the prevalence at the study site. The wetland was studied and sampled from September 2015 to October 2020. In each visit, fresh fecal samples of wild birds were collected for IAV detection by real-time RT-PCR. Furthermore, a count of wild birds present at the site was performed and environmental variables, such as temperature, rainfall, vegetation coverage (Normalized Difference Vegetation Index—NDVI), and water body size were determined. A generalized linear mixed model (GLMM) was built to assess the association between AIV prevalence and explanatory variables. Influenza positive samples were sequenced, and the host species was determined by barcoding. Of the 4349 samples screened during the study period, overall prevalence in the wetland was 2.07% (95% CI: 1.68 to 2.55) and monthly prevalence of AIV ranged widely from 0% to 8.6%. Several hemagglutinin (HA) and neuraminidase (NA) subtypes were identified, and 10 viruses were isolated and sequenced, including low pathogenic H5, H7, and H9 strains. In addition, several reservoir species were recognized (both migratory and resident birds), including the newly identified host Chilean flamingo (Phoenicopterus chilensis). Regarding environmental variables, prevalence of AIV was positively associated with NDVI (OR = 3.65, p < 0.05) and with the abundance of migratory birds (OR = 3.57, p < 0.05). These results emphasize the importance of the Lluta wetland as a gateway to Chile for viruses that come from the Northern Hemisphere and contribute to the understanding of AIV ecological drivers.

1. Introduction

Considering that all influenza A virus (IAV) subtypes have the potential to contribute to the emergence of a pandemic strain, knowledge of infection dynamics and surveillance in wild birds is central to prevent the disease in humans and food-producing animals [1,2]. Wild waterfowl, particularly the Anseriformes and Charadriiformes orders, are considered main reservoirs of IAV in its low pathogenic form (LPAI) [1,3]. These birds have the potential to spread IAV when they migrate within and between continents, representing a risk for the emergence of highly pathogenic avian influenza (HPAI) outbreaks in domestic birds [4,5]. However, while much is known about IAV prevalence in Eurasian and North American wild birds [5,6,7,8], there remains a lack of information on IAV prevalence in wild birds from South America [9,10,11,12,13].
Chile has a coastline of more than 4000 km, along which there are hundreds of wetlands that serve as breeding grounds for resident birds and wintering areas for the many migratory species that arrive each year on the Pacific and Atlantic flyways [1,14]. Recent studies carried out by our research group have demonstrated the presence of a wide diversity of IAV subtypes in wild birds including low pathogenic H5 and H7 strains [9]. Sporadic outbreaks of LPAI have also been found in commercial farms and a widespread circulation of IAV has been found in backyard production systems (BPSs), which could be the result of spillover events from wild birds [15,16,17].
Along the Chilean coastline, there is a wide diversity of climates, which determines the environmental characteristics of each wetland. For example, the Lluta River wetland in Arica and Parinacota Region has a desert climate (<1 mm yearly rainfall); the El Yali wetland in the Valparaíso Region (central zone) has a climate comprising a cloudy steppe and warm temperate steppe with winter rainfall,; and the Itata River in Biobio Region (southern zone) has a warm temperate climate with dry and rainy seasons. This environmental variety of Chilean wetlands and the size of the country may affect the structure and composition of bird assemblages; however, this is still poorly studied in Chile, especially in the northern zone [18,19].
The coastal wetlands of northern Chile, being important sources of water in an extremely arid area, are ecosystems that shelter and serve as refuge for many species of wild birds [19]. The Lluta River wetland is the northernmost coastal wetland in Chile and belongs to the biogeographical province of the Coastal Desert. Its main water source is the Lluta River, the only exoreic river in the province, which experiences its most important floods during the summer season (December to March) due to the “highland winter” rains. The average annual temperature fluctuates between 16° and 22 °C and rainfall is almost non-existent throughout the year, reaching an annual maximum of 1.5 mm [19,20].
This wetland represents the main habitat of the local avifauna in the Atacama Desert and the gateway to the country for many migratory species that use the Pacific route [19,21]. The movements of those species may play an important role in the maintenance and spread of North American lineage viruses; therefore, this wetland may be a possible point of confluence where North and South American IAVs intermix [9], thereby representing a priority site for avian influenza virus surveillance in Chile.
However, to understand the role of wild birds in the spread and persistence of IAV, surveillance efforts must focus not only on the isolation and characterization of circulating viruses at sites where these birds concentrate; it is also critical that environmental factors that may impact IAV prevalence in these ecosystems be considered [22]. These factors include abiotic elements, such as rainfall and temperature; biotic factors, such as bird species present and amount of vegetation; and anthropogenic factors, such as land use and land cover types [22,23,24,25,26]. However, to date, studies that include ecological aspects of the virus and that explore the interaction between IAV prevalence in wild birds and the environment remain scarce [12,22,24,25,27]. So far, the only study performed in Chile that explores ecological and environmental factors related to the prevalence of IAV in wild birds demonstrated that the prevalence of IAV in central Chilean wetlands is not constant over the year, being higher in the summer and fall months. Environmental factors driving this prevalence included minimum temperature, Normalized Difference Vegetation Index (NDVI), and size of the wetland water body [12].
The aim of this study was to determine the prevalence of IAV in the Lluta River wetland between September 2015 and October 2020, identify subtype diversity, and evaluate ecological and environmental factors that drive these variables at the study site.

2. Materials and Methods

2.1. Study Area

The Lluta wetland (S 18°24′55″–W 70°19′20″) geographically belongs to the Lluta River watershed, the only exoreic river in the province of Arica in the Region of Arica and Parinacota. It is located on the coastal border, 10 km north of downtown Arica and 8 km from the border with Peru. It has an area of approximately 300 ha, of which only 31 ha is protected under the denomination of “Municipal Nature Reserve”. For birds, it is the most important wetland in northern Chile, hosting more than 150 species (both resident and migratory), which represent 32% of the bird species described for the country [19,20,21]. In 2010, it was designated a shorebird reserve by the Western Hemisphere Shorebird Reserve Network (WHSRN), a conservation initiative for the protection of the ecological integrity of critical habitats for shorebirds throughout the Americas [20].
This wetland corresponds to the first stopover point for many migratory species that arrive in the country along the Pacific migratory route, such as Franklin’s gull (Leucophaeus pipixcan), Whimbrel (Numenius phaeopus), Least sandpiper (Calidris minutilla), Baird’s sandpiper (Calidris bairdii), Lesser Yellowlegs (Tringa flavipes), and Greater Yellowlegs (Tringa melanoleuca), which congregate in high densities during the spring and summer months [14,19,21].

2.2. Sampling and Sample Size

Sampling was carried out between September 2015 and October 2020, with a total of 19 sampling campaigns performed during that period. Samples were collected at irregular intervals. In each visit, fresh fecal samples of wild birds were collected for IAV detection. To minimize the probability of sampling the same individual’s feces twice, samples were collected uniformly by line transects throughout the area where a flock was previously observed. Limited samples were taken per each flock, with several sampling points in the site [12]. Samples were collected using single-use sterile flocked swabs (Copan®, Brescia, Italy) and stored in tubes containing 1 mL of universal transport medium (UTM, Copan®, Brescia, Italy). Samples were maintained at 4 °C until arrival at the laboratory of Veterinary Epidemiology of the School of Veterinary Sciences of the University of Chile. Samples were stored at −80 °C until processing. A minimum sample size of 178 samples per sampling occasion was calculated and estimated from the sample size calculation formula for detecting the presence of disease for finite populations (Equation (1)) [28]. A population size of 1000 fresh feces at the site was considered. In addition, based on results from a previous study [9], a minimum expected prevalence of 1.5% with a significance level of 5% was considered.
Equation (1):
n = ( 1 − ( α ) 1 D ) ( N − D − 1 2 ) ,
where:
  • n = required sampling size;
  • N = population size;
  • α = 1-confidence level;
  • D = estimated minimum number of positive samples.

2.3. Influenza A virus Detection

The collected samples were processed individually, and RNA extraction and real-time RT-PCR (RT-qPCR) analysis were performed at the Laboratory of Veterinary Epidemiology. Briefly, viral RNA extraction was performed on 50 μL of swab sample using the Ambion MagMax-96 AI/ND viral isolation kit (Life Technologies Corporation, Grand Island, NY, USA). Sample screening was performed using RT-qPCR (Stratagene mx3000 p Santa Clara, CA, USA) with the TaqMan Fast Virus 1-Step Master Mix (Applied Biosystems, Foster City, CA, USA). Specific primer/probes for the influenza M gene were used for the RT-qPCR reaction as previously described [29]. Samples with a cycle threshold value (Ct) ≤ 38 were considered positive [30]. Swab samples that tested positive by RT-qPCR with a Ct ≤ 35 were inoculated in embryonated chicken eggs to virus isolation as previously described [31]. Host identification was attempted for IAV-positive fecal samples using primers designed to amplify a segment of the mitochondrial cytochrome-oxidase I as described elsewhere [32].

2.4. Sequencing and Phylogenic Analysis

Sample sequencing was carried out at the St. Jude Children’s Hospital Hartwell Center on the Illumina MiSeq sequencer platform (Illumina, San Diego, CA, USA). Libraries were prepared using the Nextera XT DNA-Seq library preparation kit, as previously described [33]. Reads were de novo assembled using the SAPDes package [34]. Reference sequences for the phylogenetic analysis were obtained from the Influenza Virus Database [35]. Representative sequences from North America, Europe, and Asia were included in each tree. Briefly, sequence editing was performed using BioEdit version 7.2.5 [36] and sequence alignment was performed with MUSCLE version 3.8.3 [37]. Phylogenetic inference was undertaken for the HAs of the H5, H7, and H9 isolates by maximum likelihood, with a general time-reversible model and gamma-distributed rate variation among sites. Maximum likelihood analysis was performed using RAxML version 8.0 [38]. To provide statistical robustness of each node, 1000 bootstrap resampling processes were performed. Trees were visualized on FigTree (v.1.4.3). Genbank accession numbers are KX185901, KX185918, MH499035, MH498968, MH498978, MH499057, MH499142, MK164009, OP888556, and OQ820949.

2.5. Ecological and Environmental Variables

The ecological and environmental variables evaluated in this study were grouped into three categories (Table 1).
Wild bird community: At each visit and prior to sampling, a count of wild birds present at the site was performed during morning hours. A point-counting approach with several experienced observers was used [12]. These data were used to estimate total bird abundance, species richness, and abundance of migratory birds present at the time of sampling.
Landscape variables: The vegetative cover of the wetland and the size of the water body was measured for each sampling month, using images of the LANDSAT 7-ETM and 8-OLI satellites. In addition, to assess the cumulative effect of landscape characteristics and IAV prevalence at the site, these variables were also collected one month, two months, and three months before the sampling month [12]. Satellite images were downloaded from the website of the United States Geological Survey (USGS EarthExplorer at https://earthexplorer.usgs.gov/) and processed using ENVI® software version 4.7. The Normalized Difference Vegetation Index (NDVI) was used as an indicator of the vegetation cover; to delineate water bodies, the Modified Normalized Difference Water Index (MNDWI) was used [39,40].
Meteorological data: Monthly average temperatures, humidity, and cumulative rainfall were measured for each sampling month. To investigate the cumulative effect of rainfall, humidity, and temperature on IAV prevalence, these variables were also collected one month, two months, and three months beforehand [12]. The information was obtained from the Lluta Bajo weather station, which belongs to Agrometeorological Network (Agromet) of the Agricultural Research Institute (INIA) (https://agrometeorologia.cl/).

2.6. Statistical Analysis

To determine the association between the prevalence of IAV in the wetland and the ecological and environmental explanatory variables defined in Table 1, a GLMM was used. The response variable was defined as the number of IAV-positive samples obtained in each sampling in relation to sample size (prevalence). The sampling occasion was entered as a random effect in the model [12]. The unconditional association between IAV prevalence and each of the recorded explanatory variables was assessed in the first bivariate GLMM (Supplemental Table S1). Variables associated with the outcome at a liberal p-value of <0.15 were selected for inclusion in the multivariable model.
The linearity of continuous explanatory variables against the log odds of IAV positivity was assessed visually. Non-linear variables were categorized using the median. A forward stepwise inclusion of variables, guided by the Akaike Information Criterion (AIC), was performed to build the final multivariable mixed model. All analyses were performed using R statistical software using the “glmer” function of the “lme4” package [41]. The significance level was set at 5%.

3. Results

3.1. Influenza Prevalence and Species Richness

Between September 2015 and October 2020, a total of 4349 fecal samples were collected, of which 90 (2.07%, 95% CI: 1.68 to 2.55) were positive for the influenza virus M gene according to RT-qPCR. Thirty-five of these samples were inoculated on eggs from which the isolates were obtained (Supplemental Table S2).
Over the five years of sampling, the prevalence of IAV varied widely (between 0% and 8.6%), showing apparent peaks during the months when wild migratory birds were present (Figure 1). However, because sampling was carried out at irregular intervals, it was not possible to determine seasonality in this prevalence.
Bird numbers fluctuated widely among sampling occasions over the study period (Figure 2A), where a greater number of individuals were observed during the summer and early autumn months, mainly due to the arrival of migratory birds at the site, highlighting species such as Franklin’s gull (Leucophaeus pipixcan), Whimbrel (Numenius phaeopus), and many species of shorebirds.
Species richness observed during the study months ranged from a minimum of 17 to a maximum of 44, showing an increase during the spring and summer months with the arrival of migratory species from the Northern Hemisphere (Figure 2B). During the study period, a total of 71 species were identified in the wetland, of which 41 are considered resident, 21 migratory, and 9 accidental (Supplemental Table S3).

3.2. Environmental Variables

Meteorological variables (precipitation, humidity, maximum and minimum temperatures) did not show large variations during the months of study (Supplemental Table S4) and none of them showed significance with the prevalence of IAV in the wetland.
Regarding environmental variables, the amount of water in the wetland increased between January and March and reached its lowest level between October and December, while the vegetation cover (measured by NDVI) was higher between February and March and reached its lowest level in the months of October and November. Figure 3 shows the months in which the greatest variations in NDVI (lowest and highest) occurred during the study period.
The final multivariable model showed that IAV prevalence was positively associated with NDVI for the month of sampling and with the abundance of migratory birds (Table 2). The months in which 113 or more individuals belonging to migratory species were counted were associated with higher prevalence (OR = 3.57, p < 0.05). With respect to vegetation, the months in which NDVI values were equal to or greater than 0.27 were associated with higher prevalence (OR = 3.65, p < 0.05).
None of the variables included in the “Meteorological data” category were retained in the final multivariable model.

3.3. Subtype Diversity and Sequence Analysis

A wide diversity of hemagglutinin (HA) and neuraminidase (NA) subtypes were identified, and 10 viruses were isolated, including low pathogenic H5 and H7 strains (Table 3). Six reservoir species were recognized, including Larus dominicanus (Kelp gull, n = 2), Haematopus palliatus (American oystercatcher, n = 3), Haematopus ater (Blackish oystercatcher, n = 1), Larus pipixcan (Franklin’s gull, n = 2), Pluvialis squatarola (Grey plover, n = 1), and newly identified host Phoenicopterus chilensis (Chilean flamingo, n = 1). All obtained viruses were LPAI according to HA cleavage sequence analysis. Results of the phylogenetic analysis performed on H5, H7, and H9 hemagglutinins (Figure 4), revealed that the H5 hemagglutinin was most like other H5 hemagglutinins obtained in the central region of Chile between 2015 and 2017. Conversely, the H7 hemagglutinin belonged to the North American clade of viruses, with similar sequences obtained from birds belonging to the Mississippi and Pacific flyways. Additionally, H9 viruses did not form a distinct clade within the global H9 phylogeny but were closely related to each other and to other isolates obtained in North America.

4. Discussion

This study provides information on the prevalence of IAV and the environmental factors that could influence its occurrence in a wetland of great importance for IAV surveillance in Chile. In addition, this wetland has been a key site for surveillance of the current HPAI H5N1 2.3.4.4b outbreak that is affecting several countries in South America, and is the site where the first positive case was detected in Chile in December 2022 [42,43]. Therefore, our results provide important considerations for structuring surveillance and early warning actions in a wetland that is the gateway to the viruses from North America.
The overall IAV prevalence identified at the Lluta wetland during the months sampled was 2.07% (95% CI: 1.68–2.55), which is lower than that recorded in sites in the central zone of Chile during the same study period (4.28%) [12]. This could be due to the difference between the bird assemblages in both areas. While in Lluta wetland the predominant species correspond to shorebirds, in the wetlands of central Chile there are many resident ducks that have already been identified as primary hosts of IAV viruses [9] and may have an important role in the perpetuation of the virus throughout the year [12]. In another previously published study that included sampling in the Lluta wetland, the prevalence found was also higher (5.15%); however, the aforementioned prevalence corresponds to a single sampling event [9].
Regarding subtype diversity and origins, this is the first time that an H7 virus originating from North America has been found in the northern part of Chile, and its presence may have significant implications. Unlike the South American clade of H7 viruses, which has only been associated with HPAI in South America once [44], North American strains have been responsible for numerous outbreaks of HPAI [45,46]. On the other hand, the H5 and H9 HAs appear to be of South American origin, although further genomic surveillance is necessary to definitively establish their clade formation.
With respect to the ecological variables evaluated in this study, only the abundance of migratory birds was positively related to the prevalence of IAV at the site. None of the other variables in the “bird community” category (species richness and total bird abundance) showed statistical association with prevalence. This is probably because bird abundance also increased in months where prevalence was low (autumn and winter), due to large numbers of birds that are considered resident at the site, such as the Grey gull (Leucophaeus modestus). On the other hand, although species richness has been considered a factor influencing the prevalence of the virus in studies carried out in Spain [25,27], other studies around the world have found that bird density or the abundance of migratory species have a greater influence on the prevalence of AIV than bird richness [22,47,48], which is in agreement with the present results obtained in Lluta wetland.
The number of migratory birds in the Lluta River wetland increases considerably between October and March, with high densities at the site, and then declines to almost zero during the rest of the year. However, there are variations in the abundance of migratory birds detected in the same months between different years (Figure 2A). This could be explained because migratory patterns are driven by environmental conditions (such as temperature, food availability, and landscape structure) [49,50], which may vary from year to year or be influenced by climate change [51] or by climatic phenomena such as the El Niño–Southern Oscillation (ENSO), the main driver of interannual climate extremes in South America [52].
The species of migratory birds identified belong to the families Laridae, Scolopacidae, and Charadriidae, which have been identified as important reservoirs of IAV [1,3]. Franklin’s gull (Leucophaeus pipixcan), Whimbrels (Numenius phaeopus), and Sanderling (Calidris alba) were the species with the highest number of individuals during the sampling period (Supplemental Table S1). In addition, it is important to highlight that three viral isolates in our study were found in migratory species: two in Franklin’s gulls (H13N9) and one H9N7 subtype in a Grey plover (Pluvialis squatarola). Recent studies in Chile and Peru have already described the infection of IAV in Franklin’s gulls [9,53,54], Grey plover [9], and Whimbrels [9,10,55]. In one of these studies, which determined a wide diversity of IAVs in wetlands in Chile, it was also found that the viruses isolated in Lluta were the most genetically diverse, with segments of viruses from North American and South American lineages [9], confirming the thesis of the importance of migratory birds in the dissemination of IAV and in the emergence of new subtypes.
With respect to the meteorological variables evaluated in this study, none were significantly associated with IAV positivity to the virus in the wetland, although climatic variables such as temperature and precipitation have been described as important drivers of IAV in wild birds in other regions of the world [24,26,56,57]. This could be explained mainly because climatic conditions in the wetland are stable throughout the year. Precipitation in the area is almost non-existent and the average annual temperature is around 18 °C [21]. The same is observed in wetlands in tropical Africa, where climatic factors were poorly related to IAV prevalence in wild birds, but host ecological factors (such as bird density in wetlands and the arrival of migrants from Eurasia) played a much more important role [22]. This highlights the relevance of conducting studies that include environmental factors in a local context.
Regarding landscape variables, vegetative coverage in the sampling month (measured by NDVI) was positively associated with the prevalence of IAV in the wetland. During the study period, NDVI values fluctuated between 0.18 and 0.31, which is related to the presence of shrub and herbaceous vegetation [58] consisting mainly of communities of horsetail, chilca shrubs, salt grass, and reed [19]. The months in which NDVI values were equal to or greater than 0.27 were associated with higher prevalence. These results are similar to those found in Europe, where vegetation surrounding wetlands has been recognized as an important driver in the presentation dynamics of the IAV [25,26], because wild birds are impacted by the availability of food resources and shelter provided by a wetland. Higher NDVI may indicate the presence of better vegetative food resources and shelter, favoring a higher concentration of birds and therefore a higher risk of pathogen transmission [25,59]. In another study conducted in the Poyang Lake wetland in China, it was determined that there is a positive correlation coefficient between NDVI and the number of birds present in the wetland, especially migratory birds [59]. In the Lluta wetland, NDVI reaches its maximum between October and March. These months coincide with the arrival of migratory birds from the Northern Hemisphere and with a greater growth of vegetation in spring, which extends into the summer mainly due to the increase in the flow of the Lluta River due to the “highland winter” rains (December to March).
With respect to water size, this variable was not significant in the prevalence of AIV in the Lluta wetland. These results were different from those found in central Chile where the extreme variations in the water body size experienced by the wetlands during summer months directly affect the congregation of birds and therefore the dynamics of infection of the IAV [12]. During the study period, the Lluta wetland did not show significant variations in the surface area flooded by water. Only during January 2018 was a larger area of water observed in the wetland compared to the other months sampled. However, because the area of the wetland that is flooded with water is relatively small (during the sampling months it fluctuated between 0.5 and 2 ha), the concentration of hosts in those areas, especially during the times of migrant arrivals, generates high densities of individuals favoring the increase in prevalence at the site due to density-dependent transmission of the agent. This is similar to what is observed in wetlands in tropical Africa [60].
With respect to the IAV host species identified, in addition to the two migratory species mentioned above, IAV isolates were found in American and blackish oystercatcher (Haematopus palliates, Haematopus ater), Kelp gull (Larus dominicanus), and a Chilean flamingo (Phoenicopterus chilensis). Although the first three species have already been described as hosts in other studies in Chile and South America [9,10,11,61,62], this is the first time that IAV has been isolated in a Chilean flamingo. This species corresponds to a neotropical migrant that is widely distributed throughout South America, even reaching the Atlantic coast [63]. Therefore, it may be an important species for the spread of IAV between South American countries.
Although our results contribute to understanding the ecological drivers of IAV in the Lluta River wetland, the study had some limitations. Because sampling was conducted at irregular intervals, it was not possible to determine a seasonal pattern in this prevalence. Sampling was carried out in this way, because the Lluta wetland is in a remote area, where, for logistical reasons, it was not possible to carry out a more systematic sampling. In addition, the last year of the study (2020) coincided with the SARS-CoV-2 pandemic. Therefore, we performed only one sampling during that year.
However, due to the epidemiological importance of this wetland for IAV surveillance in Chile, it would be beneficial to extend this long-term longitudinal study. Ideally, the study would include a greater number of samples at regular intervals to determine whether there is a seasonal pattern in the prevalence of the virus at the site. In addition, this would allow us to understand the difference in prevalence found in the same months between different study years and would allow a better evaluation of the risk associated with the transmission of the virus from wildlife reservoirs to domestic birds. This would also improve surveillance programs by directing sampling efforts towards seasons that are more favorable for the maintenance and transmission of the virus.

5. Conclusions

Our results demonstrate that there is a wide circulation of IAV in wild birds of the Lluta River wetland and a wide range of different influenza subtypes present in the bird populations, including ones considered as “high-risk”. Ten influenza viruses were isolated and sequenced, and this is the first time that an H7 virus having a North American origin was found in the northern part of Chile.
The prevalence of IAV was positively associated with NDVI and with the abundance of migratory birds in the wetland. However, additional long-term studies are needed to include other host-related variables (such as age structure of the population and density of individuals) and to evaluate whether there is a seasonal pattern of prevalence in this wetland.
Our results constitute a unique contribution to the understanding of the ecological drivers that may modulate the occurrence of IAV in northern Chile, providing important considerations for the global surveillance of IAV in wild birds and highlighting the importance of the Lluta wetland as a gateway to Chile for viruses that come from the Northern Hemisphere.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15061241/s1. Supplemental Table S1: Bivariable associations between IAV prevalence and explanatory variables; Supplemental Table S2: Sampling effort detail. Total number of samples collected in the study period, number of positive samples, samples inoculated in eggs and successfully sequenced; Supplemental Table S3. Species identified in the Lluta River wetland between September 2015 and October 2020; Supplemental Table S4: Ranges of meteorological and landscape variables during the sampling months.

Author Contributions

Conceptualization, S.R., C.M., S.S.-C., C.H.-W. and P.J.-B.; Methodology, P.G., C.B., M.A.G., C.R., C.O., K.O., C.M., P.F., B.S., C.H.-W. and P.J.-B.; Validation, C.M.; Formal analysis, S.R. and P.J.-B.; Investigation, S.R. and K.O.; Resources, P.G., C.B., M.A.G., C.R., C.O., K.O., P.F. and B.S.; Data curation, B.S.; Writing – original draft, S.R., C.H.-W. and P.J.-B.; Writing – review & editing, S.R., C.M., S.S.-C., C.H.-W. and P.J.-B.; Supervision, C.H.-W. and P.J.-B.; Project administration, P.J.-B.; Funding acquisition, S.S.-C., C.H.-W. and P.J.-B. All authors have read and agreed to the published version of the manuscript.

Funding

This research project was financed by the National Institute of Allergy and Infectious Diseases grant HHSN272201400006C and 75N93021C00016 to SSC., Fondecyt grant 1191747 to CHW, and Fondecyt grant 11190755 to PJB.

Institutional Review Board Statement

The animal study protocol was approved by Institutional Animal Care and Use Committee (CICUA) of the University of Chile 19265-VET-UCH.

Informed Consent Statement

Not appliable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We would like to thank the Municipality of Arica for giving us unrestricted access to the Lluta wetland and for logistical support during sampling.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Olsen, B.; Munster, V.J.; Wallensten, A.; Waldenström, J.; Osterhaus, A.D.; Fouchier, R.A. Global patterns of influenza A virus in wild birds. Science 2006, 312, 384–388. [Google Scholar] [CrossRef] [PubMed]
  2. Webster, R.; Hulse, D. Microbial adaptation and change: Avian influenza. Rev. Sci. Technol.-Off. Int. Epizoot. 2004, 23, 453–466. [Google Scholar] [CrossRef] [PubMed]
  3. Webster, R.G.; Bean, W.J.; Gorman, O.T.; Chambers, T.M.; Kawaoka, Y. Evolution and ecology of influenza A viruses. Microbiol. Rev. 1992, 56, 152–179. [Google Scholar] [CrossRef]
  4. Caliendo, V.; Lewis, N.; Pohlmann, A.; Baillie, S.; Banyard, A.; Beer, M.; Brown, I.; Fouchier, R.; Hansen, R.; Lameris, T. Transatlantic spread of highly pathogenic avian influenza H5N1 by wild birds from Europe to North America in 2021. Sci. Rep. 2022, 12, 11729. [Google Scholar] [CrossRef] [PubMed]
  5. Zhang, G.; Li, B.; Raghwani, J.; Vrancken, B.; Jia, R.; Hill, S.C.; Fournié, G.; Cheng, Y.; Yang, Q.; Wang, Y. Bidirectional movement of emerging H5N8 avian influenza viruses between Europe and Asia via migratory birds since early 2020. Mol. Biol. Evol. 2023, 40, msad019. [Google Scholar] [CrossRef]
  6. Gorsich, E.E.; Webb, C.T.; Merton, A.A.; Hoeting, J.A.; Miller, R.S.; Farnsworth, M.L.; Swafford, S.R.; DeLiberto, T.J.; Pedersen, K.; Franklin, A.B. Continental-scale dynamics of avian influenza in US waterfowl are driven by demography, migration, and temperature. Ecol. Appl. 2021, 31, e2245. [Google Scholar] [CrossRef] [PubMed]
  7. Hoye, B.J.; Munster, V.J.; Nishiura, H.; Klaassen, M.; Fouchier, R.A. Surveillance of wild birds for avian influenza virus. Emerg. Infect. Dis. 2010, 16, 1827. [Google Scholar] [CrossRef]
  8. Verhagen, J.H.; Fouchier, R.A.; Lewis, N. Highly pathogenic avian influenza viruses at the wild–domestic bird interface in Europe: Future directions for research and surveillance. Viruses 2021, 13, 212. [Google Scholar] [CrossRef]
  9. Jiménez-Bluhm, P.; Karlsson, E.A.; Freiden, P.; Sharp, B.; Di Pillo, F.; Osorio, J.E.; Hamilton-West, C.; Schultz-Cherry, S. Wild birds in Chile Harbor diverse avian influenza A viruses. Emerg. Microbes Infect. 2018, 7, 1–4. [Google Scholar] [CrossRef]
  10. Nelson, M.I.; Pollett, S.; Ghersi, B.; Silva, M.; Simons, M.P.; Icochea, E.; Gonzalez, A.E.; Segovia, K.; Kasper, M.R.; Montgomery, J.M. The genetic diversity of influenza A viruses in wild birds in Peru. PLoS ONE 2016, 11, e0146059. [Google Scholar]
  11. Rimondi, A.; Gonzalez-Reiche, A.S.; Olivera, V.S.; Decarre, J.; Castresana, G.J.; Romano, M.; Nelson, M.I.; van Bakel, H.; Pereda, A.J.; Ferreri, L. Evidence of a fixed internal gene constellation in influenza A viruses isolated from wild birds in Argentina (2006–2016). Emerg. Microbes Infect. 2018, 7, 1–13. [Google Scholar] [CrossRef]
  12. Ruiz, S.; Jimenez-Bluhm, P.; Di Pillo, F.; Baumberger, C.; Galdames, P.; Marambio, V.; Salazar, C.; Mattar, C.; Sanhueza, J.; Schultz-Cherry, S. Temporal dynamics and the influence of environmental variables on the prevalence of avian influenza virus in main wetlands in central Chile. Transbound. Emerg. Dis. 2021, 68, 1601–1614. [Google Scholar] [CrossRef]
  13. Xu, K.; Ferreri, L.; Rimondi, A.; Olivera, V.; Romano, M.; Ferreyra, H.; Rago, V.; Uhart, M.; Chen, H.; Sutton, T. Isolation and characterization of an H9N2 influenza virus isolated in Argentina. Virus Res. 2012, 168, 41–47. [Google Scholar] [CrossRef] [PubMed]
  14. García-Walther, J.; Senner, N.; Norambuena, H.; Schmitt, F. Atlas De Las Aves Playeras de Chile: Sitios Importantes Para Su Conservación; Universidad Santo Tomás: Santiago, Chile, 2017. [Google Scholar]
  15. Bravo-Vasquez, N.; Di Pillo, F.; Lazo, A.; Jiménez-Bluhm, P.; Schultz-Cherry, S.; Hamilton-West, C. Presence of influenza viruses in backyard poultry and swine in El Yali wetland, Chile. Prev. Vet. Med. 2016, 134, 211–215. [Google Scholar] [CrossRef] [PubMed]
  16. Jimenez-Bluhm, P.; Bravo-Vasquez, N.; Torchetti, M.K.; Killian, M.L.; Livingston, B.; Herrera, J.; Fuentes, M.; Schultz-Cherry, S.; Hamilton-West, C. Low pathogenic avian influenza (H7N6) virus causing an outbreak in commercial Turkey farms in Chile. Emerg. Microbes Infect. 2019, 8, 479–485. [Google Scholar] [CrossRef]
  17. Jimenez-Bluhm, P.; Di Pillo, F.; Bahl, J.; Osorio, J.; Schultz-Cherry, S.; Hamilton-West, C. Circulation of influenza in backyard productive systems in central Chile and evidence of spillover from wild birds. Prev. Vet. Med. 2018, 153, 1–6. [Google Scholar] [CrossRef] [PubMed]
  18. Muñoz-Pedreros, A.; Merino, C. Diversity of aquatic bird species in a wetland complex in southern Chile. J. Nat. Hist. 2014, 48, 1453–1465. [Google Scholar] [CrossRef]
  19. Sielfeld, W.; Peredo, R.; Fuentes, R.; Malinarich, V.; Olivares, F. Coastal Wetlands of Northern Chile. In The Ecology and Natural History of Chilean Saltmarshes; Springer: Berlin/Heidelberg, Germany, 2017; pp. 105–167. [Google Scholar]
  20. Navarro, N.; Abad, M.; Bonnail, E.; Izquierdo, T. The Arid Coastal Wetlands of Northern Chile: Towards an Integrated Management of Highly Threatened Systems. J. Mar. Sci. Eng. 2021, 9, 948. [Google Scholar] [CrossRef]
  21. Peredo, R. La Desembocadura Del Río Lluta: Un Humedal Para Las Aves, En El Desierto Costero De Chile. 2007. Available online: https://biblioteca.cehum.org/ (accessed on 9 December 2022).
  22. Gaidet, N.; Caron, A.; Cappelle, J.; Cumming, G.S.; Balança, G.; Hammoumi, S.; Cattoli, G.; Abolnik, C.; Servan de Almeida, R.; Gil, P. Understanding the ecological drivers of avian influenza virus infection in wildfowl: A continental-scale study across Africa. Proc. R. Soc. B Biol. Sci. 2012, 279, 1131–1141. [Google Scholar] [CrossRef]
  23. Cumming, G.S.; Abolnik, C.; Caron, A.; Gaidet, N.; Grewar, J.; Hellard, E.; Henry, D.A.; Reynolds, C. A social–ecological approach to landscape epidemiology: Geographic variation and avian influenza. Landsc. Ecol. 2015, 30, 963–985. [Google Scholar] [CrossRef]
  24. Ferenczi, M.; Beckmann, C.; Warner, S.; Loyn, R.; O’riley, K.; Wang, X.; Klaassen, M. Avian influenza infection dynamics under variable climatic conditions, viral prevalence is rainfall driven in waterfowl from temperate, south-east Australia. Vet. Res. 2016, 47, 23. [Google Scholar] [CrossRef] [PubMed]
  25. Pérez-Ramírez, E.; Acevedo, P.; Allepuz, A.; Gerrikagoitia, X.; Alba, A.; Busquets, N.; Díaz-Sánchez, S.; Álvarez, V.; Abad, F.X.; Barral, M. Ecological factors driving avian influenza virus dynamics in Spanish wetland ecosystems. PLoS ONE 2012, 7, e46418. [Google Scholar] [CrossRef] [PubMed]
  26. Si, Y.; Wang, T.; Skidmore, A.K.; de Boer, W.F.; Li, L.; Prins, H.H. Environmental factors influencing the spread of the highly pathogenic avian influenza H5N1 virus in wild birds in Europe. Ecol. Soc. 2010, 15, 26. [Google Scholar] [CrossRef]
  27. Torrontegi, O.; Alvarez, V.; Acevedo, P.; Gerrikagoitia, X.; Höfle, U.; Barral, M. Long-term avian influenza virus epidemiology in a small Spanish wetland ecosystem is driven by the breeding Anseriformes community. Vet. Res. 2019, 50, 4. [Google Scholar] [CrossRef] [PubMed]
  28. Dohoo, I.R.; Martin, W.; Stryhn, H.E. Veterinary Epidemiologic Research; University of Prince Edward Island: Charlottetown, PE, Canada, 2003. [Google Scholar]
  29. WHO. CDC protocol of realtime RTPCR for swine influenza A (H1N1). In CDC Protocol of Realtime RTPCR for Swine Influenza A (H1N1); World Health Organization (WHO): Geneva, Switzerland, 2009. [Google Scholar]
  30. Shu, B.; Wu, K.-H.; Emery, S.; Villanueva, J.; Johnson, R.; Guthrie, E.; Berman, L.; Warnes, C.; Barnes, N.; Klimov, A. Design and performance of the CDC real-time reverse transcriptase PCR swine flu panel for detection of 2009 A (H1N1) pandemic influenza virus. J. Clin. Microbiol. 2011, 49, 2614–2619. [Google Scholar] [CrossRef] [PubMed]
  31. Moresco, K.A.; Stallknecht, D.E.; Swayne, D.E. Evaluation and attempted optimization of avian embryos and cell culture methods for efficient isolation and propagation of low pathogenicity avian influenza viruses. Avian Dis. 2010, 54, 622–626. [Google Scholar] [CrossRef]
  32. Cheung, P.P.; Leung, Y.C.; Chow, C.-K.; Ng, C.-F.; Tsang, C.-L.; Wu, Y.-O.; Ma, S.-K.; Sia, S.-F.; Guan, Y.; Peiris, J.M. Identifying the species-origin of faecal droppings used for avian influenza virus surveillance in wild-birds. J. Clin. Virol. 2009, 46, 90–93. [Google Scholar] [CrossRef]
  33. Kaplan, B.S.; Russier, M.; Jeevan, T.; Marathe, B.; Govorkova, E.A.; Russell, C.J.; Kim-Torchetti, M.; Choi, Y.K.; Brown, I.; Saito, T. Novel highly pathogenic avian A (H5N2) and A (H5N8) influenza viruses of clade 2.3. 4.4 from North America have limited capacity for replication and transmission in mammals. Msphere 2016, 1, e00003-16. [Google Scholar] [CrossRef]
  34. Prjibelski, A.; Antipov, D.; Meleshko, D.; Lapidus, A.; Korobeynikov, A. Using SPAdes de novo assembler. Curr. Protoc. Bioinform. 2020, 70, e102. [Google Scholar] [CrossRef]
  35. Bao, Y.; Bolotov, P.; Dernovoy, D.; Kiryutin, B.; Zaslavsky, L.; Tatusova, T.; Ostell, J.; Lipman, D. The influenza virus resource at the National Center for Biotechnology Information. J. Virol. 2008, 82, 596–601. [Google Scholar] [CrossRef]
  36. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. In Nucleic Acids Symposium Series; Oxford University Press: Oxford, UK; pp. 95–98.
  37. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed]
  38. Jones, P.; Binns, D.; Chang, H.-Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G. InterProScan 5: Genome-scale protein function classification. Bioinformatics 2014, 30, 1236–1240. [Google Scholar] [CrossRef] [PubMed]
  39. Rouse, J.W., Jr.; Haas, R.H.; Schell, J.; Deering, D. Monitoring the Vernal Advancement and Retrogradation (Green Wave Effect) of Natural Vegetation; National Space Agency: Greenbelt, MD, USA, 1973. [Google Scholar]
  40. Xu, H. Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery. Int. J. Remote Sens. 2006, 27, 3025–3033. [Google Scholar] [CrossRef]
  41. Team, R.C. R: A Language and Environment for Statistical Computing. 2018. Available online: http://www.r-project.org/index.html (accessed on 7 December 2022).
  42. Empres-i, F. Diseases. 2022. Available online: https://empres-i.apps.fao.org/diseases (accessed on 25 November 2022).
  43. Wahis. Sistema Mundial de Información Zoosanitaria. 2022. Available online: https://wahis.woah.org/#/home (accessed on 25 November 2022).
  44. Jones, Y.L.; Swayne, D.E. Comparative pathobiology of low and high pathogenicity H7N3 Chilean avian influenza viruses in chickens. Avian Dis. 2004, 48, 119–128. [Google Scholar] [CrossRef]
  45. Youk, S.; Leyson, C.; Killian, M.L.; Torchetti, M.K.; Lee, D.-h.; Suarez, D.L.; Pantin-Jackwood, M.J. Evolution of the North American Lineage H7 Avian Influenza Viruses in Association with H7 Virus’s Introduction to Poultry. J. Virol. 2022, 96, e00278-22. [Google Scholar] [CrossRef]
  46. Dietze, K.; Graaf, A.; Homeier-Bachmann, T.; Grund, C.; Forth, L.; Pohlmann, A.; Jeske, C.; Wintermann, M.; Beer, M.; Conraths, F.J. From low to high pathogenicity—Characterization of H7N7 avian influenza viruses in two epidemiologically linked outbreaks. Transbound. Emerg. Dis. 2018, 65, 1576–1587. [Google Scholar] [CrossRef]
  47. Caron, A.; Abolnik, C.; Mundava, J.; Gaidet, N.; Burger, C.E.; Mochotlhoane, B.; Bruinzeel, L.; Chiweshe, N.; de Garine-Wichatitsky, M.; Cumming, G.S. Persistence of low pathogenic avian influenza virus in waterfowl in a Southern African ecosystem. EcoHealth 2011, 8, 109–115. [Google Scholar] [CrossRef]
  48. Van Dijk, J.G.; Hoye, B.J.; Verhagen, J.H.; Nolet, B.A.; Fouchier, R.A.; Klaassen, M. Juveniles and migrants as drivers for seasonal epizootics of avian influenza virus. J. Anim. Ecol. 2014, 83, 266–275. [Google Scholar] [CrossRef]
  49. Kouba, M.; Bartoš, L.; Tomášek, V.; Popelková, A.; Šťastný, K.; Zárybnická, M. Home range size of Tengmalm’s owl during breeding in Central Europe is determined by prey abundance. PLoS ONE 2017, 12, e0177314. [Google Scholar] [CrossRef]
  50. Tucker, M.A.; Alexandrou, O.; Bierregaard, R.O., Jr.; Bildstein, K.L.; Böhning-Gaese, K.; Bracis, C.; Brzorad, J.N.; Buechley, E.R.; Cabot, D.; Calabrese, J.M. Large birds travel farther in homogeneous environments. Glob. Ecol. Biogeogr. 2019, 28, 576–587. [Google Scholar] [CrossRef]
  51. Canavan, B.C. Opening Pandora’s Box at the roof of the world: Landscape, climate and avian influenza (H5N1). Acta Trop. 2019, 196, 93–101. [Google Scholar] [CrossRef]
  52. Cai, W.; McPhaden, M.J.; Grimm, A.M.; Rodrigues, R.R.; Taschetto, A.S.; Garreaud, R.D.; Dewitte, B.; Poveda, G.; Ham, Y.-G.; Santoso, A. Climate impacts of the El Niño–southern oscillation on South America. Nat. Rev. Earth Environ. 2020, 1, 215–231. [Google Scholar] [CrossRef]
  53. Di Pillo, F.; Baumberger, C.; Salazar, C.; Galdames, P.; Ruiz, S.; Sharp, B.; Freiden, P.; Tan, S.; Schultz-Cherry, S.; Hamilton-West, C. Novel Low Pathogenic Avian Influenza H6N1 in Backyard Chicken in Easter Island (Rapa Nui), Chilean Polynesia. Viruses 2022, 14, 718. [Google Scholar] [CrossRef] [PubMed]
  54. Mathieu, C.; Moreno, V.; Pedersen, J.; Jeria, J.; Agredo, M.; Gutiérrez, C.; García, A.; Vásquez, M.; Avalos, P.; Retamal, P. Avian influenza in wild birds from Chile, 2007–2009. Virus Res. 2015, 199, 42–45. [Google Scholar] [CrossRef] [PubMed]
  55. Ghersi, B.M.; Blazes, D.L.; Icochea, E.; Gonzalez, R.I.; Kochel, T.; Tinoco, Y.; Sovero, M.M.; Lindstrom, S.; Shu, B.; Klimov, A. Avian influenza in wild birds, central coast of Peru. Emerg. Infect. Dis. 2009, 15, 935. [Google Scholar] [CrossRef]
  56. Farnsworth, M.L.; Miller, R.S.; Pedersen, K.; Lutman, M.W.; Swafford, S.R.; Riggs, P.D.; Webb, C.T. Environmental and demographic determinants of avian influenza viruses in waterfowl across the contiguous United States. PLoS ONE 2012, 7, e32729. [Google Scholar]
  57. Fuller, T.; Ducatez, M.; Njabo, K.; Couacy-Hymann, E.; Chasar, T.; Aplogan, G.; Lao, S.; Awoume, F.; Tehou, A.; Langeois, Q. Seasonality dynamics of avian influenza occurrences in Central and West Africa. bioRxiv 2014. [Google Scholar] [CrossRef]
  58. Ashok, A.; Rani, H.P.; Jayakumar, K. Monitoring of dynamic wetland changes using NDVI and NDWI based landsat imagery. Remote Sens. Appl. Soc. Environ. 2021, 23, 100547. [Google Scholar] [CrossRef]
  59. Wu, X.; Lv, M.; Jin, Z.; Michishita, R.; Chen, J.; Tian, H.; Tu, X.; Zhao, H.; Niu, Z.; Chen, X. Normalized difference vegetation index dynamic and spatiotemporal distribution of migratory birds in the Poyang Lake wetland, China. Ecol. Indic. 2014, 47, 219–230. [Google Scholar] [CrossRef]
  60. Gaidet, N. Ecology of avian influenza virus in wild birds in tropical Africa. Avian Dis. 2016, 60, 296–301. [Google Scholar] [CrossRef]
  61. Mathieu, C.; Moreno, V.; Retamal, P.; Gonzalez, A.; Rivera, A.; Fuller, J.; Jara, C.; Lecocq, C.; Rojas, M.; García, A. Pandemic (H1N1) 2009 in breeding turkeys, Valparaiso, Chile. Emerg. Infect. Dis. 2010, 16, 709. [Google Scholar] [CrossRef] [PubMed]
  62. Pereda, A.J.; Uhart, M.; Perez, A.A.; Zaccagnini, M.E.; La Sala, L.; Decarre, J.; Goijman, A.; Solari, L.; Suarez, R.; Craig, M.I. Avian influenza virus isolated in wild waterfowl in Argentina: Evidence of a potentially unique phylogenetic lineage in South America. Virology 2008, 378, 363–370. [Google Scholar] [CrossRef] [PubMed]
  63. Sullivan, B.L.; Wood, C.L.; Iliff, M.J.; Bonney, R.E.; Fink, D.; Kelling, S. eBird: A citizen-based bird observation network in the biological sciences. Biol. Conserv. 2009, 142, 2282–2292. [Google Scholar] [CrossRef]
Viruses 15 01241 g001 550
Figure 1. Prevalence (95% CI) of influenza A virus in the Lluta wetland during the sampling months.
Figure 1. Prevalence (95% CI) of influenza A virus in the Lluta wetland during the sampling months.
Viruses 15 01241 g001
Viruses 15 01241 g002 550
Figure 2. Total number of birds counted (A) and species richness (B) identified during the sampling occasions (month–year) in the Lluta River wetland.
Figure 2. Total number of birds counted (A) and species richness (B) identified during the sampling occasions (month–year) in the Lluta River wetland.
Viruses 15 01241 g002
Viruses 15 01241 g003 550
Figure 3. NDVI images in Lluta wetland in the months of highest (March 2018) and lowest (October 2016) vegetation cover values during the study period.
Figure 3. NDVI images in Lluta wetland in the months of highest (March 2018) and lowest (October 2016) vegetation cover values during the study period.
Viruses 15 01241 g003
Viruses 15 01241 g004 550
Figure 4. Maximum likelihood phylogenetic analysis of the hemagglutinins H5, H7, and H9. North American clade in blue, Eurasian clade in purple, and South American clade in green. Study sequences in red and italic. Bootstrap support > 70% indicated.
Figure 4. Maximum likelihood phylogenetic analysis of the hemagglutinins H5, H7, and H9. North American clade in blue, Eurasian clade in purple, and South American clade in green. Study sequences in red and italic. Bootstrap support > 70% indicated.
Viruses 15 01241 g004
Table
Table 1. Definition of ecological and environmental variables.
Table 1. Definition of ecological and environmental variables.
VariablesDefinition
Wild bird communityTotal abundanceTotal number of birds present at the time of sampling
Species richnessNumber of species present at the time of sampling
Abundance of migrantsNumber of migratory birds present at the time of sampling
LandscapeVegetation coverageMean NDVI for the month of sampling, one month before, two months before, and three months before sampling.
Water body size (km2)Water body size for the month of sampling, one month before, two months before, and three months before sampling.
Meteorological dataMaximum monthly temperature (°C)Monthly mean of maximum daily temperature for the month of sampling, one month before, two months before, and three months before sampling.
Minimum monthly temperature (°C)Monthly mean of minimum daily temperature for the month of sampling, one month before, two months before, and three months before sampling.
Total monthly rainfall (mm)Total rainfall at the month of sampling, one month before, two months before, and three months before sampling.
Humidity (%)Relative air humidity at the month of sampling, one month before, two months before, and three months before sampling.
Table
Table 2. Multivariable model results of the association between IAV prevalence and environmental variables.
Table 2. Multivariable model results of the association between IAV prevalence and environmental variables.
VariablesCategoriesEstimatep-ValueOR95% CI
(Intercept) −5.5872<0.0010.003(0.001–0.01)
NDVILow (<0.27)Reference
High (≥0.27)1.29490.002303.65(1.58–8.39)
Abundance of MigrantsLow (<113)Reference
High (≥113)1.27270.007333.57(1.41–9.05)
IAVs isolated from wild birds in and hosts species.
Table
Table 3. IAV strain name by year, subtypes, host species, and GenBank accession of IAVs identified in the Lluta River wetland between September 2015 and October 2020.
Table 3. IAV strain name by year, subtypes, host species, and GenBank accession of IAVs identified in the Lluta River wetland between September 2015 and October 2020.
Strain NameYearIsolate SubtypeHost SpeciesOrderAccession Number
A/American oystercatcher/Chile/C1307/20152015H9N2Haematopus palliatusCharadriiformesKX185901
A/Grey plover/Chile/C1313/20152015H9N7Pluvialis squatarolaCharadriiformesKX185918
A/American oystercatcher/Chile/C17359/20162016H3N8Haematopus palliatusCharadriiformesMH499035
A/American oystercatcher/Chile/C17400/20162016H3N8Haematopus palliatusCharadriiformesMH498968
A/Franklin’s gull/Chile/C17421/20162016H13N9Larus pipixcanCharadriiformesMH498978
A/Franklin’s gull/Chile/C17422/20162016H13N9Larus pipixcanCharadriiformesMH499057
A/Kelp gull/Chile/C27733/20172017H13N8Larus dominicanusCharadriiformesMH499142
A/Kelp gull/Chile/C35523/20172017H5N3Larus dominicanusCharadriiformesMK164009
A/blackish oystercatcher/Chile/C40194/20182018H13N2Haematopus aterCharadriiformesOP888556
A/Chilean flamingo/Chile/C52796/20202020H7N9Phoenicopterus chilensisPhoenicopteriformesOQ820949
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

Ruiz, S.; Galdames, P.; Baumberger, C.; Gonzalez, M.A.; Rojas, C.; Oyarzun, C.; Orozco, K.; Mattar, C.; Freiden, P.; Sharp, B.; et al. Remote Sensing and Ecological Variables Related to Influenza A Prevalence and Subtype Diversity in Wild Birds in the Lluta Wetland of Northern Chile. Viruses 2023, 15, 1241. https://doi.org/10.3390/v15061241

AMA Style

Ruiz S, Galdames P, Baumberger C, Gonzalez MA, Rojas C, Oyarzun C, Orozco K, Mattar C, Freiden P, Sharp B, et al. Remote Sensing and Ecological Variables Related to Influenza A Prevalence and Subtype Diversity in Wild Birds in the Lluta Wetland of Northern Chile. Viruses. 2023; 15(6):1241. https://doi.org/10.3390/v15061241

Chicago/Turabian Style

Ruiz, Soledad, Pablo Galdames, Cecilia Baumberger, Maria Antonieta Gonzalez, Camila Rojas, Cristobal Oyarzun, Katherinne Orozco, Cristian Mattar, Pamela Freiden, Bridgette Sharp, and et al. 2023. "Remote Sensing and Ecological Variables Related to Influenza A Prevalence and Subtype Diversity in Wild Birds in the Lluta Wetland of Northern Chile" Viruses 15, no. 6: 1241. https://doi.org/10.3390/v15061241

APA Style

Ruiz, S., Galdames, P., Baumberger, C., Gonzalez, M. A., Rojas, C., Oyarzun, C., Orozco, K., Mattar, C., Freiden, P., Sharp, B., Schultz-Cherry, S., Hamilton-West, C., & Jimenez-Bluhm, P. (2023). Remote Sensing and Ecological Variables Related to Influenza A Prevalence and Subtype Diversity in Wild Birds in the Lluta Wetland of Northern Chile. Viruses, 15(6), 1241. https://doi.org/10.3390/v15061241

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