Exploring virulence in Mycobacterium bovis: clues from comparative genomics and perspectives for the future
Irish Veterinary Journal volume 76, Article number: 26 (2023)
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.
Mycobacterium bovis is a member of the Mycobacterium tuberculosis complex (MTBC), the group of related pathogens that causes tuberculosis (TB) in mammals . A key feature of the MTBC pathogens is their adaptation to sustain in specific host populations: by this, we mean the ability of these pathogens to maintain cycles of infection, disease, and transmission within their preferred hosts. Mycobacterium bovis is seen as the archetypal animal-adapted member of the MTBC, able to sustain across cattle, badgers, deer, possums, etc. This contrasts with Mycobacterium tuberculosis which is the prototypic human TB pathogen but with apparently limited ability to sustain in domesticated animals or wildlife. The high degree of genetic identity between M. bovis and M. tuberculosis therefore offers the possibility of using comparative genomic approaches to identify the molecular basis of how M. bovis has evolved into such an effective, and problematic, animal pathogen. Through the identification of such pathogen adaptations, and understanding their mechanisms of action, we hope to better thwart M. bovis and ultimately accelerate eradication of bovine TB.
It was Theobald Smith in 1896 who first described morphological, culture and pathological distinctions between TB bacilli isolated from humans and animals, leading to his suggestion of a distinct bovine TB bacillus [2, 3]. Through this work Smith showed how inoculation of cattle with bacilli isolated from either human sputum, or from cattle, led to distinct pathological outcomes in the infected cattle, with the bovine-derived bacilli showing greater virulence in cattle than human sputum-derived bacilli. The announcement in 1901 by Robert Koch at the 1901 British Congress on Tuberculosis that the bovine tubercle bacillus presented a negligible risk to humans, and that he did “not deem it advisable to take any measures against it”  sparked a major controversy. Multiple research studies were established that would subsequently show the risk of infection of the bovine bacillus to humans. A British Royal Commission established in 1901 to examine the evidence concluded in 1911 that humans “must therefore be added to the list of animals notably susceptible to bovine tubercle bacilli” . While the bovine bacillus was often seen as a subtype of M. tuberculosis and hence referred to as M. tuberculosis var bovis or M. tuberculosis subsp bovis, the designation Mycobacterium bovis was established in 1970 .
It is worth returning to the point above relating to the early comparative virulence work on human and bovine bacilli. In addition to Smith’s work showing the attenuation of human-derived bacilli for cattle as compared to bovine varieties, von Behring also promoted the use of human bacilli as a vaccine against TB in cattle, termed ‘bovo-vaccine’. This was based on von Behring’s own work showing the attenuation of human bacilli for cattle [7, 8]. However, subsequent work was to show that inoculation of cattle with bovo-vaccine led to shedding of live human tubercle bacilli in milk , while others reported variable protection afforded by the vaccine , and the use of bovo-vaccine was ultimately abandoned. As described by Palmer and Waters  in their review of the US bovine TB eradication program, reports of variable virulence levels of human bacilli for cattle were often not clear cut. In his 1901 report as Chair of the Bureau of Animal Industry, Salmon  referred to work from Chauveau showing that “human tubercular virus (sic) acts on the bovine species exactly like the tubercular virus which comes from the bovine species itself”. In this same report Salmon also described work from Ravanel who fed 4 calves with human sputum, and that “postmortem examination proved that all had become affected with tuberculosis, the lesions in two being quite extensive.” Ravenel had also obtained bacilli from a child who died of TB meningitis; using this for infection studies he found that “bacilli from this material have proved very virulent for bovine animals” . However, the exact designation of bacilli used in these early studies is often unclear, inoculation dose and routes varied, and immune status of the cattle used in the infection studies not stated.
To readdress virulence of M. bovis and M. tuberculosis in cattle, Whelan et al.  undertook experimental infections using M. tuberculosis H37Rv and M. bovis AF2122/97, which at the time were the best characterised genome sequenced isolates of the respective species. In their work, cattle infected with M. tuberculosis H37Rv and M. bovis AF2122/97 both gave positive skin test and interferon gamma tests, but cattle infected with M. bovis AF2122/97 showed extensive pathology while those infected with M. tuberculosis H37Rv showed no obvious pathological lesions. Of course, it could be argued that M. tuberculosis H37Rv is a laboratory-adapted strain; it was derived by Steenken and colleagues in the 1930s by in vitro passage of an M. tuberculosis strain termed ‘H37’ which itself had been isolated by Baldwin in 1905 . Hence the apparent attenuation of M. tuberculosis H37Rv for cattle could be an artefact of its in vitro adaptation and not be representative of M. tuberculosis in general. To further explore the virulence of M. tuberculosis in experimental cattle infections, a strain of M. tuberculosis termed BTB1558, that had been isolated from the lesion of a Zebu bull in Ethiopia, was assessed in head-to-head comparisons with M. bovis AF2122/97 and M. tuberculosis H37Rv . In this latter study, Villarreal-Ramos and colleagues found that M. bovis AF2122/97 generated more pathology in infected animals than either M. tuberculosis H37Rv or the M. tuberculosis BTB1558 strain. Hence both experimental infection studies have shown that the M. tuberculosis strains used showed reduced virulence as compared to M. bovis AF2122/97 in cattle.
Before we draw definitive conclusions from these experiments on the attenuation of M. tuberculosis in the bovine host per se, it should be noted that both M. tuberculosis strains used in the studies of Whelan et al.  and Villarreal-Ramos et al.  were from M. tuberculosis lineage 4, the so-called Euro-American lineage. There are nine currently defined lineages of M. tuberculosis, with many of these showing distinct geographical distributions [16, 17]. Reports of the isolation of M. tuberculosis from lineage 1 , 3 and 4  from cattle in India suggests that conclusions derived from the study of M. tuberculosis BTB1558 and H37Rv may not apply across all M. tuberculosis lineages for cattle. Also, the immune status of an animal could modulate its susceptibility to infection and disease regardless of the genetic background of the infecting strain. Furthermore, virulence is not a ‘linear’ concept, and pathogen success is not necessarily defined by increased virulence for the host population to which the pathogen is adapted. Nevertheless, M. tuberculosis H37Rv and M. bovis AF2122/97 can provide ‘extremes’ in terms of the outcome of infection in cattle, with comparative analysis of these strains providing a route to define the molecular basis of M. bovis adaptation to cattle. Furthermore, the isolation of Mycobacterium orygis from cattle in South Asia [20,21,22] suggests that this MTBC species has also evolved a successful strategy to sustain in cattle populations , suggesting that comparative analysis using this species will also prove fruitful.
Mycobacterium bovis genomics
The landmark publication of the genome sequence of M. tuberculosis H37Rv in 1998  heralded a new era in our ability to explore the internal workings of mycobacterial pathogens. The M. bovis sequencing project came to fruition a few years later with the publication in 2003 of the complete genome sequence of M. bovis AF2122/97, a strain which had been isolated in Cornwall, UK in 1997 . These genome sequences were elucidated using Sanger sequencing technology. The advent of ‘next generation’ sequencing approaches in the mid -2000s meant that mycobacterial pathogens can now be rapidly sequenced at a fraction of the cost of the original genome sequencing projects. The exquisite resolution and throughput provided by next generation sequencing has transformed the ability to trace the evolution of M. bovis  and track transmission across livestock and wildlife . As other presentations at the M. bovis 2022 conference showcased the application of pathogen genome sequencing technology for epidemiological studies, we shall not deal with it in this review. Rather we will focus on what genomics has taught us about the functional capacity of the M. bovis genome.
As the first approach to the analysis of what makes M. bovis ‘tick’ as a successful animal pathogen, comparative analysis with the genome of M. tuberculosis strains, in particular M. tuberculosis H37Rv, was pursued . This analysis revealed > 99% identity at the nucleotide level between the bacilli, as well as confirming the presence of several deletions from the M. bovis genome as compared to M. tuberculosis [25, 28]. Notably the RD4 locus was deleted in all isolates of M. bovis studied, hence identifying it as a useful target locus in developing simple PCR-based screens to differentiate M. bovis from other members of the MTBC [29, 30]. Another immediate finding of the genomic studies was that there were no unique genes per se in the genome of M. bovis compared to M. tuberculosis . This latter finding indicated that the adaptation(s) of M. bovis to sustain in animal populations were not due to the presence of a particular gene that was lacking in human adapted bacilli. Instead, differential expression of genes across the human and bovine TB bacilli appeared a more likely basis for their unique host adaptations. Studies identifying differentially expressed genes and proteins between the animal- and human-adapted bacilli were therefore pursued as a route to defining their functional impact on host tropism.
An example of the impact of single nucleotide polymorphisms (SNPs) between the genomes of M. bovis and M. tuberculosis was the identification of a SNP in the pykA gene encoding pyruvate kinase, a key gene in the metabolism of carbohydrates that converts phosphoenolpyruvate to pyruvate which feeds into the TCA cycle. The M. bovis SNP leads to the inactivation of pyruvate kinase due to substitution of the glutamic acid residue at position 220, located at the enzyme’s active site, with aspartic acid . The lack of a functional pyruvate kinase means M. bovis cannot use carbohydrates as a sole carbon source, an insight that provided an explanation for the long-held observation that M. bovis is unable to grow on glycerol as a sole carbon source. Studies have shown that inactivation of the pykA of M. tuberculosis resulted in catabolism of fatty acids for energy production, whereas complementation of M. bovis with a functional pykA led to carbon conservation through gluconeogenesis and the glyoxylate cycle . Whether the inactivation of the M. bovis pykA is deleterious for the bacillus, or offers some advantage via metabolic remodelling, is unknown. Interestingly, analysis of an expanding M. bovis clone in the UK (termed ‘type 17’) revealed a deletion that altered carbon metabolism such that propionate derived from fatty acid metabolism was shuttled to generate pyruvate, hence perhaps compensating for the inactivation of PykA .
With the advent of ‘Next Generation’ sequencing technologies M. bovis AF2122/97 was re-sequenced using both Illumina and PacBio approaches . A surprising finding from this resequencing was the presence of an approximately 3.1 kb region of DNA which had been missed in the original M. bovis AF2122/97 sequencing project. This 3.1 kb region corresponded to the RD900 locus which had been described by Bentley and colleagues as present in Mycobacterium africanum but deleted from many other members of the MTBC, including M. bovis . However, the resequencing and de novo assembly showed that this locus was in fact present in the genome of M. bovis AF2122/97. RD900 contains genes encoding a potential transporter and two protein kinases (PknH1 and PknH2); these latter proteins are through to transduce extracellular signals and trigger signalling cascades. The genome of M. tuberculosis H37Rv lacks the RD900 locus and has a single PknH . Subsequent work by Mata and colleagues  showed that expression of the M. tuberculosis PknH orthologue in M. bovis led to altered expression of 154 genes in M. bovis, indicating that the M. tuberculosis and M. bovis PknH protein kinases have distinct functional activities.
To update the functional annotation of the M. bovis genome, Malone et al.  used a combination of quantitative transcriptomics and proteomics to identify the most highly expressed genes and most abundant proteins in M. bovis AF2122/97 and compared this to parallel datasets from M. tuberculosis H37Rv. This revealed 77 genes upregulated in M. bovis AF2122/97 at both the RNA and protein level as compared to M. tuberculosis H37Rv . To explore the mechanistic basis for differential gene expression analysis across the two species, a transcription factor enrichment analysis revealed significant association of a core set of transcription factors with the differentially expressed genes. Among these transcription factors was SigK and PhoPR which we discuss in greater detail below.
A further update to the M. bovis AF2122/97 genome annotation was performed by Farrell and colleagues . Their analysis focused on genes that had originally been described as ‘hypothetical’, or of unknown function. These latter designations represented a limitation of the initial genome annotation, namely a reliance on ascribing potential gene function based on information available from other organisms. Due to the paucity of comparative data available at the original time of genome annotation in the early 2000s, almost 30% of coding sequences in a genome would have no known function and comparison to other organisms would reveal limited information. However, with the exponential increase in bacterial genomes sequenced and analysed using ‘omic methodologies over the past 20 years, plus the advent of more sophisticated computational tools for comparative analysis of orthologous genes, many M. bovis hypothetical proteins could now be placed into potential functional groupings. Through this analysis, Farrell et al.  updated functional information for over 600 M. bovis genes and updated 69 gene names. The current M. bovis AF2122/97 genome and functional annotation is available via the DDBJ/ENA/GenBank accession number LT708304.
As highlighted above, the transcriptomic and proteomic information layered onto the genome sequence identified key differentially regulated loci between M. bovis AF2122/97 and M. tuberculosis H37Rv. Given the divergent interaction of these species with cattle, the genes and proteins identified as differentially expressed provide points a focus for research efforts that seek to refine our understanding of M. bovis virulence and host adaptation. In the following sections we focus on regions that show differential regulation in M. bovis and summarise what ongoing studies are revealing about the potential functional impacts of these differences.
The SigK regulon
SigK is an extra cytoplasmic function (ECF) sigma factor which, as the name suggests, regulates responses to extracellular conditions. Behr and colleagues were the first to show that the SigK regulon was dysregulated between M. bovis and M. tuberculosis [39, 40]. They showed that the anti-sigma factor regulator of SigK, RskA, harboured a mutation in M. bovis that inactivated its function, leading to constitutive expression of the SigK regulon in M. bovis [39, 41]. Amongst the genes under the control of SigK are mpb70 and mpb83, encoding the proteins MPB70 and MPB83, respectively (Fig. 1). The ‘MPB’ nomenclature was introduced by Nagai et al.  with, for example, MPB70 standing for “a mycobacterin which is a protein fraction from M. bovis BCG as a band with a relative mobility of 0.70”. While MPB70 is abundant in culture filtrates of M. bovis it was not detected within sonicated extracts from the bacteria prewashed with aqueous buffers, indicating efficient secretion of the protein . In contrast, MPB83 was detected in sonicated lysates of washed bacilli, and flow cytometry of M. bovis BCG Tokyo with a monoclonal antibody against MPB83 showed that MPB83 was anchored at the bacterial surface .
It is noteworthy that upregulated expression of the SigK regulon was also seen in M. orygis, but by an independent mutation . This is intriguing because M. orygis has been suggested to be another bovine-adapted member of the MTBC [20, 21, 23], implying a potentially selective advantage to upregulation of the SigK regulon in pathogens that sustain in bovine hosts. Furthermore, it was demonstrated  that RskA is not a simple negative regulator that sequesters its associated sigma factor SigK and releases it upon environmental stimulation. Instead, RskA possesses a dual function of activator-inhibitor of SigK. Both the M. bovis and M. orygis mutated forms of RskA do not result in increased SigK activity simply due to loss of the RskA inhibitory function, but rather by enhancement of the RskA activator function, hence leading to increased SigK activity .
The genes under the control of SigK are restricted to two loci on the genome [39, 40, 45]. One region is the mpb70/83 locus that encodes: MPB83, a homolog of MPB70 that is anchored to the membrane and post-translationally modified by acylation and glycosylation; DipZ, an integral membrane protein; MPB70, a secreted protein; and Mb2901, a possible transmembrane protein. The second region is the SigK locus itself, which comprises: Mb0457c and Mb0456c, proteins of unknown function; UfaA1, a phospholipid synthase; Mb0453c, again of unknown function; and finally, SigK and RskA . While ECF sigma factors are known to control relatively small regulons, the SigK regulon is remarkably small; in comparison to the 10 genes encompassed by the SigK regulon, SigH regulates the expression of 41 genes , while SigC controls the expression of 18 genes .
The best characterised member of the SigK regulon is MPB70. The structure of mature MPB70 was determined by nuclear magnetic resonance (NMR) spectroscopy in 2003 , revealing that MPB70 folds into a globular tertiary structure made of a seven stranded β-barrel adjoined to eight α-helices that pack on one side of the barrel. This study unveiled significant similarity between the backbone topology of residues 27 to 159 of MPB70 (i.e., 133 amino acids of the total 163 amino acid sequence of MPB70) and the tertiary structures of the third and fourth FAS1 (Fasciclin-like) domains of the Drosophila melanogaster Fasciclin I protein. This latter protein was the first for which the structure of a FAS1 domain was determined, using X-ray crystallography, and therefore provided the structural prototype for the FAS1 domain . In general, FAS1 domains are made of approximately 140 amino acids and are present in extracellular proteins, either secreted or bound to the plasma membrane, where they are implicated in cell-to-cell, cell-to-extracellular matrix, and extracellular matrix structure adhesion , or in homophilic domain-to-domain trans-interactions between FAS1 proteins to mediate cell-to-cell interactions . The FAS1 domain is an ancient structural motif, as it is shown to be extensively represented in animal, bacteria, and plant proteins [52, 53]. FAS1 proteins were found in both eubacteria and archaea, suggesting that the domain originated prior to the last universal common ancestor, making it an ancient adhesive extracellular domain . While FAS1 proteins are widely represented throughout various—if not all—biological groups, they are not ubiquitous (especially in microorganisms), suggesting that the FAS1 domain takes part in specialized cellular interactions .
A range of human proteins contain FAS1 domains, including periostin, TGFBI/βig-h3, stabilin-1 and stabilin-2. The latter two stabilin proteins are homologous transmembrane proteins with C-terminal cytosolic tails . The stabilins contain a total of seven FAS1 domains, composed of three dual FAS1 domain tandems and one isolated FAS1 domain, each separated by a varying number of Epidermal Growth Factor (EGF) repeats and an X-link domain, with these latter domains thought to be involved in phosphatidylserine and hyaluronic acid binding, respectively [50, 55]. Both stabilin-1 and stabilin-2 are expressed in various cell types, including tissue macrophages. Recombinant mammalian cell lines that had been engineered to express high levels of stabilin-2 were shown to become sticky and aggregate, with the FAS1 domains of stabilin-2 found to be responsible for homophilic cell–cell adhesion . This work also revealed an intriguing characteristic of MPB70. Briefly, Park et al.  performed a series of aggregation assays with a mouse fibroblast cell line expressing (or not) stabilin-2 and firstly demonstrated that stabilin-2 mediated cell aggregation in suspension. Then, to find out if the FAS1 domains of stabilin-2 in particular were responsible for cell aggregation, the authors used recombinant MPB70 – used as a model FAS1 protein carrying a unique FAS1 domain—as a further element in their aggregation assay. The authors expected that the FAS1 domains of MPB70 would inhibit cell aggregation by interfering with homophilic interactions of stabilin-2. However, the opposite happened: when MPB70 was present in the media it further enhanced the aggregation of the cells. This suggested a role of MPB70 in cellular aggregation.
MPB70: host interactions
Queval et al.  explored the interaction of M. bovis and M. tuberculosis with human and bovine macrophages. They observed that the in vitro infection of bovine blood monocyte-derived macrophages (MDMs) with M. bovis, but not M. tuberculosis, promoted the formation of multinucleated cells (MNC) by membrane fusion. Furthermore, they only obtained the multinucleation phenotype when using bovine MDMs; infection of human MDMs with M. bovis did not trigger multinucleation, suggesting a species-specific interaction. Multinucleation is a complex and multi-step process which, amongst other determinants, requires the interaction and fusion of cellular membranes. As seen previously, stabilin-2 facilitated cellular membrane interaction leading to aggregation and played a key role during myoblast multinucleation [51, 57, 58]. The observation by Park and colleagues that MPB70 enhanced cell aggregation led to the hypothesis that the differential phenotype of multinucleation of bovine MDMs could be a consequence of the differential profile of expression of MPB70 by the two mycobacterial strains (as M. bovis has constitutive expression of MPB70, while MPB70 expression is induced in M. tuberculosis after phagocytosis).
To address this hypothesis, Queval et al.  assessed whether a lack of MPB70 in M. bovis would prevent the macrophage multinucleation phenotype they had seen. An in vitro infection of bovine MDMs was therefore performed with using an M. bovis AF2122/97 wild type strain, an M. bovis AF2122/97 Δmpb70 knockout mutant, and the complemented M. bovis AF2122/97 Δmpb70 knockout with synthesis of MPB70 restored. At 24 h post infection lower numbers of MNC were obtained after infection of bovine MDMs with the Δmpb70 knockout strain, similar to numbers obtained in the non-infected cells. They furthermore demonstrated that infection of bovine macrophages with the complemented M. bovis mpb70 knockout strain led to similar MNC levels to those observed after infection with M. bovis wild type. These results therefore supported the hypothesis that MPB70 played a role in the phenotype of multinucleation of bovine MDMs.
The observation of M. bovis-specific interaction with bovine immune cells, dependent on MPB70, are intriguing. However, many questions remain. For example, Queval et al.  studied cellular phenotypes over 24 h; it is possible that the MPB70 effect may be temporal, specific to early time points in vitro, and that later time points may reveal multinucleation for both species. At the cellular level, we do not know what the macrophage proteins are that interact with MPB70. While stabilin-2 or other FAS-1 containing proteins may be candidates, this needs to be clarified. The in vitro observation of multinucleation was made with bovine blood MDMs; however, in vivo, tissue-resident alveolar macrophages are of a different origin and are functionally distinct to recruited, interstitial macrophages. It would therefore be interesting to see what, if any, differential phenotype would be observed in bovine alveolar macrophages infected with the M. bovis wild type, the mpb70 knockout mutant, and complemented mpb70 mutant. Experimental infection of cattle with the M. bovis AF2122/97 Δmpb70 knockout would also reveal how the phenotypes observed in vitro translate to the in vivo situation. Indeed, it remains to be established if MPB70 is a lone player in host interactions, or if other genes in the SigK regulon also play a role.
Another system that showed differential expression between M. bovis and M. tuberculosis was the PhoPR system (Fig. 1). PhoPR is a master regulator of virulence in M. tuberculosis. It is a two-component system, made up of a sensor kinase (PhoR) that detects changes in environmental stimuli, and a response regulator (PhoP). When activated by PhoR, PhoP forms a dimeric structure via its N-terminal region, and its C-terminal binding domain directly binds target sequences across more than 30 M. tuberculosis genes to regulate their expression [59, 60].
A major virulence system that is controlled by PhoPR is the ESX-1 secretion system which exports, among other factors, the dominant T-cell antigens ESAT-6 (EsxA) and CFP-10 (EsxB), and the ESX-1 secretion-associated proteins (Esp) EspA and EspC . The latter proteins are encoded in an operon, espACD, which is located distally from the main ESX-1 locus. Secretion from the ESX-1 locus is co-dependent on both the ESX-1 and espACD loci, with PhoP binding upstream of espACD to regulate ESX-1 secretion . The PhoPR system also regulates the synthesis of many lipids restricted to pathogenic members of the MTBC such as diacyltrehaloses (DAT), polyacyltrehaloses (PAT), and sulfolipids (SL) [62, 63]. The major binding site of the PhoP regulator is the small non-coding RNA (ncRNA) encoded by mcr7 that modulates the activity of the Twin Arginine Translocation (TAT) secretion apparatus. In PhoP mutants mcr7 expression is abolished, which in turn affects the secretion of TAT substrates, such as the Ag85 complex . The PhoPR system is hence a critical system for the control of a range of virulence functions in M. tuberculosis.
Given the central role of PhoPR in controlling virulence gene expression, it was a surprise when Gonzalo-Asensio et al.  reported defects in the PhoPR system in M. bovis. Specifically, they showed that M. bovis PhoR harbours a mutation at codon 71 that replaces a glycine with an isoleucine (Gly71Ile), leading to a defective PhoPR system. The apparent loss of a functional PhoPR in M. bovis compared to M. tuberculosis explains some of the known differences across the bacilli; for example, the absence of sulfolipid in M. bovis compared to M. tuberculosis. However, if the PhoPR system controls the secretion of ESAT-6 and CFP-10, antigens that are known to be secreted by M. bovis, how does the bovine pathogen still secrete these proteins? To explain this Gonzalo-Asensio and colleagues showed that deletion of the RD8 region from M. bovis removed espACD upstream control regions which, in combination with SNPs in the espACD promoter, allow espACD expression and hence ESX-1 secretion . This again shows how study of genetic differences between M. bovis and M. tuberculosis can reveal the molecular basis for functional differences between the bacilli.
While the PhoPR system shows differences between M. bovis and M. tuberculosis, it is not the case that PhoPR is non-functional in M. bovis. Garcia et al.  showed that inactivation of PhoP in M. bovis led to decreased ability of the mutant to block phagosome maturation after engulfment by macrophages. Through microarray analysis of global gene expression, they also showed that the M. bovis phoP mutant had altered expression of 70 genes as compared to the M. bovis wild type, including genes involved in sulfolipid biosynthesis and transport, oxidative stress responses, and ESX-1 secretion. The M. bovis phoP mutant was also more sensitive to acid stress than the wild type, a stress condition that mimics the in vivo environment. Subsequent work by Garcia et al.  revealed that PhoP controls the production and secretion of ammonia from M. bovis that buffers against acid stress. Our own work (Urtasun-Elizari and Gordon, unpublished) has also found that inactivation of the phoPR locus in M. bovis leads to altered expression of multiple genes, again indicating that the locus is functional in M. bovis. However, the M. bovis PhoR orthologue may respond to distinct signals to that of M. tuberculosis.
An organism’s genome is often referred to as its ‘genetic blueprint’, as if it were a detailed plan to build the organism from scratch. However, the analogy is flawed; the genome is better viewed as a parts list, and a partial one at that, with no instructions on how the various parts come together to build a functioning organism. Through analysis of the M. bovis genome we aim to translate the genome information into an understanding of what makes the pathogen ‘tick’ in terms of its ability to infect and transmit across animal populations. As described herein, insights gleaned to date include an understanding of how evolution has shaped the genome of M. bovis to be distinct from that of the human-adapted pathogen M. tuberculosis, altering the expression of multiple genes that are involved in host–pathogen interaction. These findings have provided targets for us to focus our efforts on elucidating the mechanistic basis of pathogenesis. Such ‘omics studies parallel the use of genomics to define M. bovis population structure and transmission dynamics, to identify novel antigens, and to develop rapid nucleic acid diagnostics. Genomics is now a cornerstone of M. bovis research. We look forward to the next international M. bovis meeting to see genomics research translated to improved bovine TB control.
Availability of data and materials
Mycobacterium tuberculosis Complex
Single nucleotide polymorphism
Malone KM, Gordon SV. Mycobacterium tuberculosis Complex Members Adapted to Wild and Domestic Animals. Adv Exp Med Biol. 2017;1019:135–54. https://doi.org/10.1007/978-3-319-64371-7_7.
Smith T. A Comparative Study of Bovine Tubercle Bacilli and of Human Bacilli from Sputum. J Exp Med. 1898;3(4–5):451–511. https://doi.org/10.1084/jem.3.4-5.451.
Smith T. Two varieties of the tubercle bacillus from mammals. Trans Assn Am Physicians. 1896;11:75–93.
Koch R. An Address on the Fight against Tuberculosis in the Light of the Experience that has been Gained in the Successful Combat of other Infectious Diseases. Br Med J. 1901;2(2117):189–93. https://doi.org/10.1136/bmj.2.2117.189.
Francis J. The work of the British Royal Commission on Tuberculosis, 1901–1911. Tubercle. 1959;40(2):124–32. https://doi.org/10.1016/s0041-3879(59)80153-7.
Karlson AG, Lessel EF. Mycobactrium bovis nom. nov. Int J Syst Evol Microbiol. 1970;20(3):273–82.
Behring E. Nobel Lecture: Serum therapy in therapeutics and medical science. 1901.
Morris CD. Dr Von Behring’s Bovovaccine as an immunising virus. Vet J. 1900;62(3):154–6.
Griffith AS. Human tubercle bacilli in the milk of a vaccinated cow. J Pathol Bacteriol. 1912;17(3):323–8.
Prausnitz C. Bacteriological Notes. J R Inst Public Health. 1907;15(11):681–3.
Palmer MV, Waters WR. Bovine tuberculosis and the establishment of an eradication program in the United States: role of veterinarians. Vet Med Int. 2011;2011:816345. https://doi.org/10.4061/2011/816345.
Salmon DE: Relation of Bovine Tuberculosis to Public Health. Bureau of Animal Industry Bulletin No. 33. In. Edited by Agriculture USDo, vol. 33. Washington: Goverment Printing Office; 1901.
Whelan AO, Coad M, Cockle PJ, Hewinson G, Vordermeier M, Gordon SV. Revisiting host preference in the Mycobacterium tuberculosis complex: experimental infection shows M. tuberculosis H37Rv to be avirulent in cattle. PLoS One. 2010;5(1):e8527. https://doi.org/10.1371/journal.pone.0008527.
Steenken W, Oatway WH, Petroff SA. Biological Studies of the Tubercle Bacillus : Iii. Dissociation and Pathogenicity of the R and S Variants of the Human Tubercle Bacillus (H(37)). J Exp Med. 1934;60(4):515–40. https://doi.org/10.1084/jem.60.4.515.
Villarreal-Ramos B, Berg S, Whelan A, Holbert S, Carreras F, Salguero FJ, et al. Experimental infection of cattle with Mycobacterium tuberculosis isolates shows the attenuation of the human tubercle bacillus for cattle. Sci Rep. 2018;8(1):894. https://doi.org/10.1038/s41598-017-18575-5.
Ngabonziza JCS, Loiseau C, Marceau M, Jouet A, Menardo F, Tzfadia O, et al. A sister lineage of the Mycobacterium tuberculosis complex discovered in the African Great Lakes region. Nat Commun. 2020;11(1):2917. https://doi.org/10.1038/s41467-020-16626-6.
Coscolla M, Gagneux S, Menardo F, Loiseau C, Ruiz-Rodriguez P, Borrell S, et al. Phylogenomics of Mycobacterium africanum reveals a new lineage and a complex evolutionary history. Microb Genom. 2021;7(2); https://doi.org/10.1099/mgen.0.000477.
Palaniyandi K, Kumar N, Veerasamy M, Kabir Refaya A, Dolla C, Balaji S, et al. Isolation and comparative genomics of Mycobacterium tuberculosis isolates from cattle and their attendants in South India. Sci Rep. 2019;9(1):17892. https://doi.org/10.1038/s41598-019-54268-x.
Mukherjee F, Bahekar VS, Pasha SY, Kannan P, Prasad A, Rana SK, et al. Isolation and analysis of the molecular epidemiology and zoonotic significance of Mycobacterium tuberculosis in domestic and wildlife ruminants from three states in India. Rev Sci Tech. 2018;37(3):999–1012. https://doi.org/10.20506/rst.37.3.2902.
Duffy SC, Srinivasan S, Schilling MA, Stuber T, Danchuk SN, Michael JS, et al. Reconsidering Mycobacterium bovis as a proxy for zoonotic tuberculosis: a molecular epidemiological surveillance study. Lancet Microbe. 2020;1(2):e66–73. https://doi.org/10.1016/S2666-5247(20)30038-0.
Rahim Z, Thapa J, Fukushima Y, van der Zanden AGM, Gordon SV, Suzuki Y, et al. Tuberculosis Caused by Mycobacterium orygis in Dairy Cattle and Captured Monkeys in Bangladesh: a New Scenario of Tuberculosis in South Asia. Transbound Emerg Dis. 2017;64(6):1965–9. https://doi.org/10.1111/tbed.12596.
Refaya AK, Kumar N, Raj D, Veerasamy M, Balaji S, Shanmugam S, et al. Whole-Genome Sequencing of a Mycobacterium orygis Strain Isolated from Cattle in Chennai, India. Microbiol Resour Announc. 2019;8(40); https://doi.org/10.1128/MRA.01080-19.
Brites D, Loiseau C, Menardo F, Borrell S, Boniotti MB, Warren R, et al. A New Phylogenetic Framework for the Animal-Adapted Mycobacterium tuberculosis Complex. Front Microbiol. 2018;9:2820. https://doi.org/10.3389/fmicb.2018.02820.
Cole ST, Brosch R, Parkhill J, Garnier T, Churcher C, Harris D, et al. Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature. 1998;393(6685):537–44. https://doi.org/10.1038/31159.
Garnier T, Eiglmeier K, Camus JC, Medina N, Mansoor H, Pryor M, et al. The complete genome sequence of Mycobacterium bovis. Proc Natl Acad Sci U S A. 2003;100(13):7877–82. https://doi.org/10.1073/pnas.1130426100.
Loiseau C, Menardo F, Aseffa A, Hailu E, Gumi B, Ameni G, et al. An African origin for Mycobacterium bovis. Evol Med Public Health. 2020;2020(1):49–59. https://doi.org/10.1093/emph/eoaa005.
Kao RR, Price-Carter M, Robbe-Austerman S. Use of genomics to track bovine tuberculosis transmission. Rev Sci Tech. 2016;35(1):241–58. https://doi.org/10.20506/rst.35.1.2430.
Gordon SV, Brosch R, Billault A, Garnier T, Eiglmeier K, Cole ST. Identification of variable regions in the genomes of tubercle bacilli using bacterial artificial chromosome arrays. Mol Microbiol. 1999;32(3):643–55. https://doi.org/10.1046/j.1365-2958.1999.01383.x.
Taylor GM, Worth DR, Palmer S, Jahans K, Hewinson RG. Rapid detection of Mycobacterium bovis DNA in cattle lymph nodes with visible lesions using PCR. BMC Vet Res. 2007;3:12. https://doi.org/10.1186/1746-6148-3-12.
Warren RM, van GeyPittius NC, Barnard M, Hesseling A, Engelke E, de Kock M, et al. Differentiation of Mycobacterium tuberculosis complex by PCR amplification of genomic regions of difference. Int J Tuberc Lung Dis. 2006;10(7):818–22.
Keating LA, Wheeler PR, Mansoor H, Inwald JK, Dale J, Hewinson RG, et al. The pyruvate requirement of some members of the Mycobacterium tuberculosis complex is due to an inactive pyruvate kinase: implications for in vivo growth. Mol Microbiol. 2005;56(1):163–74. https://doi.org/10.1111/j.1365-2958.2005.04524.x.
Chavadi S, Wooff E, Coldham NG, Sritharan M, Hewinson RG, Gordon SV, et al. Global effects of inactivation of the pyruvate kinase gene in the Mycobacterium tuberculosis complex. J Bacteriol. 2009;191(24):7545–53. https://doi.org/10.1128/JB.00619-09.
Wheeler PR, Brosch R, Coldham NG, Inwald JK, Hewinson RG, Gordon SV. Functional analysis of a clonal deletion in an epidemic strain of Mycobacterium bovis reveals a role in lipid metabolism. Microbiology (Reading). 2008;154(Pt 12):3731–42. https://doi.org/10.1099/mic.0.2008/022269-0.
Malone KM, Farrell D, Stuber TP, Schubert OT, Aebersold R, Robbe-Austerman S, et al. Updated Reference Genome Sequence and Annotation of Mycobacterium bovis AF2122/97. Genome Announc. 2017;5(14): https://doi.org/10.1128/genomeA.00157-17.
Bentley SD, Comas I, Bryant JM, Walker D, Smith NH, Harris SR, et al. The genome of Mycobacterium africanum West African 2 reveals a lineage-specific locus and genome erosion common to the M. tuberculosis complex. PLoS Negl Trop Dis. 2012;6(2):e1552. https://doi.org/10.1371/journal.pntd.0001552.
Mata E, Farrell D, Ma R, Uranga S, Gomez AB, Monzon M, et al. Independent genomic polymorphisms in the PknH serine threonine kinase locus during evolution of the Mycobacterium tuberculosis Complex affect virulence and host preference. PLoS Pathog. 2020;16(12):e1009061. https://doi.org/10.1371/journal.ppat.1009061.
Malone KM, Rue-Albrecht K, Magee DA, Conlon K, Schubert OT, Nalpas NC, et al. Comparative 'omics analyses differentiate Mycobacterium tuberculosis and Mycobacterium bovis and reveal distinct macrophage responses to infection with the human and bovine tubercle bacilli. Microb Genom. 2018;4(3): https://doi.org/10.1099/mgen.0.000163.
Farrell D, Crispell J, Gordon SV. Updated functional annotation of the Mycobacterium bovis AF2122/97 reference genome. Access Microbiol. 2020;2(7):acmi000129. https://doi.org/10.1099/acmi.0.000129.
Charlet D, Mostowy S, Alexander D, Sit L, Wiker HG, Behr MA. Reduced expression of antigenic proteins MPB70 and MPB83 in Mycobacterium bovis BCG strains due to a start codon mutation in sigK. Mol Microbiol. 2005;56(5):1302–13. https://doi.org/10.1111/j.1365-2958.2005.04618.x.
Said-Salim B, Mostowy S, Kristof AS, Behr MA. Mutations in Mycobacterium tuberculosis Rv0444c, the gene encoding anti-SigK, explain high level expression of MPB70 and MPB83 in Mycobacterium bovis. Mol Microbiol. 2006;62(5):1251–63. https://doi.org/10.1111/j.1365-2958.2006.05455.x.
Veyrier FJ, Nieves C, Lefrancois LH, Trigui H, Vincent AT, Behr MA. RskA Is a Dual Function Activator-Inhibitor That Controls SigK Activity Across Distinct Bacterial Genera. Front Microbiol. 2020;11:558166. https://doi.org/10.3389/fmicb.2020.558166.
Nagai S, Matsumoto J, Nagasuga T. Specific skin-reactive protein from culture filtrate of Mycobacterium bovis BCG. Infect Immun. 1981;31(3):1152–60. https://doi.org/10.1128/iai.31.3.1152-1160.1981.
Wiker HG, Harboe M, Nagai S. A localization index for distinction between extracellular and intracellular antigens of Mycobacterium tuberculosis. J Gen Microbiol. 1991;137(4):875–84. https://doi.org/10.1099/00221287-137-4-875.
Harboe M, Wiker HG, Ulvund G, Lund-Pedersen B, Andersen AB, Hewinson RG, et al. MPB70 and MPB83 as indicators of protein localization in mycobacterial cells. Infect Immun. 1998;66(1):289–96. https://doi.org/10.1128/IAI.66.1.289-296.1998.
Veyrier F, Said-Salim B, Behr MA. Evolution of the mycobacterial SigK regulon. J Bacteriol. 2008;190(6):1891–9. https://doi.org/10.1128/JB.01452-07.
Sharp JD, Singh AK, Park ST, Lyubetskaya A, Peterson MW, Gomes AL, et al. Comprehensive Definition of the SigH Regulon of Mycobacterium tuberculosis Reveals Transcriptional Control of Diverse Stress Responses. PLoS One. 2016;11(3):e0152145. https://doi.org/10.1371/journal.pone.0152145.
Sun R, Converse PJ, Ko C, Tyagi S, Morrison NE, Bishai WR. Mycobacterium tuberculosis ECF sigma factor sigC is required for lethality in mice and for the conditional expression of a defined gene set. Mol Microbiol. 2004;52(1):25–38. https://doi.org/10.1111/j.1365-2958.2003.03958.x.
Carr MD, Bloemink MJ, Dentten E, Whelan AO, Gordon SV, Kelly G, et al. Solution structure of the Mycobacterium tuberculosis complex protein MPB70: from tuberculosis pathogenesis to inherited human corneal desease. J Biol Chem. 2003;278(44):43736–43. https://doi.org/10.1074/jbc.M307235200.
Clout NJ, Tisi D, Hohenester E. Novel fold revealed by the structure of a FAS1 domain pair from the insect cell adhesion molecule fasciclin I. Structure. 2003;11(2):197–203. https://doi.org/10.1016/s0969-2126(03)00002-9.
Seifert GJ. Fascinating Fasciclins: A Surprisingly Widespread Family of Proteins that Mediate Interactions between the Cell Exterior and the Cell Surface. Int J Mol Sci. 2018;19(6): https://doi.org/10.3390/ijms19061628.
Park SY, Jung MY, Kim IS. Stabilin-2 mediates homophilic cell-cell interactions via its FAS1 domains. FEBS Lett. 2009;583(8):1375–80. https://doi.org/10.1016/j.febslet.2009.03.046.
Huber O, Sumper M. Algal-CAMs: isoforms of a cell adhesion molecule in embryos of the alga Volvox with homology to Drosophila fasciclin I. EMBO J. 1994;13(18):4212–22. https://doi.org/10.1002/j.1460-2075.1994.tb06741.x.
Moody RG, Williamson MP. Structure and function of a bacterial Fasciclin I Domain Protein elucidates function of related cell adhesion proteins such as TGFBIp and periostin. FEBS Open Bio. 2013;3:71–7. https://doi.org/10.1016/j.fob.2013.01.001.
Burroughs AM, Balaji S, Iyer LM, Aravind L. Small but versatile: the extraordinary functional and structural diversity of the beta-grasp fold. Biol Direct. 2007;2:18. https://doi.org/10.1186/1745-6150-2-18.
Politz O, Gratchev A, McCourt PA, Schledzewski K, Guillot P, Johansson S, et al. Stabilin-1 and -2 constitute a novel family of fasciclin-like hyaluronan receptor homologues. Biochem J. 2002;362(Pt 1):155–64. https://doi.org/10.1042/0264-6021:3620155.
Queval CJ, Fearns A, Botella L, Smyth A, Schnettger L, Mitermite M, et al. Macrophage-specific responses to human- and animal-adapted tubercle bacilli reveal pathogen and host factors driving multinucleated cell formation. PLoS Pathog. 2021;17(3):e1009410. https://doi.org/10.1371/journal.ppat.1009410.
Park SY, Yun Y, Lim JS, Kim MJ, Kim SY, Kim JE, et al. Stabilin-2 modulates the efficiency of myoblast fusion during myogenic differentiation and muscle regeneration. Nat Commun. 2016;7:10871. https://doi.org/10.1038/ncomms10871.
Park SY, Kim SY, Jung MY, Bae DJ, Kim IS. Epidermal growth factor-like domain repeat of stabilin-2 recognizes phosphatidylserine during cell corpse clearance. Mol Cell Biol. 2008;28(17):5288–98. https://doi.org/10.1128/MCB.01993-07.
Menon S, Wang S. Structure of the response regulator PhoP from Mycobacterium tuberculosis reveals a dimer through the receiver domain. Biochemistry. 2011;50(26):5948–57. https://doi.org/10.1021/bi2005575.
Solans L, Gonzalo-Asensio J, Sala C, Benjak A, Uplekar S, Rougemont J, et al. The PhoP-dependent ncRNA Mcr7 modulates the TAT secretion system in Mycobacterium tuberculosis. PLoS Pathog. 2014;10(5):e1004183. https://doi.org/10.1371/journal.ppat.1004183.
Groschel MI, Sayes F, Simeone R, Majlessi L, Brosch R. ESX secretion systems: mycobacterial evolution to counter host immunity. Nat Rev Microbiol. 2016;14(11):677–91. https://doi.org/10.1038/nrmicro.2016.131.
Perez E, Samper S, Bordas Y, Guilhot C, Gicquel B, Martin C. An essential role for phoP in Mycobacterium tuberculosis virulence. Mol Microbiol. 2001;41(1):179–87. https://doi.org/10.1046/j.1365-2958.2001.02500.x.
Gonzalo Asensio J, Maia C, Ferrer NL, Barilone N, Laval F, Soto CY, et al. The virulence-associated two-component PhoP-PhoR system controls the biosynthesis of polyketide-derived lipids in Mycobacterium tuberculosis. J Biol Chem. 2006;281(3):1313–6. https://doi.org/10.1074/jbc.C500388200.
Gonzalo-Asensio J, Malaga W, Pawlik A, Astarie-Dequeker C, Passemar C, Moreau F, et al. Evolutionary history of tuberculosis shaped by conserved mutations in the PhoPR virulence regulator. Proc Natl Acad Sci U S A. 2014;111(31):11491–6. https://doi.org/10.1073/pnas.1406693111.
Garcia EA, Blanco FC, Bigi MM, Vazquez CL, Forrellad MA, Rocha RV, et al. Characterization of the two component regulatory system PhoPR in Mycobacterium bovis. Vet Microbiol. 2018;222:30–8. https://doi.org/10.1016/j.vetmic.2018.06.016.
Garcia EA, Blanco FC, Klepp LI, Pazos A, McNeil MR, Jackson M, et al. Role of PhoPR in the response to stress of Mycobacterium bovis. Comp Immunol Microbiol Infect Dis. 2021;74: 101593. https://doi.org/10.1016/j.cimid.2020.101593.
We thank all colleagues whose research has informed and guided the ideas described in this manuscript.
About this supplement
This article has been published as part of Irish Veterinary Journal Volume 76 Supplement 1, 2023: Seventh International Conference on Mycobacterium bovis. The full contents of the supplement are available online at https://irishvetjournal.biomedcentral.com/articles/supplements/volume-76-supplement-1.
Work described in this manuscript was supported by the Wellcome Trust (PhD studentship 109166/Z/15/A to MM), Science Foundation Ireland (SFI/15/IA/3154 to SG and JMUE), Department of Agriculture Food and the Marine (DAFM) award 2019R404 (BTBGenIE to SG and DF), and the China Scholarship Council (CSC to RM). Funding bodies had no role in the conceptualization, design, data collection, analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
Consent for publication
None of the authors have any competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Mitermite, M., Elizari, J.M.U., Ma, R. et al. Exploring virulence in Mycobacterium bovis: clues from comparative genomics and perspectives for the future. Ir Vet J 76 (Suppl 1), 26 (2023). https://doi.org/10.1186/s13620-023-00257-6