scholarly journals Mycobacterium tuberculosis complex lineage 5 exhibits high levels of within-lineage genomic diversity and differing gene content compared to the type strain H37Rv

2021 ◽  
Vol 7 (7) ◽  
Author(s):  
C. N'Dira Sanoussi ◽  
Mireia Coscolla ◽  
Boatema Ofori-Anyinam ◽  
Isaac Darko Otchere ◽  
Martin Antonio ◽  
...  

Pathogens of the Mycobacterium tuberculosis complex (MTBC) are considered to be monomorphic, with little gene content variation between strains. Nevertheless, several genotypic and phenotypic factors separate strains of the different MTBC lineages (L), especially L5 and L6 (traditionally termed Mycobacterium africanum ) strains, from each other. However, this genome variability and gene content, especially of L5 strains, has not been fully explored and may be important for pathobiology and current approaches for genomic analysis of MTBC strains, including transmission studies. By comparing the genomes of 355 L5 clinical strains (including 3 complete genomes and 352 Illumina whole-genome sequenced isolates) to each other and to H37Rv, we identified multiple genes that were differentially present or absent between H37Rv and L5 strains. Additionally, considerable gene content variability was found across L5 strains, including a split in the L5.3 sub-lineage into L5.3.1 and L5.3.2. These gene content differences had a small knock-on effect on transmission cluster estimation, with clustering rates influenced by the selected reference genome, and with potential overestimation of recent transmission when using H37Rv as the reference genome. We conclude that full capture of the gene diversity, especially high-resolution outbreak analysis, requires a variation of the single H37Rv-centric reference genome mapping approach currently used in most whole-genome sequencing data analysis pipelines. Moreover, the high within-lineage gene content variability suggests that the pan-genome of M. tuberculosis is at least several kilobases larger than previously thought, implying that a concatenated or reference-free genome assembly (de novo) approach may be needed for particular questions.

2020 ◽  
Author(s):  
C. N’Dira Sanoussi ◽  
Mireia Coscolla ◽  
Boatema Ofori-Anyinam ◽  
Isaac Darko Otchere ◽  
Martin Antonio ◽  
...  

AbstractPathogens of the Mycobacterium tuberculosis complex (MTBC) are considered monomorphic, with little gene content variation between strains. Nevertheless, several genotypic and phenotypic factors separate the different MTBC lineages (L), especially L5 and L6 (traditionally termed Mycobacterium africanum), from each other. However, genome variability and gene content especially of L5 and L6 strains have not been fully explored and may be potentially important for pathobiology and current approaches for genomic analysis of MTBC isolates, including transmission studies.We compared the genomes of 358 L5 clinical isolates (including 3 completed genomes and 355 Illumina WGS (whole genome sequenced) isolates) to the L5 complete genomes and H37Rv, and identified multiple genes differentially present or absent between H37Rv and L5 strains. Additionally, considerable gene content variability was found across L5 strains, including a split in the L5.3 sublineage into L5.3.1 and L5.3.2. These gene content differences had a small knock on effect on transmission cluster estimation, with clustering rates influenced by the selection of reference genome, and with potential over-estimation of recent transmission when using H37Rv as the reference genome.Our data show that the use of H37Rv as reference genome results in missing SNPs in genes unique for L5 strains. This potentially leads to an underestimation of the diversity present in the genome of L5 strains and in turn affects the transmission clustering rates. As such, a full capture of the gene diversity, especially for high resolution outbreak analysis, requires a variation of the single H37Rv-centric reference genome mapping approach currently used in most WGS data analysis pipelines. Moreover, the high within-lineage gene content variability suggests that the pan-genome of M. tuberculosis is at least several kilobases larger than previously thought, implying a concatenated or reference-free genome assembly (de novo) approach may be needed for particular questions.Data summarySequence data for the Illumina dataset are available at European Genome-phenome Archive (EGA; https://www.ebi.ac.uk/ega/) under the study accession numbers PRJEB38317 and PRJEB38656. Individual runs accession numbers are indicated in Table S8.PacBio raw reads for the L5 Benin genome are available on the ENA accession SAME3170744. The assembled L5 Benin genome is available on NCBI with accession PRJNA641267. To ensure naming conventions of the genes in the three L5 genomes can be followed, we have uploaded these annotated GFF files to figshare at https://doi.org/10.6084/m9.figshare.12911849.v1.Custom python scripts used in this analysis can be found at https://github.com/conmeehan/pathophy.


2021 ◽  
Vol 7 (7) ◽  
Author(s):  
Casper Jamin ◽  
Sien De Koster ◽  
Stefanie van Koeveringe ◽  
Dieter De Coninck ◽  
Klaas Mensaert ◽  
...  

Whole-genome sequencing (WGS) is becoming the de facto standard for bacterial typing and outbreak surveillance of resistant bacterial pathogens. However, interoperability for WGS of bacterial outbreaks is poorly understood. We hypothesized that harmonization of WGS for outbreak surveillance is achievable through the use of identical protocols for both data generation and data analysis. A set of 30 bacterial isolates, comprising of various species belonging to the Enterobacteriaceae family and Enterococcus genera, were selected and sequenced using the same protocol on the Illumina MiSeq platform in each individual centre. All generated sequencing data were analysed by one centre using BioNumerics (6.7.3) for (i) genotyping origin of replications and antimicrobial resistance genes, (ii) core-genome multi-locus sequence typing (cgMLST) for Escherichia coli and Klebsiella pneumoniae and whole-genome multi-locus sequencing typing (wgMLST) for all species. Additionally, a split k-mer analysis was performed to determine the number of SNPs between samples. A precision of 99.0% and an accuracy of 99.2% was achieved for genotyping. Based on cgMLST, a discrepant allele was called only in 2/27 and 3/15 comparisons between two genomes, for E. coli and K. pneumoniae, respectively. Based on wgMLST, the number of discrepant alleles ranged from 0 to 7 (average 1.6). For SNPs, this ranged from 0 to 11 SNPs (average 3.4). Furthermore, we demonstrate that using different de novo assemblers to analyse the same dataset introduces up to 150 SNPs, which surpasses most thresholds for bacterial outbreaks. This shows the importance of harmonization of data-processing surveillance of bacterial outbreaks. In summary, multi-centre WGS for bacterial surveillance is achievable, but only if protocols are harmonized.


Microbiology ◽  
2021 ◽  
Vol 167 (9) ◽  
Author(s):  
Nirbhay Singh ◽  
Anu Chauhan ◽  
Ram Kumar ◽  
Sudheer Kumar Singh

Branched-chain amino acids (BCAAs) are essential amino acids, but their biosynthetic pathway is absent in mammals. Ketol-acid reductoisomerase (IlvC) is a BCAA biosynthetic enzyme that is coded by Rv3001c in Mycobacterium tuberculosis H37Rv (Mtb-Rv) and MRA_3031 in M. tuberculosis H37Ra (Mtb-Ra). IlvCs are essential in Mtb-Rv as well as in Escherichia coli . Compared to wild-type and IlvC-complemented Mtb-Ra strains, IlvC knockdown strain showed reduced survival at low pH and under low pH+starvation stress conditions. Further, increased expression of IlvC was observed under low pH and starvation stress conditions. Confirmation of a role for IlvC in pH and starvation stress was achieved by developing E. coli BL21(DE3) IlvC knockout, which was defective for growth in M9 minimal medium, but growth could be rescued by isoleucine and valine supplementation. Growth was also restored by complementing with over-expressing constructs of Mtb-Ra and E. coli IlvCs. The E. coli knockout also had a survival deficit at pH=5.5 and 4.5 and was more susceptible to killing at pH=3.0. The biochemical characterization of Mtb-Ra and E. coli IlvCs confirmed that both have NADPH-dependent activity. In conclusion, this study demonstrates the functional complementation of E. coli IlvC by Mtb-Ra IlvC and also suggests that IlvC has a role in tolerance to low pH and starvation stress.


Author(s):  
Hisami Kobayashi ◽  
Yasuhiro Tanizawa ◽  
Mitsuo Sakamoto ◽  
Moriya Ohkuma ◽  
Masanori Tohno

The taxonomic status of the species Clostridium methoxybenzovorans was assessed. The 16S rRNA gene sequence, whole-genome sequence and phenotypic characterizations suggested that the type strain deposited in the American Type Culture Collection ( C. methoxybenzovorans ATCC 700855T) is a member of the species Eubacterium callanderi . Hence, C. methoxybenzovorans ATCC 700855T cannot be used as a reference for taxonomic study. The type strain deposited in the German Collection of Microorganism and Cell Cultures GmbH (DSM 12182T) is no longer listed in its online catalogue. Also, both the 16S rRNA gene and the whole-genome sequences of the original strain SR3T showed high sequence identity with those of Lacrimispora indolis (recently reclassified from Clostridium indolis ) as the most closely related species. Analysis of the two genomes showed average nucleotide identity based on blast and digital DNA–DNA hybridization values of 98.3 and 87.9 %, respectively. Based on these results, C. methoxybenzovorans SR3T was considered to be a member of L. indolis .


2013 ◽  
Vol 63 (Pt_10) ◽  
pp. 3920-3926 ◽  
Author(s):  
Julia S. Bennett ◽  
Keith A. Jolley ◽  
Martin C. J. Maiden

Phylogenies generated from whole genome sequence (WGS) data provide definitive means of bacterial isolate characterization for typing and taxonomy. The species status of strains recently defined with conventional taxonomic approaches as representing Neisseria oralis was examined by the analysis of sequences derived from WGS data, specifically: (i) 53 Neisseria ribosomal protein subunit (rps) genes (ribosomal multi-locus sequence typing, rMLST); and (ii) 246 Neisseria core genes (core genome MLST, cgMLST). These data were compared with phylogenies derived from 16S and 23S rRNA gene sequences, demonstrating that the N. oralis strains were monophyletic with strains described previously as representing ‘ Neisseria mucosa var. heidelbergensis’ and that this group was of equivalent taxonomic status to other well-described species of the genus Neisseria . Phylogenetic analyses also indicated that Neisseria sicca and Neisseria macacae should be considered the same species as Neisseria mucosa and that Neisseria flavescens should be considered the same species as Neisseria subflava . Analyses using rMLST showed that some strains currently defined as belonging to the genus Neisseria were more closely related to species belonging to other genera within the family; however, whole genome analysis of a more comprehensive selection of strains from within the family Neisseriaceae would be necessary to confirm this. We suggest that strains previously identified as representing ‘ N. mucosa var. heidelbergensis’ and deposited in culture collections should be renamed N. oralis . Finally, one of the strains of N. oralis was able to ferment lactose, due to the presence of β-galactosidase and lactose permease genes, a characteristic previously thought to be unique to Neisseria lactamica , which therefore cannot be thought of as diagnostic for this species; however, the rMLST and cgMLST analyses confirm that N. oralis is most closely related to N. mucosa .


2020 ◽  
Vol 70 (4) ◽  
pp. 2602-2610 ◽  
Author(s):  
Shih-Yao Lin ◽  
Wen-Ming Chen ◽  
Asif Hameed ◽  
Guan-Hua Huang ◽  
Mei-Hua Hung ◽  
...  

A novel Gram-stain-positive, aerobic, motile with peritrichous flagella, rod-shaped bacterium, designated CC-MHH1044T, was isolated from a preserved vegetable sample. A polyphasic taxonomic approach was applied to the isolates in order to clarify its taxonomic position. Growth of the strain CC-MHH1044T occurred at 15–50 °C (optimum, 30 °C), pH 6.0–8.0 (optimum, pH 7.0) and with 0–2.0 % (w/v) NaCl (optimum, 1 %, w/v). The genome of strain CC-MHH1044T consisted of 8.5 Mb and the genomic DNA G+C content was 58.5 mol%. Comparison of the 16S rRNA gene sequences showed that CC-MHH1044T belonged to the genus Cohnella and showed a close relationship with the type strains of Cohnella damuensis (96.2 %) and Cohnella panacarvi (95.9 %), and lower sequence similarity to other species. Average nucleotide identity values calculated from whole-genome sequencing data proved that CC-MHH1044T represents a distinct Cohnella species. The dominant cellular fatty acids (>5 %) included iso-C14 : 0(7.4 %), iso-C15 : 0 (6.4 %), anteiso-C15 : 0(40.3 %), C16 : 0 (6.6 %) and iso-C16 : 0 (27.0 %). The polar lipid profile consisted of diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, four unidentified aminophospholipids, one unidentified phospholipid and glycolipid. The major polyamine was spermidine. The predominant isoprenoid quinone was menaqinone 7 (MK-7). Based on its distinct phylogenetic, phenotypic and chemotaxonomic traits, together with results of comparative 16S rRNA gene sequence, average nucleotide identity and digital DNA–DNA hybridization analyses, we conclude that strain CC-MHH1044T represents a novel member of the genus Cohnella , for which the name Cohnella fermenti sp. nov. is proposed. The type strain is CC-MHH1044T (=BCRC 81147T=JCM 32834T).


Microbiology ◽  
2021 ◽  
Author(s):  
Céline Rens ◽  
Joseph D. Chao ◽  
Danielle L. Sexton ◽  
Elitza I. Tocheva ◽  
Yossef Av-Gay

The success of Mycobacterium tuberculosis as a pathogen is well established: tuberculosis is the leading cause of death by a single infectious agent worldwide. The threat of multi- and extensively drug-resistant bacteria has renewed global concerns about this pathogen and understanding its virulence strategies will be essential in the fight against tuberculosis. The current review will focus on phthiocerol dimycocerosates (PDIMs), a long-known and well-studied group of complex lipids found in the M. tuberculosis cell envelope. Numerous studies show a role for PDIMs in several key steps of M. tuberculosis pathogenesis, with recent studies highlighting its involvement in bacterial virulence, in association with the ESX-1 secretion system. Yet, the mechanisms by which PDIMs help M. tuberculosis to control macrophage phagocytosis, inhibit phagosome acidification and modulate host innate immunity, remain to be fully elucidated.


2021 ◽  
Vol 7 (5) ◽  
Author(s):  
Antoni P. A. Hendrickx ◽  
Fabian Landman ◽  
Angela de Haan ◽  
Sandra Witteveen ◽  
Marga G. van Santen-Verheuvel ◽  
...  

Carbapenem-hydrolysing enzymes belonging to the OXA-48-like group are encoded by bla OXA-48-like alleles and are abundant among Enterobacterales in the Netherlands. Therefore, the objective here was to investigate the characteristics, gene content and diversity of the bla OXA-48-like carrying plasmids and chromosomes of Escherichia coli and Klebsiella pneumoniae collected in the Dutch national surveillance from 2014 to 2019 in comparison with genome sequences from 29 countries. A combination of short-read genome sequencing with long-read sequencing enabled the reconstruction of 47 and 132 complete bla OXA-48-like plasmids for E. coli and K. pneumoniae , respectively. Seven distinct plasmid groups designated as pOXA-48-1 to pOXA-48-5, pOXA-181 and pOXA-232 were identified in the Netherlands which were similar to internationally reported plasmids obtained from countries from North and South America, Europe, Asia and Oceania. The seven plasmid groups varied in size, G+C content, presence of antibiotic resistance genes, replicon family and gene content. The pOXA-48-1 to pOXA-48-5 plasmids were variable, and the pOXA-181 and pOXA-232 plasmids were conserved. The pOXA-48-1, pOXA-48-2, pOXA-48-3 and pOXA-48-5 groups contained a putative conjugation system, but this was absent in the pOXA-48-4, pOXA-181 and pOXA-232 plasmid groups. pOXA-48 plasmids contained the PemI antitoxin, while the pOXA-181 and pOXA-232 plasmids did not. Furthermore, the pOXA-181 plasmids carried a virB2-virB3-virB9-virB10-virB11 type IV secretion system, while the pOXA-48 plasmids and pOXA-232 lacked this system. A group of non-related pOXA-48 plasmids from the Netherlands contained different resistance genes, non-IncL-type replicons or no replicons. Whole genome multilocus sequence typing revealed that the bla OXA-48-like plasmids were found in a wide variety of genetic backgrounds in contrast to chromosomally encoded bla OXA-48-like alleles. Chromosomally localized bla OXA-48 and bla OXA-244 alleles were located on genetic elements of variable sizes and comprised regions of pOXA-48 plasmids. The bla OXA-48-like genetic element was flanked by a direct repeat upstream of IS1R, and was found at multiple locations in the chromosomes of E. coli . Lastly, K. pneumoniae isolates carrying bla OXA-48 or bla OXA-232 were mostly resistant for meropenem, whereas E. coli bla OXA-48, bla OXA-181 and chromosomal bla OXA-48 or bla OXA-244 isolates were mostly sensitive. In conclusion, the overall bla OXA-48-like plasmid population in the Netherlands is conserved and similar to that reported for other countries, confirming global dissemination of bla OXA-48-like plasmids. Variations in size, presence of antibiotic resistance genes and gene content impacted pOXA-48, pOXA-181 and pOXA-232 plasmid architecture.


2021 ◽  
Vol 7 (11) ◽  
Author(s):  
Isabelle Bernaquez ◽  
Christiane Gaudreau ◽  
Pierre A. Pilon ◽  
Sadjia Bekal

Many public health laboratories across the world have implemented whole-genome sequencing (WGS) for the surveillance and outbreak detection of foodborne pathogens. PulseNet-affiliated laboratories have determined that most single-strain foodborne outbreaks are contained within 0–10 multi-locus sequence typing (MLST)-based allele differences and/or core genome single-nucleotide variants (SNVs). In addition to being a food- and travel-associated outbreak pathogen, most Shigella spp. cases occur through continuous person-to-person transmission, predominantly involving men who have sex with men (MSM), leading to long-term and recurrent outbreaks. Continuous transmission patterns coupled to genetic evolution under antibiotic treatment pressure require an assessment of existing WGS-based subtyping methods and interpretation criteria for cluster inclusion/exclusion. An evaluation of 4 WGS-based subtyping methods [SNVPhyl, coreMLST, core genome MLST (cgMLST) and whole-genome MLST (wgMLST)] was performed on 9 foodborne-, travel- and MSM-related retrospective outbreaks from a collection of 91 Shigella flexneri and 232  Shigella sonnei isolates to determine the methods’ epidemiological concordance, discriminatory power, robustness and ability to generate stable interpretation criteria. The discriminatory powers were ranked as follows: coreMLST<SNVPhyl<cgMLST<wgMLST (range: 0.970–1.000). The genetic differences observed for non-MSM-related Shigella spp. outbreaks respect the standard 0–10 allele/SNV guideline; however, mobile genetic element (MGE)-encoded loci caused inflated genetic variation and discrepant phylogenies for prolonged MSM-related S. sonnei outbreaks via wgMLST. The S. sonnei correlation coefficients of wgMLST were also the lowest at 0.680, 0.703 and 0.712 for SNVPhyl, coreMLST and cgMLST, respectively. Plasmid maintenance, mobilization and conjugation-associated genes were found to be the main source of genetic distance inflation in addition to prophage-related genes. Duplicated alleles arising from the repeated nature of IS elements were also responsible for many false cg/wgMLST differences. The coreMLST approach was shown to be the most robust, followed by SNVPhyl and wgMLST for inter-laboratory comparability. Our results highlight the need for validating species-specific subtyping methods based on microbial genome plasticity and outbreak dynamics in addition to the importance of filtering confounding MGEs for cluster detection.


2019 ◽  
Vol 5 (7) ◽  
Author(s):  
Charles H. D. Williamson ◽  
Nathan E. Stone ◽  
Amalee E. Nunnally ◽  
Heidie M. Hornstra ◽  
David M. Wagner ◽  
...  

Clostridioides difficile is a ubiquitous, diarrhoeagenic pathogen often associated with healthcare-acquired infections that can cause a range of symptoms from mild, self-limiting disease to toxic megacolon and death. Since the early 2000s, a large proportion of C. difficile cases have been attributed to the ribotype 027 (RT027) lineage, which is associated with sequence type 1 (ST1) in the C. difficile multilocus sequence typing scheme. The spread of ST1 has been attributed, in part, to resistance to fluoroquinolones used to treat unrelated infections, which creates conditions ideal for C. difficile colonization and proliferation. In this study, we analysed 27 isolates from a healthcare network in northern Arizona, USA, and 1352 publicly available ST1 genomes to place locally sampled isolates into a global context. Whole genome, single nucleotide polymorphism analysis demonstrated that at least six separate introductions of ST1 were observed in healthcare facilities in northern Arizona over an 18-month sampling period. A reconstruction of transmission networks identified potential nosocomial transmission of isolates, which were only identified via whole genome sequence analysis. Antibiotic resistance heterogeneity was observed among ST1 genomes, including variability in resistance profiles among locally sampled ST1 isolates. To investigate why ST1 genomes are so common globally and in northern Arizona, we compared all high-quality C. difficile genomes and identified that ST1 genomes have gained and lost a number of genomic regions compared to all other C. difficile genomes; analyses of other toxigenic C. difficile sequence types demonstrate that this loss may be anomalous and could be related to niche specialization. These results suggest that a combination of antimicrobial resistance and gain and loss of specific genes may explain the prominent association of this sequence type with C. difficile infection cases worldwide. The degree of genetic variability in ST1 suggests that classifying all ST1 genomes into a quinolone-resistant hypervirulent clone category may not be appropriate. Whole genome sequencing of clinical C. difficile isolates provides a high-resolution surveillance strategy for monitoring persistence and transmission of C. difficile and for assessing the performance of infection prevention and control strategies.


Sign in / Sign up

Export Citation Format

Share Document