HLA-DR Expression in Natural Killer Cells Marks Distinct Functional States, Depending on Cell Differentiation Stage.

HLA-DR-positive NK cells, found in both healthy individuals and patients with different inflammatory diseases, are characterized as activated cells. However, data on their capacity for IFNγ production or cytotoxic response vary between studies. Thus, more precise investigation is needed of the mechanisms related to the induction of HLA-DR expression in NK cells, their associations with NK cell differentiation stage, and functional or metabolic state. In this work, HLA-DR-expressing NK cell subsets were investigated using transcriptomic analysis, metabolic activity assays, and analysis of intercellular signaling cascades. We demonstrated that HLA-DR+CD56bright NK cells were characterized by a proliferative phenotype, while HLA-DR+CD56dim NK cells exhibited features of adaptive cells and loss of inhibitory receptors with increased expression of MHC class II trans-activator CIITA. The activated state of HLA-DR-expressing NK cells was confirmed by higher levels of ATP and mitochondrial mass observed in this subset compared to HLA-DR- cells, both ex vivo and after stimulation in culture. We showed that HLA-DR expression in NK cells in vitro can be induced both through stimulation by exogenous IL-2 and IL-21, as well as through auto-stimulation by NK-cell-produced IFNγ. At the intracellular level, HLA-DR expression depended on the activation of STAT3- and ERK1/2-mediated pathways, with subsequent activation of isoform 3 of the transcription factor CIITA. The obtained results broaden the knowledge about HLA-DR-positive NK cell appearance, diversity, and functions, which might be useful in terms of understanding the role of this subset in innate immunity and assessing their possible implications in NK cell-based therapy.

Aberrant adaptive immune response underlies genetic susceptibility to tuberculosis.

Mycobacterium tuberculosis (Mtb) remains a major threat worldwide, although only a fraction of infected individuals develops tuberculosis (TB). TB susceptibility is shaped by multiple genetic factors, and we performed comparative immunological analysis of two mouse strains to uncover relevant mechanisms underlying susceptibility and resistance. C57BL/6 mice are relatively TB-resistant, whereas I/St mice are prone to develop severe TB, partly due to the MHC-II allelic variant that shapes suboptimal CD4+ T cell receptor repertoire. We investigated the repertoires of lung-infiltrating helper T cells and B cells at the progressed stage in both strains. We found that lung CD4+ T cell repertoires of infected C57BL/6 but not I/St mice contained convergent TCR clusters with functionally confirmed Mtb specificity. Transcriptomic analysis revealed a more prominent Th1 signature in C57BL/6, and expression of pro-inflammatory IL-16 in I/St lung-infiltrating helper T cells. The two strains also showed distinct Th2 signatures. Furthermore, the humoral response of I/St mice was delayed, less focused, and dominated by IgG/IgM isotypes, whereas C57BL/6 mice generated more Mtb antigen-focused IgA response. We conclude that the inability of I/St mice to produce a timely and efficient anti-Mtb adaptive immune responses arises from a suboptimal helper T cell landscape that also impacts the humoral response, leading to diffuse inflammation and severe disease.

Targeted checkpoint control of B cells undergoing positive selection in germinal centers by follicular regulatory T cells.

Follicular regulatory T cells (Tfr) can play opposite roles in the regulation of germinal center (GC) responses. Depending on the studies, Tfr suppress or support GC and B cell affinity maturation. However, which factors determine positive vs. negative effects of Tfr on the GC B cell is unclear. In this study, we show that GC centrocytes that express MYC up-regulate expression of CCL3 chemokine that is needed for both the positive and negative regulation of GC B cells by Tfr. B cell-intrinsic expression of CCL3 contributes to Tfr-dependent positive selection of foreign Ag-specific GC B cells. At the same time, expression of CCL3 is critical for direct Tfr-mediated suppression of GC B cells that acquire cognate to Tfr nuclear proteins. Our study suggests that CCR5 and CCR1 receptors promote Tfr migration to CCL3 and highlights Ccr5 expression on the Tfr subset that expresses Il10. Based on our findings and previous studies, we suggest a model of chemotactically targeted checkpoint control of B cells undergoing positive selection in GCs by Tfr, where Tfr directly probe and license foreign antigen-specific B cells to complete their positive selection in GCs but, at the same time, suppress GC B cells that present self-antigens cognate to Tfr.

Data on the temporal changes in soil properties and microbiome composition after a jet-fuel contamination during the pot and field experiments.

The soil response to a jet-fuel contamination is uncertain. In this article, original data on the influence of a jet-fuel spillage on the topsoil properties are presented. The data set is obtained during a one-year long pot and field experiments with Dystric Arenosols, Fibric Histosols and Albic Luvisols. Kerosene loads were 1, 5, 10, 25 and 100 g/kg. The data set includes information about temporal changes in kerosene concentration; physicochemical properties, such as рН, moisture, cation exchange capacity, content of soil organic matter, available P and K, exchangeable NH4 +, and water-soluble NO3 -; and biological properties, such as biological consumption of oxygen, and cellulolytic activity. Also, we provide sequencing data on variable regions of 16S ribosomal RNA of microbial communities from the respective soil samples.

The Influence of Kerosene on Microbiomes of Diverse Soils.

One of the most important challenges for soil science is to determine the limits for the sustainable functioning of contaminated ecosystems. The response of soil microbiomes to kerosene pollution is still poorly understood. Here, we model the impact of kerosene leakage on the composition of the topsoil microbiome in pot and field experiments with different loads of added kerosene (loads up to 100 g/kg; retention time up to 360 days). At four time points we measured kerosene concentration and sequenced variable regions of 16S ribosomal RNA in the microbial communities. Mainly alkaline Dystric Arenosols with low content of available phosphorus and soil organic matter had an increased fraction of Actinobacteriota, Firmicutes, Nitrospirota, Planctomycetota, and, to a lesser extent, Acidobacteriota and Verrucomicobacteriota. In contrast, in highly acidic Fibric Histosols, rich in soil organic matter and available phosphorus, the fraction of Acidobacteriota was higher, while the fraction of Actinobacteriota was lower. Albic Luvisols occupied an intermediate position in terms of both physicochemical properties and microbiome composition. The microbiomes of different soils show similar response to equal kerosene loads. In highly contaminated soils, the proportion of anaerobic bacteria-metabolizing hydrocarbons increased, whereas the proportion of aerobic bacteria decreased. During the field experiment, the soil microbiome recovered much faster than in the pot experiments, possibly due to migration of microorganisms from the polluted area. The microbial community of Fibric Histosols recovered in 6 months after kerosene had been loaded, while microbiomes of Dystric Arenosols and Albic Luvisols did not restore even after a year.

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.

Naïve Regulatory T Cell Subset Is Altered in X-Linked Agammaglobulinemia.

The interplay between T- and B-cell compartments during naïve, effector and memory T cell maturation is critical for a balanced immune response. Primary B-cell immunodeficiency arising from X-linked agammaglobulinemia (XLA) offers a model to explore B cell impact on T cell subsets, starting from the thymic selection. Here we investigated characteristics of naïve and effector T cell subsets in XLA patients, revealing prominent alterations in the corresponding T-cell receptor (TCR) repertoires. We observed immunosenescence in terms of decreased diversity of naïve CD4+ and CD8+ TCR repertoires in XLA donors. The most substantial alterations were found within naïve CD4+ subsets, and we have investigated these in greater detail. In particular, increased clonality and convergence, along with shorter CDR3 regions, suggested narrower focused antigen-specific maturation of thymus-derived naïve Treg (CD4+CD45RA+CD27+CD25+) in the absence of B cells - normally presenting diverse self and commensal antigens. The naïve Treg proportion among naïve CD4 T cells was decreased in XLA patients, supporting the concept of impaired thymic naïve Treg selection. Furthermore, the naïve Treg subset showed prominent differences at the transcriptome level, including increased expression of genes specific for antigen-presenting and myeloid cells. Altogether, our findings suggest active B cell involvement in CD4 T cell subsets maturation, including B cell-dependent expansion of the naïve Treg TCR repertoire that enables better control of self-reactive T cells.

Functionally specialized human CD4+ T-cell subsets express physicochemically distinct TCRs.

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.

Genome-Wide Transcription Start Site Mapping and Promoter Assignments to a Sigma Factor in the Human Enteropathogen Clostridioides difficile.

The emerging human enteropathogen Clostridioides difficile is the main cause of diarrhea associated with antibiotherapy. Regulatory pathways underlying the adaptive responses remain understudied and the global view of C. difficile promoter structure is still missing. In the genome of C. difficile 630, 22 genes encoding sigma factors are present suggesting a complex pattern of transcription in this bacterium. We present here the first transcriptional map of the C. difficile genome resulting from the identification of transcriptional start sites (TSS), promoter motifs and operon structures. By 5'-end RNA-seq approach, we mapped more than 1000 TSS upstream of genes. In addition to these primary TSS, this analysis revealed complex structure of transcriptional units such as alternative and internal promoters, potential RNA processing events and 5' untranslated regions. By following an in silico iterative strategy that used as an input previously published consensus sequences and transcriptomic analysis, we identified candidate promoters upstream of most of protein-coding and non-coding RNAs genes. This strategy also led to refine consensus sequences of promoters recognized by major sigma factors of C. difficile. Detailed analysis focuses on the transcription in the pathogenicity locus and regulatory genes, as well as regulons of transition phase and sporulation sigma factors as important components of C. difficile regulatory network governing toxin gene expression and spore formation. Among the still uncharacterized regulons of the major sigma factors of C. difficile, we defined the SigL regulon by combining transcriptome and in silico analyses. We showed that the SigL regulon is largely involved in amino-acid degradation, a metabolism crucial for C. difficile gut colonization. Finally, we combined our TSS mapping, in silico identification of promoters and RNA-seq data to improve gene annotation and to suggest operon organization in C. difficile. These data will considerably improve our knowledge of global regulatory circuits controlling gene expression in C. difficile and will serve as a useful rich resource for scientific community both for the detailed analysis of specific genes and systems biology studies.

MHC-II alleles shape the CDR3 repertoires of conventional and regulatory naïve CD4+ T cells.

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 (Tconv) and naïve regulatory CD4+ T (Treg) cells. Compared with tuberculosis-resistant C57BL/6 (H2-Ab) mice, the tuberculosis-susceptible H2-Aj mice had fewer CD4+ T cells of both subsets in the thymus. In the periphery, this deficiency was only apparent for Tconv and was compensated for by peripheral reconstitution for Treg We show that H2-Aj 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-Aj and H2-Ab 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.

Chlamydia pan-genomic analysis reveals balance between host adaptation and selective pressure to genome reduction.

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.

Micro-evolution of three Streptococcus species: selection, antigenic variation, and horizontal gene inflow.

BACKGROUND: The genus Streptococcus comprises pathogens that strongly influence the health of humans and animals. Genome sequencing of multiple Streptococcus strains demonstrated high variability in gene content and order even in closely related strains of the same species and created a newly emerged object for genomic analysis, the pan-genome. Here we analysed the genome evolution of 25 strains of Streptococcus suis, 50 strains of Streptococcus pyogenes and 28 strains of Streptococcus pneumoniae. RESULTS: Fractions 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. We have identified the set of genes with phylogenies inconsistent with species and non-conserved location in the chromosome; these genes are rare in at least one species and have likely experienced recent horizontal transfer between species. 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. Mapping the rearrangements to the phylogenetic tree revealed large parallel inversions in all species. A parallel inversion of length 15 kB with breakpoints formed by genes encoding surface antigen proteins PhtD and PhtB in S. pneumoniae leads to replacement of gene fragments that likely indicates the action of an antigen variation mechanism. CONCLUSIONS: Members of 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, integrated analysis of all aspects of genome evolution is important for the identification of potential pathogens and design of drugs and vaccines.

Microbiomes of gall-inducing copepod crustaceans from the corals Stylophora pistillata (Scleractinia) and Gorgonia ventalina (Alcyonacea).

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.

The genes of the sulphoquinovose catabolism in Escherichia coli are also associated with a previously unknown pathway of lactose degradation.

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.

Sugar Lego: gene composition of bacterial carbohydrate metabolism genomic loci.

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).

Genome-wide analysis of cell type-specific gene transcription during spore formation in Clostridium difficile.

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.