BCG vaccination alters the epigenetic landscape of progenitor cells in human bone marrow to influence innate immune responses.

Although the Bacille-Calmette-Guérin (BCG) vaccine is used to prevent tuberculosis, it also offers protection against a diverse range of non-mycobacterial infections. However, the underlying protective mechanisms in humans are not yet fully understood. Here, we surveyed at single-cell resolution the gene expression and chromatin landscape of human bone marrow, aspirated before and 90 days after BCG vaccination or placebo. We showed that BCG alters both the gene expression and epigenetic profiles of human hematopoietic stem and progenitor cells (HSPCs). Changes in gene expression occurred primarily within uncommitted stem cells. By contrast, changes in chromatin accessibility were most prevalent within differentiated progenitor cells at sites influenced by Kruppel-like factor (KLF) and early growth response (EGR) transcription factors and were highly correlated (r > 0.8) with the interleukin (IL)-1ß secretion capacity of paired peripheral blood mononuclear cells (PBMCs). Our findings shed light on BCG vaccination's profound and lasting effects on HSPCs and its influence on innate immune responses and trained immunity.

Model-based impact evaluation of new tuberculosis vaccines in aging populations under different modeling scenarios: the case of China.

Introduction: The slow descent in TB burden, the COVID-19 pandemic, along with the rise of multidrug-resistant strains of Mycobacterium tuberculosis, seriously threaten TB control and the goals of the End TB strategy. To fight back, several vaccine candidates are under development, with some of them undergoing the phases 2B and 3 of the development pipeline. The impact of these vaccines on the general population needs to be addressed using disease-transmission models, and, in a country like China, which last year ranked third in number of cases worldwide, and where the population is aging at a fast pace, the impact of TB vaccination campaigns may depend heavily upon the age of targeted populations, the mechanistic descriptions of the TB vaccines and the coupling between TB dynamics and demographic evolution. Methods: In this work, we studied the potential impact of a new TB vaccine in China targeting adolescents (15-19 y.o.) or older adults (60-64 y.o.), according to varying vaccine descriptions that represent reasonable mechanisms of action leading to prevention of disease, or prevention of recurrence, each of them targetting specific routes to TB disease. To measure the influence of the description of the coupling between transmission dynamics and aging in TB transmission models, we explored two different approaches to compute the evolution of the contact matrices, which relate to the spreading among different age strata. Results: Our findings highlight the dependence of model-based impact estimates on vaccine profiles and the chosen modeling approach for describing the evolution of contact matrices. Our results also show, in line with previous modeling works, that older adult vaccination is a suitable option in China to reduce the incidence of TB as long as the vaccine is able to protect already exposed individuals. Discussion: This study underscores the importance of considering vaccine characteristics and demographic dynamics in shaping TB control strategies. In this sense, older adult vaccination emerges as a promising avenue for mitigating TB transmission in China but also remarks the need for tailored intervention strategies aligned with demographic trends.

Maintenance of cell wall remodeling and vesicle production are connected in Mycobacterium tuberculosis.

Pathogenic and nonpathogenic mycobacteria secrete extracellular vesicles (EVs) under various conditions. EVs produced by Mycobacterium tuberculosis ( Mtb ) have raised significant interest for their potential in cell communication, nutrient acquisition, and immune evasion. However, the relevance of vesicle secretion during tuberculosis infection remains unknown due to the limited understanding of mycobacterial vesicle biogenesis. We have previously shown that a transposon mutant in the LCP-related gene virR ( virR mut ) manifested a strong attenuated phenotype during experimental macrophage and murine infections, concomitant to enhanced vesicle release. In this study, we aimed to understand the role of VirR in the vesicle production process in Mtb . We employ genetic, transcriptional, proteomics, ultrastructural and biochemical methods to investigate the underlying processes explaining the enhanced vesiculogenesis phenomenon observed in the virR mut . Our results establish that VirR is critical to sustain proper cell permeability via regulation of cell envelope remodeling possibly through the interaction with similar cell envelope proteins, which control the link between peptidoglycan and arabinogalactan. These findings advance our understanding of mycobacterial extracellular vesicle biogenesis and suggest that these set of proteins could be attractive targets for therapeutic intervention.

BCG vaccination impacts the epigenetic landscape of progenitor cells in human bone marrow.

While the Bacille-Calmette-Guérin (BCG) vaccine is used to prevent tuberculosis, it also offers protection against a diverse range of non-mycobacterial infections. However, the underlying protective mechanisms in humans are not yet fully understood. Here, we surveyed at single-cell resolution the gene expression and chromatin landscape of human bone marrow, aspirated before and 90 days after BCG vaccination or placebo administration. We show that BCG vaccination significantly alters both the gene expression and epigenetic profiles of human hematopoietic stem and progenitor cells (HSPCs). Changes in gene expression occur primarily on the most uncommitted stem cells and are reflective of a persistent myeloid bias. In contrast, BCG-induced changes in chromatin accessibility are most prevalent within differentiated progenitor cells at sites influenced by Kruppel-like factor (KLF)/SP and EGR transcription factors (TFs). These TFs are also activated in the most uncommitted stem cells, indicating that activated TFs, which drive persistent changes in HSC gene expression, likely also drive chromatin dynamics appearing within downstream progenitor cells. This perspective contests the prevailing notion that epigenetic modifications linked to innate immune memory transfer directly from stem cells to their differentiated derivatives. Finally, we show that alterations in gene expression and chromatin accessibility in HSPCs due to BCG vaccination were highly correlated (r>0.8) with the IL-1ß secretion capacity of paired PBMCs upon secondary immune challenge. Overall, our findings shed light on BCG vaccination's profound and lasting effects on HSPCs and its influence on innate immune responses.

Mycobacterium tuberculosis infection drives differential responses in the bone marrow hematopoietic stem and progenitor cells.

Hematopoietic stem and progenitor cells (HSPCs) play a vital role in the host response to infection through the rapid and robust production of mature immune cells. These HSPC responses can be influenced, directly and indirectly, by pathogens as well. Infection with Mycobacterium tuberculosis (Mtb) can drive lymphopoiesis through modulation of type I interferon (IFN) signaling. We have previously found that the presence of a drug resistance (DR)-conferring mutation in Mtb drives altered host-pathogen interactions and heightened type I IFN production in vitro. But the impacts of this DR mutation on in vivo host responses to Mtb infection, particularly the hematopoietic compartment, remain unexplored. Using a mouse model, we show that, while drug-sensitive Mtb infection induces expansion of HSPC subsets and a skew toward lymphopoiesis, DR Mtb infection fails to induce an expansion of these subsets and an accumulation of mature granulocytes in the bone marrow. Using single-cell RNA sequencing, we show that the HSCs from DR Mtb-infected mice fail to upregulate pathways related to cytokine signaling across all profiled HSC subsets. Collectively, our studies report a novel finding of a chronic infection that fails to induce a potent hematopoietic response that can be further investigated to understand pathogen-host interaction at the level of hematopoiesis.

Corrigendum: Iron deprivation enhances transcriptional responses to in vitro growth arrest of Mycobacterium tuberculosis.

[This corrects the article DOI: 10.3389/fmicb.2022.956602.].

Iron deprivation enhances transcriptional responses to in vitro growth arrest of Mycobacterium tuberculosis.

The establishment of Mycobacterium tuberculosis (Mtb) long-term infection in vivo depends on several factors, one of which is the availability of key nutrients such as iron (Fe). The relation between Fe deprivation inside and outside the granuloma, and the capacity of Mtb to accumulate lipids and persist in the absence of growth is not well understood. In this context, current knowledge of how Mtb modifies its lipid composition in response to growth arrest, depending on iron availability, is scarce. To shed light on these matters, in this work we compare genome-wide transcriptomic and lipidomic profiles of Mtb at exponential and stationary growth phases using cultures with glycerol as a carbon source, in the presence or absence of iron. As a result, we found that transcriptomic responses to growth arrest, considered as the transition from exponential to stationary phase, are iron dependent for as many as 714 genes (iron-growth interaction contrast, FDR <0.05), and that, in a majority of these genes, iron deprivation enhances the magnitude of the transcriptional responses to growth arrest in either direction. On the one hand, genes whose upregulation upon growth arrest is enhanced by iron deprivation were enriched in functional terms related to homeostasis of ion metals, and responses to several stressful cues considered cardinal features of the intracellular environment. On the other hand, genes showing negative responses to growth arrest that are stronger in iron-poor medium were enriched in energy production processes (TCA cycle, NADH dehydrogenation and cellular respiration), and key controllers of ribosomal activity shut-down, such as the T/A system mazE6/F6. Despite of these findings, a main component of the cell envelope, lipid phthiocerol dimycocerosate (PDIM), was not detected in the stationary phase regardless of iron availability, suggesting that lipid changes during Mtb adaptation to non-dividing phenotypes appear to be iron-independent. Taken together, our results indicate that environmental iron levels act as a key modulator of the intensity of the transcriptional adaptations that take place in the bacterium upon its transition between dividing and dormant-like phenotypes in vitro.

Modeling the impact of COVID-19 on future tuberculosis burden.

Background: The ongoing COVID-19 pandemic has greatly disrupted our everyday life, forcing the adoption of non-pharmaceutical interventions in many countries and putting public health services and healthcare systems worldwide under stress. These circumstances are leading to unintended effects such as the increase in the burden of other diseases. Methods: Here, using a data-driven epidemiological model for tuberculosis (TB) spreading, we describe the expected rise in TB incidence and mortality if COVID-associated changes in TB notification are sustained and attributable entirely to disrupted diagnosis and treatment adherence. Results: Our calculations show that the reduction in diagnosis of new TB cases due to the COVID-19 pandemic could result in 228k (CI 187-276) excess deaths in India, 111k (CI 93-134) in Indonesia, 27k (CI 21-33) in Pakistan, and 12k (CI 9-18) in Kenya. Conclusions: We show that it is possible to reverse these excess deaths by increasing the pre-covid diagnosis capabilities from 15 to 50% for 2 to 4 years. This would prevent almost all TB-related excess mortality that could be caused by the COVID-19 pandemic if no additional preventative measures are introduced. Our work therefore provides guidelines for mitigating the impact of COVID-19 on tuberculosis epidemic in the years to come.

The COVID-19 pandemic has disrupted everyday life and put public health services and healthcare systems worldwide under stress. This has compromised the ability to control other diseases such as Malaria, Cancer and Tuberculosis. In this work we predict the rise in Tuberculosis occurrence and mortality when healthcare systems are impacted and diagnosis capabilities blocked in 4 countries where TB is prevalent. Our calculations show that an increase in new TB cases due to the COVID-19 pandemic could result in almost 400,000 additional deaths from TB in India, Indonesia, Pakistan and Kenya. We also show that increased diagnosis capabilities after the pandemic could reduce the additional deaths from TB resulting from the COVID-19 pandemic impact.

Human alveolar macrophage metabolism is compromised during Mycobacterium tuberculosis infection.

Pulmonary macrophages have two distinct ontogenies: long-lived embryonically-seeded alveolar macrophages (AM) and bone marrow-derived macrophages (BMDM). Here, we show that after infection with a virulent strain of Mycobacterium tuberculosis (H37Rv), primary murine AM exhibit a unique transcriptomic signature characterized by metabolic reprogramming distinct from conventional BMDM. In contrast to BMDM, AM failed to shift from oxidative phosphorylation (OXPHOS) to glycolysis and consequently were unable to control infection with an avirulent strain (H37Ra). Importantly, healthy human AM infected with H37Ra equally demonstrated diminished energetics, recapitulating our observation in the murine model system. However, the results from seahorse showed that the shift towards glycolysis in both AM and BMDM was inhibited by H37Rv. We further demonstrated that pharmacological (e.g. metformin or the iron chelator desferrioxamine) reprogramming of AM towards glycolysis reduced necrosis and enhanced AM capacity to control H37Rv growth. Together, our results indicate that the unique bioenergetics of AM renders these cells a perfect target for Mtb survival and that metabolic reprogramming may be a viable host targeted therapy against TB.

Alveolar macrophages from persons living with HIV show impaired epigenetic response to Mycobacterium tuberculosis.

Persons living with HIV (PLWH) are at increased risk of tuberculosis (TB). HIV-associated TB is often the result of recent infection with Mycobacterium tuberculosis (M. tuberculosis) followed by rapid progression to disease. Alveolar macrophages (AMs) are the first cells of the innate immune system that engage M. tuberculosis, but how HIV and antiretroviral therapy (ART) affect the anti-mycobacterial response of AMs is not known. To investigate the impact of HIV and ART on the transcriptomic and epigenetic response of AMs to M. tuberculosis, we obtained AMs by bronchoalveolar lavage from 20 PLWH receiving ART, 16 control subjects who were HIV-free (HC), and 14 subjects who received ART as preexposure prophylaxis (PrEP) to prevent HIV infection. Following in vitro challenge with M. tuberculosis, AMs from each group displayed overlapping but distinct profiles of significantly up- and downregulated genes in response to M. tuberculosis. Comparatively, AMs isolated from both PLWH and PrEP subjects presented a substantially weaker transcriptional response. In addition, AMs from HC subjects challenged with M. tuberculosis responded with pronounced chromatin accessibility changes while AMs obtained from PLWH and PrEP subjects displayed no significant changes in their chromatin state. Collectively, these results revealed a stronger adverse effect of ART than HIV on the epigenetic landscape and transcriptional responsiveness of AMs.

M. tuberculosis Reprograms Hematopoietic Stem Cells to Limit Myelopoiesis and Impair Trained Immunity.

A greater understanding of hematopoietic stem cell (HSC) regulation is required for dissecting protective versus detrimental immunity to pathogens that cause chronic infections such as Mycobacterium tuberculosis (Mtb). We have shown that systemic administration of Bacille Calmette-Guérin (BCG) or ß-glucan reprograms HSCs in the bone marrow (BM) via a type II interferon (IFN-II) or interleukin-1 (IL1) response, respectively, which confers protective trained immunity against Mtb. Here, we demonstrate that, unlike BCG or ß-glucan, Mtb reprograms HSCs via an IFN-I response that suppresses myelopoiesis and impairs development of protective trained immunity to Mtb. Mechanistically, IFN-I signaling dysregulates iron metabolism, depolarizes mitochondrial membrane potential, and induces cell death specifically in myeloid progenitors. Additionally, activation of the IFN-I/iron axis in HSCs impairs trained immunity to Mtb infection. These results identify an unanticipated immune evasion strategy of Mtb in the BM that controls the magnitude and intrinsic anti-microbial capacity of innate immunity to infection.

Social history and exposure to pathogen signals modulate social status effects on gene regulation in rhesus macaques.

Social experience is an important predictor of disease susceptibility and survival in humans and other social mammals. Chronic social stress is thought to generate a proinflammatory state characterized by elevated antibacterial defenses and reduced investment in antiviral defense. Here we manipulated long-term social status in female rhesus macaques to show that social subordination alters the gene expression response to ex vivo bacterial and viral challenge. As predicted by current models, bacterial lipopolysaccharide polarizes the immune response such that low status corresponds to higher expression of genes in NF-κB-dependent proinflammatory pathways and lower expression of genes involved in the antiviral response and type I IFN signaling. Counter to predictions, however, low status drives more exaggerated expression of both NF-κB- and IFN-associated genes after cells are exposed to the viral mimic Gardiquimod. Status-driven gene expression patterns are linked not only to social status at the time of sampling, but also to social history (i.e., past social status), especially in unstimulated cells. However, for a subset of genes, we observed interaction effects in which females who fell in rank were more strongly affected by current social status than those who climbed the social hierarchy. Taken together, our results indicate that the effects of social status on immune cell gene expression depend on pathogen exposure, pathogen type, and social history-in support of social experience-mediated biological embedding in adulthood, even in the conventionally memory-less innate immune system.

Bridging the gap between efficacy trials and model-based impact evaluation for new tuberculosis vaccines.

In Tuberculosis (TB), given the complexity of its transmission dynamics, observations of reduced epidemiological risk associated with preventive interventions can be difficult to translate into mechanistic interpretations. Specifically, in clinical trials of vaccine efficacy, a readout of protection against TB disease can be mapped to multiple dynamical mechanisms, an issue that has been overlooked so far. Here, we describe this limitation and its effect on model-based evaluations of vaccine impact. Furthermore, we propose a methodology to analyze efficacy trials that circumvents it, leveraging a combination of compartmental models and stochastic simulations. Using our approach, we can disentangle the different possible mechanisms of action underlying vaccine protection effects against TB, conditioned to trial design, size, and duration. Our results unlock a deeper interpretation of the data emanating from efficacy trials of TB vaccines, which renders them more interpretable in terms of transmission models and translates into explicit recommendations for vaccine developers.

Natural selection contributed to immunological differences between hunter-gatherers and agriculturalists.

The shift from a hunter-gatherer to an agricultural mode of subsistence is believed to have been associated with profound changes in the burden and diversity of pathogens across human populations. Yet, the extent to which the advent of agriculture affected the evolution of the human immune system remains unknown. Here we present a comparative study of variation in the transcriptional responses of peripheral blood mononuclear cells to bacterial and viral stimuli between Batwa rainforest hunter-gatherers and Bakiga agriculturalists from Uganda. We observed increased divergence between hunter-gatherers and agriculturalists in the early transcriptional response to viruses compared with that for bacterial stimuli. We demonstrate that a significant fraction of these transcriptional differences are under genetic control and we show that positive natural selection has helped to shape population differences in immune regulation. Across the set of genetic variants underlying inter-population immune-response differences, however, the signatures of positive selection were disproportionately observed in the rainforest hunter-gatherers. This result is counter to expectations on the basis of the popularized notion that shifts in pathogen exposure due to the advent of agriculture imposed radically heightened selective pressures in agriculturalist populations.

Efficient and Robust NK-Cell Transduction With Baboon Envelope Pseudotyped Lentivector.

NK-cell resistance to transduction is a major technical hurdle for developing NK-cell immunotherapy. By using Baboon envelope pseudotyped lentiviral vectors (BaEV-LVs) encoding eGFP, we obtained a transduction rate of 23.0 ± 6.6% (mean ± SD) in freshly-isolated human NK-cells (FI-NK) and 83.4 ± 10.1% (mean ± SD) in NK-cells obtained from the NK-cell Activation and Expansion System (NKAES), with a sustained transgene expression for at least 21 days. BaEV-LVs outperformed Vesicular Stomatitis Virus type-G (VSV-G)-, RD114- and Measles Virus (MV)- pseudotyped LVs (p < 0.0001). mRNA expression of both BaEV receptors, ASCT1 and ASCT2, was detected in FI-NK and NKAES, with higher expression in NKAES. Transduction with BaEV-LVs encoding for CAR-CD22 resulted in robust CAR-expression on 38.3 ± 23.8% (mean ± SD) of NKAES cells, leading to specific killing of NK-resistant pre-B-ALL-RS4;11 cell line. Using a larger vector encoding a dual CD19/CD22-CAR, we were able to transduce and re-expand dual-CAR-expressing NKAES, even with lower viral titer. These dual-CAR-NK efficiently killed both CD19KO- and CD22KO-RS4;11 cells. Our results suggest that BaEV-LVs may efficiently enable NK-cell biological studies and translation of NK-cell-based immunotherapy to the clinic.

Social status alters chromatin accessibility and the gene regulatory response to glucocorticoid stimulation in rhesus macaques.

Low social status is an important predictor of disease susceptibility and mortality risk in humans and other social mammals. These effects are thought to stem in part from dysregulation of the glucocorticoid (GC)-mediated stress response. However, the molecular mechanisms that connect low social status and GC dysregulation to downstream health outcomes remain elusive. Here, we used an in vitro GC challenge to investigate the consequences of experimentally manipulated social status (i.e., dominance rank) for immune cell gene regulation in female rhesus macaques, using paired control and GC-treated peripheral blood mononuclear cell samples. We show that social status not only influences immune cell gene expression but also chromatin accessibility at hundreds of regions in the genome. Social status effects on gene expression were less pronounced following GC treatment than under control conditions. In contrast, social status effects on chromatin accessibility were stable across conditions, resulting in an attenuated relationship between social status, chromatin accessibility, and gene expression after GC exposure. Regions that were more accessible in high-status animals and regions that become more accessible following GC treatment were enriched for a highly concordant set of transcription factor binding motifs, including motifs for the GC receptor cofactor AP-1. Together, our findings support the hypothesis that social status alters the dynamics of GC-mediated gene regulation and identify chromatin accessibility as a mechanism involved in social stress-driven GC resistance. More broadly, they emphasize the context-dependent nature of social status effects on gene regulation and implicate epigenetic remodeling of chromatin accessibility as a contributing factor.

Projecting social contact matrices to different demographic structures.

The modeling of large-scale communicable epidemics has greatly benefited in the last years from the increasing availability of highly detailed data. Particullarly, in order to achieve quantitative descriptions of the evolution of epidemics, contact networks and mixing patterns are key. These heterogeneous patterns depend on several factors such as location, socioeconomic conditions, time, and age. This last factor has been shown to encapsulate a large fraction of the observed inter-individual variation in contact patterns, an observation validated by different measurements of age-dependent contact matrices. Recently, several works have studied how to project those matrices to areas where empirical data are not available. However, the dependence of contact matrices on demographic structures and their time evolution has been largely neglected. In this work, we tackle the problem of how to transform an empirical contact matrix that has been obtained for a given demographic structure into a different contact matrix that is compatible with a different demography. The methodology discussed here allows to extrapolate a contact structure measured in a particular area to any other whose demographic structure is known, as well as to obtain the time evolution of contact matrices as a function of the demographic dynamics of the populations they refer to. To quantify the effect of considering time-dynamics of contact patterns on disease modeling, we implemented a Susceptible-Exposed-Infected-Recovered (SEIR) model on 16 different countries and regions and evaluated the impact of neglecting the temporal evolution of mixing patterns. Our results show that simulated disease incidence rates, both at the aggregated and age-specific levels, are significantly dependent on contact structures variation driven by demographic evolution. The present work opens the path to eliminate technical biases from model-based impact evaluations of future epidemic threats and warns against the use of contact matrices to model diseases without correcting for demographic evolution or geographic variations.

Genetic and evolutionary determinants of human population variation in immune responses.

Humans display remarkable immune response variation when exposed to identical immune challenges. However, our understanding of the genetic, evolutionary, and environmental factors that impact this inter-individual and inter-population immune response heterogeneity is still in its early days. In this review, we discuss three fundamental questions concerning the recent evolution of the human immune system: the degree to which individuals from different populations vary in their innate immune responses, the genetic variants accounting for such differences, and the evolutionary mechanisms that led to the establishment of these variants in modern human populations. We also discuss how past selective events might have contributed to the uneven distribution of immune-related disorders across populations.

Data-driven model for the assessment of Mycobacterium tuberculosis transmission in evolving demographic structures.

In the case of tuberculosis (TB), the capabilities of epidemic models to produce quantitatively robust forecasts are limited by multiple hindrances. Among these, understanding the complex relationship between disease epidemiology and populations' age structure has been highlighted as one of the most relevant. TB dynamics depends on age in multiple ways, some of which are traditionally simplified in the literature. That is the case of the heterogeneities in contact intensity among different age strata that are common to all airborne diseases, but still typically neglected in the TB case. Furthermore, while demographic structures of many countries are rapidly aging, demographic dynamics are pervasively ignored when modeling TB spreading. In this work, we present a TB transmission model that incorporates country-specific demographic prospects and empirical contact data around a data-driven description of TB dynamics. Using our model, we find that the inclusion of demographic dynamics is followed by an increase in the burden levels predicted for the next decades in the areas of the world that are most hit by the disease today. Similarly, we show that considering realistic patterns of contacts among individuals in different age strata reshapes the transmission patterns reproduced by the models, a result with potential implications for the design of age-focused epidemiological interventions.