Phylogenetic analysis, molecular changes, and adaptation to chickens of Mexican lineage H5N2 low-pathogenic avian influenza viruses from 1994 to 2019.

The Mexican lineage H5N2 low pathogenic avian influenza viruses (LPAIVs) were first detected in 1994 and mutated to highly pathogenic avian influenza viruses (HPAIVs) in 1994-1995 causing widespread outbreaks in poultry. By using vaccination and other control measures, the HPAIVs were eradicated but the LPAIVs continued circulating in Mexico and spread to several other countries. To get better resolution of the phylogenetics of this virus, the full genome sequences of 44 H5N2 LPAIVs isolated from 1994 to 2011, and 6 detected in 2017 and 2019, were analysed. Phylogenetic incongruence demonstrated genetic reassortment between two separate groups of the Mexican lineage H5N2 viruses between 2005 and 2010. Moreover, the recent H5N2 viruses reassorted with previously unidentified avian influenza viruses. Bayesian phylogeographic results suggested that mechanical transmission involving human activity is the most probable cause of the virus spillover to Central American, Caribbean, and East Asian countries. Increased infectivity and transmission of a 2011 H5N2 LPAIV in chickens compared to a 1994 virus demonstrates improved adaptation to chickens, while low virus shedding, and limited contact transmission was observed in mallards with the same 2011 virus. The sporadic increase in basic amino acids in the HA cleavage site, changes in potential N-glycosylation sites in the HA, and truncations of PB1-F2 should be further examined in relation to the increased infectivity and transmission in poultry. The genetic changes that occur as this lineage of H5N2 LPAIVs continues circulating in poultry is concerning not only because of the effect of these changes on vaccination efficacy, but also because of the potential of the viruses to mutate to the highly pathogenic form. Continued vigilance and surveillance efforts, and the pathogenic and genetic characterization of circulating viruses, are required for the effective control of this virus.

The Pathobiology of H7N3 Low and High Pathogenicity Avian Influenza Viruses from the United States Outbreak in 2020 Differs between Turkeys and Chickens.

An outbreak caused by H7N3 low pathogenicity avian influenza virus (LPAIV) occurred in commercial turkey farms in the states of North Carolina (NC) and South Carolina (SC), United States in March of 2020. Subsequently, H7N3 high pathogenicity avian influenza virus (HPAIV) was detected on a turkey farm in SC. The infectivity, transmissibility, and pathogenicity of the H7N3 HPAIV and two LPAIV isolates, including one with a deletion in the neuraminidase (NA) protein stalk, were studied in turkeys and chickens. High infectivity [<2 log10 50% bird infectious dose (BID50)] and transmission to birds exposed by direct contact were observed with the HPAIV in turkeys. In contrast, the HPAIV dose to infect chickens was higher than for turkeys (3.7 log10 BID50), and no transmission was observed. Similarly, higher infectivity (<2-2.5 log10 BID50) and transmissibility were observed with the H7N3 LPAIVs in turkeys compared to chickens, which required higher virus doses to become infected (5.4-5.7 log10 BID50). The LPAIV with the NA stalk deletion was more infectious in turkeys but did not have enhanced infectivity in chickens. These results show clear differences in the pathobiology of AIVs in turkeys and chickens and corroborate the high susceptibility of turkeys to both LPAIV and HPAIV infections.

Surveillance and Genetic Characterization of Virulent Newcastle Disease Virus Subgenotype V.3 in Indigenous Chickens from Backyard Poultry Farms and Live Bird Markets in Kenya.

Kenyan poultry consists of ~80% free-range indigenous chickens kept in small flocks (~30 birds) on backyard poultry farms (BPFs) and they are traded via live bird markets (LBMs). Newcastle disease virus (NDV) was detected in samples collected from chickens, wild farm birds, and other domestic poultry species during a 2017-2018 survey conducted at 66 BPFs and 21 LBMs in nine Kenyan counties. NDV nucleic acids were detected by rRT-PCR L-test in 39.5% (641/1621) of 1621 analyzed samples, of which 9.67% (62/641) were NDV-positive by both the L-test and a fusion-test designed to identify the virulent virus, with a majority being at LBMs (64.5%; 40/62) compared to BPFs (25.5%; 22/62). Virus isolation and next-generation sequencing (NGS) on a subset of samples resulted in 32 complete NDV genome sequences with 95.8-100% nucleotide identities amongst themselves and 95.7-98.2% identity with other east African isolates from 2010-2016. These isolates were classified as a new sub-genotype, V.3, and shared 86.5-88.9% and 88.5-91.8% nucleotide identities with subgenotypes V.1 and V.2 viruses, respectively. The putative fusion protein cleavage site (113R-Q-K-R↓F 117) in all 32 isolates, and a 1.86 ICPI score of an isolate from a BPF chicken that had clinical signs consistent with Newcastle disease, confirmed the high virulence of the NDVs. Compared to genotypes V and VI viruses, the attachment (HN) protein of 18 of the 32 vNDVs had amino acid substitutions in the antigenic sites. A time-scaled phylogeographic analysis suggests a west-to-east dispersal of the NDVs via the live chicken trade, but the virus origins remain unconfirmed due to scarcity of continuous and systematic surveillance data. This study reveals the widespread prevalence of vNDVs in Kenyan backyard poultry, the central role of LBMs in the dispersal and possibly generation of new virus variants, and the need for robust molecular epidemiological surveillance in poultry and non-poultry avian species.

Whole gut microbiome composition of damselfish and cardinalfish before and after reef settlement.

The Pomacentridae (damselfish) and Apogonidae (cardinalfish) are among the most common fish families on coral reefs and in the aquarium trade. Members of both families undergo a pelagic larvae phase prior to settlement on the reef, where adults play key roles in benthic habitat structuring and trophic interactions. Fish-associated microbial communities (microbiomes) significantly influence fish health and ecology, yet little is known of how microbiomes change with life stage. We quantified the taxonomic (16S rRNA gene) composition of whole gut microbiomes from ten species of damselfish and two species of cardinalfish from Lizard Island, Australia, focusing specifically on comparisons between pelagic larvae prior to settlement on the reef versus post-settlement juvenile and adult individuals. On average, microbiome phylogenetic diversity increased from pre- to post-settlement, and was unrelated to the microbial composition in the surrounding water column. However, this trend varied among species, suggesting stochasticity in fish microbiome assembly. Pre-settlement fish were enriched with bacteria of the Endozoicomonaceae, Shewanellaceae, and Fusobacteriaceae, whereas settled fish harbored higher abundances of Vibrionaceae and Pasteurellaceae. Several individual operational taxonomic units, including ones related to Vibrio harveyi, Shewanella sp., and uncultured Endozoicomonas bacteria, were shared between both pre and post-settlement stages and may be of central importance in the intestinal niche across development. Richness of the core microbiome shared among pre-settlement fish was comparable to that of settled individuals, suggesting that changes in diversity with adulthood are due to the acquisition or loss of host-specific microbes. These results identify a key transition in microbiome structure across host life stage, suggesting changes in the functional contribution of microbiomes over development in two ecologically dominant reef fish families.

Standard filtration practices may significantly distort planktonic microbial diversity estimates.

Fractionation of biomass by filtration is a standard method for sampling planktonic microbes. It is unclear how the taxonomic composition of filtered biomass changes depending on sample volume. Using seawater from a marine oxygen minimum zone, we quantified the 16S rRNA gene composition of biomass on a prefilter (1.6 µm pore-size) and a downstream 0.2 µm filter over sample volumes from 0.05 to 5 L. Significant community shifts occurred in both filter fractions, and were most dramatic in the prefilter community. Sequences matching Vibrionales decreased from ~40 to 60% of prefilter datasets at low volumes (0.05-0.5 L) to less than 5% at higher volumes, while groups such at the Chromatiales and Thiohalorhabdales followed opposite trends, increasing from minor representation to become the dominant taxa at higher volumes. Groups often associated with marine particles, including members of the Deltaproteobacteria, Planctomycetes, and Bacteroidetes, were among those showing the greatest increase with volume (4 to 27-fold). Taxon richness (97% similarity clusters) also varied significantly with volume, and in opposing directions depending on filter fraction, highlighting potential biases in community complexity estimates. These data raise concerns for studies using filter fractionation for quantitative comparisons of aquatic microbial diversity, for example between free-living and particle-associated communities.

Microbial eukaryote diversity in the marine oxygen minimum zone off northern Chile.

Molecular surveys are revealing diverse eukaryotic assemblages in oxygen-limited ocean waters. These communities may play pivotal ecological roles through autotrophy, feeding, and a wide range of symbiotic associations with prokaryotes. We used 18S rRNA gene sequencing to provide the first snapshot of pelagic microeukaryotic community structure in two cellular size fractions (0.2-1.6 µm, >1.6 µm) from seven depths through the anoxic oxygen minimum zone (OMZ) off northern Chile. Sequencing of >154,000 amplicons revealed contrasting patterns of phylogenetic diversity across size fractions and depths. Protist and total eukaryote diversity in the >1.6 µm fraction peaked at the chlorophyll maximum in the upper photic zone before declining by ~50% in the OMZ. In contrast, diversity in the 0.2-1.6 µm fraction, though also elevated in the upper photic zone, increased four-fold from the lower oxycline to a maximum at the anoxic OMZ core. Dinoflagellates of the Dinophyceae and endosymbiotic Syndiniales clades dominated the protist assemblage at all depths (~40-70% of sequences). Other protist groups varied with depth, with the anoxic zone community of the larger size fraction enriched in euglenozoan flagellates and acantharean radiolarians (up to 18 and 40% of all sequences, respectively). The OMZ 0.2-1.6 µm fraction was dominated (11-99%) by Syndiniales, which exhibited depth-specific variation in composition and total richness despite uniform oxygen conditions. Metazoan sequences, though confined primarily to the 1.6 µm fraction above the OMZ, were also detected within the anoxic zone where groups such as copepods increased in abundance relative to the oxycline and upper OMZ. These data, compared to those from other low-oxygen sites, reveal variation in OMZ microeukaryote composition, helping to identify clades with potential adaptations to oxygen-depletion.

Metagenomic analysis of size-fractionated picoplankton in a marine oxygen minimum zone.

Marine oxygen minimum zones (OMZs) support diverse microbial communities with roles in major elemental cycles. It is unclear how the taxonomic composition and metabolism of OMZ microorganisms vary between particle-associated and free-living size fractions. We used amplicon (16S rRNA gene) and shotgun metagenome sequencing to compare microbial communities from large (>1.6 µm) and small (0.2-1.6 µm) filter size fractions along a depth gradient in the OMZ off Chile. Despite steep vertical redox gradients, size fraction was a significantly stronger predictor of community composition compared to depth. Phylogenetic diversity showed contrasting patterns, decreasing towards the anoxic OMZ core in the small size fraction, but exhibiting maximal values at these depths within the larger size fraction. Fraction-specific distributions were evident for key OMZ taxa, including anammox planctomycetes, whose coding sequences were enriched up to threefold in the 0.2-1.6 µm community. Functional gene composition also differed between fractions, with the >1.6 µm community significantly enriched in genes mediating social interactions, including motility, adhesion, cell-to-cell transfer, antibiotic resistance and mobile element activity. Prokaryotic transposase genes were three to six fold more abundant in this fraction, comprising up to 2% of protein-coding sequences, suggesting that particle surfaces may act as hotbeds for transposition-based genome changes in marine microbes. Genes for nitric and nitrous oxide reduction were also more abundant (three to seven fold) in the larger size fraction, suggesting microniche partitioning of key denitrification steps. These results highlight an important role for surface attachment in shaping community metabolic potential and genome content in OMZ microorganisms.