Search

Supplementary table from Accounting for B-cell Behavior and Sampling Bias Predicts Anti–PD-L1 Response in Bladder Cancer

Cancer immunotherapy is predominantly based on T cell–centric approaches. At the same time, the adaptive immune response in the tumor environment also includes clonally produced immunoglobulins and clonal effector/memory B cells that participate in antigen-specific decisions through their interactions with T cells. Here, we investigated the role of infiltrating B cells in bladder cancer via patient dataset analysis of intratumoral immunoglobulin repertoires. We showed that the IgG1/IgA ratio is a prognostic indicator for several subtypes of bladder cancer and for the whole IMVigor210 anti–PD-L1 immunotherapy study cohort. A high IgG1/IgA ratio associated with the prominence of a cytotoxic gene signature, T-cell receptor signaling, and IL21-mediated signaling. Immunoglobulin repertoire analysis indicated that effector B-cell function, rather than clonally produced antibodies, was involved in antitumor responses. From the T-cell side, we normalized a cytotoxic signature against the extent of immune cell infiltration to neutralize the artificial sampling-based variability in immune gene expression. Resulting metrics reflected proportion of cytotoxic cells among tumor-infiltrating immune cells and improved prediction of anti–PD-L1 responses. At the same time, the IgG1/IgA ratio remained an independent prognostic factor. Integration of the B-cell, natural killer cell, and T-cell signatures allowed for the most accurate prediction of anti–PD-L1 therapy responses. On the basis of these findings, we developed a predictor called PRedIctive MolecUlar Signature (PRIMUS), which outperformed PD-L1 expression scores and known gene signatures. Overall, PRIMUS allows for reliable identification of responders among patients with muscle-invasive urothelial carcinoma, including the subcohort with the low-infiltrated “desert” tumor phenotype.

Supplementary Data from Accounting for B-cell Behavior and Sampling Bias Predicts Anti–PD-L1 Response in Bladder Cancer

T cell maturation and activation depend upon T cell receptor (TCR) interactions with a wide variety of antigenic peptides displayed in a given major histocompatibility complex (MHC) context. Complementarity-determining region 3 (CDR3) is the most variable part of the TCRα and -β chains, which govern interactions with peptide–MHC complexes. However, it remains unclear how the CDR3 landscape is shaped by individual MHC context during thymic selection of naïve T cells. We established two mouse strains carrying distinct allelic variants of H2-A and analyzed thymic and peripheral production and TCR repertoires of naïve conventional CD4 + T (T conv ) and naïve regulatory CD4 + T (T reg ) cells. Compared with tuberculosis-resistant C57BL/6 (H2-A b ) mice, the tuberculosis-susceptible H2-A j mice had fewer CD4 + T cells of both subsets in the thymus. In the periphery, this deficiency was only apparent for T conv and was compensated for by peripheral reconstitution for T reg . We show that H2-A j favors selection of a narrower and more convergent repertoire with more hydrophobic and strongly interacting amino acid residues in the middle of CDR3α and CDR3β, suggesting more stringent selection against a narrower peptide–MHC-II context. H2-A j and H2-A b mice have prominent reciprocal differences in CDR3α and CDR3β features, probably reflecting distinct modes of TCR fitting to MHC-II variants. These data reveal the mechanics and extent of how MHC-II shapes the naïve CD4 + T cell CDR3 landscape, which essentially defines adaptive response to infections and self-antigens.

The purpose of this work is to adapt and apply Next Generation Sequencing methods in oil and gas well field studies. Relatively recent NGS methods provide a description of a geological formation by analyzing millions of DNA sequences and represent an entirely new way to obtain information about oil and gas reservoirs and the composition of their fluids, which could significantly change the approach to exploration and field development. We present the results of pilot work to determine the inflow profile in a horizontal well based on DNA markers. The technology is based on the comparison of bacterial DNA from drill cuttings obtained while drilling with DNA from microorganisms of fluids obtained during production at the wellhead. Because of their high selectivity, individual microbes live only under certain conditions (salinity, oil saturation, temperature) and can be used as unique natural biomarkers. The comparison of DNA samples of drilling cutting and produced fluid allows for identification of the segment of the horizontal well from which the main flow comes, as well as identifying the type of incoming fluid (water, oil, gas) without stopping the operation process and without conducting expensive downhole operations. As a result of these studies, the microbial communities of the oil-bearing sands and formation fluids of the Cretaceous deposits (group BS) in Western Siberia were identified, and the relative numerical ratio of microorganisms in the formations was determined. It was shown that the microbiome diversity changes with depth, and depends on the lithological composition, and sequencing data obtained from cuttings samples correlate with data from wellhead samples of produced fluid. Thus, the practical applicability of DNA sequencing for solving field problems in oil and gas field development, in particular for determining the inflow profile in horizontal wells, was confirmed.

Cancer immunotherapy is predominantly based on T cell–centric approaches. At the same time, the adaptive immune response in the tumor environment also includes clonally produced immunoglobulins and clonal effector/memory B cells that participate in antigen-specific decisions through their interactions with T cells. Here, we investigated the role of infiltrating B cells in bladder cancer via patient dataset analysis of intratumoral immunoglobulin repertoires. We showed that the IgG1/IgA ratio is a prognostic indicator for several subtypes of bladder cancer and for the whole IMVigor210 anti–PD-L1 immunotherapy study cohort. A high IgG1/IgA ratio associated with the prominence of a cytotoxic gene signature, T-cell receptor signaling, and IL21-mediated signaling. Immunoglobulin repertoire analysis indicated that effector B-cell function, rather than clonally produced antibodies, was involved in antitumor responses. From the T-cell side, we normalized a cytotoxic signature against the extent of immune cell infiltration to neutralize the artificial sampling-based variability in immune gene expression. Resulting metrics reflected proportion of cytotoxic cells among tumor-infiltrating immune cells and improved prediction of anti–PD-L1 responses. At the same time, the IgG1/IgA ratio remained an independent prognostic factor. Integration of the B-cell, natural killer cell, and T-cell signatures allowed for the most accurate prediction of anti–PD-L1 therapy responses. On the basis of these findings, we developed a predictor called PRedIctive MolecUlar Signature (PRIMUS), which outperformed PD-L1 expression scores and known gene signatures. Overall, PRIMUS allows for reliable identification of responders among patients with muscle-invasive urothelial carcinoma, including the subcohort with the low-infiltrated “desert” tumor phenotype.

Background Chlamydia are ancient intracellular pathogens with reduced, though strikingly conserved genome. Despite their parasitic lifestyle and isolated intracellular environment, these bacteria managed to avoid accumulation of deleterious mutations leading to subsequent genome degradation characteristic for many parasitic bacteria. Results We report pan-genomic analysis of sixteen species from genus Chlamydia including identification and functional annotation of orthologous genes, and characterization of gene gains, losses, and rearrangements. We demonstrate the overall genome stability of these bacteria as indicated by a large fraction of common genes with conserved genomic locations. On the other hand, extreme evolvability is confined to several paralogous gene families such as polymorphic membrane proteins and phospholipase D, and likely is caused by the pressure from the host immune system. Conclusions This combination of a large, conserved core genome and a small, evolvable periphery likely reflect the balance between the selective pressure towards genome reduction and the need to adapt to escape from the host immunity.

Tumor-infiltrating B cells and intratumorally-produced immunoglobulins (IG) play important roles in the tumor microenvironment and response to immunotherapy1–5. IgG antibodies produced by intratumoral B cells may drive antibody-dependent cellular cytotoxicity (ADCC) and enhance antigen presentation by dendritic cells6–8. Furthermore, B cells are efficient antigen-specific antigen presenters that can essentially modulate the behaviour of helper T cells9–11.Here we investigated the role of intratumoral IG isotype and clonality in bladder cancer. Our results show that the IgG1/IgA ratio offers a strong and independent prognostic indicator for the Basal squamous molecular subtype and for the whole ImVigor210 cohort in anti-PD-L1 immunotherapy. Our findings also indicate that effector B cell functions, rather than clonally-produced antibodies, are involved in the antitumor response. High IgG1/IgA ratio was associated with relative abundance of cytotoxic genes and prominence of the IL-21/IL-21R axis suggesting importance of T cell/B cell interaction.We integrated the B, NK, and T cell components, employing immFocus-like normalization to account for the stochastic nature of tumor tissue sampling. Using a random forest model with nested cross-validation, we developed a tumor RNA-Seq-based predictor of anti-PD-L1 therapy response in muscle-invasive urothelial carcinoma. The resulting PRIMUS (PRedIctive MolecUlar Signature) predictor achieves superior sensitivity compared to PD-L1 expression scores or existing gene signatures, allowing for reliable identification of responders even within the desert patient subcohort analyzed as a hold out set.

The organizational integrity of the adaptive immune system is determined by functionally discrete subsets of CD4+ T cells, but it has remained unclear to what extent lineage choice is influenced by clonotypically expressed T-cell receptors (TCRs). To address this issue, we used a high-throughput approach to profile the αβ TCR repertoires of human naive and effector/memory CD4+ T-cell subsets, irrespective of antigen specificity. Highly conserved physicochemical and recombinatorial features were encoded on a subset-specific basis in the effector/memory compartment. Clonal tracking further identified forbidden and permitted transition pathways, mapping effector/memory subsets related by interconversion or ontogeny. Public sequences were largely confined to particular effector/memory subsets, including regulatory T cells (Tregs), which also displayed hardwired repertoire features in the naive compartment. Accordingly, these cumulative repertoire portraits establish a link between clonotype fate decisions in the complex world of CD4+ T cells and the intrinsic properties of somatically rearranged TCRs.

Comparative genomics analysis of conserved gene cassettes demonstrated resemblance between a recently described cassette of genes involved in sulphoquinovose degradation in Escherichia coli K-12 MG1655 and a Bacilli cassette linked with lactose degradation. Six genes from both cassettes had similar functions related to carbohydrate metabolism, namely, hydrolase, aldolase, kinase, isomerase, transporter, and transcription factor. The Escherichia coli sulphoglycolysis cassette was thus predicted to be associated with lactose degradation. This prediction was confirmed experimentally: expression of genes coding for aldolase (yihT), isomerase (yihS), and kinase (yihV) was dramatically increased during growth on lactose. These genes were previously shown to be activated during growth on sulphoquinovose, so our observation may indicate multi-functional capabilities of the respective proteins. Transcription starts for yihT, yihV and yihW were mapped in silico, in vitro and in vivo. Out of three promoters for yihT, one was active only during growth on lactose. We further showed that switches in yihT transcription are controlled by YihW, a DeoR-family transcription factor in the Escherichia coli cassette. YihW acted as a carbon source-dependent dual regulator involved in sustaining the baseline growth in the absence of lac-operon, with function either complementary, or opposite to a global regulator of carbohydrate metabolism, cAMP-CRP.

Background Bacterial carbohydrate metabolism is extremely diverse, since carbohydrates serve as a major energy source and are involved in a variety of cellular processes. Bacterial genes belonging to same metabolic pathway are often co-localized in the chromosome, but it is not a strict rule. Gene co-localization in linked to co-evolution and co-regulation. This study focuses on a large-scale analysis of bacterial genomic loci related to the carbohydrate metabolism. Results We demonstrate that only 53% of 148,000 studied genes from over six hundred bacterial genomes are co-localized in bacterial genomes with other carbohydrate metabolism genes, which points to a significant role of singleton genes. Co-localized genes form cassettes, ranging in size from two to fifteen genes. Two major factors influencing the cassette-forming tendency are gene function and bacterial phylogeny. We have obtained a comprehensive picture of co-localization preferences of genes for nineteen major carbohydrate metabolism functional classes, over two hundred gene orthologous clusters, and thirty bacterial classes, and characterized the cassette variety in size and content among different species, highlighting a significant role of short cassettes. The preference towards co-localization of carbohydrate metabolism genes varies between 40 and 76% for bacterial taxa. Analysis of frequently co-localized genes yielded forty-five significant pairwise links between genes belonging to different functional classes. The number of such links per class range from zero to eight, demonstrating varying preferences of respective genes towards a specific chromosomal neighborhood. Genes from eleven functional classes tend to co-localize with genes from the same class, indicating an important role of clustering of genes with similar functions. At that, in most cases such co-localization does not originate from local duplication events. Conclusions Overall, we describe a complex web formed by evolutionary relationships of bacterial carbohydrate metabolism genes, manifested as co-localization patterns. Reviewers This article was reviewed by Daria V. Dibrova (A.N. Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Moscow, Russia), nominated by Armen Mulkidjanian (University of Osnabrück, Germany), Igor Rogozin (NCBI, NLM, NIH, USA) and Yuri Wolf (NCBI, NLM, NIH, USA). Electronic supplementary material The online version of this article (10.1186/s13062-017-0200-7) contains supplementary material, which is available to authorized users.

BackgroundGenome sequencing of multiple strains demonstrated high variability in gene content even in closely related strains of the same species and created a newly emerged object for genomic analysis, the pan-genome, that is, the complete set of genes observed in a given species or a higher level taxon. Here we analysed the pan-genome structure and the genome evolution of 25 strains of Streptococcus suis, 50 strains of Streptococcus pyogenes and 28 strains of Streptococcus pneumoniae.ResultsFractions of the pan-genome, unique, periphery, and universal genes differ in size, functional composition, the level of nucleotide substitutions, and predisposition to horizontal gene transfer and genomic rearrangements. The density of substitutions in intergenic regions appears to be correlated with selection acting on adjacent genes, implying that more conserved genes tend to have more conserved regulatory regions. The total pan-genome of the genus is open, but only due to strain-specific genes, whereas other pan-genome fractions reach saturation. The strain-specific fraction is enriched with mobile elements and hypothetical proteins, but also contains a number of candidate virulence-related genes, so it may have a strong impact on adaptability and pathogenicity.About 7% of single-copy periphery genes have been found in different syntenic regions. More than a half of these genes are rare in all Streptococcus species; others are rare in at least one species. We have identified the set of genes with phylogenies inconsistent with species and non-conserved location in the chromosome; these genes are candidates for horizontal transfer between species.An inversion of length 15 kB found in four independent branches of S. pneumoniae has breakpoints formed by genes encoding a surface antigen protein (PhtD). The observed parallelism may indicate the action of an antigen variation mechanism.ConclusionsMembers of the genus Streptococcus have a highly dynamic, open pan-genome, that potentially confers them with the ability to adapt to changing environmental conditions, i.e. antibiotic resistance or transmission between different hosts. Hence, understanding of genome evolution is important for the identification of potential pathogens and design of drugs and vaccines.

Corals harbor complex and diverse microbial communities that strongly impact host fitness and resistance to diseases, but these microbes themselves can be influenced by stresses, like those caused by the presence of macroscopic symbionts. In addition to directly influencing the host, symbionts may transmit pathogenic microbial communities. We analyzed two coral gall-forming copepod systems by using 16S rRNA gene metagenomic sequencing: (1) the sea fan Gorgonia ventalina with copepods of the genus Sphaerippe from the Caribbean and (2) the scleractinian coral Stylophora pistillata with copepods of the genus Spaniomolgus from the Saudi Arabian part of the Red Sea. We show that bacterial communities in these two systems were substantially different with Actinobacteria, Alphaproteobacteria, and Betaproteobacteria more prevalent in samples from Gorgonia ventalina, and Gammaproteobacteria in Stylophora pistillata. In Stylophora pistillata, normal coral microbiomes were enriched with the common coral symbiont Endozoicomonas and some unclassified bacteria, while copepod and gall-tissue microbiomes were highly enriched with the family ME2 (Oceanospirillales) or Rhodobacteraceae. In Gorgonia ventalina, no bacterial group had significantly different prevalence in the normal coral tissues, copepods, and injured tissues. The total microbiome composition of polyps injured by copepods was different. Contrary to our expectations, the microbial community composition of the injured gall tissues was not directly affected by the microbiome of the gall-forming symbiont copepods.

Clostridium difficile, a Gram positive, anaerobic, spore-forming bacterium is an emergent pathogen and the most common cause of nosocomial diarrhea. Although transmission of C. difficile is mediated by contamination of the gut by spores, the regulatory cascade controlling spore formation remains poorly characterized. During Bacillus subtilis sporulation, a cascade of four sigma factors, σF and σG in the forespore and σE and σK in the mother cell governs compartment-specific gene expression. In this work, we combined genome wide transcriptional analyses and promoter mapping to define the C. difficile σF, σE, σG and σK regulons. We identified about 225 genes under the control of these sigma factors: 25 in the σF regulon, 97 σE-dependent genes, 50 σG-governed genes and 56 genes under σK control. A significant fraction of genes in each regulon is of unknown function but new candidates for spore coat proteins could be proposed as being synthesized under σE or σK control and detected in a previously published spore proteome. SpoIIID of C. difficile also plays a pivotal role in the mother cell line of expression repressing the transcription of many members of the σE regulon and activating sigK expression. Global analysis of developmental gene expression under the control of these sigma factors revealed deviations from the B. subtilis model regarding the communication between mother cell and forespore in C. difficile. We showed that the expression of the σE regulon in the mother cell was not strictly under the control of σF despite the fact that the forespore product SpoIIR was required for the processing of pro-σE. In addition, the σK regulon was not controlled by σG in C. difficile in agreement with the lack of pro-σK processing. This work is one key step to obtain new insights about the diversity and evolution of the sporulation process among Firmicutes., Author Summary Clostridium difficile, a major cause of antibiotic-associated diarrhea, produces resistant spores that facilitate the persistence of this bacterium in the environment including hospitals. Its transmission is mediated by contamination of gut by spores. Understanding how this complex developmental process is regulated is fundamental to decipher C. difficile transmission and pathogenesis. The regulatory cascade controlling sporulation that involves four sigma factors, σF and σG in the forespore and σE and σK in the mother cell remains poorly characterized in C. difficile. By combining transcriptome analysis and promoter mapping, we identified genes expressed under the specific control of each sigma factor. Among sporulation-controlled proteins detected in spore, we can propose candidates for new spore coat proteins important for spore resistance. We also showed differences in the intercompartment communication between forespore and mother cell in C. difficile compared to the Bacillus subtilis model. In C. difficile, we observed that the activation of the σE regulon was partially independent of σF and that the σK regulon was not controlled by σG. Our finding suggests that the C. difficile sporulation process might be more ancestral compared to that of B. subtilis.