Mortality in Sea Lions is associated with the introduction of the H5N1 clade 2.3.4.4b virus in Brazil, October 2023: Whole genome sequencing and phylogenetic analysis

doi:10.21203/rs.3.rs-3793926/v1

Download PDF

Research Article

Mortality in Sea Lions is associated with the introduction of the H5N1 clade 2.3.4.4b virus in Brazil, October 2023: Whole genome sequencing and phylogenetic analysis

https://doi.org/10.21203/rs.3.rs-3793926/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Clade 2.3.4.4b highly pathogenic avian influenza (HPAI) H5N1 virus was detected in the South American sea lions found dead in Santa Catarina, Brazil, in October 2023. Whole genome sequencing and comparative phylogenetic analysis were conducted to investigate the origin, genetic diversity, and zoonotic potentials of the H5N1 viruses. The H5N1 viruses belonged to the genotype B3.2 of clade 2.3.4.4b H5N1 virus, which was identified in North America and disseminated to South America. They have acquired new amino acid substitutions related to mammalian host affinity. Our study provides insights into the genetic landscape of HPAI H5N1 viruses in Brazil, highlighting the continuous evolutionary processes contributing to their possible adaptation to mammalian hosts.

High Pathogenicity Avian Influenza

Whole Genome Sequencing

Marine mammal

South America

Virus

Public Health

The recent panzootic of the clade 2.3.4.4 b highly pathogenic avian influenza (HPAI) A (H5N1) virus caused massive outbreaks including wild birds, poultry, and mammalian species all around the globe during 2022 and 2023[1]. Particularly, mass mortality of coastal sea birds [2–4] and marine mammals caused by the HPAI H5N1 virus were reported from South America including Peru and Chile [4, 5]. In Brazil, the virus was identified for the first time in May 2023 in sea birds in the Southeast region and in the state of Santa Catarina in June 2023. Increase in detections of H5N1 virus in land and marine mammals are of wildlife and public health concern, as the viruses can rapidly evolve and adapt to mammals [3]. Here, we report the detection of clade 2.3.4.4 b HPAI H5N1 virus from South American sea lions (Otaria flavescens) in Santa Catarina, Brazil, during October of 2023. We conducted whole genome sequencing (WGS) and comparative phylogenetic analysis, to investigate the molecular epidemiology, genetic diversity, and mammalian adaptation markers of the H5N1 viruses.

A sudden surge of mortality in South American sea lions (Otaria flavescens) was recorded in October in the State of Santa Catarina, Brazil when compared to that of the previous years (Fig. 1). Most of the carcasses were at advanced decomposition without apparent gross lesion (SIMBA, 2023). To investigate if this unusual mortality was related to H5N1 virus, tissue samples of the digestive and respiratory tract were collected from October 10th to 15th, 2023 and shipped to the laboratory. Additionally, to investigate the potential transmission among viruses infecting birds and marine mammals in the same geographic area, we collected cloacal swab samples from 190 seabirds (Fig. 1). Birds were found weakened on the beaches and taken for treatment in the Santa Catarina veterinary care network, stabilization, and rehabilitation centers. They were from five different orders and 17 species as follows: 101 Charadriiformes (Larus dominicanus, Rynchops niger, Sterna hirundinacea, Sterna hirundo and Thalasseus acuflavidus), 75 Sphenisciformes (Spheniscus magellanicus), 23 Suliformes (Fregata magnificens, Phalacrocorax brasilianus and Sula leucogaster), eight Procellariiformes (Thalassarche melanophris, Pterodroma sp., Puffinus puffinus and Procellaria aequinoctialis), and two Pelecaniformes (Phimosus infuscatus). Samples were collected by the Santos Basin Beach Monitoring Project, PMP-BS, under license ABIO 240/2015. This project has been continuously monitoring nearly the entire coastline of Santa Catarina state since August 2015, and the analyzed samples were from June 20th to December 4th, 2023. Total RNA was extracted using a MagMax core nucleic acid purification kit (ThermoFisher, California, USA) and tested for RT-qPCR targeting the matrix gene [6]. Data on the positive birds for HPAI H5N1 in the area was also collected from public data provided by the Brazilian Ministry of Agriculture (https://mapa-indicadores.agricultura.gov.br/publico/extensions/SRN/SRN.html).

For each sample, amplification of all 8 target genes was conducted using the SuperScript IV Reverse Transcriptase kit (Thermo Fisher Scientific, USA) with the universal primers as previously described [7]. For next generation sequencing (NGS), library preparation was done using the Nextera XT DNA kit (Illumina) and run on NextSeq 550 (Illumina, San Diego, California, USA) with NextSeq 500/550 Mid Output Kit v2 (300-cycles). Genome assembly was performed and visualized using the map to reference tool in Geneious Prime software (https://www.geneious.com). Assembled consensus genome sequences were uploaded in the EpiFlu database of Global Initiative on Sharing All Influenza Data (GISAID: https://www.gisaid.org), under the accession numbers: (EPI_ISL_18669979 and EPI_ISL_18669996)

The complete coding sequences of the viruses were submitted to the GISAID BLAST database (https://gisaid.org/) for the query of top 250 BLAST results of all eight genes Top 250 BLAST result sequences were downloaded and aligned using the MAFFT versions 7.520 (https://mafft.cbrc.jp/alignment/software/) and 100% homology sequences were eliminated using the ElimDupes (https://www.hiv.lanl.gov/content/sequence/elimdupesv2/elimdupes.html). Maximum-likelihood (ML) trees were constructed using the RaxML v8.0 software using rapid bootstrapping replicate of 1000 trees to search for best-scoring ML tree with general time-reversible + Gamma nucleotide substitution model. ML trees were visualized using Interactive Tree of Life v6 (https://itol.embl.de/, accessed on 11 December 2023)[8]. Genotype of the viruses were determined by examining ML phylogeny topology and query result using the GenoFlu (https://github.com/USDA-VS/GenoFLU) [9]. Molecular markers associated with mammalian adaptation were screened using the FluSurver mutations app (http://flusurver.bii.a-star.edu.sg) and previous publications [10]. For the estimation of the time to the most recent ancestor (tMRCA) harboring the mutation, concatenated complete genome sequences of South American H5N1 viruses of were used in combination with Bayesian phylogenetic analysis using BEAST package 1.10.4 as previously described [11, 12]. Five out of eight samples from the respiratory tract of four sea lions were tested positive by the RT-qPCR test using the AIV matrix gene as the target. Among the positive samples, we sequenced two samples that had a cycle threshold (Ct) value ranging below 25. We obtained complete coding genome sequences of two viruses and named as A/South American sea lion/Brazil/OF-358R/2023 and A/South American sea lion/Brazil/OF-359R/2023 (hereafter, OF-358R/23 and OF-359R/23). Only one sample from the digestive tract (animal OF-359) was positive by the RT- qPCR with a Ct value of 33. None of the asymptomatic birds collected by our group tested positive for the RT- qPCR specific for AIV. Ten outbreaks in birds and one in sea lions were reported in Santa Catarina from September to October by the Official Veterinary Service (OVS) of the Ministry of Agriculture (Fig. 2). Maximum-likelihood (ML) phylogenies constructed using the internal genes revealed that the OF-358R/23 and OF-359R/23 viruses clustered with the previously reported South American clade 2.3.4.4 b H5N1 HPAIVs (Supplementary Figure S1). Genotype BLAST search using the GenoFlu also identified the eight genes as genotype B3.2 with less than 1% nucleotide difference [9]. No evidence of genetic reassortment was observed in these viruses. Hemagglutinin (HA) cleavage site of the viruses possessed a HPAI motif, PLREKRRKR/GLF, which is identical to the recent clade 2.3.4.4 b HPAI H5N1 viruses found in South America.

Amino acid substitutions related to mammalian adaptation were observed in the genomes of OF-358R/23 and OF-359R/23 viruses. Notably, amino acid substitutions of Q591K and D701N in PB2 gene, R57Q in PA gene, and V226T in NS were newly acquired when compared to the HPAI H5N1 isolates previously reported from wild birds in Brazil (Table 1) [13, 14]. These amino acid substitutions found in the H5N1 viruses are known to increase virulence, transmissibility, or replication in mammalian host species.

OF-358R/23 and OF-359R/23 viruses clustered with the human isolate originating from Chile and sea lion isolates originating from Peru and Chile (Fig. 3 and Supplementary Figs. 1 and 2). For better resolution of the estimation of tMRCA of the unique amino acid substitutions, concatenated genome sequences of all eight influenza genes were used. R57Q in PA gene first arose around September 23rd, 2022 (95% HPD: September 5th, 2022 to October 11th, 2022). D701N in PB2 gene most likely occurred twice in a separate instance, the first occurrence was around December 26th, 2022 (95% HPD: December 8th, 2022 to January 13th, 2023), the second occurrence was also in sea lion isolate from Peru around February 4th, 2023 (95% HPD: January 30th, 2023 to February 7th, 2023) (Fig. 3). V226T in NS gene was first observed during the same period as the first appearance of D701N in PB2. Q591K in PB2 emerged around January 27th, 2023 (95% HPD: January 9th, 2023 to February 13th, 2023) (Fig. 3).

Here, we sequenced and analyzed the complete genome sequences of clade 2.3.4.4 b HPAI H5N1 viruses detected from carcasses of South American sea lions. Phylogenetic analysis of the whole genome sequences indicated no evidence of reassortment from the genotype B3.2 initially introduced to South America. Most probable transmission route of the H5N1 virus is from North America to South America by wild birds, subsequently disseminated to marine mammals. Previous studies have suggested the possibility of establishment and transmission among marine mammals along South America American coastlines [3, 5]. It is unclear if the OF-358R/23 and OF-359R/23 viruses were originated from marine mammals or wild birds since there was no direct ancestral genome sequences available in public databases.

OF-358R/23 and OF-359R/23 showed close genetic relationship with the human isolate from Chile and Sea lion isolate from Peru, which also harbored amino acid substitutions such as Q591K, D701N, and V226T. Q57R emerged in the late 2022, at least a year before the detection of OF-358R/23 and OF-359R/23 (Fig. 3). As previously described, our analysis supported that D701N have emerged in two separate instances of sea lion virus clusters [3]. Convergent mutations like D701N highlight the need for the close monitoring of mutations arising in the viruses. Bayesian inference using the concatenated complete genome sequences revealed that the mutations in OF-358R/23 and OF-359R/23 were sustained after their emergence and has been circulating since the late 2022 and early 2023. The accumulation of these mutations may have aided in the adaptation of the H5N1 viruses to transmit among sea lion population along the coast of South America and increase their potential for zoonotic transmission as observed in Chile.

Comparative genome analysis with the recent mammalian origin viruses from South America revealed that the amino acid substitutions of Q591K and D704N in PB2, R57Q in PA, and V226T in NS were not universally found in all South American mammalian viruses, indicating that mutations associated with mammalian affinity may be accumulating. Increased detection of HPAI in marine mammals and accumulation of such mammalian molecular markers raise concern as increased affinity to mammalian host species pose as a significant threat to both mammalian wildlife and public health. Enhanced surveillance and continued monitoring effort is needed to elucidate the evidence of transmission among marine mammals and be wary of the potential spillover into human population.

Ethical approval and consent to participate

The authors confirm that the ethical policies of the journal, as noted on the journal’s author guidelines page, have been adhered to and the appropriate ethical review committee approval has been received. All applicable international, national, and/or institutional guidelines for the care and use of animals were followed.

Samples were collected under the scope of The Santos Basin Beach Monitoring Project (PMP/BS) conducted by PETROBRAS, as part of the environmental licensing process conducted by the Brazilian Institute of the Environment and Renewable Natural Resources (IBAMA) of the Brazilian Ministry of Environment, under license ABIO Nº 640/2015. Sample collections were authorized by the Ethics Committee on the Use of Animals of the University of São Paulo (CEUA/FZEA-USP: 1128141120), and National System for the Management of Genetic Heritage and Associated Traditional Knowledge (SISGEN A028CF5).

Consent for publication

Not applicable

Ethical Guidelines

Not applicable

Availability of data and materials

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: Global Initiative on Sharing All Influenza Data (GISAID: https://www.gisaid.org) EpiFlu database, under the accession numbers: EPI_ISL: 18669979 and 18669996

Funding

This work is funded by the Ministry of Science, Technology, and Innovation (MCTI-Brazil) and by the National Council for Scientific and Technological Development (CNPq: 403761/2020-4 and 400172/2022-4, 405786/2022-0, and INCT-One: 405786/2022-0). DP, ACA, FF, AG, HLF, and CWA are recipients of the CNPq scholarship. LMNS is the recipient of the CAPES scholarship. We also thank the National Network for Virus Surveillance in Wild Animals (PREVIR Network).

A collaborative research project under joint funding from FAPESP (The São Paulo Research Foundation) and NRF (National Research Foundation of Korea). The PMP-BS is funded by PETROBRAS (contract n. 5900.0111281.19.2).

Competing interest

The authors declare that they have no competing interests.

Authors’ Contribution

Conceptualization and data analysis: AYC, ACA, LMNS, ASA, EM, CKMK, RM, JVV, ASB, JLD, PVC CWA, D-HL, and HLF. Methodology, data curation, and writing—original draft preparation: AYC, D-HL, ACA, LMNS, and HLF. Sample preparation: ACA, AYC, LMNS, TCC, GCS, ASA, EM, CKMK, RM, JVV, GGL, ASB, JLD, PVC, DA, FM, MM, EGD, SP, SHL, CWA, ELD, JRRP. Writing—review and editing: AYC, D-HL, and HLF. Supervision and funding acquisition: D-HL, SHL, ELD, JRRP, CWA, and HLF. All authors have read and agreed to the published version of the manuscript.

Acknowledgements

The authors gratefully acknowledge Julia Cristina Benassi Bueno, and Alessandra Greatti for their technical assistance. We thank the support of the Brazilian Post and Telegraph Company (Correios) for transporting our fieldwork supplies and samples.

Samples were collected under the scope of The Santos Basin Beach Monitoring Project (PMP/BS) conducted by Petrobrás, licensed by the Brazilian Institute of the Environment and Renewable Natural Resources (IBAMA) of the Brazilian Ministry of Environment under ABIO Nº 640/2015.

The authors acknowledge all data contributions, i.e, the authors and their originating laboratories responsible for obtaining the specimens, and their submitting laboratories for generating the genetic sequences and metadata and sharing via the GISAID Initiative, on which this research is based.

Wille M, Barr IG: Resurgence of avian influenza virus. Science 2022, 376(6592):459-460.
Authority EFS, Prevention ECfD, Control, Influenza EURLfA, Adlhoch C, Fusaro A, Gonzales JL, Kuiken T, Mirinavičiūtė G, Niqueux É et al: Avian influenza overview June–September 2023. EFSA Journal 2023, 21(10):e08328.
Leguia M, Garcia-Glaessner A, Muñoz-Saavedra B, Juarez D, Barrera P, Calvo-Mac C, Jara J, Silva W, Ploog K, Amaro L et al: Highly pathogenic avian influenza A (H5N1) in marine mammals and seabirds in Peru. Nature Communications 2023, 14(1):5489.
Ariyama N, Pardo-Roa C, Muñoz G, Aguayo C, Ávila C, Mathieu C, Almonacid L, Medina R, Brito B, Johow M et al: Highly Pathogenic Avian Influenza A(H5N1) Clade 2.3.4.4b Virus in Wild Birds, Chile. Emerging Infectious Disease journal 2023, 29(9).
Gamarra-Toledo V, Plaza P, Gutiérrez R, Inga-Diaz G, Saravia-Guevara P, Pereyra-Meza O, Coronado-Flores E, Calderón-Cerrón A, Quiroz-Jiménez G, Martinez P et al: Mass Mortality of Sea Lions Caused by Highly Pathogenic Avian Influenza A(H5N1) Virus. Emerging Infectious Disease journal 2023, 29(12):2553.
Spackman E, Senne DA, Bulaga LL, Myers TJ, Perdue ML, Garber LP, Lohman K, Daum LT, Suarez DL: Development of real-time RT-PCR for the detection of avian influenza virus. Avian Diseases 2003, 47(3 Suppl):1079-1082.
Lee DH: Complete Genome Sequencing of Influenza A Viruses Using Next-Generation Sequencing. Methods Mol Biol 2020, 2123:69-79.
Letunic I, Bork P: Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Research 2021, 49(W1):W293-W296.
Youk S, Torchetti MK, Lantz K, Lenoch JB, Killian ML, Leyson C, Bevins SN, Dilione K, Ip HS, Stallknecht DE et al: H5N1 highly pathogenic avian influenza clade 2.3.4.4b in wild and domestic birds: Introductions into the United States and reassortments, December 2021–April 2022. Virology 2023, 587:109860.
Suttie A, Karlsson EA, Deng Y-M, Horm SV, Yann S, Tok S, Sorn S, Holl D, Tum S, Hurt AC et al: Influenza A(H5N1) viruses with A(H9N2) single gene (matrix or PB1) reassortment isolated from Cambodian live bird markets. Virology 2018, 523:22-26.
Suchard MA, Lemey P, Baele G, Ayres DL, Drummond AJ, Rambaut A: Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol 2018, 4(1):vey016.
Lee DH, Torchetti MK, Hicks J, Killian ML, Bahl J, Pantin-Jackwood M, Swayne DE: Transmission Dynamics of Highly Pathogenic Avian Influenza Virus A(H5Nx) Clade 2.3.4.4, North America, 2014-2015. Emerg Infect Dis 2018, 24(10):1840-1848.
Jackson D, Hossain MJ, Hickman D, Perez DR, Lamb RA: A new influenza virus virulence determinant: the NS1 protein four C-terminal residues modulate pathogenicity. Proc Natl Acad Sci U S A 2008, 105(11):4381-4386.
Suttie A, Deng YM, Greenhill AR, Dussart P, Horwood PF, Karlsson EA: Inventory of molecular markers affecting biological characteristics of avian influenza A viruses. Virus Genes 2019, 55(6):739-768.

Table 1 is available in the Supplementary Files section.

No competing interests reported.

Download PDF

Editorial decision: Revision requested
01 Jan, 2024
Editor assigned by journal
01 Jan, 2024
Submission checks completed at journal
25 Dec, 2023
First submitted to journal
22 Dec, 2023

You are reading this latest preprint version

Mortality in Sea Lions is associated with the introduction of the H5N1 clade 2.3.4.4b virus in Brazil, October 2023: Whole genome sequencing and phylogenetic analysis

Status:

Version 1

Abstract

Figures

Introduction, Methods, and Results

Discussion

Declarations

References

Table

Additional Declarations

Supplementary Files

Status:

Version 1