Risk assessment of a highly pathogenic H5N1 influenza virus from mink.

Outbreaks of highly pathogenic H5N1 clade 2.3.4.4b viruses in farmed mink and seals combined with isolated human infections suggest these viruses pose a pandemic threat. To assess this threat, using the ferret model, we show an H5N1 isolate derived from mink transmits by direct contact to 75% of exposed ferrets and, in airborne transmission studies, the virus transmits to 37.5% of contacts. Sequence analyses show no mutations were associated with transmission. The H5N1 virus also has a low infectious dose and remains virulent at low doses. This isolate carries the adaptive mutation, PB2 T271A, and reversing this mutation reduces mortality and airborne transmission. This is the first report of a H5N1 clade 2.3.4.4b virus exhibiting direct contact and airborne transmissibility in ferrets. These data indicate heightened pandemic potential of the panzootic H5N1 viruses and emphasize the need for continued efforts to control outbreaks and monitor viral evolution.

Genetic drift and purifying selection shape within-host influenza A virus populations during natural swine infections.

Patterns of within-host influenza A virus (IAV) diversity and evolution have been described in natural human infections, but these patterns remain poorly characterized in non-human hosts. Elucidating these dynamics is important to better understand IAV biology and the evolutionary processes that govern spillover into humans. Here, we sampled an IAV outbreak in pigs during a week-long county fair to characterize viral diversity and evolution in this important reservoir host. Nasal wipes were collected on a daily basis from all pigs present at the fair, yielding up to 421 samples per day. Subtyping of PCR-positive samples revealed the co-circulation of H1N1 and H3N2 subtype swine IAVs. PCR-positive samples with robust Ct values were deep-sequenced, yielding 506 sequenced samples from a total of 253 pigs. Based on higher-depth re-sequenced data from a subset of these initially sequenced samples (260 samples from 168 pigs), we characterized patterns of within-host IAV genetic diversity and evolution. We find that IAV genetic diversity in single-subtype infected pigs is low, with the majority of intrahost Single Nucleotide Variants (iSNVs) present at frequencies of <10%. The ratio of the number of nonsynonymous to the number of synonymous iSNVs is significantly lower than under the neutral expectation, indicating that purifying selection shapes patterns of within-host viral diversity in swine. The dynamic turnover of iSNVs and their pronounced frequency changes further indicate that genetic drift also plays an important role in shaping IAV populations within pigs. Taken together, our results highlight similarities in patterns of IAV genetic diversity and evolution between humans and swine, including the role of stochastic processes in shaping within-host IAV dynamics.

Genetic drift and purifying selection shape within-host influenza A virus populations during natural swine infections.

Patterns of within-host influenza A virus (IAV) diversity and evolution have been described in natural human infections, but these patterns remain poorly characterized in non-human hosts. Elucidating these dynamics is important to better understand IAV biology and the evolutionary processes that govern spillover into humans. Here, we sampled an IAV outbreak in pigs during a week-long county fair to characterize viral diversity and evolution in this important reservoir host. Nasal wipes were collected on a daily basis from all pigs present at the fair, yielding up to 421 samples per day. Subtyping of PCR-positive samples revealed the co-circulation of H1N1 and H3N2 subtype IAVs. PCR-positive samples with robust Ct values were deep-sequenced, yielding 506 sequenced samples from a total of 253 pigs. Based on higher-depth re-sequenced data from a subset of these initially sequenced samples (260 samples from 168 pigs), we characterized patterns of within-host IAV genetic diversity and evolution. We find that IAV genetic diversity in single-subtype infected pigs is low, with the majority of intra-host single nucleotide variants (iSNVs) present at frequencies of <10%. The ratio of the number of nonsynonymous to the number of synonymous iSNVs is significantly lower than under the neutral expectation, indicating that purifying selection shapes patterns of within-host viral diversity in swine. The dynamic turnover of iSNVs and their pronounced frequency changes further indicate that genetic drift also plays an important role in shaping IAV populations within pigs. Taken together, our results highlight similarities in patterns of IAV genetic diversity and evolution between humans and swine, including the role of stochastic processes in shaping within-host IAV dynamics.

Viral expansion after transfer is a primary driver of influenza A virus transmission bottlenecks.

For many viruses, narrow bottlenecks acting during transmission sharply reduce genetic diversity in a recipient host relative to the donor. Since genetic diversity represents adaptive potential, such losses of diversity are though to limit the opportunity for viral populations to undergo antigenic change and other adaptive processes. Thus, a detailed picture of evolutionary dynamics during transmission is critical to understanding the forces driving viral evolution at an epidemiologic scale. To advance this understanding, we used a novel barcoded virus library and a guinea pig model of transmission to decipher where in the transmission process diversity is lost for influenza A viruses. In inoculated guinea pigs, we show that a high level of viral genetic diversity is maintained across time. Continuity in the barcodes detected furthermore indicates that stochastic effects are not pronounced within inoculated hosts. Importantly, in both aerosol-exposed and direct contact-exposed animals, we observed many barcodes at the earliest time point(s) positive for infectious virus, indicating robust transfer of diversity through the environment. This high viral diversity is short-lived, however, with a sharp decline seen 1-2 days after initiation of infection. Although major losses of diversity at transmission are well described for influenza A virus, our data indicate that events that occur following viral transfer and during the earliest stages of natural infection have a predominant role in this process. This finding suggests that immune selection may have greater opportunity to operate during influenza A transmission than previously recognized.

Intestinal barrier dysfunction in murine sickle cell disease is associated with small intestine neutrophilic inflammation, oxidative stress, and dysbiosis.

The intestinal microbiome has emerged as a potential contributor to the severity of sickle cell disease (SCD). We sought to determine whether SCD mice exhibit intestinal barrier dysfunction, inflammation, and dysbiosis. Using the Townes humanized sickle cell mouse model, we found a 3-fold increase in intestinal permeability as assessed via FITC-dextran (4 kDa) assay in SS (SCD) mice compared to AA (wild type) mice (n = 4, p < 0.05). This was associated with 25 to 50% decreases in claudin-1, 3, and 15 and zonula occludens-1 gene expression (n = 8-10, p < 0.05) in the small intestine. Increased Ly6G staining demonstrated more neutrophils in the SS small intestine (3-fold, n = 5, p < 0.05) associated with increased expression of TNFα, IL-17A, CXCL1, and CD68 (2.5 to 5-fold, n = 7-10, p < 0.05). In addition, we observed 30 to 55% decreases in superoxide dismutase-1, glutathione peroxidase-1, and catalase antioxidant enzyme expression (n = 7-8, p < 0.05) concomitant to an increase in superoxide (2-fold, n = 4, p < 0.05). Importantly, all significant observations of a leaky gut phenotype and inflammation were limited to the small intestine and not observed in the colon. Finally, characterization of the composition of the microbiome within the small intestine revealed dysbiosis in SS mice compared to their AA littermates with 47 phyla to species-level significant alterations in amplicon sequence variants. We conclude that the intestinal barrier is compromised in SCD, associated with decreased gene expression of tight junction proteins, enhanced inflammation, oxidative stress, and gut microbiome dysbiosis, all specific to the small intestine.

Host Association and Spatial Proximity Shape but Do Not Constrain Population Structure in the Mutualistic Symbiont Xenorhabdus bovienii.

To what extent are generalist species cohesive evolutionary units rather than a compilation of recently diverged lineages? We examine this question in the context of host specificity and geographic structure in the insect pathogen and nematode mutualist Xenorhabdus bovienii. This bacterial species partners with multiple nematode species across two clades in the genus Steinernema. We sequenced the genomes of 42 X. bovienii strains isolated from four different nematode species and three field sites within a 240-km2 region and compared them to globally available reference genomes. We hypothesized that X. bovienii would comprise several host-specific lineages, such that bacterial and nematode phylogenies would be largely congruent. Alternatively, we hypothesized that spatial proximity might be a dominant signal, as increasing geographic distance might lower shared selective pressures and opportunities for gene flow. We found partial support for both hypotheses. Isolates clustered largely by nematode host species but did not strictly match the nematode phylogeny, indicating that shifts in symbiont associations across nematode species and clades have occurred. Furthermore, both genetic similarity and gene flow decreased with geographic distance across nematode species, suggesting differentiation and constraints on gene flow across both factors, although no absolute barriers to gene flow were observed across the regional isolates. Several genes associated with biotic interactions were found to be undergoing selective sweeps within this regional population. The interactions included several insect toxins and genes implicated in microbial competition. Thus, gene flow maintains cohesiveness across host associations in this symbiont and may facilitate adaptive responses to a multipartite selective environment. IMPORTANCE Microbial populations and species are notoriously hard to delineate. We used a population genomics approach to examine the population structure and the spatial scale of gene flow in Xenorhabdus bovienii, an intriguing species that is both a specialized mutualistic symbiont of nematodes and a broadly virulent insect pathogen. We found a strong signature of nematode host association, as well as evidence for gene flow connecting isolates associated with different nematode host species and collected from distinct study sites. Furthermore, we saw signatures of selective sweeps for genes involved with nematode host associations, insect pathogenicity, and microbial competition. Thus, X. bovienii exemplifies the growing consensus that recombination not only maintains cohesion but can also allow the spread of niche-beneficial alleles.

A conserved Bacteroidetes antigen induces anti-inflammatory intestinal T lymphocytes.

The microbiome contributes to the development and maturation of the immune system. In response to commensal bacteria, intestinal CD4+ T lymphocytes differentiate into functional subtypes with regulatory or effector functions. The development of small intestine intraepithelial lymphocytes that coexpress CD4 and CD8αα homodimers (CD4IELs) depends on the microbiota. However, the identity of the microbial antigens recognized by CD4+ T cells that can differentiate into CD4IELs remains unknown. We identified ß-hexosaminidase, a conserved enzyme across commensals of the Bacteroidetes phylum, as a driver of CD4IEL differentiation. In a mouse model of colitis, ß-hexosaminidase-specific lymphocytes protected against intestinal inflammation. Thus, T cells of a single specificity can recognize a variety of abundant commensals and elicit a regulatory immune response at the intestinal mucosa.

Microbial metabolite delta-valerobetaine is a diet-dependent obesogen.

Obesity and obesity-related metabolic disorders are linked to the intestinal microbiome. However, the causality of changes in the microbiome-host interaction affecting energy metabolism remains controversial. Here, we show the microbiome-derived metabolite δ-valerobetaine (VB) is a diet-dependent obesogen that is increased with phenotypic obesity and is correlated with visceral adipose tissue mass in humans. VB is absent in germ-free mice and their mitochondria but present in ex-germ-free conventionalized mice and their mitochondria. Mechanistic studies in vivo and in vitro show VB is produced by diverse bacterial species and inhibits mitochondrial fatty acid oxidation through decreasing cellular carnitine and mitochondrial long-chain acyl-coenzyme As. VB administration to germ-free and conventional mice increases visceral fat mass and exacerbates hepatic steatosis with a western diet but not control diet. Thus, VB provides a molecular target to understand and potentially manage microbiome-host symbiosis or dysbiosis in diet-dependent obesity.

Rapid evolutionary turnover of mobile genetic elements drives bacterial resistance to phages.

Although it is generally accepted that phages drive bacterial evolution, how these dynamics play out in the wild remains poorly understood. We found that susceptibility to viral killing in marine Vibrio is mediated by large and highly diverse mobile genetic elements. These phage defense elements display exceedingly fast evolutionary turnover, resulting in differential phage susceptibility among clonal bacterial strains while phage receptors remain invariant. Protection is cumulative, and a single bacterial genome can harbor 6 to 12 defense elements, accounting for more than 90% of the flexible genome among close relatives. The rapid turnover of these elements decouples phage resistance from other genomic features. Thus, resistance to phages in the wild follows evolutionary trajectories alternative to those predicted from laboratory-based evolutionary experiments.

Recombinant SARS-CoV-2 genomes are currently circulating at low levels.

Viral recombination can generate novel genotypes with unique phenotypic characteristics, including transmissibility and virulence. Although the capacity for recombination among betacoronaviruses is well documented, there is limited evidence of recombination between SARS-CoV-2 strains. By identifying the mutations that primarily determine SARS-CoV-2 clade structure, we developed a lightweight approach for detecting recombinant genomes. Among the over 537,000 genomes queried, we detect 1175 putative recombinants that contain multiple mutational markers from distinct clades. Additional phylogenetic analysis and the observed co-circulation of predicted parent clades in the geographic regions of exposure further support the feasibility of recombination in these detected cases. An analysis of these detected cases did not reveal any evidence for recombination hotspots in the SARS-CoV-2 genome. Although most recombinant genotypes were detected a limited number of times, at least two recombinants are now widely transmitted. Recombinant genomes were also found to contain substitutions of concern for elevated transmissibility and lower vaccine efficacy, including D614G, N501Y, E484K, and L452R. Adjusting for an unequal probability of detecting recombinants derived from different parent clades, and for geographic variation in clade abundance, we estimate that at most 5% of circulating viruses in the USA and UK are recombinant. While the phenotypic characterization of detected recombinants was beyond the scope of our analysis, the identification of transmitted recombinants involving substitutions of concern underscores the need to sustain efforts to monitor the emergence of new genotypes generated through recombination.

Recombinant SARS-CoV-2 genomes circulated at low levels over the first year of the pandemic.

Viral recombination can generate novel genotypes with unique phenotypic characteristics, including transmissibility and virulence. Although the capacity for recombination among betacoronaviruses is well documented, recombination between strains of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has not been characterized in detail. Here, we present a lightweight approach for detecting genomes that are potentially recombinant. This approach relies on identifying the mutations that primarily determine SARS-CoV-2 clade structure and then screening genomes for ones that contain multiple mutational markers from distinct clades. Among the over 537,000 genomes queried that were deposited on GISAID.org prior to 16 February 2021, we detected 1,175 potential recombinant sequences. Using a highly conservative criteria to exclude sequences that may have originated through de novo mutation, we find that at least 30 per cent (n = 358) are likely of recombinant origin. An analysis of deep-sequencing data for these putative recombinants, where available, indicated that the majority are high quality. Additional phylogenetic analysis and the observed co-circulation of predicted parent clades in the geographic regions of exposure further support the feasibility of recombination in this subset of potential recombinants. An analysis of these genomes did not reveal evidence for recombination hotspots in the SARS-CoV-2 genome. While most of the putative recombinant sequences we detected were genetic singletons, a small number of genetically identical or highly similar recombinant sequences were identified in the same geographic region, indicative of locally circulating lineages. Recombinant genomes were also found to have originated from parental lineages with substitutions of concern, including D614G, N501Y, E484K, and L452R. Adjusting for an unequal probability of detecting recombinants derived from different parent clades and for geographic variation in clade abundance, we estimate that at most 0.2-2.5 per cent of circulating viruses in the USA and UK are recombinant. Our identification of a small number of putative recombinants within the first year of SARS-CoV-2 circulation underscores the need to sustain efforts to monitor the emergence of new genotypes generated through recombination.

Topological analysis reveals state transitions in human gut and marine bacterial communities.

Microbiome dynamics influence the health and functioning of human physiology and the environment and are driven in part by interactions between large numbers of microbial taxa, making large-scale prediction and modeling a challenge. Here, using topological data analysis, we identify states and dynamical features relevant to macroscopic processes. We show that gut disease processes and marine geochemical events are associated with transitions between community states, defined as topological features of the data density. We find a reproducible two-state succession during recovery from cholera in the gut microbiomes of multiple patients, evidence of dynamic stability in the gut microbiome of a healthy human after experiencing diarrhea during travel, and periodic state transitions in a marine Prochlorococcus community driven by water column cycling. Our approach bridges small-scale fluctuations in microbiome composition and large-scale changes in phenotype without details of underlying mechanisms, and provides an assessment of microbiome stability and its relation to human and environmental health.

A Human Microbiota-Associated Murine Model for Assessing the Impact of the Vaginal Microbiota on Pregnancy Outcomes.

Disease states are often linked to large scale changes in microbial community structure that obscure the contributions of individual microbes to disease. Establishing a mechanistic understanding of how microbial community structure contribute to certain diseases, however, remains elusive thereby limiting our ability to develop successful microbiome-based therapeutics. Human microbiota-associated (HMA) mice have emerged as a powerful approach for directly testing the influence of microbial communities on host health and disease, with the transfer of disease phenotypes from humans to germ-free recipient mice widely reported. We developed a HMA mouse model of the human vaginal microbiota to interrogate the effects of Bacterial Vaginosis (BV) on pregnancy outcomes. We collected vaginal swabs from 19 pregnant African American women with and without BV (diagnosed per Nugent score) to colonize female germ-free mice and measure its impact on birth outcomes. There was considerable variability in the microbes that colonized each mouse, with no association to the BV status of the microbiota donor. Although some of the women in the study had adverse birth outcomes, the vaginal microbiota was not predictive of adverse birth outcomes in mice. However, elevated levels of pro-inflammatory cytokines in the uterus of HMA mice were detected during pregnancy. Together, these data outline the potential uses and limitations of HMA mice to elucidate the influence of the vaginal microbiota on health and disease.

How can microbial population genomics inform community ecology?

Populations are fundamental units of ecology and evolution, but can we define them for bacteria and archaea in a biologically meaningful way? Here, we review why population structure is difficult to recognize in microbes and how recent advances in measuring contemporary gene flow allow us to identify clearly delineated populations among collections of closely related genomes. Such structure can arise from preferential gene flow caused by coexistence and genetic similarity, defining populations based on biological mechanisms. We show that such gene flow units are sufficiently genetically isolated for specific adaptations to spread, making them also ecological units that are differentially adapted compared to their closest relatives. We discuss the implications of these observations for measuring bacterial and archaeal diversity in the environment. We show that operational taxonomic units defined by 16S rRNA gene sequencing have woefully poor resolution for ecologically defined populations and propose monophyletic clusters of nearly identical ribosomal protein genes as an alternative measure for population mapping in community ecological studies employing metagenomics. These population-based approaches have the potential to provide much-needed clarity in interpreting the vast microbial diversity in human and environmental microbiomes. This article is part of the theme issue 'Conceptual challenges in microbial community ecology'.

Diarrhoeal events can trigger long-term Clostridium difficile colonization with recurrent blooms.

Although Clostridium difficile is widely considered an antibiotic- and hospital-associated pathogen, recent evidence indicates that this is an insufficient depiction of the risks and reservoirs. A common thread that links all major risk factors of infection is their association with gastrointestinal disturbances, but this relationship to C. difficile colonization has never been tested directly. Here, we show that disturbances caused by diarrhoeal events trigger susceptibility to C. difficile colonization. Using survey data of the human gut microbiome, we detected C. difficile colonization and blooms in people recovering from food poisoning and Vibrio cholerae infections. Carriers remained colonized for year-long time scales and experienced highly variable patterns of C. difficile abundance, where increased shedding over short periods of 1-2 d interrupted week-long periods in which C. difficile was undetectable. Given that short shedding events were often linked to gastrointestinal disturbances, our results help explain why C. difficile is frequently detected as a co-infecting pathogen in patients with diarrhoea. To directly test the impact of diarrhoea on susceptibility to colonization, we developed a mouse model of variable disturbance intensity, which allowed us to monitor colonization in the absence of disease. As mice exposed to avirulent C. difficile spores ingested increasing quantities of laxatives, more individuals experienced C. difficile blooms. Our results indicate that the likelihood of colonization is highest in the days immediately following acute disturbances, suggesting that this could be an important window during which transmission could be interrupted and the incidence of infection lowered.

A Reverse Ecology Approach Based on a Biological Definition of Microbial Populations.

Delineating ecologically meaningful populations among microbes is important for identifying their roles in environmental and host-associated microbiomes. Here, we introduce a metric of recent gene flow, which when applied to co-existing microbes, identifies congruent genetic and ecological units separated by strong gene flow discontinuities from their next of kin. We then develop a pipeline to identify genome regions within these units that show differential adaptation and allow mapping of populations onto environmental variables or host associations. Using this reverse ecology approach, we show that the human commensal bacterium Ruminococcus gnavus breaks up into sharply delineated populations that show different associations with health and disease. Defining populations by recent gene flow in this way will facilitate the analysis of bacterial and archaeal genomes using ecological and evolutionary theory developed for plants and animals, thus allowing for testing unifying principles across all biology.

Viruses of the Nahant Collection, characterization of 251 marine Vibrionaceae viruses.

Viruses are highly discriminating in their interactions with host cells and are thought to play a major role in maintaining diversity of environmental microbes. However, large-scale ecological and genomic studies of co-occurring virus-host pairs, required to characterize the mechanistic and genomic foundations of virus-host interactions, are lacking. Here, we present the largest dataset of cultivated and sequenced co-occurring virus-host pairs that captures ecologically representative fine-scale diversity. Using the ubiquitous and ecologically diverse marine Vibrionaceae as a host platform, we isolate and sequence 251 dsDNA viruses and their hosts from three time points within a 93-day time-series study. The virus collection includes representatives of the three Caudovirales tailed virus morphotypes, a novel family of nontailed viruses, and the smallest (10,046 bp) and largest (348,911 bp) Vibrio virus genomes described. We provide general characterization and annotation of the viruses and describe read-mapping protocols to standardize genome presentation. The rich ecological and genomic contextualization of hosts and viruses make the Nahant Collection a unique platform for high-resolution studies of environmental virus-host infection networks.

A major lineage of non-tailed dsDNA viruses as unrecognized killers of marine bacteria.

The most abundant viruses on Earth are thought to be double-stranded DNA (dsDNA) viruses that infect bacteria. However, tailed bacterial dsDNA viruses (Caudovirales), which dominate sequence and culture collections, are not representative of the environmental diversity of viruses. In fact, non-tailed viruses often dominate ocean samples numerically, raising the fundamental question of the nature of these viruses. Here we characterize a group of marine dsDNA non-tailed viruses with short 10-kb genomes isolated during a study that quantified the diversity of viruses infecting Vibrionaceae bacteria. These viruses, which we propose to name the Autolykiviridae, represent a novel family within the ancient lineage of double jelly roll (DJR) capsid viruses. Ecologically, members of the Autolykiviridae have a broad host range, killing on average 34 hosts in four Vibrio species, in contrast to tailed viruses which kill on average only two hosts in one species. Biochemical and physical characterization of autolykiviruses reveals multiple virion features that cause systematic loss of DJR viruses in sequencing and culture-based studies, and we describe simple procedural adjustments to recover them. We identify DJR viruses in the genomes of diverse major bacterial and archaeal phyla, and in marine water column and sediment metagenomes, and find that their diversity greatly exceeds the diversity that is currently captured by the three recognized families of such viruses. Overall, these data suggest that viruses of the non-tailed dsDNA DJR lineage are important but often overlooked predators of bacteria and archaea that impose fundamentally different predation and gene transfer regimes on microbial systems than on tailed viruses, which form the basis of all environmental models of bacteria-virus interactions.

Population genomics of the symbiotic plasmids of sympatric nitrogen-fixing Rhizobium species associated with Phaseolus vulgaris.

Cultivated common beans are the primary protein source for millions of people around the world who subsist on low-input agriculture, enabled by the symbiotic N2 -fixation these legumes perform in association with rhizobia. Within a single agricultural plot, multiple Rhizobium species can nodulate bean roots, but it is unclear how genetically isolated these species remain in sympatry. To better understand this issue, we sequenced and compared the genomes of 33 strains isolated from the rhizosphere and root nodules of a particular bean variety grown in the same agricultural plot. We found that the Rhizobium species we observed coexist with low genetic recombination across their core genomes. Accessory plasmids thought to be necessary for the saprophytic lifestyle in soil show similar levels of genetic isolation, but with higher rates of recombination than the chromosomes. However, the symbiotic plasmids are extremely similar, with high rates of recombination and do not appear to have co-evolved with the chromosome or accessory plasmids. Therefore, while Rhizobium species are genetically isolated units within the microbial community, a common symbiotic plasmid allows all Rhizobium species to engage in symbiosis with the same host in a single agricultural plot.