Functional analysis of a clonal deletion in an epidemic strain of Mycobacterium bovis reveals a role in lipid metabolism

Enhanced analysis of the genomic diversity of Mycobacterium bovis in Great Britain to aid control of bovine tuberculosis

Bovine tuberculosis (bTB) is an endemic disease in Great Britain (GB) that affects mainly cattle but also other livestock and wild mammal species, leading to significant economic and social impact. Traditional genotyping of Mycobacterium bovis (M. bovis) isolates, which cause bTB, had been used routinely since the late 1990s as the main resource of genetic information in GB to describe their population and to understand their epidemiology. Since 2017, whole-genome sequencing (WGS) has been implemented on M. bovis isolates collected during routine surveillance. In this study, we analysed genome sequences from 3,052 M. bovis isolates from across GB to characterise their diversity and population structure in more detail. Our findings show that the M. bovis population in GB, based on WGS, is more diverse than previously indicated by traditional genotyping and can be divided into seven major clades, with one of them subdivided further into 29 clades that differ from each other by at least 70 single-nucleotide polymorphisms (SNPs). Based on the observed phylogenetic structure, we present a SNP-based classification system that replaces the genotype scheme that had been used until recently in GB. The predicted function and associated processes of the genes harbouring these SNPs are discussed with potential implications for phenotypic/functional differences between the identified clades. At the local scale, we show that WGS provides greater discriminatory power and that it can reveal the origin of infection and associated risk pathways even in areas of high bTB prevalence. The difficulty in determining transmission pathways due to the limited discrimination of isolates by traditional typing methods has compromised bTB control, as without such information it is harder to determine the relative efficacy of potential intervention measures. This study demonstrates that the higher resolution provided by WGS data can improve determination of infection sources and transmission pathways, provide important insights that will inform and shape bTB control policies in GB, as well as improve farm specific advice on interventions that are likely to be effective.

Gut microbiome is associated with personality traits of free-ranging Tibetan macaques (Macaca thibetana)

Recent studies have emphasized that there is a strong link between the gut microbiome and the brain that affects social behavior and personality in animals. However, the interface between personality and the gut microbiome in wild primates remains poorly understood. Here, we used high-throughput sequencing and ethological methods in primate behavioral ecology to investigate the relationship between gut microbiome and personality in Tibetan macaques (Macaca thibetana). The behavioral assessment results indicated three personality dimensions including socialization, shyness, and anxiety. There was significant variation in alpha diversity only for shyness, with a significantly lower alpha diversity indices (including Shannon, Chao1, and PD) for bold individuals than for shy individuals. Using regression models to control for possible confounding factors, we found that the relative abundance of three genera, Akkermansia, Dialister, and Asteroleplasma, was significantly and positively correlated with the sociability scores in the macaques. In addition, Oscillospiraceae exhibited a positive correlation with scores for Shy Dimension. Furthermore, we found that the predicted functional genes for propionate and pyruvate, porphyrin and chlorophyll metabolic pathways related to animal behavior, were significant enriched in shyness group. We propose that the gut microbiome may play an important role in the formation of personality of Tibetan macaques.

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.

Lack of methoxy-mycolates characterizes the geographically restricted lineage 7 of Mycobacterium tuberculosis complex

Lineage 7 (L7) emerged in the phylogeny of the Mycobacterium tuberculosis complex (MTBC) subsequent to the branching of 'ancient' lineage 1 and prior to the Eurasian dispersal of 'modern' lineages 2, 3 and 4. In contrast to the major MTBC lineages, the current epidemiology suggests that prevalence of L7 is highly confined to the Ethiopian population, or when identified outside of Ethiopia, it has mainly been in patients of Ethiopian origin. To search for microbiological factors that may contribute to its restricted distribution, we compared the genome of L7 to the genomes of globally dispersed MTBC lineages. The frequency of predicted functional mutations in L7 was similar to that documented in other lineages. These include mutations characteristic of modern lineages - such as constitutive expression of nitrate reductase - as well as mutations in the VirS locus that are commonly found in ancient lineages. We also identified and characterized multiple lineage-specific mutations in L7 in biosynthesis pathways of cell wall lipids, including confirmed deficiency of methoxy-mycolic acids due to a stop-gain mutation in the mmaA3 gene that encodes a methoxy-mycolic acid synthase. We show that the abolished biosynthesis of methoxy-mycolates of L7 alters the cell structure and colony morphology on selected growth media and impacts biofilm formation. The loss of these mycolic acid moieties may change the host-pathogen dynamic for L7 isolates, explaining the limited geographical distribution of L7 and contributing to further understanding the spread of MTBC lineages across the globe.

Mycobacterial MCE proteins as transporters that control lipid homeostasis of the cell wall

Mammalian cell entry (mce) genes are not only present in genomes of pathogenic mycobacteria, including Mycobacterium tuberculosis (the causative agent of tuberculosis), but also in saprophytic and opportunistic mycobacterial species. MCE are conserved cell-wall proteins encoded by mce operons, which maintain an identical structure in all mycobacteria: two yrbE genes (A and B) followed by six mce genes (A, B, C, D, E and F). Although these proteins are known to participate in the virulence of pathogenic mycobacteria, the presence of the operons in nonpathogenic mycobacteria and other bacteria indicates that they play another role apart from host cell invasion. In this respect, more recent studies suggest that they are functionally similar to ABC transporters and form part of lipid transporters in Actinobacteria. To date, most reviews on mce operons in the literature discuss their role in virulence. However, according to data from transcriptional studies, mce genes, particularly the mce1 and mce4 operons, modify their expression according to the carbon source and upon hypoxia, starvation, surface stress and oxidative stress; which suggests a role of MCE proteins in the response of Mycobacteria to external stressors. In addition to these data, this review also summarizes the studies demonstrating the role of MCE proteins as lipid transporters as well as the relevance of their transport function in the interaction of pathogenic Mycobacteria with the hosts. Altogether, the evidence to date would indicate that MCE proteins participate in the response to the stress conditions that mycobacteria encounter during infection, by participating in the cell wall remodelling and possibly contributing to lipid homeostasis.

Correlations between drug resistance of Beijing/W lineage clinical isolates of Mycobacterium tuberculosis and sublineages: a 2009-2013 prospective study in Xinjiang province, China

BACKGROUND

The prevalence of drug-resistant tuberculosis (TB) in Xinjiang is higher than in other regions of China, and Beijing/W lineage Mycobacterium tuberculosis (MTB) is the dominant strain of MTB in Xinjiang. However, information on multidrug-resistant (MDR) and extensively drug-resistant (XDR) TB, particularly the correlation between MDR and the Beijing/W lineage and the correlation between drug resistance and the Beijing/W sublineage strains, is limited.

MATERIAL/METHODS

We conducted a prospective study to describe the prevalence of MDR/XDR TB, Beijing/W lineage and sublineage strains in Xinjiang in China from 2009 to 2013. All MTB underwent drug susceptibility testing to the first- and second-line anti-tuberculosis drugs. The Beijing/W lineages and sublineages were detected by large-sequence polymorphisms with polymerase chain reaction.

RESULTS

A total of 410 clinical isolates were identified. The overall percentage of MDR and XDR cases in Xinjiang was 13.2% (54/410) and 13.0% (7/54), respectively. Overall, 9.8% (14/143) of the Beijing lineage MTB were MDR patients, and 15.6% (40/257) of the Non-Beijing lineage MTB were MDR patients. In the 143 Beijing MTB lineages, 11.2% isolates were in sublineage 105, 15.4% isolates were in sublineage 207, 69.2% isolates were in sublineage 181, and 4.2% isolates were in sublineage 150. None of the isolates were detected in sublineage 142. Significant differences between the Beijing/W and non-Beijing/W strains were observed regarding INH and EMB resistance, respectively.

CONCLUSIONS

The prevalence of the MDR TB in Xinjiang remains high and imposes challenges for TB control. Four Beijing/W sublineage isolates were observed in Xinjiang. There was no correlation between MDR and the Beijing/W lineage and no correlation between drug resistance and the Beijing/W sublineage strains. Surveillance of the clinical isolates of MTB is recommended to strengthen the identification of MDR/XDR TB and sublineages of the Beijing/W strains.

Pathogen roid rage: cholesterol utilization by Mycobacterium tuberculosis

The ability of science and medicine to control the pathogen Mycobacterium tuberculosis (Mtb) requires an understanding of the complex host environment within which it resides. Pathological and biological evidence overwhelmingly demonstrate how the mammalian steroid cholesterol is present throughout the course of infection. Better understanding Mtb requires a more complete understanding of how it utilizes molecules like cholesterol in this environment to sustain the infection of the host. Cholesterol uptake, catabolism and broader utilization are important for maintenance of the pathogen in the host and it has been experimentally validated to contribute to virulence and pathogenesis. Cholesterol is catabolized by at least three distinct sub-pathways, two for the ring system and one for the side chain, yielding dozens of steroid intermediates with varying biochemical properties. Our ability to control this worldwide infectious agent requires a greater knowledge of how Mtb uses cholesterol to its advantage throughout the course of infection. Herein, the current state of knowledge of cholesterol metabolism by Mtb is reviewed from a biochemical perspective with a focus on the metabolic genes and pathways responsible for cholesterol steroid catabolism.

Genome-level analyses of Mycobacterium bovis lineages reveal the role of SNPs and antisense transcription in differential gene expression

Background

Bovine tuberculosis (bTB) is a disease with major implications for animal welfare and productivity, as well as having the potential for zoonotic transmission. In Great Britain (GB) alone, controlling bTB costs in the region of £100 million annually, with the current control scheme seemingly unable to stop the inexorable spread of infection. One aspect that may be driving the epidemic is evolution of the causative pathogen, Mycobacterium bovis. To understand the underlying genetic changes that may be responsible for this evolution, we performed a comprehensive genome-level analyses of 4 M. bovis strains that encompass the main molecular types of the pathogen circulating in GB.

Results

We have used a combination of genome sequencing, transcriptome analyses, and recombinant DNA technology to define genetic differences across the major M. bovis lineages circulating in GB that may give rise to phenotypic differences of practical importance. The genomes of three M. bovis field isolates were sequenced using Illumina sequencing technology and strain specific differences in gene expression were measured during in vitro growth and in ex vivo bovine alveolar macrophages using a whole genome amplicon microarray and a whole genome tiled oligonucleotide microarray. SNP/small base pair insertion and deletions and gene expression data were overlaid onto the genomic sequence of the fully sequenced strain of M. bovis 2122/97 to link observed strain specific genomic differences with differences in RNA expression.

Conclusions

We show that while these strains show extensive similarities in their genetic make-up and gene expression profiles, they exhibit distinct expression of a subset of genes. We provide genomic, transcriptomic and functional data to show that synonymous point mutations (sSNPs) on the coding strand can lead to the expression of antisense transcripts on the opposing strand, a finding with implications for how we define a 'silent’ nucleotide change. Furthermore, we show that transcriptomic data based solely on amplicon arrays can generate spurious results in terms of gene expression profiles due to hybridisation of antisense transcripts. Overall our data suggest that subtle genetic differences, such as sSNPS, may have important consequences for gene expression and subsequent phenotype.

Shrinking the FadE proteome of Mycobacterium tuberculosis: insights into cholesterol metabolism through identification of an α2β2 heterotetrameric acyl coenzyme A dehydrogenase family

The ability of the pathogen Mycobacterium tuberculosis to metabolize steroids like cholesterol and the roles that these compounds play in the virulence and pathogenesis of this organism are increasingly evident. Here, we demonstrate through experiments and bioinformatic analysis the existence of an architecturally distinct subfamily of acyl coenzyme A (acyl-CoA) dehydrogenase (ACAD) enzymes that are α2β2 heterotetramers with two active sites. These enzymes are encoded by two adjacent ACAD (fadE) genes that are regulated by cholesterol. FadE26-FadE27 catalyzes the dehydrogenation of 3β-hydroxy-chol-5-en-24-oyl-CoA, an analog of the 5-carbon side chain cholesterol degradation intermediate. Genes encoding the α2β2 heterotetrameric ACAD structures are present in multiple regions of the M. tuberculosis genome, and subsets of these genes are regulated by four different transcriptional repressors or activators: KstR1 (also known as KstR), KstR2, Mce3R, and SigE. Homologous ACAD gene pairs are found in other Actinobacteria, as well as Proteobacteria. Their structures and genomic locations suggest that the α2β2 heterotetrameric structural motif has evolved to enable catalysis of dehydrogenation of steroid- or polycyclic-CoA substrates and that they function in four subpathways of cholesterol metabolism.

High throughput phenotypic analysis of Mycobacterium tuberculosis and Mycobacterium bovis strains' metabolism using biolog phenotype microarrays

Tuberculosis is a major human and animal disease of major importance worldwide. Genetically, the closely related strains within the Mycobacterium tuberculosis complex which cause disease are well-characterized but there is an urgent need better to understand their phenotypes. To search rapidly for metabolic differences, a working method using Biolog Phenotype MicroArray analysis was developed. Of 380 substrates surveyed, 71 permitted tetrazolium dye reduction, the readout over 7 days in the method. By looking for ≥5-fold differences in dye reduction, 12 substrates differentiated M. tuberculosis H37Rv and Mycobacterium bovis AF2122/97. H37Rv and a Beijing strain of M. tuberculosis could also be distinguished in this way, as could field strains of M. bovis; even pairs of strains within one spoligotype could be distinguished by 2 to 3 substrates. Cluster analysis gave three clear groups: H37Rv, Beijing, and all the M. bovis strains. The substrates used agreed well with prior knowledge, though an unexpected finding that AF2122/97 gave greater dye reduction than H37Rv with hexoses was investigated further, in culture flasks, revealing that hexoses and Tween 80 were synergistic for growth and used simultaneously rather than in a diauxic fashion. Potential new substrates for growth media were revealed, too, most promisingly N-acetyl glucosamine. Osmotic and pH arrays divided the mycobacteria into two groups with different salt tolerance, though in contrast to the substrate arrays the groups did not entirely correlate with taxonomic differences. More interestingly, these arrays suggested differences between the amines used by the M. tuberculosis complex and enteric bacteria in acid tolerance, with some hydrophobic amino acids being highly effective. In contrast, γ-aminobutyrate, used in the enteric bacteria, had no effect in the mycobacteria. This study proved principle that Phenotype MicroArrays can be used with slow-growing pathogenic mycobacteria and already has generated interesting data worthy of further investigation.

The cytotoxic early protein 77 of mycobacteriophage L5 interacts with MSMEG_3532, an L-serine dehydratase of Mycobacterium smegmatis

Mycobacteriophage L5 is a temperate phage infecting a broad range of mycobacterial species. Upon induction of lytic growth, L5 rapidly switches off host protein synthesis. We have recently identified the mycobacteriophage L5 early protein gp77 as a host shut-off protein that acts growth inhibitory in the mycobacterial host when expressed through the corresponding phage promoter. Here we present data showing that this purified phage protein of unknown function specifically binds to protein MSMEG_3532 when incubated with cell lysates of Mycobacterium smegmatis. This interaction was confirmed by pull-down assays using purified MSMEG_3532 as bait which co-purified with gp77. The amino acid sequence of MSMEG_3532 is nearly identical to that of threonine dehydratases, serine dehydratases and an L-threo-3-hydroxyaspartate dehydratase. An enzymatic assay identified this host protein as a pyridoxal-5'-phosphate-dependent L-serine dehydratase (SdhA) which converts L-serine to pyruvate. This is the first biochemical characterization of a SdhA derived from mycobacteria. Though the addition of purified gp77 to the established in vitro assay had no influence on SdhA activity at a saturating L-serine concentration, the specific interaction of phage protein and dehydratase in vivo may well have a role in altering the amino acid pool or the products of amino acid metabolism in favour of phage maturation.

Global effects of inactivation of the pyruvate kinase gene in the Mycobacterium tuberculosis complex

To better understand the global effects of "natural" lesions in genes involved in the pyruvate metabolism in Mycobacterium bovis, null mutations were made in the Mycobacterium tuberculosis H37Rv ald and pykA genes to mimic the M. bovis situation. Like M. bovis, the M. tuberculosis DeltapykA mutant yielded dysgonic colonies on solid medium lacking pyruvate, whereas colony morphology was eugonic on pyruvate-containing medium. Global effects of the loss of the pykA gene, possibly underlying colony morphology, were investigated by using proteomics on cultures grown in the same conditions. The levels of Icd2 increased and those of Icl and PckA decreased in the DeltapykA knockout. Proteomics suggested that the synthesis of enzymes involved in fatty acid and lipid biosynthesis were decreased, whereas those involved in beta-oxidation were increased in the M. tuberculosis DeltapykA mutant, as confirmed by direct assays for these activities. Thus, the loss of pykA from M. tuberculosis results in fatty acids being used principally for energy production, in contrast to the situation in the host when carbon from fatty acids is conserved through the glyoxylate cycle and gluconeogenesis; when an active pykA gene was introduced into M. bovis, the opposite effects occurred. Proteins involved in oxidative stress-AhpC, KatG, and SodA-showed increased synthesis in the DeltapykA mutant, and iron-regulated proteins were also affected. Ald levels were decreased in the DeltapykA knockout, explaining why an M. tuberculosis DeltapykA Deltaald double mutant showed little additional phenotypic effect. Overall, these data show that the loss of the pykA gene has powerful, global effects on proteins associated with central metabolism.