Comparative Genomics of a Bovine Mycobacterium tuberculosis Isolate and Other Strains Reveals Its Potential Mechanism of Bovine Adaptation

Genomic Characterization of IS6110 Insertions in Mycobacterium orygis

Mycobacterium orygis, a subspecies of the Mycobacterium tuberculosis complex (MTBC), has emerged as a significant concern in the context of One Health, with implications for zoonosis or zooanthroponosis or both. MTBC strains are characterized by the unique insertion element IS6110, which is widely used as a diagnostic marker. IS6110 transposition drives genetic modifications in MTBC, imparting genome plasticity and profound biological consequences. While IS6110 insertions are customarily found in the MTBC genomes, the evolutionary trajectory of strains seems to correlate with the number of IS6110 copies, indicating enhanced adaptability with increasing copy numbers. Here, we present a comprehensive analysis of IS6110 insertions in the M. orygis genome, utilizing ISMapper, and elucidate their genetic consequences in promoting successful host adaptation. Our study encompasses a panel of 67 paired-end reads, comprising 11 isolates from our laboratory and 56 sequences downloaded from public databases. Among these sequences, 91% exhibited high-copy, 4.5% low-copy, and 4.5% lacked IS6110 insertions. We identified 255 insertion loci, including 141 intragenic and 114 intergenic insertions. Most of these loci were either unique or shared among a limited number of isolates, potentially influencing strain behavior. Furthermore, we conducted gene ontology and pathway analysis, using eggNOG-mapper 5.0, on the protein sequences disrupted by IS6110 insertions, revealing 63 genes involved in diverse functions of Gene Ontology and 45 genes participating in various KEGG pathways. Our findings offer novel insights into IS6110 insertions, their preferential insertion regions, and their impact on metabolic processes and pathways, providing valuable knowledge on the genetic changes underpinning IS6110 transposition in M. orygis.

Identification of Differential Circular RNA Expression Profiles and Functional Networks in Human Macrophages Induced by Virulent and Avirulent Mycobacterium tuberculosis Strains

Circular RNAs (circRNAs) are noncoding RNAs with diverse functions. However, most Mycobacterium tuberculosis (M.tb)-related circRNAs remain undiscovered. In this study, we infected THP-1 cells with virulent and avirulent M.tb strains and then sequenced the cellular circRNAs. Bioinformatic analysis predicted 58,009 circRNAs in all the cells. In total, 2035 differentially expressed circRNAs were identified between the M.tb-infected and uninfected THP-1 cells and 1258 circRNAs were identified in the virulent and avirulent M.tb strains. Further, the top 10 circRNAs were confirmed by Sanger sequencing, among which four circRNAs, namely circSOD2, circCHSY1, circTNFRSF21, and circDHTKD1, which were highly differentially expressed in infected cells compared with those in uninfected cells, were further confirmed by ring formation, specific primers, and RNase R digestion. Next, circRNA-miRNA-mRNA subnetworks were constructed, such as circDHTKD1/miR-660-3p/IL-12B axis. Some of the individual downstream genes, such as miR-660-3p and IL-12B, were previously reported to be associated with cellular defense against pathological processes induced by M.tb infection. Because macrophages are important immune cells and the major host cells of M.tb, these findings provide novel ideas for exploring the M.tb pathogenesis and host defense by focusing on the regulation of circRNAs during M.tb infection.

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.

The ctpF Gene Encoding a Calcium P-Type ATPase of the Plasma Membrane Contributes to Full Virulence of Mycobacterium tuberculosis

Identification of alternative attenuation targets of Mycobacterium tuberculosis (Mtb) is pivotal for designing new candidates for live attenuated anti-tuberculosis (TB) vaccines. In this context, the CtpF P-type ATPase of Mtb is an interesting target; specifically, this plasma membrane enzyme is involved in calcium transporting and response to oxidative stress. We found that a mutant of MtbH37Rv lacking ctpF expression (MtbΔctpF) displayed impaired proliferation in mouse alveolar macrophages (MH-S) during in vitro infection. Further, the levels of tumor necrosis factor and interferon-gamma in MH-S cells infected with MtbΔctpF were similar to those of cells infected with the parental strain, suggesting preservation of the immunogenic capacity. In addition, BALB/c mice infected with Mtb∆ctpF showed median survival times of 84 days, while mice infected with MtbH37Rv survived 59 days, suggesting reduced virulence of the mutant strain. Interestingly, the expression levels of ctpF in a mouse model of latent TB were significantly higher than in a mouse model of progressive TB, indicating that ctpF is involved in Mtb persistence in the dormancy state. Finally, the possibility of complementary mechanisms that counteract deficiencies in Ca²⁺ transport mediated by P-type ATPases is suggested. Altogether, our results demonstrate that CtpF could be a potential target for Mtb attenuation.

Comparative Analysis on Proteomics Profiles of Intracellular and Extracellular M.tb and BCG From Infected Human Macrophages

Tuberculosis is the second cause in infectious diseases leading to human death. Understanding the virulence mechanism is inevitable if the disease needs to be fully cured. Therefore, this study aimed to reveal this mechanism by comparing proteomic profiles of intracellular and extracellular virulent strain M.tb and bacille Calmette–Guérin (BCG) from infected THP-1cells. First, M.tb and BCG infected THP-1 at MOI 10:1. Twelve hours postinfection, intracellular bacteria of M.tb and BCG were collected, whereas the two bacilli cultured in 7H9 broth media were used as the control. Then four groups of bacilli were subjected to proteomic analysis, and differential proteomic profiles between M.tb and BCG were comparatively analyzed with bioinformatics tools. As a result, we identified a total of 1,557 proteins. Further, they were divided into four groups for comparison of M.tb versus BCG under 7H9 culture (shorten as out), M.tb in (intracellular) versus M.tb out, BCG in versus BCG out and M.tb in versus BCG in. Between M.tb in versus BCG in, a total of 211 differentially expressed proteins were found. Eight proteins like ESAT-6 distributed in six RDs and some known proteins related to virulence. Besides, five uncharacterized proteins were differentially expressed. Further analysis revealed enriched pathways were associated with glyoxylate and dicarboxylate metabolism pathways. In M.tb out versus BCG out, a total of 144 differential proteins were identified and mainly involved in metabolism pathways. Then, 121 differential proteins in the group of M.tb in versus M.tb out were enriched in ribosome and oxidative phosphorylation related to adaptation to the host environment. The group of BCG in versus BCG out shared the same trend of different pathways to the M.tb in versus M.tb out. Finally, 42 proteins were identified to be up-regulated only in intracellular M.tb including eight RD proteins, whereas 22 up-regulated uniquely in intracellular BCG. Besides, only two proteins (Pks13 and Rv1405c) were commonly up-regulated in intracellular M.tb and BCG. Further, some unknown proteins were uniquely up-regulated in the intracellular M.tb and BCG. These findings provide valuable data for further exploration of molecular mechanism for M.tb virulence and BCG immune response.

Identification of Unique Key miRNAs, TFs, and mRNAs in Virulent MTB Infection Macrophages by Network Analysis

Although Mycobacterium tuberculosis (MTB) has existed for thousands of years, its immune escape mechanism remains obscure. Increasing evidence signifies that microRNAs (miRNAs) play pivotal roles in the progression of tuberculosis (TB). RNA sequencing was used to sequence miRNAs in human acute monocytic leukemia cells (THP-1) infected by the virulent MTB-1458 strain and the avirulent vaccine strain Mycobacterium bovis Bacillus Calmette-Guérin (BCG). Sets of differentially expressed miRNAs (DE-miRNAs) between MTB-1458/BCG-infected groups and uninfected groups were identified, among which 18 were differentially expressed only in the MTB-1458-infected THP-1 group. Then, 13 transcription factors (TFs) and 81 target genes of these 18 DE-miRNAs were matched. Gene Ontology classification as well as Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis showed that the candidate targets were predominantly involved in apoptotic-associated and interferon-γ-mediated signaling pathways. A TF-miRNA-mRNA interaction network was constructed to analyze the relationships among these 18 DE-miRNAs and their targets and TFs, as well as display the hub miRNAs, TFs, and target genes. Considering the degrees from network analysis and the reported functions, this study focused on the BHLHE40-miR-378d-BHLHE40 regulation axis and confirmed that BHLHE40 was a target of miR-378d. This cross-talk among DE-miRNAs, mRNAs, and TFs might be an important feature in TB, and the findings merited further study and provided new insights into immune defense and evasion underlying host-pathogen interactions.

Small RNA Profiling in Mycobacterium Provides Insights Into Stress Adaptability

Mycobacteria encounter a number of environmental changes during infection and respond using different mechanisms. Small RNA (sRNA) is a post-transcriptionally regulatory system for gene functions and has been investigated in many other bacteria. This study used Mycobacterium tuberculosis and Mycobacterium bovis Bacillus Calmette-Guérin (BCG) infection models and sequenced whole bacterial RNAs before and after host cell infection. A comparison of differentially expressed sRNAs using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) and target prediction was carried out. Six pathogenically relevant stress conditions, growth rate, and morphology were used to screen and identify sRNAs. From these data, a subset of sRNAs was differentially expressed in multiple infection groups and stress conditions. Many were found associated with lipid metabolism. Among them, ncBCG427 was significantly downregulated when BCG entered into macrophages and was associated with increased biofilm formation. The reduction of virulence possibility depends on regulating lipid metabolism.

MicroRNA-18b-5p Downregulation Favors Mycobacterium tuberculosis Clearance in Macrophages via HIF-1α by Promoting an Inflammatory Response

The modulation of the interaction between macrophages and Mycobacterium tuberculosis (M.tb) through microRNA during M.tb infection is increasingly capturing the attention of researchers. However, the potential role of microRNA-18b-5p (miR-18b) is not elucidated yet. In this study, miR-18b was found to be downregulated in M.tb-infected macrophage cell lines (THP-1 and RAW264.7) in time- and dose-dependent manners. Furthermore, when the miR-18b mimic and inhibitor and small interfering RNA hypoxia-inducible factor 1α (si-HIF-1α) were transfected into the macrophages separately or in combination, it was found that miR-18b targeted hypoxia-inducible factor 1α (HIF-1α). During M.tb infection, the decrease in the expression of miR-18b facilitated HIF-1α expression, which led to the increased production of pro-inflammatory cytokines, such as IL-6, resulting in decreased bacterial survival in the host cells. Moreover, the phosphorylation of p38 MAPK and NF-κB p65 was activated by the miR-18b inhibitor. Our findings expand the current understanding of the M.tb-cell interaction mechanism and provide a potential target to control M.tb infection.

Genetic differentiation of Mycobacterium bovis and Mycobacterium tuberculosis isolated from cattle and human sources in, Egypt (Suez Canal area)

Bovine tuberculosis is a devastating illness in cattle and it has the ability to transmit causing severe troubles in human. Mycobacterium bovis (M. bovis) infection in human indeed becomes increasingly critical especially in developing countries. Early diagnosis is very important to control and limit its spreading. The aim of this study is to examine the genetic differentiation and possibilities of transmission between cattle and human. Lymph node and sputum samples were collected from cattle and patients showing tuberculin test positive; respectively for phenotypic identification and for molecular examination by detection of IS6110 and oxyR genes which are specific for MTC and M. bovis; respectively. The phenotypic identification of sputum samples showed 80 % positive by both stain and culture, while, lymph nodes revealed 66 % and 84 % positive by stain and culture method; respectively. Alignment of oxyR gene sequences of M. tuberculosis and M. bovis was used as a feature for differentiation between the 2 genes in these two genetically closely similar microorganisms showed 99 % identities between the 2 genes. Alignment and phylogenetic analysis of Mpb70 gene sequences from animal and human origin showed very high relatedness (99.32 %) to each other confirming that the zoonotic transmission is most probably occurred.

TRIM25 upregulation by Mycobacterium tuberculosis infection promotes intracellular survival of M.tb in RAW264.7 cells

Tripartite motif 25 (TRIM25) is a TRIM family member which is involved in innate immunity. However, its role in the modulation of host defense against Mycobacterium tuberculosis (M.tb) infection has not been investigated. Therefore, this study aimed to demonstrate the significance of TRIM25 in the regulation of macrophage responses to M.tb infection. TRIM25 was found to be significantly overexpressed (3.476-fold) in peripheral blood mononuclear cells (PBMCs) of 67 patients with pulmonary tuberculosis compared with 48 healthy controls. TRIM25 expression was enhanced following M.tb infection of RAW264.7 cells, a macrophage cell line. Overexpression of TRIM25 in M.tb-infected RAW264.7 cells led to a significant increase in phosphorylated p38 levels; however, the production of IL-6, IL-1β, and TNF-α were significantly reduced. Finally, M.tb intracellular survival increased by 90% at 12 h post-infection (PI) (p < 0.01). To validate the previous results, TRIM25 levels in M.tb-infected RAW264.7 macrophages were down-regulated using small interfering RNA (siRNA). Therefore, it was concluded that TRIM25 promotes intracellular survival of M.tb in RAW264.7 cells, likely by enhancing p38 pathways and thereby inhibiting the production of proinflammatory cytokines. These results contribute to the further understanding of the host defense against M.tb infection.

Comparative Study of the Molecular Basis of Pathogenicity of M. bovis Strains in a Mouse Model

It is widely accepted that different strains of Mycobacterium tuberculosis have variable degrees of pathogenicity and induce different immune responses in infected hosts. Similarly, different strains of Mycobacterium bovis have been identified but there is a lack of information regarding the degree of pathogenicity of these strains and their ability to provoke host immune responses. Therefore, in the current study, we used a mouse model to evaluate various factors involved in the severity of disease progression and the induction of immune responses by two strains of M. bovis isolated from cattle. Mice were infected with both strains of M. bovis at different colony-forming unit (CFU) via inhalation. Gross and histological findings revealed more severe lesions in the lung and spleen of mice infected with M. bovis N strain than those infected with M. bovis C68004 strain. In addition, high levels of interferon-γ (IFN-γ), interleukin-17 (IL-17), and IL-22 production were observed in the serum samples of mice infected with M. bovis N strain. Comparative genomic analysis showed the existence of 750 single nucleotide polymorphisms and 145 small insertions/deletions between the two strains. After matching with the Virulence Factors Database, mutations were found in 29 genes, which relate to 17 virulence factors. Moreover, we found an increased number of virulent factors in M. bovis N strain as compared to M. bovis C68004 strain. Taken together, our data reveal that variation in the level of pathogenicity is due to the mutation in the virulence factors of M. bovis N strain. Therefore, a better understanding of the mechanisms of mutation in the virulence factors will ultimately contribute to the development of new strategies for the control of M. bovis infection.