Comparative Genomics of Field Isolates of Mycobacterium bovis and M. caprae Provides Evidence for Possible Correlates with Bacterial Viability and Virulence

Mycobacteria of the Mycobacterium tuberculosis complex (MTBC) greatly affect humans and animals worldwide. The life cycle of mycobacteria is complex and the mechanisms resulting in pathogen infection and survival in host cells are not fully understood. Recently, comparative genomics analyses have provided new insights into the evolution and adaptation of the MTBC to survive inside the host. However, most of this information has been obtained using M. tuberculosis but not other members of the MTBC such as M. bovis and M. caprae. In this study, the genome of three M. bovis (MB1, MB3, MB4) and one M. caprae (MB2) field isolates with different lesion score, prevalence and host distribution phenotypes were sequenced. Genome sequence information was used for whole-genome and protein-targeted comparative genomics analysis with the aim of finding correlates with phenotypic variation with potential implications for tuberculosis (TB) disease risk assessment and control. At the whole-genome level the results of the first comparative genomics study of field isolates of M. bovis including M. caprae showed that as previously reported for M. tuberculosis, sequential chromosomal nucleotide substitutions were the main driver of the M. bovis genome evolution. The phylogenetic analysis provided a strong support for the M. bovis/M. caprae clade, but supported M. caprae as a separate species. The comparison of the MB1 and MB4 isolates revealed differences in genome sequence, including gene families that are important for bacterial infection and transmission, thus highlighting differences with functional implications between isolates otherwise classified with the same spoligotype. Strategic protein-targeted analysis using the ESX or type VII secretion system, proteins linking stress response with lipid metabolism, host T cell epitopes of mycobacteria, antigens and peptidoglycan assembly protein identified new genetic markers and candidate vaccine antigens that warrant further study to develop tools to evaluate risks for TB disease caused by M. bovis/M.caprae and for TB control in humans and animals.

Complete Genome Sequences of Field Isolates of Mycobacterium bovis and Mycobacterium caprae

Here we report the complete genome sequences of field isolates of Mycobacterium bovis and the related mycobacterial species, Mycobacterium caprae. The genomes of three M. bovis (MB1, MB3, MB4) and one M. caprae (MB2) field isolates with different virulence, prevalence, and host distribution phenotypes were sequenced.

Tonsils of the soft palate do not mediate the response of pigs to oral vaccination with heat-inactivated Mycobacterium bovis

Mycobacterium bovis causes animal tuberculosis (TB) in cattle, humans, and other mammalian species, including pigs. The goal of this study was to experimentally assess the responses of pigs with and without a history of tonsillectomy to oral vaccination with heat-inactivated M. bovis and challenge with a virulent M. bovis field strain, to compare pig and wild boar responses using the same vaccination model as previously used in the Eurasian wild boar (Sus scrofa), to evaluate the use of several enzyme-linked immunosorbent assays (ELISAs) and lateral flow tests for in vivo TB diagnosis in pigs, and to verify if these tests are influenced by oral vaccination with inactivated M. bovis. At necropsy, the lesion and culture scores were 20% to 43% higher in the controls than those in the vaccinated pigs. Massive M. bovis growth from thoracic tissue samples was observed in 4 out of 9 controls but in none of the 10 vaccinated pigs. No effect of the presence or absence of tonsils was observed on these scores, suggesting that tonsils are not involved in the protective response to this vaccine in pigs. The serum antibody levels increased significantly only after challenge. At necropsy, the estimated sensitivities of the ELISAs and dual path platform (DPP) assays ranged from 89% to 94%. In the oral mucosa, no differences in gene expression were observed in the control group between the pigs with and without tonsils. In the vaccinated group, the mRNA levels for chemokine (C-C motif) receptor 7 (CCR7), interferon beta (IFN-β), and methylmalonyl coenzyme A mutase (MUT) were higher in pigs with tonsils. Complement component 3 mRNA levels in peripheral blood mononuclear cells (PBMC) increased with vaccination and decreased after M. bovis challenge. This information is relevant for pig production in regions that are endemic for M. bovis and for TB vaccine research.

Oral vaccination with heat inactivated Mycobacterium bovis activates the complement system to protect against tuberculosis

Tuberculosis (TB) remains a pandemic affecting billions of people worldwide, thus stressing the need for new vaccines. Defining the correlates of vaccine protection is essential to achieve this goal. In this study, we used the wild boar model for mycobacterial infection and TB to characterize the protective mechanisms elicited by a new heat inactivated Mycobacterium bovis vaccine (IV). Oral vaccination with the IV resulted in significantly lower culture and lesion scores, particularly in the thorax, suggesting that the IV might provide a novel vaccine for TB control with special impact on the prevention of pulmonary disease, which is one of the limitations of current vaccines. Oral vaccination with the IV induced an adaptive antibody response and activation of the innate immune response including the complement component C3 and inflammasome. Mycobacterial DNA/RNA was not involved in inflammasome activation but increased C3 production by a still unknown mechanism. The results also suggested a protective mechanism mediated by the activation of IFN-γ producing CD8+ T cells by MHC I antigen presenting dendritic cells (DCs) in response to vaccination with the IV, without a clear role for Th1 CD4+ T cells. These results support a role for DCs in triggering the immune response to the IV through a mechanism similar to the phagocyte response to PAMPs with a central role for C3 in protection against mycobacterial infection. Higher C3 levels may allow increased opsonophagocytosis and effective bacterial clearance, while interfering with CR3-mediated opsonic and nonopsonic phagocytosis of mycobacteria, a process that could be enhanced by specific antibodies against mycobacterial proteins induced by vaccination with the IV. These results suggest that the IV acts through novel mechanisms to protect against TB in wild boar.

Oral re-vaccination of Eurasian wild boar with Mycobacterium bovis BCG yields a strong protective response against challenge with a field strain

BACKGROUND Field vaccination trials with Mycobacterium bovis BCG, an attenuated mutant of M. bovis, are ongoing in Spain, where the Eurasian wild boar (Sus scrofa) is regarded as the main driver of animal tuberculosis (TB). The oral baiting strategy consists in deploying vaccine baits twice each summer, in order to gain access to a high proportion of wild boar piglets. The aim of this study was to assess the response of wild boar to re-vaccination with BCG and to subsequent challenge with an M. bovis field strain. RESULTS BCG re-vaccinated wild boar showed reductions of 75.8% in lesion score and 66.9% in culture score, as compared to unvaccinated controls. Only one of nine vaccinated wild boar had a culture-confirmed lung infection, as compared to seven of eight controls. Serum antibody levels were highly variable and did not differ significantly between BCG re-vaccinated wild boar and controls. Gamma IFN levels differed significantly between BCG re-vaccinated wild boar and controls. The mRNA levels for IL-1b, C3 and MUT were significantly higher in vaccinated wild boar when compared to controls after vaccination and decreased after mycobacterial challenge. CONCLUSIONS Oral re-vaccination of wild boar with BCG yields a strong protective response against challenge with a field strain. Moreover, re-vaccination of wild boar with BCG is not counterproductive. These findings are relevant given that re-vaccination is likely to happen under real (field) conditions.

Assessment of an oral Mycobacterium bovis BCG vaccine and an inactivated M. bovis preparation for wild boar in terms of adverse reactions, vaccine strain survival, and uptake by nontarget species

Wildlife vaccination is increasingly being considered as an option for tuberculosis control. We combined data from laboratory trials and an ongoing field trial to assess the risk of an oral Mycobacterium bovis BCG vaccine and a prototype heat-inactivated Mycobacterium bovis preparation for Eurasian wild boar (Sus scrofa). We studied adverse reactions, BCG survival, BCG excretion, and bait uptake by nontarget species. No adverse reactions were observed after administration of BCG (n = 27) or inactivated M. bovis (n = 21). BCG was not found at necropsy (175 to 300 days postvaccination [n = 27]). No BCG excretion was detected in fecal samples (n = 162) or in urine or nasal, oral, or fecal swab samples at 258 days postvaccination (n = 29). In the field, we found no evidence of loss of BCG viability in baits collected after 36 h (temperature range, 11°C to 41°C). Camera trapping showed that wild boar (39%) and birds (56%) were the most frequent visitors to bait stations (selective feeders). Wild boar activity patterns were nocturnal, while diurnal activities were recorded for all bird species. We found large proportions of chewed capsules (29%) (likely ingestion of the vaccine) and lost baits (39%) (presumably consumed), and the proportion of chewed capsules showed a positive correlation with the presence of wild boar. Both results suggest proper bait consumption (68%). These results indicate that BCG vaccination in wild boar is safe and that, while bait consumption by other species is possible, this can be minimized by using selective cages and strict timing of bait deployment.

Splitting of a prevalent Mycobacterium bovis spoligotype by variable-number tandem-repeat typing reveals high heterogeneity in an evolving clonal group

Mycobacterium bovis populations in countries with persistent bovine tuberculosis usually show a prevalent spoligotype with a wide geographical distribution. This study applied mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing to a random panel of 115 M. bovis isolates that are representative of the most frequent spoligotype in the Iberian Peninsula, SB0121. VNTR typing targeted nine loci: ETR-A (alias VNTR2165), ETR-B (VNTR2461), ETR-D (MIRU4, VNTR580), ETR-E (MIRU31, VNTR3192), MIRU26 (VNTR2996), QUB11a (VNTR2163a), QUB11b (VNTR2163b), QUB26 (VNTR4052), and QUB3232 (VNTR3232). We found a high degree of diversity among the studied isolates (discriminatory index [D] = 0.9856), which were split into 65 different MIRU-VNTR types. An alternative short-format MIRU-VNTR typing targeting only the four loci with the highest variability values was found to offer an equivalent discriminatory index. Minimum spanning trees using the MIRU-VNTR data showed the hypothetical evolution of an apparent clonal group. MIRU-VNTR analysis was also applied to the isolates of 176 animals from 15 farms infected by M. bovis SB0121; in 10 farms, the analysis revealed the coexistence of two to five different MIRU types differing in one to six loci, which highlights the frequency of undetected heterogeneity.

An inter-laboratory validation of a real time PCR assay to measure host excretion of bacterial pathogens, particularly of Mycobacterium bovis

Advances in the diagnosis of Mycobacterium bovis infection in wildlife hosts may benefit the development of sustainable approaches to the management of bovine tuberculosis in cattle. In the present study, three laboratories from two different countries participated in a validation trial to evaluate the reliability and reproducibility of a real time PCR assay in the detection and quantification of M. bovis from environmental samples. The sample panels consisted of negative badger faeces spiked with a dilution series of M. bovis BCG Pasteur and of field samples of faeces from badgers of unknown infection status taken from badger latrines in areas with high and low incidence of bovine TB (bTB) in cattle. Samples were tested with a previously optimised methodology. The experimental design involved rigorous testing which highlighted a number of potential pitfalls in the analysis of environmental samples using real time PCR. Despite minor variation between operators and laboratories, the validation study demonstrated good concordance between the three laboratories: on the spiked panels, the test showed high levels of agreement in terms of positive/negative detection, with high specificity (100%) and high sensitivity (97%) at levels of 10(5) cells g(-1) and above. Quantitative analysis of the data revealed low variability in recovery of BCG cells between laboratories and operators. On the field samples, the test showed high reproducibility both in terms of positive/negative detection and in the number of cells detected, despite low numbers of samples identified as positive by any laboratory. Use of a parallel PCR inhibition control assay revealed negligible PCR-interfering chemicals co-extracted with the DNA. This is the first example of a multi-laboratory validation of a real time PCR assay for the detection of mycobacteria in environmental samples. Field studies are now required to determine how best to apply the assay for population-level bTB surveillance in wildlife.

Limitations of spoligotyping and variable-number tandem-repeat typing for molecular tracing of Mycobacterium bovis in a high-diversity setting

This study describes the attempt to trace the first Mycobacterium bovis outbreak in alpacas (Lama pacos) in Spain by spoligotyping and variable-number tandem-repeat (VNTR) analysis. Due to high genotype diversity, no matching source was identified, but local expansion of a clonal group was found and its significance for molecular tracing is discussed.

African 2, a clonal complex of Mycobacterium bovis epidemiologically important in East Africa

We have identified a clonal complex of Mycobacterium bovis isolated at high frequency from cattle in Uganda, Burundi, Tanzania, and Ethiopia. We have named this related group of M. bovis strains the African 2 (Af2) clonal complex of M. bovis. Af2 strains are defined by a specific chromosomal deletion (RDAf2) and can be identified by the absence of spacers 3 to 7 in their spoligotype patterns. Deletion analysis of M. bovis isolates from Algeria, Mali, Chad, Nigeria, Cameroon, South Africa, and Mozambique did not identify any strains of the Af2 clonal complex, suggesting that this clonal complex of M. bovis is localized in East Africa. The specific spoligotype pattern of the Af2 clonal complex was rarely identified among isolates from outside Africa, and the few isolates that were found and tested were intact at the RDAf2 locus. We conclude that the Af2 clonal complex is localized to cattle in East Africa. We found that strains of the Af2 clonal complex of M. bovis have, in general, four or more copies of the insertion sequence IS6110, in contrast to the majority of M. bovis strains isolated from cattle, which are thought to carry only one or a few copies.

Tuberculosis in alpacas (Lama pacos) caused by Mycobacterium bovis

We report three cases of tuberculosis in alpacas from Spain caused by Mycobacterium bovis. The animals revealed two different lesional patterns. Mycobacterial culture and PCR assay yielded positive results for M. bovis. Molecular typing of the isolates identified spoligotype SB0295 and identical variable-number tandem repeat (VNTR) allele sizes.

Rapid identification and differentiation of Mycobacterium avium subspecies paratuberculosis types by use of real-time PCR and high-resolution melt analysis of the MAP1506 locus

High-resolution melt (HRM) analysis can identify sequence polymorphisms by comparing the melting curves of amplicons generated by real-time PCR amplification. We describe the application of this technique to identify Mycobacterium avium subspecies paratuberculosis types I, II, and III. The HRM approach was based on type-specific nucleotide sequences in MAP1506, a member of the PPE (proline-proline-glutamic acid) gene family.

Occurrence of Mycobacterium avium subspecies paratuberculosis across host species and European countries with evidence for transmission between wildlife and domestic ruminants

BACKGROUND Mycobacterium avium subspecies paratuberculosis (Map) causes an infectious chronic enteritis (paratuberculosis or Johne's disease) principally of ruminants. The epidemiology of Map is poorly understood, particularly with respect to the role of wildlife reservoirs and the controversial issue of zoonotic potential (Crohn's disease). Genotypic discrimination of Map isolates is pivotal to descriptive epidemiology and resolving these issues. This study was undertaken to determine the genetic diversity of Map, enhance our understanding of the host range and distribution and assess the potential for interspecies transmission. RESULTS 164 Map isolates from seven European countries representing 19 different host species were genotyped by standardized IS900--restriction fragment length polymorphism (IS900-RFLP), pulsed-field gel electrophoresis (PFGE), amplified fragment length polymorphisms (AFLP) and mycobacterial interspersed repeat unit-variable number tandem repeat (MIRU-VNTR) analyses. Six PstI and 17 BstEII IS900-RFLP, 31 multiplex [SnaBI-SpeI] PFGE profiles and 23 MIRU-VNTR profiles were detected. AFLP gave insufficient discrimination of isolates for meaningful genetic analysis. Point estimates for Simpson's index of diversity calculated for the individual typing techniques were in the range of 0.636 to 0.664 but a combination of all three methods increased the discriminating power to 0.879, sufficient for investigating transmission dynamics. Two predominant strain types were detected across Europe with all three typing techniques. Evidence for interspecies transmission between wildlife and domestic ruminants on the same property was demonstrated in four cases, between wildlife species on the same property in two cases and between different species of domestic livestock on one property. CONCLUSION The results of this study showed that it is necessary to use multiple genotyping techniques targeting different sources of genetic variation to obtain the level of discrimination necessary to investigate transmission dynamics and trace the source of Map infections. Furthermore, the combination of genotyping techniques may depend on the geographical location of the population to be tested. Identical genotypes were obtained from Map isolated from different host species co-habiting on the same property strongly suggesting that interspecies transmission occurs. Interspecies transmission of Map between wildlife species and domestic livestock on the same property provides further evidence to support a role for wildlife reservoirs of infection.

Single nucleotide polymorphisms in the IS900 sequence of Mycobacterium avium subsp. paratuberculosis are strain type specific

Insertion sequence IS900 is used as a target for the identification of Mycobacterium avium subsp. paratuberculosis. Previous reports have revealed single nucleotide polymorphisms within IS900. This study, which analyzed the IS900 sequences of a panel of isolates representing M. avium subsp. paratuberculosis strain types I, II, and III, revealed conserved type-specific polymorphisms that could be utilized as a tool for diagnostic and epidemiological purposes.

Discovery of stable and variable differences in the Mycobacterium avium subsp. paratuberculosis type I, II, and III genomes by pan-genome microarray analysis

Mycobacterium avium subsp. paratuberculosis is an important animal pathogen widely disseminated in the environment that has also been associated with Crohn's disease in humans. Three M. avium subsp. paratuberculosis genomotypes are recognized, but genomic differences have not been fully described. To further investigate these potential differences, a 60-mer oligonucleotide microarray (designated the MAPAC array), based on the combined genomes of M. avium subsp. paratuberculosis (strain K-10) and Mycobacterium avium subsp. hominissuis (strain 104), was designed and validated. By use of a test panel of defined M. avium subsp. paratuberculosis strains, the MAPAC array was able to identify a set of large sequence polymorphisms (LSPs) diagnostic for each of the three major M. avium subsp. paratuberculosis types. M. avium subsp. paratuberculosis type II strains contained a smaller genomic complement than M. avium subsp. paratuberculosis type I and M. avium subsp. paratuberculosis type III genomotypes, which included a set of genomic regions also found in M. avium subsp. hominissuis 104. Specific PCRs for genes within LSPs that differentiated M. avium subsp. paratuberculosis types were devised and shown to accurately screen a panel (n = 78) of M. avium subsp. paratuberculosis strains. Analysis of insertion/deletion region INDEL12 showed deletion events causing a reduction in the complement of mycobacterial cell entry genes in M. avium subsp. paratuberculosis type II strains and significantly altering the coding of a major immunologic protein (MPT64) associated with persistence and granuloma formation. Analysis of MAPAC data also identified signal variations in several genomic regions, termed variable genomic islands (vGIs), suggestive of transient duplication/deletion events. vGIs contained significantly low GC% and were immediately flanked by insertion sequences, integrases, or short inverted repeat sequences. Quantitative PCR demonstrated that variation in vGI signals could be associated with colony growth rate and morphology.