Contrasting effects of host species and phylogenetic diversity on the occurrence of HPAI H5N1 in European wild birds

Strong and consistent effects of waterbird composition on HPAI H5 occurrences across Europe

Since 2014, highly pathogenic avian influenza (HPAI) H5 viruses of clade 2.3.4.4 have been dominating the outbreaks across Europe, causing massive deaths among poultry and wild birds. However, the factors shaping these broad-scale outbreak patterns, especially those related to waterbird community composition, remain unclear. In particular, we do not know whether these risk factors differ from those of other H5 clades. Addressing this knowledge gap is important for predicting and preventing future HPAI outbreaks. Using extensive waterbird survey datasets from about 6883 sites, we here explored the effect of waterbird community composition on HPAI H5Nx (clade 2.3.4.4) spatial patterns in the 2016/2017 and 2020/2021 epidemics in Europe, and compared it with the 2005/2006 HPAI H5N1 (clade 2.2) epidemic. We showed that HPAI H5 occurrences in wild birds in the three epidemics were strongly associated with very similar waterbird community attributes, which suggested that, in nature, similar interspecific transmission processes operate between the HPAI H5 subtypes or clades. Importantly, community phylogenetic diversity consistently showed a negative association with H5 occurrence in all three epidemics, suggesting a dilution effect of phylogenetic diversity. In contrast, waterbird community variables showed much weaker associations with HPAI H5Nx occurrence in poultry. Our results demonstrate that models based on previous epidemics can predict future HPAI H5 patterns in wild birds, implying that it is important to include waterbird community factors in future HPAI studies to predict outbreaks and improve surveillance activities.

Spatio-temporal dynamics and drivers of highly pathogenic avian influenza H5N1 in Chile

Introduction

Highly pathogenic avian influenza A H5N1 clade 2.3.4.4b (hereafter H5N1) is causing vast impacts on biodiversity and poultry around the globe. In Chile, lethal H5N1 cases have been reported in a wide range of wild bird species, marine mammals, backyard and industrial poultry, and humans. This study describes the spatio-temporal patterns of the current epizootic of H5N1 in Chile and test drivers that could be associated with outbreak occurrence.

Methods

We used H5N1 cases reported by the Chilean National Animal Health Authority from 5 December 2022 to 5 April 2023. These included wild bird cases confirmed through an avian influenza-specific real-time reverse transcription PCR assay (RT-qPCR), obtained from passive and active surveillance. Data were analyzed to detect the presence of H5N1 clusters under space-time permutation probability modeling, the association of H5N1 with distance and days since the first outbreak through linear regression, and the correlation of H5N1 presence with a number of ecological and anthropogenic variables using general linear modeling.

Results

From 445 H5N1 identified outbreaks involving 613 individual cases in wild birds, a consistent wave-like spread of H5N1 from north to south was identified, which may help predict hotspots of outbreak risk. For instance, seven statistically significant clusters were identified in central and northern Chile, where poultry production and wildlife mortality are concentrated. The presence of outbreaks was correlated with landscape-scale variables, notably temperature range, bird richness, and human footprint.

Discussion

In less than a year, H5N1 has been associated with the unusual mortality of >100,000 individuals of wild animals in Chile, mainly coastal birds and marine mammals. It is urgent that scientists, the poultry sector, local communities, and national health authorities co-design and implement science-based measures from a One Health perspective to avoid further H5N1 spillover from wildlife to domestic animals and humans, including rapid removal and proper disposal of wild dead animals and the closure of public areas (e.g., beaches) reporting high wildlife mortalities.

The more, the healthier: Tree diversity reduces forest pests and pathogens

How frequently, and under what conditions, biodiversity reduces disease through "dilution effects" has been a subject of ongoing research. A new study of forest pests in PLOS Biology provides strong evidence for their generality.

Effect of epidemic diseases on wild animal conservation

Under the background of global species extinction, the impact of epidemic diseases on wild animal protection is increasingly prominent. Here, we review and synthesize the literature on this topic, and discuss the relationship between diseases and biodiversity. Diseases usually reduce species diversity by decreasing or extinction of species populations, but also accelerate species evolution and promote species diversity. At the same time, species diversity can regulate disease outbreaks through dilution or amplification effects. The synergistic effect of human activities and global change is emphasized, which further aggravates the complex relationship between biodiversity and diseases. Finally, we emphasize the importance of active surveillance of wild animal diseases, which can protect wild animals from potential diseases, maintain population size and genetic variation, and reduce the damage of diseases to the balance of the whole ecosystem and human health. Therefore, we suggest that a background survey of wild animal populations and their pathogens should be carried out to assess the impact of potential outbreaks on the population or species level. The mechanism of dilution and amplification effect between species diversity and diseases of wild animals should be further studied to provide a theoretical basis and technical support for human intervention measures to change biodiversity. Most importantly, we should closely combine the protection of wild animals with the establishment of an active surveillance, prevention, and control system for wild animal epidemics, in an effort to achieve a win-win situation between wild animal protection and disease control.

The episodic resurgence of highly pathogenic avian influenza H5 virus

Highly pathogenic avian influenza (HPAI) H5N1 activity has intensified globally since 2021, increasingly causing mass mortality in wild birds and poultry and incidental infections in mammals1-3. However, the ecological and virological properties that underscore future mitigation strategies still remain unclear. Using epidemiological, spatial and genomic approaches, we demonstrate changes in the origins of resurgent HPAI H5 and reveal significant shifts in virus ecology and evolution. Outbreak data show key resurgent events in 2016-2017 and 2020-2021, contributing to the emergence and panzootic spread of H5N1 in 2021-2022. Genomic analysis reveals that the 2016-2017 epizootics originated in Asia, where HPAI H5 reservoirs are endemic. In 2020-2021, 2.3.4.4b H5N8 viruses emerged in African poultry, featuring mutations altering HA structure and receptor binding. In 2021-2022, a new H5N1 virus evolved through reassortment in wild birds in Europe, undergoing further reassortment with low-pathogenic avian influenza in wild and domestic birds during global dissemination. These results highlight a shift in the HPAI H5 epicentre beyond Asia and indicate that increasing persistence of HPAI H5 in wild birds is facilitating geographic and host range expansion, accelerating dispersion velocity and increasing reassortment potential. As earlier outbreaks of H5N1 and H5N8 were caused by more stable genomic constellations, these recent changes reflect adaptation across the domestic-bird-wild-bird interface. Elimination strategies in domestic birds therefore remain a high priority to limit future epizootics.

Functional traits explain waterbirds' host status, subtype richness, and community-level infection risk for avian influenza

Species functional traits can influence pathogen transmission processes, and consequently affect species' host status, pathogen diversity, and community-level infection risk. We here investigated, for 143 European waterbird species, effects of functional traits on host status and pathogen diversity (subtype richness) for avian influenza virus at species level. We then explored the association between functional diversity and HPAI H5Nx occurrence at the community level for 2016/17 and 2021/22 epidemics in Europe. We found that both host status and subtype richness were shaped by several traits, such as diet guild and dispersal ability, and that the community-weighted means of these traits were also correlated with community-level risk of H5Nx occurrence. Moreover, functional divergence was negatively associated with H5Nx occurrence, indicating that functional diversity can reduce infection risk. Our findings highlight the value of integrating trait-based ecology into the framework of diversity-disease relationship, and provide new insights for HPAI prediction and prevention.

Recurring Trans-Atlantic Incursion of Clade 2.3.4.4b H5N1 Viruses by Long Distance Migratory Birds from Northern Europe to Canada in 2022/2023

In December 2022 and January 2023, we isolated clade 2.3.4.4b H5N1 high-pathogenicity avian influenza (HPAI) viruses from six American crows (Corvus brachyrhynchos) from Prince Edward Island and a red fox (Vulpes vulpes) from Newfoundland, Canada. Using full-genome sequencing and phylogenetic analysis, these viruses were found to fall into two distinct phylogenetic clusters: one group containing H5N1 viruses that had been circulating in North and South America since late 2021, and the other one containing European H5N1 viruses reported in late 2022. The transatlantic re-introduction for the second time by pelagic/Icelandic bird migration via the same route used during the 2021 incursion of Eurasian origin H5N1 viruses into North America demonstrates that migratory birds continue to be the driving force for transcontinental dissemination of the virus. This new detection further demonstrates the continual long-term threat of H5N1 viruses for poultry and mammals and the subsequent impact on various wild bird populations wherever these viruses emerge. The continual emergence of clade 2.3.4.4b H5Nx viruses requires vigilant surveillance in wild birds, particularly in areas of the Americas, which lie within the migratory corridors for long-distance migratory birds originating from Europe and Asia. Although H5Nx viruses have been detected at higher rates in North America since 2021, a bidirectional flow of H5Nx genes of American origin viruses to Europe has never been reported. In the future, coordinated and systematic surveillance programs for HPAI viruses need to be launched between European and North American agencies.

Risk factors for Lyme disease: A scale-dependent effect of host species diversity and a consistent negative effect of host phylogenetic diversity

Biodiversity can influence disease risk. One example of a diversity-disease relationship is the dilution effect, which suggests higher host species diversity (often indexed by species richness) reduces disease risk. While numerous studies support the dilution effect, its generality remains controversial. Most studies of diversity-disease relationships have overlooked the potential importance of phylogenetic diversity. Furthermore, most studies have tested diversity-disease relationships at one spatial scale, even though such relationships are likely scale dependent. Using Lyme disease as a model system, we investigated the effects of host species richness and phylogenetic relatedness on the number of reported Lyme disease cases in humans in the U.S.A. at two spatial scales (the county level and the state level) using piecewise structural equation modelling. We also accounted for relevant climatic and habitat-related factors and tested their correlations with the number of Lyme disease cases. We found that species assemblages with more related species (i.e., host species in the order Rodentia) were associated with more Lyme disease cases in humans. Host species richness correlated negatively with the number of Lyme disease cases at the state level (i.e., a dilution effect), a pattern that might be explained by the higher number of reservoir-incompetent species at high levels of species richness at this larger spatial scale. In contrast, a positive correlation was found between species richness and the number of Lyme disease cases at the county level, where a higher proportion of rodent species was associated with higher levels of species richness, potentially amplifying the disease risk. Our results highlight that analyse at a single spatial scale can miss some impacts of biodiversity on human health. Thus, multi-scale analyses with consideration of host phylogenetic diversity are critical for improving our understanding of diversity-disease relationships.

Wild Bird Densities and Landscape Variables Predict Spatial Patterns in HPAI Outbreak Risk across The Netherlands

Highly pathogenic avian influenza viruses’ (HPAIVs) transmission from wild birds to poultry occurs globally, threatening animal and public health. To predict the HPAI outbreak risk in relation to wild bird densities and land cover variables, we performed a case-control study of 26 HPAI outbreaks (cases) on Dutch poultry farms, each matched with four comparable controls. We trained machine learning classifiers to predict outbreak risk with predictors analyzed at different spatial scales. Of the 20 best explaining predictors, 17 consisted of densities of water-associated bird species, 2 of birds of prey, and 1 represented the surrounding landscape, i.e., agricultural cover. The spatial distribution of mallard (Anas platyrhynchos) contributed most to risk prediction, followed by mute swan (Cygnus olor), common kestrel (Falco tinnunculus) and brant goose (Branta bernicla). The model successfully distinguished cases from controls, with an area under the receiver operating characteristic curve of 0.92, indicating accurate prediction of HPAI outbreak risk despite the limited numbers of cases. Different classification algorithms led to similar predictions, demonstrating robustness of the risk maps. These analyses and risk maps facilitate insights into the role of wild bird species and support prioritization of areas for surveillance, biosecurity measures and establishments of new poultry farms to reduce HPAI outbreak risks.

Host diversity and behavior determine patterns of interspecies transmission and geographic diffusion of avian influenza A subtypes among North American wild reservoir species

Wild birds can carry avian influenza viruses (AIV), including those with pandemic or panzootic potential, long distances. Even though AIV has a broad host range, few studies account for host diversity when estimating AIV spread. We analyzed AIV genomic sequences from North American wild birds, including 303 newly sequenced isolates, to estimate interspecies and geographic viral transition patterns among multiple co-circulating subtypes. Our results show high transition rates within Anseriformes and Charadriiformes, but limited transitions between these orders. Patterns of transition between species were positively associated with breeding habitat range overlap, and negatively associated with host genetic distance. Distance between regions (negative correlation) and summer temperature at origin (positive correlation) were strong predictors of transition between locations. Taken together, this study demonstrates that host diversity and ecology can determine evolutionary processes that underlie AIV natural history and spread. Understanding these processes can provide important insights for effective control of AIV.

A high-quality genome and comparison of short- versus long-read transcriptome of the palaearctic duck Aythya fuligula (tufted duck)

BACKGROUND

The tufted duck is a non-model organism that experiences high mortality in highly pathogenic avian influenza outbreaks. It belongs to the same bird family (Anatidae) as the mallard, one of the best-studied natural hosts of low-pathogenic avian influenza viruses. Studies in non-model bird species are crucial to disentangle the role of the host response in avian influenza virus infection in the natural reservoir. Such endeavour requires a high-quality genome assembly and transcriptome.

FINDINGS

This study presents the first high-quality, chromosome-level reference genome assembly of the tufted duck using the Vertebrate Genomes Project pipeline. We sequenced RNA (complementary DNA) from brain, ileum, lung, ovary, spleen, and testis using Illumina short-read and Pacific Biosciences long-read sequencing platforms, which were used for annotation. We found 34 autosomes plus Z and W sex chromosomes in the curated genome assembly, with 99.6% of the sequence assigned to chromosomes. Functional annotation revealed 14,099 protein-coding genes that generate 111,934 transcripts, which implies a mean of 7.9 isoforms per gene. We also identified 246 small RNA families.

CONCLUSIONS

This annotated genome contributes to continuing research into the host response in avian influenza virus infections in a natural reservoir. Our findings from a comparison between short-read and long-read reference transcriptomics contribute to a deeper understanding of these competing options. In this study, both technologies complemented each other. We expect this annotation to be a foundation for further comparative and evolutionary genomic studies, including many waterfowl relatives with differing susceptibilities to avian influenza viruses.

Spatiotemporal risk assessment for avian influenza outbreak based on the dynamics of habitat suitability for wild birds

Highly pathogenic avian influenza (HPAI) has predominantly damaged the poultry industry worldwide. The fundamental prevention and control strategy for HPAI includes early detection and timely intervention enforcement through a systematic surveillance system for wild birds based on the ecological understanding of the dynamics of wild birds' movements. Our study aimed to develop a spatiotemporal risk assessment model for avian influenza (AI) infection in wild birds to empower surveillance information for a contingency strategy. For this purpose, first, we predicted the monthly habitat suitability of seven waterfowl species, using 227,671 Global Positioning System (GPS) tracking records of 562 birds from 2014 to 2018 in the Republic of Korea (ROK). Then, that predicted habitat suitability and 421 coordinates of AI detection sites in wild birds were used to build the risk assessment model. Subsequently, we compared the monthly predicted risk of avian influenza virus (AIv) identification in wild birds between case and non-case poultry farms with HPAI H5N6 outbreak in the ROK between 2016 and 2017. The results reported considerable variation of monthly habitat suitability of seven waterfowls and the impact of predicting AI occurrences in wild birds. The high habitat suitability for spot-billed ducks (contribution rate in November = 40.9%) and mallards (contribution rate in January = 34.3%) significantly contributed to predicting the average risk of AIv identification in wild birds, with high predictive performance [the monthly mean of area under the curve (AUC) = 0.978]. Moreover, our model showed that the averaged risk of identification AI in wild birds was significantly higher in HPAI infected premises, with infected domestic duck holdings exhibiting a significantly higher risk than the chicken farms in November. This study suggests that animal health authority establishes a risk-based HPAI surveillance system grounded on the ecological nature of wild birds to improve the effectiveness of prevention and preparedness of emerging epidemics.

Mammal assemblage composition predicts global patterns in emerging infectious disease risk

As a source of emerging infectious diseases, wildlife assemblages (and related spatial patterns) must be quantitatively assessed to help identify high-risk locations. Previous assessments have largely focussed on the distributions of individual species; however, transmission dynamics are expected to depend on assemblage composition. Moreover, disease-diversity relationships have mainly been studied in the context of species loss, but assemblage composition and disease risk (e.g. infection prevalence in wildlife assemblages) can change without extinction. Based on the predicted distributions and abundances of 4466 mammal species, we estimated global patterns of disease risk through the calculation of the community-level basic reproductive ratio R0, an index of invasion potential, persistence, and maximum prevalence of a pathogen in a wildlife assemblage. For density-dependent diseases, we found that, in addition to tropical areas which are commonly viewed as infectious disease hotspots, northern temperate latitudes included high-risk areas. We also forecasted the effects of climate change and habitat loss from 2015 to 2035. Over this period, many local assemblages showed no net loss of species richness, but the assemblage composition (i.e. the mix of species and their abundances) changed considerably. Simultaneously, most areas experienced a decreased risk of density-dependent diseases but an increased risk of frequency-dependent diseases. We further explored the factors driving these changes in disease risk. Our results suggest that biodiversity and changes therein jointly influence disease risk. Understanding these changes and their drivers and ultimately identifying emerging infectious disease hotspots can help health officials prioritize resource distribution.

Effects of host extinction and vector preferences on vector-borne disease risk in phylogenetically structured host-hector communities

Anthropogenic disturbance impacts the phylogenetic composition and diversity of ecological communities. While changes in diversity are known to dramatically change species interactions and alter disease dynamics, the effects of phylogenetic changes in host and vector communities on disease have been relatively poorly studied. Using a theoretical model, we investigated how phylogeny and extinction influence network structural characteristics relevant to disease transmission in disturbed environments. We modelled a multi-host, multi-vector community as a bipartite ecological network, where nodes represent host and vector species and edges represent connections among them through vector feeding, and we simulated vector preferences and threat status on host and parasite phylogenies. We then simulated loss of hosts, including phylogenetically clustered losses, to investigate how extinction influences network structure. We compared effects of phylogeny and extinction to those of host specificity, which we predicted to strongly increase network modularity and reduce disease prevalence. The simulations revealed that extinction often increased modularity, with higher modularity as species loss increased, although not as much as increasing host specificity did. These results suggest that extinction itself, all else being equal, may reduce disease prevalence in disturbed communities. However, in real communities, systematic patterns in species loss (e.g. favoring high competence species) or changes in abundance may counteract these effects. Unexpectedly, we found that effects of phylogenetic signal in host and vector traits were relatively weak, and only important when phylogenetic signal of host and vector traits were similar, or when these traits both varied.

Effects of migration network configuration and migration synchrony on infection prevalence in geese

Migration can influence dynamics of pathogen-host interactions. However, it is not clearly known how migration pattern, in terms of the configuration of the migration network and the synchrony of migration, affects infection prevalence. We therefore applied a discrete-time SIR model, integrating environmental transmission and migration, to various migration networks, including networks with serial, parallel, or both serial and parallel stopover sites, and with various levels of migration synchrony. We applied the model to the infection of avian influenza virus in a migratory geese population. In a network with only serial stopover sites, increasing the number of stopover sites reduced infection prevalence, because with every new stopover site, the amount of virus in the environment was lower than that in the previous stopover site, thereby reducing the exposure of the migratory population. In a network with parallel stopover sites, both increasing the number and earlier appearance of the stopover sites led to an earlier peak of infection prevalence in the migratory population, because the migratory population is exposed to larger total amount of virus in the environment, speeding-up the infection accumulation. Furthermore, higher migration synchrony reduced the average number of cumulative infection, because the majority of the population can fly to a new stopover site where the amount of virus is still relatively low and has not been increased due to virus shedding of infected birds. Our simulations indicate that a migration pattern with multiple serial stopover sites and with highly synchronized migration reduces the infection prevalence.