Functional analysis of the Mycobacterium bovis AF2122/97 PhoPR system.

The PhoPR system is a master regulator in Mycobacterium tuberculosis. A key difference between M. tuberculosis and Mycobacterium bovis is a G71I substitution in the M. bovis PhoR orthologue. Functional studies of the M. bovis PhoPR system have generated conflicting findings, with some research suggesting that the M. bovis PhoPR is defective while others indicate it is functional. We sought to revisit the functionality of the M. bovis PhoPR system. To address this, we constructed a phoPR mutant in the reference strain M. bovis AF2122/97. We employed a combination of growth assays and transcriptomics analyses to assess the phenotype of the mutant vs wild type and complemented strains. We found that the M. bovis AF2122/97 ΔphoPR mutant showed a growth defect on solid and liquid media compared to the wild type and complemented strains. The transcriptome of the M. bovis AF2122/97 ΔphoPR mutant was also altered as compared to wild type, including differential expression of genes involved in lipid metabolism and secretion. Our work provides further insight into the activity of PhoPR in M. bovis and underlines the importance of the PhoPR system as a master regulator of gene expression in the Mycobacterium tuberculosis complex.

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.

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

[This corrects the article DOI: 10.3389/fmicb.2022.817952.].

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

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

Evidence for local and international spread of Mycobacterium avium subspecies paratuberculosis through whole genome sequencing of isolates from the island of Ireland.

We describe application of whole genome sequencing (WGS) to a collection of 197 Mycobacterium avium subsp paratuberculosis (MAP) isolates gathered from 122 cattle herds across 27 counties of the island of Ireland. We compare WGS to MAP diversity quantified using mycobacterial interspersed random unit - variable number tandem repeats (MIRU-VNTR). While MIRU-VNTR showed only two major types, WGS could split the 197 isolates into eight major groups. We also found six isolates corresponding to INMV 13, a novel MIRU-VNTR type for Ireland. Evidence for dispersal of MAP across Ireland via cattle movement could be discerned from the data, with mixed infections present in several herds. Furthermore, comparisons of MAP WGS data from Ireland to data from Great Britain and continental Europe revealed many instances of close genetic similarity and hence evidence for international transmission of infection. BEAST MASCOT structured coalescent analyses, with relaxed and strict molecular clocks, estimated the substitution rate to be 0.10-0.13 SNPs/site/year and disclosed greater transitions per lineage per year from Europe to Ireland, indicating transmission into Ireland. Our work therefore reveals new insight into the seeding of MAP infection across Ireland, highlighting how WGS can inform policy formulation to ultimately control MAP transmission at local, national and international scales.

Performance and Agreement Between WGS Variant Calling Pipelines Used for Bovine Tuberculosis Control: Toward International Standardization.

Whole genome sequencing (WGS) and allied variant calling pipelines are a valuable tool for the control and eradication of infectious diseases, since they allow the assessment of the genetic relatedness of strains of animal pathogens. In the context of the control of tuberculosis (TB) in livestock, mainly caused by Mycobacterium bovis, these tools offer a high-resolution alternative to traditional molecular methods in the study of herd breakdown events. However, despite the increased use and efforts in the standardization of WGS methods in human tuberculosis around the world, the application of these WGS-enabled approaches to control TB in livestock is still in early development. Our study pursued an initial evaluation of the performance and agreement of four publicly available pipelines for the analysis of M. bovis WGS data (vSNP, SNiPgenie, BovTB, and MTBseq) on a set of simulated Illumina reads generated from a real-world setting with high TB prevalence in cattle and wildlife in the Republic of Ireland. The overall performance of the evaluated pipelines was high, with recall and precision rates above 99% once repeat-rich and problematic regions were removed from the analyses. In addition, when the same filters were applied, distances between inferred phylogenetic trees were similar and pairwise comparison revealed that most of the differences were due to the positioning of polytomies. Hence, under the studied conditions, all pipelines offer similar performance for variant calling to underpin real-world studies of M. bovis transmission dynamics.

epitopepredict: a tool for integrated MHC binding prediction.

A key step in the cellular adaptive immune response is the presentation of antigens to T cells. Computational prediction of T cell epitopes has many applications in vaccine design and immuno-diagnostics. This is the basis of immunoinformatics, which allows in silico screening of peptides before experiments are performed. With the availability of whole genomes for many microbial species it is now feasible to computationally screen whole proteomes for candidate peptides. epitopepredict is a programmatic framework and command line tool designed to aid this process. It provides access to multiple binding prediction algorithms under a single interface and scales for whole genomes using multiple target MHC alleles. A web interface is provided to assist visualization and filtering of the results. The software is freely available under an open-source license from https://github.com/dmnfarrell/epitopepredict.

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.

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.

Identification of microRNAs in bovine faeces and their potential as biomarkers of Johne's Disease.

Extracellular microRNAs (miRNAs) are detectable in the peripheral blood and have been touted as potential biomarkers for a range of maladies. The presence and biomarker potential of miRNAs in other biofluids has been less thoroughly explored, particularly in the veterinary realm. Faecal miRNAs are a case in point; while they have been identified largely in rodents and humans, they have not been reported in cattle but may have prognostic or diagnostic value for Johne's Disease (JD) in cattle, a chronic granulomatous inflammation of the ileum caused by Mycobacterium avium subspecies paratuberculosis (MAP). The aim of this study was thus to characterise the bovine faecal miRNome and to determine the utility of these transcripts as biomarkers for JD. Real-time PCR arrays consisting of 752 miRNA targets, optimised for detection of human miRNA, were used to screen RNA purified from faecal samples obtained from confirmed JD clinical cases vs. healthy controls. Two hundred and fifty-eight miRNAs were detected in bovine faeces, three of which are potentially novel orthologs of known human miRNAs. Differential abundance of three miRNA was evident in animals with clinical JD as compared to healthy controls. Our study has therefore identified a variety of miRNAs in bovine faeces and has demonstrated their utility in differentiating healthy animals from those with late-stage JD, providing potential biomarkers for MAP infection and disease progression.

Draft Genome Sequences of Three Novel Acinetobacter Isolates from an Irish Commercial Pig Farm.

Acinetobacter species are important in the emergence and spread of antimicrobial resistance (AMR), which threatens human and animal health worldwide. Here, we present the draft genome sequences of three Acinetobacter species strains (RF14B, RF15A, and RF15B) isolated from pig feces and the floor of a pig hospital pen in Ireland.

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

We report here an update to the reference genome sequence of the bovine tuberculosis bacillus Mycobacterium bovis AF2122/97, generated using an integrative multiomics approach. The update includes 42 new coding sequences (CDSs), 14 modified annotations, 26 single-nucleotide polymorphism (SNP) corrections, and disclosure that the RD900 locus, previously described as absent from the genome, is in fact present.

Integrated computational prediction and experimental validation identifies promiscuous T cell epitopes in the proteome of Mycobacterium bovis.

The discovery of novel antigens is an essential requirement in devising new diagnostics or vaccines for use in control programmes against human tuberculosis (TB) and bovine tuberculosis (bTB). Identification of potential epitopes recognised by CD4+ T cells requires prediction of peptide binding to MHC class-II, an obligatory prerequisite for T cell recognition. To comprehensively prioritise potential MHC-II-binding epitopes from Mycobacterium bovis, the agent of bTB and zoonotic TB in humans, we integrated three binding prediction methods with the M. bovisproteome using a subset of human HLA alleles to approximate the binding of epitope-containing peptides to the bovine MHC class II molecule BoLA-DRB3. Two parallel strategies were then applied to filter the resulting set of binders: identification of the top-scoring binders or clusters of binders. Our approach was tested experimentally by assessing the capacity of predicted promiscuous peptides to drive interferon-γ secretion from T cells of M. bovis infected cattle. Thus, 376 20-mer peptides, were synthesised (270 predicted epitopes, 94 random peptides with low predictive scores and 12 positive controls of known epitopes). The results of this validation demonstrated significant enrichment (>24 %) of promiscuously recognised peptides predicted in our selection strategies, compared with randomly selected peptides with low prediction scores. Our strategy offers a general approach to the identification of promiscuous epitopes tailored to target populations where there is limited knowledge of MHC allelic diversity.

Analysis of Biobanked Serum from a Mycobacterium avium subsp paratuberculosis Bovine Infection Model Confirms the Remarkable Stability of Circulating miRNA Profiles and Defines a Bovine Serum miRNA Repertoire.

Johne's Disease (JD) is a chronic enteritis of ruminants caused by Mycobacterium avium subspecies paratuberculosis (MAP). Current disease control strategies are hampered by the lack of sensitive and specific diagnostic modalities. Therefore, novel diagnostic and prognostic tools are needed, and circulating microRNAs (miRNAs) may hold potential in this area. The aims of this study were twofold: (i) to address the stability of miRNA in bovine sera from biobanked samples, and (ii) to assess the potential of miRNAs as biomarkers for JD disease progression. To address these aims we used bovine sera from an experimental MAP infection model that had been stored at -20°C for over a decade, allowing us to also assess the stability of miRNA profiles in biobanked serum samples through comparison with fresh sera. Approximately 100-200 intact miRNAs were identified in each sample with 83 of these being consistently detected across all 57 samples. The miRNA profile of the biobanked sera stored at -20°C for over 10 years was highly similar to the profile of <1 year-old sera stored at -80°C, with an overlap of 73 shared miRNAs. IsomiR analysis also indicated a distinct bovine serum-specific isomiR profile as compared to previously reported bovine macrophage miRNA profiles. To explore the prognostic potential of miRNA profiles cattle defined as seropositive for anti-MAP antibodies (n = 5) were compared against seronegative cattle (n = 7). No significant differential expressed miRNAs were detected at either the early (6 months) or late (43, 46 and 49 months) intervals (FDR≤0.05, fold-change≥1.5) across seropositive or seronegative animals. However, comparing pre-infection sera to the early and late time-points identified increased miR-29a and miR-92b abundance (2-fold) that may be due to blood-cell population changes over time (P<0.001). In conclusion our study has demonstrated that bovine circulating miRNAs retain their integrity under long-term sub-optimal storage temperatures opening the way for increased miRNA analyses from biobanked samples for a range of infectious and non-infectious diseases.

Epitopemap: a web application for integrated whole proteome epitope prediction.

BACKGROUND: Predictions of MHC binding affinity are commonly used in immunoinformatics for T cell epitope prediction. There are multiple available methods, some of which provide web access. However there is currently no convenient way to access the results from multiple methods at the same time or to execute predictions for an entire proteome at once. RESULTS: We designed a web application that allows integration of multiple epitope prediction methods for any number of proteins in a genome. The tool is a front-end for various freely available methods. Features include visualisation of results from multiple predictors within proteins in one plot, genome-wide analysis and estimates of epitope conservation. CONCLUSIONS: We present a self contained web application, Epitopemap, for calculating and viewing epitope predictions with multiple methods. The tool is easy to use and will assist in computational screening of viral or bacterial genomes.

The Identification of Circulating MiRNA in Bovine Serum and Their Potential as Novel Biomarkers of Early Mycobacterium avium subsp paratuberculosis Infection.

Mycobacterium avium subspecies paratuberculosis (MAP) is the aetiological agent of Johne's disease (JD), a chronic enteritis in ruminants that causes substantial economic loses to agriculture worldwide. Current diagnostic assays are hampered by low sensitivity and specificity that seriously complicate disease control; a new generation of diagnostic and prognostic assays are therefore urgently needed. Circulating microRNAs (miRNAs) have been shown to have significant potential as novel biomarkers for a range of human diseases, but their potential application in the veterinary sphere has been less well characterised. The aim of this study was therefore to apply RNA-sequencing approaches to serum from an experimental JD infection model as a route to identify novel diagnostic and prognostic miRNA biomarkers. Sera from experimental MAP-challenged calves (n = 6) and age-matched controls (n = 6) were used. We identified a subset of known miRNAs from bovine serum across all samples, with approximately 90 being at potentially functional abundance levels. The majority of known bovine miRNAs displayed multiple isomiRs that differed from the canonical sequences. Thirty novel miRNAs were identified after filtering and were found within sera from all animals tested. No significant differential miRNA expression was detected when comparing sera from MAP-challenged animals to their age-matched controls at six-month's post-infection. However, comparing sera from pre-infection bleeds to six-month's post-infection across all 12 animals did identify increased miR-205 (2-fold) and decreased miR-432 (2-fold) within both challenged and control groups, which suggests changes in circulating miRNA profiles due to ageing or development (P<0.00001). In conclusion our study has identified a range of novel miRNA in bovine serum, and shown the utility of small RNA sequencing approaches to explore the potential of miRNA as novel biomarkers for infectious disease in cattle.

Protein dielectric constants determined from NMR chemical shift perturbations.

Understanding the connection between protein structure and function requires a quantitative understanding of electrostatic effects. Structure-based electrostatic calculations are essential for this purpose, but their use has been limited by a long-standing discussion on which value to use for the dielectric constants (ε(eff) and ε(p)) required in Coulombic and Poisson-Boltzmann models. The currently used values for ε(eff) and ε(p) are essentially empirical parameters calibrated against thermodynamic properties that are indirect measurements of protein electric fields. We determine optimal values for ε(eff) and ε(p) by measuring protein electric fields in solution using direct detection of NMR chemical shift perturbations (CSPs). We measured CSPs in 14 proteins to get a broad and general characterization of electric fields. Coulomb's law reproduces the measured CSPs optimally with a protein dielectric constant (ε(eff)) from 3 to 13, with an optimal value across all proteins of 6.5. However, when the water-protein interface is treated with finite difference Poisson-Boltzmann calculations, the optimal protein dielectric constant (ε(p)) ranged from 2 to 5 with an optimum of 3. It is striking how similar this value is to the dielectric constant of 2-4 measured for protein powders and how different it is from the ε(p) of 6-20 used in models based on the Poisson-Boltzmann equation when calculating thermodynamic parameters. Because the value of ε(p) = 3 is obtained by analysis of NMR chemical shift perturbations instead of thermodynamic parameters such as pK(a) values, it is likely to describe only the electric field and thus represent a more general, intrinsic, and transferable ε(p) common to most folded proteins.

DataPipeline: automated importing and fitting of large amounts of biophysical data.

Raw data from experiments across the biological sciences comes in a large variety of text formats. In small or medium sized laboratories researchers often use an assorted collection of software to interpret, fit, and visualize their data. The spreadsheet is commonly the core component of such a workflow. The limitations of such programs for large amounts of heterogeneous data can be frustrating. We report the construction of DataPipeline, a desktop and command-line application that automates the tasks of importing, fitting, and plotting of text-based data. The software is designed to simplify the process of importing text data from various sources using simple configuration files to describe raw file formats. Once imported, curve fitting can be performed using custom fitting models designed by the user inside the application. Fitted parameters can be grouped together as new datasets to be fitted to other models and experimental uncertainties propagated to give error estimates. This software will be useful for processing of data from high through-put biological experiments or for rapid visualization of pilot data without the need for a chain of different programs to carry out each step. DataPipeline and source code is available under an open source license. The software can be freely downloaded at http://code.google.com/p/peat/downloads/list.

Toward fast determination of protein stability maps: experimental and theoretical analysis of mutants of a Nocardiopsis prasina serine protease.

The stability of serine proteases is of major importance for their application in industrial processes. Here we study the determinants of the stability of a Nocardiopsis prasina serine protease using fast residual activity assays, a feature classification algorithm, and structure-based energy calculation algorithms for 121 micropurified mutant enzyme clones containing multiple point mutations. Using a multivariate regression analysis, we deconvolute the data for the mutant clones and find that mutations of residues Asn47 and Pro124 are deleterious to the stability of the enzyme. Both of these residues are situated in loops that are known to be important for the stability of the highly homologous α-lytic protease. Structure-based energy calculations with PEATSA give a good general agreement with the trend of experimentally measured values but also identify a number of clones that the algorithm fails to predict correctly. We discuss the significance of the results in relation to the structure and function of closely related proteases, comment on the optimal experimental design when performing high-throughput experiments for characterizing the determinants of protein stability, and discuss the performance of structure-based energy calculations with complex data sets such as the one presented here.