Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Nature Communications volume 14, Article number: 5489 (2023)
4711
120
Metrics details
Highly pathogenic avian influenza (HPAI) A/H5N1 viruses (lineage 2.3.4.4b) are rapidly invading the Americas, threatening wildlife, poultry, and potentially evolving into the next global pandemic. In November 2022 HPAI arrived in Peru, triggering massive pelican and sea lion die-offs. We report genomic characterization of HPAI/H5N1 in five species of marine mammals and seabirds (dolphins, sea lions, sanderlings, pelicans and cormorants). Peruvian viruses belong to lineage 2.3.4.4b, but they are 4:4 reassortants where 4 genomic segments (PA, HA, NA and MP) position within the Eurasian lineage that initially entered North America from Eurasia, while the other 4 genomic segments (PB2, PB1, NP and NS) position within the American lineage (clade C) that circulated in North America. These viruses are rapidly accruing mutations, including mutations of concern, that warrant further examination and highlight an urgent need for active local surveillance to manage outbreaks and limit spillover into other species, including humans.
The recent emergence of highly pathogenic avian influenza (HPAI) H5N1 viruses in mammals and birds in the Americas1,2 presents a severe threat to wild and endangered species, to poultry production3,4, and to public health when the virus spills over into humans3,5. The H5N1 (clade 2.3.4.4b) virus arrived in North America in late 2021 from Eurasia and spread across the continent in wild birds, spilling over into poultry farms and infecting an alarming number of wild terrestrial mammals, including fox, skunk, bear, bobcat, and raccoon6,7,8. In October 2022, an apparent mink-to-mink transmission of clade 2.3.4.4b H5N1 viruses in Spain9 further heightened concern that the avian virus was adapting to mammals and that an H5N1 global pandemic in humans could be approaching.
In November 2022, Peruvian pelicans (Pelecanus thagus) along the coast and in offshore islands began experiencing a mass die-off10,11. HPAI A/H5N1 was identified in samples collected from dead birds, resulting in the country declaring a sanitary alert10. This was followed by several spillover events into other domestic and wild birds, including zoo animals and wild raptors4,11. By the beginning of 2023, the Peru outbreak had spread to marine mammals, particularly affecting the South American sea lion (Otaria flavescens), which also began to experience a mass die-off12.
The Pacific coast of Peru hosts a rich biodiversity of marine mammals and seabirds13. The ecosystem is home to large populations of South American sea lion, “guano birds” like the Guanay cormorant (Phalacrocorax bougainvillii), the Peruvian booby (Sula variegata) and the Peruvian pelican, and endangered birds like the Humboldt penguin (Spheniscus humboldti). The region also serves as stopover points and feeding grounds for diverse avian migratory species, including the Franklin’s gull (Leucophaeus pipixcan) and several species of sandpipers (Calidris spp.)14,15 that travel south from the northern hemisphere during the boreal winter. Endemic avian influenza A viruses circulate in birds in Peru16, including unusual reassortant viruses with combinations of genes from the dominant “American” lineage that is widespread in North America and the “South American” lineage that is commonly found in waterfowl in Argentina and Chile17,18. However, HPAI viruses of Eurasian origin that caused a major US outbreak in 2014–201519 never reached Peru and, to our knowledge, the 2022–2023 HPAI outbreak in South America represents the first incursion of HPAI in the region.
Here we report detection, genomic characterisation, phylogenetic analysis, and mutation analysis of HPAI A/H5N1 viruses identified in marine mammals (sea lion and common dolphin) and seabirds (pelican, cormorant and sanderling) sampled along the coast of Peru since November 2022.
In response to a multi-species outbreak starting in November 2022, we collected 69 swabs of external orifices and internal organ tissues from seven species of marine mammals and seabirds (common dolphin (Delphinus delphis), South American sea lion (Otaria flavescens), sanderling (Calidris alba), Peruvian pelican (Pelecanus thagus), Guanay cormorant (Phalacrocorax bougainvillii), Peruvian gull (Larus belcheri) and Humboldt penguin (Spheniscus humboldti)) in regions representing the northern (Piura), central (Lima) and southern (Arequipa and Tacna) coast of Peru (Supplementary Table 1, Fig. 1A). All animals sampled were either deceased or manifested clear signs of disease, including respiratory, digestive, and/or neurological symptoms indicative of acute encephalitis. Seabirds exhibited disorientation, ataxia, circling, nystagmus, torticollis, congested conjunctivae and dyspnea; sea lions exhibited disorientation, ataxia, circling, “stargazing” posture, copious nasal discharge, sialorrhea, dyspnea and seizures; the dolphin was found recently deceased.
A Map of Peru showing collection sites representative of the northern, central and southern regions. Mass die-offs of both birds and mammals occurred along the entire coast of Peru, but map only shows positive locations of animals samples in this study. B Photographic record of animals sampled, including common dolphin, South American sea lion, sanderling, Guanay cormorant, Peruvian booby and Peruvian pelican.
To date, we have tested 69 samples from a total of 28 individuals by RT-qPCR, confirming 11 individuals (1 dolphin, 4 sea lions and 6 seabirds), plus 1 pooled sample from 5 sea lions, as positive for influenza A. All samples were also tested for coronaviruses, alphaviruses, bunyaviruses and flaviviruses using pan-PCR assays, and all samples were negative for these tests. Influenza A positives were subjected to NGS for subtyping and to generate full genomes for phylogenetic and mutational analysis. We confirmed 11 individuals as positive for HPAI A/H5N1 (1/1 dolphin, 3/4 individual sea lions, 1/1 pooled sample from 5 sea lions and 6/6 seabirds). It was not possible to type the remaining positive sea lion, as the sample contained very low viral loads or extensive signs of nucleic acid degradation (Supplementary Table 1), likely associated with extensively decomposed tissues in deceased animals (Fig. 1B). Of the 11 individuals that yielded quality sequence data, we generated complete sequences for most genomic regions in most samples (Table 1 and Supplementary Table 1). All sequences have been deposited in GenBank (Accession Numbers OQ550419-OQ550478, and OQ925704-OQ925729).
To further characterise our Peruvian HPAI A/H5N1 isolates we carried out phylogenetic analyses for each genome segment separately, using all HPAI A/H5N1 reference sequences available globally for avian and mammalian H5 viruses submitted to GISAID since January 1, 2021, including the original A/H5N1 goose/Guangdong strain identified in 199620. Our analysis confirms that all samples sequenced are closely related and cluster within the HPAI A/H5N1 lineage 2.3.4.4b, and further, that the eight Peruvian HPAI A/H5N1 viruses are reassortants (Figs. 2 and 3). Peruvian viruses are positioned in the Eurasian lineage in trees inferred for the PA, HA, NA and MP segments, and cluster within the large clade of North American viruses that were first introduced from Eurasia to the North American’s Atlantic flyway in late 2021 (Fig. 2). In contrast, Peruvian viruses are positioned in the American lineage for trees inferred for the PB2, PB1, NP and NS segments, which is evidence of a 4:4 reassortment event. Most of the H5N1 viruses observed since the summer of 2022 in the Americas are reassortants containing various combinations of segments from the Eurasian and American lineages (Fig. 3). The Eurasian lineage H5N1 virus that originally invaded North America in 2021 has reassorted multiple times with the endemic American lineage since arriving in North America, as evidenced by multiple reassortant clades positioned within the American lineage (Fig. 2). The Peru viruses are positioned in clade C on the PB2, PB1, NP and NS trees and their 4:4 reassortment pattern is referred to as “R6” (Fig. 3).
Phylogenetic trees for AIVs collected globally and submitted to GISAID since January 1, 2021, inferred using ML methods. Branches shaded by AIV lineage; tips shaded by AIV subtype and location. H5N1 viruses from the 2021–2023 Western hemisphere outbreak are shaded yellow, South American H5N1 viruses (except Peru) are shaded green, and Peruvian H5N1 viruses are shaded pink. Cartoon animals denote notable outbreaks in mink in Spain, sea lions in Peru, and seals in Maine, USA. A detailed tree is provided for the PB2 segment, including bootstrap values, labels, and a box around the R6 reassortant clade that is presented in greater detail in Fig. 3A. Smaller trees are provided for the other 7 segments. Raw tree files for each genome segment are available at GitHub (https://github.com/mostmarmot/Peru_AIV/).
Each genotype box represents one of eight genome segments (1–8), shaded by lineage: grey = Eurasian lineage; pink = American lineage (clade A); blue = American lineage (clade B); green = American lineage (clade C); orange = American lineage (clade D).
R6 reassortants were first detected in North American poultry in March 2022 (e.g., EPI_ISL_11971482|A/turkey/South Dakota/22-008485-002/2022) (Fig. 4A), spread regionally across the central United States during the spring and summer of 2022 (Fig. 4B), and became the dominant genotype among sequenced H5N1 viruses in the Americas during the autumn of 2022 (Fig. 4C). The population genetics of H5N1 in the Americas underwent a shift during 2022, as reassortants displaced the original non-reassortant Eurasian H5N1 virus, which has not been detected in the Western hemisphere since June 2022. All currently detected H5N1 reassortants in the Western hemisphere still retain Eurasian HA, NA, and MP segments, but their PB2 and NP segments belong to the American lineage. Both Eurasian and American lineages continue to co-circulate in the Americas for the PB1, PA, and NS segments in various reassortant backgrounds (Fig. 3). Viruses isolated from poultry in Colombia and Ecuador21 in November 2022 (e.g., A/chicken/Ecuador/02/2022) have the same R6 reassortant genotype as the Peru viruses that were also first detected in November 2022. However, the South America H5N1 viruses are not monophyletic (Fig. 4A) and represent at least four independent viral introductions from North American birds migrating to South American countries during fall migration (1 introduction to Ecuador, 1 introduction to Chile/Peru (see below), and 2 introductions to Colombia). The H5N1 viruses isolated from pelicans in Venezuela in November 2022 have a different reassortant genotype (R7) that is primarily found in the eastern US, including in Florida, and represents a fifth independent H5N1introduction from North to South America during the fall of 2022 (Fig. 3).
A Time-scaled maximum clade credibility (MCC) tree inferred for the PB2 segment using Bayesian approaches for 98 representative H5N1 viruses (R6 genotype) belonging to the reassortant clade C in the American lineage (see Fig. 2). Branches shaded by location. South American virus names shaded by country or Peruvian region. Posterior probabilities provided for key nodes. Estimated timing of R6 virus entry into Peru/Chile is provided. Diamonds indicate the 3 viruses with PB2:D701N mutations. Raw MCC tree available on GitHub (https://github.com/mostmarmot/Peru_AIV/). B Discrete phylogeographic transition history of R6 viruses at four time points. Red lines with arrows indicate interhemispheric transitions from North America to South America. C Proportion of H5N1 viruses sequenced in the Americas between 2021–2023 that belong to different genotypes (Fig. 2) over time.
H5N1 viruses from Peru and Chile with the R6 genotype form a single clade, representing the only phylogenetic evidence of H5N1 spread between two South American countries to date (Fig. 4A). The Peru-Chile clade appears to descend from a single viral introduction from North America, estimated to have occurred during October 2022 (95% HPD, September 05, 2022 to October 23, 2022). Whether the introduction from North America first arrived in Peru or Chile is difficult to infer from the tree, since the posterior probabilities are similar for both locations at this node (0.41 for Peru; 0.47 for Chile). Within the Peru-Chile clade, the H5N1 viruses collected from northern Peru (Piura region) and central Peru (Lima region) cluster together, consistent with transmission in Peru between avian and mammalian species (pelican, cormorant, sandpiper, dolphin and sea lion). However, two H5N1 viruses collected from sea lions in southern Peru (Arequipa region) during the latter phases of the outbreak cluster with Chilean viruses, reflecting instances of viral gene flow across Peru’s southern border. Both of these viruses have the PB2 D701N mutation that is associated with enhanced mammalian transmission in mammals (see below), however, these Arequipa viruses were collected in sea lions a month apart (February 7, 2023 vs. March 6, 2023) and they do not cluster together on the same phylogenetic trees. One of the Arequipa viruses (A/South American sea lion/PeruAQP-SER00R/2023) with the D701N mutation clusters on the PB2 tree with the H5N1 virus isolated from a human in Chile22 a few weeks later (March 24, 2023) that also has the D701N mutation (Figs. 4A and 5), but this pattern does not hold for trees inferred using other genome segments. Taken together, these data provide evidence that the D701N mutation is emerging repeatedly in H5N1 viruses that infect mammalian hosts along South America’s Pacific coast, however, these data cannot yet ascertain whether viruses with this mutation are transmitting within mammals.
Additional reference sequences used are available through GenBank and GISAID. A total of 22 variable sites across genomic segments were identified relative to the original A/H5N1 goose/Guangdong reference from 1996, the A/Vietnam/1203/2004 reference used to annotate amino acid positions in the CDC inventory23, plus two additional more recent references representing R6 reassortants from 2022 identified in birds (chicken/wyoming/2022) and mammals (skunk/washington/22-019274-001/2022) (shaded in grey). Mutation PB2 D701N (shown with an *) has been previously linked to mammalian host adaptation and enhanced transmission24,25. The remaining 21/22 sites have not been previously characterised, and 4 may warrant special consideration (shown in bold) as they are concentrated in recent mammalian samples (Supplement Table 2), including 2 (PB2 D701N and PA M86I) that later show up in the genome sequenced from the human case in Chile (A/Chile/25946/2023)22. Dots represent amino acids identical to those present in annotation reference A/Goose/Guangdong/1/96; colours have been randomly assigned to amino acids simply to aid in identification of differences; X indicates there is no genomic information available for that position; – – indicates a deletion in the sequence.
We also performed a detailed SNP and mutational analysis to identify amino acid changes potentially linked to increased virulence, transmission, or mammalian host adaptation, and to assess if we could identify specific differences between host species (mammals vs. birds). We focused on variable sites that were different from both the original A/H5N1 goose/Guangdong strain from 199620 and the A/Vietnam/1203/2004 strain used to annotate amino acid positions in the CDC inventory23. Based on these criteria, we identified more than 70 variable sites spread across all genomic regions (Supplementary Table 2), including 21 that have been previously reported as linked to specific phenotypes, such as altered polymerase activity and replication efficiency (usually enhanced), increased virus binding to α2–3 and α2–6, enhanced transmission, and increased virulence and pathogenicity, including in mammals23,24. However, given the increased age of the two reference sequences used to define variable sites, which would make them unsuitable to identify potentially relevant recent changes, we also included in the analysis two additional references representing R6 reassortants from 2022 identified in birds (A/chicken/Wyoming/22-032216-001/2022) and mammals (A/skunk/Washington/22-019274-001/2022). Based on this more stringent comparison, we identified 23 variable positions in Peruvian viruses that contain mutations that are not present in either the Guangdong/Vietnam reference strains from 1996/2004, or the R6 reassortants from 2022 (Fig. 5, Supplementary Table 2). Within this group, 15 mutations are present in all birds and mammals from Peru (3 cases) or present in single individuals from Peru (12 cases). The remaining 8 sites are of particularly high interest because they repeat two or more times. These sites include 1 in PB2 (D701N), 1 in PB1 (S515A), 3 in PA (R57Q, T85V, M86I), 1 in HA (H355R, according to the H5 numbering scheme), 1 in NP (Y289F) and 1 in NA (A81I), and half of them cluster solely in mammalian viruses from Peru. Particularly noteworthy is the PB2 D701N mutation, which appears in sea lion samples collected during the late phases of the outbreak in Peru in February and March, and subsequently appears in the virus sequenced from the human case reported in Chile22. This mutation, along with two other PB2 mutations (E627K, K702R) not present in our data, have been specifically linked to mammalian host adaptation and enhanced transmission previously24,25. In addition, despite not finding PB2 E627K and K702R, we have identified two additional mutations (PA M861I and NS1 D26K) in mammalian samples from Peru that are also present in the human case from Chile22, raising the possibility that the viruses are changing in a host-specific manner24,25 that supports mammalian host adaptation. The sea lion viruses contain additional unique mutations, such as PB1 L378M and S515A, that may warrant further observation, as these are also showing up in Chilean genomes, mostly from birds. Sampling of additional individuals will be needed to fully assess the significance of these uncharacterised variable sites, especially in terms of their pathogenicity to mammals. However, it is clear that most are restricted to sequences reported in Latin America (Fig. 5, Supplementary Table 2), which indicates that the virus is indeed changing, and perhaps even adapting to mammals, as it travels south from the northern hemisphere.
For decades, South America’s relative geographic isolation has largely insulated its fragile coastal ecosystems and poultry industry from the HPAI outbreaks that periodically ravage US and Mexican farms to the north. But the newest variant of HPAI A/H5N1, clade 2.3.4.4b, is spreading faster, causing mass mortality in wildlife, infecting mammals, and invading countries like Peru that had remained HPAI-free for decades. The arrival of HPAI in regions with less experience managing highly pathogenic viruses in wildlife and poultry is highly concerning. These events should prompt immediate cross-sectoral capacity strengthening and coordinated response activities throughout the region. Here, we rapidly established new surveillance partnerships between government and academia to respond to mass mortality events involving Peruvian pelicans and South American sea lions. We confirmed the presence of HPAI A/H5N1 clade 2.3.4.4b in both pelicans and sea lions, as well as in Guanay cormorants, sanderlings and dolphins, and further surmised that HPAI A/H5N1 is the likely causative agent of the mass wildlife die-offs. We suspect that direct HPAI transmission between sea lions could be occurring, rather than independent spillovers into sea lions from avian sources, but additional sequence data and analysis will be required to further characterise mammal-to-mammal transmission. We report more than 70 variable sites, many of which have been previously linked to altered polymerase activity and replication efficiency, increased virus binding to α2–3 and α2–6, enhanced transmission, and increased virulence and pathogenicity, including in mammals23,24,26, but we specifically focus on 23 sites that are different even from R6 resassortants that circulated in 2022. Within this subset we are particularly concerned by the presence of PB2 D701N in 2 sea lion samples, and in a human case reported in Chile22, as this mutation has been specifically linked to mammalian host adaptation and enhanced transmission24,25. The same concern applies to PA M861I and NS1 D26K, which are present in mammalian samples from Peru and in the human case from Chile22. Two other mutations of concern (PB2 E627K and K702R), linked to mammalian host adaptation and enhanced transmission24,25, as well as the PB2 T271A mutation observed in the mink outbreak in Spain9, were not present in our samples. Our analysis has focused on recent non-synonymous mutations, however, other mutations, including synonymous substitutions should be further examined, as these can change viral RNA structure and splicing, which in turn can potentially impact viral pathogenesis and fitness.
Our phylogenetic analysis supports a single introduction of 2.3.4.4b into Peru from North America, presumably through the movements of migratory wild birds that travel south during the boreal winter, setting the stage for infection of local sea birds that share habitats with marine mammals. There are multiple possible transmission routes for local transmission among species that involve direct contact or indirect environmental transmission. For one, seabirds share feeding spaces with both sea lions and dolphins, providing ample opportunities for direct contact between animals at sea27,28,29,30. Direct contact also occurs on islands, islets, and guano headlands, especially in protected areas where large and dense breeding colonies of sea lions and seabirds cohabitate, and where indirect transmission is also possible on land and via guano runoff into the surrounding waters27,28,30,31. Another scenario for transmission involves carnivory and scavenging of infected animal carcasses by marine and terrestrial carnivores, as well as by raptors, gulls, and other scavenger birds2,32,33. Fishing docks, where fishermen often dispose of waste by dumping it at sea, attract seabirds, sea lions, marine otters and others that come to feed. Many docks along the Peruvian coast also function as tourist attractions, where seabirds and sea lions are purposely fed to create photo opportunities, building large congregations of wild animals that also increase opportunities for contact with humans. The confirmed presence of HPAI A/H5N1 in 2 species of resident guano seabirds, the Peruvian pelican and Guanay cormorant, provides another potential route for future transmission to humans34,35, as guano is widely used to fertilise crops. Finally, the Peruvian desert coast is home to large poultry operations with millions of chickens adjacent to infected wild and migratory birds, placing this production industry at risk and furthering the potential for animal and human contact with circulating HPAI lineages in the region.
There are outstanding questions about which migratory bird species are involved in the long-distance dissemination of HPAI from North to South America, possibly by way of Central America. We detected clade 2.3.4.4b in a migratory sanderling (Calidris alba) that would have arrived in Peru after breeding in the Canadian arctic. However, Calidris spp. are an unlikely conduit for HPAIV A/H5N1 because experimental inoculations result in death or disease within 5 to 11 days of inoculation36,37. Given the unlikelihood of a successful long-distance migration for a clinically infected bird, we suspect the sanderling was infected locally. Our phylogenetic analysis supports multiple independent introductions of HPAI from North America into South American countries for which sequence data was available at the time of this study, including Peru, Ecuador, and Venezuela. This contrasts with the single introduction of HPAI from Eurasia to North America observed earlier in 2021. Although North America is the primary source of HPAI for South America’s initial HPAI outbreaks, South American countries are likely to become more important sources for regional HPAI outbreaks as the virus establishes locally. After observing HPAI A/H5N1 reassort repeatedly within North American viruses, it is possible that the virus will continue to evolve in South America by mutation and reassortment with the genetically distinct South American AIV lineage that is commonly detected in Argentina17 and Chile18. For this reason, and given the uncharacterised variable sites reported here, there is an urgent need to establish pipelines for efficient real-time genomic sequencing of HPAI to track viral evolution and spread across Peru and other countries in South America, as well as funding to support characterisation of possible new mutations.
The impact of HPAI A/H5N1 on the morbidity and mortality of endangered species like the Andean condor (Vultur gryphus)38, the Humboldt penguin (Spheniscus humboldti)39, the marine otter (Lontra felina)40, and many others is concerning. The Peruvian coast is one of the few places in the world where scavenger condors feed on dead marine animals41,42, putting them at risk of infection if they consume contaminated carcasses. The recent infection and death of several critically endangered California condors (Gymnogyps californianus) in the United States highlights these risks43. The Humboldt penguin lives in large colonies and shares space with guano birds and sea lions, and although the single individual tested here was not positive for HPAI A/H5N1, Chile is observing record mortality in this species, now reaching 10% of their population, with four confirmed positive to HPAI A/H5N144. Similarly, the marine otter40 is an aquatic mustelid that inhabits the same rocky shores of Peru and Chile inhabited by sea lions. HPAI A/H5N1 has been confirmed previously in otters2 and recently in two animals from Chile44, and although none were tested here, we have anecdotal knowledge of otter carcasses washing onto Peruvian shores. Fortunately, marine otters do not live in large groups45, which might limit intraspecies contagion. However, direct mammal to mammal transmission has been suggested as a possible explanation for an outbreak in a Spanish farm among breeding mink, another mustelid, though evidence of mammal-to-mammal transmission was inconclusive9. Various cases of HPAI A/H1 in mammals, like the dolphin case reported here and others in Chile44, along with the massive South American sea lion die-off in Peru supports direct mammal to mammal transmission as a possible viral dissemination route, but confirmation will require further investigations with larger samples sizes, deeper genomic analyses and epidemiologic data that are often difficult to come by in outbreak settings. Peruvian pelicans have also suffered massive die-offs at the beginning of the outbreak in Peru, later followed by Guanay cormorants and Peruvian boobies. Peruvian pelicans are considered near-threatened worldwide46, but in Peru they are classified as endangered due to recent large population declines resulting from severe El Niño events associated with global warming and overfishing of their main source of food, the Peruvian anchovy15. For these reasons, efforts are urgently needed to assess the impact of ongoing mass mortalities on Peruvian marine sea birds and mammals.
Finally, an even larger concern is the possibility of spillover into human populations, as has been already documented22,47,48, followed by massive human-to-human transmission. Previous human cases have resulted in fatalities49, which has led the WHO to declare that the current zoonotic threat from HPAI A/H5N1 remains elevated and that member states should remain vigilant and consider mitigation steps to reduce human exposure50. In Peru, the outbreak occurred along the Pacific coast and during the austral summer, when many people go to the beach. It is not uncommon for beachgoers (and their pets) to interact with sick and disoriented animals without any knowledge of the risks, or for free-roaming dogs in rural and semi-rural coastal areas to encounter sick or dead animals as they scavenge for food. This has led government authorities to relocate live animals that show up in places where they do not belong, or to euthanize sick individuals and appropriately dispose of their carcasses. However, both the pelican and sea lion die-offs have been so massive that it has been very challenging for the authorities to respond in a timely manner. The recent confirmed HPAI H5N1 human case in Chile22 has been considered an environmental exposure due to an abundance of dead wildlife in the area close to the patient’s residence, as the human virus sequenced was identical to local wild bird strains. Heightened public awareness campaigns are needed, including educating the public to avoid contact with infected animals51. Animal workers, particularly municipal personnel tasked with cleaning duties, need additional training in the proper use of personal protective equipment, and on management and disposal of infected carcasses51,52. People in contact with sick and dead animals infected with HPAI A/H5N1 are at risk of infection, and human cases could be missed in the absence of active and obvious human-to-human transmission. Syndromic surveillance based on PCR-based diagnostics and/or serology of suspected cases in people working in wildlife and poultry outbreak response is highly recommended for early detection of zoonotic transmission. Subsequent serologic studies of people in close contact with infected animals, especially in outbreak settings, can also inform on the extent of zoonotic spillover and direct future surveillance at important animal-human interfaces. This approach promotes opportunities to strengthen ties among governmental institutions at all levels, expert researchers within academia, and the larger scientific community operationalizing a One Health approach, which are essential for epidemic and pandemic preparedness. The work reported here exemplifies such collaboration, and although many of the variable genomic sites identified in Peruvian H5N1 viruses will need further validation with both downstream in-vitro and in-vivo studies, this study of HPAI A/H5N1 in marine birds and mammals in Peru represents a significant contribution to the understanding of a highly pathogenic virus with serious pandemic potential for humans.
Starting on November 22, 2022, samples were collected opportunistically from both live and dead animals by trained veterinarians from the Peruvian wildlife service (Servicio Nacional Forestal y de Fauna Silvestre, SERFOR) and the NGO Wildlife Conservation Society Perú, under permit 2023-009134, using full personal protective equipment and following standard protocols for cleaning, disinfection, and disposal of hazardous waste52. Some severely clinically ill animals were humanely euthanized53. Briefly, animals were chemically immobilised by intramuscular administration of ketamine and xylazine, followed by intravenous or intracardiac injection of T-61TM (MSD Animal Health USA)53. Samples and tissues were collected into cryovials containing 0.5 mL of DNA/RNA Shield (Zymo R1200 125) or viral transport media (VTM), and transported at 4 °C within 1–10 days of collection to the laboratory for testing. Nucleic acids were extracted using Quick-DNA/RNA Viral Extraction Kits (Zymo D7021) and tested for influenza A by RT-qPCR using CDC protocols54 on a Bio-Rad CFX96 instrument.
Samples positive for influenza A by RT-qPCR were subtyped using a combination of directed amplification with universal primers targeting conserved genomic regions55,56,57, followed by next-generation sequencing (NGS). All primer sequences used are provided in Supplementary Table 3. Briefly, RNA samples were reverse transcribed using Superscript IV (Invitrogen 18090050) and amplified using Q5 High-fidelity DNA polymerase (NEB M0491L). Amplification products were prepared into barcoded sequencing libraries using DNA Prep Kits (Illumina 20060059) and Nextera DNA UD Indexes (Illumina 20027215). The resulting libraries were quality controlled using High Sensitivity DNA kits (Agilent 5067–4626) on a Bioanalyzer 2100 instrument. Libraries were normalised to 4 nM each, pooled, re-quantified using Qubit 1x dsDNA HS Kits (Invitrogen Q33230), and sequenced using High Output Sequencing Kits (Illumina FC-420-1003) on an Illumina MiniSeq instrument.
Illumina paired-end raw sequence data was pre-processed to trim sequencing adaptors and filter out low quality/low complexity reads (Phred scores <Q20, 75 bp minimum length) using Geneious Prime 2023.0.4 and BBDuk58. Pre-processed reads were then filtered by reference-mapping to various HA (H1, H2, H3, H5, H7) and NA (N1, N2, N3, N5, N7) references (Accession Numbers: NC_026433.1, NC_007374.1, NC_007366.1, NC_007362.1, NC_026425.1, AF144304.1, NC_007382.1, OP806485, MF046172.2, OP723829.1), and to single references in the case of all other segments (Accession Numbers MT624412.1, MN254461.1, KY635563.1, MT825070.1, MW110227.1, MT982385.1). Filtered reads were then re-assembled de-novo using SPAdes59 to generate complete genomes whenever possible, and to further confirm subtyping by BLAST. All sequences have been deposited in GenBank (Accession Numbers OQ550419-OQ550478, and OQ925704-OQ925729).
HPAI H5N1 reference sequences were obtained from GenBank and GISAID, and together with the sequences generated here, they were concatenated and aligned using MAFFT60. Maximum likelihood (ML) trees were prepared with FastTree61 incorporating a general time-reversible (GTR) model of nucleotide substitution with gamma-distributed rate variation among sites. Trees were annotated using the WHO/OIE/FAO nomenclature for highly pathogenic avian influenza A H562. To place the Peruvian viruses in a global context, we downloaded (on February 14, 2023) an additional background dataset of influenza A virus genomes that included all sequences from avian and mammalian H5 viruses submitted to GISAID since January 1, 2021. To keep the dataset up to date, additional South American sequences were added on April 21, 2023. Partial sequences were excluded. Phylogenetic relationships were inferred for each of the eight genome segments using the ML methods available in IQ-Tree 263 with a GTR model and a gamma distribution as described above. Due to the size of the dataset, we used the high-performance computational capabilities of the Biowulf Linux cluster at the National Institutes of Health (http://biowulf.nih.gov). To assess the robustness of each node, a bootstrap resampling process was performed with 1000 replicates. Finally, mutation analysis was done using the CDC H5N1 genetic changes inventory for SNP analysis and various other previously published mutations of concern23,24.
To examine the evolution of the clade of reassortant R6 viruses in Peru and other locations in greater detail, we performed a time-scaled Bayesian analysis on a dataset of 98 PB2 sequences from the R6 clade. The dataset included all South American sequences, all mammalian sequences, and a representative subsample of North American avian sequences, including one representative virus from each cluster of viruses from the same US or Canadian state that were sampled from the same avian species on the same date. We used the Markov chain Monte Carlo (MCMC) method available in the BEAST package, v1.10.464, again using the NIH Biowulf Linux cluster, and a Bayesian non-parametric (SkyGrid) demographic model65, with a GTR model of nucleotide substitution with gamma-distributed rate variation among sites. Each tip was assigned a location state and we performed a phylogeographic discrete trait analysis66 to examine the routes of viral spatial dissemination between North America and South America, within South America, and within Peru. The MCMC chain was run separately four times for each dataset using the BEAGLE 367 library to improve computational performance, until all parameters reached convergence, as assessed visually using Tracer v.1.7.268. At least 10% of the chain was removed as burn-in, and runs for the same dataset were combined using LogCombiner v1.10.4. A MCC tree was summarised using TreeAnnotator v.1.10.4 and visualised in FigTree v1.4.4. SPREAD 469 was used to visualise the spatial-temporal reconstruction. XMLs, MCC and ML trees, and GISAID acknowledgement tables are available in a GitHub repository (https://github.com/mostmarmot/Peru_AIV/).
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
We gratefully acknowledge the authors and both originating and submitting laboratories of the sequences from GISAID’s EpiFlu™ Database on which this research is based. The GISAID server was unable to download a large table of 4000+ bird viruses used in our analyses, but tables for all South American and mammalian viruses used are provided in the supplemental materials. We also provide GenBank accession numbers with hyperlinks for all the sequences generated as part of this study in Supplementary Table 4. In addition, XMLs, MCC and ML trees, and GISAID acknowledgement tables are also available in a GitHub repository (https://github.com/mostmarmot/Peru_AIV/).
For NGS data analysis we used a suite of bioinformatic tools, including BBDuk and SPAdes, that are contained within the Geneious Prime 2023.0.4 package. In addition, XMLs, MCC and ML trees, and GISAID acknowledgement tables are available in a GitHub repository (https://github.com/mostmarmot/Peru_AIV/).
Ramey, A. M. et al. Highly pathogenic avian influenza is an emerging disease threat to wild birds in North America. J. Wildlife Manag. 86, e22171 (2022).
Article Google Scholar
Vreman, S. et al. Zoonotic mutation of highly pathogenic avian influenza H5N1 virus identified in the brain of multiple wild carnivore species. Pathogens 12, 168 (2023).
Article CAS PubMed PubMed Central Google Scholar
Alders, R., Awuni, J. A., Bagnol, B., Farrell, P. & de Haan, N. Impact of avian influenza on village poultry production globally. Ecohealth 11, 63–72 (2014).
PubMed Google Scholar
Ministerio de Salud del Peru, M. Alerta Epidemiológica: Seguimiento de la Epizootia de Influenza Aviar, Tipo A, Subtipo H5 en aves silvestres y aves de traspatio en el país. Código: AE 002-2023, https://www.dge.gob.pe/epipublic/uploads/alertas/alertas_20232_31_094534.pdf (2023).
Kalthoff, D., Globig, A. & Beer, M. (Highly pathogenic) avian influenza as a zoonotic agent. Vet. Microbiol. 140, 237–245 (2010).
Article PubMed Google Scholar
Caliendo, V. et al. Transatlantic spread of highly pathogenic avian influenza H5N1 by wild birds from Europe to North America in 2021. Sci. Rep. 12, 11729–11729 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Elsmo, E. J. et al. Pathology of natural infection with highly pathogenic avian influenza virus (H5N1) clade 2.3.4.4b in wild terrestrial mammals in the United States in 2022. bioRxiv, https://doi.org/10.1101/2023.03.10.532068 (2023).
Alkie, T. N. et al. A threat from both sides: multiple introductions of genetically distinct H5 HPAI viruses into Canada via both East Asia-Australasia/Pacific and Atlantic flyways. Virus Evol. 8, https://doi.org/10.1093/ve/veac077 (2022).
Agüero, M. et al. Highly pathogenic avian influenza A(H5N1) virus infection in farmed minks, Spain, October 2022. Euro. Surveill. 28, https://doi.org/10.2807/1560-7917.Es.2023.28.3.2300001 (2023).
Ministerio de Salud del Peru, M. Alerta Epidemiológica: Epizootia de Influenza Aviar, Tipo A, Subtipo H5 en aves silvestres en el litoral del país. Código: AE 028-2022, https://www.dge.gob.pe/epipublic/uploads/alertas/alertas_202228_24_145625.pdf (2022).
Ministerio de Salud del Peru, M. Sala de influenza aviar, https://www.dge.gob.pe/influenza-aviar-ah5/ (2023).
Servicio Nacional de Áreas Naturales Protegidas por el Estado, S. Sernanp despliega protocolo de monitoreo ante casos de aves y lobos marinos afectados por influenza aviar en áreas naturales protegidas, https://www.gob.pe/institucion/sernanp/noticias/697084-sernanp-despliega-protocolo-de-monitoreo-ante-casos-de-aves-y-lobos-marinos-afectados-por-influenza-aviar-en-areas-naturales-protegidas (2023).
Miloslavich, P. et al. Marine biodiversity in the Atlantic and Pacific coasts of South America: knowledge and gaps. PLoS One 6, e14631 (2011).
Article ADS CAS PubMed PubMed Central Google Scholar
Schulenberg, T. S., Stotz, D. F., Lane, D. F., O’Neill, J. P. & Parker III, T. A. Aves de Perú. 1 (Princeton University Press, 2010).
Servicio Nacional Forestal y de Fauna Silvestre, S. Libro Rojo de la Fauna Silvestre Amenazada del Perú, 1st edn (2018).
Nelson, M. I. et al. The genetic diversity of influenza A viruses in wild birds in Peru. PLoS One 11, e0146059 (2016).
Article PubMed PubMed Central Google Scholar
Rimondi, A. et al. Evidence of a fixed internal gene constellation in influenza A viruses isolated from wild birds in Argentina (2006–2016). Emerg. Microbes Infect. 7, 194 (2018).
Article PubMed PubMed Central Google Scholar
Jiménez-Bluhm, P. et al. Wild birds in Chile Harbor diverse avian influenza A viruses. Emerg. Microbes Infect. 7, 44 (2018).
Article PubMed PubMed Central Google Scholar
Lee, D. H. et al. Highly pathogenic Avian Influenza Viruses and generation of novel reassortants, United States, 2014–2015. Emerg. Infect. Dis. 22, 1283–1285 (2016).
Article CAS PubMed PubMed Central Google Scholar
Harfoot, R. & Webby, R. J. H5 influenza, a global update. J. Microbiol. 55, 196–203 (2017).
Article PubMed Google Scholar
Bruno, A. et al. Phylogenetic analysis reveals that the H5N1 avian influenza A outbreak in poultry in Ecuador in November 2022 is associated with the highly pathogenic clade 2.3.4.4b. Int. J. Infect. Dis. 133, 27–30 (2023).
Article CAS PubMed Google Scholar
World Health Organization, W. Human Infection caused by Avian Influenza A (H5N1) – Chile, https://www.who.int/emergencies/disease-outbreak-news/item/2023-DON461 (2023).
Centers for Disease Control Prevention, C. H5N1 genetic changes inventory: a tool for international surveillance | Avian Influenza (Flu), https://www.cdc.gov/flu/avianflu/h5n1-genetic-changes.htm (2015).
Suttie, A. et al. Inventory of molecular markers affecting biological characteristics of avian influenza A viruses. Virus Genes 55, 739–768 (2019).
Article CAS PubMed PubMed Central Google Scholar
Bordes, L. et al. Highly pathogenic avian influenza H5N1 virus infections in wild red foxes (Vulpes vulpes) show neurotropism and adaptive virus mutations. Microbiol. Spectrum 11, e02867–02822 (2023).
Article Google Scholar
Bevins, S. N. et al. Intercontinental movement of highly pathogenic avian influenza A(H5N1) Clade 2.3.4.4 virus to the United States, 2021. Emerg. Infect. Dis. 28, 1006–1011 (2022).
Article CAS PubMed PubMed Central Google Scholar
Crespo, E. A., de Oliveira, L. R. & Sepúlveda, M. in Ecology and conservation of pinnipeds in Latin America (eds G. Heckel & Y. Schramm) 93–118 (Springer International Publishing, 2021).
Reperant, L. A., Rimmelzwaan, G. F. & Kuiken, T. Avian influenza viruses in mammals. Rev. Sci. Tech. 28, 137–159 (2009).
Article CAS PubMed Google Scholar
Thiebault, A., Semeria, M., Lett, C. & Tremblay, Y. How to capture fish in a school? Effect of successive predator attacks on seabird feeding success. J. Anim. Ecol. 85, 157–167 (2016).
Article PubMed Google Scholar
Weimerskirch, H., Bertrand, S., Silva, J., Bost, C. & Peraltilla, S. Foraging in Guanay cormorant and Peruvian booby, the major guano-producing seabirds in the Humboldt Current System. Marine Ecol. Progr. Ser. 458, 231–245 (2012).
Article ADS Google Scholar
Martin, G., Becker, D. J. & Plowright, R. K. Environmental persistence of influenza H5N1 is driven by temperature and salinity: insights from a bayesian meta-analysis. Front. Ecol. Evol. 6, https://doi.org/10.3389/fevo.2018.00131 (2018).
US Department of Agriculture, A. 2022–2023 detections of highly pathogenic avian influenza in mammals, https://www.aphis.usda.gov/aphis/ourfocus/animalhealth/animal-disease-information/avian/avian-influenza/hpai-2022/2022-hpai-mammals (2023).
Reperant, L. A. et al. Highly pathogenic avian influenza virus (H5N1) infection in red foxes fed infected bird carcasses. Emerg. Infect. Dis. 14, 1835–1841 (2008).
Article PubMed PubMed Central Google Scholar
Kurmi, B. et al. Survivability of highly pathogenic avian influenza H5N1 virus in poultry faeces at different temperatures. Indian J. Virol. 24, 272–277 (2013).
Article PubMed PubMed Central Google Scholar
Kandun, I. N. et al. Chicken faeces garden fertilizer: possible source of human avian influenza H5N1 infection. Zoonoses Public Health 57, 285–290 (2010).
Article CAS PubMed Google Scholar
Hall, J. S. et al. Experimental challenge and pathology of highly pathogenic avian influenza virus H5N1 in dunlin (Calidris alpina), an intercontinental migrant shorebird species. Influenza Other Respir. Viruses 5, 365–372 (2011).
Article PubMed PubMed Central Google Scholar
Reperant, L. A. et al. Highly pathogenic avian influenza virus H5N1 infection in a long-distance migrant shorebird under migratory and non-migratory states. PLOS One 6, e27814 (2011).
Article ADS CAS PubMed PubMed Central Google Scholar
International Union for Conservation of Nature, I. Vultur gryphus, https://www.iucnredlist.org/species/22697641/181325230 (2020).
International Union for Conservation of Nature, I. Spheniscus humboldti, https://www.iucnredlist.org/species/22697817/182714418 (2020).
International Union for Conservation of Nature, I. Lontra felina, https://www.iucnredlist.org/species/12303/215395045 (2022).
Wallace, M. P. & Temple, S. A. Competitive interactions within and between species in a guild of avian scavengers. The Auk 104, 290–295 (1987).
Article Google Scholar
Stucchi, M. Presencia del Cóndor Andino (Vultur gryphus) en la costa peruana. Boletín de Lima 155, 124–133 (2013).
Google Scholar
National Parks Service, N. Highly pathogenic avian influenza confirmed as cause of three california condor mortalities in Arizona, https://www.nps.gov/orgs/1207/highly-pathogenic-avian-influenza-confirmed-as-cause-of-three-california-condor-mortalities-in-arizona.htm (2023).
Servicio Nacional de Pesca y Acuicultura, S. Influenza aviar, http://www.sernapesca.cl/influenza-aviar (2023).
Valqui, J. The marine otter Lontra felina (Molina, 1782): A review of its present status and implications for future conservation. Mamm. Biol. 77, 75–83, https://doi.org/10.1016/j.mambio.2011.08.004 (2012).
Article Google Scholar
International Union for Conservation of Nature, I. Pelecanus thagus, https://www.iucnredlist.org/species/22697619/132596827 (2018).
World Health Organization, W. Human infection with avian influenza A(H5) viruses – Avian Influenza Weekly Update Number 881, https://www.who.int/docs/default-source/wpro—documents/emergency/surveillance/avian-influenza/ai_20230203.pdf (2023).
Centers for Disease Control Prevention, C. U.S. case of human avian influenza A(H5) virus reported, https://www.cdc.gov/media/releases/2022/s0428-avian-flu.html (2022).
Lai, S. et al. Global epidemiology of avian influenza A H5N1 virus infection in humans, 1997–2015: a systematic review of individual case data. Lancet Infect. Dis. 16, e108–e118 (2016).
Article PubMed PubMed Central Google Scholar
World Health Organization, W. Assessment of risk associated with recent influenza A(H5N1) clade 2.3.4.4b viruses, https://www.who.int/publications/m/item/assessment-of-risk-associated-with-recent-influenza-a(h5n1)-clade-2.3.4.4b-viruses (2022).
Geraci, J. R. & Lounsbury, V. J. Marine mammals ashore: a field guide for strandings (National Aquarium in Baltimore, 2005).
Miller, L. P., Miknis, R. A. & Flory, G. A. Carcass management guidelines: effective disposal of animal carcasses and contaminated materials on small to medium-sized farms, Vol. 23 (Food & Agriculture Org., 2020).
Underwood, W. & Anthony, R. AV MA guidelines for the euthanasia of animals: 2020 edition. Retrieved on March 2013, 2020–2021 (SERFOR (Servicio Nacional Forestal y de Fauna Silvestre), 2020).
Shu, B. et al. Multiplex real-time reverse transcription PCR for influenza A virus, influenza B virus, and severe acute respiratory syndrome coronavirus 2. Emerg. Infect. Dis. 27, 1821–1830 (2021).
Article CAS PubMed PubMed Central Google Scholar
Hoffmann, E., Stech, J., Guan, Y., Webster, R. G. & Perez, D. R. Universal primer set for the full-length amplification of all influenza A viruses. Arch. Virol. 146, 2275–2289 (2001).
Article CAS PubMed Google Scholar
Zhou, B. et al. Single-reaction genomic amplification accelerates sequencing and vaccine production for classical and swine origin human influenza A viruses. J. Virol. 83, 10309–10313 (2009).
Article CAS PubMed PubMed Central Google Scholar
Hurtado, R. et al. Molecular characterization of subtype H11N9 avian influenza virus isolated from shorebirds in Brazil. PLoS One 10, e0145627 (2015).
Article PubMed PubMed Central Google Scholar
Bushnell, B. BBMap short read aligner, and other bioinformatic tools, https://sourceforge.net/projects/bbmap/ (2022).
Bankevich, A. et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477 (2012).
Article MathSciNet CAS PubMed PubMed Central Google Scholar
Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780 (2013).
Article CAS PubMed PubMed Central Google Scholar
Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol. Biol. Evol. 26, 1641–1650 (2009).
Article CAS PubMed PubMed Central Google Scholar
Smith, G. J. & Donis, R. O. Nomenclature updates resulting from the evolution of avian influenza A(H5) virus clades 2.1.3.2a, 2.2.1, and 2.3.4 during 2013–2014. Influenza Other Respir. Viruses 9, 271–276 (2015).
Article PubMed PubMed Central Google Scholar
Nguyen, L. T., Schmidt, H. A., von Haeseler, A. & Minh, B. Q. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274 (2015).
Article CAS PubMed Google Scholar
Suchard, M. A. et al. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 4, vey016 (2018).
Article PubMed PubMed Central Google Scholar
Gill, M. S. et al. Improving Bayesian population dynamics inference: a coalescent-based model for multiple loci. Mol. Biol. Evol. 30, 713–724 (2013).
Article CAS PubMed Google Scholar
Lemey, P., Rambaut, A., Drummond, A. J. & Suchard, M. A. Bayesian phylogeography finds its roots. PLoS Comput. Biol. 5, e1000520 (2009).
Article ADS MathSciNet PubMed PubMed Central Google Scholar
Ayres, D. L. et al. BEAGLE 3: improved performance, scaling, and usability for a high-performance computing library for statistical phylogenetics. Syst. Biol. 68, 1052–1061 (2019).
Article PubMed PubMed Central Google Scholar
Rambaut, A., Drummond, A. J., Xie, D., Baele, G. & Suchard, M. A. Posterior summarization in bayesian phylogenetics using tracer 1.7. Syst. Biol. 67, 901–904 (2018).
Article CAS PubMed PubMed Central Google Scholar
Nahata, K. D. et al. SPREAD 4: online visualisation of pathogen phylogeographic reconstructions. Virus Evol. 8, veac088 (2022).
Article PubMed PubMed Central Google Scholar
Download references
Authors are employees of the Peruvian Ministerio de Desarrollo Agrario y Riego (MIDAGRI), of the Pontificia Universidad Catolica del Peru (PUCP), of the Wildlife Conservation Society, of the University of California – Davis, and of the US Government. This work was prepared as part of their official duties, with additional support from the National Institute of Allergy and Infectious Diseases at the National Institutes of Health under award U01AI151814 (M.L., A.G.G., B.M.S., D.J., P.B., C.K.J., C.C.M. and M.M.U.); from the Gordon and Betty Moore Foundation through Grant GBMF10809 to WCS Peru (P.C.C.); from the Intramural Research Program of the US National Library of Medicine at the NIH and the Centers of Excellence for Influenza Research and Response, National Institute of Allergy and Infectious Diseases, National Institutes of Health (NIH), Department of Health and Human Services, under contract 75N93021C00014 (M.I.N.). The funders had no role in study design, sample collection, data collection and analysis, decision to publish, or preparation of the article. The authors thank Sherilym Castillo, Esthefany Ruíz, Pedro Ramírez, Max Guerra, Antero Martínez, Wendy Rojas, Renato Colán, José Cerón, and Karina Espinoza from SERFOR for support during field work; Julio Reyes and Joe Macalupú from IMARPE for support with dolphin species identification and sampling permissions, respectively; Mariana Peralta from ATFFS-Arequipa, as well as Ferdinand Ortiz Mamani, Daniel Alama and Paul Cieza from SENASA for support with sample conservation and transport from field to lab; Miryam Quevedo from UNMSM for providing equipment for fieldwork; and Rosa Vento and Jorge Martínez from WCS for general logistical support.
Laboratorio de Genómica, Pontificia Universidad Católica del Perú (PUCP), Lima, Peru
Mariana Leguia, Alejandra Garcia-Glaessner, Breno Muñoz-Saavedra, Diana Juarez & Patricia Barrera
EpiCenter for Emerging Infectious Disease Intelligence, Centers for Research in Emerging Infectious Diseases, Lima, Peru
Mariana Leguia, Alejandra Garcia-Glaessner, Breno Muñoz-Saavedra, Diana Juarez, Patricia Barrera, Carlos Calvo-Mac, Christine K. Johnson & Marcela M. Uhart
Servicio Nacional Forestal y de Fauna Silvestre (SERFOR), Ministerio de Desarrollo Agrario y Riego (MIDAGRI) del Perú, Lima, Peru
Javier Jara, Walter Silva, Karl Ploog, Lady Amaro & Jesus Lescano
Wildlife Conservation Society (WCS) – Perú, Lima, Peru
Paulo Colchao-Claux
One Health Institute, School of Veterinary Medicine, University of California, Davis, CA, USA
Christine K. Johnson & Marcela M. Uhart
National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health (NIH), Bethesda, MD, USA
Martha I. Nelson
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Study Design: M.L. and J.L. Field work and sample collection: J.J., W.S., K.P., L.A., P.C.C. and J.L. Lab work and sample processing: A.G.G., D.J., P.B. Bioinformatics processing and data analysis: M.L., A.G.G., B.M.S., D.J., P.B., M.I.N. Manuscript preparation: M.L., A.G.G., C.C.M., L.A., C.K.J., M.M.U., M.I.N. Manuscript editing: all authors. Review and approval of final manuscript: all authors.
Correspondence to Mariana Leguia.
The authors declare no competing interests.
Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Leguia, M., Garcia-Glaessner, A., Muñoz-Saavedra, B. et al. Highly pathogenic avian influenza A (H5N1) in marine mammals and seabirds in Peru. Nat Commun 14, 5489 (2023). https://doi.org/10.1038/s41467-023-41182-0
Download citation
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-023-41182-0
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Collection
Focus
Advertisement
© 2023 Springer Nature Limited
Sign up for the Nature Briefing: Translational Research newsletter — top stories in biotechnology, drug discovery and pharma.