Batf3+ DCs and the 4-1BB/4-1BBL axis are required at the effector phase in the tumor microenvironment for PD-1/PD-L1 blockade efficacy.

The cellular source of positive signals that reinvigorate T cells within the tumor microenvironment (TME) for the therapeutic efficacy of programmed death-1 (PD-1)/programmed death-ligand 1 (PD-L1) blockade has not been clearly defined. We now show that Batf3-lineage dendritic cells (DCs) are essential in this process. Flow cytometric analysis, gene-targeted mice, and blocking antibody studies revealed that 4-1BBL is a major positive co-stimulatory signal provided by these DCs within the TME that translates to CD8+ T cell functional reinvigoration and tumor regression. Immunofluorescence and spatial transcriptomics on human tumor samples revealed clustering of Batf3+ DCs and CD8+ T cells, which correlates with anti-PD-1 efficacy. In addition, proximity to Batf3+ DCs within the TME is associated with CD8+ T cell transcriptional states linked to anti-PD-1 response. Our results demonstrate that Batf3+ DCs within the TME are critical for PD-1/PD-L1 blockade efficacy and indicate a major role for the 4-1BB/4-1BB ligand (4-1BBL) axis during this process.

Checkpoint Blockade-Induced Dermatitis and Colitis Are Dominated by Tissue-Resident Memory T Cells and Th1/Tc1 Cytokines.

Immune checkpoint blockade is therapeutically successful for many patients across multiple cancer types. However, immune-related adverse events (irAE) frequently occur and can sometimes be life threatening. It is critical to understand the immunologic mechanisms of irAEs with the goal of finding novel treatment targets. Herein, we report our analysis of tissues from patients with irAE dermatitis using multiparameter immunofluorescence (IF), spatial transcriptomics, and RNA in situ hybridization (RISH). Skin psoriasis cases were studied as a comparison, as a known Th17-driven disease, and colitis was investigated as a comparison. IF analysis revealed that CD4+ and CD8+ tissue-resident memory T (TRM) cells were preferentially expanded in the inflamed portion of skin in cutaneous irAEs compared with healthy skin controls. Spatial transcriptomics allowed us to focus on areas containing TRM cells to discern functional phenotype and revealed expression of Th1-associated genes in irAEs, compared with Th17-asociated genes in psoriasis. Expression of PD-1, CTLA-4, LAG-3, and other inhibitory receptors was observed in irAE cases. RISH technology combined with IF confirmed expression of IFNγ, CXCL9, CXCL10, and TNFα in irAE dermatitis, as well as IFNγ within TRM cells specifically. The Th1-skewed phenotype was confirmed in irAE colitis cases compared with healthy colon.

Rephine.r: a pipeline for correcting gene calls and clusters to improve phage pangenomes and phylogenies.

BACKGROUND: A pangenome is the collection of all genes found in a set of related genomes. For microbes, these genomes are often different strains of the same species, and the pangenome offers a means to compare gene content variation with differences in phenotypes, ecology, and phylogenetic relatedness. Though most frequently applied to bacteria, there is growing interest in adapting pangenome analysis to bacteriophages. However, working with phage genomes presents new challenges. First, most phage families are under-sampled, and homologous genes in related viruses can be difficult to identify. Second, homing endonucleases and intron-like sequences may be present, resulting in fragmented gene calls. Each of these issues can reduce the accuracy of standard pangenome analysis tools. METHODS: We developed an R pipeline called Rephine.r that takes as input the gene clusters produced by an initial pangenomics workflow. Rephine.r then proceeds in two primary steps. First, it identifies three common causes of fragmented gene calls: (1) indels creating early stop codons and new start codons; (2) interruption by a selfish genetic element; and (3) splitting at the ends of the reported genome. Fragmented genes are then fused to create new sequence alignments. In tandem, Rephine.r searches for distant homologs separated into different gene families using Hidden Markov Models. Significant hits are used to merge families into larger clusters. A final round of fragment identification is then run, and results may be used to infer single-copy core genomes and phylogenetic trees. RESULTS: We applied Rephine.r to three well-studied phage groups: the Tevenvirinae (e.g., T4), the Studiervirinae (e.g., T7), and the Pbunaviruses (e.g., PB1). In each case, Rephine.r recovered additional members of the single-copy core genome and increased the overall bootstrap support of the phylogeny. The Rephine.r pipeline is provided through GitHub (https://www.github.com/coevoeco/Rephine.r) as a single script for automated analysis and with utility functions to assist in building single-copy core genomes and predicting the sources of fragmented genes.

UPΦ phages, a new group of filamentous phages found in several members of Enterobacteriales.

Filamentous phages establish chronic infections in their bacterial hosts, and new phages are secreted by infected bacteria for multiple generations, typically without causing host death. Often, these viruses integrate in their host's genome by co-opting the host's XerCD recombinase system. In several cases, these viruses also encode genes that increase bacterial virulence in plants and animals. Here, we describe a new filamentous phage, UPϕ901, which we originally found integrated in a clinical isolate of Escherichia coli from urine. UPϕ901 and closely related phages can be found in published genomes of over 200 other bacteria, including strains of Citrobacter koseri, Salmonella enterica, Yersinia enterocolitica, and Klebsiella pneumoniae. Its closest relatives are consistently found in urine or in the blood and feces of patients with urinary tract infections. More distant relatives can be found in isolates from other environments, including sewage, water, soil, and contaminated food. Each of these phages, which we collectively call 'UPϕ viruses', also harbors two or more novel genes of unknown function.

Comparative Genomic Study of Lactobacillus jensenii and the Newly Defined Lactobacillus mulieris Species Identifies Species-Specific Functionality.

Lactobacilli are dominant members of the "healthy" female urogenital microbiota. One of these species, Lactobacillus jensenii, is routinely identified in the urinary microbiota of women both with and without urinary tract symptoms. In March 2020, the new bacterial species Lactobacillus mulieris was introduced, and phylogenetic and average nucleotide identity analysis identified eight L. jensenii strains that should be classified as members of the L. mulieris species. This prompted our phylogenomic study of all publicly available L. jensenii and L. mulieris genome sequences. While there is little variation in the 16S rRNA gene sequences, the core genome shows a clear distinction between genomes of the two species. We find eight additional strains of the species L. mulieris among these genomes. Furthermore, one strain, currently classified as L. mulieris UMB7784, is distinct from both L. jensenii and L. mulieris strains. As part of our comparative genomic study, we also investigated the genetic content that distinguishes these two species. Unique to the L. jensenii genomes are several genes related to catabolism of disaccharides. In contrast, L. mulieris genomes encode several cell surface and secreted proteins that are not found within the L. jensenii genomes. These L. jensenii-specific and L. mulieris-specific loci provide insight into phenotypic differences of these two species.IMPORTANCELactobacillus species play a key role in the health of the urinary tract. For instance, Lactobacillus crispatus and L. jensenii have been found to inhibit uropathogenic Escherichia coli growth. While L. crispatus is typically found only within the microbiota of women without lower urinary tract symptoms (LUTS), L. jensenii has been found in the microbiota of women both with and without LUTS. With the recent introduction of the new species Lactobacillus mulieris, several strains of L. jensenii were reclassified as L. mulieris based upon gene marker and average nucleotide identity. We took a phylogenomic and comparative genomic approach to ascertain the genetic determinants of these two species. Looking at a larger data set, we identified additional L. mulieris strains, including one distinct from other members of the species-L. mulieris UMB7784. Furthermore, we identified unique loci in each species that may have clinical implications.

Introducing Lu-1, a Novel Lactobacillus jensenii Phage Abundant in the Urogenital Tract.

Bacteriophages (phages) play a key role in shaping microbial communities, including those of the human body. Phages are abundant members of the urogenital tract, most often persisting through the lysogenic life cycle as prophages integrated within the genomes of their bacterial hosts. While numerous studies of the urogenital microbiota have focused on the most abundant bacterial member of this niche-Lactobacillus species-very little is known about Lactobacillus phages. Focusing on Lactobacillus jensenii strains from the urinary tract, we identified numerous prophages related to the previously characterized Lv-1 phage from a vaginal L. jensenii strain. Furthermore, we identified a new L. jensenii phage, Lu-1. Evidence suggests that both phages are abundant within the urogenital tract. CRISPR spacer sequences matching to Lv-1 and Lu-1 prophages were identified. While first detected in urinary isolates, the Lu-1 phage was also discovered in L. jensenii isolates from vaginal and perineal swabs, and both phages were found in metagenomic data sets. The prevalence of these phages in the isolates suggests that both phages are active members of the urogenital microbiota.

Genomic Survey of E. coli From the Bladders of Women With and Without Lower Urinary Tract Symptoms.

Urinary tract infections (UTIs) are one of the most common human bacterial infections. While UTIs are commonly associated with colonization by Escherichia coli, members of this species also have been found within the bladder of individuals with no lower urinary tract symptoms (no LUTS), also known as asymptomatic bacteriuria. Prior studies have found that both uropathogenic E. coli (UPEC) strains and E. coli isolates that are not associated with UTIs encode for virulence factors. Thus, the reason(s) why E. coli sometimes causes UTI-like symptoms remain(s) elusive. In this study, the genomes of 66 E. coli isolates from adult female bladders were sequenced. These isolates were collected from four cohorts, including women: (1) without lower urinary tract symptoms, (2) overactive bladder symptoms, (3) urgency urinary incontinence, and (4) a clinical diagnosis of UTI. Comparative genomic analyses were conducted, including core and accessory genome analyses, virulence and motility gene analyses, and antibiotic resistance prediction and testing. We found that the genomic content of these 66 E. coli isolates does not correspond with the participant's symptom status. We thus looked beyond the E. coli genomes to the composition of the entire urobiome and found that the presence of E. coli alone was not sufficient to distinguish between the urobiomes of individuals with UTI and those with no LUTS. Because E. coli presence, abundance, and genomic content appear to be weak predictors of UTI status, we hypothesize that UTI symptoms associated with detection of E. coli are more likely the result of urobiome composition.

Gene Co-occurrence Networks Reflect Bacteriophage Ecology and Evolution.

Bacteriophages are the most abundant and diverse biological entities on the planet, and new phage genomes are being discovered at a rapid pace. As more phage genomes are published, new methods are needed for placing these genomes in an ecological and evolutionary context. Phages are difficult to study by phylogenetic methods, because they exchange genes regularly, and no single gene is conserved across all phages. Here, we demonstrate how gene-level networks can provide a high-resolution view of phage genetic diversity and offer a novel perspective on virus ecology. We focus our analyses on virus host range and show how network topology corresponds to host relatedness, how to find groups of genes with the strongest host-specific signatures, and how this perspective can complement phage host prediction tools. We discuss extensions of gene network analysis to predicting the emergence of phages on new hosts, as well as applications to features of phage biology beyond host range.IMPORTANCE Bacteriophages (phages) are viruses that infect bacteria, and they are critical drivers of bacterial evolution and community structure. It is generally difficult to study phages by using tree-based methods, because gene exchange is common, and no single gene is shared among all phages. Instead, networks offer a means to compare phages while placing them in a broader ecological and evolutionary context. In this work, we build a network that summarizes gene sharing across phages and test how a key constraint on phage ecology, host range, corresponds to the structure of the network. We find that the network reflects the relatedness among phage hosts, and phages with genes that are closer in the network are likelier to infect similar hosts. This approach can also be used to identify genes that affect host range, and we discuss possible extensions to analyze other aspects of viral ecology.

Evolution of mutualism from parasitism in experimental virus populations.

While theory suggests conditions under which mutualism may evolve from parasitism, few studies have observed this transition empirically. Previously, we evolved Escherichia coli and the filamentous bacteriophage M13 in 96-well microplates, an environment in which the ancestral phage increased the growth rate and yield of the ancestral bacteria. In the majority of populations, mutualism was maintained or even enhanced between phages and coevolving bacteria; however, these same phages evolved traits that harmed the ancestral E. coli genotype. Here, we set out to determine if mutualism could evolve from this new parasitic interaction. To do so, we chose six evolved phage populations from the original experiment and used them to establish new infections of the ancestral bacteria. After 20 passages, mutualism evolved in almost all replicates, with the remainder growing commensally. Many phage populations also evolved to benefit both their local, evolving bacteria and the ancestral bacteria, though these phages were less beneficial to their co-occurring hosts than phages that harm the ancestral bacteria. These results demonstrate the rapid recovery of mutualism from parasitism, and we discuss how our findings relate to the evolution of phages that enhance the virulence of bacterial pathogens.

Adding biotic complexity alters the metabolic benefits of mutualism.

Mutualism is ubiquitous in nature and plays an integral role in most communities. To predict the eco-evolutionary dynamics of mutualism it is critical to extend classic pair-wise analysis to include additional species. We investigated the effect of adding a third species to a pair-wise mutualism in a spatially structured environment. We tested the hypotheses that selection for costly excretions in a focal population (i) decreases when an exploiter is added (ii) increases when a third mutualist is added relative to the pair-wise scenario. We assayed the selection acting on Salmonella enterica when it exchanges methionine for carbon in an obligate mutualism with an auxotrophic Escherichia coli. A third bacterium, Methylobacterium extorquens, was then added and acted either as an exploiter of the carbon or third obligate mutualist depending on the nitrogen source. In the tripartite mutualism M. extorquens provided nitrogen to the other species. Contrary to our expectations, adding an exploiter increased selection for methionine excretion in S. enterica. Conversely, selection for cooperation was lower in the tripartite mutualism relative to the pair-wise system. Genome-scale metabolic models helped identify the mechanisms underlying these changes in selection. Our results highlight the utility of connecting metabolic mechanisms and eco-evolutionary dynamics.

Evolution of parasitism and mutualism between filamentous phage M13 and Escherichia coli.

Background. How host-symbiont interactions coevolve between mutualism and parasitism depends on the ecology of the system and on the genetic and physiological constraints of the organisms involved. Theory often predicts that greater reliance on horizontal transmission favors increased costs of infection and may result in more virulent parasites or less beneficial mutualists. We set out to understand transitions between parasitism and mutualism by evolving the filamentous bacteriophage M13 and its host Escherichia coli. Results. The effect of phage M13 on bacterial fitness depends on the growth environment, and initial assays revealed that infected bacteria reproduce faster and to higher density than uninfected bacteria in 96-well microplates. These data suggested that M13 is, in fact, a facultative mutualist of E. coli. We then allowed E. coli and M13 to evolve in replicated environments, which varied in the relative opportunity for horizontal and vertical transmission of phage in order to assess the evolutionary stability of this mutualism. After 20 experimental passages, infected bacteria from treatments with both vertical and horizontal transmission of phage had evolved the fastest growth rates. At the same time, phage from these treatments no longer benefited the ancestral bacteria. Conclusions. These data suggest a positive correlation between the positive effects of M13 on E. coli hosts from the same culture and the negative effects of the same phage toward the ancestral bacterial genotype. The results also expose flaws in applying concepts from the virulence-transmission tradeoff hypothesis to mutualism evolution. We discuss the data in the context of more recent theory on how horizontal transmission affects mutualisms and explore how these effects influence phages encoding virulence factors in pathogenic bacteria.

Freshwater Metaviromics and Bacteriophages: A Current Assessment of the State of the Art in Relation to Bioinformatic Challenges.

Advances in bioinformatics and sequencing technologies have allowed for the analysis of complex microbial communities at an unprecedented rate. While much focus is often placed on the cellular members of these communities, viruses play a pivotal role, particularly bacteria-infecting viruses (bacteriophages); phages mediate global biogeochemical processes and drive microbial evolution through bacterial grazing and horizontal gene transfer. Despite their importance and ubiquity in nature, very little is known about the diversity and structure of viral communities. Though the need for culture-based methods for viral identification has been somewhat circumvented through metagenomic techniques, the analysis of metaviromic data is marred with many unique issues. In this review, we examine the current bioinformatic approaches for metavirome analyses and the inherent challenges facing the field as illustrated by the ongoing efforts in the exploration of freshwater phage populations.

Genomes of Gardnerella Strains Reveal an Abundance of Prophages within the Bladder Microbiome.

Bacterial surveys of the vaginal and bladder human microbiota have revealed an abundance of many similar bacterial taxa. As the bladder was once thought to be sterile, the complex interactions between microbes within the bladder have yet to be characterized. To initiate this process, we have begun sequencing isolates, including the clinically relevant genus Gardnerella. Herein, we present the genomic sequences of four Gardnerella strains isolated from the bladders of women with symptoms of urgency urinary incontinence; these are the first Gardnerella genomes produced from this niche. Congruent to genomic characterization of Gardnerella isolates from the reproductive tract, isolates from the bladder reveal a large pangenome, as well as evidence of high frequency horizontal gene transfer. Prophage gene sequences were found to be abundant amongst the strains isolated from the bladder, as well as amongst publicly available Gardnerella genomes from the vagina and endometrium, motivating an in depth examination of these sequences. Amongst the 39 Gardnerella strains examined here, there were more than 400 annotated prophage gene sequences that we could cluster into 95 homologous groups; 49 of these groups were unique to a single strain. While many of these prophages exhibited no sequence similarity to any lytic phage genome, estimation of the rate of phage acquisition suggests both vertical and horizontal acquisition. Furthermore, bioinformatic evidence indicates that prophage acquisition is ongoing within both vaginal and bladder Gardnerella populations. The abundance of prophage sequences within the strains examined here suggests that phages could play an important role in the species' evolutionary history and in its interactions within the complex communities found in the female urinary and reproductive tracts.

The impact of transmission mode on the evolution of benefits provided by microbial symbionts.

While past work has often examined the effects of transmission mode on virulence evolution in parasites, few studies have explored the impact of horizontal transmission on the evolution of benefits conferred by a symbiont to its host. Here, we identify three mechanisms that create a positive covariance between horizontal transmission and symbiont-provided benefits: pleiotropy within the symbiont genome, partner choice by the host, and consumption of host waste by-products by symbionts. We modify a susceptible-infected model to incorporate the details of each mechanism and examine the evolution of symbiont benefits given variation in either the immigration rate of susceptible hosts or the rate of successful vertical transmission. We find conditions for each case under which greater opportunity for horizontal transmission (higher migration rate) favors the evolution of mutualism. Further, we find the surprising result that vertical transmission can inhibit the evolution of benefits provided by symbionts to hosts when horizontal transmission and symbiont-provided benefits are positively correlated. These predictions may apply to a number of natural systems, and the results may explain why many mutualisms that rely on partner choice often lack a mechanism for vertical transmission.