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.

Abstract IA20: 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.

Citation Format: Robin Reschke, Jason W Shapiro, Jovian Yu, Sherin J Rouhani, Daniel J Olson, Yuanyuan Zha, Thomas F Gajewski. Checkpoint blockade–induced dermatitis and colitis are dominated by tissue-resident memory T cells and Th1/Tc1 cytokines [abstract]. In: Proceedings of the AACR Special Conference: Tumor Immunology and Immunotherapy; 2022 Oct 21-24; Boston, MA. Philadelphia (PA): AACR; Cancer Immunol Res 2022;10(12 Suppl):Abstract nr IA20.

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.

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.

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.

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.

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.

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.

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.

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.

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.