Genetic barriers more than environmental associations explain Serratia marcescens population structure.

Bacterial species often comprise well-separated lineages, likely emerged and maintained by genetic isolation and/or ecological divergence. How these two evolutionary actors interact in the shaping of bacterial population structure is currently not fully understood. In this study, we investigate the genetic and ecological drivers underlying the evolution of Serratia marcescens, an opportunistic pathogen with high genomic flexibility and able to colonise diverse environments. Comparative genomic analyses reveal a population structure composed of five deeply-demarcated genetic clusters with open pan-genome but limited inter-cluster gene flow, partially explained by Restriction-Modification (R-M) systems incompatibility. Furthermore, a large-scale research on hundred-thousands metagenomic datasets reveals only a partial habitat separation of the clusters. Globally, two clusters only show a separate gene composition coherent with ecological adaptations. These results suggest that genetic isolation has preceded ecological adaptations in the shaping of the species diversity, an evolutionary scenario coherent with the Evolutionary Extended Synthesis.

Resistance patterns and transmission of mono- and polyresistant TB: clinical impact of WGS.

Objectives: Rapidly diagnosing drug-resistant TB is crucial for improving treatment and transmission control. WGS is becoming increasingly accessible and has added value to the diagnosis and treatment of TB. The aim of the study was to perform WGS to determine the rate of false-positive results of phenotypic drug susceptibility testing (pDST) and characterize the molecular mechanisms of resistance and transmission of mono- and polyresistant Mycobacterium (M.) tuberculosis. Methods: WGS was performed on 53 monoresistant and 25 polyresistant M. tuberculosis isolates characterized by pDST. Sequencing data were bioinformatically processed to infer mutations encoding resistance and determine the origin of resistance and phylogenetic relationship between isolates studied. Results: The data showed the variable sensitivity and specificity of WGS in comparison with pDST as the gold standard: isoniazid 92.7% and 92.3%; streptomycin 41.9% and 100.0%; pyrazinamide 15% and 94.8%; and ethambutol 75.0% and 98.6%, respectively. We found novel mutations encoding resistance to streptomycin (in gidB) and pyrazinamide (in kefB). Most isolates belonged to lineage 4 (80.1%) and the overall clustering rate was 11.5%. We observed lineage-specific gene variations encoding resistance to streptomycin and pyrazinamide. Conclusions: This study highlights the clinical potential of WGS in ruling out false-positive drug resistance following phenotypic or genetic drug testing, and recommend this technology together with the WHO catalogue in designing an optimal individualized treatment regimen and preventing the development of MDR TB. Our results suggest that resistance is primarily developed through spontaneous mutations or selective pressure.

Comprehensive Drug Resistance Characterization of Pulmonary Tuberculosis in Algeria: Insights on Mycobacterium tuberculosis Strains by Whole-Genome Sequencing.

In this study, we aimed to characterize drug-resistant strains by whole-genome sequencing (WGS), to describe the spreading lineages and the history of transmission. Drug susceptibility testing was performed by 96-well broth microdilution plates. The genomic DNA was extracted and purified; libraries were prepared and run on the Illumina NextSeq500 System. Among 82 isolates, 21 tuberculosis (TB) isolates (25.6%) were drug resistant, including 10 MDR and 4 pre-extensively drug-resistant (XDR)-TB. The mutation Ser315Thr in the katG gene was confirmed in 15 isolates. In rpoB, Ser450Leu and His445Asp mutations were the most common. Asp94Asn and Ala90Val mutations were reported in gyrA. The LAM family, the most TB drug resistant, was widely predominant in the north and the T sublineage in the south of the country. This study provides the first insight on TB drug resistance using WGS in Algeria and clearly describes the first pre-XDR-TB cases and lineage distribution across the country.

Bedaquiline and clofazimine resistance in Mycobacterium tuberculosis: an in-vitro and in-silico data analysis.

BACKGROUND: Bedaquiline is a core drug for the treatment of multidrug-resistant tuberculosis; however, the understanding of resistance mechanisms is poor, which is hampering rapid molecular diagnostics. Some bedaquiline-resistant mutants are also cross-resistant to clofazimine. To decipher bedaquiline and clofazimine resistance determinants, we combined experimental evolution, protein modelling, genome sequencing, and phenotypic data. METHODS: For this in-vitro and in-silico data analysis, we used a novel in-vitro evolutionary model using subinhibitory drug concentrations to select bedaquiline-resistant and clofazimine-resistant mutants. We determined bedaquiline and clofazimine minimum inhibitory concentrations and did Illumina and PacBio sequencing to characterise selected mutants and establish a mutation catalogue. This catalogue also includes phenotypic and genotypic data of a global collection of more than 14 000 clinical Mycobacterium tuberculosis complex isolates, and publicly available data. We investigated variants implicated in bedaquiline resistance by protein modelling and dynamic simulations. FINDINGS: We discerned 265 genomic variants implicated in bedaquiline resistance, with 250 (94%) variants affecting the transcriptional repressor (Rv0678) of the MmpS5-MmpL5 efflux system. We identified 40 new variants in vitro, and a new bedaquiline resistance mechanism caused by a large-scale genomic rearrangement. Additionally, we identified in vitro 15 (7%) of 208 mutations found in clinical bedaquiline-resistant isolates. From our in-vitro work, we detected 14 (16%) of 88 mutations so far identified as being associated with clofazimine resistance and also seen in clinically resistant strains, and catalogued 35 new mutations. Structural modelling of Rv0678 showed four major mechanisms of bedaquiline resistance: impaired DNA binding, reduction in protein stability, disruption of protein dimerisation, and alteration in affinity for its fatty acid ligand. INTERPRETATION: Our findings advance the understanding of drug resistance mechanisms in M tuberculosis complex strains. We have established an extended mutation catalogue, comprising variants implicated in resistance and susceptibility to bedaquiline and clofazimine. Our data emphasise that genotypic testing can delineate clinical isolates with borderline phenotypes, which is essential for the design of effective treatments. FUNDING: Leibniz ScienceCampus Evolutionary Medicine of the Lung, Deutsche Forschungsgemeinschaft, Research Training Group 2501 TransEvo, Rhodes Trust, Stanford University Medical Scientist Training Program, National Institute for Health and Care Research Oxford Biomedical Research Centre, Oxford University Hospitals NHS Foundation Trust, Bill & Melinda Gates Foundation, Wellcome Trust, and Marie Sklodowska-Curie Actions.

Advantages of long- and short-reads sequencing for the hybrid investigation of the Mycobacterium tuberculosis genome.

Introduction: In the fight to limit the global spread of antibiotic resistance, computational challenges associated with sequencing technology can impact the accuracy of downstream analysis, including drug resistance identification, transmission, and genome resolution. About 10% of Mycobacterium tuberculosis (MTB) genome is constituted by the PE/PPE family, a GC-rich repetitive genome region. Although sequencing using short read technology is widely used, it is well recognized its limit in the PE/PPE regions due to the unambiguously mapping process onto the reference genome. The aim of this study was to compare the performances of short-reads (SRS), long-reads (LRS) and hybrid-reads (HYBR) based analysis over different common investigative tasks: genome coverage estimation, variant calling and cluster analysis, drug resistance detection and de novo assembly. Methods: For the study 13 model MTB clinical isolates were sequenced with both SRS and LRS. HYBR were produced correcting the long reads with the short reads. The fastq from the three approaches were then processed using a customized version of MTBseq for genome coverage estimation and variant calling and using two different assemblers for de novo assembly evaluation. Results: Estimation of genome coverage performances showed lower 8X breadth coverage for SRS respect to LRS and HYBR: considering the PE/PPE genes, SRS showed low results for the PE_PGRS family, while obtained acceptable coverage in PE and PPE genes; LRS and HYBR reached optimal coverages in PE/PPE genes. For variant calling HYBR showed the highest resolution, detecting the highest percentage of uniquely identified mutations compared to LRS and SRS. All three approaches agreed on the identification of two major clusters, with HYBR identifying an higher number of SNPs between the two clusters. Comparing the quality of the assemblies, HYBR and LRS obtained better results than SRS. Discussion: In conclusion, depending on the aim of the investigation, both SRS and LRS present complementary advantages and limitations implying that for a full resolution of MTB genomes, where all the mentioned analyses and both technologies are needed, the use of the HYBR approach represents a valid option and a well-rounded strategy.

Whole genome sequencing of Bulgarian rifampicin resistant Mycobacterium tuberculosis strains.

INTRODUCTION: The transmission of drug-resistant tuberculosis is one of the greatest challenges facing the global tuberculosis control. AIM: The aim of the study was to investigate the resent transmission of rifampicin resistant tuberculosis in Bulgaria and to describe the mutations related to the antimicrobials' resistance using whole genome sequencing. MATERIALS AND METHODS: As part of an ECDC funded pilot study for evaluation of the systematic use of whole genome sequencing (WGS) of Mycobacterium tuberculosis (MTB) surveillance (EUSeqMyTB), Bulgaria provided 65 rifampicin resistant isolates over a three years' timeframe (2017-2019) representing 87.5% of the notified rifampicin resistant cases. Drug resistance prediction and relatedness analysis of the resistant isolates was performed in collaboration with San Raffaele Scientific Institute, Milan, Italy. RESULTS: Almost all of the isolates were identified as Euro-American lineage (96.9%); 18.5% of the isolates were found to be resistant to fluoroquinolones, but no mutations conferring resistance to bedaquiline or linezolid could be identified. Less than half (43.3%) of the isolates were clustered (<5 SNPs distance) into a total of seven national SNP-based clusters, while a total of six isolates were found to be part of different cross-border clusters. All clustered cases originated from Bulgaria. CONCLUSIONS: WGS has proven to be a reliable tool for surveillance and tracing of recent transmission of tuberculosis and has the potential for resistance prediction for most of the antituberculosis drugs.

Whole genome sequencing of multidrug-resistant Mycobacterium tuberculosis isolates collected in the Czech Republic, 2005-2020.

The emergence and spread of resistant tuberculosis (TB) pose a threat to public health, so it is necessary to diagnose the drug-resistant forms in a clinically short time frame and closely monitor their transmission. In this study, we carried out a first whole genome sequencing (WGS)-based analysis of multidrug resistant (MDR) M. tuberculosis strains to explore the phylogenetic lineages diversity, drug resistance mechanisms, and ongoing transmission chains within the country. In total, 65 isolates phenotypically resistant to at least rifampicin and isoniazid collected in the Czech Republic in 2005-2020 were enrolled for further analysis. The agreement of the results obtained by WGS with phenotypic drug susceptibility testing (pDST) in the determination of resistance to isoniazid, rifampicin, pyrazinamide, streptomycin, second-line injectables and fluoroquinolones was more than 80%. Phylogenetic analysis of WGS data revealed that the majority of MDR M. tuberculosis isolates were the Beijing lineage 2.2.1 (n = 46/65; 70.8%), while the remaining strains belonged to Euro-American lineage. Cluster analysis with a predefined cut-off distance of less than 12 single nucleotide polymorphisms between isolates showed 19 isolates in 6 clusters (clustering rate 29.2%), located mainly in the region of the capital city of Prague. This study highlights the utility of WGS as a high-resolution approach in the diagnosis, characterization of resistance patterns, and molecular-epidemiological analysis of resistant TB in the country.

Rapid SARS-CoV-2 Intra-Host and Within-Household Emergence of Novel Haplotypes.

In February 2020, the municipality of Vo', a small town near Padua (Italy) was quarantined due to the first coronavirus disease 19 (COVID-19)-related death detected in Italy. To investigate the viral prevalence and clinical features, the entire population was swab tested in two sequential surveys. Here we report the analysis of 87 viral genomes, which revealed that the unique ancestor haplotype introduced in Vo' belongs to lineage B, carrying the mutations G11083T and G26144T. The viral sequences allowed us to investigate the viral evolution while being transmitted within and across households and the effectiveness of the non-pharmaceutical interventions implemented in Vo'. We report, for the first time, evidence that novel viral haplotypes can naturally arise intra-host within an interval as short as two weeks, in approximately 30% of the infected individuals, regardless of symptom severity or immune system deficiencies. Moreover, both phylogenetic and minimum spanning network analyses converge on the hypothesis that the viral sequences evolved from a unique common ancestor haplotype that was carried by an index case. The lockdown extinguished both the viral spread and the emergence of new variants.

Whole genome sequencing of isoniazid monoresistant clinical isolates of Mycobacterium tuberculosis reveals novel genetic polymorphisms.

In an attempt to uncover genotypic indicators for isoniazid (INH) resistance in M. tuberculosis, in addition to the canonical mutations in genes associated with INH resistance, including katG, inhA and fabG promoter; we analyzed, two INH monoresistant isolates, ASTS24/13 (INHR1) and SHR1/14 (INHR2). Targeted Sanger sequencing detected a canonical mutation at katG315 only in INHR2. Infection of THP-1 cells and exposure to antituberculosis drugs led to two-fold increase in the minimum inhibitory concentration of INH in INHR2. Whole genome sequences revealed that INHR1 and INHR2 belonged to Delhi Central Asian Strain and East African Indian lineages, respectively. The sequences were compared with INH susceptible isolates with the same lineage as the INH monoresistant strains. INHR1 had a novel unique mutation STOP420Trp in the efflux pump gene Rv0849, while INHR2 had a novel mutation Arg579Ser in efflux pump gene mmpL5. Comparison of lipid associated genes showed novel mutations in INHR1 in fadE16, fadD3 and fbpD; while INHR2 had mutations in fadE1, Rv0145, Rv1425, fadD9 and mmaA3. Both isolates also demonstrated novel mutations in cell wall associated genes. Our study highlights the importance of searching for alternate mechanisms of INH resistance that may contribute to the development of more comprehensive diagnostic tools.

Detection of NDM-1/5 and OXA-48 co-producing extensively drug-resistant hypervirulent Klebsiella pneumoniae in Northern Italy.

OBJECTIVES: Using a hybrid long-read sequencing approach, we aimed to fully characterise four extensively-drug resistant (XDR) hypervirulent Klebsiella pneumoniae isolates, one of which represented the first strain isolated in Italy co-expressing NDM-1/5 and OXA-48 carbapenemases. METHODS: Whole-genome sequencing was performed using Illumina and Oxford Nanopore Technology platforms. An assembly pipeline was used to recover the structures both of the chromosome and plasmids. RESULTS: Multilocus sequence typing (MLST) showed that these strains belonged to high-risk sequence types (STs) not commonly circulating in Italy (ST383, ST147 and ST15). The hybrid sequencing approach allowed to characterise three multidrug resistance plasmids, which demonstrated high homology with previously sequenced plasmids, that were simultaneously detected in one ST383 strain carrying, respectively, blaNDM-1, blaNDM-5 and blaOXA-48. CONCLUSION: This is the first report in Italy of new hypervirulent XDR K. pneumoniae clones characterised by co-production of OXA-48, NDM-1 and NDM-5. The discovery of new high-risk clones harbouring multiple mobile elements is a growing problem that poses a great challenge for public health.

Anti-tuberculosis drug resistance in Slovakia, 2018-2019: The first whole-genome epidemiological study.

OBJECTIVE: The resistance of Mycobacterium (M.) tuberculosis to antituberculosis drugs poses a major threat to global public health. Whole genome sequencing (WGS) is an increasingly preferred method in the diagnostics and monitoring of the transmission dynamics of resistant forms of tuberculosis (TB). The aim of the study was to, for the first time, use the sequencing-based analysis to study the transmission and resistance patterns of a systematic and recent collection of extensively drug resistant (XDR) and multidrug resistant tuberculosis (MDR-TB) isolates and to expand our knowledge about drug resistant (DR) TB epidemiological dynamics in Slovakia. DESIGN: A total of 495 patients with pulmonary TB, who were referred to National Reference Laboratory for Mycobacteriology (Vysné Hágy, Slovakia) in the years 2018-2019, were studied. Out of the total of 495 patients, 4 XDR-TB (0.8%) and 8 (1.6%) MDR-TB isolates were identified by conventional drug susceptibility testing on Löwenstein-Jensen solid medium and subjected to whole genome sequencing. Sequencing data were evaluated for molecular-epidemiological analysis and identification of resistance patterns. RESULTS: Phylogenetic and cluster analysis showed extensive recent transmission events and the predominance of Euro-American lineage 4.7 in Slovakia. However, phylogenetic analysis revealed the circulation of several lineages that originally occurred in Eastern European countries. Resistance patterns for first- and second-line antituberculosis drugs characterized by whole genome sequencing were in high concordance with the results of phenotypic drug susceptibility testing. CONCLUSION: Forty percent of at least MDR-TB isolates were not genetically linked, indicating that appropriate measures should be taken to monitor and prevent the spread of drug-resistant tuberculosis within the country as well as in other regions.

Enhanced Molecular Dynamics Method to Efficiently Increase the Discrimination Capability of Computational Protein-Protein Docking.

Protein-protein docking typically consists of the generation of putative binding conformations, which are subsequently ranked by fast heuristic scoring functions. The simplicity of these functions allows for computational efficiency but has severe repercussions on their discrimination capabilities. In this work, we show the effectiveness of suitable descriptors calculated along short scaled molecular dynamics runs in recognizing the nearest-native bound conformation among a set of putative structures generated by the HADDOCK tool for eight protein-protein systems.

Role of Epistasis in Amikacin, Kanamycin, Bedaquiline, and Clofazimine Resistance in Mycobacterium tuberculosis Complex.

Antibiotic resistance among bacterial pathogens poses a major global health threat. Mycobacterium tuberculosis complex (MTBC) is estimated to have the highest resistance rates of any pathogen globally. Given the low growth rate and the need for a biosafety level 3 laboratory, the only realistic avenue to scale up drug susceptibility testing (DST) for this pathogen is to rely on genotypic techniques. This raises the fundamental question of whether a mutation is a reliable surrogate for phenotypic resistance or whether the presence of a second mutation can completely counteract its effect, resulting in major diagnostic errors (i.e., systematic false resistance results). To date, such epistatic interactions have only been reported for streptomycin that is now rarely used. By analyzing more than 31,000 MTBC genomes, we demonstrated that the eis C-14T promoter mutation, which is interrogated by several genotypic DST assays endorsed by the World Health Organization, cannot confer resistance to amikacin and kanamycin if it coincides with loss-of-function (LoF) mutations in the coding region of eis. To our knowledge, this represents the first definitive example of antibiotic reversion in MTBC. Moreover, we raise the possibility that mmpR (Rv0678) mutations are not valid markers of resistance to bedaquiline and clofazimine if these coincide with an LoF mutation in the efflux pump encoded by mmpS5 (Rv0677c) and mmpL5 (Rv0676c).

A Simplified Amino Acidic Alphabet to Unveil the T-Cells Receptors Antigens: A Computational Perspective.

The exposure to pathogens triggers the activation of adaptive immune responses through antigens bound to surface receptors of antigen presenting cells (APCs). T cell receptors (TCR) are responsible for initiating the immune response through their physical direct interaction with antigen-bound receptors on the APCs surface. The study of T cell interactions with antigens is considered of crucial importance for the comprehension of the role of immune responses in cancer growth and for the subsequent design of immunomodulating anticancer drugs. RNA sequencing experiments performed on T cells represented a major breakthrough for this branch of experimental molecular biology. Apart from the gene expression levels, the hypervariable CDR3α/ß sequences of the TCR loops can now be easily determined and modelled in the three dimensions, being the portions of TCR mainly responsible for the interaction with APC receptors. The most direct experimental method for the investigation of antigens would be based on peptide libraries, but their huge combinatorial nature, size, cost, and the difficulty of experimental fine tuning makes this approach complicated time consuming, and costly. We have implemented in silico methodology with the aim of moving from CDR3α/ß sequences to a library of potentially antigenic peptides that can be used in immunologically oriented experiments to study T cells' reactivity. To reduce the size of the library, we have verified the reproducibility of experimental benchmarks using the permutation of only six residues that can be considered representative of all ensembles of 20 natural amino acids. Such a simplified alphabet is able to correctly find the poses and chemical nature of original antigens within a small subset of ligands of potential interest. The newly generated library would have the advantage of leading to potentially antigenic ligands that would contribute to a better understanding of the chemical nature of TCR-antigen interactions. This step is crucial in the design of immunomodulators targeted towards T-cells response as well as in understanding the first principles of an immune response in several diseases, from cancer to autoimmune disorders.

Tuning Local Hydration Enables a Deeper Understanding of Protein-Ligand Binding: The PP1-Src Kinase Case.

Water plays a key role in biomolecular recognition and binding. Despite the development of several computational and experimental approaches, it is still challenging to comprehensively characterize water-mediated effects on the binding process. Here, we investigate how water affects the binding of Src kinase to one of its inhibitors, PP1. Src kinase is a target for treating several diseases, including cancer. We use biased molecular dynamics simulations, where the hydration of predetermined regions is tuned at will. This computational technique efficiently accelerates the SRC-PP1 binding simulation and allows us to identify several key and yet unexplored aspects of the solvent's role. This study provides a further perspective on the binding phenomenon, which may advance the current drug design approaches for the development of new kinase inhibitors.

Application of Targeted Next-Generation Sequencing Assay on a Portable Sequencing Platform for Culture-Free Detection of Drug-Resistant Tuberculosis from Clinical Samples.

Targeted next-generation sequencing (tNGS) has emerged as a comprehensive alternative to existing methods for drug susceptibility testing (DST) of Mycobacterium tuberculosis from patient sputum samples for clinical diagnosis of drug-resistant tuberculosis (DR-TB). However, the complexity of sequencing platforms has limited their uptake in low-resource settings. The goal of this study was to evaluate the use of the tNGS-based DST solution Genoscreen Deeplex Myc-TB, for use on the compact, low-cost Oxford Nanopore Technologies MinION sequencer. One hundred four DNA samples extracted from smear-positive sputum sediments, previously sequenced using the Deeplex assay on an Illumina MiniSeq, were resequenced on MinION after applying a custom library preparation. MinION read quality, mapping statistics, and variant calling were computed using an in-house pipeline and compared to the reference MiniSeq data. The average percentage of MinION reads mapped to an H37RV reference genome was 90.8%, versus 99.5% on MiniSeq. The mean depths of coverage were 4,151× and 4,177× on MinION and MiniSeq, respectively, with heterogeneous distribution across targeted genes. Composite reference coverage breadth was >99% for both platforms. We observed full concordance between technologies in reporting the clinically relevant drug-resistant markers, including full gene deletions. In conclusion, we demonstrated that the workflow and sequencing data obtained from Deeplex on MinION are comparable to those for the MiniSeq, despite the higher raw error rates on MinION, with the added advantage of MinION's portability, versatility, and low capital costs. Targeted NGS on MinION is a promising DST solution for rapidly providing clinically relevant data to manage complex DR-TB cases.

NanoShaper-VMD interface: computing and visualizing surfaces, pockets and channels in molecular systems.

SUMMARY: NanoShaper is a program specifically aiming the construction and analysis of the molecular surface of nanoscopic systems. It uses ray-casting for parallelism and it performs analytical computations whenever possible to maximize robustness and accuracy of the approach. Among the other features, NanoShaper provides volume, surface area, including that of internal cavities, for any considered molecular system. It identifies pockets via a very intuitive definition based on the concept of probe radius, intrinsic to the definition of the solvent excluded surface. We show here that, with a suitable choice of the parameters, the same approach can also permit the visualisation of molecular channels. NanoShaper has now been interfaced with the widely used molecular visualization software VMD, further enriching its already well furnished toolset. AVAILABILITY AND IMPLEMENTATION: VMD is available at http://www.ks.uiuc.edu/Research/vmd/. NanoShaper, its documentation, tutorials and supporting programs are available at http://concept.iit.it/downloads. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.