Data-driven model for the assessment of Mycobacterium tuberculosis transmission in evolving demographic structures.

In the case of tuberculosis (TB), the capabilities of epidemic models to produce quantitatively robust forecasts are limited by multiple hindrances. Among these, understanding the complex relationship between disease epidemiology and populations' age structure has been highlighted as one of the most relevant. TB dynamics depends on age in multiple ways, some of which are traditionally simplified in the literature. That is the case of the heterogeneities in contact intensity among different age strata that are common to all airborne diseases, but still typically neglected in the TB case. Furthermore, while demographic structures of many countries are rapidly aging, demographic dynamics are pervasively ignored when modeling TB spreading. In this work, we present a TB transmission model that incorporates country-specific demographic prospects and empirical contact data around a data-driven description of TB dynamics. Using our model, we find that the inclusion of demographic dynamics is followed by an increase in the burden levels predicted for the next decades in the areas of the world that are most hit by the disease today. Similarly, we show that considering realistic patterns of contacts among individuals in different age strata reshapes the transmission patterns reproduced by the models, a result with potential implications for the design of age-focused epidemiological interventions.

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.

Rapid diagnosis of tuberculosis. Detection of drug resistance mechanisms. / Diagnóstico rápido de la tuberculosis. Detección de mecanismos de resistencia.

Tuberculosis is still a serious public health problem, with 10.8 million new cases and 1.8 million deaths worldwide in 2015. The diversity among members of the Mycobacterium tuberculosis complex, the causal agent of tuberculosis, is conducive to the design of different methods for rapid diagnosis. Mutations in the genes involved in resistance mechanisms enable the bacteria to elude the treatment. We have reviewed the methods for the rapid diagnosis of M. tuberculosis complex and the detection of susceptibility to drugs, both of which are necessary to prevent the onset of new resistance and to establish early, appropriate treatment.

In Vivo IS6110 Profile Changes in a Mycobacterium tuberculosis Strain as Determined by Tracking over 14 Years.

Transposition and homologous recombination of IS6110 appear in Mycobacterium tuberculosis along in vivo sequential infections. These events were checked in different clones of a successful strain, M. tuberculosis Zaragoza, with the focus on a variant in which integration of a copy of IS6110 in the origin of replication (oriC) region occurred.

A specific polymorphism in Mycobacterium tuberculosis H37Rv causes differential ESAT-6 expression and identifies WhiB6 as a novel ESX-1 component.

The ESX-1 secreted virulence factor ESAT-6 is one of the major and most well-studied virulence factors of Mycobacterium tuberculosis, given that its inactivation severely attenuates virulent mycobacteria. In this work, we show that clinical isolates of M. tuberculosis produce and secrete larger amounts of ESAT-6 than the widely used M. tuberculosis H37Rv laboratory strain. A search for the genetic polymorphisms underlying this observation showed that whiB6 (rv3862c), a gene upstream of the ESX-1 genetic locus that has not previously been found to be implicated in the regulation of the ESX-1 secretory apparatus, presents a unique single nucleotide insertion in its promoter region in strains H37Rv and H37Ra. This polymorphism is not present in any of the other publicly available M. tuberculosis complex genomes or in any of the 76 clinical M. tuberculosis isolates analyzed in our laboratory. We demonstrate that in consequence, the virulence master regulator PhoP downregulates whiB6 expression in H37Rv, while it upregulates its expression in clinical strains. Importantly, reintroduction of the wild-type (WT) copy of whiB6 in H37Rv restored ESAT-6 production and secretion to the level of clinical strains. Hence, we provide clear evidence that in M. tuberculosis--with the exception of the H37Rv strain--ESX-1 expression is regulated by WhiB6 as part of the PhoP regulon, which adds another level of complexity to the regulation of ESAT-6 secretion with a potential role in virulence adaptation.

Proposal of a consensus set of hypervariable mycobacterial interspersed repetitive-unit-variable-number tandem-repeat loci for subtyping of Mycobacterium tuberculosis Beijing isolates.

Mycobacterium tuberculosis Beijing strains represent targets of special importance for molecular surveillance of tuberculosis (TB), especially because they are associated with spread of multidrug resistance in some world regions. Standard 24-locus mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing lacks resolution power for accurately discriminating closely related clones that often compose Beijing strain populations. Therefore, we evaluated a set of 7 additional, hypervariable MIRU-VNTR loci for better resolution and tracing of such strains, using a collection of 535 Beijing isolates from six world regions where these strains are known to be prevalent. The typeability and interlaboratory reproducibility of these hypervariable loci were lower than those of the 24 standard loci. Three loci (2163a, 3155, and 3336) were excluded because of their redundant variability and/or more frequent noninterpretable results compared to the 4 other markers. The use of the remaining 4-locus set (1982, 3232, 3820, and 4120) increased the number of types by 52% (from 223 to 340) and reduced the clustering rate from 58.3 to 36.6%, when combined with the use of the standard 24-locus set. Known major clonal complexes/24-locus-based clusters were all subdivided, although the degree of subdivision varied depending on the complex. Only five single-locus variations were detected among the hypervariable loci of an additional panel of 92 isolates, representing 15 years of clonal spread of a single Beijing strain in a geographically restricted setting. On this calibrated basis, we propose this 4-locus set as a consensus for subtyping Beijing clonal complexes and clusters, after standard typing.

Single nucleotide polymorphism (SNP) analysis used for the phylogeny of the Mycobacterium tuberculosis complex based on a pyrosequencing assay.

BACKGROUND: Different polymorphisms have been described as markers to classify the lineages of the Mycobacterium tuberculosis complex. The analysis of nine single nucleotide polymorphisms (SNPs) was used to describe seven SNPs cluster groups (SCGs). We attempted to classify those strains that could not been categorized into lineages by the genotyping methods used in the routine testing. RESULTS: The M. tuberculosis complex isolates collected in 2010 in our region were analysed. A new method based on multiplex-PCRs and pyrosequencing to analyse these SNPs was designed. For the pyrosequencing assay nine SNPs that defined the seven SCGs were selected from the literature: 1977, 74092, 105139, 232574, 311613, 913274, 2460626, 3352929 and gyrA95. In addition, SNPs in katG(463), mgtC(182), Ag85C(103) and RD(Rio) deletion were detected. CONCLUSIONS: This work has permitted to achieve a better classification of Aragonian strains into SCGs and in some cases, to assign strains to its certain lineage. Besides, the description of a new pattern shared by two isolates "SCG-6c" reinforces the interest of SNPs to follow the evolution of M. tuberculosis complex.

Role of wild bird and rodents in the epidemiology of subclinical salmonellosis in finishing pigs.

Wild birds and rodents may play an important role in the dynamics of subclinical pig salmonellosis, either as the introducers of the bacteria into the farm or as receptors of an infection already established in the farm. We tried to gain further insight into the epidemiology of this infection by studying the phenotypic (i.e., serotype and antimicrobial resistance patterns) and molecular characteristics of Salmonella strains isolated from samples collected from pigs and wildlife captured in the vicinity of pig farms. Salmonella-positive pig fecal samples were identified in 56.1% of the 41 farms investigated. Birds shedding Salmonella spp. were detected in 21.4% of the farms despite the low numbers of birds captured in many farms. Most Salmonella isolates from birds (74%) did not show any antimicrobial resistance (AR) pattern and belonged to phage types rarely seen in the pig population (U310, DT56, DT137, DT164), supporting the likely avian source of infection for most birds. The proportion of farms showing Salmonella-infected rodents was higher (46.2%), with Salmonella isolates showing a high homology with those likely originated from pigs. Salmonella-positive environmental samples were found in >50% of the farms, and the characteristics of these Salmonella strains supported the idea of pigs as a major source of Salmonella contamination of the farm environment. Dissemination of Salmonella in pig farms from areas of high Salmonella prevalence appeared to depend to some extent upon rodents and wild birds present in the farm, but the role of rodents in its maintenance seemed to be somewhat more relevant than that of birds. In conclusion, activities aimed at reducing the contact of these wild species with pigs will probably assist in the control of pig salmonellosis. Strict hygienic measures should be considered in areas of high prevalence of infection to lower the high load of environmental contamination.

Mapping IS6110 in high-copy number Mycobacterium tuberculosis strains shows specific insertion points in the Beijing genotype.

BACKGROUND: Mycobacterium tuberculosis Beijing strains are characterized by a large number of IS6110 copies, suggesting the potential implication of this element in the virulence and capacity for rapid dissemination characteristic of this family. This work studies the insetion points of IS6110 in high-copy clinical isolates specifically focusing on the Beijing genotype. RESULTS: In the present work we mapped the insertion points of IS6110 in all the Beijing strains available in the literature and in the DNA sequence databases. We generated a representative primer collection of the IS6110 locations, which was used to analyse 61 high-copy clinical isolates. A total of 440 points of insertion were identified and analysis of their flanking regions determined the exact location, the direct repeats (DRs), the orientation and the distance to neighboring genes of each copy of IS6110. We identified specific points of insertion in Beijing strains that enabled us to obtain a dendrogram that groups the Beijing genotype. CONCLUSIONS: This work presents a detailed analysis of locations of IS6110 in high-copy clinical isolates, showing points of insertion present with high frequency in the Beijing family and absent in other strains.

Analysis of mutations in streptomycin-resistant strains reveals a simple and reliable genetic marker for identification of the Mycobacterium tuberculosis Beijing genotype.

The Mycobacterium tuberculosis pandemic is a major health problem, further complicated by an increasing incidence of drug-resistant isolates and the existence of highly transmissible strains, such as those in the Beijing family. Streptomycin (STR)-resistant M. tuberculosis clinical isolates have been analyzed to look for mutations in the rpsL, rrs, and gidB genes. In addition, the Rv1258c gene, which encodes Tap, an efflux pump that transports STR, has been sequenced. Mutations affecting codons 43 and 88 of the rpsL gene were found in 44.4% of the strains, and 16.7% of the strains carried mutations in the rrs gene, both of which probably contribute to STR resistance. Many strains presented with mutations in the gidB gene, but the implication of those mutations in STR resistance remains unclear. Interestingly, a cytosine nucleotide insertion between positions 580 and 581 (denominated Tap(580)) in the Rv1258c gene has been found in all Beijing isolates included in this study, suggesting that it might be a novel polymorphism specific to the Beijing family of M. tuberculosis. A simple and fast restriction fragment length polymorphism (RFLP)-PCR method for detecting the Tap(580) insertion has been developed and used to screen a collection of 220 DNA samples obtained from cultures of M. tuberculosis isolates and 30 respiratory specimens. In all cases, the Beijing and non-Beijing representative samples were identified correctly. Tap(580) is a novel polymorphism specific to the highly transmissible Beijing family, which allows for fast detection of these strains even at the very early stages of infection.

Global study of IS6110 in a successful Mycobacterium tuberculosis strain: clues for deciphering its behavior and for its rapid detection.

The Mycobacterium tuberculosis insertion sequence IS6110, besides being a very useful tool in molecular epidemiology, seems to have an impact on the biology of bacilli. In the present work, we mapped the 12 points of insertion of IS6110 in the genome of a successful strain named M. tuberculosis Zaragoza (which has been referred to as the MTZ strain). This strain, belonging to principal genetic group 3, caused a large unsuspected tuberculosis outbreak involving 85 patients in Zaragoza, Spain, in 2001 to 2004. The mapping of the points of insertion of IS6110 in the genome of the Zaragoza strain offers clues for a better understanding of the adaptability and virulence of M. tuberculosis. Surprisingly, the presence of one copy of IS6110 was found in Rv2286c, as was recently described for a successful Beijing sublineage. As a result of this analysis, a rapid method for detecting this particular M. tuberculosis strain has been designed.

Transcriptomic profile of the most successful Mycobacterium tuberculosis strain in Aragon, the MtZ strain, during exponential and stationary growth phases.

IMPORTANCE: Aragon Community suffered, during the first years of the beginning of this century, a large outbreak caused by the MtZ strain, producing more than 240 cases to date. MtZ strain and the outbreak have been previously studied from an epidemiological and molecular point of view. In this work, we analyzed the transcriptomic profile of the strain for better understanding of its success among our population. We have discovered that MtZ has some upregulated virulence pathways, such as the ESX-1 system, the cholesterol degradation pathway or the peptidoglycan biosynthesis. Interestingly, MtZ has downregulated the uptake of iron. Another special feature of MtZ strain is the interruption of desA3 gene, essential for producing oleic acid. Although the strain takes a long time to grow in the initial culture media, eventually it is able to reach normal optical densities, suggestive of the presence of another route for obtaining oleic acid in Mycobacterium tuberculosis.

Rapid Identification of Lineage and Drug Resistance in Clinical Samples of Mycobacterium tuberculosis.

BACKGROUND: Mycobacterium tuberculosis is a slow-growing bacterium, which could delay its diagnosis and, therefore, promote the spread of the disease. Whole-genome sequencing allows us to obtain the complete drug-resistance profile of the strain; however, bacterial cultivation of clinical samples, along with complex processing, is required. METHODS: In this work, we explore AmpliSeq, an amplicon-based enrichment method for preparing libraries for targeted next-generation sequencing, to identify lineage and drug resistance directly from clinical samples. RESULTS: In our study, 111 clinical samples were tested. The lineage was identified in 100% of the culture-derived samples (52/52), in 95% of the smear (BK)-positive clinical samples (38/40) and in 42.1% of the BK-negative clinical samples (8/19). The drug-resistance profile was accurately identified in all but 11 samples, in which some phenotypic and genotypic discrepancies were found. In this respect, our panels were not exact in the detection of streptomycin resistance for isolates derived from clinical samples, as an extremely high number of SNPs in the rrs and rrl genes were detected due to cross-contamination. CONCLUSION: This technique has demonstrated high sensitivity in obtaining the drug-resistance profile of the isolates, as even those samples with DNA concentrations below the detection limit of Qubit produced a result. AmpliSeq technology is cheaper than whole-genome sequencing, easy to perform by laboratory technicians and applicable to any microorganism using the Ion Torrent platform.

Deciphering the Tangible Spatio-Temporal Spread of a 25-Year Tuberculosis Outbreak Boosted by Social Determinants.

Outbreak strains of Mycobacterium tuberculosis are promising candidates as targets in the search for intrinsic determinants of transmissibility, as they are responsible for many cases with sustained transmission; however, the use of low-resolution typing methods and restricted geographical investigations represent flaws in assessing the success of long-lived outbreak strains. We can now address the nature of outbreak strains by combining large genomic data sets and phylodynamic approaches. We retrospectively sequenced the whole genome of representative samples assigned to an outbreak circulating in the Canary Islands (the GC strain) since 1993, which accounts for ~20% of local tuberculosis cases. We selected a panel of specific single nucleotide polymorphism (SNP) markers for an in-silico search for additional outbreak-related sequences within publicly available tuberculosis genomic data. Using this information, we inferred the origin, spread, and epidemiological parameters of the GC strain. Our approach allowed us to accurately trace the historical and more recent dispersion of the GC strain. We provide evidence of a highly successful nature within the Canarian archipelago but limited expansion abroad. Estimation of epidemiological parameters from genomic data disagree with a distinctive biology of the GC strain. With the increasing availability of genomic data allowing for the accurate inference of strain spread and critical epidemiological parameters, we can now revisit the link between Mycobacterium tuberculosis genotypes and transmission, as is routinely carried out for SARS-CoV-2 variants of concern. We demonstrate that social determinants rather than intrinsically higher bacterial transmissibility better explain the success of the GC strain. Importantly, our approach can be used to trace and characterize strains of interest worldwide. IMPORTANCE Infectious disease outbreaks represent a significant problem for public health. Tracing outbreak expansion and understanding the main factors behind emergence and persistence remain critical to effective disease control. Our study allows researchers and public health authorities to use Whole-Genome Sequencing-based methods to trace outbreaks, and shows how available epidemiological information helps to evaluate the factors underpinning outbreak persistence. Taking advantage of all the freely available information placed in public repositories, researchers can accurately establish the expansion of an outbreak beyond original boundaries, and determine the potential risk of a strain to inform health authorities which, in turn, can define target strategies to mitigate expansion and persistence. Finally, we show the need to evaluate strain transmissibility in different geographic contexts to unequivocally associate spread to local or pathogenic factors, an important lesson taken from genomic surveillance of SARS-CoV-2.

IS-seq: a novel high throughput survey of in vivo IS6110 transposition in multiple Mycobacterium tuberculosis genomes.

BACKGROUND: The insertion element IS6110 is one of the main sources of genomic variability in Mycobacterium tuberculosis, the etiological agent of human tuberculosis. Although IS 6110 has been used extensively as an epidemiological marker, the identification of the precise chromosomal insertion sites has been limited by technical challenges. Here, we present IS-seq, a novel method that combines high-throughput sequencing using Illumina technology with efficient combinatorial sample multiplexing to simultaneously probe 519 clinical isolates, identifying almost all the flanking regions of the element in a single experiment. RESULTS: We identified a total of 6,976 IS6110 flanking regions on the different isolates. When validated using reference strains, the method had 100% specificity and 98% positive predictive value. The insertions mapped to both coding and non-coding regions, and in some cases interrupted genes thought to be essential for virulence or in vitro growth. Strains were classified into families using insertion sites, and high agreement with previous studies was observed. CONCLUSIONS: This high-throughput IS-seq method, which can also be used to map insertions in other organisms, extends previous surveys of in vivo interrupted loci and provides a baseline for probing the consequences of disruptions in M. tuberculosis strains.

Rapid test for identification of a highly transmissible Mycobacterium tuberculosis Beijing strain of sub-Saharan origin.

The development of a rapid test to identify Mycobacterium tuberculosis Beijing isolates and specifically strain GC1237, coming from a sub-Saharan country, is needed due to its alarming wide spread on Gran Canaria Island (Spain). A rapid test that detects IS6110 present between dnaA and dnaN in the Beijing strains and in a specific site for GC1237 (Rv2180c) has been developed. This test would be a useful tool in the surveillance of subsequent cases.

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.

Estimation of the mutation rate of Mycobacterium tuberculosis in cases with recurrent tuberculosis using whole genome sequencing.

The study of tuberculosis latency is problematic due to the difficulty of isolating the bacteria in the dormancy state. Despite this, several in vivo approaches have been taken to mimic the latency process. Our group has studied the evolution of the bacteria in 18 cases of recurrent tuberculosis. We found that HIV positive patients develop recurrent tuberculosis earlier, generally in the first two years (p value = 0.041). The genome of the 36 Mycobacterium tuberculosis paired isolates (first and relapsed isolates) showed that none of the SNPs found within each pair was observed more than once, indicating that they were not directly related to the recurrence process. Moreover, some IS6110 movements were found in the paired isolates, indicating the presence of different clones within the patient. Finally, our results suggest that the mutation rate remains constant during all the period as no correlation was found between the number of SNPs and the time to relapse.

Analysis of the twenty-six largest outbreaks of tuberculosis in Aragon using whole-genome sequencing for surveillance purposes.

The incidence of tuberculosis in Aragon, Spain, is around ten cases per 100,000 inhabitants. Since 2004, a molecular surveillance protocol has been carried out; therefore, all M. tuberculosis strains are genotyped. Recently, whole-genome sequencing has been implemented for relevant isolates. The aim of this work is to characterise at the molecular level the causative strains of the 26 largest outbreaks of the community (including ten or more cases), genotyped by IS6110-RFLP and causing 26% of tuberculosis cases. To achieve this objective, two or three isolates of each IS6110-cluster belonging to different years were selected for sequencing. We found that strains of lineages L4.8, L4.3 and L4.1.2 were the most frequent. The threshold of 12 SNPs as the maximum distance for confirming the belonging to an outbreak was met for 18 of the 26 IS6110-clusters. Four pairs of isolates with more than 90 SNPs were identified as not belonging to the same strain, and four other pairs were kept in doubt as the number of SNPs was close to 12, between 14 and 35. The study of Regions of Difference revealed that they are lineage conserved. Moreover, we could analyse the IS6110 locations for all genome-sequenced isolates, finding some frequent locations in isolates belonging to the same lineage and certain IS6110 movements between the paired isolates. In the vast majority, these movements were not captured by the IS6110-RFLP pattern. After classifying the genes containing SNP by their functional category, we could confirm that the number of SNPs detected in genes considered as virulence factors and the number of cases the strain produced were not related, suggesting that a particular SNP is more relevant than the number. The characteristics found in the most successful strains in our community could be useful for other researchers in epidemiology, virulence and pathogenesis.

Successful control of Serratia marcescens outbreak in a neonatal unit of a tertiary-care hospital in Spain.

OBJECTIVE: Serratia marcescens is a Gram-negative bacterium that is found in hospital environments and commonly associated with outbreaks in neonatal units. One S. marcescens isolate was detected from a bloodstream culture from a neonate in our hospital that was followed by an outbreak. The aim of this study was to describe the molecular epidemiology of a S. marcescens outbreak in the neonatal unit. METHODS: In order to investigate the outbreak, weekly surveillance rectal swabs were submitted for culture from all patients admitted in this unit from August to September 2018. Environmental samples were obtained from potential sources in September 2018. Typing of isolates was performed by pulsed field gel electrophoresis (PFGE). In addition, we studied the in vitro activity of chlorhexidine against S. marcescens. RESULTS: During this period, 146 infants were hospitalised in our neonatal unit, of which 16 patients had a S. marcescens-positive sample. A total of 36 environmental surveillance samples were collected, and one sample from a stethoscope from an incubator of a colonized baby was positive for S. marcescens. All the 18 isolates, including the isolate from the stethoscope, belonged to a single PFGE cluster. We found that very low concentrations of chlorhexidine, even with application times close to 0 achieved significant reductions in the amount of S. marcescens. CONCLUSION: A unique clone of S. marcescens caused this outbreak, including isolates from patients and from one stethoscope. The outbreak was controlled with the early implementation of specific control measures.