Comparative 'omics analyses differentiate Mycobacterium tuberculosis and Mycobacterium bovis and reveal distinct macrophage responses to infection with the human and bovine tubercle bacilli

Multilevel Proteomics Reveals Epigenetic Signatures in BCG-Mediated Macrophage Activation

The bacillus Calmette-Guérin BCG vaccine (Mycobacterium bovis) is primarily used to prevent tuberculosis (TB) infections but has wide-ranging immunogenic effects. One of its most notable properties is its ability to induce trained immunity, a memory-like response in innate immune cells such as macrophages. Through targeted analyses of well-established histone marks, prior research has shown that these changes are generated through epigenetic modification. Mass spectrometry-based proteomic approaches provide a way to globally profile various aspects of the proteome, providing data to further identify unexplored mechanisms of BCG-mediated immunomodulation. Here we use multi-level proteomics (total, histone, and phospho to identify networks and potential mechanisms that mediate BCG-induced immunomodulation in macrophages. Histone-focused proteomics and total proteomics were performed at the University of Cape Town (data available via ProteomeXchange with identifier PXD051187), while phosphoproteomics data was retrieved from the ProteomeXchange Repository (identifier PXD013171). We identify several epigenetic mechanisms that may drive BCG-induced training phenotypes. Evidence across the proteomics and histone-focused proteomics data set pair 6 epigenetic effectors (NuA4, NuRD, NSL, Sin3A, SIRT2, SIRT6) and their substrates.

Preferential differential gene expression within the WC1.1+ γδ T cell compartment in cattle naturally infected with Mycobacterium bovis

Bovine tuberculosis (bTB), caused by infection with Mycobacterium bovis, continues to cause significant issues for the global agriculture industry as well as for human health. An incomplete understanding of the host immune response contributes to the challenges of control and eradication of this zoonotic disease. In this study, high-throughput bulk RNA sequencing (RNA-seq) was used to characterise differential gene expression in γδ T cells - a subgroup of T cells that bridge innate and adaptive immunity and have known anti-mycobacterial response mechanisms. γδ T cell subsets are classified based on expression of a pathogen-recognition receptor known as Workshop Cluster 1 (WC1) and we hypothesised that bTB disease may alter the phenotype and function of specific γδ T cell subsets. Peripheral blood was collected from naturally M. bovis-infected (positive for single intradermal comparative tuberculin test (SICTT) and IFN-γ ELISA) and age- and sex-matched, non-infected control Holstein-Friesian cattle. γδ T subsets were isolated using fluorescence activated cell sorting (n = 10-12 per group) and high-quality RNA extracted from each purified lymphocyte subset (WC1.1+, WC1.2+, WC1- and γδ-) was used to generate transcriptomes using bulk RNA-seq (n = 6 per group, representing a total of 48 RNA-seq libraries). Relatively low numbers of differentially expressed genes (DEGs) were observed between most cell subsets; however, 189 genes were significantly differentially expressed in the M. bovis-infected compared to the control groups for the WC1.1+ γδ T cell compartment (absolute log2 FC ≥ 1.5 and FDR P adj. ≤ 0.1). The majority of these DEGs (168) were significantly increased in expression in cells from the bTB+ cattle and included genes encoding transcription factors (TBX21 and EOMES), chemokine receptors (CCR5 and CCR7), granzymes (GZMA, GZMM, and GZMH) and multiple killer cell immunoglobulin-like receptor (KIR) proteins indicating cytotoxic functions. Biological pathway overrepresentation analysis revealed enrichment of genes with multiple immune functions including cell activation, proliferation, chemotaxis, and cytotoxicity of lymphocytes. In conclusion, γδ T cells have important inflammatory and regulatory functions in cattle, and we provide evidence for preferential differential activation of the WC1.1+ specific subset in cattle naturally infected with M. bovis.

Global-scale GWAS associates a subset of SNPs with animal-adapted variants in M. tuberculosis complex

BACKGROUND

While Mycobacterium tuberculosis complex (MTBC) variants are clonal, variant tuberculosis is a human-adapted pathogen, and variant bovis infects many hosts. Despite nucleotide identity between MTBC variants exceeding 99.95%, it remains unclear what drives these differences. Markers of adaptation into variants were sought by bacterial genome-wide association study of single nucleotide polymorphisms extracted from 6,362 MTBC members from varied hosts and countries.

RESULTS

The search identified 120 genetic loci associated with MTBC variant classification and certain hosts. In many cases, these changes are uniformly fixed in certain variants while absent in others in this dataset, providing good discriminatory power in distinguishing variants by polymorphisms. Multiple changes were seen in genes for cholesterol and fatty acid metabolism, pathways previously proposed to be important for host adaptation, including Mce4F (part of the fundamental cholesterol intake Mce4 pathway), 4 FadD and FadE genes (playing roles in cholesterol and fatty acid utilization), and other targets like Rv3548c and PTPB, genes shown essential for growth on cholesterol by transposon studies.

CONCLUSIONS

These findings provide a robust set of genetic loci associated with the split of variant bovis and variant tuberculosis, and suggest that adaptation to new hosts could involve adjustments in uptake and catabolism of cholesterol and fatty acids, like the proposed specialization to different populations in MTB lineages by alterations to host lipid composition. Future studies are required to elucidate how the associations between cholesterol profiles and pathogen utilization differences between hosts and MTBC variants, as well as the investigation of uncharacterized genes discovered in this study. This information will likely provide an understanding on the diversification of MBO away from humans and specialization towards a broad host range.

Exploring virulence in Mycobacterium bovis: clues from comparative genomics and perspectives for the future

Here we provide a summary of a plenary lecture delivered on Mycobacterium bovis, the bovine TB bacillus, at the M. bovis 2022 meeting held in Galway, Ireland, in June 2022. We focus on the analysis of genetic differences between M. bovis and the human pathogen Mycobacterium tuberculosis as a route to gain knowledge on what makes M. bovis function as an animal pathogen. We provide a brief historical background around M. bovis and comparative virulence experiments with M. tuberculosis, before moving to what we have learned from the studies of the M. bovis genome sequence. We discuss the need to translate knowledge on the molecular basis of virulence in M. bovis into improved control of bovine tuberculosis.

Animal models for COVID-19 and tuberculosis

Respiratory infections cause tremendous morbidity and mortality worldwide. Amongst these diseases, tuberculosis (TB), a bacterial illness caused by Mycobacterium tuberculosis which often affects the lung, and coronavirus disease 2019 (COVID-19) caused by the Severe Acute Respiratory Syndrome Coronavirus type 2 (SARS-CoV-2), stand out as major drivers of epidemics of global concern. Despite their unrelated etiology and distinct pathology, these infections affect the same vital organ and share immunopathogenesis traits and an imperative demand to model the diseases at their various progression stages and localizations. Due to the clinical spectrum and heterogeneity of both diseases experimental infections were pursued in a variety of animal models. We summarize mammalian models employed in TB and COVID-19 experimental investigations, highlighting the diversity of rodent models and species peculiarities for each infection. We discuss the utility of non-human primates for translational research and emphasize on the benefits of non-conventional experimental models such as livestock. We epitomize advances facilitated by animal models with regard to understanding disease pathophysiology and immune responses. Finally, we highlight research areas necessitating optimized models and advocate that research of pulmonary infectious diseases could benefit from cross-fertilization between studies of apparently unrelated diseases, such as TB and COVID-19.

In-depth systems biological evaluation of bovine alveolar macrophages suggests novel insights into molecular mechanisms underlying Mycobacterium bovis infection

Objective

Bovine tuberculosis (bTB) is a chronic respiratory infectious disease of domestic livestock caused by intracellular Mycobacterium bovis infection, which causes ~$3 billion in annual losses to global agriculture. Providing novel tools for bTB managements requires a comprehensive understanding of the molecular regulatory mechanisms underlying the M. bovis infection. Nevertheless, a combination of different bioinformatics and systems biology methods was used in this study in order to clearly understand the molecular regulatory mechanisms of bTB, especially the immunomodulatory mechanisms of M. bovis infection.

Methods

RNA-seq data were retrieved and processed from 78 (39 non-infected control vs. 39 M. bovis-infected samples) bovine alveolar macrophages (bAMs). Next, weighted gene co-expression network analysis (WGCNA) was performed to identify the co-expression modules in non-infected control bAMs as reference set. The WGCNA module preservation approach was then used to identify non-preserved modules between non-infected controls and M. bovis-infected samples (test set). Additionally, functional enrichment analysis was used to investigate the biological behavior of the non-preserved modules and to identify bTB-specific non-preserved modules. Co-expressed hub genes were identified based on module membership (MM) criteria of WGCNA in the non-preserved modules and then integrated with protein-protein interaction (PPI) networks to identify co-expressed hub genes/transcription factors (TFs) with the highest maximal clique centrality (MCC) score (hub-central genes).

Results

As result, WGCNA analysis led to the identification of 21 modules in the non-infected control bAMs (reference set), among which the topological properties of 14 modules were altered in the M. bovis-infected bAMs (test set). Interestingly, 7 of the 14 non-preserved modules were directly related to the molecular mechanisms underlying the host immune response, immunosuppressive mechanisms of M. bovis, and bTB development. Moreover, among the co-expressed hub genes and TFs of the bTB-specific non-preserved modules, 260 genes/TFs had double centrality in both co-expression and PPI networks and played a crucial role in bAMs-M. bovis interactions. Some of these hub-central genes/TFs, including PSMC4, SRC, BCL2L1, VPS11, MDM2, IRF1, CDKN1A, NLRP3, TLR2, MMP9, ZAP70, LCK, TNF, CCL4, MMP1, CTLA4, ITK, IL6, IL1A, IL1B, CCL20, CD3E, NFKB1, EDN1, STAT1, TIMP1, PTGS2, TNFAIP3, BIRC3, MAPK8, VEGFA, VPS18, ICAM1, TBK1, CTSS, IL10, ACAA1, VPS33B, and HIF1A, had potential targets for inducing immunomodulatory mechanisms by M. bovis to evade the host defense response.

Conclusion

The present study provides an in-depth insight into the molecular regulatory mechanisms behind M. bovis infection through biological investigation of the candidate non-preserved modules directly related to bTB development. Furthermore, several hub-central genes/TFs were identified that were significant in determining the fate of M. bovis infection and could be promising targets for developing novel anti-bTB therapies and diagnosis strategies.

The rate and role of pseudogenes of the Mycobacterium tuberculosis complex

Whole-genome sequence analyses have significantly contributed to the understanding of virulence and evolution of the Mycobacterium tuberculosis complex (MTBC), the causative pathogens of tuberculosis. Most MTBC evolutionary studies are focused on single nucleotide polymorphisms and deletions, but rare studies have evaluated gene content, whereas none has comprehensively evaluated pseudogenes. Accordingly, we describe an extensive study focused on quantifying and predicting possible functions of MTBC and Mycobacterium canettii pseudogenes. Using NCBI's PGAP-detected pseudogenes, we analysed 25 837 pseudogenes from 158 MTBC and M. canetii strains and combined transcriptomics and proteomics of M. tuberculosis H37Rv to gain insights about pseudogenes' expression. Our results indicate significant variability concerning rate and conservancy of in silico predicted pseudogenes among different ecotypes and lineages of tuberculous mycobacteria and pseudogenization of important virulence factors and genes of the metabolism and antimicrobial resistance/tolerance. We show that in silico predicted pseudogenes contribute considerably to MTBC genetic diversity at the population level. Moreover, the transcription machinery of M. tuberculosis can fully transcribe most pseudogenes, indicating intact promoters and recent pseudogene evolutionary emergence. Proteomics of M. tuberculosis and close evaluation of mutational lesions driving pseudogenization suggest that few in silico predicted pseudogenes are likely capable of neofunctionalization, nonsense mutation reversal, or phase variation, contradicting the classical definition of pseudogenes. Such findings indicate that genome annotation should be accompanied by proteomics and protein function assays to improve its accuracy. While indels and insertion sequences are the main drivers of the observed mutational lesions in these species, population bottlenecks and genetic drift are likely the evolutionary processes acting on pseudogenes' emergence over time. Our findings unveil a new perspective on MTBC's evolution and genetic diversity.

Evidence, Challenges, and Knowledge Gaps Regarding Latent Tuberculosis in Animals

Mycobacterium bovis and other Mycobacterium tuberculosis complex (MTBC) pathogens that cause domestic animal and wildlife tuberculosis have received considerably less attention than M. tuberculosis, the primary cause of human tuberculosis (TB). Human TB studies have shown that different stages of infection can exist, driven by host-pathogen interactions. This results in the emergence of heterogeneous subpopulations of mycobacteria in different phenotypic states, which range from actively replicating (AR) cells to viable but slowly or non-replicating (VBNR), viable but non-culturable (VBNC), and dormant mycobacteria. The VBNR, VBNC, and dormant subpopulations are believed to underlie latent tuberculosis (LTB) in humans; however, it is unclear if a similar phenomenon could be happening in animals. This review discusses the evidence, challenges, and knowledge gaps regarding LTB in animals, and possible host-pathogen differences in the MTBC strains M. tuberculosis and M. bovis during infection. We further consider models that might be adapted from human TB research to investigate how the different phenotypic states of bacteria could influence TB stages in animals. In addition, we explore potential host biomarkers and mycobacterial changes in the DosR regulon, transcriptional sigma factors, and resuscitation-promoting factors that may influence the development of LTB.

Defining the Genes Required for Survival of Mycobacterium bovis in the Bovine Host Offers Novel Insights into the Genetic Basis of Survival of Pathogenic Mycobacteria

Tuberculosis has severe impacts on both humans and animals. Understanding the genetic basis of survival of both Mycobacterium tuberculosis, the human-adapted species, and Mycobacterium bovis, the animal-adapted species, is crucial to deciphering the biology of both pathogens. There are several studies that identify the genes required for survival of M. tuberculosis in vivo using mouse models; however, there are currently no studies probing the genetic basis of survival of M. bovis in vivo. In this study, we utilize transposon insertion sequencing in M. bovis AF2122/97 to determine the genes required for survival in cattle. We identify genes encoding established mycobacterial virulence functions such as the ESX-1 secretion system, phthiocerol dimycocerosate (PDIM) synthesis, mycobactin synthesis, and cholesterol catabolism that are required in vivo. We show that, as in M. tuberculosis H37Rv, phoPR is required by M. bovis AF2122/97 in vivo despite the known defect in signaling through this system. Comparison to studies performed in species that are able to use carbohydrates as an energy source, such as M. bovis BCG and M. tuberculosis, suggests that there are differences in the requirement for genes involved in cholesterol import (mce4 operon) and oxidation (hsd). We report a good correlation with existing mycobacterial virulence functions but also find several novel virulence factors, including genes involved in protein mannosylation, aspartate metabolism, and glycerol-phosphate metabolism. These findings further extend our knowledge of the genetic basis of survival in vivo in bacteria that cause tuberculosis and provide insight for the development of novel diagnostics and therapeutics. IMPORTANCE This is the first report of the genetic requirements of an animal-adapted member of the Mycobacterium tuberculosis complex (MTBC) in a natural host. M. bovis has devastating impacts on cattle, and bovine tuberculosis is a considerable economic, animal welfare, and public health concern. The data highlight the importance of mycobacterial cholesterol catabolism and identify several new virulence factors. Additionally, the work informs the development of novel differential diagnostics and therapeutics for TB in both human and animal populations.

Rational development of mycobacteria cell factory for advancing the steroid biomanufacturing

Steroidal resource occupies a vital proportion in the pharmaceutical industry attributing to their important therapeutic effects on fertility, anti-inflammatory and antiviral activities. Currently, microbial transformation from phytosterol has become the dominant strategy of steroidal drug intermediate synthesis that bypasses the traditional chemical route. Mycobacterium sp. serve as the main industrial microbial strains that are capable of introducing selective functional modifications of steroidal intermediate, which has become an indispensable platform for steroid biomanufacturing. By reviewing the progress in past two decades, the present paper concentrates mainly on the microbial rational modification aspects that include metabolic pathway editing, key enzymes engineering, material transport pathway reinforcement, toxic metabolic intermediates removal and byproduct reconciliation. In addition, progress on omics analysis and direct genetic manipulation are summarized and classified that may help reform the industrial hosts with more efficiency. The paper provides an insightful present for steroid biomanufacturing especially on the current trends and prospects of mycobacteria.

The TbD1 Locus Mediates a Hypoxia-Induced Copper Response in Mycobacterium bovis

The Mycobacterium tuberculosis complex (MTBC) contains the causative agents of tuberculosis (TB) in mammals. The archetypal members of the MTBC, Mycobacterium tuberculosis and Mycobacterium bovis, cause human tuberculosis and bovine tuberculosis, respectively. Although M. tuberculosis and M. bovis share over 99.9% genome identity, they show distinct host adaptation for humans and animals; hence, while the molecular basis of host adaptation is encoded in their genomes, the mechanistic basis of host tropism is still unclear. Exploration of the in vitro phenotypic consequences of known genetic difference between M. bovis and M. tuberculosis offers one route to explore genotype–phenotype links that may play a role in host adaptation. The TbD1 (“Mycobacterium tuberculosis deletion 1 region”) locus encompasses the mmpS6 and mmpL6 genes. TbD1 is absent in M. tuberculosis “modern” lineages (Lineages 2, 3, and 4) but present in “ancestral” M. tuberculosis (Lineages 1 and 7), Mycobacterium africanum lineages (Lineages 5 and 6), newly identified M. tuberculosis lineages (Lineages 8 and 9), and animal adapted strains, such as M. bovis. The function of TbD1 has previously been investigated in M. tuberculosis, where conflicting data has emerged on the role of TbD1 in sensitivity to oxidative stress, while the underlying mechanistic basis of such a phenotype is unclear. In this study, we aimed to shed further light on the role of the TbD1 locus by exploring its function in M. bovis. Toward this, we constructed an M. bovis TbD1 knockout (ΔTbD1) strain and conducted comparative transcriptomics to define global gene expression profiles of M. bovis wild-type (WT) and the ΔTbD1 strains under in vitro culture conditions (rolling and standing cultures). This analysis revealed differential induction of a hypoxia-driven copper response in WT and ΔTbD1 strains. In vitro phenotypic assays demonstrated that the deletion of TbD1 sensitized M. bovis to H2O2 and hypoxia-specific copper toxicity. Our study provides new information on the function of the TbD1 locus in M. bovis and its role in stress responses in the MTBC.

Probing Differences in Gene Essentiality Between the Human and Animal Adapted Lineages of the Mycobacterium tuberculosis Complex Using TnSeq

Members of the Mycobacterium tuberculosis complex (MTBC) show distinct host adaptations, preferences and phenotypes despite being >99% identical at the nucleic acid level. Previous studies have explored gene expression changes between the members, however few studies have probed differences in gene essentiality. To better understand the functional impacts of the nucleic acid differences between Mycobacterium bovis and Mycobacterium tuberculosis, we used the Mycomar T7 phagemid delivery system to generate whole genome transposon libraries in laboratory strains of both species and compared the essentiality status of genes during growth under identical in vitro conditions. Libraries contained insertions in 54% of possible TA sites in M. bovis and 40% of those present in M. tuberculosis, achieving similar saturation levels to those previously reported for the MTBC. The distributions of essentiality across the functional categories were similar in both species. 527 genes were found to be essential in M. bovis whereas 477 genes were essential in M. tuberculosis and 370 essential genes were common in both species. CRISPRi was successfully utilised in both species to determine the impacts of silencing genes including wag31, a gene involved in peptidoglycan synthesis and Rv2182c/Mb2204c, a gene involved in glycerophospholipid metabolism. We observed species specific differences in the response to gene silencing, with the inhibition of expression of Mb2204c in M. bovis showing significantly less growth impact than silencing its orthologue (Rv2182c) in M. tuberculosis. Given that glycerophospholipid metabolism is a validated pathway for antimicrobials, our observations suggest that target vulnerability in the animal adapted lineages cannot be assumed to be the same as the human counterpart. This is of relevance for zoonotic tuberculosis as it implies that the development of antimicrobials targeting the human adapted lineage might not necessarily be effective against the animal adapted lineage. The generation of a transposon library and the first reported utilisation of CRISPRi in M. bovis will enable the use of these tools to further probe the genetic basis of survival under disease relevant conditions.

Challenges in defining the functional, non-coding, expressed genome of members of the Mycobacterium tuberculosis complex

A definitive transcriptome atlas for the non-coding expressed elements of the members of the Mycobacterium tuberculosis complex (MTBC) does not exist. Incomplete lists of non-coding transcripts can be obtained for some of the reference genomes (e.g., M. tuberculosis H37Rv) but to what extent these transcripts have homologues in closely related species or even strains is not clear. This has implications for the analysis of transcriptomic data; non-coding parts of the transcriptome are often ignored in the absence of formal, reliable annotation. Here, we review the state of our knowledge of non-coding RNAs in pathogenic mycobacteria, emphasizing the disparities in the information included in commonly used databases. We then proceed to review ways of combining computational solutions for predicting the non-coding transcriptome with experiments that can help refine and confirm these predictions.

Comparative transcriptomic analysis of THP-1-derived macrophages infected with Mycobacterium tuberculosis H37Rv, H37Ra and BCG

Tuberculosis (TB) remains a worldwide healthcare concern, and the exploration of the host‐pathogen interaction is essential to develop therapeutic modalities and strategies to control Mycobacterium tuberculosis (M.tb). In this study, RNA sequencing (transcriptome sequencing) was employed to investigate the global transcriptome changes in the macrophages during the different strains of M.tb infection. THP‐1 cells derived from macrophages were exposed to the virulent M.tb strain H37Rv (Rv) or the avirulent M.tb strain H37Ra (Ra), and the M.tb BCG vaccine strain was used as a control. The cDNA libraries were prepared from M.tb‐infected macrophages and then sequenced. To assess the transcriptional differences between the expressed genes, the bioinformatics analysis was performed using a standard pipeline of quality control, reference mapping, differential expression analysis, protein‐protein interaction (PPI) networks, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Q‐PCR and Western blot assays were also performed to validate the data. Our findings indicated that, when compared to BCG or M.tb H37Ra infection, the transcriptome analysis identified 66 differentially expressed genes in the M.tb H37Rv‐infected macrophages, out of which 36 genes were up‐regulated, and 30 genes were down‐regulated. The up‐regulated genes were associated with immune response regulation, chemokine secretion, and leucocyte chemotaxis. In contrast, the down‐regulated genes were associated with amino acid biosynthetic and energy metabolism, connective tissue development and extracellular matrix organization. The Q‐PCR and Western blot assays confirmed increased expression of pro‐inflammatory factors, altered energy metabolic processes, enhanced activation of pro‐inflammatory signalling pathways and increased pyroptosis in H37Rv‐infected macrophage. Overall, our RNA sequencing‐based transcriptome study successfully identified a comprehensive, in‐depth gene expression/regulation profile in M.tb‐infected macrophages. The results demonstrated that virulent M.tb strain H37Rv infection triggers a more severe inflammatory immune response associated with increased tissue damage, which helps in understanding the host‐pathogen interaction dynamics and pathogenesis features in different strains of M.tb infection.

Effect of co-positivity for brucellosis and tuberculosis on milk yield and fertility of Holstein cows

This study aimed to determine whether cows detected as tuberculosis (bTB) reactors and seropositive to brucellosis (bBR), as well as co-positive to bBR and bTB (bBR-bTB) and with a complete lactation before slaughter, were associated with reduced milk production and fertility. A total of 8068 productive and reproductive records of high-yielding Holstein cows from a single large dairy herd with a high prevalence of bTB and bBR were collected from 2012 to 2015. Lactation derived either from calving (n = 6019) or hormonally induced lactation (n = 2049), and all cows received growth hormone throughout lactation. For cows not induced into lactation, pregnancy rate to first service for healthy cows (C; 26.6%) was higher (P < 0.01) than bBR (15.2%), bTB (15.8%), and bBR-bTB (1.3%) cows. For induced cows, pregnancy rate to first service did not differ significantly among C, bBR, and bTB (14.5-17.3%) cows, but the percentage success of first service was extremely low (1.3%; P < 0.01) in bBR-bTB cows. Services per pregnancy (only pregnant cows) were lowest for C (3.3 ± 2.9; P < 0.01) and highest (6.4 ± 3.4) for bBR-bTB non-induced cows. This variable was lowest for C (2.9 ± 2.5; P < 0.01) and highest for bBR-bTB non-induced cows (6.3 ± 3.1). Pregnancy rate to all services did not differed for C (79.5%), bBR (76.7%), and bTB (75.9%) but was lower (58.9%; P < 0.01) for bBR-bTB non-induced cows. For induced cows this variable was highest for bBR (53.3%) and lowest for bBR-bTB (34.1%; P < 0.01) non-induced cows. 305-d milk production was increased by 4%, and total milk yield by 7% in TB-positive cows compared to that of the negative cows non-induced hormonally into lactation. This study showed the negative impact of the co-positivity for bTB and bBR on the reproductive efficiency of Holstein cows, although positive bTB and bBR tests enhanced milk yield.

Mycobacterial Infection of Precision-Cut Lung Slices Reveals Type 1 Interferon Pathway Is Locally Induced by Mycobacterium bovis but Not M. tuberculosis in a Cattle Breed

Tuberculosis exacts a terrible toll on human and animal health. While Mycobacterium tuberculosis (Mtb) is restricted to humans, Mycobacterium bovis (Mb) is present in a large range of mammalian hosts. In cattle, bovine TB (bTB) is a noticeable disease responsible for important economic losses in developed countries and underestimated zoonosis in the developing world. Early interactions that take place between mycobacteria and the lung tissue early after aerosol infection govern the outcome of the disease. In cattle, these early steps remain poorly characterized. The precision-cut lung slice (PCLS) model preserves the structure and cell diversity of the lung. We developed this model in cattle in order to study the early lung response to mycobacterial infection. In situ imaging of PCLS infected with fluorescent Mb revealed bacilli in the alveolar compartment, in adjacent or inside alveolar macrophages, and in close contact with pneumocytes. We analyzed the global transcriptional lung inflammation signature following infection of PCLS with Mb and Mtb in two French beef breeds: Blonde d'Aquitaine and Charolaise. Whereas, lungs from the Blonde d'Aquitaine produced high levels of mediators of neutrophil and monocyte recruitment in response to infection, such signatures were not observed in the Charolaise in our study. In the Blonde d'Aquitaine lung, whereas the inflammatory response was highly induced by two Mb strains, AF2122 isolated from cattle in the UK and Mb3601 circulating in France, the response against two Mtb strains, H37Rv, the reference laboratory strain, and BTB1558, isolated from zebu in Ethiopia, was very low. Strikingly, the type I interferon pathway was only induced by Mb but not Mtb strains, indicating that this pathway may be involved in mycobacterial virulence and host tropism. Hence, the PCLS model in cattle is a valuable tool to deepen our understanding of early interactions between lung host cells and mycobacteria. It revealed striking differences between cattle breeds and mycobacterial strains. This model could help in deciphering biomarkers of resistance vs. susceptibility to bTB in cattle as such information is still critically needed for bovine genetic selection programs and would greatly help the global effort to eradicate bTB.

RNA-Seq Transcriptome Analysis of Peripheral Blood From Cattle Infected With Mycobacterium bovis Across an Experimental Time Course

Bovine tuberculosis, caused by infection with members of the Mycobacterium tuberculosis complex, particularly Mycobacterium bovis, is a major endemic disease affecting cattle populations worldwide, despite the implementation of stringent surveillance and control programs in many countries. The development of high-throughput functional genomics technologies, including RNA sequencing, has enabled detailed analysis of the host transcriptome to M. bovis infection, particularly at the macrophage and peripheral blood level. In the present study, we have analysed the transcriptome of bovine whole peripheral blood samples collected at −1 week pre-infection and +1, +2, +6, +10, and +12 weeks post-infection time points. Differentially expressed genes were catalogued and evaluated at each post-infection time point relative to the −1 week pre-infection time point and used for the identification of putative candidate host transcriptional biomarkers for M. bovis infection. Differentially expressed gene sets were also used for examination of cellular pathways associated with the host response to M. bovis infection, construction of de novo gene interaction networks enriched for host differentially expressed genes, and time-series analyses to identify functionally important groups of genes displaying similar patterns of expression across the infection time course. A notable outcome of these analyses was identification of a 19-gene transcriptional biosignature of infection consisting of genes increased in expression across the time course from +1 week to +12 weeks post-infection.

Integrative genomics of the mammalian alveolar macrophage response to intracellular mycobacteria

Background

Bovine TB (bTB), caused by infection with Mycobacterium bovis, is a major endemic disease affecting global cattle production. The key innate immune cell that first encounters the pathogen is the alveolar macrophage, previously shown to be substantially reprogrammed during intracellular infection by the pathogen. Here we use differential expression, and correlation- and interaction-based network approaches to analyse the host response to infection with M. bovis at the transcriptome level to identify core infection response pathways and gene modules. These outputs were then integrated with genome-wide association study (GWAS) data sets to enhance detection of genomic variants for susceptibility/resistance to M. bovis infection.

Results

The host gene expression data consisted of RNA-seq data from bovine alveolar macrophages (bAM) infected with M. bovis at 24 and 48 h post-infection (hpi) compared to non-infected control bAM. These RNA-seq data were analysed using three distinct computational pipelines to produce six separate gene sets: 1) DE genes filtered using stringent fold-change and P-value thresholds (DEG-24: 378 genes, DEG-48: 390 genes); 2) genes obtained from expression correlation networks (CON-24: 460 genes, CON-48: 416 genes); and 3) genes obtained from differential expression networks (DEN-24: 339 genes, DEN-48: 495 genes). These six gene sets were integrated with three bTB breed GWAS data sets by employing a new genomics data integration tool—gwinteR. Using GWAS summary statistics, this methodology enabled detection of 36, 102 and 921 prioritised SNPs for Charolais, Limousin and Holstein-Friesian, respectively.

Conclusions

The results from the three parallel analyses showed that the three computational approaches could identify genes significantly enriched for SNPs associated with susceptibility/resistance to M. bovis infection. Results indicate distinct and significant overlap in SNP discovery, demonstrating that network-based integration of biologically relevant transcriptomics data can leverage substantial additional information from GWAS data sets. These analyses also demonstrated significant differences among breeds, with the Holstein-Friesian breed GWAS proving most useful for prioritising SNPS through data integration. Because the functional genomics data were generated using bAM from this population, this suggests that the genomic architecture of bTB resilience traits may be more breed-specific than previously assumed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07643-w.

Macrophage-specific responses to human- and animal-adapted tubercle bacilli reveal pathogen and host factors driving multinucleated cell formation

The Mycobacterium tuberculosis complex (MTBC) is a group of related pathogens that cause tuberculosis (TB) in mammals. MTBC species are distinguished by their ability to sustain in distinct host populations. While Mycobacterium bovis (Mbv) sustains transmission cycles in cattle and wild animals and causes zoonotic TB, M. tuberculosis (Mtb) affects human populations and seldom causes disease in cattle. The host and pathogen determinants underlying host tropism between MTBC species are still unknown. Macrophages are the main host cell that encounters mycobacteria upon initial infection, and we hypothesised that early interactions between the macrophage and mycobacteria influence species-specific disease outcome. To identify factors that contribute to host tropism, we analysed blood-derived primary human and bovine macrophages (hMϕ or bMϕ, respectively) infected with Mbv and Mtb. We show that Mbv and Mtb reside in different cellular compartments and differentially replicate in hMϕ whereas both Mbv and Mtb efficiently replicate in bMϕ. Specifically, we show that out of the four infection combinations, only the infection of bMϕ with Mbv promoted the formation of multinucleated giant cells (MNGCs), a hallmark of tuberculous granulomas. Mechanistically, we demonstrate that both MPB70 from Mbv and extracellular vesicles released by Mbv-infected bMϕ promote macrophage multinucleation. Importantly, we extended our in vitro studies to show that granulomas from Mbv-infected but not Mtb-infected cattle contained higher numbers of MNGCs. Our findings implicate MNGC formation in the contrasting pathology between Mtb and Mbv for the bovine host and identify MPB70 from Mbv and extracellular vesicles from bMϕ as mediators of this process.

The identification of host and pathogen factors contributing to host-pathogen interaction is crucial to understand the pathogenesis and dissemination of tuberculosis. This is particularly the case in deciphering the mechanistic basis for host-tropism across the MTBC. Here, we show that in vitro, M. bovis but not M. tuberculosis induces multinucleated cell formation in bovine macrophages. We identified host and pathogen mechanistic drivers of multinucleated cell formation: MPB70 as the M. bovis factor and bovine macrophage extracellular vesicles. Using a cattle experimental infection model, we confirmed differential multinucleated cell formation in vivo. Thus, we have identified host and pathogen factors that contribute to host tropism in human/bovine tuberculosis. Additionally, this work provides an explanation for the long-standing association of multinucleated cells with tuberculosis pathogenesis.

Weighted Gene Co-Expression Network Analysis Identifies Key Modules and Hub Genes Associated with Mycobacterial Infection of Human Macrophages

Tuberculosis (TB) is still a leading cause of death worldwide. Treatments remain unsatisfactory due to an incomplete understanding of the underlying host–pathogen interactions during infection. In the present study, weighted gene co-expression network analysis (WGCNA) was conducted to identify key macrophage modules and hub genes associated with mycobacterial infection. WGCNA was performed combining our own transcriptomic results using Mycobacterium aurum-infected human monocytic macrophages (THP1) with publicly accessible datasets obtained from three types of macrophages infected with seven different mycobacterial strains in various one-to-one combinations. A hierarchical clustering tree of 11,533 genes was built from 198 samples, and 47 distinct modules were revealed. We identified a module, consisting of 226 genes, which represented the common response of host macrophages to different mycobacterial infections that showed significant enrichment in innate immune stimulation, bacterial pattern recognition, and leukocyte chemotaxis. Moreover, by network analysis applied to the 74 genes with the best correlation with mycobacteria infection, we identified the top 10 hub-connecting genes: NAMPT, IRAK2, SOCS3, PTGS2, CCL20, IL1B, ZC3H12A, ABTB2, GFPT2, and ELOVL7. Interestingly, apart from the well-known Toll-like receptor and inflammation-associated genes, other genes may serve as novel TB diagnosis markers and potential therapeutic targets.

The genetic proteome: Using genetics to inform the proteome of mycobacterial pathogens

Mycobacterial pathogens pose a sustained threat to human health. There is a critical need for new diagnostics, therapeutics, and vaccines targeting both tuberculous and nontuberculous mycobacterial species. Understanding the basic mechanisms used by diverse mycobacterial species to cause disease will facilitate efforts to design new approaches toward detection, treatment, and prevention of mycobacterial disease. Molecular, genetic, and biochemical approaches have been widely employed to define fundamental aspects of mycobacterial physiology and virulence. The recent expansion of genetic tools in mycobacteria has further increased the accessibility of forward genetic approaches. Proteomics has also emerged as a powerful approach to further our understanding of diverse mycobacterial species. Detection of large numbers of proteins and their modifications from complex mixtures of mycobacterial proteins is now routine, with efforts of quantification of these datasets becoming more robust. In this review, we discuss the “genetic proteome,” how the power of genetics, molecular biology, and biochemistry informs and amplifies the quality of subsequent analytical approaches and maximizes the potential of hypothesis-driven mycobacterial research. Published proteomics datasets can be used for hypothesis generation and effective post hoc supplementation to experimental data. Overall, we highlight how the integration of proteomics, genetic, molecular, and biochemical approaches can be employed successfully to define fundamental aspects of mycobacterial pathobiology.

Independent genomic polymorphisms in the PknH serine threonine kinase locus during evolution of the Mycobacterium tuberculosis Complex affect virulence and host preference

Species belonging to the Mycobacterium tuberculosis Complex (MTBC) show more than 99% genetic identity but exhibit distinct host preference and virulence. The molecular genetic changes that underly host specificity and infection phenotype within MTBC members have not been fully elucidated. Here, we analysed RD900 genomic region across MTBC members using whole genome sequences from 60 different MTBC strains so as to determine its role in the context of MTBC evolutionary history. The RD900 region comprises two homologous genes, pknH1 and pknH2, encoding a serine/threonine protein kinase PknH flanking the tbd2 gene. Our analysis revealed that RD900 has been independently lost in different MTBC lineages and different strains, resulting in the generation of a single pknH gene. Importantly, all the analysed M. bovis and M. caprae strains carry a conserved deletion within a proline rich-region of pknH, independent of the presence or absence of RD900. We hypothesized that deletion of pknH proline rich-region in M. bovis may affect PknH function, having a potential role in its virulence and evolutionary adaptation. To explore this hypothesis, we constructed two M. bovis ‘knock-in’ strains containing the M. tuberculosis pknH gene. Evaluation of their virulence phenotype in mice revealed a reduced virulence of both M. bovis knock-in strains compared to the wild type, suggesting that PknH plays an important role in the differential virulence phenotype of M. bovis vs M. tuberculosis.

Tuberculosis is caused in humans and animals by organisms from the Mycobacterium tuberculosis Complex (MTBC), that share more than 99% genetic identity but exhibit distinct host preference and virulence. While Mycobacterium tuberculosis is the main causative agent of human TB, Mycobacterium bovis is responsible for bovine TB disease, that exacts a tremendous economic burden worldwide, as well as being a zoonotic threat. Unlike the human restriction of M. tuberculosis, M. bovis has a broader host range and it has been found to be more virulent than M. tuberculosis in different animal models. However, the molecular basis for host preference and virulence divergence between M. tuberculosis and M. bovis is not fully elucidated. Here we study the genetic variations of the genomic region RD900 in the context of MTBC phylogeny. RD900 contains two genes encoding orthologues of the serine/threonine kinase PknH, which is linked to the regulation of several bacterial processes including virulence. We found that M. bovis pknH genes show a conserved deletion that is not present in M. tuberculosis strains, and we evaluated the potential impact of these variations in the regulation of M. bovis vs M. tuberculosis virulence through the construction and in vivo characterization of M. bovis pknH mutant strains.

Evidence for Anti-Viral Effects of Complete Freund's Adjuvant in the Mouse Model of Enterovirus Infection

Group B coxsackieviruses (CVBs) belonging to the genus, Enterovirus and contain six serotypes that induce various diseases, whose occurrence may involve the mediation of more than one serotype. We recently identified immunogenic epitopes within coxsackieviruses B3 (CVB3) viral protein 1 that induce anti-viral T cell responses in mouse models of CVB infections. In our investigations to determine the protective responses of the viral epitopes, we unexpectedly noted that animals immunized with complete Freund’s adjuvant (CFA) alone and later challenged with CVB3 were completely protected against myocarditis. Similarly, the pancreatitis-inducing ability of CVB3 was remarkably reduced to only 10% in the CFA group as opposed to 73.3% in the control group that received no CFA. Additionally, no mortalities were noted in the CFA group, whereas 40% of control animals died during the course of 21 days post-infection with CVB3. Taken together, our data suggest that the adjuvant effects of CFA may be sufficient for protection against CVB infections. These observations may provide new insights into our understanding of the occurrence of viral infections.

Complete Genome Sequences of Mycobacterium tuberculosis Isolates Subjected to 200 Continuous Passages

We report 6 draft genome sequences corresponding to Mycobacterium tuberculosis H37Rv and M. tuberculosis DR689, a Beijing isolate, plus their counterparts subjected to 200 continuous passages in Middlebrook 7H9 broth, either alone or with ox bile.

We report 6 draft genome sequences corresponding to Mycobacterium tuberculosis H37Rv and M. tuberculosis DR689, a Beijing isolate, plus their counterparts subjected to 200 continuous passages in Middlebrook 7H9 broth, either alone or with ox bile.

Does Mycobacterium bovis persist in cattle in a non-replicative latent state as Mycobacterium tuberculosis in human beings?

Members of the Mycobacterium tuberculosis complex (MTBC) are responsible for tuberculosis in several mammals. In this complex, Mycobacterium tuberculosis and Mycobacterium bovis, which are closely related, show host preference for humans and cattle, respectively. Although human and bovine tuberculosis are clinically similar, M. tuberculosis mostly causes latent infection in humans, whereas M. bovis frequently leads to an acute infection in cattle. This review attempts to connect the pathology in experimental animal models as well as the cellular responses to M. bovis and M. tuberculosis regarding the differences in protein expression and regulatory mechanisms of both pathogens that could explain their apparent divergent latency behaviour. The occurrence of latent bovine tuberculosis (bTB) would represent a serious complication for the eradication of the disease in cattle, with the risk of onward transmission to humans. Thus, understanding the physiological events that may lead to the state of latency in bTB could assist in the development of appropriate prevention and control tools.

Human Alveolar and Splenic Macrophage Populations Display a Distinct Transcriptomic Response to Infection With Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) infects alveolar macrophages (AMs), causing pulmonary tuberculosis (PTB), the most common form of the disease. Less frequently, Mtb is disseminated to many other organs and tissues, resulting in different extrapulmonary forms of TB. Nevertheless, very few studies have addressed the global mRNA response of human AMs, particularly from humans with the active form of the disease. Strikingly, almost no studies have addressed the response of human extrapulmonary macrophages to Mtb infection. In this pilot study, using microarray technology, we examined the transcriptomic ex vivo response of AMs from PTB patients (AMTBs) and AMs from control subjects (AMCTs) infected with two clinical isolates of Mtb. Furthermore, we also studied the infection response of human splenic macrophages (SMs) to Mtb isolates, as a model for extrapulmonary infection, and compared the transcriptomic response between AMs and SMs. Our results showed a striking difference in global mRNA profiles in response to infection between AMs and SMs, implicating a tissue-specific macrophage response to Mtb.

Alveolar Macrophage Chromatin Is Modified to Orchestrate Host Response to Mycobacterium bovis Infection

Bovine tuberculosis is caused by infection with Mycobacterium bovis, which can also cause disease in a range of other mammals, including humans. Alveolar macrophages are the key immune effector cells that first encounter M. bovis and how the macrophage epigenome responds to mycobacterial pathogens is currently not well understood. Here, we have used chromatin immunoprecipitation sequencing (ChIP-seq), RNA-seq and miRNA-seq to examine the effect of M. bovis infection on the bovine alveolar macrophage (bAM) epigenome. We show that H3K4me3 is more prevalent, at a genome-wide level, in chromatin from M. bovis-infected bAM compared to control non-infected bAM; this was particularly evident at the transcriptional start sites of genes that determine programmed macrophage responses to mycobacterial infection (e.g. M1/M2 macrophage polarisation). This pattern was also supported by the distribution of RNA Polymerase II (Pol II) ChIP-seq results, which highlighted significantly increased transcriptional activity at genes demarcated by permissive chromatin. Identification of these genes enabled integration of high-density genome-wide association study (GWAS) data, which revealed genomic regions associated with resilience to infection with M. bovis in cattle. Through integration of these data, we show that bAM transcriptional reprogramming occurs through differential distribution of H3K4me3 and Pol II at key immune genes. Furthermore, this subset of genes can be used to prioritise genomic variants from a relevant GWAS data set.

Transcriptomic Profiles of Monocyte-Derived Macrophages in Response to Escherichia coli is Associated with the Host Genetics

Reactive Nitrogen Species (RNS) are a group of bactericidal molecules produced by macrophages in response to pathogens in a process called oxidative burst. Nitric oxide (NO^-) is a member of RNS produced from arginine by inducible Nitric Oxide Synthase (iNOS) enzyme. The activity of iNOS and production of NO^- by macrophages following stimulation is one of the indicators of macrophage polarization towards M1/proinflammatory. Production of NO^- by bovine monocyte-derived macrophage (MDM) and mouse peritoneal macrophages has been shown to be strongly associated with host genetic with the heritability of 0.776 in bovine MDM and 0.8 in mouse peritoneal macrophages. However, the mechanism of genetic regulation of macrophage response has remained less explored. In the current study, the transcriptome of bovine MDMs was compared between two extreme phenotypes that had been classified as high and low responder based on NO^- production. The results showed that 179 and 392 genes were differentially expressed (DE) between high and low responder groups at 3 and 18 hours after exposure to Escherichia coli, respectively. A set of 11 Transcription Factors (TFs) (STAT1, IRF7, SPI1, STAT4, IRF1, HIF1A, FOXO3, REL, NFAT5, HIC1, and IRF4) at 3 hours and a set of 13 TFs (STAT1, IRF1, HIF1A, STAT4, ATF4, TP63, EGR1, CDKN2A, RBL1, E2F1, PRDM1, GATA3, and IRF4) at 18 hours after exposure to E. coli were identified to be differentially regulated between the high and low responder phenotypes. These TFs were found to be divided into two clusters of inflammatory- and hypoxia-related TFs. Functional analysis revealed that some key canonical pathways such as phagocytosis, chemotaxis, antigen presentation, and cell-to-cell signalling are enriched among the over-expressed genes by high responder phenotype. Based on the results of this study, it was inferred that the functional characteristics of bovine MDMs are associated with NO-based classification. Since NO^- production is strongly associated with host genetics, this study for the first time shows the distinct proinflammatory profiles of macrophages are controlled by the natural genetic polymorphism in an outbred population. In addition, the results suggest that genetics can be considered as a new dimension in the current model of macrophage polarization which is currently described by the combination of stimulants, only.

The Core Proteome of Biofilm-Grown Clinical Pseudomonas aeruginosa Isolates

Comparative genomics has greatly facilitated the identification of shared as well as unique features among individual cells or tissues, and thus offers the potential to find disease markers. While proteomics is recognized for its potential to generate quantitative maps of protein expression, comparative proteomics in bacteria has been largely restricted to the comparison of single cell lines or mutant strains. In this study, we used a data independent acquisition (DIA) technique, which enables global protein quantification of large sample cohorts, to record the proteome profiles of overall 27 whole genome sequenced and transcriptionally profiled clinical isolates of the opportunistic pathogen Pseudomonas aeruginosa. Analysis of the proteome profiles across the 27 clinical isolates grown under planktonic and biofilm growth conditions led to the identification of a core biofilm-associated protein profile. Furthermore, we found that protein-to-mRNA ratios between different P. aeruginosa strains are well correlated, indicating conserved patterns of post-transcriptional regulation. Uncovering core regulatory pathways, which drive biofilm formation and associated antibiotic tolerance in bacterial pathogens, promise to give clues to interactions between bacterial species and their environment and could provide useful targets for new clinical interventions to combat biofilm-associated infections.

Genome-wide mutational biases fuel transcriptional diversity in the Mycobacterium tuberculosis complex

The Mycobacterium tuberculosis complex (MTBC) members display different host-specificities and virulence phenotypes. Here, we have performed a comprehensive RNAseq and methylome analysis of the main clades of the MTBC and discovered unique transcriptional profiles. The majority of genes differentially expressed between the clades encode proteins involved in host interaction and metabolic functions. A significant fraction of changes in gene expression can be explained by positive selection on single mutations that either create or disrupt transcriptional start sites (TSS). Furthermore, we show that clinical strains have different methyltransferases inactivated and thus different methylation patterns. Under the tested conditions, differential methylation has a minor direct role on transcriptomic differences between strains. However, disruption of a methyltransferase in one clinical strain revealed important expression differences suggesting indirect mechanisms of expression regulation. Our study demonstrates that variation in transcriptional profiles are mainly due to TSS mutations and have likely evolved due to differences in host characteristics.

Genomic determinants of speciation and spread of the Mycobacterium tuberculosis complex

Models on how bacterial lineages differentiate increase our understanding of early bacterial speciation events and the genetic loci involved. Here, we analyze the population genomics events leading to the emergence of the tuberculosis pathogen. The emergence is characterized by a combination of recombination events involving core pathogenesis functions and purifying selection on early diverging loci. We identify the phoR gene, the sensor kinase of a two-component system involved in virulence, as a key functional player subject to pervasive positive selection after the divergence of the Mycobacterium tuberculosis complex from its ancestor. Previous evidence showed that phoR mutations played a central role in the adaptation of the pathogen to different host species. Now, we show that phoR mutations have been under selection during the early spread of human tuberculosis, during later expansions, and in ongoing transmission events. Our results show that linking pathogen evolution across evolutionary and epidemiological time scales points to past and present virulence determinants.

The bovine alveolar macrophage DNA methylome is resilient to infection with Mycobacterium bovis

DNA methylation is pivotal in orchestrating gene expression patterns in various mammalian biological processes. Perturbation of the bovine alveolar macrophage (bAM) transcriptome, due to Mycobacterium bovis (M. bovis) infection, has been well documented; however, the impact of this intracellular pathogen on the bAM epigenome has not been determined. Here, whole genome bisulfite sequencing (WGBS) was used to assess the effect of M. bovis infection on the bAM DNA methylome. The methylomes of bAM infected with M. bovis were compared to those of non-infected bAM 24 hours post-infection (hpi). No differences in DNA methylation (CpG or non-CpG) were observed. Analysis of DNA methylation at proximal promoter regions uncovered >250 genes harbouring intermediately methylated (IM) promoters (average methylation of 33–66%). Gene ontology analysis, focusing on genes with low, intermediate or highly methylated promoters, revealed that genes with IM promoters were enriched for immune-related GO categories; this enrichment was not observed for genes in the high or low methylation groups. Targeted analysis of genes in the IM category confirmed the WGBS observation. This study is the first in cattle examining genome-wide DNA methylation at single nucleotide resolution in an important bovine cellular host-pathogen interaction model, providing evidence for IM promoter methylation in bAM.

PPE37 Is Essential for Mycobacterium tuberculosis Heme-Iron Acquisition (HIA), and a Defective PPE37 in Mycobacterium bovis BCG Prevents HIA

Mycobacterium tuberculosis, one of the world's leading causes of death, must acquire nutrients, such as iron, from the host to multiply and cause disease. Iron is an essential metal and M. tuberculosis possesses two different systems to acquire iron from its environment: siderophore-mediated iron acquisition (SMIA) and heme-iron acquisition (HIA), involving uptake and degradation of heme to release ferrous iron. We have discovered that Mycobacterium bovis BCG, the tuberculosis vaccine strain, is severely deficient in HIA, and we exploited this phenotypic difference between BCG and M. tuberculosis to identify genes involved in HIA by complementing BCG's defect with a fosmid library. We identified ppe37, an iron-regulated PPE family gene, as being essential for HIA. BCG complemented with M. tuberculosisppe37 exhibits HIA as efficient as that of M. tuberculosis, achieving robust growth with <0.2 µM hemin. Conversely, deletion of ppe37 from M. tuberculosis results in a strain severely attenuated in HIA, with a phenotype nearly identical to that of BCG, requiring a 200-fold higher concentration of hemin to achieve growth equivalent to that of its parental strain. A nine-amino-acid deletion near the N terminus of BCG PPE37 (amino acids 31 to 39 of the M. tuberculosis PPE37 protein) underlies BCG's profound defect in HIA. Significant genetic variability exists in ppe37 genes across different M. tuberculosis strains, with more than 60% of sequences from completely sequenced M. tuberculosis genomes having mutations that result in altered PPE37 proteins; furthermore, these altered PPE37 proteins are nonfunctional in HIA. Our findings should allow delineation of the relative roles of HIA and SMIA in M. tuberculosis pathogenesis.