Epidemiological and clinical insights into the enterovirus D68 upsurge in Europe 2021/22 and the emergence of novel B3-derived lineages, ENPEN multicentre study.

Enterovirus D68 (EV-D68) infections are associated with severe respiratory disease and acute flaccid myelitis (AFM). The European Non-Polio Enterovirus Network (ENPEN) aimed to investigate the epidemiological and genetic characteristics of EV-D68 and its clinical impact during the fall-winter season of 2021/22. From 19 European countries, 58 institutes reported 10,481 (6.8%) EV-positive samples of which 1,004 (9.6%) were identified as EV-D68 (852 respiratory samples). Clinical data was reported for 969 cases. 78.9% of infections were reported in children (0-5 years); 37.9% of cases were hospitalised. Acute respiratory distress was commonly noted (93.1%) followed by fever (49.4%). Neurological problems were observed in 6.4% of cases with six reported with AFM. Phylodynamic/Nextstrain and phylogenetic analyses based on 694 sequences showed the emergence of two novel B3-derived lineages, with no regional clustering. In conclusion, we describe a large-scale EV-D68 European upsurge with severe clinical impact and the emergence of B3-derived lineages.

Host DNA depletion can increase the sensitivity of Mycobacterium spp. detection through shotgun metagenomics in sputum.

Identification and phenotypic drug-susceptibility testing for mycobacteria are time-consuming and challenging but essential for managing mycobacterial infections. Next-generation sequencing (NGS) technologies can increase diagnostic speed and quality, but standardization is still lacking for many aspects (e.g., unbiased extraction, host depletion, bioinformatic analysis). Targeted PCR approaches directly on sample material are limited by the number of targets that can be included. Unbiased shotgun metagenomics on direct material is hampered by the massive amount of host DNA, which should be removed to improve the microbial detection sensitivity. For this reason, we developed a method for NGS-based diagnosis of mycobacteria directly from patient material. As a model, we used the non-tuberculous mycobacterium (NTM) Mycobacterium abscessus. We first compared the efficiency of three different DNA extraction kits for isolating DNA (quality and concentration). The two most efficient kits were then used in a follow-up study using artificial sputum. Finally, one extraction kit was selected and further evaluated for DNA isolation from a patients' sputum mixture spiked with M. abscessus at three concentrations (final concentrations 108, 107, 106 CFU/ml). The spiked sputum samples were processed with and without saponin treatment (ST) in combination with DNAse treatment prior to bacterial DNA extraction to evaluate the recovery of bacteria and depletion of host DNA by PCR and Illumina sequencing. While Ct values of the qPCR targeting mycobacterial ITS DNA remained rather stable, Ct values in the qPCR targeting the human ß-actin gene increased by five Ct values in ST samples. In subsequent Illumina sequencing, a decrease of 89% of reads mapped to the human genome was observed in ST samples. The percentage of reads mapped to M. abscessus (108 CFU/ml) increased by 89%, and the sequencing depth increased two times when undergoing ST. In conclusion, the sensitivity of M. abscessus detection in artificial sputum was increased using a saponin pre-treatment step. The saponin followed by the DNase I treatment approach could be efficiently applied to detect and characterize mycobacterial infections, including tuberculosis, directly from sputum.

Genetic profiling of Mycobacterium tuberculosis revealed "modern" Beijing strains linked to MDR-TB from Southwestern Colombia.

Beijing strains of Mycobacterium tuberculosis (lineage 2) have been associated with drug-resistance and transmission of tuberculosis worldwide. Most of the Beijing strains identified in the Colombian Pacific coast have exhibited a multidrug resistant (MDR) phenotype. We sought to evaluate the clonality and sublineage of Beijing strains circulating in Southwestern Colombia. Thirty-seven Beijing strains were identified through spoligotyping out of 311 clinical isolates collected in 9 years from 2002-2010. Further analysis by MIRU-VNTR 24 loci was conducted for the Beijing strains. For sublineage classification, deletions of RD105, RD207, and RD131 and point mutations at fbpB, mutT2, and acs were evaluated. Drug-resistance associated mutations to first- and second-line anti-TB drugs were also evaluated. Additionally, two Beijing strains were Illumina-whole genome sequenced (one MDR and one drug-susceptible). Among the 37 Beijing strains characterized, 36 belonged to the SIT190 type from which 28 were MDR, four pre-extensively drug resistant (XDR) TB, and four XDR-TB. The remaining strain was SIT1 and drug susceptible. MIRU-VNTR analysis allowed the identification of three Beijing clusters and two unique strains. Beijing strains were confirmed as "modern" sublineage. The mutations rpoB S531L and katG S315T were the most common among MDR strains. Moreover, the two strains evaluated by whole genome sequencing (WGS) shared most of the genetic features with the sublineage 2.2.1 "modern" Beijing previously characterized from Asian strains. WGS analysis of the MDR strain revealed the presence of eight SNPs previously reported in other MDR "Beijing-like" strains from Colombia. The presence of "modern" Beijing strains in Southwestern Colombia, most of them with MDR phenotype, suggests a different origin of this M. tuberculosis sublineage compared to other Beijing strains found in neighboring South American countries. This work may serve as a genetic baseline to study the evolution and spread of M. tuberculosis Beijing strains in Colombia, which play an important role in the propagation of MDR-TB.

Pyrazinamide resistance-conferring mutations in pncA and the transmission of multidrug resistant TB in Georgia.

BACKGROUND: The ongoing epidemic of multidrug-resistant tuberculosis (MDR-TB) in Georgia highlights the need for more effective control strategies. A new regimen to treat MDR-TB that includes pyrazinamide (PZA) is currently being evaluated and PZA resistance status will largely influence the success of current and future treatment strategies. PZA susceptibility testing was not routinely performed at the National Reference Laboratory (NRL) in Tbilisi between 2010 and September 2015. We here provide a first insight into the prevalence of PZA resistant TB in this region. METHODS: Phenotypic susceptibility to PZA was determined in a convenience collection of well-characterised TB patient isolates collected at the NRL in Tbilisi between 2012 and 2013. In addition, the pncA gene was sequenced and whole genome sequencing was performed on two isolates. RESULTS: Out of 57 isolates tested 33 (57.9%) showed phenotypic drug resistance to PZA and had a single pncA mutation. All of these 33 isolates were MDR-TB strains. pncA mutations were absent in all but one of the 24 PZA susceptible isolate. In total we found 18 polymorphisms in the pncA gene. From the two major MDR-TB clusters represented (94-32 and 100-32), 10 of 15, 67.0% and 13 of 14, 93.0% strains, respectively were PZA resistant. We also identified a member of the potentially highly transmissive clade A strain carrying the characteristic I6L substitution in PncA. Another strain with the same MLVA type as the clade A strain acquired a different mutation in pncA and was genetically more distantly related suggesting that different branches of this particular lineage have been introduced into this region. CONCLUSION: In this high MDR-TB setting more than half of the tested MDR-TB isolates were resistant to PZA. As PZA is part of current and planned MDR-TB treatment regimens this is alarming and deserves the attention of health authorities. Based on our typing and sequence analysis results we conclude that PZA resistance is the result of primary transmission as well as acquisition within the patient and recommend prospective genotyping and PZA resistance testing in high MDR-TB settings. This is of utmost importance in order to preserve bacterial susceptibility to PZA to help protect (new) second line drugs in PZA containing regimens.

Quick and cheap MIRU-VNTR typing of Mycobacterium tuberculosis species complex using duplex PCR.

While minisatellites are usually typed using capillary sequencers or qiaplex systems in developed countries, many low-resource regions cannot afford it. We propose an optimized agarose gel electrophoresis method to genotype Mycobacterium tuberculosis species complex minisatellites in their standardized format (24 MIRU-VNTR). It is based on duplex PCRs combining VNTR loci harboring distinct amplicon sizes whatever the repetition number of each locus. This method performs well both on DNA extracts of good quality and on thermolysates while reducing workload and reagents costs.

Mechanisms of Phenotypic Rifampicin Tolerance in Mycobacterium tuberculosis Beijing Genotype Strain B0/W148 Revealed by Proteomics.

The "successful" Russian clone B0/W148 of Mycobacterium tuberculosis Beijing is well-known for its capacity to develop antibiotic resistance. During treatment, resistant mutants can occur that have inheritable resistance to specific antibiotics. Next to mutations, M. tuberculosis has several mechanisms that increase their tolerance to a variety of antibiotics. Insights in the phenotypic mechanisms that contribute to drug tolerance will increase our understanding of how antibiotic resistance develops in M. tuberculosis. In this study, we examined the (phospho)proteome dynamics in M. tuberculosis Beijing strain B0/W148 when exposed to a high dose of rifampicin; one of the most potent first-line antibiotics. A total of 2,534 proteins and 191 phosphorylation sites were identified, and revealed the differential regulation of DosR regulon proteins, which are necessary for the development of a dormant phenotype that is less susceptible to antibiotics. By examining independent phenotypic markers of dormancy, we show that persisters of in vitro rifampicin exposure entered a metabolically hypoactive state, which yields rifampicin and other antibiotics largely ineffective. These new insights in the role of protein regulation and post-translational modifications during the initial phase of rifampicin treatment reveal a shortcoming in the antituberculosis regimen that is administered to 8-9 million individuals annually.

Predominance of modern Mycobacterium tuberculosis strains and active transmission of Beijing sublineage in Jayapura, Indonesia Papua.

Mycobacterium tuberculosis genotype distribution is different between West and Central Indonesia, but there are no data on the most Eastern part, Papua. We aimed to identify the predominant genotypes of M. tuberculosis responsible for tuberculosis in coastal Papua, their transmission, and the association with patient characteristics. A total of 199 M. tuberculosis isolates were collected. Spoligotyping was applied to describe the population structure of M. tuberculosis, lineage identification was performed using a combination of lineage-specific markers, and genotypic clusters were identified using a combination of 24-locus-MIRU-VNTR and spoligotyping. A high degree of genetic diversity was observed among isolates based on their spoligopatterns. Strains from modern lineage 4 made up almost half of strains (46.9%), being more abundant than the ancient lineage 1 (33.7%), and modern lineage 2 (19.4%). Thirty-five percent of strains belonged to genotypic clusters, especially strains in the Beijing genotype. Previous TB treatment and mutations associated with drug resistance were more common in patients infected with strains of the Beijing genotype. Papua shows a different distribution of M. tuberculosis genotypes compared to other parts of Indonesia. Clustering and drug resistance of modern strains recently introduced to Papua may contribute to the high tuberculosis burden in this region.

Genomic characterization of Mycobacterium tuberculosis lineage 7 and a proposed name: 'Aethiops vetus'.

Lineage 7 of the Mycobacterium tuberculosis complex has recently been identified among strains originating from Ethiopia. Using different DNA typing techniques, this study provides additional information on the genetic heterogeneity of five lineage 7 strains collected in the Amhara Region of Ethiopia. It also confirms the phylogenetic positioning of these strains between the ancient lineage 1 and TbD1-deleted, modern lineages 2, 3 and 4 of Mycobacterium tuberculosis. Four newly identified large sequence polymorphisms characteristic of the Amhara Region lineage 7 strains are described. While lineage 7 strains have been previously identified in the Woldiya area, we show that lineage 7 strains circulate in other parts of the Amhara Region and also among foreign-born individuals from Eritrea and Somalia in The Netherlands. For ease of documenting future identification of these strains in other geographical locations and recognizing the place of origin, we propose to assign lineage 7 strains the lineage name 'Aethiops vetus'.

Optimizing multiplex SNP-based data analysis for genotyping of Mycobacterium tuberculosis isolates.

BACKGROUND: Multiplex ligation-dependent probe amplification (MLPA) is a powerful tool to identify genomic polymorphisms. We have previously developed a single nucleotide polymorphism (SNP) and large sequence polymorphisms (LSP)-based MLPA assay using a read out on a liquid bead array to screen for 47 genetic markers in the Mycobacterium tuberculosis genome. In our assay we obtain information regarding the Mycobacterium tuberculosis lineage and drug resistance simultaneously. Previously we called the presence or absence of a genotypic marker based on a threshold signal level. Here we present a more elaborate data analysis method to standardize and streamline the interpretation of data generated by MLPA. The new data analysis method also identifies intermediate signals in addition to classification of signals as positive and negative. Intermediate calls can be informative with respect to identifying the simultaneous presence of sensitive and resistant alleles or infection with multiple different Mycobacterium tuberculosis strains. RESULTS: To validate our analysis method 100 DNA isolates of Mycobacterium tuberculosis extracted from cultured patient material collected at the National TB Reference Laboratory of the National Center for Tuberculosis and Lung Diseases in Tbilisi, Republic of Georgia were tested by MLPA. The data generated were interpreted blindly and then compared to results obtained by reference methods. MLPA profiles containing intermediate calls are flagged for expert review whereas the majority of profiles, not containing intermediate calls, were called automatically. No intermediate signals were identified in 74/100 isolates and in the remaining 26 isolates at least one genetic marker produced an intermediate signal. CONCLUSION: Based on excellent agreement with the reference methods we conclude that the new data analysis method performed well. The streamlined data processing and standardized data interpretation allows the comparison of the Mycobacterium tuberculosis MLPA results between different experiments. All together this will facilitate the implementation of the MLPA assay in different settings.

Optimization of standard in-house 24-locus variable-number tandem-repeat typing for Mycobacterium tuberculosis and its direct application to clinical material.

Variable-number tandem-repeat (VNTR) typing with a panel of 24 loci is the current gold standard in the molecular typing of Mycobacterium tuberculosis complex isolates. However, because of technical problems, a part of the loci often cannot be amplified by multiplex PCRs. Therefore, a considerable number of single-locus PCRs have to be performed for the loci with missing results, which impairs the laboratory work flow. Therefore, the original in-house method described by Supply et al. in 2006 was reevaluated. We modified seven primers and the PCR master mixture and obtained a strongly optimized in-house 24-locus VNTR typing method. The percentage of instantly complete 24-locus VNTR patterns detected in the routine flow of typing activities increased to 84.7% from the 72.3% obtained with the typing conducted with the commercially available Genoscreen MIRU-VNTR typing kit. The analytical sensitivity of the optimized in-house method was assessed by serial dilutions of M. tuberculosis in bronchoalveolar lavage fluid. A 1:10 dilution of the different strains tested was the lowest dilution for the detection of a complete 24-locus VNTR pattern. The optimized in-house 24-locus VNTR typing method will reduce the turnaround time of typing significantly and also the financial burden of these activities.

Comparison of 14 molecular assays for detection of Mycobacterium tuberculosis complex in bronchoalveolar lavage fluid.

We compared 14 molecular assays for their ability to detect the Mycobacterium tuberculosis complex in bronchoalveolar lavage fluid samples. Three approaches were followed. First, by using DNA from Mycobacterium bovis BCG, we determined the detection limits of the assays using routine molecular methods. Second, in order to determine the analytical sensitivities of the assays, we added one of four M. tuberculosis isolates with various numbers of the insertion sequence IS6110 to N-acetyl-l-cysteine (NALC)-NaOH-treated bronchoalveolar lavage fluid samples in dilutions of 1:10 to 1:10,000,000. Third, intertest variabilities were measured and defined by the standard deviations for the quantitation cycle (Cq) values of three positive test results per dilution per assay. The 14 assays tested had similar analytical sensitivities, except for GeneXpert, which had an analytical sensitivity that was 10- to 100-fold lower than that of the other assays. The MP MTB/NTM test and the in-house TaqMan-10 revealed the best performances for the detection limit and had the highest analytical sensitivities. Most of the tests performed well regarding detection limit and analytical sensitivity for the detection of the M. tuberculosis complex in serial dilutions, and the differences were small. The MP MTB/NTM and the in-house TaqMan-10 assays revealed the best, and GeneXpert the worst, overall performances.

Clustering of tuberculosis cases based on variable-number tandem-repeat typing in relation to the population structure of Mycobacterium tuberculosis in the Netherlands.

The population structure of 3,776 Mycobacterium tuberculosis isolates was determined using variable-number tandem-repeat (VNTR) typing. The degree of clonality was so high that a more relaxed definition of clustering cannot be applied. Among recent immigrants with non-Euro-American isolates, transmission is overestimated if based on identical VNTR patterns.

Inferring patient to patient transmission of Mycobacterium tuberculosis from whole genome sequencing data.

BACKGROUND: Mycobacterium tuberculosis is characterised by limited genomic diversity, which makes the application of whole genome sequencing particularly attractive for clinical and epidemiological investigation. However, in order to confidently infer transmission events, an accurate knowledge of the rate of change in the genome over relevant timescales is required. METHODS: We attempted to estimate a molecular clock by sequencing 199 isolates from epidemiologically linked tuberculosis cases, collected in the Netherlands spanning almost 16 years. RESULTS: Multiple analyses support an average mutation rate of ~0.3 SNPs per genome per year. However, all analyses revealed a very high degree of variation around this mean, making the confirmation of links proposed by epidemiology, and inference of novel links, difficult. Despite this, in some cases, the phylogenetic context of other strains provided evidence supporting the confident exclusion of previously inferred epidemiological links. CONCLUSIONS: This in-depth analysis of the molecular clock revealed that it is slow and variable over short time scales, which limits its usefulness in transmission studies. However, the superior resolution of whole genome sequencing can provide the phylogenetic context to allow the confident exclusion of possible transmission events previously inferred via traditional DNA fingerprinting techniques and epidemiological cluster investigation. Despite the slow generation of variation even at the whole genome level we conclude that the investigation of tuberculosis transmission will benefit greatly from routine whole genome sequencing.

Comparative study of IS6110 restriction fragment length polymorphism and variable-number tandem-repeat typing of Mycobacterium tuberculosis isolates in the Netherlands, based on a 5-year nationwide survey.

In order to switch from IS6110 and polymorphic GC-rich repetitive sequence (PGRS) restriction fragment length polymorphism (RFLP) to 24-locus variable-number tandem-repeat (VNTR) typing of Mycobacterium tuberculosis complex isolates in the national tuberculosis control program in The Netherlands, a detailed evaluation on discriminatory power and agreement with findings in a cluster investigation was performed on 3,975 tuberculosis cases during the period of 2004 to 2008. The level of discrimination of the two typing methods did not differ substantially: RFLP typing yielded 2,733 distinct patterns compared to 2,607 in VNTR typing. The global concordance, defined as isolates labeled unique or identically distributed in clusters by both methods, amounted to 78.5% (n = 3,123). Of the remaining 855 cases, 12% (n = 479) of the cases were clustered only by VNTR, 7.7% (n = 305) only by RFLP typing, and 1.8% (n = 71) revealed different cluster compositions in the two approaches. A cluster investigation was performed for 87% (n = 1,462) of the cases clustered by RFLP. For the 740 cases with confirmed or presumed epidemiological links, 92% were concordant with VNTR typing. In contrast, only 64% of the 722 cases without an epidemiological link but clustered by RFLP typing were also clustered by VNTR typing. We conclude that VNTR typing has a discriminatory power equal to IS6110 RFLP typing but is in better agreement with findings in a cluster investigation performed on an RFLP-clustering-based cluster investigation. Both aspects make VNTR typing a suitable method for tuberculosis surveillance systems.

First worldwide proficiency study on variable-number tandem-repeat typing of Mycobacterium tuberculosis complex strains.

Although variable-number tandem-repeat (VNTR) typing has gained recognition as the new standard for the DNA fingerprinting of Mycobacterium tuberculosis complex (MTBC) isolates, external quality control programs have not yet been developed. Therefore, we organized the first multicenter proficiency study on 24-locus VNTR typing. Sets of 30 DNAs of MTBC strains, including 10 duplicate DNA samples, were distributed among 37 participating laboratories in 30 different countries worldwide. Twenty-four laboratories used an in-house-adapted method with fragment sizing by gel electrophoresis or an automated DNA analyzer, nine laboratories used a commercially available kit, and four laboratories used other methods. The intra- and interlaboratory reproducibilities of VNTR typing varied from 0% to 100%, with averages of 72% and 60%, respectively. Twenty of the 37 laboratories failed to amplify particular VNTR loci; if these missing results were ignored, the number of laboratories with 100% interlaboratory reproducibility increased from 1 to 5. The average interlaboratory reproducibility of VNTR typing using a commercial kit was better (88%) than that of in-house-adapted methods using a DNA analyzer (70%) or gel electrophoresis (50%). Eleven laboratories using in-house-adapted manual typing or automated typing scored inter- and intralaboratory reproducibilities of 80% or higher, which suggests that these approaches can be used in a reliable way. In conclusion, this first multicenter study has documented the worldwide quality of VNTR typing of MTBC strains and highlights the importance of international quality control to improve genotyping in the future.

Rapidly growing nontuberculous mycobacteria cultured from home tap and shower water.

Tap and shower water at two locations in the Netherlands was examined for the presence of rapidly growing nontuberculous mycobacteria. Cultures yielded Mycobacterium peregrinum, M. salmoniphilum, M. llatzerense, M. septicum, and three potentially novel species, a distribution different from that in clinical samples.

Characterization of Mycobacterium africanum subtype I among cows in a dairy farm in Bangladesh using spoligotyping.

Acid-fast bacilli (AFB) were detected in the autopsy lung tissue homogenate samples of four cows (variety Frisian cross) in a dairy farm in Bangladesh. Histopathological examination of the lung tissue demonstrated prominent granulomas, caseating necrosis and calcification indicative of tuberculosis (TB) infection. Mycobacteria could not be cultured from the tissue homogenate samples by Lowenstein-Jensen based conventional culture method though AFB were evident by Ziehl-Neelsen (ZN) staining of the smears of tissue homogenate and in paraffin embedded tissue slices. Spoligotyping performed on DNA extracts of paraffin embedded lung tissue samples confirmed the AFB as a member of the M. tuberculosis complex (MTBC) with a pattern assigned to M. africanum subtype I. This characterization by spoligotyping was confirmed by subjecting M. africanum subtype I isolates from other parts of the world to an alternative identification method based on DNA polymorphism in the gyrB gene (Hain Life Science, GmbH, Nehren, Germany). Since M. africanum is believed to be a human pathogen, general infection in cattle may be a public health threat. The presence of these bacteria in the animal reservoir most likely originated from a caretaker.

M. tuberculosis genotypic diversity and drug susceptibility pattern in HIV-infected and non-HIV-infected patients in northern Tanzania.

BACKGROUND: Tuberculosis (TB) is a major health problem and HIV is the major cause of the increase in TB. Sub-Saharan Africa is endemic for both TB and HIV infection. Determination of the prevalence of M. tuberculosis strains and their drug susceptibility is important for TB control.TB positive culture, BAL fluid or sputum samples from 130 patients were collected and genotyped. The spoligotypes were correlated with anti-tuberculous drug susceptibility in HIV-infected and non-HIV patients from Tanzania. RESULTS: One-third of patients were TB/HIV co-infected. Forty-seven spoligotypes were identified. Fourteen isolates (10.8%) had new and unique spoligotypes while 116 isolates (89.2%) belonged to 33 known spoligotypes. The major spoligotypes contained nine clusters: CAS1-Kili 30.0%, LAM11- ZWE 14.6%, ND 9.2%, EAI 6.2%, Beijing 5.4%, T-undefined 4.6%, CAS1-Delhi 3.8%, T1 3.8% and LAM9 3.8%. Twelve (10.8%) of the 111 phenotypically tested strains were resistant to anti-TB drugs. Eight (7.2%) were monoresistant strains: 7 to isoniazid (INH) and one to streptomycin. Four strains (3.5%) were resistant to multiple drugs: one (0.9%) was resistant to INH and streptomycin and the other three (2.7%) were MDR strains: one was resistant to INH, rifampicin and ethambutol and two were resistant to all four anti-TB drugs. Mutation in the katG gene codon 315 and the rpoB hotspot region showed a low and high sensitivity, respectively, as predictor of phenotypic drug resistance. CONCLUSION: CAS1-Kili and LAM11-ZWE were the most common families. Strains of the Beijing family and CAS1-Kili were not or least often associated with resistance, respectively. HIV status was not associated with spoligotypes, resistance or previous TB treatment.