Alignment of multiple complete genomes suggests that gene rearrangements may contribute towards the speciation of Mycobacteria

Characterisation of IS1311 in Mycobacterium avium subspecies paratuberculosis genomes: Typing, continental clustering, microbial evolution and host adaptation

Johne's disease (JD), caused by Mycobacterium avium subspecies paratuberculosis (MAP) is a global burden for livestock producers and has an association with Crohn's disease in humans. Within MAP there are two major lineages, S/Type I/TypeIII and C/Type II, that vary in phenotype including culturability, host preference and virulence. These lineages have been identified using the IS1311 element, which contains a conserved, single nucleotide polymorphism. IS1311 and the closely related IS1245 element belong to the IS256 family of insertion sequences, are dispersed throughout M. avium taxa but remain poorly characterised. To investigate the distribution and diversity of IS1311 in MAP, 805 MAP genomes were collated from public databases. IS1245 was absent, while IS1311 sequence, copy number and insertion loci were conserved between MAP S lineages and varied within the MAP C lineage. One locus was specific to the S strains, which contained nine IS1311 copies. In contrast, C strains contained either seven or eight IS1311 loci. Most insertion loci were associated with the boundaries of homologous regions that had undergone genome rearrangement between the MAP lineages, suggesting that this sequence may be a driver of recombination. Phylogenomic geographic clustering of MAP subtypes was demonstrated for the first time, at continental scale, and indicated that there may have been recent MAP transmission between Europe and North America, in contrast to Australia where importation of live ruminants is generally prohibited. This investigation confirmed the utility of IS1311 typing in epidemiological studies and resolved anomalies in past studies. The results shed light on potential mechanisms of niche/host adaptation, virulence of MAP and global transmission dynamics.

Genomic structural plasticity of rodent-associated Bartonella in nature

Rodent‐associated Bartonella species have shown a remarkable genetic diversity and pathogenic potential. To further explore the extent of the natural intraspecific genomic variation and its potential role as an evolutionary driver, we focused on a single genetically diverse Bartonella species, Bartonella krasnovii, which circulates among gerbils and their associated fleas. Twenty genomes from 16 different B. krasnovii genotypes were fully characterized through a genome sequencing assay (using short and long read sequencing), pulse field gel electrophoresis (PFGE), and PCR validation. Genomic analyses were performed in comparison to the B. krasnovii strain OE 1–1. While, single nucleotide polymorphism represented only a 0.3% of the genome variation, structural diversity was identified in these genomes, with an average of 51 ± 24 structural variation (SV) events per genome. Interestingly, a large proportion of the SVs (>40%) was associated with prophages. Further analyses revealed that most of the SVs, and prophage insertions were found at the chromosome replication termination site (ter), suggesting this site as a plastic zone of the B. krasnovii chromosome. Accordingly, six genomes were found to be unbalanced, and essential genes near the ter showed a shift between the leading and lagging strands, revealing the SV effect on these genomes. In summary, our findings demonstrate the extensive genomic diversity harbored by wild B. krasnovii strains and suggests that its diversification is initially promoted by structural changes, probably driven by phages. These events may constantly feed the system with novel genotypes that ultimately lead to inter‐ and intraspecies competition and adaptation.

Tuberculosis caused by Mycobacterium africanum: Knowns and unknowns

Tuberculosis (TB), one of the deadliest threats to human health, is mainly caused by 2 highly related and human-adapted bacteria broadly known as Mycobacterium tuberculosis and Mycobacterium africanum. Whereas M. tuberculosis is widely spread, M. africanum is restricted to West Africa, where it remains a significant cause of tuberculosis. Although several differences have been identified between these 2 pathogens, M. africanum remains a lot less studied than M. tuberculosis. Here, we discuss the genetic, phenotypic, and clinical similarities and differences between strains of M. tuberculosis and M. africanum. We also discuss our current knowledge on the immune response to M. africanum and how it possibly articulates with distinct disease progression and with the geographical restriction attributed to this pathogen. Understanding the functional impact of the diversity existing in TB-causing bacteria, as well as incorporating this diversity in TB research, will contribute to the development of better, more specific approaches to tackle TB.

Prevalence, Risk Factors, and Diagnostic Efficacy of Bovine Tuberculosis in Slaughtered Animals at the Chiang Mai Municipal Abattoir, Thailand

This study aimed to (1) investigate the prevalence of bovine tuberculosis (bTB) in slaughtered animals at the Chiang Mai Municipal abattoir in Chiang Mai, Thailand; (2) identify animal-level risk factors for bTB at the abattoir; and (3) evaluate the performance of techniques for bTB detection at the abattoir. From April 2020 to March 2021, 161 animals registered for slaughter were randomly selected for the study. Animal data including age, sex, species, body condition scores, and origins of the animals were collected. Meat inspection was performed by a trained meat inspector. Tissue samples of the lung, liver, and lymph nodes were collected for histopathological diagnosis and polymerase chain reaction (PCR) detection of Mycobacteria and specifically Mycobacterium bovis. The prevalence of bTB during meat inspection and PCR was calculated separately. Animal-level factors affecting bTB were determined using multivariate logistic regression analysis. The performance of meat inspection and PCR was evaluated using a Bayesian approach. The prevalence of bTB was 12.4% (20/161) and 34.8% (56/161) when the disease was diagnosed using meat inspection and PCR, respectively. Buffaloes had a significantly higher risk of being identified as bTB-positive using PCR compared to beef cattle (odds ratio = 2.19; confidence interval = 1.11–4.30). The median of posterior estimates of sensitivity (Se) and specificity (Sp) to detect bTB using meat inspection were 20.8% [95% posterior probability interval (PPI) = 9.1–36.5%] and 87.8% (95% PPI = 79.6–95.4%), respectively. The medians of the posterior estimates of Se and Sp for PCR were 88.6% (95% PPI = 70.5–98.3%) and 94.4% (95% PPI = 84.7–98.8%), respectively. These findings demonstrate that bTB is highly prevalent among slaughtered animals. PCR can be used as an ancillary test for bTB surveillance at abattoirs in Thailand.

Human-to-Cattle Mycobacterium tuberculosis Complex Transmission in the United States

The Mycobacterium tuberculosis complex (MTBC) species includes both M. tuberculosis, the primary cause of human tuberculosis (TB), and M. bovis, the primary cause of bovine tuberculosis (bTB), as well as other closely related Mycobacterium species. Zoonotic transmission of M. bovis from cattle to humans was recognized more than a century ago, but transmission of MTBC species from humans to cattle is less often recognized. Within the last decade, multiple published reports from around the world describe human-to-cattle transmission of MTBC. Three probable cases of human-to-cattle MTBC transmission have occurred in the United States since 2013. In the first case, detection of active TB disease (M. bovis) in a dairy employee in North Dakota prompted testing and ultimate detection of bTB infection in the dairy herd. Whole genome sequencing (WGS) demonstrated a match between the bTB strain in the employee and an infected cow. North Dakota animal and public health officials concluded that the employee's infection was the most likely source of disease introduction in the dairy. The second case involved a Wisconsin dairy herd with an employee diagnosed with TB disease in 2015. Subsequently, the herd was tested twice with no disease detected. Three years later, a cow originating from this herd was detected with bTB at slaughter. The strain in the slaughter case matched that of the past employee based on WGS. The third case was a 4-month-old heifer calf born in New Mexico and transported to Texas. The calf was TB tested per Texas entry requirements and found to have M. tuberculosis. Humans are the suspected source of M. tuberculosis in cattle; however, public health authorities were not able to identify an infected human associated with the cattle operation. These three cases provide strong evidence of human-to-cattle transmission of MTBC organisms and highlight human infection as a potential source of introduction of MTBC into dairy herds in the United States. To better understand and address the issue, a multisectoral One Health approach is needed, where industry, public health, and animal health work together to better understand the epidemiology and identify preventive measures to protect human and animal health.

Retrospective Analysis of Archived Pyrazinamide Resistant Mycobacterium tuberculosis Complex Isolates from Uganda-Evidence of Interspecies Transmission

The contribution of Mycobacterium bovis to the proportion of tuberculosis cases in humans is unknown. A retrospective study was undertaken on archived Mycobacterium tuberculosis complex (MTBC) isolates from a reference laboratory in Uganda to identify the prevalence of human M. bovis infection. A total of 5676 isolates maintained in this repository were queried and 136 isolates were identified as pyrazinamide resistant, a hallmark phenotype of M. bovis. Of these, 1.5% (n = 2) isolates were confirmed as M. bovis by using regions of difference PCR analysis. The overall size of whole genome sequences (WGSs) of these two M. bovis isolates were ~4.272 Mb (M. bovis Bz_31150 isolated from a captive chimpanzee) and 4.17 Mb (M. bovis B2_7505 from a human patient), respectively. Alignment of these genomes against 15 MTBC genome sequences revealed 7248 single nucleotide polumorphisms (SNPs). Theses SNPs were used for phylogenetic analysis that indicated a strong relationship between M. bovis and the chimpanzee isolate (Bz_31150) while the other M. bovis genome from the human patient (B2_7505) analyzed did not cluster with any M. bovis or M. tuberculosis strains. WGS analysis also revealed multidrug resistance genotypes; these genomes revealed pncA mutations at positions H57D in Bz_31150 and B2_7505. Phenotypically, B2_7505 was an extensively drug-resistant strain and this was confirmed by the presence of mutations in the major resistance-associated proteins for all anti-tuberculosis (TB) drugs, including isoniazid (KatG (S315T) and InhA (S94A)), fluoroquinolones (S95T), streptomycin (rrs (R309C)), and rifampin (D435Y, a rare but disputed mutation in rpoB). The presence of these mutations exclusively in the human M. bovis isolate suggested that these occurred after transmission from cattle. Genome analysis in this study identified M. bovis in humans and great apes, suggesting possible transmission from domesticated ruminants in the area due to a dynamic and changing interface, which has created opportunity for exposure and transmission.

Whole-genome sequence analysis of the Mycobacterium avium complex and proposal of the transfer of Mycobacterium yongonense to Mycobacterium intracellulare subsp. yongonense subsp. nov

Bacterial whole-genome sequences contain informative features of their evolutionary pathways. Comparison of whole-genome sequences have become the method of choice for classification of prokaryotes, thus allowing the identification of bacteria from an evolutionary perspective, and providing data to resolve some current controversies. Currently, controversy exists about the assignment of members of the Mycobacterium avium complex, as is for the cases of Mycobacterium yongonense and 'Mycobacterium indicus pranii'. These two mycobacteria, closely related to Mycobacterium intracellulare on the basis of standard phenotypic and single gene-sequences comparisons, were not considered a member of such species on the basis on some particular differences displayed by a single strain. Whole-genome sequence comparison procedures, namely the average nucleotide identity and the genome distance, showed that those two mycobacteria should be considered members of the species M. intracellulare. The results were confirmed with other whole-genome comparison supplementary methods. According to the data provided, Mycobacterium yongonense and 'Mycobacterium indicus pranii' should be considered and renamed and included as members of M. intracellulare. This study highlights the problems caused when a novel species is accepted on the basis of a single strain, as was the case for M. yongonense. Based mainly on whole-genome sequence analysis, we conclude that M. yongonense should be reclassified as a subspecies of Mycobacterium intracellulareas Mycobacterium intracellularesubsp. yongonense and 'Mycobacterium indicus pranii' classified in the same subspecies as the type strain of Mycobacterium intracellulare and classified as Mycobacterium intracellularesubsp. intracellulare.

Phylogenomic analysis of the species of the Mycobacterium tuberculosis complex demonstrates that Mycobacterium africanum, Mycobacterium bovis, Mycobacterium caprae, Mycobacterium microti and Mycobacterium pinnipedii are later heterotypic synonyms of Mycobacterium tuberculosis

The species within the Mycobacterium tuberculosis Complex (MTBC) have undergone numerous taxonomic and nomenclatural changes, leaving the true structure of the MTBC in doubt. We used next-generation sequencing (NGS), digital DNA-DNA hybridization (dDDH), and average nucleotide identity (ANI) to investigate the relationship between these species. The type strains of Mycobacterium africanum, Mycobacterium bovis, Mycobacterium caprae, Mycobacterium microti and Mycobacterium pinnipedii were sequenced via NGS. Pairwise dDDH and ANI comparisons between these, previously sequenced MTBC type strain genomes (including 'Mycobacterium canettii', 'Mycobacterium mungi' and 'Mycobacterium orygis') and M. tuberculosis H37Rv^T were performed. Further, all available genome sequences in GenBank for species in or putatively in the MTBC were compared to H37Rv^T. Pairwise results indicated that all of the type strains of the species are extremely closely related to each other (dDDH: 91.2-99.2 %, ANI: 99.21-99.92 %), greatly exceeding the respective species delineation thresholds, thus indicating that they belong to the same species. Results from the GenBank genomes indicate that all the strains examined are within the circumscription of H37Rv^T (dDDH: 83.5-100 %). We, therefore, formally propose a union of the species of the MTBC as M. tuberculosis. M. africanum, M. bovis, M. caprae, M. microti and M. pinnipedii are reclassified as later heterotypic synonyms of M. tuberculosis. 'M. canettii', 'M. mungi', and 'M. orygis' are classified as strains of the species M. tuberculosis. We further recommend use of the infrasubspecific term 'variant' ('var.') and infrasubspecific designations that generally retain the historical nomenclature associated with the groups or otherwise convey such characteristics, e.g. M. tuberculosis var. bovis.

Gene and whole genome analyses reveal that the mycobacterial strain JS623 is not a member of the species Mycobacterium smegmatis

Unexpected differences were found between the genome of strain JS623, used in bioremediation studies, and the genome of strain mc²155, a model organism for investigating basic biology of mycobacteria. Both strains are currently assigned in the databases to the species Mycobacterium smegmatis and, consequently, the environmental isolate JS623 is increasingly included as a representative of that species in comparative genome‐based approaches aiming at identifying distinctive traits of the different members of the genus Mycobacterium. We applied traditional molecular taxonomic procedures – inference of single and concatenated gene trees – to re‐evaluate the membership of both strains to the same species, adopting the latest accepted cut‐off values for species delimitation. Additionally, modern whole genome‐based in silico methods where performed in a comprehensive molecular phylogenetic analysis of JS623 and other members of the genus Mycobacterium. These analyses showed that all relevant genome parameters of JS623 clearly separate this strain from M. smegmatis. The strain JS623 should be corrected as Mycobacterium sp. not only in the literature but, even more importantly, in the database entries, as inclusion of the genome wrongly attributed to the M. smegmatis species in comparative studies will result in misleading conclusions.

PolyTB: a genomic variation map for Mycobacterium tuberculosis

Tuberculosis (TB) caused by Mycobacterium tuberculosis (Mtb) is the second major cause of death from an infectious disease worldwide. Recent advances in DNA sequencing are leading to the ability to generate whole genome information in clinical isolates of M. tuberculosis complex (MTBC). The identification of informative genetic variants such as phylogenetic markers and those associated with drug resistance or virulence will help barcode Mtb in the context of epidemiological, diagnostic and clinical studies. Mtb genomic datasets are increasingly available as raw sequences, which are potentially difficult and computer intensive to process, and compare across studies. Here we have processed the raw sequence data (>1500 isolates, eight studies) to compile a catalogue of SNPs (n = 74,039, 63% non-synonymous, 51.1% in more than one isolate, i.e. non-private), small indels (n = 4810) and larger structural variants (n = 800). We have developed the PolyTB web-based tool (http://pathogenseq.lshtm.ac.uk/polytb) to visualise the resulting variation and important meta-data (e.g. in silico inferred strain-types, location) within geographical map and phylogenetic views. This resource will allow researchers to identify polymorphisms within candidate genes of interest, as well as examine the genomic diversity and distribution of strains. PolyTB source code is freely available to researchers wishing to develop similar tools for their pathogen of interest.

Big data and clinicians: a review on the state of the science

BACKGROUND

In the past few decades, medically related data collection saw a huge increase, referred to as big data. These huge datasets bring challenges in storage, processing, and analysis. In clinical medicine, big data is expected to play an important role in identifying causality of patient symptoms, in predicting hazards of disease incidence or reoccurrence, and in improving primary-care quality.

OBJECTIVE

The objective of this review was to provide an overview of the features of clinical big data, describe a few commonly employed computational algorithms, statistical methods, and software toolkits for data manipulation and analysis, and discuss the challenges and limitations in this realm.

METHODS

We conducted a literature review to identify studies on big data in medicine, especially clinical medicine. We used different combinations of keywords to search PubMed, Science Direct, Web of Knowledge, and Google Scholar for literature of interest from the past 10 years.

RESULTS

This paper reviewed studies that analyzed clinical big data and discussed issues related to storage and analysis of this type of data.

CONCLUSIONS

Big data is becoming a common feature of biological and clinical studies. Researchers who use clinical big data face multiple challenges, and the data itself has limitations. It is imperative that methodologies for data analysis keep pace with our ability to collect and store data.

Unusual large-scale chromosomal rearrangements in Mycobacterium tuberculosis Beijing B0/W148 cluster isolates

The Mycobacterium tuberculosis (MTB) Beijing family isolates are geographically widespread, and there are examples of Beijing isolates that are hypervirulent and associated with drug resistance. One-fourth of Beijing genotype isolates found in Russia belong to the B0/W148 group. The aim of the present study was to investigate features of these endemic strains on a genomic level. Four Russian clinical isolates of this group were sequenced, and the data obtained was compared with published sequences of various MTB strain genomes, including genome of strain W-148 of the same B0/W148 group. The comparison of the W-148 and H37Rv genomes revealed two independent inversions of large segments of the chromosome. The same inversions were found in one of the studied strains after deep sequencing using both the fragment and mate-paired libraries. Additionally, inversions were confirmed by RFLP hybridization analysis. The discovered rearrangements were verified by PCR in all four newly sequenced strains in the study and in four additional strains of the same Beijing B0/W148 group. The other 32 MTB strains from different phylogenetic lineages were tested and revealed no inversions. We suggest that the initial largest inversion changed the orientation of the three megabase (Mb) segment of the chromosome, and the second one occurred in the previously inverted region and partly restored the orientation of the 2.1 Mb inner segment of the region. This is another remarkable example of genomic rearrangements in the MTB in addition to the recently published of large-scale duplications. The described cases suggest that large-scale genomic rearrangements in the currently circulating MTB isolates may occur more frequently than previously considered, and we hope that further studies will help to determine the exact mechanism of such events.

Computational genomics-proteomics and Phylogeny analysis of twenty one mycobacterial genomes (Tuberculosis & non Tuberculosis strains)

Background

The genus Mycobacterium comprises different species, among them the most contagious and infectious bacteria. The members of the complex Mycobacterium tuberculosis are the most virulent microorganisms that have killed human and other mammals since millennia. Additionally, with the many different mycobacterial sequences available, there is a crucial need for the visualization and the simplification of their data. In this present study, we aim to highlight a comparative genome, proteome and phylogeny analysis between twenty-one mycobacterial (Tuberculosis and non tuberculosis) strains using a set of computational and bioinformatics tools (Pan and Core genome plotting, BLAST matrix and phylogeny analysis).

Results

Considerably the result of pan and core genome Plotting demonstrated that less than 1250 Mycobacterium gene families are conserved across all species, and a total set of about 20,000 gene families within the Mycobacterium pan-genome of twenty one mycobacterial genomes.

Viewing the BLAST matrix a high similarity was found among the species of the complex Mycobacterium tuberculosis and less conservation is found with other slow growing pathogenic mycobacteria.

Phylogeny analysis based on both protein conservation, as well as rRNA clearly resolve known relationships between slow growing mycobacteria.

Conclusion

Mycobacteria include important pathogenic species for human and animals and the Mycobacterium tuberculosis complex is the most cause of death of the humankind. The comparative genome analysis could provide a new insight for better controlling and preventing these diseases.