Analysis of IS6110 insertion sites provide a glimpse into genome evolution of Mycobacterium tuberculosis

Detection of Mycobacterium tuberculosis from tongue swabs using sonication and sequence-specific hybridization capture

Tongue swabs hold promise as a non-invasive sample for diagnosing tuberculosis (TB). However, their utility as replacements for sputum has been limited by their varied diagnostic performance in PCR assays compared to sputum. The use of silica-based DNA extraction methods may limit sensitivity due to incomplete lysis of Mycobacterium tuberculosis (MTB) cells and co-extraction of non-target nucleic acid, which may inhibit PCR. Specificity may also be compromised because these methods are labor-intensive and prone to cross-contamination. To address these limitations, we developed a sample preparation method that combines sonication for MTB lysis and a sequence-specific MTB DNA capture method using hybridization probes immobilized on magnetic beads. In spiked tongue swabs, our hybridization capture method demonstrated a 100-fold increase in MTB DNA yield over silica-based Qiagen DNA extraction and ethanol precipitation. In a study conducted on clinical samples from South Africa, our protocol had 74% (70/94) sensitivity and 98% (41/42) specificity for detecting active pulmonary TB with sputum Xpert MTB/RIF Ultra as the reference standard. While hybridization capture did not show improved sensitivity over Qiagen DNA extraction and ethanol precipitation, it demonstrated better specificity than previously reported methods and was easier to perform. With integration into point-of-care platforms, these strategies have the potential to help enable rapid non-sputum-based TB diagnosis across key underserved patient populations.

Detection of Mycobacterium tuberculosis complex in saliva by quantitative PCR: A potential alternative specimen for pulmonary tuberculosis diagnosis

BACKGROUND

Current tuberculosis (TB) diagnostic tests primarily rely on sputum samples, yet many TB patients cannot produce sputum. This study explored whether saliva could be used instead of sputum to diagnose pulmonary TB (PTB).

METHOD

The study included 32 patients with confirmed PTB and 30 patients with other respiratory diseases (ORD). Saliva from all study participants was subjected to quantitative (qPCR) assays targeting the IS1081 gene for detection of M. tuberculosis complex DNA.

RESULTS

The sensitivity of saliva IS1081 qPCR was 65.6 % (95 % CI 48.4-80.2 %) with positive results for 21/32 PTB cases, while the specificity was 96.7 % (95 % CI 85.9-99.6 %) with negative results for 29/30 participants with ORD. Sensitivity improved to 72.4 % (95 % CI 54.6-86.0 %) when sputum-Xpert was used as the reference standard, while remaining similar at 65.5 % (95 % CI 47.4-80.7 %) when culture was used as the reference standard. In receiver operating characteristic (ROC) curve analysis, the area under the curve (AUC) for saliva IS1081 qPCR was 82.5 % (95 % CI 71.7-93.3 %).

CONCLUSION

Saliva testing offers a promising alternative to sputum for TB diagnosis among confirmed PTB cases. Larger multicenter studies, encompassing diverse clinical TB characteristics, are needed to provide improved estimates of diagnostic sensitivity and specificity.

Delineating the Acquired Genetic Diversity and Multidrug Resistance in Alcaligenes from Poultry Farms and Nearby Soil

Alcaligenes faecalis is one of the most important and clinically significant environmental pathogens, increasing in importance due to its isolation from soil and nosocomial environments. The Gram-negative soil bacterium is associated with skin endocarditis, bacteremia, dysentery, meningitis, endophthalmitis, urinary tract infections, and pneumonia in patients. With emerging antibiotic resistance in A. faecalis, it has become crucial to understand the origin of such resistance genes within this clinically significant environmental and gut bacterium. In this research, we studied the impact of antibiotic overuse in poultry and its effect on developing resistance in A. faecalis. We sampled soil and faecal materials from five poultry farms, performed whole genome sequencing & analysis and identified four strains of A. faecalis. Furthermore, we characterized the genes in the genomic islands of A. faecalis isolates. We found four multidrug-resistant A. faecalis strains that showed resistance against vancomycin (MIC >1000 μg/ml), ceftazidime (50 μg/ml), colistin (50 μg/ml) and ciprofloxacin (50 μg/ml). From whole genome comparative analysis, we found more than 180 resistance genes compared to the reference sequence. Parts of our assembled contigs were found to be similar to different bacteria which included pbp1A and pbp2 imparting resistance to amoxicillin originally a part of Helicobacter and Bordetella pertussis. We also found the Mycobacterial insertion element IS6110 in the genomic islands of all four genomes. This prominent insertion element can be transferred and induce resistance to other bacterial genomes. The results thus are crucial in understanding the transfer of resistance genes in the environment and can help in developing regimes for antibiotic use in the food and poultry industry.

Diagnostic accuracy of a sequence-specific Mtb-DNA hybridization assay in urine: a case-control study including subclinical TB cases

Tuberculosis (TB) caused by Mycobacterium tuberculosis (Mtb) remains one of the deadliest infectious diseases globally. Timely diagnosis is a key step in the management of TB patients and in the prevention of further transmission events. Current diagnostic tools are limited in these regards. There is an urgent need for new accurate non-sputum-based diagnostic tools for the detection of symptomatic as well as subclinical TB. In this study, we recruited 52 symptomatic TB patients (sputum Xpert MTB/RIF positive) and 58 household contacts to assess the accuracy of a sequence-specific hybridization assay that detects the presence of Mtb cell-free DNA in urine. Using sputum Xpert MTB/RIF as a reference test, the magnetic bead-capture assay could discriminate active TB from healthy household contacts with an overall sensitivity of 72.1% [confidence interval (CI) 0.59-0.86] and specificity of 95.5% (CI 0.90-1.02) with a positive predictive value of 93.9% and negative predictive value of 78.2%. The detection of Mtb-specific DNA in urine suggested four asymptomatic TB infection cases that were confirmed in all instances either by concomitant Xpert MTB/RIF sputum testing or by follow-up investigation raising the specificity of the index test to 100%. We conclude that sequence-specific hybridization assays on urine specimens hold promise as non-invasive tests for the detection of subclinical TB.

IMPORTANCE

There is an urgent need for a non-sputum-based diagnostic tool allowing sensitive and specific detection of all forms of tuberculosis (TB) infections. In that context, we performed a case-control study to assess the accuracy of a molecular detection method enabling the identification of cell-free DNA from Mycobacterium tuberculosis that is shed in the urine of tuberculosis patients. We present accuracy data that would fulfill the target product profile for a non-sputum test. In addition, recent epidemiological data suggested that up to 50% of individuals secreting live bacilli do not present with symptoms at the time of screening. We report, here, that the investigated index test could also detect instances of asymptomatic TB infections among household contacts.

Genomic Characterization of IS6110 Insertions in Mycobacterium orygis

Mycobacterium orygis, a subspecies of the Mycobacterium tuberculosis complex (MTBC), has emerged as a significant concern in the context of One Health, with implications for zoonosis or zooanthroponosis or both. MTBC strains are characterized by the unique insertion element IS6110, which is widely used as a diagnostic marker. IS6110 transposition drives genetic modifications in MTBC, imparting genome plasticity and profound biological consequences. While IS6110 insertions are customarily found in the MTBC genomes, the evolutionary trajectory of strains seems to correlate with the number of IS6110 copies, indicating enhanced adaptability with increasing copy numbers. Here, we present a comprehensive analysis of IS6110 insertions in the M. orygis genome, utilizing ISMapper, and elucidate their genetic consequences in promoting successful host adaptation. Our study encompasses a panel of 67 paired-end reads, comprising 11 isolates from our laboratory and 56 sequences downloaded from public databases. Among these sequences, 91% exhibited high-copy, 4.5% low-copy, and 4.5% lacked IS6110 insertions. We identified 255 insertion loci, including 141 intragenic and 114 intergenic insertions. Most of these loci were either unique or shared among a limited number of isolates, potentially influencing strain behavior. Furthermore, we conducted gene ontology and pathway analysis, using eggNOG-mapper 5.0, on the protein sequences disrupted by IS6110 insertions, revealing 63 genes involved in diverse functions of Gene Ontology and 45 genes participating in various KEGG pathways. Our findings offer novel insights into IS6110 insertions, their preferential insertion regions, and their impact on metabolic processes and pathways, providing valuable knowledge on the genetic changes underpinning IS6110 transposition in M. orygis.

Direct TAMRA-dUTP labeling of M. tuberculosis genes using loop-mediated isothermal amplification (LAMP)

Fluorescent molecule-based direct labeling of amplified DNA is a sensitive method employed across diverse DNA detection and diagnostics systems. However, using pre-labeled primers only allows for the attachment of a single fluorophore to each DNA strand and any modifications of the system are less flexible, requiring new sets of primers. As an alternative, direct labeling of amplified products with modified nucleotides is available, but still poorly characterized. To address these limitations, we sought a direct and adaptable approach to label amplicons produced through Loop-mediated isothermal amplification (LAMP), using labeled nucleotides (dUTPs) rather than primers. The focus of this study was the development and examination of a direct labeling technique of specific genes, including those associated with drug resistance in Mycobacterium tuberculosis. We used 5-(3-Aminoallyl)-2'-deoxyuridine-5'triphosphate, tagged with 5/6-TAMRA (TAMRA-dUTP) for labeling LAMP amplicons during the amplification process and characterized amplification and incorporation efficiency. The optimal TAMRA-dUTP concentration was first determined based on amplification efficiency (0.5% to total dNTPs). Higher concentrations of modified nucleotides reduced or completely inhibited the amplification yield. Target size also showed to be determinant to the success of amplification, as longer sequences showed lower amplification rates, thus less TAMRA incorporated amplicons. Finally, we were able to successfully amplify all four M. tuberculosis target genes using LAMP and TAMRA-modified dUTPs.

Accuracy of Xpert® MTB/RIF Ultra test for posterior oropharyngeal saliva for the diagnosis of paucibacillary pulmonary tuberculosis: a prospective multicenter study

BACKGROUND

Posterior oropharyngeal saliva (POS) is increasingly recognized as an alternative specimen for detecting respiratory pathogens. The accuracy of Xpert® MTB/RIF Ultra (X-Ultra), when performed on POS obtained from patients with paucibacillary pulmonary tuberculosis (TB) is unclear.

METHODS

We consecutively recruited adults with symptoms suggestive of pulmonary TB who were negative by both smear microscopy and Xpert MTB/RIF (X-Classic). Each participant was required to provide one bronchoalveolar lavage fluid (BALF) and one POS specimen, respectively. Diagnostic performances of X-Ultra and X-Classic on POS were compared against clinical and mycobacterial reference standards.

FINDINGS

686 participants meeting inclusion criteria were consecutively enrolled into the study. The overall diagnostic sensitivities of X-Ultra and X-Classic on POS samples were 78.9% [95% confidence interval (CI): 72.8-83.8] and 56.4% (95% CI: 49.7-62.9), respectively; the specificities were 96.6% (95% CI: 94.3-98.1) for X-Ultra and 97.6 (95CI: 95.5-98.8) for X-Classic in POS specimens. Notably, the sensitivity of X-Ultra on POS was as sensitive as X-Classic on BALF against microbiological reference standard (78.9% VS 73.1%). Against clinical diagnosis as a reference standard, the sensitivities of X-Ultra and X-Classic on POS were 55.9% (95% CI: 50.5-61.2; 193/345) and 40.0% (95% CI: 34.8-45.4; 138/345), respectively. The risk of negative results with POS was dramatically increased with decreasing bacterial loads.

CONCLUSIONS

The testing of POS using X-Ultra shows promise as a tool to identify patients with paucibacillary TB. Considering that bronchoscopy is a semi-invasive procedure, POS testing ahead of bronchoscopy, may decrease the need for bronchoscopic procedures, and the cost of care.

Spoligotyping of Mycobacterium tuberculosis - Comparing in vitro and in silico approaches

Spoligotyping is one of the molecular typing methods widely used for exploring the genetic variety of Mycobacterium tuberculosis. The aim of this study was to compare the spoligoprofiles of M. tuberculosis clinical isolates, obtained using in vitro and in silico approaches. The study included 230 M. tuberculosis isolates, recovered from Poland and Lithuania between 2018 and 2021. Spoligotyping in vitro was performed with a commercially available kit. Whole genome sequencing (WGS) was done with Illumina NovaSeq 6000 sequencer. Spoligotype International Types (SITs) were assigned according to the SITVIT2 database or using three different in silico tools, and based on WGS data, namely SpoTyping, SpolPred, and lorikeet. Upon in vitro spoligotyping, the isolates produced 65 different spoligotypes. Spoligotypes inferred from the WGS data were congruent with in vitro generated patterns in 81.7% (188/230) for lorikeet and 81.3% (187/230) for SpolPred and SpoTyping. Spacers 18 and 31 produced the highest ratio of discrepant results between in vitro and in silico approaches, with their signals discordantly assigned for 15 (6.5%) and 9 (3.9%) isolates, respectively. All three in silico approaches used were similarly efficient for M. tuberculosis spoligotype prediction. However, only SpoTyping could predict spoligotypes without a need for manual curation. Thus, we consider it as the most accurate tool. Its use is further advocated by the shortest time of analysis. A relatively high (ca. 20%) discordance between in vitro and in silico spoligotyping results was observed. While we discourage comparing conventional spoligotyping with in silico equivalents, we advise the use of the latter, as it improves the accuracy of spoligopatterns, and thus depicts the relatedness between the isolates more reliably.

Recent progress in the genotyping of bovine tuberculosis and its rapid diagnosis via nanoparticle-based electrochemical biosensors

Bovine tuberculosis (bTB) is considered a worldwide infectious zoonotic disease. Mycobacterium bovis causes bTB disease. It is one of the Mycobacterium tuberculosis complex (MTBC) members. MTBC is a clonal complex of close relatives with approximately 99.95% similarity. M. bovis is a spillover pathogen that can transmit from animals to humans and rarely from humans to animals with contact. Genotyping techniques are important to discriminate and differentiate between MTBC species. Spoligotyping and mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) are widely used but they have some limitations. As an alternative, whole genome sequencing approaches have been utilized due to their high-resolution power. They are employed in typing M. bovis and explain the evolutionary and phylogenetic relationships between isolates. The control of bTB disease has attracted a large amount of attention. Rapid and proper diagnosis is necessary for monitoring the disease as an initial step for its control and treatment. Nanotechnology has a potential impact on the rapid diagnosis and treatment of bTB through the use of nanocarrier and metal nanoparticles (NPs). Special attention has been paid to voltammetric and impedimetric electrochemical strategies as facile, sensitive, and selective methods for the efficient detection of tuberculosis. The efficacy of these sensors is enhanced in the presence of NPs, which act as recognition and/or redox probes. Gold, silver, copper, cobalt, graphene, and magnetic NPs, as well as polypyrrole nanowires and multiwalled carbon nanotubes have been employed for detecting tuberculosis. Overall, NP-based electrochemical sensors represent a promising tool for the diagnosis of bTB.

Non-invasive specimen collections for Mycobacterium tuberculosis detection in free-ranging long-tailed macaques (Macaca fascicularis)

Surveillance of infectious diseases in free-ranging or wild animals has been widely conducted in many habitat-range countries after the COVID-19 episode. Thailand is located in the center of the distribution range of long-tailed macaques (Macaca fascicularis; Mf) where the animals have both frequent human contact and a high prevalence of human tuberculosis. For the large-scale detection of Mycobacterium tuberculosis complex (MTBC) using IS6110-nested PCR in free-ranging Mf, non-invasive sampling was developed using oral (via rope bait) and fecal (direct swabs of fresh feces) specimen collection. Firstly, the MTBC-IS6110-nested PCR was validated in non-invasively collected specimens, in terms of its specificity and sensitivity, and then compared with those of the invasively collected oral and rectal swabs in 24 captive MTBC-suspected Mf. After validation, these methods were applied to survey for the prevalence of shed MTBC (MTBCS) in four previously reported MTBC-infected populations. A total of 173 baited rope specimens and 204 freshly defecated excretions were collected. The limit of detection of the IS6110-nested PCR technique was 10 fg/μL and the 181-bp PCR amplicon showed 100% sequence similarity with the MTB H37Rv genome sequence. Comparing the MTBCS detection between the invasive and non-invasive collected specimens in captive suspected Mf revealed a significant correlation between the two types of oral specimens (oral swabs and baited ropes; n = 24, r2 = 1, p-value < 0.001), but fresh fecal swabs showed higher MTBCS frequencies than the rectal swabs. Moreover, the proportion of MTBCS-positive free-ranging Mf were significantly higher in the fresh fecal swabs (8.82%; 95% CI; 4.9-12.7%) than in the baited ropes (5.20%; 95% CI; 1.9-8.5%). This result indicates that oral sampling via baited ropes and fecal sampling via defecated excretion swabs can serve as ancillary specimens for MTBCS detection in free-ranging non-human primates.

Genome-wide tiled detection of circulating Mycobacterium tuberculosis cell-free DNA using Cas13

Detection of microbial cell-free DNA (cfDNA) circulating in the bloodstream has emerged as a promising new approach for diagnosing infection. Microbial diagnostics based on cfDNA require assays that can detect rare and highly fragmented pathogen nucleic acids. We now report WATSON (Whole-genome Assay using Tiled Surveillance Of Nucleic acids), a method to detect low amounts of pathogen cfDNA that couples pooled amplification of genomic targets tiled across the genome with pooled CRISPR/Cas13-based detection of these targets. We demonstrate that this strategy of tiling improves cfDNA detection compared to amplification and detection of a single targeted locus. WATSON can detect cfDNA from Mycobacterium tuberculosis in plasma of patients with active pulmonary tuberculosis, a disease that urgently needs accurate, minimally-invasive, field-deployable diagnostics. We thus demonstrate the potential for translating WATSON to a lateral flow platform. WATSON demonstrates the ability to capitalize on the strengths of targeting microbial cfDNA to address the need for point-of-care diagnostic tests for infectious diseases.

Prevalence of Mycobacterium tuberculosis Complex among Wild Rhesus Macaques and 2 Subspecies of Long-Tailed Macaques, Thailand, 2018-2022

We identified tuberculosis in 1,836 macaques from 6 wild rhesus (Macaca mulatta), 23 common long-tailed (M. fascicularis fascicularis), and 6 Burmese long-tailed (M. fascicularis aurea) macaque populations in Thailand. We captured, anesthetized, and collected throat, buccal, and rectal swab specimens from the macaques. We screened swabs for Mycobacterium tuberculosis complex (MTBC) using insertion sequence 6110-specific nested PCR. We found higher MTBC prevalence at both population and individual levels among M. mulatta than M. fascicularis fascicularis macaques; all 3 M. fascicularis aurea macaque populations were positive for tuberculosis. We found that throat swab specimens provided the best sample medium for detecting MTBC. Our results showed no difference in MTBC prevalence between male and female animals, but a higher percentage of adults were infected than subadults and juveniles. Although we detected no association between frequency of human-macaque interaction and MTBC prevalence, bidirectional zoonotic transmission should be considered a possible public health concern.

Selection of IS6110 conserved regions for the detection of Mycobacterium tuberculosis using qPCR and LAMP

IS6110 insertion sequence is a frequently used target for Mycobacterium tuberculosis detection. However, its sequence variability is studied insufficiently. We aimed to identify the most conservative and variable regions in IS6110 sequences and develop qPCR and LAMP oligonucleotide sets for the conservative regions. Using in-house Python scripts, 3609 M. tuberculosis genome sequences from the NCBI database were aligned; conservative regions were identified to design oligonucleotide sets. IS6110 fragments located within the 31-231 bp region were the most conservative and represented in genomes and were used to design qPCR and LAMP oligonucleotides. The in silico sensitivity of the qPCR oligonucleotides on the whole genome set was 99.1% and 98.4%. For the LAMP primers developed, the sensitivity was 96.9%. For qPCR, the limit of detection with 95% confidence (LoD_95%) was four IS6110 copies per reaction, with LoD_90% being 200 BCG cells per ml of artificial sputum. For LAMP, LoD_95% was 16 copies per reaction, with LoD_90% being 400 Mycobacterium bovis Bacille Calmette-Guerin (BCG) cells per ml of artificial sputum. We have demonstrated the IS6110 sequence variability and designed highly sensitive qPCR and LAMP oligonucleotides to detect M. tuberculosis.

Mycobacterium tuberculosis Lineages Associated with Mutations and Drug Resistance in Isolates from India

Current knowledge on resistance-conferring determinants in Mycobacterium tuberculosis is biased toward globally dominant lineages 2 and 4. In contrast, lineages 1 and 3 are predominant in India. In this study, we performed whole-genome sequencing of 498 MDR M. tuberculosis isolates from India to determine the prevalence of drug resistance mutations and to understand the genomic diversity. A retrospective collection of 498 M. tuberculosis isolates submitted to the National Institute for Research in Tuberculosis for phenotypic susceptibility testing between 2014 to 2016 were sequenced. Genotypic resistance prediction was performed using known resistance-conferring determinants. Genotypic and phenotypic results for 12 antituberculosis drugs were compared, and sequence data were explored to characterize lineages and their association with drug resistance. Four lineages were identified although lineage 1 predominated (43%). The sensitivity of prediction for isoniazid and rifampicin was 92% and 98%, respectively. We observed lineage-specific variations in the proportion of isolates with resistance-conferring mutations, with drug resistance more common in lineages 2 and 3. Disputed mutations (codons 430, 435, 445, and 452) in the rpoB gene were more common in isolates other than lineage 2. Phylogenetic analysis and pairwise SNP difference revealed high genetic relatedness of lineage 2 isolates. WGS based resistance prediction has huge potential, but knowledge of regional and national diversity is essential to achieve high accuracy for resistance prediction. IMPORTANCE Current knowledge on resistance-conferring determinants in Mycobacterium tuberculosis is biased toward globally dominant lineages 2 and 4. In contrast, lineages 1 and 3 are predominant in India. We performed whole-genome sequencing of 498 MDR M. tuberculosis isolates from India to determine the prevalence of drug resistance mutations and to understand genomic diversity. Four lineages were identified although lineage 1 predominated (43%). The sensitivity of prediction for isoniazid and rifampicin was 92% and 98%, respectively. We observed lineage-specific variations in the proportion of isolates with resistance-conferring mutations, with drug resistance more common in lineages 2 and 3. Disputed mutations (codons 430, 435, 445, and 452) in the rpoB gene were more common in isolates other than lineage 2. Phylogenetic analysis and pairwise SNP difference revealed high genetic relatedness of lineage 2 isolates. WGS based resistance prediction has huge potential, but knowledge of regional and national diversity is essential to achieve high accuracy for resistance prediction.

IS6110 Copy Number in Multi-Host Mycobacterium bovis Strains Circulating in Bovine Tuberculosis Endemic French Regions

IS6110 is an insertion sequence found in the Mycobacterium tuberculosis complex, to which Mycobacterium bovis belongs, which can play a role in genome plasticity and in bacterial evolution. In this study, the abundance and location of IS6110 on M. bovis genomic data of French animal field strains were studied. A first analysis was performed on a panel of 81 strains that reflect the national M. bovis population’s genetic diversity. The results show that more than one-third of them are IS6110 multicopy and that 10% have IS6110 in a high copy number (more than 6 copies). Multicopy strains are those circulating in the regions where prevalence was above the national average. Further study of 93 such strains, with an IS6110 copy number of 10-12, showed stability of IS6110 copy number and genome location over time and between host species. The correlation between M. bovis multicopy strains and high bovine tuberculosis (bTB) prevalence leads us to consider whether their epidemiological success could be partly due to genetic changes originated by IS6110 transposition.

In-depth Analysis of IS6110 Genomic Variability in the Mycobacterium tuberculosis Complex

The insertion sequence (IS) 6110 is a repetitive mobile element specific for the Mycobacterium tuberculosis complex (MTBC) used for years to diagnose and genotype this pathogen. It contains the overlapping reading frames orfA and orfB that encode a transposase. Its genetic variability is difficult to study because multiple copies are present in the genome. IS6110 is randomly located, nevertheless some preferential locations have been reported, which could be related to the behaviour of the strains. The aim of this work was to determine the intra- and inter-strain genetic conservation of this element in the MTBC. For this purpose, we analysed 158 sequences of IS6110 copies from 55 strains. Eighty-four copies were from 17 strains for which we knew all the locations in their genome. In addition, we studied 74 IS6110 copies in 38 different MTBC strains in which the location was characteristic of different families including Haarlem, LAM, S, and L6 strains. We observed mutation in 13.3% of the copies studied and we found 10 IS6110 variants in 21 copies belonging to 16 strains. The high copy number strains showed 6.2% of their IS6110 copies mutated, in contrast with the 31.1% in the low-copy-number strains. The apparently more ancient copy localised in the DR region was that with more variant copies, probably because this was the most studied location. Notably, all Haarlem and X family strains studied have an IS6110 in Rv0403c, suggesting a common origin for both families. Nevertheless, we detected a variant specific for the X family that would have occurred in this location after the phylogenetic separation. This variant does not prevent transposition although it may occur at a lower frequency, as X strains remain with low copy number (LCN) of IS6110.

Beyond to the Stable: Role of the Insertion Sequences as Epidemiological Descriptors in Corynebacterium striatum

In recent years, epidemiological studies of infectious agents have focused mainly on the pathogen and stable components of its genome. The use of these stable components makes it possible to know the evolutionary or epidemiological relationships of the isolates of a particular pathogen. Under this approach, focused on the pathogen, the identification of resistance genes is a complementary stage of a bacterial characterization process or an appendix of its epidemiological characterization, neglecting its genetic components’ acquisition or dispersal mechanisms. Today we know that a large part of antibiotic resistance is associated with mobile elements. Corynebacterium striatum, a bacterium from the normal skin microbiota, is also an opportunistic pathogen. In recent years, reports of infections and nosocomial outbreaks caused by antimicrobial multidrug-resistant C. striatum strains have been increasing worldwide. Despite the different existing mobile genomic elements, there is evidence that acquired resistance genes are coupled to insertion sequences in C. striatum. This perspective article reviews the insertion sequences linked to resistance genes, their relationship to evolutionary lineages, epidemiological characteristics, and the niches the strains inhabit. Finally, we evaluate the potential of the insertion sequences for their application as a descriptor of epidemiological scenarios, allowing us to anticipate the emergence of multidrug-resistant lineages.

Past and Present Approaches to Diagnosis of Active Pulmonary Tuberculosis

Tuberculosis disease continues to contribute to the mortality burden globally. Due to the several shortcomings of the available diagnostic methods, tuberculosis disease continues to spread. The difficulty to obtain sputum among the very ill patients and the children also affects the quick diagnosis of tuberculosis disease. These challenges warrant investigating different sample types that can provide results in a short time. Highlighted in this review are the approved pulmonary tuberculosis diagnostic methods and ongoing research to improve its diagnosis. We used the PRISMA guidelines for systematic reviews to search for studies that met the selection criteria for this review. In this review we found out that enormous biosignature research is ongoing to identify host biomarkers that can be used as predictors of active PTB disease. On top of this, more research was also being done to improve already existing diagnostic tests. Host markers required more optimization for use in different settings given their varying sensitivity and specificity in PTB endemic and non-endemic settings.

Targeted RNA-Seq Reveals the M. tuberculosis Transcriptome from an In Vivo Infection Model

Simple Summary

High-throughput sequencing techniques such as RNA-seq allow a more detailed characterization of the gene expression profile during in vivo infections. However, using this strategy for intracellular pathogens such as Mycobacterium tuberculosis (Mtb) entails technical limitations. Some authors have resorted to flow cytometers to separate infected cells or significantly increase sequencing depth to obtain pathogens’ gene expression. However, these options carry additional expenses in specialized equipment. We propose an experimental protocol based on differential cell lysis and a probe-based ribosomal depletion to determine the gene expression of Mtb and its host during in vivo infection. This method allowed us to increase the number of observed expressed genes from 13 using a traditional RNA-seq approach to 702. In addition, we observed the expression of genes essential for establishing the infection, codifying proteins such as PE-PGRS, lipoproteins lppN and LpqH, and three ncRNAs (small RNA MTS2823, transfer-messenger RNA RF00023, and ribozyme RF00010). We believe our method represents a valuable alternative to current RNA-seq approaches to study host–pathogen interactions and will help explore host–pathogen mechanisms in tuberculosis and other similar models of intracellular infections.

Abstract

The study of host-pathogen interactions using in vivo models with intracellular pathogens like Mycobacterium tuberculosis (Mtb) entails technical limitations, such as: (i) Selecting an efficient differential lysis system to enrich the pathogen cells; (ii) obtaining sufficient high-quality RNA; and (iii) achieving an efficient rRNA depletion. Thus, some authors had used flow cytometers to separate infected cells or significantly increase the sequencing depth of host–pathogen RNA libraries to observe the pathogens’ gene expression. However, these options carry additional expenses in specialized equipment typically not available for all laboratories. Here, we propose an experimental protocol involving differential cell lysis and a probe-based ribosomal depletion to determine the gene expression of Mtb and its host during in vivo infection. This method increased the number of observed pathogen-expressed genes from 13 using the traditional RNA-seq approach to 702. After eliminating rRNA reads, we observed that 61.59% of Mtb sequences represented 702 genes, while 38.41% represented intergenic regions. Some of the most expressed genes codified for IS1081 (Rv2512c) transposase and eight PE-PGRS members, such as PGRS49 and PGRS50. As expected, a critical percent of the expressed genes codified for secreted proteins essential for infection, such as PE68, lppN, and LpqH. Moreover, three Mtb ncRNAs were highly expressed (small RNA MTS2823, transfer-messenger RNA RF00023, and ribozyme RF00010). Many of the host-expressed genes were related to the inflammation process and the expression of surfactant proteins such as the Sftpa and Sftpc, known to bind Mtb to alveolar macrophages and mi638, a microRNA with no previous associations with pulmonary diseases. The main objective of this study is to present the method, and a general catalog of the Mtb expressed genes at one point of the in vivo infection. We believe our method represents a different approach to the existing ones to study host–pathogen interactions in tuberculosis and other similar intracellular infections, without the necessity of specialized equipment.

Intragenic Distribution of IS6110 in Clinical Mycobacterium tuberculosis Strains: Bioinformatic Evidence for Gene Disruption Leading to Underdiagnosed Antibiotic Resistance

Antibiotic resistance is a global challenge for tuberculosis control, and accelerating its diagnosis is critical for therapy decisions and controlling transmission. Genotype-based molecular diagnostics now play an increasing role in accelerating the detection of such antibiotic resistance, but their accuracy depends on the instructed detection of genetic variations. Genetic mobile elements such as IS6110 are established sources of genetic variation in Mycobacterium tuberculosis, but their implication in clinical antibiotic resistance has thus far been unclear. Here, we describe the discovery of an intragenic IS6110 insertion into Rv0678 that caused antibiotic resistance in an in vitro-selected M. tuberculosis isolate. The subsequent development of bioinformatics scripts allowed genome-wide analysis of intragenic IS6110 insertions causing gene disruptions in 6,426 clinical M. tuberculosis strains. This analysis identified 10,070 intragenic IS6110 insertions distributed among 333 different genes. Focusing on genes whose disruption leads to antibiotic resistance, 12 clinical isolates were identified with high confidence to be resistant to bedaquiline, clofazimine, pyrazinamide, ethionamide, and para-aminosalicylic acid because of an IS6110-mediated gene disruption event. A number of these IS6110-mediated resistant strains had identical genomic distributions of IS6110 elements and likely represent transmission events of a single resistant isolate. These data provide strong evidence that IS6110-mediated gene disruption is a clinically relevant mechanism of antibiotic resistance in M. tuberculosis that should be considered for molecular diagnostics. Concomitantly, this analysis provides a list of 333 IS6110-disrupted genes in clinical tuberculosis isolates that can be deemed nonessential for human infection.

IMPORTANCE To help control the spread of drug-resistant tuberculosis and to guide treatment choices, it is important that rapid and accurate molecular diagnostic tools are used. Current molecular diagnostic tools detect the most common antibiotic-resistance-conferring mutations in the form of single nucleotide changes, small deletions, or insertions. Mobile genetic elements, named IS6110, are also known to move within the M. tuberculosis genome and cause significant genetic variations, although the role of this variation in clinical drug resistance remains unclear. In this work, we show that both in vitro and in data analyzed from 6,426 clinical M. tuberculosis strains, IS6110 elements are found that disrupt specific genes essential for the function of a number of pivotal antituberculosis drugs. By providing ample evidence of clinically relevant IS6110-mediated drug resistance, we believe that this shows that this form of genetic variation must not be overlooked in molecular diagnostics of drug resistance.

Phenotypic and genotypic features of the Mycobacterium tuberculosis lineage 1 subgroup in central Vietnam

Mycobacterium tuberculosis (Mtb) has different features depending on different geographic areas. We collected Mtb strains from patients with smear-positive pulmonary tuberculosis in Da Nang, central Vietnam. Using a whole genome sequencing platform, including genome assembly complemented by long-read-sequencing data, genomic characteristics were studied. Of 181 Mtb isolates, predominant Vietnamese EAI4_VNM and EAI4-like spoligotypes (31.5%), ZERO strains (5.0%), and part of EAI5 (11.1%) were included in a lineage-1 (L1) sublineage, i.e., L1.1.1.1. These strains were found less often in younger people, and they genetically clustered less frequently than other modern strains. Patients infected with ZERO strains demonstrated less lung infiltration. A region in RD2bcg spanning six loci, i.e., PE_PGRS35, cfp21, Rv1985c, Rv1986, Rv1987, and erm(37), was deleted in EAI4_VNM, EAI4-like, and ZERO strains, whereas another 118 bp deletion in furA was specific only to ZERO strains. L1.1.1.1-sublineage-specific deletions in PE_PGRS4 and PE_PGRS22 were also identified. RD900, seen in ancestral lineages, was present in majority of the L1 members. All strains without IS6110 (5.0%) had the ZERO spoligo-pattern. Distinctive features of the ancestral L1 strains provide a basis for investigation of the modern versus ancestral Mtb lineages and allow consideration of countermeasures against this heterogeneous pathogen.

A case of childhood tuberculosis from late mediaeval Somerset, England

BACKGROUND

The remains of a 3-5 year-old child from the late mediaeval cemetery serving the Priory of St. Peter and St. Paul, Taunton, Somerset, UK was the subject of an aDNA study.

OBJECTIVE

The aim was to distinguish between two differential diagnoses suggested by earlier osteological examination of the remains; either tuberculosis or Langerhans cell histiocytosis.

FINDINGS

The remains tested positive for MTB complex markers, corroborating this diagnosis reached on osteological grounds. Based on positivity for the mtp40 element and a deletion in the pks15/1 locus, we conclude that infection was due to a strain of the human pathogen M.tuberculosis belonging to lineage 4. Although DNA recovered from the case was heavily fragmented, sex determination by amelogenin PCR suggested these are the remains of a young male child. The findings are discussed considering additions to the literature since the original report.

CONCLUSIONS

Descriptions of tuberculosis in children from this period are rare and burial Sk2077 represents the first UK example of a pre-adolescent individual to have a molecular diagnosis combined with osteological pathology. This provides an important reference of childhood tuberculosis and insight into the likely presence of tuberculosis in the mediaeval adult population served by this cemetery.

CRISPRbuilder-TB: "CRISPR-builder for tuberculosis". Exhaustive reconstruction of the CRISPR locus in mycobacterium tuberculosis complex using SRA

Mycobacterium tuberculosis complex (MTC) CRISPR locus diversity has long been studied solely investigating the presence/absence of a known set of spacers. Unveiling the genetic mechanisms of its evolution requires a more exhaustive reconstruction in a large amount of representative strains.

In this article, we point out and resolve, with a new pipeline, the problem of CRISPR reconstruction based directly on short read sequences in M. tuberculosis. We first show that the process we set up, that we coin as “CRISPRbuilder-TB” (https://github.com/cguyeux/CRISPRbuilder-TB), allows an efficient reconstruction of simulated or real CRISPRs, even when including complex evolutionary steps like the insertions of mobile elements. Compared to more generalist tools, the whole process is much more precise and robust, and requires only minimal manual investigation. Second, we show that more than 1/3 of the currently complete genomes available for this complex in the public databases contain largely erroneous CRISPR loci. Third, we highlight how both the classical experimental in vitro approach and the basic in silico spoligotyping provided by existing analytic tools miss a whole diversity of this locus in MTC, by not capturing duplications, spacer and direct repeats variants, and IS6110 insertion locations. This description is extended in a second article that describes MTC-CRISPR diversity and suggests general rules for its evolution. This work opens perspectives for an in-depth exploration of M. tuberculosis CRISPR loci diversity and of mechanisms involved in its evolution and its functionality, as well as its adaptation to other CRISPR locus-harboring bacterial species.

In this article, we tackle the bioinformatical issue of the reconstruction of the Mycobacterium tuberculosis complex CRISPR locus using short read sequences without requiring genome assembly. We first show that many complete genomes, as found in public databases and often reconstructed by de novo assemblies, often contain errors on this locus as well as on other repeated sequences. We provide an in-depth description of our new method, designated as ‘CRISPRbuilder-TB’, and we show that our method provides much more exhaustive and reliable information (on DR variants, spacer diversity, global structure) than Crass and CRISPR_detector. The new and unsuspected genomic diversity we detected is described in a companion paper. Scripts are available to adapt the tool to other species.

Unexpected diversity of CRISPR unveils some evolutionary patterns of repeated sequences in Mycobacterium tuberculosis

BACKGROUND

Diversity of the CRISPR locus of Mycobacterium tuberculosis complex has been studied since 1997 for molecular epidemiology purposes. By targeting solely the 43 spacers present in the two first sequenced genomes (H37Rv and BCG), it gave a biased idea of CRISPR diversity and ignored diversity in the neighbouring cas-genes.

RESULTS

We set up tailored pipelines to explore the diversity of CRISPR-cas locus in Short Reads. We analyzed data from a representative set of 198 clinical isolates as evidenced by well-characterized SNPs. We found a relatively low diversity in terms of spacers: we recovered only the 68 spacers that had been described in 2000. We found no partial or global inversions in the sequences, letting always the Direct Variant Repeats (DVR) in the same order. In contrast, we found an unexpected diversity in the form of: SNPs in spacers and in Direct Repeats, duplications of various length, and insertions at various locations of the IS6110 insertion sequence, as well as blocks of DVR deletions. The diversity was in part specific to lineages. When reconstructing evolutionary steps of the locus, we found no evidence for SNP reversal. DVR deletions were linked to recombination between IS6110 insertions or between Direct Repeats.

CONCLUSION

This work definitively shows that CRISPR locus of M. tuberculosis did not evolve by classical CRISPR adaptation (incorporation of new spacers) since the last most recent common ancestor of virulent lineages. The evolutionary mechanisms that we discovered could be involved in bacterial adaptation but in a way that remains to be identified.

Improvement of Mycobacterium tuberculosis detection in sputum using DNA extracted by sonication

Tuberculosis (TB) is one of the infectious diseases with high mortality in the world. DNA amplification techniques have been used to overcome barriers to the diagnosis of this disease. However, the success of these methodologies is highly dependent on the DNA obtained from the sample. This study was carried out to verify whether the DNA extracted by sonication (in house method) could yield suitable DNA for amplification by real-time PCR (qPCR). Sixty sputum samples were submitted to DNA extraction using sonication compared to a commercial method (Detect-TB kit, Labtest/MG-Brazil). All DNA samples were amplified by qPCR for IS6110 region (IS6110-qPCR/SYBR Green assay). Out of 60 samples, 40 were positive for TB; of these, all had positive results when extracted by sonication (100%) and 80% when extracted by the commercial method. The limit of detection (LOD) of Mycobacterium tuberculosis (H37Rv strain) by qPCR was 14 CFU/mL when the DNA was extracted by sonication, compared to countless colonies when extracted by commercial kit. In conclusion, the sonication protocol (without purification step) proved to be a simple, fast, and suitable method for obtaining DNA for use in qPCR from sputum samples.

Mycobacterium tuberculosis complex in wildlife: Review of current applications of antemortem and postmortem diagnosis

Tuberculosis (TB) is a chronic inflammatory and zoonotic disease caused by Mycobacterium tuberculosis complex (MTBC) members, which affects various domestic animals, wildlife, and humans. Some wild animals serve as reservoir hosts in the transmission and epidemiology of the disease. Therefore, the monitoring and surveillance of both wild and domestic hosts are critical for prevention and control strategies. For TB diagnosis, the single intradermal tuberculin test or the single comparative intradermal tuberculin test, and the gamma-interferon test, which is regarded as an ancillary test, are used. Postmortem examination can identify granulomatous lesions compatible with a diagnosis of TB. In contrast, smears of the lesions can be stained for acid-fast bacilli, and samples of the affected organs can be subjected to histopathological analyses. Culture is the gold standard test for isolating mycobacterial bacilli because it has high sensitivity and specificity compared with other methods. Serology for antibody detection allows the testing of many samples simply, rapidly, and inexpensively, and the protocol can be standardized in different laboratories. Molecular biological analyses are also applicable to trace the epidemiology of the disease. In conclusion, reviewing the various techniques used in MTBC diagnosis can help establish guidelines for researchers when choosing a particular diagnostic method depending on the situation at hand, be it disease outbreaks in wildlife or for epidemiological studies. This is because a good understanding of various diagnostic techniques will aid in monitoring and managing emerging pandemic threats of infectious diseases from wildlife and also preventing the potential spread of zoonotic TB to livestock and humans. This review aimed to provide up-to-date information on different techniques used for diagnosing TB at the interfaces between wildlife, livestock, and humans.

Molecular diagnostics for verification of pleural tuberculosis in Morocco

Pleural tuberculosis (pTB) is a very common form of extrapulmonary tuberculosis (TB). pTB diagnostics represents a major burning challenge worldwide due to the limitations of available conventional diagnostic tools. These latter include microscopic examination of the pleural fluid for acid-fast bacilli, mycobacterial culture of pleural fluid in solid or liquid media, sputum or pleural tissue, and histopathological examination of pleural tissue; these tests have recognized limitations for clinical use. Hence, to overcome these limitations, attention has been devoted to new nucleic acid amplification (NAA) diagnostic tests such as the polymerase chain reaction (PCR) and real-time PCR (RT-PCR), owing to their accuracy, rapidity, high sensitivity and specificity. Within this context, this prospective study was conducted to evaluate the performance of molecular diagnosis methods for differentiation between tuberculosis and non-tuberculosis pleural effusions. Fifty patients with pleural effusion were enrolled in this prospective study in Rabat, Morocco. The efficacy of conventional polymerase chain reaction (PCR) in the diagnostics of tuberculous pleurisy by targeting IS6110 and mycobacterial internal transcribed spacer (MYITS) was evaluated compared to histopathologic examination and culture data. Our results showed that IS6110 PCR could “rule in” pTB, the sensitivity and specificity being 41.6% and 85.7%, respectively. Therefore, the findings confirmed that molecular tests exert a relatively high specificity in EPTB but lower sensitivity, thus a positive test is considered as a pTB case whereas negative one cannot exclude the disease. Although the study was limited by a small sample size, it adds to the body of evidence of usefulness of molecular testing as adjuncts to histopathologic examination for accurate diagnosis of pTB, to treat timely and to avoid the emergence and spread of drug resistant pTB. However, further efforts should be made to increase the sensitivity of NAA methods and to identify the best molecular targets to be useful in clinical practice.

Impact of insertion sequences on convergent evolution of Shigella species

Shigella species are specialised lineages of Escherichia coli that have converged to become human-adapted and cause dysentery by invading human gut epithelial cells. Most studies of Shigella evolution have been restricted to comparisons of single representatives of each species; and population genomic studies of individual Shigella species have focused on genomic variation caused by single nucleotide variants and ignored the contribution of insertion sequences (IS) which are highly prevalent in Shigella genomes. Here, we investigate the distribution and evolutionary dynamics of IS within populations of Shigella dysenteriae Sd1, Shigella sonnei and Shigella flexneri. We find that five IS (IS1, IS2, IS4, IS600 and IS911) have undergone expansion in all Shigella species, creating substantial strain-to-strain variation within each population and contributing to convergent patterns of functional gene loss within and between species. We find that IS expansion and genome degradation are most advanced in S. dysenteriae and least advanced in S. sonnei; and using genome-scale models of metabolism we show that Shigella species display convergent loss of core E. coli metabolic capabilities, with S. sonnei and S. flexneri following a similar trajectory of metabolic streamlining to that of S. dysenteriae. This study highlights the importance of IS to the evolution of Shigella and provides a framework for the investigation of IS dynamics and metabolic reduction in other bacterial species.

Xpert MTB/RIF Ultra versus Xpert MTB/RIF for the diagnosis of tuberculous meningitis: a prospective, randomised, diagnostic accuracy study

BACKGROUND

Xpert MTB/RIF Ultra (Xpert Ultra) might have higher sensitivity than its predecessor, Xpert MTB/RIF (Xpert), but its role in tuberculous meningitis diagnosis is uncertain. We aimed to compare Xpert Ultra with Xpert for the diagnosis of tuberculous meningitis in HIV-uninfected and HIV-infected adults.

METHODS

In this prospective, randomised, diagnostic accuracy study, adults (≥16 years) with suspected tuberculous meningitis from a single centre in Vietnam were randomly assigned to cerebrospinal fluid testing by either Xpert Ultra or Xpert at baseline and, if treated for tuberculous meningitis, after 3-4 weeks of treatment. Test performance (sensitivity, specificity, and positive and negative predictive values) was calculated for Xpert Ultra and Xpert and compared against clinical and mycobacterial culture reference standards. Analyses were done for all patients and by HIV status.

FINDINGS

Between Oct 16, 2017, and Feb 10, 2019, 205 patients were randomly assigned to Xpert Ultra (n=103) or Xpert (n=102). The sensitivities of Xpert Ultra and Xpert for tuberculous meningitis diagnosis against a reference standard of definite, probable, and possible tuberculous meningitis were 47·2% (95% CI 34·4-60·3; 25 of 53 patients) for Xpert Ultra and 39·6% (27·6-53·1; 21 of 53) for Xpert (p=0·56); specificities were 100·0% (95% CI 92·0-100·0; 44 of 44) and 100·0% (92·6-100·0; 48 of 48), respectively. In HIV-negative patients, the sensitivity of Xpert Ultra was 38·9% (24·8-55·1; 14 of 36) versus 22·9% (12·1-39·0; eight of 35) by Xpert (p=0·23). In HIV co-infected patients, the sensitivities were 64·3% (38·8-83·7; nine of 14) for Xpert Ultra and 76·9% (49·7-91·8; ten of 13) for Xpert (p=0·77). Negative predictive values were 61·1% (49·6-71·5) for Xpert Ultra and 60·0% (49·0-70·0) for Xpert. Against a reference standard of mycobacterial culture, sensitivities were 90·9% (72·2-97·5; 20 of 22 patients) for Xpert Ultra and 81·8% (61·5-92·7; 18 of 22) for Xpert (p=0·66); specificities were 93·9% (85·4-97·6; 62 of 66) and 96·9% (89·5-91·2; 63 of 65), respectively. Six (22%) of 27 patients had a positive test by Xpert Ultra after 4 weeks of treatment versus two (9%) of 22 patients by Xpert.

INTERPRETATION

Xpert Ultra was not statistically superior to Xpert for the diagnosis of tuberculous meningitis in HIV-uninfected and HIV-infected adults. A negative Xpert Ultra or Xpert test does not rule out tuberculous meningitis. New diagnostic strategies are urgently required.

FUNDING

Wellcome Trust and the Foundation for Innovative New Diagnostics.

Towards next-generation diagnostics for tuberculosis: identification of novel molecular targets by large-scale comparative genomics

MOTIVATION

Tuberculosis (TB) remains one of the main causes of death worldwide. The long and cumbersome process of culturing Mycobacterium tuberculosis complex (MTBC) bacteria has encouraged the development of specific molecular tools for detecting the pathogen. Most of these tools aim to become novel TB diagnostics, and big efforts and resources are invested in their development, looking for the endorsement of the main public health agencies. Surprisingly, no study has been conducted where the vast amount of genomic data available is used to identify the best MTBC diagnostic markers.

RESULTS

In this work, we used large-scale comparative genomics to identify 40 MTBC-specific loci. We assessed their genetic diversity and physiological features to select 30 that are good targets for diagnostic purposes. Some of these markers could be used to assess the physiological status of the bacilli. Remarkably, none of the most used MTBC markers is in our catalog. Illustrating the translational potential of our work, we develop a specific qPCR assay for quantification and identification of MTBC DNA. Our rational design of targeted molecular assays for TB could be used in many other fields of clinical and basic research.

AVAILABILITY AND IMPLEMENTATION

The database of non-tuberculous mycobacteria assemblies can be accessed at: 10.5281/zenodo.3374377.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Isolation and comparative genomics of Mycobacterium tuberculosis isolates from cattle and their attendants in South India

The major human pathogen Mycobacterium tuberculosis is rarely reported to cause disease in other animals. Cases in livestock are thought to occur through contact with infected handlers, but previous studies evaluating putative livestock-human transmission used typing techniques with limited resolution. Here, we undertook cross-sectional surveillance for tuberculosis in 271 livestock handlers and 167 cattle on three farms in Chennai, India and defined the relatedness of cultured isolates using whole genome sequencing. Humans and livestock were screened for active mycobacterial infection, and opportunistic post-mortem examination was performed on comparative intradermal test-positive cattle that died. Four cattle and 6 handlers on two farms were culture-positive for M. tuberculosis; M. bovis was not isolated. All 10 isolates (one from each case) belonged to Lineage 1. Pairwise genome comparisons of single nucleotide polymorphism (SNP) differences ranged from 1 to 600 SNPs, but 3 isolate pairs were less than 5 SNPs different. Two pairs were from handlers and the third pair were from two cattle on the same farm. The minimum pairwise SNP difference between a cattle and human isolate was >250 SNPs. Our study confirms the presence of M. tuberculosis infection in cattle in India, sequencing of which characterised relatedness between human and cattle-derived isolates.

The role of IS6110 in micro- and macroevolution of Mycobacterium tuberculosis lineage 2

The insertion sequence 6110 (IS6110) is the most studied transposable element in the Mycobacterium tuberculosis complex species. The element plays a significant role in genome plasticity of this important human pathogen, but still many causes and consequences of its transposition have not been fully studied. Here, we analyzed insertion sites for 902 Mycobacterium tuberculosis lineage 2 strains using whole-genome sequencing data. In total, 17,972 insertions were found, corresponding to 827 independent positions in the genome of the reference strain H37Rv. To trace the history of IS6110 expansion since proto-Beijing strains up to modern sublineages, we looked at the distribution of IS6110 across the genome-wide SNP-based phylogenetic tree. This analysis demonstrated a stepwise transposition of IS6110 that occurs by «copy-and-paste» mechanism. Additionally, we detected evolutionary-scale and sublineage-specific integration sites, which can be used for typing and for understanding the reasons for the success of the lineage. A significant part of such insertions affected the genes that are essential for the pathogen. Finally, we identified and confirmed deletions that occurred between differently oriented elements, which is uncommon for this family of insertion elements and appears to be another mechanism of genome variability.

First insights on the genetic diversity of MDR Mycobacterium tuberculosis in Lebanon

Background

Lebanon hosts a heterogeneous population coming from underdeveloped and developing countries, resulting in increasing incidences of tuberculosis over the past years. The genetic heterogeneity and lineages associated with tuberculosis, along with their resistance determinants have not been studied at the genomic level previously in the region.

Methods

Isolates were recovered from the American University of Beirut Medical Center (AUBMC). Antimicrobial susceptibility profiles were determined using the MGIT automated system for the first-line drugs at AUBMC, while second-line drug susceptibility was tested at Mayo Clinic Laboratories. Whole Genome Sequencing (WGS) was performed to classify mycobacterial lineages and highlight single nucleotide mutations causing resistance to both 1st line and 2nd line antimicrobials. wgSNP analysis provided insights on the phylogeny of the isolates along with spoligotyping and core genomic SNVs, IS6110 insertion sites, and variable number tandem repeats (VNTR).

Results

The analyzed isolates carry distinct resistance determinants to isoniazid, rifampicin, ethambutol, quinolones, and streptomycin. The isolates belonged to different lineages including the Euro/American lineage (Lineage 4) (53.8%), M. bovis (15.4%) and Delhi/Central Asia (Lineage 1) (15.4%), Beijing/East Asia (Lineage 2) (7.7%), and East Africa/Indian Ocean lineage (Lineage 3) (7.7%) showing great phylogenetic differences at the genomic level.

Conclusions

The population diversity in Lebanon holds an equally diverse and uncharacterized population of drug resistant mycobacteria. To achieve the WHO “END-TB” milestones of 2025 and 2035, Lebanon must decrease TB incidences by 95% in the next decade. This can only be done through WGS-based patient centered diagnosis with higher throughput and genomic resolution to improve treatment outcomes and to monitor transmission patterns.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3626-3) contains supplementary material, which is available to authorized users.

Approximate, simultaneous comparison of microbial genome architectures via syntenic anchoring of quiver representations

Motivation

A long-standing limitation in comparative genomic studies is the dependency on a reference genome, which hinders the spectrum of genetic diversity that can be identified across a population of organisms. This is especially true in the microbial world where genome architectures can significantly vary. There is therefore a need for computational methods that can simultaneously analyze the architectures of multiple genomes without introducing bias from a reference.

Results

In this article, we present Ptolemy: a novel method for studying the diversity of genome architectures-such as structural variation and pan-genomes-across a collection of microbial assemblies without the need of a reference. Ptolemy is a 'top-down' approach to compare whole genome assemblies. Genomes are represented as labeled multi-directed graphs-known as quivers-which are then merged into a single, canonical quiver by identifying 'gene anchors' via synteny analysis. The canonical quiver represents an approximate, structural alignment of all genomes in a given collection encoding structural variation across (sub-) populations within the collection. We highlight various applications of Ptolemy by analyzing structural variation and the pan-genomes of different datasets composing of Mycobacterium, Saccharomyces, Escherichia and Shigella species. Our results show that Ptolemy is flexible and can handle both conserved and highly dynamic genome architectures. Ptolemy is user-friendly-requires only FASTA-formatted assembly along with a corresponding GFF-formatted file-and resource-friendly-can align 24 genomes in ∼10 mins with four CPUs and <2 GB of RAM.

Availability and implementation

Github: https://github.com/AbeelLab/ptolemy.

Supplementary information

Supplementary data are available at Bioinformatics online.

Automated TruTip nucleic acid extraction and purification from raw sputum

Automated nucleic acid extraction from primary (raw) sputum continues to be a significant technical challenge for molecular diagnostics. In this work, we developed a prototype open-architecture, automated nucleic acid workstation that includes a mechanical homogenization and lysis function integrated with heating and TruTip purification; optimized an extraction protocol for raw sputum; and evaluated system performance on primary clinical specimens. Eight samples could be processed within 70 min. The system efficiently homogenized primary sputa and doubled nucleic acid recovery relative to an automated protocol that did not incorporate sample homogenization. Nucleic acid recovery was at least five times higher from raw sputum as compared to that of matched sediments regardless of smear or culture grade, and the automated workstation reproducibly recovered PCR-detectable DNA to at least 80 CFU mL-1 raw sputum. M. tuberculosis DNA was recovered and detected from 122/123 (99.2%) and 124/124 (100%) primary sputum and sediment extracts, respectively. There was no detectable cross-contamination across 53 automated system runs and amplification or fluorescent inhibitors (if present) were not detectable. The open fluidic architecture of the prototype automated workstation yields purified sputum DNA that can be used for any molecular diagnostic test. The ability to transfer TruTip protocols between personalized, on-demand pipetting tools and the fully automated workstation also affords public health agencies an opportunity to standardize sputum nucleic acid sample preparation procedures, reagents, and quality control across multiple levels of the health care system.

New insights into the transposition mechanisms of IS6110 and its dynamic distribution between Mycobacterium tuberculosis Complex lineages

The insertion Sequence IS6110, only present in the pathogens of the Mycobacterium tuberculosis Complex (MTBC), has been the gold-standard epidemiological marker for TB for more than 25 years, but biological implications of IS6110 transposition during MTBC adaptation to humans remain elusive. By studying 2,236 clinical isolates typed by IS6110-RFLP and covering the MTBC, we remarked a lineage-specific content of IS6110 being higher in modern globally distributed strains. Once observed the IS6110 distribution in the MTBC, we selected representative isolates and found a correlation between the normalized expression of IS6110 and its abundance in MTBC chromosomes. We also studied the molecular regulation of IS6110 transposition and we found a synergistic action of two post-transcriptional mechanisms: a -1 ribosomal frameshift and a RNA pseudoknot which interferes translation. The construction of a transcriptionally active transposase resulted in 20-fold increase of the transposition frequency. Finally, we examined transposition in M. bovis and M. tuberculosis during laboratory starvation and in a mouse infection model of TB. Our results shown a higher transposition in M. tuberculosis, that preferably happens during TB infection in mice and after one year of laboratory culture, suggesting that IS6110 transposition is dynamically adapted to the host and to adverse growth conditions.

Since the pioneering discovery of transposition by Barbara McClintock in eukaryotes and later in prokaryotes by Robert W. Hedges and Alan E. Jacob, it has become clear the key role of mobile genetics elements in chromosome remodelling, microbial evolution and host adaptation. The insertion sequence IS6110 is widely recognized for its utility in TB diagnosis and epidemiology because it is only present in the M. tuberculosis Complex (MTBC) and its transposition provides an excellent chromosomal polymorphic variability allowing the study of recent TB transmission. This inherent feature of IS6110 leads us to hypothesize that IS6110 plays a crucial role during the TB infectious cycle. However, the biological significance of IS6110 has been hindered by its almost exclusive use as an epidemiological marker. Here, we study the regulatory mechanisms and the distribution of IS6110 in the different MTBC lineages. We discuss the potential biological implications of IS6110, that is much more than an excellent TB epidemiological tool. Since IS6110 could play an important role in the adaptation of MTBC to the host, this study opens new avenues to decipher the biological roles of IS6110 in TB pathogenesis.

Classification of pathogenic microbes using a minimal set of single nucleotide polymorphisms derived from whole genome sequences

In a context specific manner, Intra-species genomic variation plays an important role in phenotypic diversity observed among pathogenic microbes. Efficient classification of these pathogens is important for diagnosis and treatment of several infectious diseases. NGS technologies have provided access to wealth of data that can be utilized to discover important markers for pathogen classification. In this paper, we described three different approaches (Jensen-Shannon divergence, random forest and Shewhart control chart) for identification of a minimal set of SNPs that can be used for classification of organisms. These methods are generic and can be implemented for analysis of any organism. We have shown usefulness of these approaches for analysis of Mycobacterium tuberculosis and Escherichia coli isolates. We were able to identify a minimal set of 18 SNPs that can be used as molecular markers for phylogroup based classification and 8 SNPs for pathogroup based classification of E. coli.

Spoligotyping and whole-genome sequencing analysis of lineage 1 strains of Mycobacterium tuberculosis in Da Nang, Vietnam

Background

Spacer oligonucleotide typing (spoligotyping), a widely used, classical genotyping method for Mycobacterium tuberculosis complex (MTBC), is a PCR-based dot-blot hybridization technique to detect the genetic diversity of the direct repeat (DR) region. Of the seven major MTBC lineages in the world, lineage 1 (Indo-Oceanic) mostly corresponds to the East African–Indian (EAI) spoligotype family in East Africa and Southeast Asia.

Objectives

We investigated the genomic features of Vietnamese lineage 1 strains, comparing spoligotype patterns using whole-genome sequencing (WGS) data.

Methods

M. tuberculosis strains isolated in Da Nang, Vietnam were subjected to conventional spoligotyping, followed by WGS analysis using a high-throughput sequencer. Vietnamese lineage 1 strains were further analyzed with other lineage 1 strains obtained from a public database.

Results

Indicating a major spoligotype in Da Nang, 86 (46.2%) of the 186 isolates belonged to the EAI family or lineage 1. Although typical EAI4-VNM strains are characterized by the deletion of spacers 26 and 27, 65 (75.6%) showed ambiguous signals on spacer 26. De novo assembly of the entire DR region and in silico spoligotyping analysis suggested the absence of spacer 26, and direct sequencing revealed that the 17^th spacer sequence not used for conventional typing, was cross-hybridized to the spacer 26 probe. Vietnamese EAI4-VNM, other EAI-like strains, and those showing a non-EAI pattern lacking many spacers formed a monophyletic group separate from other EAI families in the world.

Conclusion

Information about the alignment of spacers in the entire DR region obtained from WGS data provides a clue for the determination of experimentally ambiguous spoligo patterns. WGS data also helped to analyze the hidden relationships between apparently distinct spoligo patterns.

Comparative genomic analysis of Mycobacterium tuberculosis Beijing-like strains revealed specific genetic variations associated with virulence and drug resistance

Isolates of the Mycobacterium tuberculosis lineage 2/East-Asian are considered one of the most successful strains due to their increased pathogenicity, hyper-virulence associated with drug resistance, and high transmission. Recent studies in Colombia have shown that the Beijing-like genotype is associated with multidrug-resistance and high prevalence in the southwest of the country, but the genetic basis of its success in dissemination is unknown. In contribution to this matter, we obtained the whole sequences of six genomes of clinical isolates assigned to the Beijing-like genotype. The genomes were compared with the reference genome of M. tuberculosis H37Rv and 53 previously published M. tuberculosis genomes. We found that the six Beijing-like isolates belong to a modern Beijing sub-lineage and share specific genomic variants: i.e. deletion in the PPE8 gene, in Rv3806c (ubiA) responsible of high ethambutol resistance and in Rv3862c (whiB6) which is involved in granuloma formation and virulence, are some of them. Moreover, each isolated has exclusively single nucleotide polymorphisms (SNPs) in genes related with cell wall processes and cell metabolism. We identified polymorphisms in genes related to drug resistance that could explain the drug-resistant phenotypes found in the six isolates from Colombia. We hypothesize that changes due to these genetic variations contribute to the success of these strains. Finally, we analyzed the IS6110 insertion sequences finding very low variance between them, suggesting that SNPs is the major cause of variability found in Beijing-like strains circulating in Colombia.

SMRT genome assembly corrects reference errors, resolving the genetic basis of virulence in Mycobacterium tuberculosis

BACKGROUND

The genetic basis of virulence in Mycobacterium tuberculosis has been investigated through genome comparisons of virulent (H37Rv) and attenuated (H37Ra) sister strains. Such analysis, however, relies heavily on the accuracy of the sequences. While the H37Rv reference genome has had several corrections to date, that of H37Ra is unmodified since its original publication.

RESULTS

Here, we report the assembly and finishing of the H37Ra genome from single-molecule, real-time (SMRT) sequencing. Our assembly reveals that the number of H37Ra-specific variants is less than half of what the Sanger-based H37Ra reference sequence indicates, undermining and, in some cases, invalidating the conclusions of several studies. PE_PPE family genes, which are intractable to commonly-used sequencing platforms because of their repetitive and GC-rich nature, are overrepresented in the set of genes in which all reported H37Ra-specific variants are contradicted. Further, one of the sequencing errors in H37Ra masks a true variant in common with the clinical strain CDC1551 which, when considered in the context of previous work, corresponds to a sequencing error in the H37Rv reference genome.

CONCLUSIONS

Our results constrain the set of genomic differences possibly affecting virulence by more than half, which focuses laboratory investigation on pertinent targets and demonstrates the power of SMRT sequencing for producing high-quality reference genomes.

Structure and variation of CRISPR and CRISPR-flanking regions in deleted-direct repeat region Mycobacterium tuberculosis complex strains

Background

CRISPR and CRISPR-flanking genomic regions are important for molecular epidemiology of Mycobacterium tuberculosis complex (MTBC) strains, and potentially for adaptive immunity to phage and plasmid DNA, and endogenous roles in the bacterium. Genotyping in the Israel National Mycobacterium Reference Center Tel-Aviv of over 1500 MTBC strains from 2008–2013 showed three strains with validated negative 43-spacer spoligotypes, that is, with putatively deleted direct repeat regions (deleted-DR/CRISPR regions). Two isolates of each of three negative spoligotype MTBC (a total of 6 isolates) were subjected to Next Generation Sequencing (NGS). As positive controls, NGS was performed for three intact-DR isolates belonging to T3_Eth, the largest multiple-drug-resistant (MDR)-containing African-origin cluster in Israel. Other controls consisted of NGS reads and complete whole genome sequences from GenBank for 20 intact-DR MTBC and for 1 deleted-DR MTBC strain recognized as CAS by its defining RD deletion.

Results

NGS reads from negative spoligotype MTBC mapped to reference H37Rv NC_000962.3 suggested that the DR/CRISPR regions were completely deleted except for retention of the middle IS6110 mobile element. Clonally specific deletion of CRISPR-flanking genes also was observed, including deletion of at least cas2 and cas1 genes. Genomic RD deletions defined lineages corresponding to the major spoligotype families Beijing, EAI, and Haarlem, consistent with 24 loci MIRU-VNTR profiles. Analysis of NGS reads, and analysis of contigs obtained by manual PCR confirmed that all 43 gold standard DR/CRISPR spacers were missing in the deleted-DR genomes.

Conclusions

Although many negative spoligotype strains are recorded as spoligotype-international-type (SIT) 2669 in the SITVIT international database, this is the first time to our knowledge that it has been shown that negative spoligotype strains are found in at least 4 different 24 loci MIRU-VNTR and RD deletion families. We report for the first time negative spoligotype-associated total loss of CRISPR region spacers and repeats, with accompanying clonally specific loss of flanking genes, including at least CRISPR-associated genes cas2 and cas1. Since cas1 deleted E.coli shows increased sensitivity to DNA damage and impaired chromosomal segregation, we discussed the possibility of a similar phenotype in the deleted-DR strains and Beijing family strains as both lack the cas1 gene.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3560-6) contains supplementary material, which is available to authorized users.