PARNAS: Objectively Selecting the Most Representative Taxa on a Phylogeny

Eurasian 1C swine influenza A virus exhibits high pandemic risk traits

ABSTRACTRecent surveillance has identified an expansion of swine H1 1C influenza viruses in Eurasian swine. Since 2010, at least 21 spillover events of 1C virus into humans have been detected and three of these occurred from July to December of 2023. Pandemic risk assessment of H1 1C influenza virus revealed that individuals born after 1950 had limited cross-reactive antibodies, confirming that they are antigenically novel viruses. The 1C virus exhibited phenotypic signatures similar to the 2009 pandemic H1N1 virus, including human receptor preference, productive replication in human airway cells, and robust environmental stability. Efficient inter- and intraspecies airborne transmission using the swine and ferret models was observed, including efficient airborne transmission to ferrets with pre-existing human seasonal H1N1 immunity. Together our data suggest H1 1C influenza virus poses a relatively high pandemic risk.

Emergence and interstate spread of highly pathogenic avian influenza A(H5N1) in dairy cattle in the United States

Highly pathogenic avian influenza (HPAI) viruses cross species barriers and have the potential to cause pandemics. In North America, HPAI A(H5N1) viruses related to the goose/Guangdong 2.3.4.4b hemagglutinin phylogenetic clade have infected wild birds, poultry, and mammals. Our genomic analysis and epidemiological investigation showed that a reassortment event in wild bird populations preceded a single wild bird-to-cattle transmission episode. The movement of asymptomatic or presymptomatic cattle has likely played a role in the spread of HPAI within the United States dairy herd. Some molecular markers that may lead to changes in transmission efficiency and phenotype were detected at low frequencies. Continued transmission of H5N1 HPAI within dairy cattle increases the risk for infection and subsequent spread of the virus to human populations.

Genomic characterization of novel bat kobuviruses in Madagascar: implications for viral evolution and zoonotic risk

Kobuviruses (family Picornaviridae, genus Kobuvirus) are enteric viruses that infect a wide range of both human and animal hosts. Much of the evolutionary history of kobuviruses remains elusive, largely due to limited screening in wildlife. Bats have been implicated as major sources of virulent zoonoses, including coronaviruses, henipaviruses, and filoviruses, though much of the bat virome still remains uncharacterized. While most bat virus research has historically focused on immediately recognizable zoonotic clades (e.g. SARS-related coronaviruses), a handful of prior reports catalog kobuvirus infection in bats and posit the role of bats as potential progenitors of downstream kobuvirus evolution. As part of a multi-year study, we carried out metagenomic Next Generation Sequencing (mNGS) on fecal samples obtained from endemic, wild-caught Madagascar fruit bats to characterize potentially zoonotic viruses circulating within these populations. The wild bats of Madagascar represent diverse Asian and African phylogeographic histories, presenting a unique opportunity for viruses from disparate origins to mix, posing significant public health threats. Here, we report detection of kobuvirus RNA in Malagasy fruit bat (Eidolon dupreanum) feces and undertake phylogenetic characterization of one full genome kobuvirus sequence, which nests within the Aichivirus A clade - a kobuvirus clade known to infect a wide range of hosts including humans, rodents, canids, felids, birds, and bats. Given the propensity of kobuviruses for recombination and cross-species infection, further characterization of this clade is critical to accurate evaluation of future zoonotic threats.

Tiled Amplicon Sequencing Enables Culture-free Whole-Genome Sequencing of Pathogenic Bacteria From Clinical Specimens

Pathogen sequencing is an important tool for disease surveillance and demonstrated its high value during the COVID-19 pandemic. Viral sequencing during the pandemic allowed us to track disease spread, quickly identify new variants, and guide the development of vaccines. Tiled amplicon sequencing, in which a panel of primers is used for multiplex amplification of fragments across an entire genome, was the cornerstone of SARS-CoV-2 sequencing. The speed, reliability, and cost-effectiveness of this method led to its implementation in academic and public health laboratories across the world and adaptation to a broad range of viral pathogens. However, similar methods are not available for larger bacterial genomes, for which whole-genome sequencing typically requires in vitro culture. This increases costs, error rates and turnaround times. The need to culture poses particular problems for medically important bacteria such as Mycobacterium tuberculosis, which are slow to grow and challenging to culture. As a proof of concept, we developed two novel whole-genome amplicon panels for M. tuberculosis and Streptococcus pneumoniae. Applying our amplicon panels to clinical samples, we show the ability to classify pathogen subgroups and to reliably identify markers of drug resistance without culturing. Development of this work in clinical settings has the potential to dramatically reduce the time of diagnosis of drug resistance for multiple drugs in parallel, enabling earlier intervention for high priority pathogens.

2018-2019 human seasonal H3N2 influenza A virus spillovers into swine with demonstrated virus transmission in pigs were not sustained in the pig population

Human seasonal H3 clade 3C3a influenza A viruses (IAV) were detected four times in U.S. pigs from commercial swine farms in Michigan, Illinois, and Virginia in 2019. To evaluate the relative risk of this spillover to the pig population, whole genome sequencing and phylogenetic characterization were conducted, and the results revealed that all eight viral gene segments were closely related to 2018-2019 H3N2 human seasonal IAV. Next, a series of in vitro viral kinetics, receptor binding, and antigenic characterization studies were performed using a representative A/swine/Virginia/A02478738/2018(H3N2) (SW/VA/19) isolate. Viral replication kinetic studies of SW/VA/19 demonstrated less efficient replication curves than all 10 swine H3N2 viruses tested but higher than three human H3N2 strains. Serial passaging experiments of SW/VA/19 in swine cells did not increase virus replication, but changes at HA amino acid positions 9 and 159 occurred. In swine transmission studies, wild-type SW/VA/19 was shed in nasal secretions and transmitted to all indirect contact pigs, whereas the human seasonal strain A/Switzerland/9715293/2013(H3N2) from the same 3C3a clade failed to transmit. SW/VA/19 induced minimal macroscopic and microscopic lung lesions. Collectively, these findings demonstrate that these human seasonal H3N2 3C3a-like viruses did not require reassortment with endemic swine IAV gene segments for virus shedding and transmission in pigs. Limited detections in the U.S. pig population in the subsequent period of time suggest a yet-unknown restriction factor likely limiting the spread of these viruses in the U.S. pig population.IMPORTANCEInterspecies human-to-swine IAV transmission occurs globally and contributes to increased IAV diversity in pig populations. We present data that a swine isolate from a 2018-2019 human-to-swine transmission event was shed for multiple days in challenged and contact pigs. By characterizing this introduction through bioinformatic, molecular, and animal experimental approaches, these findings better inform animal health practices and vaccine decision-making. Since wholly human seasonal H3N2 viruses in the United States were not previously identified as being transmissible in pigs (i.e., reverse zoonosis), these findings reveal that the interspecies barriers for transmission to pigs may not require significant changes to all human seasonal H3N2, although additional changes may be required for sustained transmission in swine populations.

Progress and persistence of diseases of high consequence to livestock in the United States

The USDA/ARS-National Disease Center (NADC) will celebrate its 65th anniversary of existence in November 2026. NADC continues as one of the world's premier animal health research centers conducting basic and applied research on endemic diseases with economic impact on U.S. livestock and wildlife. This research center also supports a program studying important food safety pathogens such as Salmonella, E. coli and Campylobacter. NADC has contributed significantly to the elimination of a few diseases, notably hog cholera and milk fever, and made progress in reducing the impact of many other animal diseases through vaccines, therapies and managerial recommendations. Despite nearly 65 years of targeted research on these diseases and much progress, some of these continue to persist. The reasons for such persistence varies for each disease condition and they are often multifactorial involving host susceptibility, virulence and even environmental conditions. Individually and in aggregate, these disease conditions have a massive economic impact and can be devasting to animal producers, owners and individuals that become infected through zoonotic disease agents such as tuberculosis, leptospirosis and avian influenza. They also diminish the health, well-being and welfare of affected animals, which directly affects the food supply. The NADC is using all available technologies including genomic, biochemical, reverse genetics, and vaccine trials in the target host to combat these significant diseases. We review the progress and reasons for persistence of selected diseases and food safety pathogens as well as the progress and potential outcomes should research and programmatic plans to eliminate these disease conditions cease.

Influenza A virus shedding and reinfection during the post-weaning period in swine: longitudinal study of two nurseries

INTRODUCTION

Influenza A virus in swine (IAV-S) is common in the United States commercial swine population and has the potential for zoonotic transmission.

OBJECTIVE

To elucidate influenza shedding the domestic pig population, we evaluated two commercial swine farms in Illinois, United States, for 7 weeks. Farm 1 had a recent IAV-S outbreak. Farm 2 has had IAV-S circulating for several years.

METHODS

Forty post-weaning pigs on Farm 1 and 51 pigs from Farm 2 were individually monitored and sampled by nasal swabs for 7 weeks.

RESULTS

RT-PCR results over time showed most piglets shed in the first 2 weeks post weaning, with 91.2% shedding in week one, and 36.3% in week two. No difference in the number of pigs shedding was found between the two nurseries. Reinfection events did differ between the farms, with 30% of piglets on Farm 1 becoming reinfected, compared to 7.8% on Farm 2. In addition, whole genome sequencing of nasal swab samples from each farm showed identical viruses circulating between the initial infection and the reinfection periods. Sequencing also allowed for nucleic and amino acid mutation analysis in the circulating viruses, as well the identification of a potential reverse zoonosis event. We saw antigenic site mutations arising in some pigs and MxA resistance genes in almost all samples.

CONCLUSION

This study provided information on IAV-S circulation in nurseries to aid producers and veterinarians to screen appropriately for IAV-S, determine the duration of IAV-S shedding, and predict the occurrence of reinfection in the nursery period.

Immune History Modifies Disease Severity to HPAI H5N1 Clade 2.3.4.4b Viral Challenge

The most recent outbreak of highly pathogenic avian H5 influenza (HPAI) virus in cattle is now widespread across the U.S. with spillover events happening to other mammals, including humans. Several human cases have been reported with clinical signs ranging from conjunctivitis to respiratory illness. However, most of those infected report mild to moderate symptoms, while previously reported HPAI H5Nx infections in humans have had mortality rates upwards of 50%. We recently reported that mice with pre-existing immunity to A/Puerto Rico/08/1934 H1N1 virus were protected from lethal challenge from highly pathogenic clade 2.3.4.4b H5N1 influenza virus. Here, we demonstrate that mice infected with the 2009 pandemic H1N1 virus strain A/California/04/2009 (Cal09) or vaccinated with a live-attenuated influenza vaccine (LAIV) were moderately-to-highly protected against a lethal A/bovine/Ohio/B24OSU-439/2024 H5N1 virus challenge. We also observed that ferrets with mixed pre-existing immunity-either from LAIV vaccination and/or from Cal09 infection-showed protection against a HPAI H5N1 clade 2.3.4.4b virus isolated from a cat. Notably, this protection occurred independently of any detectable hemagglutination inhibition titers (HAIs) against the H5N1 virus. To explore factors that may contribute to protection, we conducted detailed T cell epitope mapping using previously published sequences from H1N1 strains. This analysis revealed a high conservation of amino acid sequences within the internal proteins of our bovine HPAI H5N1 virus strain. These data highlight the necessity to explore additional factors that contribute to protection against HPAI H5N1 viruses, such as memory T cell responses, in addition to HA-inhibition or neutralizing antibodies.

First detection of influenza A virus subtypes H1N1 and H3N8 in the Antarctic region: King George Island, 2023

RELEVANCE

Influenza A virus is characterized by a segmented single-stranded RNA genome. Such organization of the virus genome determines the possibility of reassortment, which can lead to the emergence of new virus variants. The main natural reservoir of most influenza A virus subtypes are wild waterfowl. Seasonal migrations gather waterfowl from all major migration routes to nesting areas near the northern and southern polar circles. This makes intercontinental spread of influenza A viruses possible. Objective ‒ to conduct molecular genetic monitoring and study the phylogenetic relationships of influenza A virus variants circulating in Antarctica in 2023.

MATERIALS AND METHODS

We studied 84 samples of biological material obtained from birds and marine mammals in April‒May 2023 in coastal areas of Antarctica. For 3 samples, sequencing was performed on the Miseq, Illumina platform and phylogenetic analysis of the obtained nucleotide sequences of the influenza A virus genomes was performed.

RESULTS

The circulation of avian influenza virus in the Antarctic region was confirmed. Heterogeneity of the pool of circulating variants of the influenza A virus (H3N8, H1N1) was revealed. Full-length genomes of the avian influenza virus were sequenced and posted in the GISAID database (EPI_ISL_19032103, 19174530, 19174467).

CONCLUSION

The study of the genetic diversity of influenza A viruses circulating in the polar regions of the Earth and the identification of the conditions for the emergence of new genetic variants is a relevant task for the development of measures to prevent biological threats.

Detection and spread of high pathogenicity avian influenza virus H5N1 in the Antarctic Region

Until recent events, the Antarctic was the only major geographical region in which high pathogenicity avian influenza virus (HPAIV) had never previously been detected. Here we report on the detection of clade 2.3.4.4b H5N1 HPAIV in the Antarctic and sub-Antarctic regions of South Georgia and the Falkland Islands, respectively. We initially detected H5N1 HPAIV in samples collected from brown skuas at Bird Island, South Georgia on 8th October 2023. Since this detection, mortalities were observed in several avian and mammalian species at multiple sites across South Georgia. Subsequent testing confirmed H5N1 HPAIV across several sampling locations in multiple avian species and two seal species. Simultaneously, we also confirmed H5N1 HPAIV in southern fulmar and black-browed albatross in the Falkland Islands. Genetic assessment of the virus indicates spread from South America, likely through movement of migratory birds. Critically, genetic assessment of sequences from mammalian species demonstrates no increased risk to human populations above that observed in other instances of mammalian infections globally. Here we describe the detection, species impact and genetic composition of the virus and propose both introductory routes and potential long-term impact on avian and mammalian species across the Antarctic region. We also speculate on the threat to specific populations following recent reports in the area.

Immunogenicity of an Inactivated Senecavirus A Vaccine with a Contemporary Brazilian Strain in Mice

Senecavirus A (SVA) is a picornavirus that is endemic in swine, causing a vesicular disease clinically indistinguishable from other vesicular diseases, like foot-and-mouth disease. The widespread viral circulation, constant evolution, and economic losses caused to the swine industry emphasize the need for measures to control the agent. In this study, we evaluated the immunogenicity of a whole-virus-inactivated vaccine using a representative contemporary Brazilian SVA strain in Balb/ByJ mice. The animals were vaccinated with two doses by an intramuscular route. The humoral response induced by the vaccination was evaluated by an in-house ELISA assay for IgG detection. The cellular response was assessed by flow cytometry after in vitro SVA stimulation in splenocyte cultures from vaccinated and non-vaccinated groups. Protection against SVA was assessed in the experimental groups following an oral challenge with the homologous virus. The vaccination induced high levels of IgG antibodies and the proliferation of CD45R/B220+sIgM+, CD3e+CD69+, and CD3e+CD4+CD44+CD62L- cells. These results indicate the immunogenicity and safety of the vaccine formulation in a murine model and the induction of humoral and cellular response against SVA.

Evaluation of the Increased Genetic Resolution and Utility for Source Tracking of a Recently Developed Method for Genotyping Cyclospora cayetanensis

Cyclospora cayetanensis is a foodborne parasite that causes cyclosporiasis, an enteric illness in humans. Genotyping methods are used to genetically discriminate between specimens from cyclosporiasis cases and can complement source attribution investigations if the method is sufficiently sensitive for application to food items. A very sensitive targeted amplicon sequencing (TAS) assay for genotyping C. cayetanensis encompassing 52 loci was recently designed. In this study, we analyzed 66 genetically diverse clinical specimens to assess the change in phylogenetic resolution between the TAS assay and a currently employed eight-marker scheme. Of the 52 markers, ≥50 were successfully haplotyped for all specimens, and these results were used to generate a hierarchical cluster dendrogram. Using a previously described statistical approach to dissect hierarchical trees, the 66 specimens resolved into 24 and 27 distinct genetic clusters for the TAS and an 8-loci scheme, respectively. Although the specimen composition of 15 clusters was identical, there were substantial differences between the two dendrograms, highlighting the importance of both inclusion of additional genome coverage and choice of loci to target for genotyping. To evaluate the ability to genetically link contaminated food samples with clinical specimens, C. cayetanensis was genotyped from DNA extracted from raspberries inoculated with fecal specimens. The contaminated raspberry samples were assigned to clusters with the corresponding clinical specimen, demonstrating the utility of the TAS assay for traceback efforts.

A Guide to Phylogenomic Inference

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

An improved framework for detecting discrete epidemiologically meaningful partitions in hierarchically clustered genetic data

Motivation

Hierarchical clustering of microbial genotypes has the limitation that hierarchical clusters are nested, where smaller groups of related isolates exist within larger groups that get progressively larger as relationships become increasingly distant. In an epidemiologic context, investigators must dissect hierarchical trees into discrete groupings that are epidemiologically meaningful. We recently described a statistical framework (Method A) for dissecting hierarchical trees that attempts to minimize investigator bias. Here, we apply a modified version of that framework (Method B) to a hierarchical tree constructed from 2111 genotypes of the foodborne parasite Cyclospora, including 639 genotypes linked to epidemiologically defined outbreaks. To evaluate Method B's performance, we examined the concordance between these epidemiologically defined groupings and the genetic partitions identified. We also used the same epidemiologic clusters to evaluate the performance of Method A, plus two tree-dissection methods (cutreeHybrid and cutreeDynamic) available within the Dynamic Tree Cut R package, in addition to the TreeCluster method and PARNAS.

Results

Compared to the other methods, Method B, TreeCluster, and PARNAS were the most accurate (99.4%) in identifying genetic groups that reflected the epidemiologic groupings, noting that TreeCluster and PARNAS performed identically on our dataset. CutreeHybrid identified groups reflecting patterns in the wider Cyclospora population structure but lacked finer, strain-level discrimination (Simpson's D: cutreeHybrid=0.785). CutreeDynamic displayed good strain discrimination (Simpson's D = 0.933), though lacked sensitivity (77%). At two different threshold/radius settings TreeCluster/PARNAS displayed similar utility to Method B. However, Method B computes a tree-dissection threshold automatically, and the threshold/radius settings used when executing TreeCluster/PARNAS here were computed using Method B. Using a TreeCluster threshold of 0.045 as recommended in the TreeCluster documentation, epidemiologic utility dropped markedly below that of Method B.

Availability and implementation

Relevant code and data are publicly available. Source code (Method B) and instructions for its use are available here: https://github.com/Joel-Barratt/Hierarchical-tree-dissection-framework.