Updated Reference Genome Sequence and Annotation of Mycobacterium bovis AF2122/97

The utility of whole-genome sequencing to identify likely transmission pairs for pathogens with slow and variable evolution

Pathogen whole-genome sequencing (WGS) has been used to track the transmission of infectious diseases in extraordinary detail, especially for pathogens that undergo fast and steady evolution, as is the case with many RNA viruses. However, for other pathogens evolution is less predictable, making interpretation of these data to inform our understanding of their epidemiology more challenging and the value of densely collected pathogen genome data uncertain. Here, we assess the utility of WGS for one such pathogen, in the "who-infected-whom" identification problem. We study samples from hosts (130 cattle, 111 badgers) with confirmed infection of M. bovis (causing bovine Tuberculosis), which has an estimated clock rate as slow as ∼0.1-1 variations per year. For each potential pathway between hosts, we calculate the relative likelihood that such a transmission event occurred. This is informed by an epidemiological model of transmission, and host life history data. By including WGS data, we shrink the number of plausible pathways significantly, relative to those deemed likely on the basis of life history data alone. Despite our uncertainty relating to the evolution of M. bovis, the WGS data are therefore a valuable adjunct to epidemiological investigations, especially for wildlife species whose life history data are sparse.

Whole-Genome sequencing in routine Mycobacterium bovis epidemiology - scoping the potential

Mycobacterium bovis the main agent of bovine tuberculosis (bTB), presents as a series of spatially-localised micro-epidemics across landscapes. Classical molecular typing methods applied to these micro-epidemics, based on genotyping a few variable loci, have significantly improved our understanding of potential epidemiological links between outbreaks. However, they have limited utility owing to low resolution. Conversely, whole-genome sequencing (WGS) provides the highest resolution data available for molecular epidemiology, producing richer outbreak tracing, insights into phylogeography and epidemic evolutionary history. We illustrate these advantages by focusing on a common single lineage of M. bovis (1.140) from Northern Ireland. Specifically, we investigate the spatial sub-structure of 20 years of herd-level multi locus VNTR analysis (MLVA) surveillance data and WGS data from a down sampled subset of isolates of this MLVA type over the same time frame. We mapped 2108 isolate locations of MLVA type 1.140 over the years 2000-2022. We also mapped the locations of 148 contemporary WGS isolates from this lineage, over a similar geographic range, stratifying by single nucleotide polymorphism (SNP) relatedness cut-offs of 15 SNPs. We determined a putative core range for the 1.140 MLVA type and SNP-defined sequence clusters using a 50 % kernel density estimate, using cattle movement data to inform on likely sources of WGS isolates found outside of core ranges. Finally, we applied Bayesian phylogenetic methods to investigate past population history and reproductive number of the 1.140 M. bovis lineage. We demonstrate that WGS SNP-defined clusters exhibit smaller core ranges than the established MLVA type - facilitating superior disease tracing. We also demonstrate the superior functionality of WGS data in determining how this lineage was disseminated across the landscape, likely via cattle movement and to infer how its effective population size and reproductive number has been in flux since its emergence. These initial findings highlight the potential of WGS data for routine monitoring of bTB outbreaks.

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.

Complete Genome Sequences of Mycobacterium bovis Strains Affiliated with Bovine Tuberculosis Outbreaks in Canada in 2016 and 2018

Mycobacterium bovis is the primary causative agent of bovine tuberculosis, a zoonotic infectious disease that presents a risk to public health, livestock, and wildlife. Here, we report complete genome sequences of two Mycobacterium bovis strains affiliated with bovine tuberculosis outbreaks in Canadian cattle farms in 2016 and 2018.

Use of Whole-Genome Sequencing to Unravel the Genetic Diversity of a Prevalent Mycobacterium bovis Spoligotype in a Multi-Host Scenario in Spain

Despite the efforts invested in the eradication of bovine tuberculosis in Spain, herd prevalence has remained constant in the country during the last 15 years (~1.5–1.9%) due to a combination of epidemiological factors impairing disease control, including between-species transmission. Here, our aim was to investigate the molecular diversity of Mycobacterium bovis isolates belonging to the highly prevalent SB0339 spoligotype in the cattle-wildlife interface in different regions of Spain using whole-genome sequencing (WGS). Genomic data of 136 M. bovis isolates recovered from different animal species (cattle, wild boar, fallow deer, and red deer) and locations between 2005 and 2018 were analyzed to investigate between- and within-species transmission, as well as within-herds. All sequenced isolates differed by 49–88 single nucleotide polymorphisms from their most recent common ancestor. Genetic heterogeneity was geographic rather than host species-specific, as isolates recovered from both cattle and wildlife from a given region were more closely related compared to isolates from the same species but geographically distant. In fact, a strong association between the geographic and the genetic distances separating pairs of M. bovis isolates was found, with a significantly stronger effect when cattle isolates were compared with wildlife or cattle-wildlife isolates in Spain. The same results were obtained in Madrid, the region with the largest number of sequenced isolates, but no differences depending on the host were observed. Within-herd genetic diversity was limited despite the considerable time elapsed between isolations. The detection of closely related strains in different hosts demonstrates the complex between-host transmission dynamics present in endemic areas in Spain. In conclusion, WGS results a valuable tool to track bTB infection at a high resolution and may contribute to achieve its eradication in Spain.

A critical evaluation of Mycobacterium bovis pangenomics, with reference to its utility in outbreak investigation

The increased accessibility of next generation sequencing has allowed enough genomes from a given bacterial species to be sequenced to describe the distribution of genes in the pangenome, without limiting analyses to genes present in reference strains. Although some taxa have thousands of whole genome sequences available on public databases, most genomes were sequenced with short read technology, resulting in incomplete assemblies. Studying pangenomes could lead to important insights into adaptation, pathogenicity, or molecular epidemiology, however given the known information loss inherent in analyzing contig-level assemblies, these inferences may be biased or inaccurate. In this study we describe the pangenome of a clonally evolving pathogen, Mycobacterium bovis , and examine the utility of gene content variation in M. bovis outbreak investigation. We constructed the M. bovis pangenome using 1463 de novo assembled genomes. We tested the assumption of strict clonal evolution by studying evidence of recombination in core genes and analyzing the distribution of accessory genes among core monophyletic groups. To determine if gene content variation could be utilized in outbreak investigation, we carefully examined accessory genes detected in a well described M. bovis outbreak in Minnesota. We found significant errors in accessory gene classification. After accounting for these errors, we show that M. bovis has a much smaller accessory genome than previously described and provide evidence supporting ongoing clonal evolution and a closed pangenome, with little gene content variation generated over outbreaks. We also identified frameshift mutations in multiple genes, including a mutation in glpK , which has recently been associated with antibiotic tolerance in Mycobacterium tuberculosis . A pangenomic approach enables a more comprehensive analysis of genome dynamics than is possible with reference-based approaches; however, without critical evaluation of accessory gene content, inferences of transmission patterns employing these loci could be misguided.

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.

Proteome characterization of the culture supernatant of Mycobacterium bovis in different growth stages

This study aimed to identify proteins secreted by Mycobacterium bovis into culture medium at different stages of bacterial growth. A field strain of M. bovis was grown in Middlebrook 7H9 media and culture supernatant was collected at three-time points representing three different phases of growth (early exponential, late exponential, and stationary phases). Supernatants were double filtered, digested by trypsin and analyzed by LC-MS/MS. The study found 15, 21, and 16 proteins in early, mid and late growth phases, respectively. In total, 22 proteins were identified, 18 of which were reported or predicted to have a cell wall or extracellular localization. To our knowledge, this is the first study to identify proteins secreted into the culture medium by a field strain of M. bovis in three different stages of growth. The dataset generated here provides candidate proteins with the potential for the development of serological diagnostic reagents or vaccine for bovine tuberculosis. Data are available via ProteomeXchange with identifier PXD017817.

•

Identifying extracellular proteins produced by M. bovis using mass spectrometry.

•

Culture supernatant proteins secreted in three different stages of bacterial growth.

•

Proteins were secreted from live bacteria and heat-killing was omitted.

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.

Mycobacterium tuberculosis complex lineage 5 exhibits high levels of within-lineage genomic diversity and differing gene content compared to the type strain H37Rv

Pathogens of the Mycobacterium tuberculosis complex (MTBC) are considered to be monomorphic, with little gene content variation between strains. Nevertheless, several genotypic and phenotypic factors separate strains of the different MTBC lineages (L), especially L5 and L6 (traditionally termed Mycobacterium africanum) strains, from each other. However, this genome variability and gene content, especially of L5 strains, has not been fully explored and may be important for pathobiology and current approaches for genomic analysis of MTBC strains, including transmission studies. By comparing the genomes of 355 L5 clinical strains (including 3 complete genomes and 352 Illumina whole-genome sequenced isolates) to each other and to H37Rv, we identified multiple genes that were differentially present or absent between H37Rv and L5 strains. Additionally, considerable gene content variability was found across L5 strains, including a split in the L5.3 sub-lineage into L5.3.1 and L5.3.2. These gene content differences had a small knock-on effect on transmission cluster estimation, with clustering rates influenced by the selected reference genome, and with potential overestimation of recent transmission when using H37Rv as the reference genome. We conclude that full capture of the gene diversity, especially high-resolution outbreak analysis, requires a variation of the single H37Rv-centric reference genome mapping approach currently used in most whole-genome sequencing data analysis pipelines. Moreover, the high within-lineage gene content variability suggests that the pan-genome of M. tuberculosis is at least several kilobases larger than previously thought, implying that a concatenated or reference-free genome assembly (de novo) approach may be needed for particular questions.

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.

ESX-1-Independent Horizontal Gene Transfer by Mycobacterium tuberculosis Complex Strains

Current models of horizontal gene transfer (HGT) in mycobacteria are based on “distributive conjugal transfer” (DCT), an HGT type described in the fast-growing, saprophytic model organism Mycobacterium smegmatis, which creates genome mosaicism in resulting strains and depends on an ESX-1 type VII secretion system. In contrast, only few data on interstrain DNA transfer are available for tuberculosis-causing mycobacteria, for which chromosomal DNA transfer between two Mycobacterium canettii strains was reported, a process which, however, was not observed for Mycobacterium tuberculosis strains. Here, we have studied a wide range of human- and animal-adapted members of the Mycobacterium tuberculosis complex (MTBC) using an optimized filter-based mating assay together with three selected strains of M. canettii that acted as DNA recipients. Unlike in previous approaches, we obtained a high yield of thousands of recombinants containing transferred chromosomal DNA fragments from various MTBC donor strains, as confirmed by whole-genome sequence analysis of 38 randomly selected clones. While the genome organizations of the obtained recombinants showed mosaicisms of donor DNA fragments randomly integrated into a recipient genome backbone, reminiscent of those described as being the result of ESX-1-mediated DCT in M. smegmatis, we observed similar transfer efficiencies when ESX-1-deficient donor and/or recipient mutants were used, arguing that in tubercle bacilli, HGT is an ESX-1-independent process. These findings provide new insights into the genetic events driving the pathoevolution of M. tuberculosis and radically change our perception of HGT in mycobacteria, particularly for those species that show recombinogenic population structures despite the natural absence of ESX-1 secretion systems.

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.

Identification of drug resistance mutations among Mycobacterium bovis lineages in the Americas

Identifying the Mycobacterium tuberculosis resistance mutation patterns is of the utmost importance to assure proper patient's management and devising of control programs aimed to limit spread of disease. Zoonotic Mycobacterium bovis infection still represents a threat to human health, particularly in dairy production regions. Routinary, molecular characterization of M. bovis is performed primarily by spoligotyping and mycobacterial interspersed repetitive units (MIRU) while next generation sequencing (NGS) approaches are often performed by reference laboratories. However, spoligotyping and MIRU methodologies lack the resolution required for the fine characterization of tuberculosis isolates, particularly in outbreak settings. In conjunction with sophisticated bioinformatic algorithms, whole genome sequencing (WGS) analysis is becoming the method of choice for advanced genetic characterization of tuberculosis isolates. WGS provides valuable information on drug resistance and compensatory mutations that other technologies cannot assess. Here, we performed an analysis of the most frequently identified mutations associated with tuberculosis drug resistance and their genetic relationship among 2,074 Mycobacterium bovis WGS recovered primarily from non-human hosts. Full-length gene sequences harboring drug resistant associated mutations and their phylogenetic relationships were analyzed. The results showed that M. bovis isolates harbor mutations conferring resistance to both first- and second-line antibiotics. Mutations conferring resistance for isoniazid, fluoroquinolones, streptomycin, and aminoglycosides were identified among animal strains. Our findings highlight the importance of molecular surveillance to monitor the emergence of mutations associated with multi and extensive drug resistance in livestock and other non-human mammals.

Pathogenomic analyses of Mycobacterium microti, an ESX-1-deleted member of the Mycobacterium tuberculosis complex causing disease in various hosts

Mycobacterium microti is an animal-adapted member of the Mycobacterium tuberculosis complex (MTBC), which was originally isolated from voles, but has more recently also been isolated from other selected mammalian hosts, including occasionally from humans. Here, we have generated and analysed the complete genome sequences of five representative vole and clinical M. microti isolates using PacBio- and Illumina-based technologies, and have tested their virulence and vaccine potential in SCID (severe combined immune deficient) mouse and/or guinea pig infection models. We show that the clinical isolates studied here cluster separately in the phylogenetic tree from vole isolates and other clades from publicly available M. microti genome sequences. These data also confirm that the vole and clinical M. microti isolates were all lacking the specific RD1mic region, which in other tubercle bacilli encodes the ESX-1 type VII secretion system. Biochemical analysis further revealed marked phenotypic differences between isolates in type VII-mediated secretion of selected PE and PPE proteins, which in part were attributed to specific genetic polymorphisms. Infection experiments in the highly susceptible SCID mouse model showed that the clinical isolates were significantly more virulent than the tested vole isolates, but still much less virulent than the M. tuberculosis H37Rv control strain. The strong attenuation of the ATCC 35872 vole isolate in immunocompromised mice, even compared to the attenuated BCG (bacillus Calmette-Guérin) vaccine, and its historic use in human vaccine trials encouraged us to test this strain's vaccine potential in a guinea pig model, where it demonstrated similar protective efficacy as a BCG control, making it a strong candidate for vaccination of immunocompromised individuals in whom BCG vaccination is contra-indicated. Overall, we provide new insights into the genomic and phenotypic variabilities and particularities of members of an understudied clade of the MTBC, which all share a recent common ancestor that is characterized by the deletion of the RD1mic region.

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.

Phylogenetic relationships investigation of Mycobacterium caprae strains from sympatric wild boar and goats based on whole genome sequencing

Tuberculosis (TB) in wildlife challenges epidemiological surveillance and disease control. An outbreak of TB was detected in a free‐ranging wild boar population of a Natural Park in Catalonia (Spain) and the outbreak investigation was conducted in the area. During the study period (2015–2020), 278 wild boars were analysed by gross pathology, histopathology, mycobacterial culture and DVR‐spoligotyping. In addition, all cattle (49) and goat (47) herds of the area were tested with tuberculin skin test. TB compatible lesions were detected in 21 wild boars, and Mycobacterium caprae was isolated in 17 of them with two different spoligotypes: SB0415 (13) and SB1908 (4). Only two goat herds showed TB positive animals that were subsequently slaughtered. M. caprae with the spoligotypes SB0416 and SB0415 were isolated from these animals. To investigate the phylogenetic relationships and the transmission chain of the outbreak, nine strains isolated from six wild boars and three goats of the study area were analysed by whole genome sequencing (WGS) followed by single nucleotide polymorphism (SNP) analysis by maximum likelihood and median‐joining network inference methods. Results indicated that infected wild boars maintained M. caprae strains circulation in their own population and have likely transmitted the infection to goats, thus acting as TB reservoirs, compromising the success of livestock TB eradication campaigns and posing a risk for public health. The results also highlighted the usefulness of WGS followed by SNP analysis in providing relevant epidemiological information when detailed contact data are missing.

Mycobacterium bovis genomics reveals transmission of infection between cattle and deer in Ireland

Control of bovine tuberculosis (bTB), caused by Mycobacterium bovis, in the Republic of Ireland costs €84 million each year. Badgers are recognized as being a wildlife source for M. bovis infection of cattle. Deer are thought to act as spillover hosts for infection; however, population density is recognized as an important driver in shifting their epidemiological role, and deer populations across the country have been increasing in density and range. County Wicklow represents one specific area in the Republic of Ireland with a high density of deer that has had consistently high bTB prevalence for over a decade, despite control operations in both cattle and badgers. Our research used whole-genome sequencing of M. bovis sourced from infected cattle, deer and badgers in County Wicklow to evaluate whether the epidemiological role of deer could have shifted from spillover host to source. Our analyses reveal that cattle and deer share highly similar M. bovis strains, suggesting that transmission between these species is occurring in the area. In addition, the high level of diversity observed in the sampled deer population suggests deer may be acting as a source of infection for local cattle populations. These findings have important implications for the control and ultimate eradication of bTB in Ireland.

Mycobacterium bovis: From Genotyping to Genome Sequencing

Mycobacterium bovis is the main pathogen of bovine, zoonotic, and wildlife tuberculosis. Despite the existence of programs for bovine tuberculosis (bTB) control in many regions, the disease remains a challenge for the veterinary and public health sectors, especially in developing countries and in high-income nations with wildlife reservoirs. Current bTB control programs are mostly based on test-and-slaughter, movement restrictions, and post-mortem inspection measures. In certain settings, contact tracing and surveillance has benefited from M. bovis genotyping techniques. More recently, whole-genome sequencing (WGS) has become the preferential technique to inform outbreak response through contact tracing and source identification for many infectious diseases. As the cost per genome decreases, the application of WGS to bTB control programs is inevitable moving forward. However, there are technical challenges in data analyses and interpretation that hinder the implementation of M. bovis WGS as a molecular epidemiology tool. Therefore, the aim of this review is to describe M. bovis genotyping techniques and discuss current standards and challenges of the use of M. bovis WGS for transmission investigation, surveillance, and global lineages distribution. We compiled a series of associated research gaps to be explored with the ultimate goal of implementing M. bovis WGS in a standardized manner in bTB control programs.

Updated functional annotation of the Mycobacterium bovis AF2122/97 reference genome

Mycobacterium bovis AF2122/97 is the reference strain for the bovine tuberculosis bacillus. Here we report an update to the M. bovis AF2122/97 genome annotation to reflect 616 new protein identifications that replace many of the old hypothetical coding sequences and proteins of unknown function in the genome. These changes integrate information from functional assignments of orthologous coding sequences in the Mycobacterium tuberculosis H37Rv genome. We have also added 69 additional new gene names.

The complete genome sequence of Mycobacterium bovis Mb3601, a SB0120 spoligotype strain representative of a new clonal group

Mycobacterium bovis strain Mb3601 was isolated from the lymph node of an infected bovine in a bovine tuberculosis highly enzoonotic area of Burgundy, France. It was selected to obtain a complete genome for a new clonal complex, mainly constituted by SB0120-spoligotype strains that we propose to name "European 3". It was recently described as "clonal group I" based on whole-genome SNP analysis of 87 French strains. Here we describe the 4,365,068 bp complete genome obtained by the combination of PacBio and Illumina technologies. This genome of 65.64% G + C content includes 4024 predicted protein-coding genes, 52 tRNA, 3 rRNA and 11 copies of IS6110.

Phylogenomic Perspective on a Unique Mycobacterium bovis Clade Dominating Bovine Tuberculosis Infections among Cattle and Buffalos in Northern Brazil

Lack of routine surveillance in countries endemic for bovine tuberculosis (TB) and limited laboratory support contributes to the inability to differentiate the Mycobacterium tuberculosis Complex species, leading to an underestimated burden of the disease. Here, Whole-Genome Sequencing of Mycobacterium bovis isolated from tissues with TB-like lesions obtained from cattle and buffalos at Marajó Island, Brazil, demonstrates that recent transmission of M. bovis is ongoing at distinct sites. Moreover, the M. bovis epidemiology in this setting is herein found to be dominated by an endemic and unique clade composed of strains evolved from a common ancestor that are now genetically differentiated from other M. bovis clades. Additionally, envisioning a rapid strain differentiation and tracing across multiple settings, 28 globally validated strain-specific SNPs were identified, three of which considered as robust markers for the M. bovis Marajó strain. In conclusion, this study contributes with data regarding the identification of a novel M. bovis phylogenetic clade responsible for ongoing transmission events in both cattle and buffalo species in Brazil, provides a framework to investigate the dissemination of this highly prevalent strain and, holds the potential to inform TB control strategies that may help to prevent the spread of bovine and zoonotic TB.

Combining genomics and epidemiology to analyse bi-directional transmission of Mycobacterium bovis in a multi-host system

Quantifying pathogen transmission in multi-host systems is difficult, as exemplified in bovine tuberculosis (bTB) systems, but is crucial for control. The agent of bTB, Mycobacterium bovis, persists in cattle populations worldwide, often where potential wildlife reservoirs exist. However, the relative contribution of different host species to bTB persistence is generally unknown. In Britain, the role of badgers in infection persistence in cattle is highly contentious, despite decades of research and control efforts. We applied Bayesian phylogenetic and machine-learning approaches to bacterial genome data to quantify the roles of badgers and cattle in M. bovis infection dynamics in the presence of data biases. Our results suggest that transmission occurs more frequently from badgers to cattle than vice versa (10.4x in the most likely model) and that within-species transmission occurs at higher rates than between-species transmission for both. If representative, our results suggest that control operations should target both cattle and badgers.

Genomic Polymorphism Associated with the Emergence of Virulent Isolates of Mycobacterium bovis in the Nile Delta

Mycobacterium bovis is responsible for bovine tuberculosis in both animals and humans. Despite being one of the most important global zoonotic disease, data related to the ecology and pathogenicity of bovine tuberculosis is scarce, especially in developing countries. In this report, we examined the dynamics of M. bovis transmission among dairy cattle in the Nile Delta of Egypt. Animals belonging to 27 herds from 7 governorates were tested by the Single Intradermal Comparative Skin Tuberculin (SICST), as a preliminary screen for the presence of bovine tuberculosis. Positive SICST reactors were identified in 3% of the animals spread among 40% of the examined herds. Post-mortem examination of slaughtered reactors confirmed the presence of both pulmonary and/or digestive forms of tuberculosis in > 50% of the examined animals. Targeted and whole-genome analysis of M. bovis isolates indicated the emergences of a predominant spoligotype (SB0268) between 2013–2015, suggesting a recent clonal spread of this isolate within the Nile Delta. Surprisingly, 2 isolates belonged to M. bovis BCG group, which are not allowed for animal vaccination in Egypt, while the rest of isolates belonged to the virulent M. bovis clonal complex European 2 present in Latin America and several European countries. Analysis of strain virulence in the murine model of tuberculosis indicated the emergence of a more virulent strain (MBE4) with a specific genotype. More analysis is needed to understand the molecular basis for successful spread of virulent isolates of bovine tuberculosis among animals and to establish genotype/phenotype association.

HomoplasyFinder: a simple tool to identify homoplasies on a phylogeny

A homoplasy is a nucleotide identity resulting from a process other than inheritance from a common ancestor. Importantly, by distorting the ancestral relationships between nucleotide sequences, homoplasies can change the structure of the phylogeny. Homoplasies can emerge naturally, especially under high selection pressures and/or high mutation rates, or be created during the generation and processing of sequencing data. Identification of homoplasies is critical, both to understand their influence on the analyses of phylogenetic data and to allow an investigation into how they arose. Here we present HomoplasyFinder, a java application that can be used as a stand-a-lone tool or within the statistical programming environment R. Within R and Java, HomoplasyFinder is shown to be able to automatically, and quickly, identify any homoplasies present in simulated and real phylogenetic data. HomoplasyFinder can easily be incorporated into existing analysis pipelines, either within or outside of R, allowing the user to quickly identify homoplasies to inform downstream analyses and interpretation.

A New Phylogenetic Framework for the Animal-Adapted Mycobacterium tuberculosis Complex

Tuberculosis (TB) affects humans and other animals and is caused by bacteria from the Mycobacterium tuberculosis complex (MTBC). Previous studies have shown that there are at least nine members of the MTBC infecting animals other than humans; these have also been referred to as ecotypes. However, the ecology and the evolution of these animal-adapted MTBC ecotypes are poorly understood. Here we screened 12,886 publicly available MTBC genomes and newly sequenced 17 animal-adapted MTBC strains, gathering a total of 529 genomes of animal-adapted MTBC strains. Phylogenomic and comparative analyses confirm that the animal-adapted MTBC members are paraphyletic with some members more closely related to the human-adapted Mycobacterium africanum Lineage 6 than to other animal-adapted strains. Furthermore, we identified four main animal-adapted MTBC clades that might correspond to four main host shifts; two of these clades are hypothesized to reflect independent cattle domestication events. Contrary to what would be expected from an obligate pathogen, MTBC nucleotide diversity was not positively correlated with host phylogenetic distances, suggesting that host tropism in the animal-adapted MTBC seems to be driven by contact rates and demographic aspects of the host population rather by than host relatedness. By combining phylogenomics with ecological data, we propose an evolutionary scenario in which the ancestor of Lineage 6 and all animal-adapted MTBC ecotypes was a generalist pathogen that subsequently adapted to different host species. This study provides a new phylogenetic framework to better understand the evolution of the different ecotypes of the MTBC and guide future work aimed at elucidating the molecular mechanisms underlying host range.

Whole Genome Sequencing for Determining the Source of Mycobacterium bovis Infections in Livestock Herds and Wildlife in New Zealand

The ability to DNA fingerprint Mycobacterium bovis isolates helped to define the role of wildlife in the persistence of bovine tuberculosis in New Zealand. DNA fingerprinting results currently help to guide wildlife control measures and also aid in tracing the source of infections that result from movement of livestock. During the last 5 years we have developed the ability to distinguish New Zealand (NZ) M. bovis isolates by comparing the sequences of whole genome sequenced (WGS) M. bovis samples. WGS provides much higher resolution than our other established typing methods and greatly improves the definition of the regional localization of NZ M. bovis types. Three outbreak investigations are described and results demonstrate how WGS analysis has led to the confirmation of epidemiological sourcing of infection, to better definition of new sources of infection by ruling out other possible sources, and has revealed probable wildlife infection in an area considered to be free of infected wildlife. The routine use of WGS analyses for sourcing new M. bovis infections will be an important component of the strategy employed to eradicate bovine TB from NZ livestock and wildlife.

Whole Genome Sequencing of Mycobacterium bovis Isolated From Livestock in the United States, 1989-2018

The United States official bovine tuberculosis (bTB) eradication program has utilized genotyping for Mycobacterium bovis isolates since 2000 and whole genome sequencing was implemented in 2013. The program has been highly successful, yet as bTB prevalence has reached historic lows, a small number of new bTB-affected cattle herds occur annually. Therefore, understanding the epidemiology of bTB transmission is critically important, in order to target limited resources for surveillance and achieve eradication. This evaluation described the diversity and epidemiology of M. bovis isolates identified in the USA livestock. Isolates from animals within the bTB endemic area of Michigan were excluded. Broad diversity was found among 1,248 isolates, collected from affected cattle and farmed cervids herds and fed cattle during 1989–2018. Nearly 70% of isolates from 109 herds/cases during 1999–2018 were European clonal complex 1 and 30% were European clonal complex 2. The sources of infection based on the herd investigation were known for 41% of herds/cases and 59% were not epidemiologically linked to another USA origin herd. Whole genome sequencing results were consistent with the investigation findings and previously unrecognized links between herds and cases were disclosed. For herds/cases with an unknown source of infection, WGS results suggested several possible sources, including undocumented cattle movement, imported cattle and humans. The use of WGS in new cases has reduced the time and costs associated with epidemiological investigations. Within herd SNP diversity was evaluated by examining 18 herds with 10 or more isolates sequenced. Forty percent of isolates had not diverged or accumulated any SNPs, and 86% of the isolates had accumulated 3 or fewer SNPs. The results of WGS does not support a bTB reservoir in USA cattle. The bTB eradication program appears to be highly effective as the majority of herds/cases in the USA are unique strains with limited herd to herd transmission.

Genome Sequences of Five Mycobacterium bovis Strains Isolated from Farmed Animals and Wildlife in Canada

Mycobacterium bovis is the causative agent of bovine tuberculosis, an infectious disease that affects both animals and humans and thus presents a risk to public health and the livestock industry. Here, we report the genome sequences of five Mycobacterium bovis strains that represent major genotype clusters observed in farmed animals and wildlife in Canada.

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

Members of the Mycobacterium tuberculosis complex (MTBC) are the causative agents of tuberculosis in a range of mammals, including humans. A key feature of MTBC pathogens is their high degree of genetic identity yet distinct host tropism. Notably, while Mycobacterium bovis is highly virulent and pathogenic for cattle, the human pathogen M. tuberculosis is attenuated in cattle. Previous research also suggests that host preference amongst MTBC members has a basis in host innate immune responses. To explore MTBC host tropism, we present in-depth profiling of the MTBC reference strains M. bovis AF2122/97 and M. tuberculosis H37Rv at both the global transcriptional and the translational level via RNA-sequencing and SWATH MS. Furthermore, a bovine alveolar macrophage infection time course model was used to investigate the shared and divergent host transcriptomic response to infection with M. tuberculosis H37Rv or M. bovis AF2122/97. Significant differential expression of virulence-associated pathways between the two bacilli was revealed, including the ESX-1 secretion system. A divergent transcriptional response was observed between M. tuberculosis H37Rv and M. bovis AF2122/97 infection of bovine alveolar macrophages, in particular cytosolic DNA-sensing pathways at 48 h post-infection, and highlights a distinct engagement of M. bovis with the bovine innate immune system. The work presented here therefore provides a basis for the identification of host innate immune mechanisms subverted by virulent host-adapted mycobacteria to promote their survival during the early stages of infection.

Experimental infection of cattle with Mycobacterium tuberculosis isolates shows the attenuation of the human tubercle bacillus for cattle

The Mycobacterium tuberculosis complex (MTBC) is the collective term given to the group of bacteria that cause tuberculosis (TB) in mammals. It has been reported that M. tuberculosis H37Rv, a standard reference MTBC strain, is attenuated in cattle compared to Mycobacterium bovis. However, as M. tuberculosis H37Rv was isolated in the early 1930s, and genetic variants are known to exist, we sought to revisit this question of attenuation of M. tuberculosis for cattle by performing a bovine experimental infection with a recent M. tuberculosis isolate. Here we report infection of cattle using M. bovis AF2122/97, M. tuberculosis H37Rv, and M. tuberculosis BTB1558, the latter isolated in 2008 during a TB surveillance project in Ethiopian cattle. We show that both M. tuberculosis strains caused reduced gross pathology and histopathology in cattle compared to M. bovis. Using M. tuberculosis H37Rv and M. bovis AF2122/97 as the extremes in terms of infection outcome, we used RNA-Seq analysis to explore differences in the peripheral response to infection as a route to identify biomarkers of progressive disease in contrast to a more quiescent, latent infection. Our work shows the attenuation of M. tuberculosis strains for cattle, and emphasizes the potential of the bovine model as a ‘One Health’ approach to inform human TB biomarker development and post-exposure vaccine development.

Genomic repeats, misassembly and reannotation: a case study with long-read resequencing of Porphyromonas gingivalis reference strains

BACKGROUND

Without knowledge of their genomic sequences, it is impossible to make functional models of the bacteria that make up human and animal microbiota. Unfortunately, the vast majority of publicly available genomes are only working drafts, an incompleteness that causes numerous problems and constitutes a major obstacle to genotypic and phenotypic interpretation. In this work, we began with an example from the class Bacteroidia in the phylum Bacteroidetes, which is preponderant among human orodigestive microbiota. We successfully identify the genetic loci responsible for assembly breaks and misassemblies and demonstrate the importance and usefulness of long-read sequencing and curated reannotation.

RESULTS

We showed that the fragmentation in Bacteroidia draft genomes assembled from massively parallel sequencing linearly correlates with genomic repeats of the same or greater size than the reads. We also demonstrated that some of these repeats, especially the long ones, correspond to misassembled loci in three reference Porphyromonas gingivalis genomes marked as circularized (thus complete or finished). We prove that even at modest coverage (30X), long-read resequencing together with PCR contiguity verification (rrn operons and an integrative and conjugative element or ICE) can be used to identify and correct the wrongly combined or assembled regions. Finally, although time-consuming and labor-intensive, consistent manual biocuration of three P. gingivalis strains allowed us to compare and correct the existing genomic annotations, resulting in a more accurate interpretation of the genomic differences among these strains.

CONCLUSIONS

In this study, we demonstrate the usefulness and importance of long-read sequencing in verifying published genomes (even when complete) and generating assemblies for new bacterial strains/species with high genomic plasticity. We also show that when combined with biological validation processes and diligent biocurated annotation, this strategy helps reduce the propagation of errors in shared databases, thus limiting false conclusions based on incomplete or misleading information.

Whole genome sequencing of the monomorphic pathogen Mycobacterium bovis reveals local differentiation of cattle clinical isolates

Background

Bovine tuberculosis (bTB) poses serious risks to animal welfare and economy, as well as to public health as a zoonosis. Its etiological agent, Mycobacterium bovis, belongs to the Mycobacterium tuberculosis complex (MTBC), a group of genetically monomorphic organisms featured by a remarkably high overall nucleotide identity (99.9%). Indeed, this characteristic is of major concern for correct typing and determination of strain-specific traits based on sequence diversity. Due to its historical economic dependence on cattle production, Uruguay is deeply affected by the prevailing incidence of Mycobacterium bovis. With the world’s highest number of cattle per human, and its intensive cattle production, Uruguay represents a particularly suited setting to evaluate genomic variability among isolates, and the diversity traits associated to this pathogen.

Results

We compared 186 genomes from MTBC strains isolated worldwide, and found a highly structured population in M. bovis. The analysis of 23 new M. bovis genomes, belonging to strains isolated in Uruguay evidenced three groups present in the country. Despite presenting an expected highly conserved genomic structure and sequence, these strains segregate into a clustered manner within the worldwide phylogeny. Analysis of the non-pe/ppe differential areas against a reference genome defined four main sources of variability, namely: regions of difference (RD), variable genes, duplications and novel genes. RDs and variant analysis segregated the strains into clusters that are concordant with their spoligotype identities. Due to its high homoplasy rate, spoligotyping failed to reflect the true genomic diversity among worldwide representative strains, however, it remains a good indicator for closely related populations.

Conclusions

This study introduces a comprehensive population structure analysis of worldwide M. bovis isolates. The incorporation and analysis of 23 novel Uruguayan M. bovis genomes, sheds light onto the genomic diversity of this pathogen, evidencing the existence of greater genetic variability among strains than previously contemplated.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4249-6) contains supplementary material, which is available to authorized users.

The C3HeB/FeJ mouse model recapitulates the hallmark of bovine tuberculosis lung lesions following Mycobacterium bovis aerogenous infection

Achieving the control of bovine tuberculosis (bTB) would require the discovery of an efficient combined immunodiagnostic and vaccine strategy. Since in vivo experiments on cattle are not ethically and economically acceptable there is a need for a cost-effective animal model capable of reproducing, as closely as possible, the physiopathology of bTB to (i) better characterize the cellular and molecular features of bTB immunopathogenesis and (ii) screen preclinical vaccine candidates. To develop such a model, we focused on the C3HeB/FeJ Kramnik’s mouse forming hypoxic, encapsulated granulomas with a caseous necrotic center following Mycobacterium tuberculosis infection. Our work represents the first investigation on C3HeB/FeJ interaction with M. bovis, the main agent of bTB. Detailed histopathological analysis of C3HeB/FeJ lung lesions development following aerogenous M. bovis infection unraveled a bimodal evolution of the pathology. The C3HeB/FeJ recapitulated all the hallmarks of classical bovine lung granulomas but also developed, to some extend, lethal necrotic large lesions characterized by high mycobacterial and neutrophil load, and an inefficient collagen-driven lesion encapsulation. Interestingly these rapidly invasive pneumonia lesions, occurring in a constant percentage of the mice, shared all features with some exacerbated lung lesions that we and others have observed in lungs of cattle naturally or experimentally infected with M. bovis. Together, our findings demonstrate the relevance of the C3HeB/FeJ mouse as a comprehensive model to study bTB immunopathology that could be used for further vaccine therapies in the future.