Towards the reconstruction of a global TB history using a new pipeline "TB-Annotator".

Mycobacterium tuberculosis complex (MTBC) has a population structure consisting of 9 human and animal lineages. The genomic diversity within these lineages is a pathogenesis factor that affects virulence, transmissibility, host response, and antibiotic resistance. Hence it is important to develop improved information systems for tracking and understanding the spreading and evolution of genomes. We present results obtained thanks to a new informatics platform for computational biology of MTBC, that uses a convenience sample from public/private SRAs, designated as TB-Annotator. Version 1 was a first interactive graphic-based web tool based on 15,901 representative genomes. Version 2, still interactive, is a more sophisticated database, developed using the Snakemake Workflow Management System (WMS) that allows an unsupervised global and scalable analysis of the content of the USA National Center for Biotechnology Information Short Read Archives database. This platform analyzes nucleotide variants, the presence/absence of genes, known regions of difference and detect new deletions, the insertion sites of mobile genetic elements, and allows phylogenetic trees to be built, imported in a graphical interface and interactively analyzed between the data and the tree. The objective of TB-Annotator is triple: detect recent epidemiological links, reconstruct distant phylogeographical histories as well as perform more complex phenotypic/genotypic Genome-Wide Association Studies (GWAS). In this paper, we compare the various taxonomic SNPs-based labels and hierarchies previously described in recent reference papers for L1, and present a comparative analysis that allows identification of alias and thus provides the basis of a future unifying naming scheme for L1 sublineages. We present a global phylogenetic tree built with RAxML-NG, and one on L2; at the time of writing, we characterized about 200 sublineages, with many new ones; a detail tree for Modern L2 and a hierarchical scheme allowing to facilitate L2 lineage assignment are also presented.

Mycobacterium tuberculosis complex drug-resistance, phylogenetics, and evolution in Nigeria: Comparison with Ghana and Cameroon.

In this article, we provide an in-depth analysis on the drug-resistance phenotypic characteristics of a cohort of 325 tuberculosis and characterize by Whole Genome Sequencing 24 isolates from Nigeria belonging to L4, L5 and L6. Our results suggest an alarming rate of drug-resistance of the L4.6.2.2 Mycobacterium tuberculosis complex (MTBC) lineage and a high diversity of L5. We compiled these new Sequence Read Archives (SRAs) to previously published ones from available Bioprojects run in Nigeria. We performed RAxML phylogenetic reconstructions of larger samples that include public NCBI SRAs from some neighboring countries (Cameroon, Ghana). To confront phylogenetic reconstruction to metadata, we used a new proprietary database named TB-Annotator. We show that L5 genomes in Northern Nigeria belong to new clades as the ones described until now and allow an update of the taxonomy of L5. In addition, we describe the L4.6.2.2 lineage in Nigeria, Cameroon and Ghana. We provide computations on the likely divergence time of L4.6.2.2 and suggest a new hypothesis concerning its origin. Finally we provide a short overview on M. bovis diversity in Nigeria. This study constitutes a baseline knowledge on the global genomic diversity, phylogeography and phylodynamics of MTBC in Nigeria, as well as on the natural history of this largely ignored but densely populated country of Africa. These results highlight the need of sequencing additional MTBC genomes in Nigeria and more generally in West-Africa, both for public health and for academic reasons. The likelihood of replacement of L5-L6 by L4.6.2.2 isolates, leave potentially little time to gather historical knowledge informative on the ancient history of tuberculosis in West-Africa.

Hi-plex deep amplicon sequencing for identification, high-resolution genotyping and multidrug resistance prediction of Mycobacterium leprae directly from patient biopsies by using Deeplex Myc-Lep.

BACKGROUND: Expansion of antimicrobial resistance monitoring and epidemiological surveillance are key components of the WHO strategy towards zero leprosy. The inability to grow Mycobacterium leprae in vitro precludes routine phenotypic drug susceptibility testing, and only limited molecular tests are available. We evaluated a culture-free targeted deep sequencing assay, for mycobacterial identification, genotyping based on 18 canonical SNPs and 11 core variable-number tandem-repeat (VNTR) markers, and detection of rifampicin, dapsone and fluoroquinolone resistance-associated mutations in rpoB/ctpC/ctpI, folP1, gyrA/gyrB, respectively, and hypermutation-associated mutations in nth. METHODS: The limit of detection (LOD) was determined using DNA of M. leprae reference strains and from 246 skin biopsies and 74 slit skin smears of leprosy patients, with genome copies quantified by RLEP qPCR. Sequencing results were evaluated versus whole genome sequencing (WGS) data of 14 strains, and versus VNTR-fragment length analysis (FLA) results of 89 clinical specimens. FINDINGS: The LOD for sequencing success ranged between 80 and 3000 genome copies, depending on the sample type. The LOD for minority variants was 10%. All SNPs detected in targets by WGS were identified except in a clinical sample where WGS revealed two dapsone resistance-conferring mutations instead of one by Deeplex Myc-Lep, due to partial duplication of the sulfamide-binding domain in folP1. SNPs detected uniquely by Deeplex Myc-Lep were missed by WGS due to insufficient coverage. Concordance with VNTR-FLA results was 99.4% (926/932 alleles). INTERPRETATION: Deeplex Myc-Lep may help improve the diagnosis and surveillance of leprosy. Gene domain duplication is an original putative drug resistance-related genetic adaptation in M. leprae. FUNDING: EDCTP2 programme supported by the European Union (grant number RIA2017NIM-1847 -PEOPLE). EDCTP, R2Stop: Effect:Hope, The Mission To End Leprosy, the Flemish Fonds Wetenschappelijk Onderzoek.

Mycobacterium chimaera Genomics With Regard to Epidemiological and Clinical Investigations Conducted for an Open Chest Postsurgical Mycobacterium chimaera Infection Outbreak.

BACKGROUND: Postsurgical infections due to Mycobacterium chimaera appeared as a novel nosocomial threat in 2015, with a worldwide outbreak due to contaminated heater-cooler units used in open chest surgery. We report the results of investigations conducted in France including whole-genome sequencing comparison of patient and heater-cooler unit isolates. METHODS: We sought M. chimaera infection cases from 2010 onwards through national epidemiological investigations in health care facilities performing cardiopulmonary bypass, together with a survey on good practices and systematic heater-cooler unit microbial analyses. Clinical and heater-cooler unit isolates were subjected to whole-genome sequencing analyzed with regard to the reference outbreak strain Zuerich-1. RESULTS: Only 2 clinical cases were shown to be related to the outbreak, although 23% (41/175) of heater-cooler units were declared positive for M. avium complex. Specific measures to prevent infection were applied in 89% (50/56) of health care facilities, although only 14% (8/56) of them followed the manufacturer maintenance recommendations. Whole-genome sequencing comparison showed that the clinical isolates and 72% (26/36) of heater-cooler unit isolates belonged to the epidemic cluster. Within clinical isolates, 5-9 nonsynonymous single nucleotide polymorphisms were observed, among which an in vivo mutation in a putative efflux pump gene was observed in a clinical isolate obtained for 1 patient on antimicrobial treatment. CONCLUSIONS: Cases of postsurgical M. chimaera infections have been declared to be rare in France, although heater-cooler units were contaminated, as in other countries. Genomic analyses confirmed the connection to the outbreak and identified specific single nucleotide polymorphisms, including 1 suggesting fitness evolution in vivo.

Genetic and metabolic signatures of Salmonella enterica subsp. enterica associated with animal sources at the pangenomic scale.

BACKGROUND: Salmonella enterica subsp. enterica is a public health issue related to food safety, and its adaptation to animal sources remains poorly described at the pangenome scale. Firstly, serovars presenting potential mono- and multi-animal sources were selected from a curated and synthetized subset of Enterobase. The corresponding sequencing reads were downloaded from the European Nucleotide Archive (ENA) providing a balanced dataset of 440 Salmonella genomes in terms of serovars and sources (i). Secondly, the coregenome variants and accessory genes were detected (ii). Thirdly, single nucleotide polymorphisms and small insertions/deletions from the coregenome, as well as the accessory genes were associated to animal sources based on a microbial Genome Wide Association Study (GWAS) integrating an advanced correction of the population structure (iii). Lastly, a Gene Ontology Enrichment Analysis (GOEA) was applied to emphasize metabolic pathways mainly impacted by the pangenomic mutations associated to animal sources (iv). RESULTS: Based on a genome dataset including Salmonella serovars from mono- and multi-animal sources (i), 19,130 accessory genes and 178,351 coregenome variants were identified (ii). Among these pangenomic mutations, 52 genomic signatures (iii) and 9 over-enriched metabolic signatures (iv) were associated to avian, bovine, swine and fish sources by GWAS and GOEA, respectively. CONCLUSIONS: Our results suggest that the genetic and metabolic determinants of Salmonella adaptation to animal sources may have been driven by the natural feeding environment of the animal, distinct livestock diets modified by human, environmental stimuli, physiological properties of the animal itself, and work habits for health protection of livestock.